Class II genes of the human major histocompatibility complex. Organization and evolutionary relationship of the DR beta genes

Unraveling the architecture of major histocompatibility complex class II haplotypes in rhesus macaques

The regions in the genome that encode components of the immune system are often featured by polymorphism, copy number variation, and segmental duplications. There is a need to thoroughly characterize these complex regions to gain insight into the impact of genomic diversity on health and disease. Here we resolve the organization of complete major histocompatibility complex (MHC) class II regions in rhesus macaques by using a long-read sequencing strategy (Oxford Nanopore Technologies) in concert with adaptive sampling. In particular, the expansion and contraction of the primate DRB-region appear to be a dynamic process that involves the rearrangement of different cassettes of paralogous genes. These chromosomal recombination events are propagated by a conserved pseudogene, DRB6, which features the integration of two retroviral elements. In contrast, the DRA locus appears to be protected from rearrangements, which may be owing to the presence of an adjacently located truncated gene segment, DRB9 With our sequencing strategy, the annotation, evolutionary conservation, and potential function of pseudogenes can be reassessed, an aspect that was neglected by most genome studies in primates. Furthermore, our approach facilitates the characterization and refinement of an animal model essential to study human biology and disease.

Evaluation of HLA typing content of next-generation sequencing datasets from family trios and individuals of arab ethnicity

Introduction: HLA typing is a critical tool in both clinical and research applications at the individual and population levels. Benchmarking studies have indicated HLA-HD as the preferred tool for accurate and comprehensive HLA allele calling. The advent of next-generation sequencing (NGS) has revolutionized genetic analysis by providing high-throughput sequencing data. This study aims to evaluate, using the HLA-HD tool, the HLA typing content of whole exome, whole genome, and HLA-targeted panel sequence data from the consanguineous population of Arab ethnicity, which has been underrepresented in prior benchmarking studies. Methods: We utilized sequence data from family trios and individuals, sequenced on one or more of the whole exome, whole genome, and HLA-targeted panel sequencing technologies. The performance and resolution across various HLA genes were evaluated. We incorporated a comparative quality control analysis, assessing the results obtained from HLA-HD by comparing them with those from the HLA-Twin tool to authenticate the accuracy of the findings. Results: Our analysis found that alleles across 29 HLA loci can be successfully and consistently typed from NGS datasets. Clinical-grade whole exome sequencing datasets achieved the highest consistency rate at three-field resolution, followed by targeted HLA panel, research-grade whole exome, and whole genome datasets. Discussion: The study catalogues HLA typing consistency across NGS datasets for a large array of HLA genes and highlights assessments regarding the feasibility of utilizing available NGS datasets in HLA allele studies. These findings underscore the reliability of HLA-HD for HLA typing in underrepresented populations and demonstrate the utility of various NGS technologies in achieving accurate HLA allele calling.

Distributions of 11-loci HLA alleles typed by amplicon-based next-generation sequencing in South Koreans

The range of HLA typing for successful hematopoietic stem cell transplantation (HSCT) is gradually expanding with the next-generation sequencing (NGS)-based improvement in its quality. However, it is influenced by the allocation of finances and laboratory conditions. HLA-A, -B, -C, -DRB1/3/4/5, -DQA1, -DQB1, -DPA1, and -DPB1 alleles were genotyped at the 3-field level by amplicon-based NGS using MiSeqDx system and compared to our previous study employing long-range PCR and NGS using TruSight HLA v2 kit, in healthy donors from South Korea. Exon 2, exons 2/3, exons 2/3/4 or 5 of 11-loci were amplified by multiplex PCR. The sequence reads of over 53 depth counts were consistently obtained in each sample exon, depending on the target exon determined to match the reference sequence contained in the IPD-IMGT/HLA Database. HLA alleles were investigated by combinations of the determined exons. A total of 18 alleles with a frequency over 10% were found at the 11 HLA loci. Three ambiguities of HLA-A, -C, and -DRB1 were resolved. We observed a total of 26 HLA-A ~ C ~ B and 6 HLA-DRB1 ~ DQA1 ~ DQB1 ~ DPA1 ~ DPB1 haplotypes having significant linkage disequilibrium between alleles at all neighboring HLA loci. This result was compatible with the previous one, using TruSight HLA v2 kit. Advantages are simple and short progress time because one plate is used for each PCR step in one PCR machine and 11-loci HLA typing is possible even if only eight samples. These data suggested that expanded 11-loci HLA typing data by amplicon-based NGS might help perform HSCT.

Functional immune diversity in reindeer reveals a high Arctic population at risk

Climate changes the geographic range of both species as well as pathogens, causing a potential increase in the vulnerability of populations or species with limited genetic diversity. With advances in high throughput sequencing (HTS) technologies, we can now define functional expressed genetic diversity of wild species at a larger scale and identify populations at risk. Previous studies have used genomic DNA to define major histocompatibility complex (MHC) class II diversity in reindeer. Varying numbers of expressed genes found in many ungulates strongly argues for using cDNA in MHC typing strategies to ensure that diversity estimates relate to functional genes. We have used available reindeer genomes to identify candidate genes and established an HTS approach to define expressed MHC class I and class II diversity. To capture a broad diversity we included samples from wild reindeer from Southern Norway, semi-domesticated reindeer from Northern Norway and reindeer from the high Artic archipelago Svalbard. Our data show a medium MHC diversity in semi-domesticated and wild Norwegian mainland reindeer, and low MHC diversity reindeer in Svalbard reindeer. The low immune diversity in Svalbard reindeer provides a potential risk if the pathogenic pressure changes in response to altered environmental conditions due to climate change, or increased human-related activity.

HLA class II peptide-binding-region analysis reveals funneling of polymorphism in action

BACKGROUND

HLA-class II proteins hold important roles in key physiological processes. The purpose of this study was to compile all class II alleles reported in human population and investigate patterns in pocket variants and their combinations, focusing on the peptide-binding region (PBR).

METHODS

For this purpose, all protein sequences of DPA1, DQA1, DPB1, DQB1 and DRB1 were selected and filtered, in order to have full PBR sequences. Proportional representation was used for pocket variants while population data were also used.

RESULTS

All pocket variants and PBR sequences were retrieved and analyzed based on the preference of amino acids and their properties in all pocket positions. The observed number of pocket variants combinations was much lower than the possible inferred, suggesting that PBR formation is under strict funneling. Also, although class II proteins are very polymorphic, in the majority of the reported alleles in all populations, a significantly less polymorphic pocket core was found.

CONCLUSIONS

Pocket variability of five HLA class II proteins was studied revealing favorable properties of each protein. The actual PBR sequences of HLA class II proteins appear to be governed by restrictions that lead to the establishment of only a fraction of the possible combinations and the polymorphism recorded is the result of intense funneling based on function.

Allele and haplotype frequencies of human leukocyte antigen-A, -B, -C, -DRB1, -DRB3/4/5, -DQA1, -DQB1, -DPA1, and -DPB1 by next generation sequencing-based typing in Koreans in South Korea

Allele frequencies and haplotype frequencies of HLA-A, -B, -C, -DRB1, -DRB3/4/5, -DQA1, -DQB1, -DPA1, and -DPB1 have been rarely reported in South Koreans using unambiguous, phase-resolved next generation DNA sequencing. In this study, HLA typing of 11 loci in 173 healthy South Koreans were performed using next generation DNA sequencing with long-range PCR, TruSight^® HLA v2 kit, Illumina MiSeqDx platform system, and Assign^™ for TruSight^™ HLA software. Haplotype frequencies were calculated using the PyPop software. Direct counting methods were used to investigate the association with DRB1 for samples with only one copy of a particular secondary DRB locus. We compared these allele types with the ambiguous allele combinations of the IPD-IMGT/HLA database. We identified 20, 40, 26, 31, 19, 16, 4, and 16 alleles of HLA-A, HLA-B, HLA-C, HLA-DRB1, HLA-DQA1, HLA-DQB1, HLA-DPA1, and HLA-DPB1, respectively. The number of HLA-DRB3/4/5 alleles was 4, 5, and 3, respectively. The haplotype frequencies of most common haplotypes were as follows: A*33:03:01-B*44:03:01-C*14:03-DRB1*13:02:01-DQB1*06:04:01-DPB1*04:01:01 (2.89%), A*33:03:01-B*44:03:01-C*14:03 (4.91%), DRB1*08:03:02-DQA1*01:03:01-DQB1*06:01:01-DPA1*02:02:02-DPB1*05:01:01 (5.41%), DRB1*04:05:01-DRB4*01:03:01 (12.72%), DQA1*01:03:01-DQB1*06:01:01 (13.01%), and DPA1*02:02:02-DPB1*05:01:01 (30.83%). In samples with only one copy of a specific secondary DRB locus, we examined its association with DRB1. We, thus, resolved 10 allele ambiguities in HLA-B, -C (each exon 2+3), -DRB1, -DQB1, -DQA1, and -DPB1 (each exon 2) of the IPD-IMGT/HLA database. Korean population was geographically close to Japanese and Han Chinese populations in the genetic distances by multidimensional scaling (MDS) plots. The information obtained by HLA typing of the 11 extended loci by next generation sequencing may be useful for more exact diagnostic tests on various transplantations and the genetic population relationship studies in South Koreans.

Similar patterns of genetic diversity and linkage disequilibrium in Western chimpanzees (Pan troglodytes verus) and humans indicate highly conserved mechanisms of MHC molecular evolution

BACKGROUND

Many species are threatened with extinction as their population sizes decrease with changing environments or face novel pathogenic threats. A reduction of genetic diversity at major histocompatibility complex (MHC) genes may have dramatic effects on populations' survival, as these genes play a key role in adaptive immunity. This might be the case for chimpanzees, the MHC genes of which reveal signatures of an ancient selective sweep likely due to a viral epidemic that reduced their population size a few million years ago. To better assess how this past event affected MHC variation in chimpanzees compared to humans, we analysed several indexes of genetic diversity and linkage disequilibrium across seven MHC genes on four cohorts of chimpanzees and we compared them to those estimated at orthologous HLA genes in a large set of human populations.

RESULTS

Interestingly, the analyses uncovered similar patterns of both molecular diversity and linkage disequilibrium across the seven MHC genes in chimpanzees and humans. Indeed, in both species the greatest allelic richness and heterozygosity were found at loci A, B, C and DRB1, the greatest nucleotide diversity at loci DRB1, DQA1 and DQB1, and both significant global linkage disequilibrium and the greatest proportions of haplotypes in linkage disequilibrium were observed at pairs DQA1 ~ DQB1, DQA1 ~ DRB1, DQB1 ~ DRB1 and B ~ C. Our results also showed that, despite some differences among loci, the levels of genetic diversity and linkage disequilibrium observed in contemporary chimpanzees were globally similar to those estimated in small isolated human populations, in contrast to significant differences compared to large populations.

CONCLUSIONS

We conclude, first, that highly conserved mechanisms shaped the diversity of orthologous MHC genes in chimpanzees and humans. Furthermore, our findings support the hypothesis that an ancient demographic decline affecting the chimpanzee populations - like that ascribed to a viral epidemic - exerted a substantial effect on the molecular diversity of their MHC genes, albeit not more pronounced than that experienced by HLA genes in human populations that underwent rapid genetic drift during humans' peopling history. We thus propose a model where chimpanzees' MHC genes regenerated molecular variation through recombination/gene conversion and/or balancing selection after the selective sweep.

Characterizing alleles with large deletions using region specific extraction

Two novel HLA class II alleles, DRB4*03:01N and DQB1*03:276N, containing large deletions were identified during routine typing. Extraction of DNA encompassing the deletions was carried out with a panel of capture oligonucleotides followed by whole genome amplification. Next generation DNA sequencing was then used to characterize the sequences. DRB4*03:01N has a 16 kilobase pair deletion stretching upstream from intron 2 toward centromeric DRB8. DQB1*03:276N has two deletions separated by 844 nucleotides. The first deletion (3.7 kilobase pairs) is upstream of intron 1 and the second deletion removes 3.3 kilobase pairs further upstream towards centromeric DQA2.

Genomic structure of the horse major histocompatibility complex class II region resolved using PacBio long-read sequencing technology

The mammalian Major Histocompatibility Complex (MHC) region contains several gene families characterized by highly polymorphic loci with extensive nucleotide diversity, copy number variation of paralogous genes, and long repetitive sequences. This structural complexity has made it difficult to construct a reliable reference sequence of the horse MHC region. In this study, we used long-read single molecule, real-time (SMRT) sequencing technology from Pacific Biosciences (PacBio) to sequence eight Bacterial Artificial Chromosome (BAC) clones spanning the horse MHC class II region. The final assembly resulted in a 1,165,328 bp continuous gap free sequence with 35 manually curated genomic loci of which 23 were considered to be functional and 12 to be pseudogenes. In comparison to the MHC class II region in other mammals, the corresponding region in horse shows extraordinary copy number variation and different relative location and directionality of the Eqca-DRB, -DQA, -DQB and -DOB loci. This is the first long-read sequence assembly of the horse MHC class II region with rigorous manual gene annotation, and it will serve as an important resource for association studies of immune-mediated equine diseases and for evolutionary analysis of genetic diversity in this region.

A quick and robust MHC typing method for free-ranging and captive primate species

Gene products of the major histocompatibility complex (MHC) of human and non-human primates play a crucial role in adaptive immunity, and most of the relevant genes not only show a high degree of variability (polymorphism) but also copy number variation (CNV) is observed. Due to this diversity, MHC proteins influence the capability of individuals to cope with various pathogens. MHC and/or MHC-linked gene products such as odorant receptor genes are thought to influence mate choice and reproductive success. Therefore, MHC typing of wild and captive primate populations is considered to be useful in conservation biology, which is, however, often hampered by the need of invasive and time-consuming methods. All intact Mhc-DRB genes in primates appear to possess a complex and highly divergent microsatellite, DRB-STR. A panel of 154 pedigreed olive baboons (Papio anubis) was examined for their DRB content by DRB-STR analysis of genomic DNA. Using the same methodology on DNA of feces samples, DRB variability of a silvery gibbon population (Hylobates moloch) (N = 24), an endangered species, could successfully be studied. In both species, length determination of the DRB-STR resulted in the definition of unique genotyping patterns that appeared to be specific for a certain chromosome. Moreover, the different STR lengths were shown to segregate with the allelic variation of the respective gene. The results obtained expand data gained previously on DRB-STR typing in macaques, great apes, and humans and strengthen the conclusion that this protocol is applicable in molecular ecology, conservation biology, and colony management, especially of endangered primate species.

Cost-efficient multiplex PCR for routine genotyping of up to nine classical HLA loci in a single analytical run of multiple samples by next generation sequencing

Background

HLA genotyping by next generation sequencing (NGS) requires three basic steps, PCR, NGS, and allele assignment. Compared to the conventional methods, such as PCR-sequence specific oligonucleotide primers (SSOP) and -sequence based typing (SBT), PCR-NGS is extremely labor intensive and time consuming. In order to simplify and accelerate the NGS-based HLA genotyping method for multiple DNA samples, we developed and evaluated four multiplex PCR methods for genotyping up to nine classical HLA loci including HLA-A, HLA-B, HLA-C, HLA-DRB1/3/4/5, HLA-DQB1, and HLA-DPB1.

Results

We developed multiplex PCR methods using newly and previously designed middle ranged PCR primer sets for genotyping different combinations of HLA loci, (1) HLA-DRB1/3/4/5, (2) HLA-DQB1 (3.8 kb to 5.3 kb), (3) HLA-A, HLA-B, HLA-C, and (4) HLA-DPB1 (4.6 kb to 7.2 kb). The primer sets were designed to genotype polymorphic exons to the field 3 level or 6-digit typing. When we evaluated the PCR method for genotyping all nine HLA loci (9LOCI) using 46 Japanese reference subjects who represented a distribution of more than 99.5% of the HLA alleles at each of the nine HLA loci, all of the 276 alleles genotyped, except for HLA-DRB3/4/5 alleles, were consistent with known alleles assigned by the conventional methods together with relevant locus balance and no excessive allelic imbalance. One multiplex PCR method (9LOCI) was able to provide precise genotyping data even when only 1 ng of genomic DNA was used for the PCR as a sample template.

Conclusions

In this study, we have demonstrated that the multiplex PCR approach for NGS-based HLA genotyping could serve as an alternative routine HLA genotyping method, possibly replacing the conventional methods by providing an accelerated yet robust amplification step. The method also could provide significant merits for clinical applications with its ability to amplify lower quantity of samples and the cost-saving factors.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1514-4) contains supplementary material, which is available to authorized users.

Characterising a microsatellite for DRB typing in Aotus vociferans and Aotus nancymaae (Platyrrhini)

Non-human primates belonging to the Aotus genus have been shown to be excellent experimental models for evaluating drugs and vaccine candidates against malaria and other human diseases. The immune system of this animal model must be characterised to assess whether the results obtained here can be extrapolated to humans. Class I and II major histocompatibility complex (MHC) proteins are amongst the most important molecules involved in response to pathogens; in spite of this, the techniques available for genotyping these molecules are usually expensive and/or time-consuming. Previous studies have reported MHC-DRB class II gene typing by microsatellite in Old World primates and humans, showing that such technique provides a fast, reliable and effective alternative to the commonly used ones. Based on this information, a microsatellite present in MHC-DRB intron 2 and its evolutionary patterns were identified in two Aotus species (A. vociferans and A. nancymaae), as well as its potential for genotyping class II MHC-DRB in these primates.

HLA-DRB1, -DRB3, -DRB4 and -DRB5 genotyping at a super-high resolution level by long range PCR and high-throughput sequencing

Super high-resolution single molecule sequence-based typing (SS-SBT) is a human leukocyte antigen (HLA) DNA typing method to the field 4 level of allelic resolution (formerly known as eight-digit typing) to efficiently detect new and null alleles without phase ambiguity by combination of long ranged polymerase chain reaction (PCR) amplification and next-generation sequencing (NGS) technologies. We previously reported the development and application of the SS-SBT method for the eight classical HLA loci, A, B, C, DRB1, DQA1, DQB1, DPA1 and DPB1. In this article, we describe the development of the SS-SBT method for three DRB1 linked loci, DRB3, DRB4 and DRB5 (DRB3/4/5) and characterization of DRB1-DRB3/4/5 haplotype structures to the field 4 level. Locus specific PCR primers for DRB3/4/5 were designed to amplify the gene regions from intron 1 to exon 6 [3' untranslated region (3'UTR)]. In total 20 DRB1 and 13 DRB3/4/5 allele sequences were determined by the SS-SBT to the field 4 level without phase ambiguity using 19 DR51, DR52 and DR53 positive genomic DNA samples obtained from Japanese. Moreover, 18 DRB1-DRB3/4/5 haplotypes were estimated to the field 4 level by the SS-SBT method in contrast to 10 haplotypes estimated by conventional methods to the field 1 level (formerly known as two digit typing). Therefore, DRB1-DRB3/4/5 haplotyping by SS-SBT is expected to provide informative data for improved HLA matching in medical research, transplantation procedures, HLA-related disease studies and human population diversity studies.

Evolution of HLA-DRB genes

The HLA region shows diversity concerning the number and content of DRB genes present per haplotype. Similar observations are made for the equivalent regions in other primate species. To elucidate the evolutionary history of the various HLA-DRB genes, a large panel of intron sequences obtained from humans, chimpanzees, rhesus macaques, and common marmosets has been subjected to phylogenetic analyses. Special attention was paid to the presence and absence of particular transposable elements and/or to their segments. The sharing of different parts of the same long interspersed nuclear element-2 (LINE2, L2) and various Alu insertions by the species studied demonstrates that one precursor gene must have been duplicated several times before the Old World monkey (OWM) and hominid (HOM) divergence. At least four ancestral DRB gene families appear to have been present before the radiation of OWM and HOM, and one of these even predates the speciation of Old and New World primates. Two of these families represent the pseudogenes DRB6/DRB2 and DRB7, which have been locked in the genomes of various primate species over long evolutionary time spans. Furthermore, all phylogenies of different intron segments show consistently that, apart from the pseudogenes, only DRB5 genes are shared by OWM and HOM, and they demonstrate the common history of certain DRB genes/lineages of humans and chimpanzees. In contrast, the evolutionary history of some other DRB loci is difficult to decipher, thus illustrating the complex history of the evolution of DRB genes due to a combination of mutations and recombination-like events. The selected approach allowed us to shed light on the ancestral DRB gene pool in primates and on the evolutionary relationship of the various HLA-DRB genes.

A multicentre study of acute kidney injury in severe sepsis and septic shock: association with inflammatory phenotype and HLA genotype

BACKGROUND

To investigate the association between severity of acute kidney injury (AKI) and outcome, systemic inflammatory phenotype and HLA genotype in severe sepsis.

METHODOLOGY/PRINCIPAL FINDINGS

Prospective multicenter observational study done in 4 intensive care units in two university hospitals. Severe sepsis and septic shock patients with at least 2 organ failures based on the SOFA score were classified: 1) "no AKI", 2) "mild AKI" (grouping stage 1 and 2 of AKIN score) and 3) "severe AKI" (stage 3 of AKIN score). Sequential measurements: The vasopressor dependency index (VDI; dose and types of drugs) to evaluate the association between hemodynamic status and the development of early AKI; plasma levels of IL-10, macrophage migration inhibitory factor (MIF), IL-6 and HLA-DR monocyte expression. Genotyping of the 13 HLA-DRB1 alleles with deduction of presence of HLA-DRB3, -DRB4 and -DRB5 genes. We used multivariate analysis with competitive risk model to study associations. Overall, 176 study patients (146 with septic shock) were classified from AKIN score as "no AKI" (n = 43), "mild AKI" (n = 74) or "severe AKI" (n = 59). The VDI did not differ between groups of AKI. After adjustment, "mild and severe AKI" were an independent risk factor for mortality (HR 2.42 95%CI[1.01-5.83], p = 0.048 and HR 1.99 95%CI[1.30-3.03], p = 0.001 respectively). "Severe AKI" had higher levels of plasma IL-10, MIF and IL-6 compared to "no AKI" and mild AKI (p<0.05 for each), with no difference in mHLA-DR at day 0. HLA-DRB genotyping showed a significantly lower proportion of 4 HLA-DRB alleles among patients requiring renal replacement therapy (RRT) (58%) than in patients with severe AKI who did not receive RRT (84%) (p = 0.004).

CONCLUSIONS

AKI severity is independently associated with mortality and plasma IL-10, MIF or IL-6 levels. Presence of 4 alleles of HLA-DRB in severe AKI patients seems associated with a lower need of RRT.

CCR5 blockade is well tolerated and induces changes in the tissue distribution of CCR5+ and CD25+ T cells in healthy, SIV-uninfected rhesus macaques

BACKGROUND

CCR5 is a main co-receptor for HIV, but also homes lymphocytes to sites of inflammation. We hypothesized that inhibition of CCR5 signaling would reduce HIV-associated chronic immune activation.

METHODS

To test this hypothesis, we administered an antagonistic anti-CCR5 monoclonal antibody (HGS101) to five uninfected rhesus macaques (RMs) and monitored lymphocyte dynamics in blood and tissue.

RESULTS

CCR5 blockade resulted in decreased levels of CCR5+ T cells in blood and, at later timepoints, in lymph nodes. Additionally, the levels of CD25+ T cells increased in lymph nodes, but decreased in blood, bone marrow, and rectal mucosa. Finally, a profile of gene expression from HGS101-treated RMs revealed a subtle, but consistent, in vivo signature of CCR5 blockade that suggests a mild immune-modulatory effect.

CONCLUSIONS

Treatment with anti-CCR5 antibody induces changes in the tissue distribution of CCR5+ and CD25+ T cells that may impact on the overall levels of immune activation during HIV and SIV infection.

The chimpanzee Mhc-DRB region revisited: gene content, polymorphism, pseudogenes, and transcripts

In humans, great apes, and different monkey species, the major histocompatibility complex (MHC) class II DRB region is known to display considerable copy number variation. The microsatellite D6S2878 has been shown to be a valuable marker for haplotyping the DR region in humans and macaque species. The present report illustrates that chimpanzee haplotypes also can be discriminated with this marker. The analyses resulted in the description of nine different region configurations, of which seven are present within the West African chimpanzee population studied. The region configurations vary in gene content from two up to five DRB genes. Subsequent cDNA sequencing increased the number of known full-length Patr-DRB sequences from 3 to 32, and shows that one to three Patr-DRB genes per haplotype apparently produce functional transcripts. This is more or less comparable to humans and rhesus macaques. Moreover, microsatellite analysis in concert with full-length DRB gene sequencing showed that the Patr-DRB*W9 and -DRB3*01/02 lineages most likely arose from a common ancestral lineage: hence, the Patr-DRB*W9 lineage was renamed to Patr-DRB3*07. Overall, the data demonstrate that the D6S2878 microsatellite marker allows fast and accurate haplotyping of the Patr-DRB region. In addition, the limited amount of allelic variation observed at the various Patr-DRB genes is in agreement with the fact that chimpanzees experienced a selective sweep that may have been caused by an ancient retroviral infection.

High resolution definition of HLA-DRB haplotypes by a simplified microsatellite typing technique

In humans, the region configurations DR1, DR8, DR51, DR52 and DR53 are known to display copy number as well as allelic variation, rendering high resolution typing of HLA-DRB haplotypes cumbersome. Advantage was taken of microsatellite D6S2878, present in all DRB genes/pseudogenes with an intact exon 2-intron 2 segment. This DRB-STR is highly polymorphic in composition and length. Recently, it was proven that all exon 2 sequences could be linked to a certain DRB-STR that segregates with the respective DRB allele. Because haplotypes show differential copy numbers and compositions of exon 2-positive DRB genes/pseudogenes, unique DRB-STR patterns could be described that appear to be specific for a particular DRB haplotype. The aim of this workshop project was to approve and to qualify this simple typing protocol in a larger panel covering different European populations. All participants succeeded in correctly defining the DRB-STR amplicons varying from 135 to 222 base pair (bp) lengths. The panel of 101 samples covered 50 DRB alleles distributed over 37 different haplotypes as defined by exon 2 sequence-based typing. These haplotypes could be refined into 105 haplotypes by DRB-STR typing. Thus, discrimination of exon 2-identical DRB alleles was feasible, as well as the exact description of three different crossing-over events that resulted in the generation of hybrid DR region configurations. This typing procedure appears to be a quick and highly robust technique that can easily be performed by different laboratories, even without experience in microsatellite typing; thus, it is suitable for a variety of researchers in diverse research areas.

Restriction fragment length polymorphism of DQ and DR class II genes of the bovine major histocompatibility complex

DQ alpha, DQ beta, DR alpha and DR beta class II genes of the bovine major histocompatibility complex (MHC) were investigated by Southern blot hybridizations using human probes. Hybridizations of these probes to genomic DNA, digested with PvuII or TaqI, revealed extensive restriction fragment length polymorphisms (RFLPs). The polymorphisms were interpreted genetically by analysing a family material, comprising five sires, 48 dams and 50 offspring, and a population sample comprising 197 breeding bulls. The analysis resolved 20 DQ alpha, 17 DQ beta, 5 DR alpha and 25 DR beta RFLP types. The segregation data were consistent with simple Mendelian inheritance of the RFLPs. The analysis of the bull sample showed that it is possible to apply the RFLP method for routine typing of class II polymorphism in population samples. The linkage disequilibrium in the DQ-DR region was found to be extremely strong as only about 20 DQ and about 30 DQ-DR haplotypes were observed despite the large number of possible haplotypes. Close linkage to the blood group locus M was also found; the M' allele occurred in strong linkage disequilibrium with the class II haplotype DQ1BDR alpha 4DR beta 1B. A population genetic analysis of the DQ data in the sample of breeding bulls revealed that the frequency of homozygotes was significantly lower than Hardy-Weinberg expectation and that the allele frequency distribution deviated significantly from the one expected for selectively neutral alleles.

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.

Comparative genetics of a highly divergent DRB microsatellite in different macaque species

The DRB region of the major histocompatibility complex (MHC) of cynomolgus and rhesus macaques is highly plastic, and extensive copy number variation together with allelic polymorphism makes it a challenging enterprise to design a typing protocol. All intact DRB genes in cynomolgus monkeys (Mafa) appear to possess a compound microsatellite, DRB-STR, in intron 2, which displays extensive length polymorphism. Therefore, this STR was studied in a large panel of animals, comprising pedigreed families as well. Sequencing analysis resulted in the detection of 60 Mafa-DRB exon 2 sequences that were unambiguously linked to the corresponding microsatellite. Its length is often allele specific and follows Mendelian segregation. In cynomolgus and rhesus macaques, the nucleotide composition of the DRB-STR is in concordance with the phylogeny of exon 2 sequences. As in humans and rhesus monkeys, this protocol detects specific combinations of different DRB-STR lengths that are unique for each haplotype. In the present panel, 22 Mafa-DRB region configurations could be defined, which exceeds the number detected in a comparable cohort of Indian rhesus macaques. The results suggest that, in cynomolgus monkeys, even more frequently than in rhesus macaques, new haplotypes are generated by recombination-like events. Although both macaque species are known to share several identical DRB exon 2 sequences, the lengths of the corresponding microsatellites often differ. Thus, this method allows not only fast and accurate DRB haplotyping but may also permit discrimination between highly related macaque species.

Impact of endogenous intronic retroviruses on major histocompatibility complex class II diversity and stability

The major histocompatibility complex (MHC) represents a multigene family that is known to display allelic and gene copy number variations. Primate species such as humans, chimpanzees (Pan troglodytes), and rhesus macaques (Macaca mulatta) show DRB region configuration polymorphism at the population level, meaning that the number and content of DRB loci may vary per haplotype. Introns of primate DRB alleles differ significantly in length due to insertions of transposable elements as long endogenous retrovirus (ERV) and human ERV (HERV) sequences in the DRB2, DRB6, and DRB7 pseudogenes. Although the integration of intronic HERVs resulted sooner or later in the inactivation of the targeted genes, the fixation of these endogenous retroviral segments over long time spans seems to have provided evolutionary advantage. Intronic HERVs may have integrated in a sense or an antisense manner. On the one hand, antisense-oriented retroelements such as HERV-K14I, observed in intron 2 of the DRB7 genes in humans and chimpanzees, seem to promote stability, as configurations/alleles containing these hits have experienced strong conservative selection during primate evolution. On the other hand, the HERVK3I present in intron 1 of all DRB2 and/or DRB6 alleles tested so far integrated in a sense orientation. The data suggest that multigenic regions in particular may benefit from sense introgressions by HERVs, as these elements seem to promote and maintain the generation of diversity, whereas these types of integrations may be lethal in monogenic systems, since they are known to influence transcript regulation negatively.

A highly divergent microsatellite facilitating fast and accurate DRB haplotyping in humans and rhesus macaques

The DRB region of the MHC in primate species is known to display abundant region configuration polymorphism with regard to the number and content of genes present per haplotype. Furthermore, depending on the species studied, the different DRB genes themselves may display varying degrees of allelic polymorphism. Because of this combination of diversity (differential gene number) and polymorphism (allelic variation), molecular typing methods for the primate DRB region are cumbersome. All intact DRB genes present in humans and rhesus macaques appear to possess, however, a complex and highly divergent microsatellite. Microsatellite analysis of a sizeable panel of outbred rhesus macaques, covering most of the known Mamu-DRB haplotypes, resulted in the definition of unique genotyping patterns that appear to be specific for a given haplotype. Subsequent examination of a representative panel of human cells illustrated that this approach also facilitates high-resolution HLA-DRB typing in an easy, quick, and reproducible fashion. The genetic composition of this complex microsatellite is shown to be in concordance with the phylogenetic relationships of various HLA-DRB and Mamu-DRB exon 2 gene/lineage sequences. Moreover, its length variability segregates with allelic variation of the respective gene. This simple protocol may find application in a variety of research avenues such as transplantation biology, disease association studies, molecular ecology, paternity testing, and forensic medicine.

The Relationship Between HLA Class II Polymorphisms and Somatic Deletions in Testicular B Cell Lymphomas of Dutch Patients

Several risk factors including immune deficiencies, infections, and autoimmune diseases have been established for non-Hodgkin's lymphoma (NHL). For diffuse large B cell lymphoma (DLBCL), the most common type of lymphoma, no risk factors have been described, which may be due to the intrinsic heterogeneity of this disorder. Previously we reported that, in contrast to nodal DLBCLs, the majority of testicular DLBCLs manifested complete loss of HLA-DR and -DQ expression associated with homozygous deletions of the corresponding genes. To determine the correlation between HLA class II polymorphisms and these lymphomas, we applied DNA typing for HLA-DRB1 and HLA-DQB1 on 50 Dutch patients with testicular and 48 with nodal DLBCL and compared the frequencies with a cohort of healthy Dutch controls. Both the patients with nodal and those with testicular DLBCL manifested significantly higher frequencies of HLA-DRB1*15 than the controls (p < 0.018, odds ratio 2.09 and p < 0.013, odds ratio 2.12, respectively). Moreover, a positive association was seen with HLA-DRB1*12 (p = 0.043, odds ratio 4.17) in the patients with testicular DLBCL, and a negative association was seen with HLA-DRB1*07 (p = 0.022, odds ratio 0.13) in the patients with nodal DLBCL. Homozygous deletions of the HLA-DR/DQ region, evaluated by interphase fluorescence in situ hybridization were seen in 20 of 48 testicular tumors. No preferential loss or retention of a particular HLA-DR or -DQ allele was seen because all alleles were at least once retained or involved in a homozygous deletion.

New strategies for efficient typing of HLA class-II loci DQB1 and DRB1 by using Pyrosequencing

The characterization of genetic risk factors for complex diseases located on chromosome-6 frequently requires human leucocyte antigen (HLA) genotyping of large patient cohorts. Currently available methods do not support high-throughput HLA typing beyond the major allele group level. We, thus, developed a high-throughput approach for the HLA-DQB1 and HLA-DRB1 loci that is based on Pyrosequencing. Pyrosequencing offers a higher degree of automation than direct sequencing or oligotyping. Using a dispensation order optimized for the particular HLA locus, rapid group typing and fine resolution can be achieved. We implemented the method for two important HLA loci--DQB1 and DRB1. The HLA-DQB1 typing method comprises the following steps: splitting the potential alleles after a generic polymerase chain reaction (PCR) amplification into groups with a first Pyrosequencing reaction and resolving the split allele groups by means of five further Pyrosequencing reactions. The HLA-DR gene family is known to be the most polymorphic one in the HLA class-II region because of a large number of DRB1 alleles. Because of this complex nature, HLA-DRB1 typing was performed by means of a combination of sequence-specific PCR typing and Pyrosequencing. HLA-DQB1 typing and HLA-DRB1 typing were performed successfully by using standard DNA samples with the help of known HLA genotypes and in a blind study by using the samples from the Deutscher Zell Austausch 2002 and 2003. The approach was optimized and was practically tested for genotyping in disease association studies. Our well-elaborated Pyrosequencing-based protocols offer a new alternative to the existing HLA class-II typing methods and represent a convenient and economic solution, a unique combination of high accuracy with high-sample throughput.

Estimation of linkage disequilibrium for loci with multiple alleles: basic approach and an application using data from bighorn sheep

The great expansion of population genetic data using molecular techniques now allows examination of the extent of linkage disequilibrium for many pairs of loci, each locus often with multiple alleles. The expectation-maximization (EM) algorithm for generating maximum likelihood estimates of gametic frequencies from multiallelic genotypic data is described and applied. The EM algorithm is used in desert bighorn sheep where the population size, and consequently the sample size, is often small. We calculated haplotype frequencies for all pairwise combinations of five major histocompatibility loci and three microsatellite loci in 14 populations; the performance of the algorithm is discussed. Disequilibrium values are calculated and tested for statistical significance. High levels of disequilibrium are found between all pairs of major histocompatibility complex (MHC) loci and between MHC and a linked microsatellite locus.

New DR5 sequences: a novel DRB1*11122 allele identified in Paiwan tribe members of Taiwan and a corrected sequence for the DRB1*1201 allele

We report herein the identification of a new DRB1 allele using sequence-based typing (SBT). This novel allele, HLA-DRB1*11122, was found in an aboriginal individual (SWP71) from the Paiwan tribe in the southern part of Taiwan. This individual was typed by SBT method as having an HLA genotype of HLA-A*24021/24021, HLA-B*4001/4002, HLA-DRB1*11122/15011, HLA-DRB3*0202, and HLA-DRB5*01011. This new allele differs from DRB1*1112 in the polymorphic exon 2 only at codon 34 (CAA-->CAG; both specify glutamine) and from DRB1*1110 in the exon 2 sequence only at codon 32 (CAT-->TAT; H32T). The most likely candidate allele which is found in the aboriginal populations of Taiwan and which may mutate into this new allele is DRB1*11011. DRB1*11122 allele differs from DRB1*11011 allele in the polymorphic exon 2 at both codon 34 (CAA-->CAG) and codon 37 (TAC-->TTC; T37F). This novel HLA-DRB1*11122 allele was also found in another aboriginal individual (SWP90) from the same Paiwan tribe. This SWP90 individual was typed by SBT method as having an HLA genotype of HLA-A*24021/24021, HLA-B*4002/5502, HLA-DRB1*11122/1201, and HLA-DRB3*01011/0202. However, the original DRB1*1201 sequence from HERLUFF was found to be erroneously reported and the corrected sequence from SWP90 is now presented herein.

High concentration of soluble HLA-DR in the synovial fluid: generation and significance in "rheumatoid-like" inflammatory joint diseases

In the search for its role in inflammatory joint diseases, soluble HLA-DR (sHLA-DR) was quantitated in 72 synovial fluids (SF) by a newly established immunoenzyme assay. Unlike other soluble receptors which accumulated only moderately (sCD25, sCD4) or negligibly (sHLA class I, sCD8) in the SF, SF sHLA-DR levels exceeded serum levels by up to 3 orders of magnitude and varied disease dependently from "control" values (traumatic synovitis and osteoarthritis: 9.9 +/- 6.1 ng/ml). Clear-cut different SF sHLA-DR values in HLA-DR-associated "rheumatoid-like" (136.5 +/- 130.0 ng/ml) vs HLA-B27-associated "spondylarthropathy-like" arthritic forms (28.4 +/- 29.1 ng/ml) were most significant comparing oligoarticular juvenile chronic arthritis type I (147.6 +/- 112.6 ng/ml) and type II (3.3 +/- 1.1 ng/ml), thus offering a new classification marker. Also ex vivo, large amounts of sHLA-DR were released spontaneously by SF mononuclear cells and found to be related to the T-cell activation state. SF sHLA-DR may be shed in large complexes or micelles, as it eluted mainly at >450 kDa on gel filtration. Western blotting revealed that the majority of SF sHLA-DR consisted of full-length alpha- and beta-chains. Minor fractions of smaller sized antigens seemed to be generated by proteolytic cleavage rather than by alternative splicing, since only minute amounts of HLA-DRB mRNA lacking the transmembrane exon could be amplified by RT-PCR. Distinct forms of high-dose sHLA-DR, able to provoke rather than to suppress T-cell responses, are discussed as contributing to some HLA-DR disease association.

The noncoding regions of HLA-DRB uncover interlineage recombinations as a mechanism of HLA diversification

The mechanisms generating new alleles at the MHC loci are still unknown in detail, and several proposals have been made to explain the extent of polymorphism. The patchwork pattern of polymorphism in the 2nd exon of HLA-DRB1 recommends this locus as a model for the study of the potential of interallelic gene conversion. In general, the inference of gene conversion-like events based exclusively on exon sequence comparisons may be misleading because the identity of the putative donor allele remains unknown. In this study, we describe five alleles of the HLA-DRB1 gene, which intron regions give evidence for interlineage recombination events either strictly located at the 2nd exon or involving the adjacent introns. Furthermore, we show that the noncoding regions provide important clues to the mechanisms of the generation of new alleles, and our results indicate that interlineage recombinations may be hidden and are perhaps more frequent than currently expected.

The insulator protein CTCF represses transcription on binding to the (gt)(22)(ga)(15) microsatellite in intron 2 of the HLA-DRB1(*)0401 gene

The insulator and transcription factor CTCF is a highly conserved 11 zinc finger protein possessing multiple specifities in DNA sequence recognition. CTCF regulates transcription of several genes, like the human oncogene c-myc or the chicken lysozyme gene by binding extremely divergent DNA sequences with different sets of its 11 zinc fingers. Recently, an insulator function was demonstrated for several CTCF binding elements. Here we show that CTCF binds to the (gt)(22)(ga)(15) microsatellite repeat A9 in intron 2 of the HLA-DRB1(*)0401 gene. Reporter gene activity is repressed by the A9 element. This repression is dependent on coexpressed CTCF and is even stronger compared with the CTCF binding site F1 of the chicken lysozyme gene, for which a silencer activity has been shown. This is the first report suggesting a function for microsatellite sequences in regulating specific gene expression.

Alu elements support independent origin of prosimian, platyrrhine, and catarrhine Mhc-DRB genes

The primate major histocompatibility complex (Mhc) genes fall into two classes and each of the classes into several families. Of the class II families, the DRB family has a long and complex evolutionary history marked by gene turnover, rearrangement, and molecular convergence. Because the history is not easily decipherable from sequences alone, Alu element insertions were used as cladistic markers to support the surmised phylogenetic relationships among the DRB genes. Intron 1 segments of 24 DRB genes from five platyrrhine species and five DRB genes from three prosimian species were amplified by PCR and cloned, and the amplification products were sequenced or PCR-typed for Alu repeats. Three Alu elements were identified in the platyrrhine and four in the prosimian DRB genes. One of the platyrrhine elements (Alu50J) is also found in the Catarrhini, whereas the other two (Alu62Sc, Alu63Sc) are restricted to the New World monkeys. Similarly, the four prosimian elements are found only in this taxon. This distribution of Alu elements is consistent with the phylogeny of the DRB genes as determined from their intron 1 sequences in an earlier and the present study. It contradicts the exon 2-based phylogeny and thus corroborates the conclusion that the evolution of DRB exon 2 sequences is, to some extent, shaped by molecular convergence. Taken together, the data indicate that each of the assemblages of DRB genes in prosimians, platyrrhines, and catarrhines is derived from a separate ancestral gene.

Phylogenetic history of hominoid DRB loci and alleles inferred from intron sequences

The evolutionary relationships among the MHC class II DRB4, DRB5 and DRB6 loci as well as the allelic lineages and alleles of the DRB1 locus were studied based on intron 1 and intron 2 sequences from humans, chimpanzee (Pan troglodytes), bonobo (Pan paniscus) and gorilla (Gorilla gorilla). The phylogenetic trees for these sequences indicate that most of the DRB1 allelic lineages predate the separation of the hominoid species studied, consistent with previous analysis of the coding sequences of these lineages. However, the intron sequence variation among alleles within DRB1 allelic lineages is very limited, consistent with the notion that the majority of the contemporary alleles have been generated within the last 250,000 years. The clustering of the DRB1 allelic lineages *08 and *12 with *03 supports a common ancestry for the DR8 and DR52 haplotypes. Similarly, the clustering of DRB1 allelic lineages *15 and *01 with the DRB3 locus is consistent with a common ancestry for the DR1 and DR51 haplotypes. Two cases of recombination around the second exon were observed: 1) the HLA-DRB6 locus appears to have been generated through a recombination between a DRB5 allele and an ancestral DRB6 allele, and 2) the gorilla sequence Gogo-DRB1 *03 appears to have been generated through a recombination between the DRB3 locus and an allele from the DRB1 *03 allelic lineage. The nucleotide substitution rate of DRB introns was estimated to 0.85-1.63 x 10(-9) per site per year, based on comparisons between the most closely related sequences from different hominoid species. This estimate is similar to the substitution rate for other intronic regions of the primate genome.

Molecular clock and recombination in primate Mhc genes

To set an accurate chronological framework to the evolution of primate class I and II genes in the major histocompatibility complex (Mhc), the rate of silent nucleotide substitutions in exons and introns is examined for various cDNA and genome sequences currently available. The rate is sensitive to the GC content and correlates negatively with increased GC biases at the third codon positions of Mhc genes. The intergenic recombination rate in the HLA region is estimated from the synonymous nucleotide differences at 37 linked loci. Any HLA subregion is recombined more or less at the ordinary rate of 1 cM per 1 Mb, although the rate may be reduced in some subregions. This information is used to discuss HLA haplotypes when they are applied to studies of human demography. The unusual polymorphism in the alpha-helix of HLA-DRB1 is also revisited in relation to intragenic recombination, but the molecular mechanism and the evolutionary cause both remain enigmatic.

Retroelements in the human MHC class II region

Molecular genetic studies of the human major histocompatibility complex (MHC) have led to the identification of more than 200 genes. Besides the large number of genes in the MHC, densely clustered areas of retroelements have been identified. These include short and long interspersed elements (SINEs and LINEs), and human endogenous retroviruses (HERVs). The presence of retroelements in the MHC provides a clear example of how these elements affect the genome plasticity of the host. Comparative analyses of these retroelements have proven highly useful in evolutionary studies of the MHC. Recently, HERV-encoded superantigens have been implicated as candidate autoimmune genes in type I diabetes and multiple sclerosis. In addition, genetic analyses have revealed that autoimmune diseases show strong associations with MHC class II genes. The intriguing correlations between retroviral encoded antigens, MHC class II genes and the development of autoimmune disease merit intense future investigations of retroelements, in particular those endogenous retroviruses located in the MHC class II region proper.

Presence of retroelements reveal the evolutionary history of the human DR haplotypes

Comparison of intron sequences has been a successful tool for drawing major conclusions about the evolutionary relationship of DRB genes. This complex family of genes is discussed in this review as well as a proposed model for the evolution of HLA-DR haplotypes. The model is based both on phylogenetic analysis of intron sequences as well as presence of ERV9 LTR elements located at identical position in intron 5 of a number of DRB genes. According to this model, two main evolutionary branches of DR haplotypes exist. The DR53 haplotype represents one branch, and the second branch contains the DR51, DR52, DR1, and DR8 haplotypes. After the divergence of the DR53 haplotype, an ERV9 LTR element was inserted in a primordial gene. Consequently, all DRB1 genes as well as the DRB3 gene within haplotypes of the second branch, contain this LTR element. In addition, conserved regulatory sequence motifs are found present within these LTR elements that might regulate DRB gene expression. Novel haplotypes are generated by recombinations and the maintenance of the DR haplotype variation as well as the frequent genetic rearrangements observed might be evolutionary advantageous.

High-resolution HLA typing for the DQB1 gene by sequence-based typing

The ideal high-resolution typing strategy for polymorphic genes is sequence-based typing. SBT of genomic DNA has been developed for the HLA class II genes DRB1, DRB3/4/5 and DPB1. For the DQB1 gene the sequence-based typing method was shown to cause a number of problems. To resolve those problems, different primers to amplify and sequence exon 2 of DQB1 were designed and tested. With several primer combinations, preferential amplification was observed in individuals heterozygous for DQB1*02/*03 and DQB1*02/*04. The preference was for DQB1*02 in many instances but could also be demonstrated for DQB1*03 or *04 and resulted occasionally in allelic drop-out. The best primer combination was selected and successfully used to type individuals heterozygous for DQB1*02, *03 and *04. To distinguish DQB1*0201 and *0202, primers for amplification and sequencing of exon 3 were developed and correct subtyping was obtained. The ambiguous typing DQB1*0301/*0302 and DQB1*0303/*0304 was resolved by allele-specific amplification and sequencing. A total of 258 individuals were fully typed for their DQB1 subtypes. All samples had been previously typed by PCR-SSP and serology. Concordant typing results were obtained for all individuals tested. The DQB1 alleles detected included *0501, *0502, *0503, *0601, *0602, *0603, *0604, *0609, *0201, *0202, *0301, *0302, *0303, *0304, *0401 and *0402. Sequence-based typing of the DQB1 gene proved a reliable typing strategy for assignment of the different DQB1 alleles after intensive selection of primers and test conditions.

The absence of DR51 in a DRB5-positive individual DR2ES is caused by a null allele (DRB5*0108N)

DR51, a protein encoded by the DRB5 gene, was shown to be present in almost all DR2-positive haplotypes. Exceptions were reported, some DR2-negative samples were shown to be DR51 positive and in a number of DR2-positive samples no DR51 antigen could be demonstrated. In some of them lack of the DRB5 gene was the cause of the absence of DR51 but in others the DRB5 gene was present without resulting in a detectable gene product. Many of these variants were studied in detail in previous international workshops. One of them was DR2ES from our laboratory. She is a DR15-positive DR51-negative individual of oriental origin with a clearly demonstrable DRB5*01 allele when typed by molecular techniques. To unravel the molecular mechanism responsible for the defect in expression, cDNA and DNA encoding the defective DRB5 allele were analyzed. Nucleotide sequence analysis of exon 2 showed no differences from the sequence of DRB5*0102. However, when exon 3 was examined a difference in length was noticed due to a deletion of 19 nucleotides between codon 161 and 168. The deletion caused a frameshift and a premature stopcodon resulting in a null allele. The same allele could be demonstrated in 6 other unrelated individuals of oriental origin as well as in 5 individuals from South Africa. The absence of the DR51 protein was explained by the presence of an alteration in the DRB5 allele resulting in a null allele. The allele has been officially named DRB5*0108N. This is the first description of a null allele of the DRB5 gene.

HLA-DRB4 genotyping by PCR-RFLP: diversity in the associations between HLA-DRB4 and DRB1 alleles

The serologically defined HLA-DR53 antigen is associated with HLA-DR4, -DR7, and -DR9 antigens, and these haplotypes contain two functional genes, DRB1 and DRB4, and two pseudogenes, DRB7 and DRB8. The DRB4 gene encodes the DR53 antigen, and has been officially recognized to contain three allelic variants (DRB4*0101, 0102, and 0103). In this study, we have established the polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) method for DRB4 genotyping and analyzed genetic polymorphism of the DRB4 gene in Japanese population. DRB4*0101, DRB4*0102, and DRB4*0103 could be observed at the frequencies of 0.5%, 1.1% and 32.7%, respectively. The same DRB1 allele does not necessarily share an identical DRB4 allele. Further, a tight linkage disequilibrium was found between DRB4*0102 and DRB1*0401 in Japanese population, whereas DRB1*0401 was associated with DRB4*0101 or *0103 in Caucasian population. These findings reveal extensive diversity of the HLA-DRB1 and -DRB4 haplotypes and may have important implications for HLA-disease associations and donor selection in unrelated transplantation.

Molecular analysis of the DLA DR region

The dog is an important model for studying organ transplantation. In order to develop a DNA-based typing system, a genomic analysis of the canine DR region of the MHC was undertaken. Southern analysis suggests the presence of two DRB genes in most dogs and all have one DRA gene. One dog out of 200 examined contains only one DRB gene. The DRA gene was mapped in one lambda phage clone. One DRB gene (DRBB1) was localized to two overlapping phage clones. Another DRB gene (DRBB2) was localized to two overlapping phage clones. Sequence analysis of the two DRB genes suggests that one gene is intact (DRBB1) and one gene is a pseudogene (DRBB2) because it lacks a splice acceptor signal for exon 3 and lacks exons 2 and exon 4. Intron 1 and 2 sequence data from DRBB1 allow amplification of the polymorphic exon 2 which, in turn, can serve as a basis for developing a typing system for DRB.

HLA-DRB intron 1 sequences: implications for the evolution of HLA-DRB genes and haplotypes

Human DRB genes encode beta chains of the major histocompatibility complex (MHC) class II molecules. Although nine DRB loci have been mapped to the short arm of chromosome 6, an individual chromosome contains only one to five loci and is classified into one of five major haplotypes. To elucidate the origin of human DRB loci and haplotypes, intron 1 sequences approximately 5000 bp in length were determined for three DRB1 alleles (DRB1*03, DRB1*04, and DRB1*15) and five DRB genes (DRB2, DRB3, DRB4, DRB5, and DRB7). The sequences were subjected to phylogenetic analyses together with previously determined intron 4 and 5 sequences. The sequences provided two sources of information: Nucleotide substitutions that could be used to construct phylogenetic trees and to estimate divergence times and a set of insertions (mostly Alu elements) that reveal the order of splitting of duplicated genes. The combined data indicate that the ancestor of the human DRB genes was HLA-DRB1*04-like and that the DRB2, DRB7, DRB5, and DRB3 genes arose from this ancestor by four rounds of duplication 58, 56, 53, and 36 million years (MY) ago, respectively. The DRB4 gene may have arisen 46 MY ago by a deletion from the DRB1 and DRB2 genes and the DRB6 gene is probably an allele at the DRB2 locus. During the course of its evolution, the DRB1*04 gene acquired an intron 1 segment (including two Alu elements) from a gene that became the ancestor of DRB1*03. The present-day HLA-DR haplotypes were derived from three principal ancestral haplotypes: DRB1-DRB2, DRB1-DRB5, and DRB1-DRB7.

Stem-loop potential in MHC genes: a new way of evaluating positive Darwinian selection?

The domains of polymorphic major histocompatibility complex (MHC) proteins which interact with peptides and T-cell receptors are considered to have been under positive evolutionary selection pressure. Evidence for this is a high ratio of non-synonymous to synonymous mutations in the corresponding genomic domains. By this criterion snake venom phospholipase A2 genes have also been under positive selection pressure. Recent studies of the latter genes indicate that positive selection has overridden an evolutionary pressure on base order which normally promotes the potential to extrude single-strand stem-loops from supercoiled duplex DNA ( fold pressure ). This has resulted in base order-dependent stem-loop potential being shifted to introns, which are highly conserved between species. It is now shown that, like snake venom phospholipase A2 genes, the domains of polymorphic MHC genes which appear to have responded to positive selection pressure have decreased base order-dependent stem-loop potential. The evolutionary pressure to generate stem-loop potential (believed to be important for recombination) has been overridden less in exons under negative purifying selection than in exons under positive Darwinian selection. Thus, base order-dependent stem-loop potential shows promise as an independent indicator of positive selection.

Evolutionary relationship of HLA-DRB genes inferred from intron sequences

The major histocompatibility complex (Mhc) consists of class I and class II genes. In the human Mhc (HLA) class II genes, nine DRB loci have been identified. To elucidate the origin of these duplicated loci and allelic divergences at the most polymorphic DRB1 locus, introns 4 and 5 as well as the 3' untranslated region (altogether approximately 1,000 base pairs) of seven HLA-DRB loci, three HLA-DRB1 alleles, and nine nonhuman primate DRB genes were examined. It is shown that there were two major diversification events in HLA-DRB genes, each involving gene duplications and allelic divergences. Approximately 50 million years (my) ago, DRB1*04 and an ancestor of the DRB1*03 cluster (DRB1*03, DRB1*15, and DRB3) diverged from each other and DRB5, DRB7, DRB8, and an ancestor of the DRB2 cluster (DRB2, DRB4, and DRB6) arose by gene duplication. Later, about 25 my ago, DRB1*15 diverged from DRB1*03, and DRB3 was duplicated from DRB1*03. Then, some 20 my ago, the lineage leading to the DRB2 cluster produced two new loci, DRB4 and DRB6. The DRB1*03 and DRB1*04 allelic lineages are extraordinarily old and have persisted longer than some duplicated genes. The orthologous relationships of DRB genes between human and Old World monkeys are apparent, but those between Catarrhini and New World monkeys are equivocal because of a rather rapid expansion and contraction of primate DRB genes by duplication and deletion.

Evolutionary relationship between human major histocompatibility complex HLA-DR haplotypes

HLA-DR haplotypes of the human major histocompatibility complex are organized in five different groups. They can be identified based on the serological specificity expressed by the polymorphic DRB1 locus and by the presence of a characteristic set of DRB genes. The nucleotide sequences of introns 4 and 5 of the two DRB genes (DRB1(*)01 and DRB6(*)01 ) from a DR1 haplotype and the three DRB genes (DRB1(*)15, DRB6(*)15 , and DRB5(*)15 ), from a DR51 haplotype were determined. This study identified endogenous retroviral long terminal repeat elements (ERV9 LTR) located at identical positions in intron 5 of the DRB1 genes in both the DR1 and DR51 haplotypes. Phylogenetic analyses revealed a close evolutionary relationship between these two haplotypes. The DRB5 gene, unique for the DR51 haplotype, may have been lost by a recent gene deletion event creating the DR1 haplotype. A model for the evolution of the human DR haplotypes involving separate duplication and contraction events is presented.

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.

Analysis of the fine specificities of sheep major histocompatibility complex class II-specific monoclonal antibodies using mouse L-cell transfectants

The fine specificities of two panels of monoclonal antibodies (mAbs) for sheep major histocompatibility complex (MHC) class II molecules were determined using five mouse L-cell transfectants, each expressing a defined sheep DQ or DR MHC class II A/B gene pair. Using the transfectants in an indirect fluorescence antibody assay, previous immunochemical characterization of the mAbs was confirmed for 16 of 23 mAbs tested. The MHC class II subtype specificity (DQ or DR) of each mAb was assigned without interference from the products of other expressed class II loci. This allowed the identification of both cross-locus specificities as well as defining fine specificities of mAbs previously only partially characterized by immunochemical techniques.