Gene duplications and sequence polymorphism of bovine class II DQB genes

The genetic diversity of bovine class II DQB genes was investigated by polymerase chain reaction amplification and DNA sequencing. The first domain exon was amplified from genomic DNA samples representing 14 class II haplotypes, defined by restriction fragment length polymorphism (RFLP) analysis. The presence of a polymorphism in the copy number of DQB genes was confirmed since two DQB sequences were isolated from certain haplotypes. Four subtypes of bovine DQB genes were found. DQB1 is the major type and was found in almost all haplotypes. DQB2 is very similar DQB1 but was found only in the duplicated haplotypes DQ9 to 12. DQB3 and DQB4 are two quite divergent genes only present in certain duplicated haplotypes. The bovine DQB complexity thus resembles that in the human DRB region. Bovine DQB genes were found to be highly polymorphic as ten DQB1 alleles and four DQB2 alleles were identified. The observed sequence polymorphism correlated well with previously defined DQB RFLPs. Bovine and human DQB alleles show striking similarities at the amino acid level. In contrast, the frequency of silent substitutions is much higher in comparisons of DQB alleles between species than within species ruling out the possibility that any of the contemporary DQB alleles have been maintained since the divergence of humans and cattle. The frequency of silent substitutions between DQB alleles was markedly lower in cattle than in humans, in agreement with a previous comparison of human and bovine DRB alleles.     

Cloning and sequence analysis of 14 DRB alleles of the bovine major histocompatibility complex by using the polymerase chain reaction

The genetic diversity in the first domain exon of a bovine class II DRB gene was investigated by PCR amplification and DNA sequencing. Genomic DNA samples representing 14 different class II haplotypes, defined by RFLP analysis, were used. The analysis revealed an extensive polymorphism and 14 alleles at a single locus, designated DRB3, were identified. Multiple amino acid substitutions were found in all pairwise comparisons of alleles; 5 to 21 substitutions in the 83 positions compared. The genetic diversity at the amino acid level found in cattle matches the one previously found in the DRB1 locus in man. The significantly higher frequency of replacement substitutions compared with the frequency of silent substitutions provides strong evidence that there is selection for genetic diversity in the bovine DRB3 first domain exon. A comparison of the DRB polymorphism in man and cattle reveals a striking similarity as regards the location of polymorphic positions in the DRB molecule and the degree of polymorphism at polymorphic positions. The majority of polymorphic positions in both species are found in the proposed antigen recognition site of the class II molecule. In addition, there are eight positions which are polymorphic in both species but have not been assigned to the antigen recognition site. The possible functional significance of the polymorphism of these latter positions is discussed.     

Duplicated DQ haplotypes increase the complexity of restriction element usage in cattle

The MHC of cattle encodes two distinct isotypes of class II molecules, DR and DQ. Unlike humans, cattle lack the DP locus and about half the common haplotypes express duplicated DQ genes. The number and frequency of DQA and DQB alleles means that most cattle are heterozygous. If inter- and/or intrahaplotype pairing of DQA and DQB molecules occurs, cattle carrying DQ-duplicated haplotypes may express more restriction elements than would be predicted by the number of expressed alleles. We are investigating whether duplicated haplotypes cause differences in immune response, particularly in terms of generating protective immunity. We have analyzed the Ag-presenting function of DQ molecules in two heterozygous animals, one of which carries a duplicated haplotype. We compared the class II isotype specificity of T cell clones recognizing a putative vaccinal peptide from foot-and-mouth disease virus (FMDV15). We show for the first time that bovine T cells can recognize Ag in the context of DQ molecules. We also present evidence that interhaplotype pairings of DQA and DQB molecules form functional restriction elements. Both animals showed distinct biases to usage of particular restriction elements. Mainly DQ-restricted clones were derived from the animal with duplicated DQ genes, whereas the majority of clones from the animal with a single DQ gene pair were DR restricted. Furthermore, haplotype bias was observed with both animals. These experiments show that understanding of class II chain pairing in addition to knowledge of the genotype may be important in vaccine design where effective epitope selection is essential.     

Sequence and PCR-RFLP analysis of 14 novel BoLA-DRB3 alleles

The genetic diversity of the bovine class IIDRB3 locus was investigated by polymerase chain reaction (PCR) amplification and DNA sequencing of the first domain exon. Studying 34 animals of various cattle breeds, 14 previously unrecognized DRB3 alleles were identified. In three alleles, amino acid substitutions were observed that had not been previously found in bovine DRB3, but occurred at the same position in bovine DQB and in the DRB alleles of other mammals. For all newly identified alleles, the restriction fragment length polymorphism (RFLP) patterns of PCR products obtained with the enzymes Rsa I, Bst YI, and Hae III were compared with patterns of 38 previously described alleles. Altogether, eleven novel PCR-RFLP types were defined. Twelve out of the 42 PCR-RFLP types identified so far were not found to be fully informative because they corresponded to more than one allelic sequence. PCR-RFLP may therefore be a rapid and useful method for DRB3 typing in cattle families, but for studies on outbred populations, sequencing and hybridization techniques are required.     

Evolution ofMHC polymorphism: Extensive sharing of polymorphic sequence motifs between human and bovine DRB alleles

The evolution ofMHC polymorphism has been studied by comparing the amino acid and nucleotide sequences of 14 bovine and 32 humanDRB alleles. The comparison revealed an extensive sharing of polymorphic sequence motifs in the two species. Almost identical sets of residues were found at several highly polymorphic amino acid positions in the putative antigen recognition site. Consequently, certain bovine alleles were found to be more similar to certain human alleles than to other bovine alleles. In contrast, the frequencies of silent nucleotide substitutions were found to be much higher in comparisons between species than within species implying that none of the human or bovine DRB alleles originated before the divergence of these distantly related species. The results suggest that the observed similarity inDRB polymorphism is due to convergent evolution and possibly the sharing of short ancestral sequence motifs. However, the relative role of the latter mechanism is difficult to assess due to the biased base composition in the first domain exon of polymorphic class 11 genes. The frequency of silent substitutions betweenDRB alleles was markedly lower in cattle than in man suggesting that theDRB diversity has evolved more rapidly in the former species.     

Characterization of Bison bison major histocompatibility complex class IIa haplotypes

American bison (Bison bison) and domestic cattle (Bos taurus and Bos indicus) evolved from a common ancestor 1–1.4 million years ago. Nevertheless, they show dramatic differences in their susceptibility to infectious diseases, including malignant catarrhal fever (MCF). Although bison are highly susceptible to ovine herpesvirus-2 (OvHV-2) associated MCF, about 20% of healthy domesticated and wild bison are positive for OvHV-2 antibody. We are interested in testing the hypothesis that, within the bison population, the polymorphism of major histocompatibility complex (MHC) class II genes influences resistance to MCF. However, since little was known about the MHC class II genes of bison, it was necessary to first characterize class II haplotypes present in Bi. bison (Bibi). Thus, the MHC class II haplotypes carried by 14 bison were characterized by the PCR-based cloning and sequencing of their DRB3, DQA, and DQB alleles. Twelve MHC class II haplotypes were identified in the 14 bison. These haplotypes comprised six previously reported and six new Bibi-DRB3 alleles, along with 11 Bibi-DQA and 10 Bibi-DQB alleles. For each bison class II allele, it was possible to identify closely related cattle sequences. The closest bison and bovine DQA, DQB, and DRB3 alleles, on average, differed by only 1.3, 3.5, and 5.8 amino acids, respectively. Furthermore, bison MHC haplotypes with both nonduplicated and duplicated DQ genes were identified; these haplotypes appear to have originated from the same ancestral haplotypes as orthologous cattle haplotypes. This study was supported by USDA-Agricultural Research Service grant CWU-5348-32000-018-00D. While working on this project, Dr. Bharat Bhushan was supported by a fellowship from the World-Bank-sponsored National Agricultural Technology Project of the Indian Council of Agricultural Research, Indian Ministry of Agriculture, New Delhi, India     

Visualizing human leukocyte antigen class II risk haplotypes in human systemic lupus erythematosus

Human leukocyte antigen (HLA) class I and class II alleles are implicated as genetic risk factors for many autoimmune diseases. However, the role of the HLA loci in human systemic lupus erythematosus (SLE) remains unclear. Using a dense map of polymorphic microsatellites across the HLA region in a large collection of families with SLE, we identified three distinct haplotypes that encompassed the class II region and exhibited transmission distortion. DRB1 and DQB1 typing of founders showed that the three haplotypes contained DRB1*1501/ DQB1*0602, DRB1*0801/ DQB1*0402, and DRB1*0301/DQB1*0201 alleles, respectively. By visualizing ancestral recombinants, we narrowed the disease-associated haplotypes containing DRB1*1501 and DRB1*0801 to an approximately 500-kb region. We conclude that HLA class II haplotypes containing DRB1 and DQB1 alleles are strong risk factors for human SLE.     

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.     

Polymorphism of the HLA class II loci in Siberian populations

The populations that colonized Siberia diverged from one another in the Paleolithic and evolved in isolation until today. These populations are therefore a rich source of information about the conditions under which the initial divergence of modern humans occurred. In the present study we used the HLA system, first, to investigate the evolution of the human major histocompatibility complex (MHC) itself, and second, to reveal the relationships among Siberian populations. We determined allelic frequencies at five HLA class II loci (DRB1, DQA1, DQB1, DPA1, and DPB1) in seven Siberian populations (Ket, Evenk, Koryak, Chukchi, Nivkh, Udege, and Siberian Eskimo) by the combination of single-stranded conformational polymorphism and DNA sequencing analysis. We then used the gene frequency data to deduce the HLA class II haplotypes and their frequencies. Despite high polymorphism at four of the five loci, no new alleles could be detected. This finding is consistent with a conserved evolution of human class II MHC genes. We found a high number of HLA class II haplotypes in Siberian populations. More haplotypes have been found in Siberia than in any other population. Some of the haplotypes are shared with non-Siberian populations, but most of them are new, and some represent “forbidden” combinations of DQA1 and DQB1 alleles. We suggest that a set of “public” haplotypes was brought to Siberia with the colonizers but that most of the new haplotypes were generated in Siberia by recombination and are part of a haplotype pool that is turning over rapidly. The allelic frequencies at the DRB1 locus divide the Siberian populations into eastern and central Siberian branches; only the former shows a clear genealogical relationship to Amerinds. Received: 18 August 1997 / Accepted: 6 October 1997     

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.     

Extensive polymorphism of the BOLA-DRB3 gene distinguished by PCR-RFLP

A polymerase chain reaction (PCR)-based method is described for typing of alleles of the bovine lymphocyte antigen (BoLA)-DRB3 gene. A total of 30 DRB3 alleles were distinguished by digestion of PCR amplification products of BoLA-DRB3 exon 2 with RsaI, BstYI and HaeIII (PCR-RFLP). All restriction fragment patterns, with the exception of one HaeIII pattern, were consistent with restriction sites that were found among 14 previously sequenced DRB3 alleles. The PCR-RFLP typing method was evaluated on 168 genomic DNA samples collected from animals of 10 cattle breeds, 48 of which were typed in the Fourth International BoLA Workshop for BoLA-DRB and -DQ by conventional restriction fragment length polymorphism (RFLP) analysis using heterologous and homologous DNA probes. Thirty-one DRB/DQ haplotypes containing 23 DRB3 alleles were identified among the 48 workshop animals analysed. Using PCR-RFLP, 11 DRB3 alleles were identified in 18 workshop animals for which DRB RFLPs were not informative. PCR-RFLP typing of additional animals revealed five new DRB3 alleles, of which three contained a putatively located three basepair deletion in the identical position as found for the sequenced allele DRB*2A. PCR-RFLP was shown to be a rapid and sensitive method for the detection of polymorphism in a functionally relevant domain of the BoLA-DRB3 gene and should be useful for studying the evolution of DRB polymorphism in cattle and other Bovidae.     

HLA class II antigens positively and negatively associated with hepatosplenic schistosomiasis in a Chinese population

To identify possible associations between host genetic factors and the onset of liver fibrosis following Schistosoma japonicum infection, the major histocompatibility class II alleles of 84 individuals living on an island (Jishan) endemic for schistosomiasis japonica in the Poyang Lake Region of Southern China were determined. Forty patients exhibiting advanced schistosomiasis, characterised by extensive liver fibrosis, and 44 age and sex-matched control subjects were assessed for the class II haplotypes HLA-DRB1 and HLA-DQB1. Two HLA-DRB1 alleles, HLA-DRB1*0901 (P=0.012) and *1302 (P=0.039), and two HLA-DQB1 alleles, HLA-DQB1*0303 (P=0.012) and *0609 (P=0.037), were found to be significantly associated with susceptibility to fibrosis. These associated DRB1 and DQB1 alleles are in very strong linkage disequilibrium, with DRB1*0901-DQB1*0303 and DRB1*1302-DQB1*0609 found as common haplotypes in this population. In contrast, the alleles HLA-DRB1*1501 (P=0.025) and HLA-DQB1*0601 (P=0.022) were found to be associated with resistance to hepatosplenic disease. Moreover, the alleles DQB1*0303 and DRB1*0901 did not increase susceptibility in the presence of DQB1*0601, indicating that DQB1*0601 is dominant over DQB1*0303 and DRB1*0901. The study has thus identified both positive and negative associations between HLA class II alleles and the risk of individuals developing moderate to severe liver fibrosis following schistosome infection.     

Mhc-DRB genes of the pigtail macaque (Macaca nemestrina): implications for the evolution of human DRB genes.

The DRB family of human class II major histocompatibility complex (Mhc) loci is unusual in that individuals differ in the number and combination of genes (haplotypes) they carry. Indications are that both the allelic and haplotype polymorphisms of the DRB loci predate speciation. Searching for the evolutionary origins of these polymorphisms, we have sequenced five DRB clones isolated from a cDNA library of a pigtail macaque (Macaca nemestrina) B lymphocyte line. The clones represent five different genes which we designate Mane-DRB*01-Mane-DRB*05. The genes appears to be approximately equidistant from each other, so that allelic relationships between them cannot be established on the basis of the sequence data alone. If positions coding for the peptide-binding region of the class II beta chains are eliminated from sequence comparisons, the Mane-DRB genes appear to be most closely related to the human (HLA) DRB1 genes of the DRw52 group. We interpret this finding to indicate that the ancestral gene of the DRw52 group of human DRB1 alleles separated from the rest of the HLA-DRB1 alleles before the separation of the Old World monkeys (Cercopithecoidea) from the apes (Hominoidea) in the early Oligocene. After this separation, the ancestral DRB1 gene of the DRw52 group duplicated in the Old World monkey lineage to give rise to genes at three loci at least, while in the ape lineage this gene may have remained single and diverged into a number of alleles instead. These findings suggest that some of the polymorphism currently present at the DRB1 locus is greater than 35 Myr old.     

Two divergent routes of evolution gave rise to the DRw13 haplotypes

The HLA class II genes and haplotypes have evolved over a long period of evolutionary time by mechanisms such as gene conversion, reciprocal recombination and point mutation. The extent of the diversity generated is most clearly evident in an analysis of the HLA class II alleles present within DRw13 haplotypes. This study uses cDNA sequencing to examine the first domains of DRB1, DRB3, DQA1, and DQB1 alleles from several American black individuals expressing seven different DRw13 haplotypes, five with undefined HLA-D specificities (i.e., not Dw18 or Dw19). Two new DRw13 alleles described in this study are the first examples of convergent evolution of DR alleles in which gene conversion has apparently combined segments of DRB1 alleles encoding DRw11 and DRw8 to generate two new DRB1 alleles, DRB1*1303 and DRB1*1304, that encode molecules bearing serologic determinants of a third allele, DRw13. These new DRw13 alleles are found embedded in haplotypes of DRw11 origin distinct from haplotypes encoding previously identified DRw13 alleles, DRB1*1301 and DRB1*1302. These data suggest that two evolutionary pathways may have given rise to two subgroups of alleles encoding molecules that share DRw13 serologic determinants yet which possess different structural and, likely, functional motifs. Reciprocal gene recombination events resulting in different DR, DRw52 and DQ allele combinations also appear to have played a crucial role in augmenting the level of diversity found in DRw13 haplotypes. Recombination has resulted in the association of one of the new DRw13 alleles with a DQw2 allele normally found associated with DR7 and the association of the DRw52c-associated DRw13 allele (DRB1*1302) with three different DQw1 alleles. The seven DRw13 haplotypes that have resulted from the effect of recombination on haplotypes formed by the two pathways of DRw13 allelic diversification have resulted in different repertoires of class II molecules and, most likely, different immune response profiles in individuals with these haplotypes.     

Characterization of swine leukocyte antigen alleles and haplotypes on a novel miniature pig line,Microminipig

Microminipigs are extremely small‐sized, novel miniature pigs that were recently developed for medical research. The inbred Microminipigs with defined swine leukocyte antigen (SLA) haplotypes are expected to be useful for allo‐ and xenotransplantation studies and also for association analyses between SLA haplotypes and immunological traits. To establish SLA‐defined Microminipig lines, we characterized the polymorphic SLA alleles for three class I (SLA‐1, SLA‐2 and SLA‐3) and two class II (SLA‐DRB1 and SLA‐DQB1) genes of 14 parental Microminipigs using a high‐resolution nucleotide sequence‐based typing method. Eleven class I and II haplotypes, including three recombinant haplotypes, were found in the offspring of the parental Microminipigs. Two class I and class II haplotypes, Hp‐31.0 (SLA‐1*1502–SLA‐3*070102–SLA‐2*1601) and Hp‐0.37 (SLA‐DRB1*0701–SLA‐DQB1*0502), are novel and have not so far been reported in other pig breeds. Crossover regions were defined by the analysis of 22 microsatellite markers within the SLA class III region of three recombinant haplotypes. The SLA allele and haplotype information of Microminipigs in this study will be useful to establish SLA homozygous lines including three recombinants for transplantation and immunological studies.     

PCR/oligonucleotide probe typing of HLA class II alleles in a Filipino population reveals an unusual distribution of HLA haplotypes.

We have analyzed the distribution of HLA class II alleles and haplotypes in a Filipino population by PCR amplification of the DRB1, DQB1, and DPB1 second-exon sequences from buccal swabs obtained from 124 family members and 53 unrelated individuals. The amplified DNA was typed by using nonradioactive sequence-specific oligonucleotide probes. Twenty-two different DRB1 alleles, including the novel Filipino *1105, and 46 different DRB1/DQB1 haplotypes, including the unusual DRB1*0405-DQB1*0503, were identified. An unusually high frequency (f = .383) of DPB1*0101, a rare allele in other Asian populations, was also observed. In addition, an unusual distribution of DRB1 alleles and haplotypes was seen in this population, with DR2 (f = .415) and DRB1*1502-DQB1*0502 (f = .233) present at high frequencies. This distribution of DRB1 alleles differs from the typical HLA population distribution, in which the allele frequencies are more evenly balanced. The distribution of HLA class II alleles and haplotypes in this Filipino population is different from that of other Asian and Pacific groups: of those populations studied to date; the Indonesian population is the most similar. DRB1*1502-DQB1*0502 was in strong linkage disequilibrium (D'' = .41) with DPB1*0101 (f = .126, for the extended haplotype), which is consistent with selection for this DR, DQ, DP haplotype being responsible for the high frequency of these three class II alleles in this population.     

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.     

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.