An extended HLA-DQ-DR haplotype rather than DRB1 alone contributes to RA predisposition

 In the present study, we tested our hypothesis on the role of a DQ-DR haplotype in rheumatoid arthritis (RA) predisposition. Using two groups of patients and controls, one from The Netherlands and one from Switzerland, we found that DQA1*0301-homozygous and DQA1*0301//DQA1*0101/04-heterozygous individuals are highly predisposed to RA in both populations, while DQA1*0101/04-homozygous are not. The DQA1*0301-DRB1*0403/06/07 and DQA1*0301-DRB1*0901 haplotypes are not associated with RA by themselves but strongly increase the risk of developing disease in DQA1*0301- and DQA1*0101/04-heterozygous. DRB1 alleles carrying the motif DERAA in their third hypervariable region, i.e., *0103, *0402, *1102, *1103, *1301, and *1302, provide a long-lasting protection against RA in DQA1*0101/04- but not in DQA1*0301-positive individuals. These data show that considering both DQ and DR gives a better distinction between patients and controls than the shared epitope hypothesis. Received: 5 March 1998 / Revised: 21 April     

HLA class II DR and DQ genotypes and haplotypes associated with rheumatic fever among a clinically homogeneous patient population of Latvian children

The HLA system is being paid more and more attention because it is very significant in polymorphous immunological reactions. Several studies have suggested that genetic susceptibility to rheumatic fever (RF) and rheumatic heart disease (RHD) is linked to HLA class II alleles. We hypothesized that HLA class II associations within RHD may be more consistent if analysed amongst patients with a relatively homogeneous clinical outcome. A total of 70 RF patients under the age of 18 years were surveyed and analysed in Latvia. HLA genotyping of DQA1, DQB1 and DRB1 was performed using PCR with amplification with sequence-specific primers. We also used results from a previous study of DQB1 and DRB1 genotyping. In the RF patients, HLA class II DQA1*0401 was found more frequently compared to DQA1*0102. In the RF homogeneous patient groups, DQA1*0402 has the highest odds ratio. This is also the case in the multivalvular lesion (MVL) group, together with DQA1*0501 and DQA1*0301. In the chorea minor patients, DQA1*0201 was often found. Significant HLA DQA1 protective genotypes were not detected, although DQA1 genotypes *0103/*0201 and *0301/*0501 were found significantly and frequently. In the distribution of HLA DRB1/DQA1 genotypes, *07/*0201 and *01/*0501 were frequently detected; these also occurred significantly often in the MVL group. The genotype *07/*0201 was frequently found in Sydenhamn's chorea patients that had also acquired RHD, but DRB1*04/DQA1*0401 was often apparent in RF patients without RHD. In the distribution of HLA DQA1/DQB1 genotypes, both in RF patients and in the homogeneous patient groups, the least frequent were *0102/*0602-8. The genotype DQA1*0501 with the DQB1 risk allele *0301 was often found in the MVL group. The genotype *0301/*0401-2 was frequently found in the RF and Sydenhamn's chorea patient groups. The haplotype *07-*0201-*0302 was frequently found in RF and homogeneous patient groups, including the MVL group. In addition, haplotypes *04-*0401-*0301 and *04-*0301-*0401-2 were frequent amongst patients with Sydenhamn's chorea. The protective alleles DQA1*0102 and DQB1*0602-8 in the haplotype DRB1*15 were less frequently found in RF patients. The results of the present study support our hypothesis and indicate that certain HLA class II haplotypes are associated with risk for or protection against RHD and that these associations are more evident in patients in clinically homogeneous groups.     

MHC binding peptides for designing of vaccines against Japanese encephalitis virus: A computational approach

Japanese encephalitis (JE), a viral disease has seen a drastic and fatal enlargement in the northern states of India in the current decade. The better and exact cure for the disease is still in waiting. For the cause an in silico strategy in the development of the peptide vaccine has been taken here for the study. A computational approach to find out the Major Histocompatibility Complex (MHC) binding peptide has been implemented. The prediction analysis identified MHC class I (using propred I) and MHC class II (using propred) binding peptides at an expectable percent predicted IC (50) threshold values. These predicted Human leukocyte antigen [HLA] allele binding peptides were further analyzed for potential conserved region using an Immune Epitope Database and Analysis Resource (IEDB). This analysis shows that HLA-DRB1*0101, HLA-DRB3*0101, HLA-DRB1*0401, HLA-DRB1*0102 and HLA-DRB1*07:01% of class II (in genotype 2) and HLA-A*0101, HLA-A*02, HLA-A*0301, HLA-A*2402, HLA-B*0702 and HLA-B*4402% of HLA I (in genotype 3) bound peptides are conserved. The predicted peptides MHC class I are ILDSNGDIIGLY, FVMDEAHFTDPA, KTRKILPQIIK, RLMSPNRVPNYNLF, APTRVVAAEMAEAL, YENVFHTLW and MHC class II molecule are TTGVYRIMARGILGT, NYNLFVMDEAHFTDP, AAAIFMTATPPGTTD, GDTTTGVYRIMARGI and FGEVGAVSL found to be top ranking with potential super antigenic property by binding to all HLA. Out of these the predicted peptide FVMDEAHFTDPA for allele HLA-A*02:01 in MHC class I and NYNLFVMDEAHFTDP for allele HLA-DRB3*01:01 in MHC class II was observed to be most potent and can be further proposed as a significant vaccine in the process. The reported results revealed that the immune-informatics techniques implemented in the development of small size peptide is useful in the development of vaccines against the Japanese encephalitis virus (JEV).     

Genetic History of Aleuts of the Commander Islands as Revealed by the Analysis of the HLA Class II Gene Variability

Variability of the HLA class II genes (alleles of the DRB1, DQA1, and DQB1 loci) was investigated in a sample of Aleuts of the Commanders (n = 31), whose ancestors inhabited the Commander Islands for many thousand years. Among 19 haplotypes revealed in the Aleuts of the Commanders, at most eight were inherited from the native inhabitants of the Commander Islands. Five of these haplotypes (DRB1*0401-DQA1*0301-DQB1*0301, DRB1*1401-DQA1*0101-DQB1*0503, DRB1*0802-DQA1*0401-DQB1*0402, DRB1*1101-DQA1*0501-DQB1*0301, and DRB1*1201-DQA1*0501-DQB1*0301) were typical of Beringian Mongoloids, i.e., Coastal Chukchi and Koryaks, as well as Siberian and Alaskan Eskimos. Genetic contribution of the immigrants to the genetic pool of the proper Aleuts constituted about 52%. Phylogenetic analysis based on Transberingian distribution of the DRB1 allele frequencies favored the hypothesis on the common origin of the Paleo-Aleuts, Paleo-Eskimos, and the Indians from the northwestern North America, whose direct ancestors survived in Beringian/southwestern Alaskan coastal refugia during the late Ice Age.     

HLA DOA1 and DOB1 loci in Honduran women with cervical dysplasia and invasive cervical carcinoma and their relationship to human papillomavirus infection.

Molecular and epidemiological studies have demonstrated that certain types of human papillomavirus (HPV), mainly HPV-16 and HPV-18, are the primary causes of cervical cancer and its precursor lesions; there is now evidence for a clear association with specific HLA class I and class II loci contributing independently to the expression of cervical cancer. Among Honduran women carcinoma of the cervix is the most common type of cancer, and infections with high-risk HPV types are highly prevalent. To study the interactive role of viral-host genetics, we performed PCR amplification of DNA and sequence-specific oligonucleotide probe typing on cervical scrapes from 49 women [24 with cervical intraepithelial neoplasia stage III or cervical cancer (severe cases) and 25 with stage I or II cervical intraepithelial neoplasia (mild cases)] and 75 control subjects to look for possible associations between HPV and HLA class II DQA1 and DQB1 alleles in the development of dysplasias and invasive cancer. This analysis revealed a predominance of HLA-DQA1*0301 among severe-case patients [relative risk (RR) = 3.45, p = 0.008), whereas DQA1*0501 was negatively associated (RR = 0.30, p = 0.03), suggesting a protective effect of this allele. HPV typing showed a decreased relative risk among the HPV-16 or HPV-18 carrying patients and other HPV-related positive patients in the presence of DQB1*0602 compared with positive control subjects (p = 0.04). No statistically significant allele frequency difference was observed between mild dysplasia cases and control subjects. The results suggest that DQA1*03011, which is in linkage desequilibrium with all HLA-DR4 alleles, confers an increased risk for severe cervical dysplasia and invasive cancer, whereas DQA1*0501, which is in several DR52 haplotypes, has a protective effect. Furthermore, specific HLA-DQB1 sequences may be important in determining the immune response to HPV peptides and may affect the risk for cervical cancer after HPV infection in mestizo Honduran women.     

HLA class II gene and haplotype diversity in the population of the island of Hvar, Croatia

The DRB1, DRB3, DRB5, DQA1 and DQB1 allele polymorphisms were analysed in 3 western and 3 eastern villages of the island of Hvar using PCR-SSOP method and 12th International Workshop primers and probes. Three DQB1 alleles (*0304, *0305, *0607) detected in the population of the island of Hvar (HP) have not yet been observed in general Croatian population (GCP). Significant differences were observed between two regions of Hvar for: a) DRB1*0701 allele (p < 0.001), b) DQA1*0201 allele (p < 0.01), and c) DRB1*0101-DQA1*0101-DQB1*0501 haplotypic association (p < 0.05). Two unusual haplotypic associations, which have not yet been described in general Croatian population (GCP), DRB1*0101-DQA1*0102-DQB1*0501 and DRB1*1501-DQA1 *0102-DQB1*0604 were observed in the population from the island of Hvar (HP). Measures of genetic kinship and genetic distances revealed isolation and clusterization which coincides with the known ethnohistorical, as well as biological and biocultural data obtained from a series of previous investigations. The five studied village subpopulations formed two clusters (East-West) to which the far eastern village (with the highest rii of 0.0407) joined later, thus indicating possible impact of historical immigrations from the mainland.     

Differences in DNA sequence specificity among MHC class II X box binding proteins

The class II (Ia) MHC Ag are integral membrane proteins whose expression is limited to specific cell types. A pair of consensus sequences, X and Y, is found upstream from all class II genes and deletion of each of these sequences eliminates expression of transfected genes. Cells that express Ia demonstrate a coordinate response to lymphokines and other stimuli. These conserved sequences might, therefore, play a role in tissue specificity or lymphokine inducibility of Ia gene expression. The X box sequence of the murine class II A alpha gene diverges much more substantially from the X consensus than does the Y box motif of this gene. We demonstrate that this X box motif is nonetheless recognized by sequence-specific DNA-binding proteins, as is the more closely conserved Y box. Gel retardation assays and DNase I footprints were compared for a panel of Ia+ and Ia- cells as well as for cells stimulated with the Ia-inducing lymphokines IL-4 and IFN-gamma. The level, retardation pattern and region of DNA contact were comparable in all instances. Thus the availability of active DNA-binding X and Y box factors cannot alone account for the regulation of A alpha expression. To test whether the same set of proteins binds all class II MHC conserved motifs, oligonucleotide probe binding and cross-competition experiments with X box sequences from A alpha, E alpha, and E beta genes were performed. These studies demonstrated A alpha, E alpha, and E beta DNA-protein complexes with unique mobilities and specificities. In addition, all three X box oligonucleotide probes generated one faint complex with an affinity profile of E beta greater than E alpha much greater than A alpha. These three complexes comigrated and thus may represent a communal binding protein. The data are most consistent with the conclusion that multiple proteins bind class II MHC X boxes. For A alpha, the predominant complexes represent different specificities from the predominant E alpha and E beta X box binding proteins.     

Isolation and characterization of the major histocompatibility complex <Emphasis Type="BoldItalic">DQA1</Emphasis> and <Emphasis Type="BoldItalic">DQA2</Emphasis> genes in gayal (<Emphasis Type="BoldItalic">Bos frontalis</Emphasis>)

The species origin of Yunnan gayal has been controversial since many years. However, few recent genetic studies have suggested that it has perhaps originated from the hybridization between male Bos frontalis and female B. taurus or B. indicus. Being an important semi-wild bovid species, this has also been listed under the red list of International Union of Conservation of Nature and Natural Resources. However, there is limited information available about the immunogenicity of this precarious species of Bos. Major histocompatibility complex (MHC) plays a pivotal role in immune response to infectious diseases in vertebrates. In the present study, we have investigated the structural and functional characteristics and possible duplication of the MHC-DQA genes in gayal (B. frontalis). Two full-length cDNA clones of the MHC-DQA genes were amplified and designated as Bofr-DQA1 (DQA*0101) and Bofr-DQA2 (DQA*2001) with GenBank accession numbers KT318732 and KT318733, respectively. A comparison between Bofr-DQA1, Bofr-DQA2 and to other MHC-DQA molecules from different animal species showed that nucleotide and encoded amino acid sequences of these two identified MHC-DQA genes have more similarity to alleles of specific DQA1 and DQA2 molecules from other Ruminantia species than to each other. The phylogenic investigation also demonstrated a large genetic distance between these two genes than to homologous from the other species. The large genetic distance between Bofr-DQA1 and Bofr-DQA2, and the presence of different bovine DQA putative motifs clarify that these sequences are nonallelic type. These results could suggest that duplication of the DQA genes has also occurred in gayal. The findings of the present study have strengthened our understanding to MHC diversity in rare ruminants and mutation of immunological functions, selective and evolutionary forces that affect MHC variation within and between species.     

HLA-DPB1 alleles in a population from South China

Using the polymerase chain reaction (PCR) and hybridization with oligonucleotide probes, we analyzed the distribution of DPB1 alleles in 99 healthy unrelated individuals from the city of Guangzhou (Canton), South China. Twelve different DPB1 alleles were found in this panel. The most common allele was DPB1*0501 (62.6%). Other major alleles detected included DPB1*02 (DPB1*0201 and DPB1*0202), DPB1*1301, DPB1*0401, and a recently described allele, designated DPB1*2101. The hybridization pattern of DPB1*2101 showed that this allele shared sequences with DPB1*0301 and DPB1*0601 in the A and F hypervariable regions, while the C, D, and E regions were identical to those of DPB1*0202. DPB1*2101 was observed in 11% of the subjects tested. It was found to be in strong linkage dis-equilibrium with DRB1*1202. In family studies, segregation of the haplotype DRB1*1202, DRB3*0301, DQA1*0601, DQB1*0301, DPB1*2101 was observed. The second exon of DPB1*2101 was sequenced from codon 8 to codon 90 and the sequence, inferred from the pattern of hybridization, was confirmed. DPB1*0301, DPB1*0402, DPB1*0101, DPB1*1401, DPB1*1901, and another recently recognized allele, now designated DPB1*2401, were detected with low frequencies. DPB1*2401 had the same hybridization pattern as DPB1*0501 except for a probe that matches codons 85–90. In this region, DPB1*2401 encoded the amino acid sequence GPMTLQ instead of EAVTLQ as in DPB1*0501.     

A distinct HLA-DRw8 haplotype characterizes patients with juvenile rheumatoid arthritis

We studied the first domain of the HLA-DRB1, HLA-DQA1, and HLA-DQB1 loci of 67 HLA-DRw8-positive Caucasians including 43 with early-onset pauciarticular juvenile rheumatoid arthritis (EOPA-JRA, alternatively known as early-onset pauciarticular juvenile chronic arthritis). Serology, restriction fragment length polymorphism (RFLP), and polymerase chain reaction (PCR) oligotyping revealed that 62, including all the EOPA-JRA patients, carried the HLA-DRB1*0801, DQA1*0401, DQB1*0402 genotype. Approximately onefifth of the controls carried atypical HLA-DRB1, HLA-DQA1, and/or HLA-DQB1 loci on their HLA-DRw8 haplotype confirmed by family studies. DNA sequences of HLA-DRB1, DQA1, and DQB1 alleles in patients and controls were identical to those previously reported. Disease association studies in 113 EOPA-JRA patients and 207 controls unselected for HLA-DRw8 revealed that the HLA-DRB1*0801, DQA1*0401, DQB1*0402 genotype was associated with a higher relative risk (RR) for disease (RR = 12.8, ² = 48.8, P < 10^–4) than was the serologically defined presence of HLA-DRw8 (RR = 8, ² = 39, P < 10^–4). Further analysis suggested that the DQ genes on HLA-DRw8 haplotypes are as likely as the DR genes to contribute to the pathogenesis of EOPA-JRA. This study increases to five the number of HLA-DR/DQ haplotypes identified in HLA-DRw8 Caucasians.The nucleotide sequence data reported in this paper have been submitted to the GenBank nucleotide sequence database and have been assigned the accession number M34308.     

Peptide motif for the rat MHC class II molecule RT1.Da: similarities to the multiple sclerosis-associated HLA-DRB1*1501 molecule

Experimental autoimmune encephalomyelitis induced with myelin proteins in DA and LEW.1AV1 rats is a model of multiple sclerosis (MS). It reproduces major aspects of this detrimental disease of the central nervous system. MS is associated with the HLA-DRB1*1501, DRB5*0101, and DQB1*0602 haplotype. DA and LEW.1AV1 rats share the RT1^av1 haplotype. So far, no MHC class II peptide motif of RT1.D^a molecules has been described. Sequence alignment of the chain of the rat MHC class II molecule RT1.D^a with human HLA class II molecules revealed strong similarity in the peptide-binding groove of RT1.D^a and HLA-DRB1*1501. According to the putative peptide-binding pockets of RT1.D^a, after comparison with the pockets of HLA-DRB1*1501, we predicted the peptide motif of RT1.D^a. To verify the predicted motif, naturally processed peptides were eluted by acidic treatment from immunoaffinity-purified RT1.D^a molecules of lymphoid tissue of DA rats and subsequently analyzed by ESI tandem mass spectrometry. In addition, we performed binding studies with combinatorial nonapeptide libraries to purified RT1.D^a molecules. Based on these studies we could define a peptide-binding motif for RT1.D^a characterized by aliphatic amino acid residues (L, I, V, M) and of F for the peptide pocket P1, aromatic residues (F, Y, W) for P4, basic residues (K, R) for P6, aliphatic residues (I, L, V) for P7, and aromatic residues (F, Y, W) and L for P9. Both methods revealed similar binding characteristics for peptides to RT1.D^a. This data will allow epitope predictions for analysis of peptides, relevant for experimental autoimmune diseases.     

Polymorphism in a T-cell receptor variable gene is associated with susceptibility to a juvenile rheumatoid arthritis subset

This report demonstrates a T-cell receptor (Tcr) restriction fragment length polymorphism, defined by a Tcrb-V6.1 gene probe and Bgl II restriction enzyme, to be absolutely correlated with allelic variation in the coding sequence of a Tcrb-V6.1 gene. A pair of non-conservative amino acid substitutions distinguish the Tcrb-V6.1 allelic variants. An association of this Tcrb-V6.1 gene allelic variant with one form of juvenile rheumatoid arthritis (JRA) was established in a cohort of 126 patients. The association was observed in patients possessing the HLA-DQA1*0101 gene. Among HLA-DQA*0101 individuals, 19 of 26 patients (73.1%) carried one particular Tcrb-V6.1 gene allele as opposed to 11 of 33 controls (33%; p<0.005). Haplotypes carrying this HLA gene have previously been shown to confer increased risk for progression of arthritis in JRA. This demonstration of a disease-associated Tcrb-V gene allelic variant has not, to our knowledge, been previously reported and supports the contribution of polymorphism in the Tcr variable region genomic repertoire to human autoimmune disease.The nucleotide sequence data reported in this paper have been submitted to the GenBank nucleotide sequence database and have been assigned the accession numbers M67511 for V6.1A and M67512 for V6.1B.     

Genetic history of Aleuts of the Komandor islands from results of analyzing variability of class II HLA genes

Variability of the HLA class II genes (alleles of the DRB1, DQA1, and DQB1 loci) was investigated in a sample of Aleuts of the Commanders (n = 31), whose ancestors inhabited the Commander Islands for many thousand years. Among 19 haplotypes revealed in Aleuts of the Commanders, at most eight were inherited from the native inhabitants of the Commander Islands. Five of these haplotypes (DRB1*0401-DQA1*0301-DQB1*0301, DRB1*1401-DQA1*0101-DQB1*0503, DRB1*0802-DQA1*0401-DQB1*0402, DRB1*1101-DQA1*0501-DQB1*0301, and DRB1*1201-DQA1*0501-DQB1*0301) were typical of Beringian Mongoloids, i.e., Coastal Chukchi and Koryaks, as well as Siberian and Alaskan Eskimos. Genetic contribution of the immigrants to the genetic pool of proper Aleuts constituted about 52%. Phylogenetic analysis based on Transberingian distribution of the DRB1 allele frequencies favored the hypothesis on the common origin of Paleo-Aleuts, Paleo-Eskimos, and the Indians from the northwestern North America, whose direct ancestors survived in Beringian/southwestern Alaskan coastal refugia during the late Ice Age.     

The HLA molecules DQA1*0501/B1*0201 and DQA1*0301/B1*0302 share an extensive overlap in peptide binding specificity

Assays to measure the binding capacity of peptides for HLA-DQA1*0501/B*0201 (DQ2.3) and DQA1*0301/B*0302 (DQ3.2) were developed using solubilized MHC molecules purified from EBV-transformed cell lines. These quantitative assays, based on the principle of the inhibition of binding of a high-affinity radiolabeled ligand, were validated by examining the binding capacity of known DQ-restricted epitopes or ligands. The availability of these assays allowed an investigation of patterns of cross-reactivity between different DQ molecules and with various common DR molecules. DQ2.3 and DQ3.2 were found to have significantly overlapping peptide binding repertoires. Specifically, of 13 peptides that bound either DQ2.3 or DQ3.2, nine (69.2%) bound both. The molecular basis of this high degree of cross-reactivity was further investigated with panels of single substitution analogs of the thyroid peroxidase 632-645Y epitope. It was found that DQ2.3 and DQ3.2 bind the same ligands by using similar anchor residues but different registers. These data suggest that in analogy to what was previously described for HLA-DR molecules, HLA-DQ supertypes characterized by largely overlapping binding repertoires can be defined. In light of the known linkage of both HLA-DQ2.3 and -DQ3.2 with insulin-dependent diabetes mellitus and celiac disease, these results might have important implications for understanding HLA class II autoimmune disease associations.     

How specific MHC class I and class II combinations affect disease resistance against infectious salmon anaemia in Atlantic salmon (Salmo salar)

The aim was to evaluate the performance of selected individual MHC class I and class II alpha (A) alleles, and combinations of these on disease resistance against infectious salmon anaemia (ISA). The material consisting of 1966 fish from seven families, contained five MHC class I alleles and four MHC class II A alleles. Which representing given class II A and class II beta (B) haplotypes, totalling 19 MHC class I and class II A genotypes. The fish were challenged with infectious salmon anaemia virus (ISAV), the virus causing ISA. Dead fish were collected daily during the challenge experiment and the survivors were collected at termination. All fish were genotyped for MHC class I and class II A. The total mortality in the material was 85.14%. For MHC class I, UBA*0201 and UBA*0301 were significantly the most resistant alleles, while UBA*0601 for class I and DAA*0301 for class II A were the significantly most susceptible alleles. The analysis of combined MHC class I and class II A genotypes detected that fish with the genotype UBA*0201/*0301;DAA*0201/*0201 were the most resistant fish with a hazard ratio (HR) at 0.750, while the fish with the genotypes UBA*0601/*0801;DAA*0501/*0501 and UBA*0201/*0301;DAA*0301/*0501 were the most susceptible fish with HR of 1.334 and 1.425. In addition, Cox regression analysis within family detected combined MHC class I and class II A genotypes that contributed significantly to resistance and susceptibility. The study confirmed the expectation of performance of individual MHC class I and class II A alleles, and also detected an effect of MHC class I and class II A in combinations.     

Structural and functional studies of trans-encoded HLA-DQ2.3 (DQA1*03:01/DQB1*02:01) protein molecule

MHC class II molecules are composed of one α-chain and one β-chain whose membrane distal interface forms the peptide binding groove. Most of the existing knowledge on MHC class II molecules comes from the cis-encoded variants where the α- and β-chain are encoded on the same chromosome. However, trans-encoded class II MHC molecules, where the α- and β-chain are encoded on opposite chromosomes, can also be expressed. We have studied the trans-encoded class II HLA molecule DQ2.3 (DQA1*03:01/DQB1*02:01) that has received particular attention as it may explain the increased risk of certain individuals to type 1 diabetes. We report the x-ray crystal structure of this HLA molecule complexed with a gluten epitope at 3.05 Å resolution. The gluten epitope, which is the only known HLA-DQ2.3-restricted epitope, is preferentially recognized in the context of the DQ2.3 molecule by T-cell clones of a DQ8/DQ2.5 heterozygous celiac disease patient. This preferential recognition can be explained by improved HLA binding as the epitope combines the peptide-binding motif of DQ2.5 (negative charge at P4) and DQ8 (negative charge at P1). The analysis of the structure of DQ2.3 together with all other available DQ crystal structures and sequences led us to categorize DQA1 and DQB1 genes into two groups where any α-chain and β-chain belonging to the same group are expected to form a stable heterodimer.     

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.