Mapping and characterization of the DQ subregion of the ovine MHC

A map of the ovine MHC class II DQ subregion has been constructed from overlapping cosmid clones. This region consists of two loci linked on a linear tract of 130 kb DNA. Each locus consists of a DQA and a DQB gene in a tail-to-tail orientation. The genes in each locus are transcribed but only those designated DQ1 express class II molecules at the surface of mouse L cells following DNA-mediated gene transfection. The DQA1 and DQB1 genes are separated by 11kb while the DQA2 and B2 genes are 25 kb apart. The loci are separated by 22 kb.     

Diversity analysis at MHC class II DQA locus in buffalo (Bubalus bubalis) indicates extensive duplication and trans-species evolution

Variation at MHC Class II-DQA locus in riverine and swamp buffaloes (Bubu) has been explored in this study. Through sequencing of buffalo DQA, 48 nucleotide variants identified from 17 individuals, reporting 42 novel alleles, including one pseudogene. Individual animal displayed two to seven variants, suggesting the presence of more than two Bubu-DQA loci, as an evidence of extensive duplication. dN values were found to be higher than dS values at peptide binding sites, separately for riverine and swamp buffaloes, indicating locus being under positive selection. Evolutionary analysis revealed numerous trans-species polymorphism with alleles from water buffalo assigned to at least three different loci (Bubu-DQA1, DQA2, DQA3). Alleles of both the sub-species intermixed within the cluster, showing convergent evolution of MHC alleles in bovines. The results thus suggest that both riverine and swamp buffaloes share con-current arrangement of DQA region, comparable to cattle in terms of copy number and population polymorphism.     

Major histocompatibility complex class II genetic variation in forest musk deer (Moschus berezovskii) in China

The major histocompatibility complex (MHC) plays an important role in the immune system of vertebrates. We used the second exon of four MHC class II genes (DRA, DQA1, DQA2 and DRB3) to assess the overall MHC variation in forest musk deer (Moschus berezovskii). We also compared the MHC variation in captive and wild populations. We observed 22 alleles at four loci (four at DRA, four at DQA1, four at DQA2 and 10 at DRB3), 15 of which were newly identified alleles. Results suggest that forest musk deer maintain relatively high MHC variation, which may result from balancing selection. Moreover, considerable diversity was observed at the DRA locus. We found a high frequency of Mobe‐DRA*02, Mobe‐DQA1*01 and Mobe‐DQA2*05 alleles, which may be important for pathogen resistance. A Ewens–Watterson test showed that the DRB3 locus in the wild population had experienced recent balancing selection. We detected a small divergence at the DRA locus, suggesting the effect of weak positive selection on the DRA gene. Alternatively, this locus may be young and not yet adapted a wide spectrum of alleles for pathogen resistance. The significant heterozygosity deficit observed at the DQA1 and DRB3 loci in the captive population and at all four loci in the wild population may be the result of a population bottleneck. Additionally, MHC genetic diversity was higher in the wild population than in the captive, suggesting that the wild population may have the ability to respond to a wider range of pathogens.     

Analysis of genetic diversity at the DQA loci in African cattle: evidence for a BoLA-DQA3 locus

 We describe the development of a polymerase chain reaction (PCR)-based approach for analysis of genetic diversity at the DQA loci in African Bos indicus and Bos taurus cattle. This approach, equally effective in European and Asian cattle breeds, detects the presence or absence of DQA1 and most duplicated DQA2 genes. Nucleotide and predicted amino acid sequence analysis of the highly polymorphic second exons, in addition to analysis of the locus-specific and relatively non-polymorphic transmembrane, cytoplasmic, and 3-prime untranslated regions, has provided evidence for considerable diversity between each of the duplicated DQA2 genes. Therefore, we propose the designation BoLA-DQA3 for the previously unpublished alleles at the second DQA2 locus. Fourteen distinct PCR restriction fragment length polymorphism (RFLP) patterns, each identifying families of alleles at three DQA loci, can be distinguished. Nucleotide sequence analysis of new PCR-RFLP patterns from 193 Kenyan Boran, Ethiopian Arsi (B. indicus), and Guinean N’Dama (B. taurus) cattle identified 13 DQA1 alleles within eight major allelic families, five DQA2 alleles within a single allelic family, and seven DQA3 alleles within three major allelic families. Received: 19 February 1997 / Revised: 28 February 1997     

Major histocompatibility complex variation and evolution at a single, expressed DQA locus in two genera of elephants

Genes of the vertebrate major histocompatibility complex (MHC) are crucial to defense against infectious disease, provide an important measure of functional genetic diversity, and have been implicated in mate choice and kin recognition. As a result, MHC loci have been characterized for a number of vertebrate species, especially mammals; however, elephants are a notable exception. Our study is the first to characterize patterns of genetic diversity and natural selection in the elephant MHC. We did so using DNA sequences from a single, expressed DQA locus in elephants. We characterized six alleles in 30 African elephants (Loxodonta africana) and four alleles in three Asian elephants (Elephas maximus). In addition, for two of the African alleles and three of the Asian alleles, we characterized complete coding sequences (exons 1–5) and nearly complete non-coding sequences (introns 2–4) for the class II DQA loci. Compared to DQA in other wild mammals, we found moderate polymorphism and allelic diversity and similar patterns of selection; patterns of non-synonymous and synonymous substitutions were consistent with balancing selection acting on the peptides involved in antigen binding in the second exon. In addition, balancing selection has led to strong trans-species allelism that has maintained multiple allelic lineages across both genera of extant elephants for at least 6 million years. We discuss our results in the context of MHC diversity in other mammals and patterns of evolution in elephants.     

Molecular characterization of swine leukocyte antigen gene diversity in purebred Pietrain pigs

The porcine major histocompatibility complex (MHC) harbors the highly polymorphic swine leukocyte antigen (SLA) class I and II gene clusters encoding glycoproteins that present antigenic peptides to T cells in the adaptive immune response. In Austria, the majority of commercial pigs are F ₂ descendants of F ₁ Large White/Landrace hybrids paired with Pietrain boars. Therefore, the repertoire of SLA alleles and haplotypes present in Pietrain pigs has an important influence on that of their descendants. In this study, we characterized the SLA class I ( SLA‐1 , SLA‐2 , SLA‐3 ) and class II ( SLA‐DRB1 , SLA‐DQB1 , SLA‐DQA ) genes of 27 purebred Pietrain pigs using a combination of the high‐resolution sequence‐based typing (SBT) method and a low‐resolution (Lr) PCR‐based method using allele‐group, sequence‐specific primers (PCR‐SSP). A total of 15 class I and 13 class II haplotypes were identified in the studied cohort. The most common SLA class I haplotype Lr‐43.0 ( SLA‐1 *11XX– SLA‐3 *04XX– SLA‐2 *04XX) was identified in 11 animals with a frequency of 20%. For SLA class II, the most prevalent haplotype, Lr‐0.14 [ SLA‐DRB1 *0901– SLA‐DQB1 *0801– SLA‐DQA *03XX], was found in 14 animals with a frequency of 26%. Two class II haplotypes, tentatively designated as Lr‐Pie‐0.1 [ SLA‐DRB1 *01XX/be01/ha04– SLA‐DQB1 *05XX– SLA ‐ DQA*blank] and Lr‐Pie‐0.2 [ SLA‐DRB1 *06XX– SLA‐DQB1 *03XX– SLA‐DQA *03XX], appeared to be novel and have never been reported so far in other pig populations. We showed that SLA genotyping using PCR‐SSP‐based assays represents a rapid and cost‐effective way to study SLA diversity in outbred commercial pigs and may facilitate the development of more effective vaccines or identification of disease‐resistant pigs in the context of SLA antigens to improve overall swine health.     

Diversity and locus specificity of chicken MHC B class I sequences

The major histocompatibility complex B (MHC B) region in a standard haplotype of Leghorn chickens contains two closely linked class I loci, B-FI and B-FIV. Few sequences of B-FI alleles are available, and therefore alleles of the two loci have not been compared with regard to sequence diversity or locus specificity. Here, we report eight new B-F alpha 1/alpha 2-coding sequences from broiler chicken MHC B haplotypes, and a unique recombinant between the two B-F loci. The new sequences were combined with existing B-F sequences from Leghorn and broiler haplotypes for analysis. On the basis of phylogenetic analysis and conserved sequence motifs, B-F sequences separated into two groups (Groups A and B), corresponding to B-FIV and B-FI locus, respectively. Every broiler haplotype had one B-F sequence in Group A and the second B-F sequence, if it existed, clustered in Group B. Group B (presumptive B-FI locus) sequences identified in broiler haplotypes resembled the human MHC class I HLA-C locus in their distinctive pattern of allelic polymorphism. Compared with B-FIV, B-FI alleles were less polymorphic and possessed a conserved locus-specific motif in the alpha1 helix, but nevertheless demonstrated evidence of diversifying selection. One B-FI alpha 1/alpha 2-coding nucleotide sequence was completely conserved in four different broiler haplotypes, but each allele differed in the exon encoding the alpha 3 domain.     

Genetic variation in the major histocompatibility complex of the European brown hare (Lepus europaeus) across distinct phylogeographic areas

The major histocompatibility complex is one of the best studied systems in vertebrates providing evidence for the long-term action of selection. Here, we examined the intra- and inter-population genetic diversity of the MHC class II DRB locus in European brown hare (Lepus europaeus) and correlated the results with genetic variability already estimated from the MHC DQA locus and from maternally (mitochondrial DNA (mtDNA)) and biparentally (allozymes, microsatellites) inherited loci. L. europaeus showed remarkable genetic polymorphism in both DQA and DRB1 loci. The Anatolian populations exhibited the highest genetic polymorphism for both loci. Balancing selection has established increased variability in the European populations despite the founder effects after the last glaciation. Different evolutionary rates were traced for DRB1 and DQA loci, as evidenced by the higher number of common DRB1 than DQA alleles and the greater differences between DRB1 alleles with common origin in comparison with DQA alleles. The high number of rare alleles with low frequencies detected implies that frequency-dependent selection drives MHC evolution in the brown hare through the advantage of rare alleles. Both loci were under the influence of positive selection within the peptide-binding region. The functional polymorphism, recorded as amino acid substitutions within the binding pockets, fell also within distinct geographic patterns, yet it was much narrower than the genetic polymorphism. We hypothesize that certain structural and functional characteristics of the binding pockets set limitations to the actual shape of genetic polymorphism in MHC.     

Allelic diversity at the Mhc-DQA locus of woodmouse populations (Apodemus sylvaticus) present in the islands and mainland of the northern Mediterranean

Aim Polymorphism at neutral markers and at MHC loci in rodent populations living on islands is generally low. The main genetic factors that may contribute to a reduced level of genetic variability are genetic drift, reduced gene flow and founder events. We investigated the pattern of polymorphism at the second exon of the Mhc‐DQA gene in island populations of Apodemus sylvaticus and in their mainland counterparts to investigate the pattern of MHC polymorphism in a phylogeographical context and to assess the impact of insularity on diversity at this locus. Location Eight north Mediterranean populations of Apodemus sylvaticus were studied, including five island populations (Majorca, Minorca, Porquerolles, Port‐Cros and Sicily) and three mainland populations. Methods cDNA sequencing and nucleotide sequences analyses. Synonymous and non‐synonymous substitutions were examined at the PBR and non‐PBR sites. The DQA allelic distribution in populations was compared with the woodmouse phylogeography. Results This study presents novel DQA alleles. High polymorphism of the DQA locus is recorded in natural populations of A. sylvaticus with 13 alleles being widely distributed irrespective of the geographical origin and palaeoclimatic history of populations. The DQA locus does not show the expected pattern for non‐synonymous substitutions at the PBR sites. However, island populations show a weak loss of polymorphism in comparison with their mainland counterparts. Main conclusions The DQA locus in the woodmouse seems to be subject to weak selection and does not allow resolution of phylogeographical relationships among European woodmouse populations. The presence of at least three alleles in island populations and the maintenance of five alleles between the two European lineages over 1.5 Myr suggest that balancing selection may act within populations, and more precisely within island populations, to maintain genetic variability. This study shows that phylogeographical studies are a prerequisite for any genetic investigation of selected genes in natural populations.     

Effect of swine leukocyte antigen‐DQA gene variation on diarrhea in Large White,Landrace, and Duroc piglets

Piglet diarrhea is one of the most common factors that affects the benefits of the swine industry. Although recent studies have shown that exon 2 of SLA‐DQA is associated with piglet resistance to diarrhea, contributions of genetic variation in the additional exon coding regions of this gene remain unclear. Here, we investigated variation in exons 1, 3 and 4 of the SLA‐DQA gene and evaluated their effects on diarrheal infection in 425 suckling piglets. No variation was identified in exon 1. In exon 3, there were eight alleles detected, generated by 14 single nucleotide polymorphisms (SNPs) and three nucleotide deletions, eight SNPs being newly identified. Four allele sequences and three SNPs were identified in exon 4, only one SNP being newly identified. Statistical analysis showed that the genotypes of exon 3 are significantly associated with piglet diarrhea; indeed, genotypes DQA*wb01/wb02 and wb04/wb05 are clearly associated with resistance to piglet diarrhea, as they have the lowest probabilities of infection (P < 0.05). However, no significant association was found between the genotypes of exon 4 and diarrhea (P > 0.05). These results provide important new information concerning the level of genetic diversity at the SLA‐DQA locus and suggest that further genetic association studies of piglet diarrhea resistance should include analyses of both exons 2 and 3 of this locus.     

Molecular characterization of swine leucocyte antigen class II genes in outbred pig populations

The highly polymorphic swine leucocyte antigen (SLA) genes are among the most important determinants of swine immune responses to disease and vaccines. Accurate and effective SLA genotyping methods are required to understand how SLA gene polymorphisms affect immunity, especially in outbred pigs with diverse genetic backgrounds. In this study, we present a simple and rapid molecular‐based typing system for characterizing SLA class II alleles of the DRB1, DQB1 and DQA loci. This system utilizes a set of 47 sequence‐specific PCR primers developed to differentiate alleles by groups that share similar sequence motifs. We applied this typing method to investigate the SLA class II diversity in four populations of outbred pigs (n = 206) and characterized a total of 19 SLA class II haplotypes, six of which were shared by at least three of the sampled pig populations. We found that Lr‐0.1 (DRB1*01XX–DQB1*01XX–DQA*01XX) was the most prevalent haplotype with a combined frequency of 16.0%, followed by Lr‐0.2 (DRB1*02XX–DQB1*02XX–DQA*02XX) with 14.6% and Lr‐0.15b (DRB1*04XX–DQB1*0202–DQA*02XX) with 14.1%. Over 70% of the pigs (n = 147) had at least one copy of one of these three haplotypes. The PCR‐based typing system described in this study demonstrates a reliable and unambiguous detection method for SLA class II alleles. It will be a valuable tool for studying the influence of SLA diversity on various immunological, pathological and physiological traits in outbred pigs.     

Structure and evolution of the promoter regions of the DQA genes

HLA-DQ antigens are unique among the class II antigens in that their chains are highly polymorphic. In the present study, we characterized the general structure of the promoter regions of the DQA genes derived from different DR haplotypes and defined their nucleotide sequence polymorphisms. The promoter of each DQA1 allele contains three sequence motifs which are not present in non-DQA related class II genes: one identical to a tumor necrosis factor (TNF) response element, one similar to an NFB binding element, and one similar to a W motif. All DQA alleles lack TATA and CCAAT boxes in the proximal promoter region but carry other sequence elements characteristic of MHC class II genes, including S, X, X2, and Y boxes, and a pyrimidine-rich tract upstream of the X box. Nucleotide sequence polymorphisms among the various DQA1 alleles were noted within the promoter region, with some of the differences mapping within, or close to, regulatory elements that are important for the expression of MHC class II genes. All DQA1 alleles carry an unrearranged, full length, Alu-Sx related repeat immediately upstream of the proximal promoter region. This repeat was not present in the DQA2 (DXA) genes analyzed, confirming that DQ locus duplication probably occurred before integration of the Alu repeat into the primordial DQA1 locus, some 31–43 million years (myr) ago. The DQA2 promoter region is highly conserved between DR4 and DR3 haplotypes, with the degree of conservation exceeding that expected from the neutral mutation rate.The nucleotide sequence data reported in this paper have been submitted to the GenBank nucleotide sequence data base and have been assigned the accession numbers M97 454-M97 464. Correspondence to: E. Morzycka-Wroblewska.     

Homozygosity at a class II MHC locus depresses female reproductive ability in European brown hares

The link between adaptive genetic variation, individual fitness and wildlife population dynamics is fundamental to the study of ecology and evolutionary biology. In this study, a Bayesian modelling approach was employed to examine whether individual variability at two major histocompatibility complex (MHC) class II loci (DQA and DRB) and eight neutral microsatellite loci explained variation in female reproductive success for wild populations of European brown hare (Lepus europaeus). We examined two aspects of reproduction: the ability to reproduce (sterility) and the number of offspring produced (fecundity). Samples were collected from eastern Austria, experiencing a sub‐continental climatic regime, and from Belgium with a more Atlantic‐influenced climate. As expected, reproductive success (both sterility and fecundity) was significantly influenced by age regardless of sampling locality. For Belgium, there was also a significant effect of DQA heterozygosity in determining whether females were able to reproduce (95% highest posterior density interval of the regression parameter [−3.64, −0.52]), but no corresponding effect was found for Austria. In neither region was reproduction significantly associated with heterozygosity at the DRB locus. DQA heterozygotes from both regions also showed a clear tendency, but not significantly so, to produce a larger number of offspring. Predictive simulations showed that, in Belgium, sub‐populations of homozygotes will have higher rates of sterile individuals and lower average offspring numbers than heterozygotes. No similar effect is predicted for Austria. The mechanism for the spatial MHC effect is likely to be connected to mate choice for increased heterozygosity or to the linkage of certain MHC alleles with lethal recessives at other loci.     

Evolution of the ovine MHC DQA region

Southern hybridisation was used to define an apparent gene duplication event at the ovine DQA2 locus. Approximately 500 sheep from five different breeds were genotyped at their DQA1 and DQA2 loci. A subset of these were selected for further characterisation. Southern hybridisation of TaqI digested DNA revealed no DQA1 region in some sheep. It was also noted in these DQA1 null animals the DQA2 specific probe hybridised to two bands. An EcoRV-RFLP designed to distinguish copy number confirmed this duplication of the DQA2 region. The results showed that the duplication was exclusively associated with the DQA1 null haplotype and occurred only in alleles DQA2-F, -G, -I and -J. Comparison with bovine MHC genes revealed that they also contained a DQA1 null haplotype and that this haplotype was associated with a putative DQA3 gene. The potential for an ovine DQA3 locus is discussed.     

High MHC gene copy number maintains diversity despite homozygosity in a Critically Endangered single‐island endemic bird,but no evidence of MHC‐based mate choice

Small population sizes can, over time, put species at risk due to the loss of genetic variation and the deleterious effects of inbreeding. Losing diversity in the major histocompatibility complex (MHC) could be particularly harmful, given its key role in the immune system. Here, we assess MHC class I (MHC‐I) diversity and its effects on mate choice and survival in the Critically Endangered Raso lark Alauda razae, a species restricted to the 7 km² islet of Raso, Cape Verde, since ~1460, whose population size has dropped as low as 20 pairs. Exhaustively genotyping 122 individuals, we find no effect of MHC‐I genotype/diversity on mate choice or survival. However, we demonstrate that MHC‐I diversity has been maintained through extreme bottlenecks by retention of a high number of gene copies (at least 14), aided by cosegregation of multiple haplotypes comprising 2–8 linked MHC‐I loci. Within‐locus homozygosity is high, contributing to low population‐wide diversity. Conversely, each individual had comparably many alleles, 6–16 (average 11), and the large and divergent haplotypes occur at high frequency in the population, resulting in high within‐individual MHC‐I diversity. This functional immune gene diversity will be of critical importance for this highly threatened species’ adaptive potential.     

Evolutionary relationships among the primate Mhc-DQA1 and DQA2 alleles

The variation of the Mhc-DQA1 and DQA2 loci of ten different primate species (hominoids and Old World monkeys) was studied in order to obtain an insight in the processes that generate polymorphism of major histocompatibility complex (Mhc) class II genes and to establish the evolutionary relationships of their alleles. To that end nucleotide sequences of 36 Mhc class II DQA1 and seven DQA2 second exons were determined and phylogenetic trees that illustrate their evolutionary relationships were constructed. We demonstrate the existence of four primate Mhc-DQA1 allele lineages, two of which probably existed before the separation of the ancestors of the hominoids and Old World monkeys (approximately 22–28 million years ago). Mhc-DQA2 sequences were found only in the hominoid species and showed little diversity. We found no evidence for a major contribution of recombinational events to the generation of allelic diversity of the primate Mhc-DQA1 locus. Instead, our data suggest that the primate Mhc-DQA1 and DQA2 loci are relatively stable entities that mutated primarily as a result of point mutations.The nucleotide sequence data reported in this paper have been submitted to the GenBank nucleotide sequence database and have been assigned the accession numbers M76186-M76229.     

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.     

Balancing selection and heterozygote advantage in major histocompatibility complex loci of the bottlenecked Finnish wolf population

Maintaining effective immune response is an essential factor in the survival of small populations. One of the most important immune gene regions is the highly polymorphic major histocompatibility complex (MHC). We investigated how a population bottleneck and recovery have influenced the diversity and selection in three MHC class II loci, DLA‐DRB1, DLA‐DQA1 and DLA‐DQB1, in the Finnish wolf population. We studied the larger Russian Karelian wolf population for comparison and used 17 microsatellite markers as reference loci. The Finnish and Karelian wolf populations did not differ substantially in their MHC diversities ( = 0.047, P = 0.377), but differed in neutral microsatellite diversities ( = 0.148, P = 0.008). MHC allele frequency distributions in the Finnish population were more even than expected under neutrality, implying balancing selection. In addition, an excess of nonsynonymous compared to synonymous polymorphisms indicated historical balancing selection. We also studied association between helminth (Trichinella spp. and Echinococcus canadensis) prevalence and MHC diversity at allele and SNP level. MHC‐heterozygous wolves were less often infected by Trichinella spp. and carriers of specific MHC alleles, SNP haplotypes and SNP alleles had less helminth infections. The associated SNP haplotypes and alleles were shared by different MHC alleles, which emphasizes the necessity of single‐nucleotide‐level association studies also in MHC. Here, we show that strong balancing selection has had similar effect on MHC diversities in the Finnish and Russian Karelian wolf populations despite significant genetic differentiation at neutral markers and small population size in the Finnish population.     

Polymorphism at the ovine major histocompatibility complex class II loci

Southern hybridization analysis of the ovine major histocompatibility complex (MHC) ( MhcOvar ) class II region, using sheep-specific probes for the DQA1, DQA2, DQB and DRA loci, has revealed extensive polymorphism. DQA1 and DQAP had eight and 16 alleles respectively, DQB had six and DRA had three alleles. Little information was derived from the DRB locus owing to extensive cross-hybridization between the DRB probe and the DQB locus. Differences in allele frequency between breeds were revealed. At the DQA1 locus a null allele (DQA1-N) was observed with a frequency of between 27% and 45%, making this the most common DQA1 allele in all breeds examined. The frequency of DQA1-N homozygotes was between 11% and 18%, raising questions as to the functional significance of the DQA1 gene. Linkage analysis between the DQA1, DQA2, DQB and DRA loci did not reveal any recombination.     

Sequence diversity of MHC class-II <Emphasis Type="BoldItalic">DRB</Emphasis> gene in gazelles (<Emphasis Type="BoldItalic">Gazella subgutturosa</Emphasis>) raised in Sanliurfa of Turkey

In this study, we aimed to assess the sequence diversity of major histocompatibility complex (MHC) class-II DRB gene at exon 2 in gazelles raised in Sanliurfa Province of Turkey. Twenty DNA samples isolated from gazelles (Gazella subgutturosa) were used for sequencing exon 2 of MHC class-II DRB gene. Target region was amplified by polymerase chain reaction (PCR) and their products were directly sequenced. Nine of these 20 samples yielded unambiguously readable sequences. Three of the nine samples were homozygotes and each showed different sequences. A 262-bp sequence obtained from the three homozygote samples were submitted to GenBank (accession numbers: KC309405, KC309406 and KC309407). Using an allele specific PCR, we detected 10 additional haplotypes. Among 13 haplotypes, 45 nucleotide positions were polymorphic and most of the polymorphic nucleotide positions localized at peptide-binding region (PBR). Rates of nonsynonymous substitutions were significantly higher than synonymous substitutions at PBR. Phylogenetic analysis of the haplotypes showed that 10 haplotypes of the gazelles were clustered together while three were clustered with ovine and bovine haplotypes. The results indicated that at least 13 haplotypes at exon 2 of MHC class-II DRB gene were showing high degree of nucleotide and amino acid diversity, and certain haplotypes of G. subgutturosa were more similar to haplotypes from sheep or cattle than to each other. Rates of synonymous and nonsynonymous substitutions suggested that positive selection was a driving force for diversity at this locus in G. subgutturosa.