The nucleotide sequence of bovine MHC class II DQB and DRB genes

The DRB3 gene of the bovine major histocompatibility complex: discovery, diversity and distribution of alleles in commercial breeds of cattle and applications for development of vaccines

The bovine Major Histocompatibility Complex (MHC), also known as the Bovine Leucocyte Antigen (BoLA) complex, is the genomic region that encodes the most important molecules for antigen presentation to initiate immune responses. The first evidence of MHC in bovines pointed to a locus containing 2 antigens, one detected by cytotoxic antiserum (MHC class I) and another studied by mixed lymphocyte culture tests (MHC class II). The most studied gene in the BoLA region is the highly polymorphic BoLA-DRB3, which encodes a β chain with a peptide groove domain involved in antigen presentation for T cells that will develop and co-stimulate cellular and humoral effector responses. BoLA-DRB3 alleles have been associated with outcomes in infectious diseases such as mastitis, trypanosomiasis, and tick loads, and with production traits. To catalog these alleles, 2 nomenclature methods were proposed, and the current use of both systems makes it difficult to list, comprehend and apply these data effectively. In this review we have organized the knowledge available in all of the reports on the frequencies of BoLA-DRB3 alleles. It covers information from studies made in at least 26 countries on more than 30 breeds; studies are lacking in countries that are important producers of cattle livestock. We highlight practical applications of BoLA studies for identification of markers associated with resistance to infectious and parasitic diseases, increased production traits and T cell epitope mapping, in addition to genetic diversity and conservation studies of commercial and creole and locally adapted breeds. Finally, we provide support for the need of studies to discover new BoLA alleles and uncover unknown roles of this locus in production traits.

Genetic Variation and Population Differentiation in the Bovine Lymphocyte Antigen DRB3.2 Locus of South African Nguni Crossbred Cattle

Simple Summary

Indigenous cattle breeds are important to their native environment as they confer significant and often unique adaptability traits. The Nguni is one such breeds that is indigenous to the Southern African region. This breed confers resistance to thermal stress and diseases, amongst other factors. The bovine major histocompatibility complex is an important region, which codes for alleles that have been associated with a plethora of diseases. In the current study, the genetic diversity within this region was assessed in Nguni crossbred cattle. This was done to detect the gene pool of the Nguni breed, and to identify genes that might be important within this breed. The populations displayed a high degree of genetic diversity, and some alleles were common throughout the populations and accounted for a significant portion of the total alleles. This high genetic diversity could account for the great adaptability of the Nguni breed to Southern Africa.

Abstract

The bovine lymphocyte antigen (BoLA-DRB3) gene is an important region that codes for glycoproteins responsible for the initiation of an immune response. BoLA-DRB3 alleles have been demonstrated to be associated with disease resistance/tolerance. Therefore, great genetic diversity is correlated with better adaptation, fitness, and robustness. The current study was conducted to assess the population genetic structure of the BoLA-DRB3 gene in Nguni crossbred cattle using polymerase chain reaction-sequence based typing (PCR-SBT). High genetic diversity was detected, with 30 alleles, 11 of which are novel to the study. Alleles DRB3*0201, DRB3*0701, DRB*0901, and DRB*1601 were present in all populations and accounted for nearly around 50% of all observed alleles. A mean genetic diversity (H_E) of 0.93 was detected. The high overall genetic diversity is possibly associated with pathogen-assisted selection and heterozygote advantage. Such high diversity might explain the hardiness of the Nguni crossbred cattle to the Southern African region. Low population genetic structure was identified (F_ST = 0.01), suggesting possible gene flow between populations and retention of similar alleles. The study was undertaken to bridge the dearth of such studies in South African breeds and it is imperative for effective sustainability of indigenous breeds and the implementation of effective breeding strategies.

Genetic assessment of BoLA-DRB3 polymorphisms by comparing Bangladesh, Ethiopian, and Korean cattle

Attributable to their major function in pathogen recognition, the use of bovine leukocyte antigens (BoLA) as disease markers in immunological traits in cattle is well established. However, limited report exists on polymorphism of the BoLA gene in zebu cattle breeds by high resolution typing methods. Thus, we used a polymerase chain reaction sequence-based typing (PCR-SBT) method to sequence exon 2 of the BoLA class II DRB3 gene from 100 animals (Boran, n = 13; Sheko, n = 20; Fogera, n = 16; Horro, n = 19), Hanwoo cattle (n = 18) and Bangladesh Red Chittagong zebu (n = 14). Out of the 59 detected alleles, 43 were already deposited under the Immuno Polymorphism Database for major histocompatibility complex (IPD-MHC) while 16 were unique to this study. Assessment of the level of genetic variability at the population and sequence levels with genetic distance in the breeds considered in this study showed that Zebu breeds had a gene diversity score greater than 0.752, nucleotide diversity score greater than 0.152, and mean number of pairwise differences higher than 14, being very comparable to those investigated for other cattle breeds. Regarding neutrality tests analyzed, we investigated that all the breeds except Hanwoo had an excess number of alleles and could be expected from a recent population expansion or genetic hitchhiking. Howbeit, the observed heterozygosity was not significantly (p < 0.05) higher than the expected heterozygosity. The Hardy Weinberg equilibrium (HWE) analysis revealed non-significant excess of heterozygote animals, indicative of plausible over-dominant selection. The pairwise FST values suggested a low genetic variation among all the breeds (FST = 0.056; p < 0.05), besides the rooting from the evolutionary or domestication history of the cattle. No detached clade was observed in the evolutionary divergence study of the BoLA-DRB3 gene, inferred from the phylogenetic tree based on the maximum likelihood model. The investigation herein indicated the clear differences in BoLA-DRB3 gene variability between African and Asian cattle breeds.

The Tip Region on VP2 Protein of Bluetongue Virus Contains Potential IL-4-Inducing Amino Acid Peptide Segments

Bluetongue is an infectious viral hemorrhagic disease of domestic and wild ruminants that has a considerable economic impact on domestic ruminants. There are currently at least 29 serotypes of bluetongue virus (BTV) in the world. Noteworthily, the pathogenesis among BTV serotypes is different, even in the same animal species. In this study, BTV2/KM/2003 and BTV12/PT/2003 were used to investigate the differential immunological effects on bovine peripheral blood mononuclear cells (PBMCs). The BTV viral load and the expression of cytokine messenger RNA (mRNA) in PBMCs were measured by fluorescence-based real-time reverse-transcription PCR (qRT-PCR). The immunofluorescence assay (IFA) was applied to detect BTV signals in monocyte-derived macrophages (MDMs). The SWISS-MODEL and IL-4pred prediction tools were used to predict the interleukin 4 (IL-4)-inducing peptides in BTV-coat protein VP2. Synthetic peptides of VP2 were used to stimulate PBMCs for IL-4-inducing capability. This study demonstrated that the cytokine profiles of BTV-induced PBMCs were significantly different between BTV2/KM/2003 and BTV12/PT/2003. BTV2 preferentially activated the T helper 2 (Th2) pathway, represented by the early induction of IL-4, and likely fed back to inhibit the innate immunity. In contrast, BTV12 preferentially activated the innate immunity, represented by the induction of tumor necrosis factor -α (TNF-α) and interleukin 1 (IL-1), with only minimal subsequent IL-4. The BTV nonstructural protein 3 antibody (anti-BTV-NS3) fluorescent signals demonstrated that monocytes in PBMCs and MDMs were the preferred targets of BTV replication. Bioinformatics analysis revealed that the capability to induce IL-4 was attributed to the tip region of the VP2 protein, wherein a higher number of predicted peptide segments on BTVs were positively correlated with the allergic reaction reported in cattle. Synthetic peptides of BTV2-VP2 induced significant IL-4 within 12-24 h post-infection (hpi) in PBMCs, whereas those of BTV12 did not, consistent with the bioinformatics prediction. Bovine PBMCs and synthetic peptides together seem to serve as a good model for pursuing the BTV-induced IL-4 activity that precedes the development of an allergic reaction, although further optimization of the protocol is warranted.

Genetic characterization of MHC class II DQB exon 2 variants in gayal (Bos frontalis)

In the present study, exon 2 of major histocompatibility complex (MHC) class II DQB gene from 39 gayals (Bos frontalis) was isolated, characterized and compared with previously reported patterns for other bovidae. It was revealed by sequence analyses that there are 36 DQB exon 2 variants among 39 gayals. These variants exhibited a high degree of nucleotide and amino acid substitutions with most amino acid variations occurring at positions forming the peptide-binding sites (PBS). The DQB loci were analysed for patterns of synonymous (d_S) and non-synonymous (d_N) substitution. The gayals were observed to be under strong balancing selection in the DQB exon 2 PBS (d_N = 0.094, P = 0.001). It appears that this variability among gayals could confer the ability to mount immune responses to a wide variety of peptides or pathogens.

Genetic variability of MHC class II DQB exon 2 alleles in yak (Bos grunniens)

The major histocompatibility class (MHC) DQ molecules are dimeric glycoproteins revealing antigen presentation to CD(4+) T cells. In the present study, the exon 2 of the MHC class II DQB gene from 32 yaks (Bos grunniens) was cloned, sequenced and compared with previously reported patterns for other bovidae. It was revealed by sequence analyses that there are 25 DQB exon 2 alleles among 32 yaks, all alleles are found to belong to DQB1 loci. These alleles exhibited a high degree of nucleotide and amino acid polymorphisms with most amino acid variations occurring at positions forming the peptide-binding sites. The DQB loci were analyzed for patterns of synonymous (d S) and non-synonymous (d N) substitution. The yak was observed to be under strong positive selection in the DQB exon 2 peptide-binding sites (d N = 0.15, P < 0.001). It appears that this variability among yaks confers the ability to mount immune responses to a wide variety of peptides or pathogens.

A high-resolution radiation hybrid map of the river buffalo major histocompatibility complex and comparison with BoLA

The major histocompatibility complex (MHC) in mammals codes for antigen-presenting proteins. For this reason, the MHC is of great importance for immune function and animal health. Previous studies revealed this gene-dense and polymorphic region in river buffalo to be on the short arm of chromosome 2, which is homologous to cattle chromosome 23. Using cattle-derived STS markers and a river buffalo radiation hybrid (RH) panel (BBURH5000 ), we generated a high-resolution RH map of the river buffalo MHC region. The buffalo MHC RH map (cR5000 ) was aligned with the cattle MHC RH map (cR12000 ) to compare gene order. The buffalo MHC had similar organization to the cattle MHC, with class II genes distributed in two segments, class IIa and class IIb. Class IIa was closely associated with the class I and class III regions, and class IIb was a separate cluster. A total of 53 markers were distributed into two linkage groups based on a two-point LOD score threshold of ≥8. The first linkage group included 32 markers from class IIa, class I and class III. The second linkage group included 21 markers from class IIb. Bacterial artificial chromosome clones for seven loci were mapped by fluorescence in situ hybridization on metaphase chromosomes using single- and double-color hybridizations. The order of cytogenetically mapped markers in the region corroborated the physical order of markers obtained from the RH map and served as anchor points to align and orient the linkage groups.

The major histocompatibility complex in bovines: a review

Productivity in dairy cattle and buffaloes depends on the genetic factors governing the production of milk and milk constituents as well as genetic factors controlling disease resistance or susceptibility. The immune system is the adaptive defense system that has evolved in vertebrates to protect them from invading pathogens and also carcinomas. It is remarkable in the sense that it is able to generate an enormous variety of cells and biomolecules which interact with each other in numerous ways to form a complex network that helps to recognize, counteract, and eliminate the apparently limitless number of foreign invading pathogens/molecules. The major histocompatibility complex which is found to occur in all mammalian species plays a central role in the development of the immune system. It is an important candidate gene involved in susceptibility/resistance to various diseases. It is associated with intercellular recognition and with self/nonself discrimination. It plays major role in determining whether transplanted tissue will be accepted as self or rejected as foreign.

Allelic diversity at MHC class II DQ loci in buffalo (Bubalus bubalis): evidence for duplication

The genetic diversity of MHC class II DQ genes was investigated in riverine buffalo (Bubalus bubalis) by PCR-RFLP and sequencing. Highly variable regions (exons 2-3) of DQ genes were amplified from 152 buffaloes and genotyped by PCR-RFLP. Alleles identified by differential restriction patterns were sequenced for the characterization. PCR-RFLP was a rapid method to discriminate between DQA1 and duplicated DQA2 genes in buffalo, however, the method appeared to be inadequate for determining the more complicated DQB genotypes. A total of 7 and 10 alleles were identified for DQA and DQB loci, respectively. Nucleotide as well as amino acid variations among DQ alleles particularly at peptide binding regions were high. Such variations were as expected higher in DQB than DQA alleles. The phylogenetic analysis for both genes revealed the grouping of alleles into two major sub-groups with higher genetic divergence. High divergence among DQ allelic families and the isolation of two diverse DQA and DQB sequences from individual samples indicated duplication of DQ loci was similar in buffalo to other ruminants.

Sequencing and genetic analysis of a bovine DQB cDNA clone

A BoLA-DQB cDNA clone (BoLA-DQ beta-1) was isolated by screening a bovine lymphoblastoid cDNA library with a HLA-DQB genomic clone. The DNA and predicted protein sequences were compared to class II sequences from cattle and other species. BoLA-DQ beta-1 has 92.0% similarity to the coding regions of two previously sequenced BoLA-DQB genomic clones and 69.6% similarity to a BoLA-DR beta pseudogene. However, the first domain encoded by BoLA-DQ beta-1 has 94 amino acids; one more than the predicted size of the products encoded by two previously sequenced bovine DQB genes (BoDQ beta-Q1 and BoDQ beta-Y1). Comparing all coding regions, BoLA-DQ beta-1 has greater nucleotide similarity to HLA-DQB sequences than to I-A beta, HLA-DRB and I-E beta sequences. Like the HLA-DQB gene product, the cytoplasmic domain of the predicted protein encoded by BoLA-DQ beta-1 is eight amino acids shorter than that of I-A beta, HLA-DRB and I-E beta molecules. Six clone-specific amino acid substitutions were identified in the beta 1 domain of BoLA-DQ beta-1, including an unusual cysteine residue at position 13 which is believed to be positioned on a beta-strand and face into the antigen recognition site. Southern blot analysis of PvuII-digested genomic DNA from a paternal half-sibling family (sire, and six dam-offspring pairs) using BoLA-DQ beta-1 as a probe, revealed five allelic PvuII RFLP patterns, including two patterns not previously described, that cosegregated with serologically-defined BoLA-A (class I) alleles. The evolution, polymorphism and function of a transcriptionally active BoLA-DQB gene can now be readily studied using this DQB cDNA clone as a source of allele and locus-specific oligonucleotide primers.

Nucleotide sequence and northern analysis of a bovine major histocompatibility class II DR beta-like cDNA

A 1.2-kb bovine DR beta-like cDNA clone (BoLA-DRB3) was isolated from a peripheral blood lymphocyte cDNA library utilizing a human DR beta cDNA as a probe. BoLA-DRB3 was found to have a high degree of nucleotide sequence similarity (96.8%) with a previously sequenced bovine DR beta-like gene (A1). It is believed that BoLA-DRB3 and A1 represent distinct alleles of one of the three bovine DR beta-like loci. Sequence comparison of BoLA-DRB3 with genes representing the other two bovine DR beta-like loci resulted in moderate degrees of sequence similarities (83.1% and 86.3%, respectively). Comparison of the relative abundance of RNA transcripts of the three bovine DR beta-like loci by Northern analysis of lymphocyte RNA indicated that BoLA-DRB3 is the most actively transcribed of the three bovine DR beta-like genes. Based on these results we suggest that of the three DR beta-like loci thus far identified in the bovine, only one is actively transcribed.

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.

Class II major histocompatibility complex genes of the sheep

The class II genes of the sheep major histocompatibility complex (MHC) have been cloned from two unrelated heterozygous sheep into cosmid vectors. By restriction mapping and hybridization with a number of class II probes of human and mouse origin, the cloned genetic material has been assigned to seven distinct alpha genes, 10 distinct beta genes and 14 beta-related sequences. It was difficult to identify homologues of specific HLA class II genes because of a tendency for the ovine genes to cross-hybridize between HLA probes representing different loci. Such cross-hybridization was especially marked among the beta genes. While DQ and DR homologues have been tentatively identified by several criteria, no genes corresponding to DP have been identified. Cosmids containing class II alpha and beta genes have been transfected into mouse LTK- cells, and surface expression of a sheep class II molecule has been obtained.

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.

Joint Report of the Fourth International Bovine Lymphocyte Antigen (BoLA) Workshop, East Lansing, Michigan, USA, 25 August 1990

Blood samples from 54 animals were exchanged between 15 laboratories in nine countries to improve and expand BoLA class I and class II typing. A total of 27 out of 33 (82%) of previously accepted BoLA-w specificities were represented within the cell panel. Seventeen new serum-defined BoLA specificities were accepted by the workshop participants, thus expanding the number of internationally recognized BoLA specificities to 50. The large number of new specificities detected resulted from the number of serological reagents used (n = 1139) and the genetic diversity of the cell panel. Confidence derived from the high percentage of agreement between the laboratories on antigen detection (97.3%; r = 0.84) permitted the removal of the workshop (w) notation from 23 BoLA-w specificities and their acceptance as full status BoLA-A antigens. Two new non-BoLA antigens were also detected, one completely included within the red blood cell factor S' (BoLy-S'), whereas a second (BoLy-w1) did not show any association with tested red blood cell factors. A comparison between serological, isoelectric focusing (IEF) and DNA typing for BoLA class II polymorphism was conducted with a subset of workshop cells. Correlation between the three methods was significant for three combinations of alleles. Three other serologically defined class II specificities were correlated with DR and/or DQ restriction fragment length polymorphism (RFLP) types, whereas six additional IEF types were correlated with DR and/or DQ RFLP types (r greater than or equal to 0.50). Several new IEF, DRB, DQA and DQB RFLP patterns were identified. In 46 animals that were typed for BoLA-DR and DQ genes by RFLP analysis, 46 different BoLA haplotypes were tentatively defined. These 46 haplotypes were distinguished by 31 serologically-defined BoLA-A alleles (and 2 'blanks'), 15 DRB RFLP types (plus up to 10 new DRB RFLP patterns) and 23 DQA-DQB haplotypes.

Expression and characterization of ovine major histocompatibility complex class II (OLA-DR) genes

Previous work made use of nucleic acid probes corresponding to different subtypes of the class II regions of the human and murine major histocompatibility complex (MHC) to isolate seven different alpha and 24 different beta genes of the ovine MHC from two cosmid libraries. In an attempt to identify pairs of alpha and beta genes capable of cell surface expression, all permutations of alpha and beta genes were in turn transfected into mouse L-cells. Two pairs of alpha and beta genes co-expressed and stable ovine MHC class II L-cell lines were developed. The expressed alpha genes had previously been defined as DR-alpha homologues (DRA) by differential Southern hybridization to human subtype specific class II probes. The expressed ovine beta genes were also assigned as ovine DR-beta homologues (DRB) on the basis of their sequence having a higher degree of similarity with human DRB than any other subtype. A total of eight out of 23 anti-sheep class II specific monoclonal antibodies were typed OLA-DR specific by FACScan analysis using the L-cell lines.

Detection of bovine leukocyte antigen DRB3 alleles as candidate markers for clinical mastitis resistance in Holstein x Zebu

Bovine leukocyte antigen DRB3 alleles from Holstein x Zebu crossbred dairy cows (n = 409) were analyzed using the PCR-RFLP technique. Exon II of DRB3 was amplified using locus-specific primers (HLO30/HLO32), followed by digestion with 3 restriction enzymes (RsaI, BstyI, and HaeIII). Forty alleles were found with frequency ranging from 0.005 to 0.139. The most frequently detected alleles of Holstein x Zebu were DRB3*16, *51, *23, *11, *8, and *1, accounting for 61.12% of the alleles in the population. Detection of candidate alleles for clinical mastitis occurrence was performed by logistic regression. It was found that percentage of Holstein fraction in crossbred cows had a nonsignificant effect (P > 0.05). However, parity had a significant effect on mastitis occurrence. In addition, DRB3*1 and *52 were the most associated with the occurrence of clinical mastitis, whereas *15, *51, and *22 were associated with resistance in crossbred populations. This is the first report of association of DRB3*15 and *51 with mastitis resistance. The association was validated by examining the candidate alleles in another commercial population. Highly susceptible (n = 43) and resistant (n = 42) groups of Holstein x Zebu cows were investigated. The result confirmed that DRB3*1 and *52 could be considered as susceptibility alleles, whereas *15, *51, and *22 could be considered as resistant alleles in Holstein x Zebu raised under tropical conditions. In addition, allele effects on 305-d milk production were estimated by BLUP. It was shown that most alleles associated with high clinical mastitis occurrence were related to increased milk yield. This study revealed that allele DRB3*10 had the greatest effect on increasing milk yield with moderate resistance to clinical mastitis, which could be used as a potential marker for selection in dairy genetic evaluation.

Short communication: Identification of allelic variation at the bovine DRA locus by polymerase chain reaction-single strand conformational polymorphism

The major histocompatibility complex DRA locus is noteworthy among the major histocompatibility complex class II loci for the little or no variation reported in many species. In cattle, DRA has not been investigated in depth, and the extent of variation at the locus is not yet known. In this study, we used PCR-single strand conformational polymorphism (SSCP) analysis to screen for potential sequence variation in the second exon of bovine DRA, which encodes the antigen-presentation groove. Four unique SSCP patterns were detected among 384 cattle from New Zealand. Sequence analysis revealed that these SSCP patterns represent 4 DRA alleles; and 3 single nucleotide polymorphisms were detected in exon 2. However, all of the single nucleotide polymorphisms were synonymous, and no amino acid change was therefore expected in this region. The polymorphism detected may be linked to variation elsewhere in the gene that affects its structure or function.

Genomic analysis of Ovis aries (Ovar) MHC class IIa loci

Determining the genomic organization of the Ovis aries (Ovar) major histocompatibility complex class IIa region is essential for future functional studies related to antigen presentation. In this study, a bacterial artificial chromosome (BAC) library of genomic DNA from peripheral blood leukocytes (PBL) of a Rambouillet ram was constructed, and BAC clone consisting of the major histocompatibility complex (MHC) class II DQB2, DQA2, DQB1, DQA1, and DRB1 loci was identified and completely sequenced. The BAC clone consists of 160,889 bp of finished sequence with the loci arranged in the following order: DQB2, DQA2, DQB1, DQA1, and DRB1 with 14.3, 25, 6.6, and 40.9 Kb spanning between the loci, respectively. All five of these loci were transcribed in the animal used to generate the MHC class II BAC clone. Repeat or retrotransposable elements along with MHC class II cis promoter elements consisting of S, X, and Y boxes were identified in the sequence. In addition, 16 non-coding conserved sequences amongst primates, carnivores, and ruminants were identified (p < 0.001). These conserved sequences include binding sites for transcription factors with known roles in immune cells, and they provide a basis for further functional investigation of the genes in this region. This is the first ruminant finished sequence of the DQB2-DRB1 region, and this sequence information will aid in whole genome and transcriptome analyses of MHC class II.

Diversity and duplication of DQB and DRB-like genes of the MHC in baleen whales (suborder: Mysticeti)

The molecular diversity and phylogenetic relationships of two class II genes of the baleen whale major histocompatibility complex were investigated and compared to toothed whales and out-groups. Amplification of the DQB exon 2 provided sequences showing high within-species and between-species nucleotide diversity and uninterrupted reading frames consistent with functional class II loci found in related mammals (e.g., ruminants). Cloning of amplified products indicated gene duplication in the humpback whale and triplication in the southern right whale, with average nucleotide diversity of 5.9 and 6.3%, respectively, for alleles of each species. Significantly higher nonsynonymous divergence at sites coding for peptide binding (32% for humpback and 40% for southern right) suggested that these loci were subject to positive (overdominant) selection. A population survey of humpback whales detected 23 alleles, differing by up to 21% of their inferred amino acid sequences. Amplification of the DRB exon 2 resulted in two groups of sequences. One was most similar to the DRB3 of the cow and present in all whales screened to date, including toothed whales. The second was most similar to the DRB2 of the cow and was found only in the bowhead and right whales. Both loci showed low diversity among species and apparent loss of function or altered function including interruption of reading frames. Finally, comparison of inferred protein sequence of the DRB3-like locus suggested convergence with the DQB, perhaps resulting from intergenic conversion or recombination.

Structure of the BoLA-DRB3 gene and promoter

The cattle major histocompatibility complex (MHC) class II DR gene product is a heterodimer encoded by the BoLA-DRA and -DRB3 genes. Several groups have isolated cDNA and genomic clones for these genes, but their full genomic organization has not been described. We used a combination of long-range polymerase chain reaction (PCR), cloning and sequencing to define the organization of the DRB3 gene on existing genomic clones and in genomic DNA. We estimate the size of the coding region to be 11.4 kbp. Sequencing of full-length PCR clones from two different haplotypes confirmed that they carried complete DRB3 genes and allowed the design of probes and primers to isolate and characterize the DRB3 promoter and 3' end. Fragments carrying the 5' end of the DRB3 gene and its promoter were identified on bacterial artificial chromosome (BAC) clones carrying the BoLA-DR genes. A 10-kbp promoter fragment was subcloned from one clone and a 1.7-kbp region including exon 1 and the promoter was sequenced. A 3-kbp fragment encoding exons 4-6 and the entire 3' untranslated region of the DRB3 gene was isolated from lambda clone A1 and sequenced. This provides us with improved characterization of the DRB3*0101 and DRB3*2002 alleles, and also subcloned 5' and 3' flanking regions of the polymorphic DRB3 gene for use in functional studies.

Genotyping BoLA-DRB3 alleles in Brazilian Dairy Gir cattle (Bos indicus) by temperature-gradient gel electrophoresis (TGGE) and direct sequencing

BoLA-DRB3 is a gene of the major histocompatibility complex (MHC) in cattle. The product of the BoLA-DRB3 gene is a beta chain of an MHC class II molecule, a glycoprotein expressed on the surface of antigen-presenting cells (APCs). Responses of CD4+ T lymphocytes to peptides are dependent on the presentation of peptide ligands bound to class II molecules on APCs. Genotyping of the BoLA-DRB3 gene is relatively complex due to the extensive polymorphism of this locus. Current techniques for assignment of genotypes are polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP), direct sequencing of PCR products, cloning-sequencing, polymerase chain reaction using sequence-specific oligonucleotide probes (PCR-SSOP), and denaturant-gradient gel electrophoresis. These techniques are time-consuming, do not discriminate all possible alleles, or are not readily reproducible. The objective of this study was to genotype BoLA-DRB3 using temperature-gradient gel electrophoresis (TGGE) to separate alleles before sequencing. PCRs using 28 DNA samples from Gir Dairy cattle (a Brazilian breed of Bos indicus) were submitted to TGGE. New PCR products were generated from separated alleles, purified, and sequenced. Allele separation was possible in 21 out of 26 heterozygote samples (81%). Results indicate that two sequence reads (forward and reverse) were sufficient for accurate genotyping of BoLA-DRB3 alleles. Separation of alleles by TGGE provides high-throughput, reliable typing of BoLA-DRB3, which is critical in disease association studies in cattle.

Pseudogenes: are they "junk" or functional DNA?

Pseudogenes have been defined as nonfunctional sequences of genomic DNA originally derived from functional genes. It is therefore assumed that all pseudogene mutations are selectively neutral and have equal probability to become fixed in the population. Rather, pseudogenes that have been suitably investigated often exhibit functional roles, such as gene expression, gene regulation, generation of genetic (antibody, antigenic, and other) diversity. Pseudogenes are involved in gene conversion or recombination with functional genes. Pseudogenes exhibit evolutionary conservation of gene sequence, reduced nucleotide variability, excess synonymous over nonsynonymous nucleotide polymorphism, and other features that are expected in genes or DNA sequences that have functional roles. We first review the Drosophila literature and then extend the discussion to the various functional features identified in the pseudogenes of other organisms. A pseudogene that has arisen by duplication or retroposition may, at first, not be subject to natural selection if the source gene remains functional. Mutant alleles that incorporate new functions may, nevertheless, be favored by natural selection and will have enhanced probability of becoming fixed in the population. We agree with the proposal that pseudogenes be considered as potogenes, i.e., DNA sequences with a potentiality for becoming new genes.

Short communication: characterization of DRB3 alleles in the MHC of Japanese shorthorn cattle by polymerase chain reaction-sequence-based typing

A study was made of exon 2 of the bovine leukocyte antigen BoLA-DRB3 gene of 176 Japanese Shorthorn cattle at six farms in Japan using polymerase chain reaction-sequence-based typing (PCR-SBT). An initial round of PCR using conserved locus-specific primers, a second round of PCR using a locus-specific primer, and at least one sequence-specific primer (SSP), followed by direct sequencing of products of PCR with SSP were conducted. Twenty-one BoLA-DRB3 alleles were identified with frequencies ranging from 0.3 to 19.6% in 176 individuals, and two of these alleles were new alleles that have not been reported previously. The three most frequently observed alleles (DRB3*1201, *0301, and *0801) accounted for 43.8% of the alleles in the population of these herds. Next, we tested the products of amplification by PCR of BoLA-DRB3 exon 2 with RsaI, BstYI, and HaeIII, and identified 18 previously described PCR-restriction fragment length polymorphism (RFLP) types. The PCR-RFLP types reflected the results of PCR-SBT exactly. Our results indicate that exon 2 of the BoLA-DRB3 gene is highly polymorphic in Japanese Shorthorn cattle.

Multi-primer target PCR for rapid identification of bovine DRB3 alleles

Multi-primer target polymerase chain reaction (MPT-PCR) is a rapid method for the identification of specific BoLA-DRB3 alleles. In a single PCR reaction, the presence of two alleles associated with increased risk, DRB3.2*23 (DRB3*2701-2703, 2705-2707) and decreased risk, DRB3.2*16 (DRB3*1501, 1502), of mastitis in Canadian Holstein can be detected. Two outer primers amplify exon 2 of DRB3. Simultaneously, two inner, allele-specific primers amplify individual alleles. Initially, 40 cows previously typed by PCR-restriction fragment length polymorphism (PCR-RFLP) were genotyped using the multi-primer approach. An additional 30 cows were first genotyped by multi-primer target PCR, then by PCR-RFLP. All animals were correctly identified and there were no false positives. This technique can readily be modified to identify other BoLA alleles of interest.

Characterization of major histocompatibility complex DRB diversity in the endemic South African antelope Damaliscus pygargus: a comparison in two subspecies with different demographic histories

Major histocompatibility complex (MHC) class II locus DRB was investigated by single-strand conformation polymorphism analysis (SSCP) and sequence analysis in the endemic South African antelope, Damaliscus pygargus, of which there are two subspecies. Greater polymorphism was found in the blesbok (D. p. phillipsi) subspecies (n = 44; 22 alleles) than in the bontebok (D. p. pygargus) subspecies (n = 45; 6 alleles). Erosion of allelic diversity in bontebok was most likely the result of two severe bottleneck events caused by hunting pressure and parasitic infection. A majority of the polymorphism observed was found within the peptide binding region (PBR) where dN/dS ratios were higher than for the non-PBR region. This, and the apparent trans-species relationship among alleles in a bovid phylogeny, suggest the evolution of diversity by heterosis or frequency-dependent selection.

Characterization of 18 new BoLA-DRB3 alleles

The second exon of the bovine MHC class II DRB3 gene was amplified by polymerase chain reaction (PCR) from DNA samples of 568 zebu Brahman cattle (Bos indicus) from Martinique (French West Indies). Cloning of these PCR products allowed the isolation of both alleles from each animal, which were characterized by the PCR-restriction fragment length polymorphism (RFLP) technique using the restriction enzymes RsaI, BstYI and HaeIII. Four new PCR-RFLP patterns were obtained by digestion with RsaI. These patterns were named 'v', 'w', 'x' and 'y' continuing the accepted nomenclature. Sequencing of each allele allowed the identification of 18 new BoLA-DRB3 exon 2 nucleotide sequences and their deduced amino acid sequences.

Cattle MHC: evolution in action?

Because major histocompatibility complex (MHC) genes play a major role in the development of acquired immune responses, it is essential to obtain comparative information on their organisation, expression and possible functional dichotomies in different species. In human, three classical, polymorphic class I genes (HLA-A, B- and -C) and four expressed A/B class II gene pairs (HLA-DM, -DP, -DQ and -DR) are each present on all haplotypes. With the exception of the HLA-DRB loci, it has been assumed that a similar rigid organisational situation exists in other mammalian species. However, extensive analysis of the bovine MHC (BoLA) at both the genomic and transcriptional levels has revealed a degree of genetic fluidity not described in other species. None of the four (or more) classical class I genes identified is consistently expressed, and haplotypes differ from one another in both the number and composition of expressed class I genes. Similarly, in the class II region, the number of DQ genes varies between haplotypes in both number and composition. These variations in both class I and II (which appear to reflect differences at the genomic level) are likely to play an important role in cattle immune responses. The observed phenotypic differences in cattle demonstrate very clearly the dynamic nature of the MHC region. This review addresses the functional impact of such variation in different breeds and populations, and its significance in terms of MHC evolution.

African buffalo maintain high genetic diversity in the major histocompatibility complex in spite of historically known population bottlenecks

Historical population collapses caused by rinderpest epidemics are hypothesized to have resulted in notable genetic losses in populations of the African buffalo. Polymorphism in the major histocompatibity complex (MHC) DRB3 gene was probed by means of restriction analysis of the sequence encoding the peptide-binding region. Nucleotide substitution patterns agreed with a positive selection acting on this fitness-relevant locus. Buffalo populations from four National Parks, situated in eastern and southern Africa, each revealed a surprisingly high allelic diversity. Current high levels of heterozygosity may be reconciled with historical bottlenecks by assuming that local extinctions were followed by fast recolonization, in accordance with the high dispersive capabilities of buffalo. The specific amplification of DRB3 alleles also enabled the assignment of individual genotypes. For each population sample a deficiency in the expected number of heterozygous animals was found. As overdominant selection on the MHC is predicted to yield an excess of heterozygous individuals, this may not be a locus-specific effect. Several other explanations are discussed, of which increased homozygosity caused by nonrandom mating of buffalo in populations seems the most probable.

A PCR-RFLP typing method for the caprine Mhc class II DRB gene

A polymerase chain reaction (PCR) technique associated with restriction fragment length polymorphism (RFLP) analysis has been developed for typing the second exon of the caprine Mhc class II DRB gene. By using a maximum of three restriction enzymes, corresponding to different combinations of RsaI, BsaI, AccI, NdeII, HaeIII and HpaII, this procedure allows to distinguish unequivocally 18 of the 22 caprine DRB alleles. Forty goats of several breeds were typed with this technique and two new RsaI restriction patterns which did not match with previously published sequence data were identified. Close associations have been found between RFLPs and amino acid substitutions at positions which are expected to be involved in the formation of the antigen-recognition site (ARS) of the DR molecule. These results suggest that PCR-RFLP may be a useful tool in typing the caprine DRB gene and in relating amino acid substitutions at the ARS of the DR molecule with disease resistance.

Demonstration of three DRB loci in a domestic horse family

Single-strand conformational polymorphism (SSCP) gel electrophoresis and DNA sequencing were used to characterize the second exon of the horse DRB homologue as well as to identify eight new DRB alleles. The SSCP gels presented a complex pattern, with phenotypes exhibiting between 4 and 13 bands. The DRB SSCP patterns were studied for two families (6 to 13 bands per pattern). For both families, the patterns showed simple Mendelian inheritance. The polymerase chain reaction products from two individuals possessing homozygous major histocompatibility complex (MHC) alleles by descent were cloned and retested on SSCP gels. All bands derived from the genomic DNA amplification could be accounted for with bands derived from the cloned DNA amplification products. The results were consistent with three DRB loci, though this number may be variable within the domestic horse population. Gene sequences were variable among the different products, and we were unable to assign locus designations for particular sequences. Amplification of cDNA library material derived from one of the individuals who is MHC homozygous by descent showed an SSCP profile suggesting that all three DRB loci are transcribed into mRNA.

Detection of a common BoLA-DRB3 deletion by sequence-specific oligonucleotide typing

The bovine MHC clas II allele BoLA-DRB3*2A has an amino acid deletion of unknown function at codon 65 in the second exon, which codes for the antigen-binding site. Sequence-specific oligonucleotides were designed based on published nucleotide sequences on BoLA-DRRB3 alleles, and used to detect this deletion in 51 Hereford cattle. Probes 65+ and 65- detect the presence or absence of codon 65 respectively. Oligonucleotide probes were labelled with Digoxigenin (DIG), hybridized to dot blots of BoLA-DRB3 exon 2 polymerase chain reaction (PCR) product, and detected by chemiluminescence. Of the 51 animals screened, two were homozygous and 11 were heterozygous for the deletion at codon 65. The methodology described here provides the necessary tools to screen rapidly for this deletion in a large number of animals in order to study its effect on antigen binding and immune response.

An amino acid sequence coded by the exon 2 of the BoLA DRB3 gene associated with a BoLA class I specificity constitutes a likely genetic marker of resistance to dermatophilosis in Brahman zebu cattle of Martinique (FWI)

One hundred and twenty-seven Brahman cattle from several locations in Martinique (FWI), reared under different environmental conditions, were followed over three years and checked for clinical signs of dermatophilosis. To confirm that these animals had been in contact with the pathogen Dermatophilus congolensis, their sera were tested by ELISA. On the basis of this epidemiological study, 12 animals were classified as resistant (seropositive without clinical signs), belonging to herds in which the prevalence of the disease ranged from 25 to nearly 98%. Eighteen animals classified as highly susceptible displayed severe characteristic skin lesions. These 30 selected animals were typed for class I antigens of the major histocompatibility complex (MHC). MHC class II genes were analyzed using the polymerase chain reaction (PCR) and restriction fragment length polymorphism (RFLP) techniques, on the exon 2 of the bovine leucocyte antigen (BoLA) DRB3 gene. Several alleles were found, according to patterns provided by the restriction enzymes used: Fnu 4HI, Dpn II, Hae III, and Rsa I. A particular sequence "EIAY" at amino acid positions 66/67/74/78 located in the antigen recognition sites (ARS) was found in the 12 animals classified as resistant, and 10 of them displayed also class I BoLA-A8 specificity. On the other hand, only 3 out of the 18 susceptible animals showed simultaneously the BoLA-DRB3 "EIAY" sequence and BoLA-A8 specificity. Interestingly, a serine residue at position 30 of the ARS was found in 8 of the susceptible animals and was completely absent from all resistant animals. Furthermore, in a same animal, the serine at position 30 and the EIAY sequence were never found simultaneously on the same haplotype. These results show a strong correlation between the resistant character to dermatophilosis and the association of MHC haplotypes: the BoLA-A8 specificity and the BoLA-DRB3 "EIAY" sequence at ARS positions 66/67/74/78 with the lack of serine in position 30. To confirm these results, family segregation studies are in progress and some interesting observations have been obtained.

Structure and organization of pig MHC class II DRB genes: evidence for genetic exchange between loci

The pig major histocompatibility complex DRB genes were studied by polymerase chain reaction (PCR) amplification of exon 2 from eight domestic pigs and two European wild boars. Sequence comparisons together with a phylogenetic analysis showed the existence of at least three DRB genes of which only one appears to be expressed. The two putative DRB pseudogenes contained deletions in exon 2, making it possible to confirm the presence of three non-allelic DRB genes by analyzing the length polymorphism of the amplified PCR products. The expressed gene shows allelic polymorphism at the same positions as in the human DRB1 gene. In addition, this pig gene shows extensive allelic polymorphism at positions 84-88, whereas, e.g., human DRB genes do not. Surprisingly, the two putative DRB pseudogenes also display a considerable amount of allelic polymorphism, albeit of a different character as compared with the expressed DRB gene. Short stretches of sequences are shared between individual alleles at different loci. These sequence similarities cannot be due to natural selection, since two of the three DRB genes involved are polymorphic pseudogenes constituting allelic series that have diverged after the inactivation event. Instead, the results indicate that the sequences have been exchanged between the DRB genes by intergenic recombination.

Identification and characterization of new BoLA-DRB3 alleles by heteroduplex analysis and direct sequencing

A sample of 52 mixed-breed dairy cattle (Holstein Friesian and Jersey) and 51 beef cattle (Hereford) from south-east Queensland was studied. The second exon of BoLA-DRB3 was amplified by polymerase chain reaction (PCR), and polymorphisms were detected by heteroduplex analysis. A large number of different heteroduplex patterns indicated extensive sequence polymorphism. Direct sequencing of PCR products from 17 homozygotes and cloning and sequencing of PCR product from two heterozygotes resulted in the identification and characterization of four novel alleles. The previously described allele BoLA-DRB3*2A is characterized by an amino acid deletion at position 65. We have identified three animals that are homozygous for this amino acid deletion, indicating that the deletion is unlikely to result in loss of function. Two of these animals had allele BoLA-DRB3*2A, and one had a novel allele with codon 65 deleted but differing from BoLA-DRB3*2A at a number of other amino acid positions. In conclusion, heteroduplex analysis allows rapid discrimination between homozygotes and heterozygotes, and enables rapid identification of new BoLA-DRB3 alleles.

Expression of multiple DQB genes in Bos indicus cattle

Polymorphism of expressed bovine major histocompatibility complex (MHC) class II DQB genes was investigated in a group of nine MHC-homozygous Kenya Boran cattle (Bos indicus). DQB second exon fragments were amplified by reverse transcription-polymerase chain reaction of total mononuclear cell RNA, cloned and sequenced. While a single DQB sequence was obtained from some animals, two DQB exon 2 sequences were found in others, implying expression of duplicated DQB genes. Two pairs of duplicated DQB genes were found in this group of homozygotes, and sequence analysis showed that both pairs contained distinct DQB alleles. One DQB duplication was observed in two related animals with the same MHC haplotype (F100 and F188), while the second was detected in three animals (F187, G166 and G277) that expressed the same class II DRB3 and DQB alleles but had different class I (BoLA-A) types, suggesting that this DR/DQ haplotype may be widespread in Boran cattle.

Functional effects of a natural polymorphism in the transcriptional regulatory sequence of HLA-DQB1

DNA sequence polymorphism in the genes encoding HLA class II proteins accounts for allelic diversity in antigen recognition and presentation and, thus, in the role of these cell surface glycoproteins as determinants of the scope of the T-cell repertoire. In addition, sequence polymorphism in the promoter-proximal transcriptional regulatory regions of these genes has been described, particularly for the HLA-DQB1 locus, where these differences may contribute to variation in locus- and allele-specific expression. In this study, we measured the effect of such regulatory sequence polymorphism on the expression of endogenous alleles of DQB1 in heterozygous cells. Quantitative reverse transcriptase-mediated PCR analysis showed that expression of the DQB1*0301 allele responded more rapidly to gamma interferon induction than that of DQB1*0302. We have analyzed functional effects of a prominent allelic polymorphism that consists of a TG dinucleotide present between the W and X1 consensus elements in the DQB1*0302 allele but missing in the DQB1*0301 allele. The dominant effect of this polymorphism was to introduce a variation in the spacing between the W and X1 elements of these two alleles. A secondary compensatory effect was specific for the TG dinucleotide itself, which was essential for the binding of a nuclear protein complex to the *0302 regulatory region immediately 5' of the X1 element. Derivatives of the DQB1 5' regulatory region were used to drive expression of the chloramphenicol acetyltransferase gene in transient transfections of human B-lymphoblastoid and gamma interferon-treated melanoma cell lines, demonstrating that the additional spacing between the W and X1 elements caused by the presence of the TG dinucleotide in the *0302 allele resulted in reduced expression compared with that driven by the *0301 fragment; this difference overshadowed an up-regulating effect on expression which corresponded to the binding of the TG-dependent nuclear protein complex. The presence of this polymorphism in multiple HLA-DQB1 alleles and in several species suggests selection for two alternative transcriptional regulatory mechanisms influencing expression of alleles of the same HLA locus.

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

The genetic diversity of the bovine class II DRB3 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 RsaI, BstYI, and HaeIII 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.

Low major histocompatibility complex class II diversity in European and North American moose

Major histocompatibility complex (MHC) genes encode cell surface proteins whose function is to bind and present intracellularly processed peptides to T lymphocytes of the immune system. Extensive MHC diversity has been documented in many species and is maintained by some form of balancing selection. We report here that both European and North American populations of moose (Alces alces) exhibit very low levels of genetic diversity at an expressed MHC class II DRB locus. The observed polymorphism was restricted to six amino acid substitutions, all in the peptide binding site, and four of these were shared between continents. The data imply that the moose have lost MHC diversity in a population bottleneck, prior to the divergence of the Old and New World subspecies. Sequence analysis of mtDNA showed that the two subspecies diverged at least 100,000 years ago. Thus, viable moose populations with very restricted MHC diversity have been maintained for a long period of time. Both positive selection for polymorphism and intraexonic recombination have contributed to the generation of MHC diversity after the putative bottleneck.