Human leukocyte antigen and interleukin 2, 10 and 12p40 cytokine responses to measles: is there evidence of the HLA effect?

HLA class I and class II associations were examined in relation to measles virus-specific cytokine responses in 339 healthy children who had received two doses of live attenuated measles vaccine. Multivariate linear regression modeling analysis revealed suggestions of associations between the expression of DPA1*0201 (p=0.03) and DPA1*0202 (p=0.09) alleles and interleukin-2 (IL-2) cytokine production (global p-value 0.06). Importantly, cytokine production and DQB1 allele associations (global p-value 0.04) revealed that the alleles with the strongest association with IL-10 secretion were DQB1*0302 (p=0.02), DQB1*0303 (p=0.07) and DQB1*0502 (p=0.06). Measles-specific IL-10 secretion associations approached significance with DRB1 and DQA1 loci (both global p-values 0.08). Specifically, suggestive associations were found between DRB1*0701 (p=0.07), DRB1*1103 (p=0.06), DRB1*1302 (p=0.08), DRB1*1303 (p=0.06), DQA1*0101 (p=0.08), and DQA1*0201 (p=0.04) alleles and measles-induced IL-10 secretion. Further, suggestive association was observed between specific DQA1*0505 (p=0.002) alleles and measles-specific IL-12p40 secretion (global p-value 0.09) indicating that cytokine responses to measles antigens are predominantly influenced by HLA class II genes. We found no associations between any of the alleles of HLA A, B, and Cw loci and cytokine secretion. These novel findings suggest that HLA class II genes may influence the level of cytokine production in the adaptive immune responses to measles vaccine.     

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.     

Genetic drift vs. natural selection in a long-term small isolated population: major histocompatibility complex class II variation in the Gulf of California endemic porpoise (Phocoena sinus)

Although many studies confirm long-term small isolated populations (e.g. island endemics) commonly sustain low neutral genetic variation as a result of genetic drift, it is less clear how selection on adaptive or detrimental genes interplay with random forces. We investigated sequence variation at two major histocompatibility complex (Mhc) class II loci on a porpoise endemic to the upper Gulf of California, México (Phocoena sinus, or vaquita). Its unique declining population is estimated around 500 individuals. Single-strand conformation polymorphism analysis revealed one putative functional allele fixed at the locus DQB (n = 25). At the DRB locus, we found two presumed functional alleles (n = 29), differing by a single nonsynonymous nucleotide substitution that could increase the stability at the dimer interface of alphabeta-heterodimers on heterozygous individuals. Identical trans-specific DQB1 and DRB1 alleles were identified between P. sinus and its closest relative, the Burmeister's porpoise (Phocoena spinipinnis). Comparison with studies on four island endemic mammals suggests fixation of one allele, due to genetic drift, commonly occurs at the DQA or DQB loci (effectively neutral). Similarly, deleterious alleles of small effect are also effectively neutral and can become fixed; a high frequency of anatomical malformations on vaquita gave empirical support to this prediction. In contrast, retention of low but functional polymorphism at the DRB locus was consistent with higher selection intensity. These observations indicated natural selection could maintain (and likely also purge) some crucial alleles even in the face of strong and prolonged genetic drift and inbreeding, suggesting long-term small populations should display low inbreeding depression. Low levels of Mhc variation warn about a high susceptibility to novel pathogens and diseases in vaquita.     

Polymorphism, recombination, and linkage disequilibrium within the HLA class II region.

Thirty-nine CEPH (Centre d'Etude du Polymorphisme Humain) families, comprised of 502 individuals, have been typed for the HLA class II genes DRB1, DQA1, DQB1, and DPB1 using nonradioactive sequence-specific oligonucleotide probes to analyze polymerase chain reaction amplified DNA. This population, which consists of 266 independent chromosomes, contains 27 DRB1, 7 DQA1, 12 DQB1, and 17 DPB1 alleles. Analysis of the distribution of allele frequencies using the homozygosity statistic, which gives an indication of past selection pressures, suggests that balancing selection has acted on the DRB1, DQA1, and DQB1 loci. The distribution of DPB1 alleles, however, suggests a different evolutionary past. Family data permits the estimation of recombination rates and the unambiguous assignment of haplotypes. No recombinants were found between DRB1, DQA1, and DQB1; however, recombinants were detected between DQB1 and DPB1, resulting in an estimated recombination fraction of greater than or equal to 0.008 +/- 0.004. Only 33 distinct DRB1-DQA1-DQB1 haplotypes were found in this population which illustrates the extreme nonrandom haplotypic association of alleles at these loci. A few of these haplotypes are unusual (previously unreported) for a Caucasian population and most likely result from past recombination events between the DR and DQ subregions. Examination of disequilibrium across the HLA region using these data and the available serologic HLA-A and HLA-B types of these samples shows that global disequilibrium between these loci declines with the recombination fraction, approaching statistic nonsignificance at the most distant interval, HLA-A to HLA-DP.DR-DQ haplotypes in linkage disequilibrium with DPB1 and B are noted and, finally, the evolutionary origin of certain class II haplotypes is addressed.     

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.     

Ovar-MHC Polymorphism in Malpura and Avikalin Sheep Vaccinated for Peste des Petits Ruminants (PPR) Virus

India harbors a vast diversity of sheep (40 breeds). The study was carried out to assess the genetic diversity of DRB1 and DQA2 locus of the ovar-MHC and their possible association with Peste des petits ruminants (PPR) virus vaccine response in Malpura and Avikalin sheep breeds maintained at an organized institute flock in the semi-arid region of India. Genetic analysis revealed the rich diversity of DRB1 locus with 23 alleles in Malpura and 21 alleles in Avikalin sheep that included 9 new alleles. DQA2 locus also had rich diversity with 19 alleles in Malpura and 20 alleles in Avikalin sheep that included 7 new alleles. At the protein level, high variability alike at the nucleotide level was observed. A marker for footrot susceptibility, DQA2*1101 was absent in both breeds. Genotypic association of DRB1 and DQA2 with PPR vaccine response was statistically non-significant. Vaccine response being a multifactorial (polygenic and influenced by environment) variable, could not show statistically significant association with MHC genotypes in the present study. However, rich genetic diversity of DRB1 and DQA2 gene reflects the importance of this locus for future selection programs.     

High similarity at three MHC loci between the baiji and finless porpoise: trans-species or convergent evolution?

Two DRA alleles and six MHC-I alleles were identified from a group of 15 baiji (Lipotes vexillifer), the most threatened cetacean in the world. Little sequence variation was detected at the DRA locus but extensive variation at the MHC-I locus. In combination with data at the DQB locus previously reported, three MHC loci exon 2 of the baiji all revealed striking similarity with those of the finless porpoise. Especially, some identical alleles shared by both species at the MHC-I and DQB loci suggested the convergent evolution as a consequence of common adaptive solutions to similar environmental pressures in the Yangtze River. As for DRA locus, the identity alleles were shared not only by baiji and finless porpoise but by some other cetacean species of the families Phocoenidae and Delphinidae, suggesting trans-species evolution on this gene.     

Nomenclature for factors of the BoLA system, 1996: report of the ISAG BoLA Nomenclature Committee

The BoLA (bovine lymphocyte antigen) Nomenclature Committee met during the 1994 and 1996 conferences of the International Society for Animal Genetics to define a sequence-based nomenclature system for genes of the BoLA system. The rules for acceptance of new sequences are described and names are assigned to the sequenced alleles of the class II genes DRA, DRB1, DRB2, DRB3, DQA, DQB, DYA, DIB, DMA and DMB . The assignment of BoLA class I sequences to loci will be considered at a later workshop when further sequencing/mapping data are available.     

Major Histocompatibility Complex (MHC) class II sequence polymorphism in long-finned pilot whale (Globicephala melas) from the North Atlantic

Determining how intra-specific genetic diversity is apportioned among natural populations is essential for detecting local adaptation and identifying populations with inherently low levels of extant diversity which may become a conservation concern. Sequence polymorphism at two adaptive loci (MHC DRA and DQB) was investigated in long-finned pilot whales (Globicephala melas) from four regions in the North Atlantic and compared with previous data from New Zealand (South Pacific). Three alleles were resolved at each locus, with trans-species allele sharing and higher levels of non-synonymous to synonymous substitution, especially in the DQB locus. Overall nucleotide diversities of 0.49?±?0.38% and 4.60?±?2.39% were identified for the DRA and DQB loci, respectively, which are relatively low for MHC loci in the North Atlantic, but comparable to levels previously described in New Zealand (South Pacific). There were significant differences in allele frequencies within the North Atlantic and between the North Atlantic and New Zealand. Patterns of diversity and divergence are consistent with the long-term effects of balancing selection operating on the MHC loci, potentially mediated through the effects of host-parasite coevolution. Differences in allele frequency may reflect variation in pathogen communities, coupled with the effects of differential drift and gene flow.     

HLA class II linkage disequilibrium and haplotype evolution in the Cayapa Indians of Ecuador.

DNA-based typing of the HLA class II loci in a sample of the Cayapa Indians of Ecuador reveals several lines of evidence that selection has operated to maintain and to diversify the existing level of polymorphism in the class II region. As has been noticed for other Native American groups, the overall level of polymorphism at the DRB1, DQA1, DQB1, and DPB1 loci is reduced relative to that found in other human populations. Nonetheless, the relative evenness in the distribution of allele frequencies at each of the four loci points to the role of balancing selection in the maintenance of the polymorphism. The DQA1 and DQB1 loci, in particular, have near-maximum departures from the neutrality model, which suggests that balancing selection has been especially strong in these cases. Several novel DQA1-DQB1 haplotypes and the discovery of a new DRB1 allele demonstrate an evolutionary tendency favoring the diversification of class II alleles and haplotypes. The recombination interval between the centromeric DPB1 locus and the other class II loci will, in the absence of other forces such as selection, reduce disequilibrium across this region. However, nearly all common alleles were found to be part of DR-DP haplotypes in strong disequilibrium, consistent with the recent action of selection acting on these haplotypes in the Cayapa.     

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.     

Expressed 2-5A synthetase genes and pseudogenes in the marine sponge Geodia barretti

In cattle, bovine leukocyte antigens (BoLAs) have been extensively used as markers for bovine diseases and immunological traits. In this study, we sequenced alleles of the BoLA class II loci, BoLA-DRB3 and BoLA-DQA1, from 650 Japanese cattle from six herds [three herds (507 animals) of Japanese Black cattle and three herds (143 animals) of Holstein cattle] using polymerase chain reaction-sequence-based typing (PCR-SBT) methods. We identified 26 previously reported distinct DRB3 alleles in the two populations: 22 in Japanese Black and 17 in Holstein. The number of DRB3 alleles detected in each herd ranged from 9 to 20. Next, we identified 15 previously reported distinct DQA1 alleles: 13 in Japanese Black and 10 in Holstein. The number of alleles in each herd ranged from 6 to 10. Thus, allelic divergence is significantly greater for DRB3 than for DQA1. A population tree on the basis of the frequencies of the DRB3 and DQA1 alleles showed that, although the genetic distance differed significantly between the two cattle breeds, it was closely related within the three herds of each breed. In addition, Wu-Kabat variability analysis indicated that the DRB3 gene was more polymorphic than the DQA1 gene in both breeds and in all herds, and that the majority of the hypervariable positions within both loci corresponded to pocket-forming residues. The DRB3 and DQA1 heterozygosity for both breeds within each herd were calculated based on the Hardy-Weinberg equilibrium. Only one Japanese Black herd showed a significant difference between the expected and observed heterozygosity at both loci. This is the first report presenting a detailed study of the allelic distribution of BoLA-DRB3 and -DQA1 genes in Japanese Black and Holstein cattle from different farms in Japan. These results may help to develop improved livestock breeding strategies in the future.     

Heterozygosity excess at the cattle DRB locus revealed by large scale genotyping of two closely linked microsatellites

A method for MHC DRB typing in cattle based on two closely linked and highly polymorphic microsatellites is described. The two microsatellites DRBP1ms and DRB3ms are located in intron 2 of the corresponding DRB gene. The very strong linkage disequilibrium between the two loci made it possible to establish DRB microsatellite haplotypes. The typing results with this method on reference samples followed closely that obtained with RFLP and direct sequence analysis of DRB3 exon 2. The method is well suited for large scale genotyping and was successfully applied for typing more than 600 unrelated animals representing 23 breeds. The data were used to test whether the observed DRB allele frequency distributions were consistent with that expected for selectively neutral alleles in populations at mutation-drift equilibrium. A significant heterozygosity excess was detected and there was an obvious trend across breeds towards a more even allele frequency distribution than expected. The deviation may be due to balancing selection acting on the DRB locus or by recent population bottlenecks.     

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.     

Single-strand conformational polymorphism and sequence polymorphism of Mhc-DRB in Latxa and Karrantzar sheep: implications for Caprinae phylogeny

Single-strand conformational polymorphism analysis and DNA sequencing were used to characterize Mhc-DRB second exon variability in the Latxa and Karrantzar breeds of sheep. The presence of more than two sequences in some animals indicates that alleles of two different loci have been amplified. Six new alleles were identified by sequencing. The allele frequency distribution of the DRB1 gene is striking, with two alleles accounting for half of the gene pool in both breeds under study. The most frequent allele in both breeds was the same (named DRB1*0702), with some specific amino acids: Tyr in position 31 and Thr in 51. A species variability analysis was also performed including the entire set of sheep DRB exon 2 sequences. Based on the patchwork patterns of different alleles, interallelic recombination appears to be playing a significant role in the generation of allelic diversity at this locus in sheep. The phylogenetic tree of all known Caprinae DRB sequences shows that certain alleles from one species are more closely related to those from other species than they are to each other. Allele DRB1*0702 merits special attention due to its high similarity to the Mufflon allele. As this is the most frequent in both breeds analyzed, one can hypothesize that in sheep, both Mufflon and Argali have had different influences depending on the sheep breed under study and that the relationship between domestic sheep and Mufflon is greater than previously thought. The data generated in this study can serve as a basis for developing a typing assay for the sheep DRB genes in the Latxa and Karrantzar populations.     

Polymorphisms in MHC-DRA and -DRB alleles of water buffalo (Bubalus bubalis) reveal different features from cattle DR alleles

Seventy-five individuals of Bubalus bubalis belonging to four different breeds, three of river buffalo and one of swamp buffalo, were studied for polymorphism in MHC DRB (Bubu-DRB) and DRA (Bubu-DRA) loci. Eight alleles of Bubu-DRB were found, and all alleles in the swamp type were shared with the three river breeds. All alleles sampled from the breed of European origin (Mediterranean) were present in breeds sampled in Brazil, thus variability of this locus may have been preserved to a great extent in the more recently founded Brazilian population. Bubu-DRB alleles contained higher proportions of synonymous vs. non-synonymous substitutions in the non-peptide-binding sites (PBS) region, in contrast to the pattern of variation found in BoLA-DRB3, the orthologous locus in cattle. This indicated that either the first domain exon (exon 2) of Bubu-DRB has not undergone as much recombination and/or gene conversion as in cattle alleles, or Bubu-DRB may be more ancient than BoLA-DRB3 alleles. Phylogenetic analysis of DRB alleles from Bubalus, Syncerus c. caffer, the Cape buffalo, and domestic cattle demonstrated transspecies polymorphism. Water buffalo contained two alleles of DRA that differed from each other in two amino acid positions, including one in the PBS (alpha22) that was also shared with Anoa depressicornis, the anoa. Discovery of variation in DRA was surprising as the first domain of DRA is a highly conserved polypeptide in mammals in general and especially in ruminants, where no other substitution in PBS was seen.     

Breadth of the CD4+ T cell response to Anaplasma marginale VirB9-1, VirB9-2 and VirB10 and MHC class II DR and DQ restriction elements

MHC class II molecules influence antigen-specific CD4+ T lymphocyte responses primed by immunization and infection. CD4+ T cell responses are important for controlling infection by many bacterial pathogens including Anaplasma marginale and are observed in cattle immunized with the protective A. marginale outer membrane (OM) vaccine. Immunogenic proteins that comprise the protective OM vaccine include type IV secretion system (T4SS) proteins VirB9-1, VirB9-2 and VirB10, candidates for inclusion in a multiepitope vaccine. Our goal was to determine the breadth of the VirB9-1, VirB9-2 and VirB10 T cell response and MHC class II restriction elements in six cattle with different MHC class II haplotypes defined by DRB3, DQA and DQB allele combinations for each animal. Overlapping peptides spanning each T4SS protein were tested in T cell proliferation assays with autologous antigen-presenting cells (APC) and artificial APC expressing combinations of bovine DR and DQ molecules. Twenty immunostimulatory peptides were identified; three representing two or more epitopes in VirB9-1, ten representing eight or more epitopes in VirB9-2 and seven representing seven or more epitopes in VirB10. Of the eight DRA/DRB3 molecules, four presented 15 peptides, which was biased as DRA/DRB3*1201 presented ten and DRA/DRB3*1101 presented four peptides. Four DQA/DQB molecules composed of two intrahaplotype and two interhaplotype pairs presented seven peptides, of which five were uniquely presented by DQ molecules. In addition, three functional mixed isotype (DQA/DRB3) restriction elements were identified. The immunogenicity and broad MHC class II presentation of multiple VirB9-1, VirB9-2 and VirB10 peptide epitopes justify their testing as a multiepitope vaccine against A. marginale.     

MICA-A5.1 allele is associated with atypical forms of celiac disease in HLA-DQ2-negative patients

We selected 38 consecutive celiac disease (CD) patients (from a group of 316 consecutive CD patients) and 91 healthy blood donors, all of whom were HLA-DQ2 (DQA1*0501/DQB1*0201) negative, and investigated the presence of the classically associated alleles HLA-DQ8 and HLA-DRB4. We also studied the distribution of MICA transmembrane alleles in the two clinical forms of the disease. For this reason, these 38 DQ2-negative patients were subdivided into two groups: 18 typical CD patients and 20 atypical CD patients. No differences were found in the distribution of the DRB4 allele between DQ2-negative patients and controls. The HLA-DQ8 heterodimer (DQA1*03xx/DQB1*0302) was increased in CD patients (29%) compared with controls (10%), but no statistical differences were found. No differences were observed in the frequency of these alleles between either group of CD DQ2-negative patients. MICA-A5.1 was increased in atypical CD patients when compared with the typical forms of disease ( P(c)=0.03) and with healthy controls (P(c)=0.002). No other MICA allele was found to be significantly increased in the groups under study. The presence of MICA-A5.1 in atypical CD DQ2-negative patients may indicate a possible role of this allele in the development of CD.     

Patterns of Adaptive and Neutral Diversity Identify the Xiaoxiangling Mountains as a Refuge for the Giant Panda

Genetic variation plays a significant role in maintaining the evolutionary potential of a species. Comparing the patterns of adaptive and neutral diversity in extant populations is useful for understanding the local adaptations of a species. In this study, we determined the fine-scale genetic structure of 6 extant populations of the giant panda (Ailuropoda melanoleuca) using mtDNA and DNA fingerprints, and then overlaid adaptive variations in 6 functional Aime-MHC class II genes (DRA, DRB3, DQA1, DQA2, DQB1, and DQB2) on this framework. We found that: (1) analysis of the mtDNA and DNA fingerprint-based networks of the 6 populations identified the independent evolutionary histories of the 2 panda subspecies; (2) the basal (ancestral) branches of the fingerprint-based Sichuan-derived network all originated from the smallest Xiaoxiangling (XXL) population, suggesting the status of a glacial refuge in XXL; (3) the MHC variations among the tested populations showed that the XXL population exhibited extraordinary high levels of MHC diversity in allelic richness, which is consistent with the diversity characteristics of a glacial refuge; (4) the phylogenetic tree showed that the basal clades of giant panda DQB sequences were all occupied by XXL-specific sequences, providing evidence for the ancestor-resembling traits of XXL. Finally, we found that the giant panda had many more DQ alleles than DR alleles (33∶13), contrary to other mammals, and that the XXL refuge showed special characteristics in the DQB loci, with 7 DQB members of 9 XXL-unique alleles. Thus, this study identified XXL as a glacial refuge, specifically harboring the most number of primitive DQB alleles.