MHC class II genes in European wolves: a comparison with dogs

The genome of the grey wolf, one of the most widely distributed land mammal species, has been subjected to both stochastic factors, including biogeographical subdivision and population fragmentation, and strong selection during the domestication of the dog. To explore the effects of drift and selection on the partitioning of MHC variation in the diversification of species, we present nine DQA, 10 DQB, and 17 DRB1 sequences of the second exon for European wolves and compare them with sequences of North American wolves and dogs. The relatively large number of class II alleles present in both European and North American wolves attests to their large historical population sizes, yet there are few alleles shared between these regions at DQB and DRB1. Similarly, the dog has an extensive array of class II MHC alleles, a consequence of a genetically diverse origin, but allelic overlap with wolves only at DQA. Although we might expect a progression from shared alleles to shared allelic lineages during differentiation, the partitioning of diversity between wolves and dogs at DQB and DRB1 differs from that at DQA. Furthermore, an extensive region of nucleotide sequence shared between DRB1 and DQB alleles and a shared motif suggests intergenic recombination may have contributed to MHC diversity in the Canidae.     

Allelic Combinations of Promoter and Exon 2 in <Emphasis Type="Italic">DQB1</Emphasis> in Dogs and Wolves

Polymorphism of PBRs of the major histocompatibility complex (MHC) genes is well recognized, but the polymorphism also extends to proximal promoter regions. Examining DQB1 variability in dogs and wolves, we identified 7 promoter variants and 13 exon 2 alleles among 89 dogs, including a previously unknown DQB1 exon 2 allele, and 8 promoter variants and 9 exon 2 alleles among 85 wolves. As expected from previous studies and from a close chromosomal location, strong linkage disequilibrium was demonstrated in both wolves and dogs by having significantly fewer promoter/exon 2 combinations than expected from simulations of randomized data sets. Interestingly, we noticed weaker haplotypic associations in dogs than in wolves. Dogs had twice as many promoter/exon 2 combinations as wolves and an almost 2-fold difference in the number of exon 2 alleles per promoter variant. This difference was not caused by an admixture of breeds in our group of dogs because the high ratio of observed to expected number of haplotypes persisted within a single dog breed, the German Shepherd. Ewens-Watterson tests indicated that both the promoter and exon 2 are under the balancing selection, and both regions appear to be more recently derived in the dog than in the wolf. Hence, although reasons for the differences are unknown, they may relate to altered selection pressure on patterns of expression. Deviations from normal MHC expression patterns have been associated with autoimmune diseases, which occur frequently in several dog breeds. Further knowledge about these deviations may help us understand the source of such diseases.     

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.     

Major histocompatibility complex variation at three class II loci in the northern elephant seal

Northern elephant seals were hunted to near extinction in the 19th century, yet have recovered remarkably and now number around 175,000. We surveyed 110 seals for single-strand conformation polymorphism (SSCP) and sequence variation at three major histocompatibility (MHC) class II loci (DQA, DQB and DRB) to evaluate the genetic consequences of the population bottleneck at these loci vs. other well-studied genes. We found very few alleles at each MHC locus, significant variation among breeding sites for the DQA locus, and linkage disequilibrium between the DQB and DRB loci. Northern elephant seals are evidently inbred, although there is as yet no evidence of correlative reductions in fitness.     

Substantial functional diversity accompanies limited major histocompatibility complex class II variability in golden jackal (Canis aureus): A comparison between two wild Canis species in Croatia

The golden jackal (Canis aureus) is one of the less studied carnivores and research on its major histocompatibility complex (MHC) variability is just at its early stages. MHC genes encode cell-surface receptors that serve to bind and present antigens to T cells, which is essential to initiating specific immunological responses in vertebrates. In this paper we present for the first time patterns of genetic diversity and natural selection on MHC class II DLA-DRB1, DQA1 and DQB1 loci in the golden jackal using samples from two geographically distinct regions in Croatia and further compare them to the values found in its congener grey wolf (Canis lupus). Diversity of golden jackals at all three loci was markedly lower than that of grey wolves (allelic richness values were 4, 2 and 3 in jackal versus 11.9, 6.6 and 10.2 in wolves for DRB1, DQA1 and DQB1, respectively) and can be attributed to a genetic drift rather than to the lack of historical positive selection. The finding of high evolutionary distances (16.3% for DRB1 and 8.5% for DQB1) and a substantial number of codons predicted to be under the influence of positive selection (11 for DRB1 and 9 for DQB1) suggests that the investigated golden jackal population still contains considerable functional diversity necessary for the presentation of varied foreign peptides. In contrast to neutral genetic variation, our results suggest that the Dalmatian population has a higher MHC diversity than the Slavonian population, casting doubt on its supposed isolation and calling for a more extensive investigation of the MHC variability of southern Balkan jackal populations.     

MHC-DRB exon 2 allele polymorphism in Indian river buffalo (Bubalus bubalis)

The polymorphism of the major histocompatibility complex (MHC) class II DRB gene of riverine buffalo (Bubalus bubalis) was studied. Second exon sequences from the buffalo DRB locus, homologous to the cattle DRB3 gene, were amplified and characterized. A combination of single strand conformation polymorphism (SSCP) and heteroduplex analysis (HA) in a non-denaturing gel was used to identify new DRB second exon sequences. SSCP, HA and finally sequencing allowed the identification of 22 MHC-DRB exon 2 alleles from 25 unrelated Indian river buffalo. These are the first river buffalo DRB second exon sequences reported. A high degree of polymorphism in the sequences encoding the peptide binding regions was observed and some amino acid substitutions were found unique to the river buffalo.     

Regression mapping of association between the human leukocyte antigen region and Graves disease

The human leukocyte antigen class II genes DRB1, DQB1, and DQA1 are associated with Graves disease (GD), but, because of strong linkage disequilibrium within this region, the primary etiological variant(s) remains unknown. In the present study, 871 patients with GD and 621 control subjects were genotyped at the DRB1, DQB1, and DQA1 loci. All three loci were associated with GD (P=1.45 x 10(-12), P=3.20 x 10(-5), and P=9.26 x 10(-12), respectively). Stepwise logistic-regression analysis showed that the association could be explained by either DRB1 or DQA1 but not by DQB1. To extend previous results, the amino acid sequence of the exon 2-encoded peptide-binding domain of DRB1 was predicted for each subject, and, by use of logistic regression, each position was analyzed for association with GD. Of 102 amino acids, 70 were uninformative; of the remaining 32 amino acids, 13 were associated with GD (P values ranged from 2.20 x 10(-4) to 1.2 x 10(-12)). The strongest association was at position beta 74. This analysis is consistent with the possibility that position beta 74 of exon 2 of the DRB1 molecule may have a specific and central role in autoantigen presentation by DRB1 to T lymphocytes. However, we cannot yet exclude a primary role for DQA1 or for other polymorphisms that affect DRB1 function or expression.     

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.     

A novel HURRAH protocol reveals high numbers of monomorphic MHC class II loci and two asymmetric multi-locus haplotypes in the Père David's deer

The Père David's deer is a highly inbred, but recovered, species, making it interesting to consider their adaptive molecular evolution from an immunological perspective. Prior to this study, genomic sequencing was the only method for isolating all functional MHC genes within a certain species. Here, we report a novel protocol for isolating MHC class II loci from a species, and its use to investigate the adaptive evolution of this endangered deer at the level of multi-locus haplotypes. This protocol was designated "HURRAH" based on its various steps and used to estimate the total number of MHC class II loci. We confirmed the validity of this novel protocol in the giant panda and then used it to examine the Père David's deer. Our results revealed that the Père David's deer possesses nine MHC class II loci and therefore has more functional MHC class II loci than the eight genome-sequenced mammals for which full MHC data are currently available. This could potentially account at least in part for the strong survival ability of this species in the face of severe bottlenecking. The results from the HURRAH protocol also revealed that: (1) All of the identified MHC class II loci were monomorphic at their antigen-binding regions, although DRA was dimorphic at its cytoplasmic tail; and (2) these genes constituted two asymmetric functional MHC class II multi-locus haplotypes: DRA1*01 ～ DRB1 ～ DRB3 ～ DQA1 ～ DQB2 (H1) and DRA1*02 ～ DRB2 ～ DRB4 ～ DQA2 ～ DQB1 (H2). The latter finding indicates that the current members of the deer species have lost the powerful ancestral MHC class II haplotypes of nine or more loci, and have instead fixed two relatively weak haplotypes containing five genes. As a result, the Père David's deer are currently at risk for increased susceptibility to infectious pathogens.     

High resolution molecular typing of HLA class II region in the population of the island of Krk, Croatia

The HLA class II alleles (DRB1, DRB3, DRB5, DQA1, and DQB1) and haplotypic associations were studied in the population of the island of Krk using the PCR-SSOP method and the 12th International Histocompatibility Workshop primers and probes. Allele and haplotypic frequencies were compared with the general Croatian population. Significant differences were observed between the population of the island of Krk and Croatians for: a) three broad specificities at DRB1 locus (DRB1*01, *15, and *07), b) one allele at DRB3 locus (DRB3*0301), c) one allele at DQA1 locus (DQA1*0201), d) one allele at DQB1 locus (DQB1*0303). Four unusual haplotypic associations, which have not yet been described in the Croatian population, DRB1*1301-DQA1*0103-DQB1*0607, DRB1*1302-DQA1*0102-DQB1*0605, DRB1*1305-DQA1*0102-DQB1*0605 and DRB1*1305-DQA1*0103-DQB1*0603 were observed in the population from the island of Krk.     

Mutations in the<Emphasis Type="Italic"> HLA</Emphasis> class II genes leading to loss of expression of HLA-DR and HLA-DQ in diffuse large B-cell lymphoma

Loss of expression of human leukocyte antigen (HLA) class II molecules on tumor cells affects the onset and modulation of the immune response through lack of activation of CD4+ T lymphocytes. Previously, we showed that the frequent loss of expression of HLA class II in diffuse large B-cell lymphoma (DLBCL) of the testis and the central nervous system (CNS) is mainly due to homozygous deletions in the HLA region on chromosome band 6p21.3. A minority of cases showed hemizygous deletions or mitotic recombination, implying that mutation of the remaining copy of the class II genes might be involved. Here, we studied three DLBCLs with loss of HLA-DQ expression for mutations in the DQB1 and DQA1 genes and three tumors with loss of HLA-DR expression for mutations in the DRB1 and DRA genes. In one case, a point mutation in exon 2 of the DQB1 gene, leading to the formation of a stop codon, was detected at position 47. In a second case, a stop codon was found at position 11 due to a deletion of 19 bp in exon 1 of the DRA gene. No mutations were found in the promoter sequences of the DRA, DQA1 and DQB1 genes. We conclude that both homozygous deletions and hemizygous deletions or mitotic recombination with mutations of the remaining allele may lead to loss of expression of the HLA class II genes, which is comparable to the mechanisms affecting HLA class I expression in solid cancers.     

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.     

Organization and polymorphism of bovine major histocompatibility complex class II genes as revealed by genomic hybridizations with bovine probes

Previous studies on restriction fragment length polymorphism of bovine major histocompatibility complex class II genes have primarily been based on the use of human probes. In the present study bovine probes for DQA, DQB, DRB and DYA were used for RFLP analysis of cattle genomic DNA digested with PvuII and TaqI. There was an excellent agreement between the RFLP results obtained with homologous and heterologous probes. Although a few 'new' restriction fragments were revealed with the bovine probes there was no discrepancy with regard to the classification of allelic types with the two types of probes. The major advantages of using bovine probes were a better hybridization signal and reduced cross-hybridization between loci. Hybridization experiments with DQA probes for the first domain exon from two different genomic clones revealed the presence of two distinct types of bovine DQA genes. Surprisingly, these probes did not cross-hybridize at high stringency, indicating that the two genes are quite divergent. Hybridization with a recently described genomic clone for a novel bovine alpha-chain gene confirmed that it corresponds to the DYA gene which had previously been identified by cross-hybridization to a human DQA probe.     

Cloning of the bovine major histocompatibility complex class II genes

Class II genes of the bovine major histocompatibility complex (MHC) have been cloned from a genomic library. The library was constructed in the bacteriophage lambda vector EMBL3 and comprises approximately 10 times the equivalent of the haploid genome. Half the library was screened with the human DQA, DQB, DRA and DRB cDNA probes. Of the 100 positively hybridizing phage clones, 37 were eventually fully characterized and mapped by means of Southern blot analysis. The exons encoding the first, second and transmembrane domain of all different A and B genes were subcloned and mapped in more detail. These analyses showed that these 37 clones were derived from five different A and 10 different B genes. The hybridization studies indicate that we have cloned and mapped two DQA genes, one DRA gene, two other A genes, four DQB genes, three DRB genes and three other B genes. Since the library was made from a heterozygous animal, this would suggest that there are at least one DQA, one DRA one other undefined A, two DQB, two DRB and one or two other undefined B genes in the haploid genome of Holstein Friesian cattle.