Simplifying genetic locus assignment of HLA-DRB genes

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.

Human Leukocyte Antigen (HLA) Testing in Pharmacogenomics

The genetic region on the short arm of chromosome 6 where the human leukocyte antigen (HLA) genes are located is the major histocompatibility complex. The genes in this region are highly polymorphic, and some loci have a high degree of homology with other genes and pseudogenes. Histocompatibility testing has traditionally been performed in the setting of transplantation and involves determining which specific alleles are present. Several HLA alleles have been associated with disease risk or increased risk of adverse drug reaction (ADR) when treated with certain medications. Testing for these applications differs from traditional histocompatibility in that the desired result is simply presence or absence of the allele of interest, rather than determining which allele is present. At present, the majority of HLA typing is done by molecular methods using commercially available kits. A subset of pharmacogenomics laboratories has developed their own methods, and in some cases, query single nucleotide variants associated with certain HLA alleles rather than directly testing for the allele. In this chapter, a brief introduction to the HLA system is provided, followed by an overview of a variety of testing technologies including those specifically used in pharmacogenomics, and the chapter concludes with details regarding specific HLA alleles associated with ADR.

State-of-the-art genome inference in the human MHC

The Major Histocompatibility Complex (MHC) on the short arm of chromosome 6 is associated with more diseases than any other region of the genome; it encodes the antigen-presenting Human Leukocyte Antigen (HLA) proteins and is one of the key immunogenetic regions of the genome. Accurate genome inference and interpretation of MHC association signals have traditionally been hampered by the region's uniquely complex features, such as high levels of polymorphism; inter-gene sequence homologies; structural variation; and long-range haplotype structures. Recent algorithmic and technological advances have, however, significantly increased the accessibility of genetic variation in the MHC; these developments include (i) accurate SNP-based HLA type imputation; (ii) genome graph approaches for variation-aware genome inference from next-generation sequencing data; (iii) long-read-based diploid de novo assembly of the MHC; (iv) cost-effective targeted MHC sequencing methods. Applied to hundreds of thousands of samples over the last years, these technologies have already enabled significant biological discoveries, for example in the field of autoimmune disease genetics. Remaining challenges concern the development of integrated methods that leverage haplotype-resolved de novo assembly of the MHC for the development of improved MHC genotyping methods for short reads and the construction of improved reference panels for SNP-based imputation. Improved genome inference in the MHC can crucially contribute to an improved genetic and functional understanding of many immune-related phenotypes and diseases.

Next-generation sequencing of 11 HLA loci in a large dengue vaccine cohort from the Philippines

HLA genotyping by next-generation sequencing (NGS) has evolved with significant advancements in the last decade. Here we describe full-length HLA genotyping of 11 loci in 612 individuals comprising a dengue vaccine cohort from Cebu province in the Philippines. The multi-locus individual tagging NGS (MIT-NGS) method that we developed initially for genotyping 4-6 loci in one MiSeq run was expanded to 11 loci including HLA-A, B, C, DPA1, DPB1, DQA1, DQB1, DRB1, and DRB3/4/5. This change did not affect the overall coverage or depth of the sequencing reads. HLA alleles with frequencies greater than 10% were A*11:01:01, A*24:02:01, A*24:07:01, A*34:01:01, B*38:02:01, B*15:35, B*35:05:01, C*07:02:01, C*04:01:01, DPA1*02:02:02, DPB1*05:01:01, DPB1*01:01:01, DQA1*01:02:01, DQA1*06:01:01, DQB1*05:02:01, DQB1*03:01:01, DRB1*15:02:01, DRB1*12:02:01, DRB3*03:01:03, DRB4*01:03:01, and DRB5*01:01:01. Improvements in sequencing library preparation provide uniform and even coverage across all exons and introns. This has led to a marked reduction in allele imbalance and dropout. Furthermore, including more loci, such as DRB3/4/5, decreases cross-mapping and incorrect allele assignment at the DRB1 locus. The increased number of loci sequenced for each sample does not reduce the number of samples that can be multiplexed on a single MiSeq run and is therefore more cost-efficient. We believe that such improvements will help HLA genotyping by NGS to gain momentum over other conventional methods by increasing confidence in the calls.

The prevalence of antibodies against the HLA-DRB3 protein in kidney transplantation and the correlation with HLA expression

Human leukocyte antigen (HLA)-DRB3 is a functional HLA class II gene, which has a limited allele diversity in the human population. Furthermore, the HLA-DRB3 gene is only present in a subset of individuals. Therefore, in organ transplantation, this HLA molecule is frequently mismatched between patient and graft donor and thus antibodies against this mismatched HLA molecule can develop. In this study, we aimed to evaluate the prevalence and reactivity of these antibodies and aimed to identify factors that underlie antibody formation against HLA-DRB3. We showed in our patient cohort that HLA-DRB3 antibodies are identified in about 7% of all patients that were screened with solid phase assays. In these assays, we observed multiple antibody reactivity patterns indicating that HLA-DRB3 harbours multiple epitopes. In those cases, where we succeeded at tracing back the induction of these antibodies to the molecular HLA typing of the immunogenic event, we noticed a different frequency of HLA-DRB1 allele groups in the donors as compared to a control group. To a certain extent this distribution (e.g. HLA-DRB1*11 individuals) could be linked to an altered expression level. However, it also appears that different HLA-DRB3 alleles (e.g. HLA-DRB3*01 group) vary in their immunogenicity without having an expression difference. In conclusion, our study provides information on the immunogenicity and reactivity patterns of antibodies against HLA-DRB3 in kidney transplantation, and it points towards the possibility of HLA expression as a factor underlying antibody formation.

Polymorphism of HLA class I and class II alleles in influenza A(H1N1)pdm09 virus infected population of Assam, Northeast India

Human leucocyte antigen (HLA) represents one of the most highly polymorphic systems which plays a central role in the immune response. Genetic polymorphism of HLA in influenza A(H1N1)pdm09 infected population may be an important factor in disease progression and severity that needs further probing. In this study, a total of 110 Influenza like illness patients were recruited from the population of Assam, Northeast India, from which 35 cases infected by A(H1N1)pdm09 viruses and 35 controls were typed for HLA-A, B and DRB1 locus by PCR-SSP method. A total of seven alleles of HLA-A, 16 alleles of HLA-B, and 11 alleles of HLA-DRB1 locus were identified. The most common alleles within each locus in cases were HLA-A*11 (85.71%, P = 0.046), HLA-B*35 (25%, P = 0.0001), and HLA-DRB1*15 (49.35%,  P = 0.133) as compared to the controls, HLA-A*11 (40.82%), HLA-B*35 (0.00%), and HLA-DRB1*15 (67.53%). The frequency of HLA-A*11 and HLA-B*35 were significantly higher in cases as compared to the controls. In DRB1 locus, HLA-DRB1*10 was significantly higher in cases (20.78%, P = 0.005) than that of controls (0.00%). Whereas, HLA-DRB1*15 showed a higher frequency in controls than in cases. In addition, HLA-DRB3*01 (P = 0.053), DRB4*01 (P = 1.000), and DRB5*01(P = 0.591) were also identified along with HLA-DRB1 haplotype. From this preliminary study, it is suspected that there may be a role of HLA-A*11, HLA-B*35 and HLA-DRB1*10 in conferring susceptibility to influenza A(H1N1)pdm09 infection in the study population. A larger extended study on HLA polymorphism may explain the association between HLA and influenza A(H1N1)pdm09 infection and provide insights for HLA restricted peptide based vaccines.

Imputing Variants in HLA-DR Beta Genes Reveals That HLA-DRB1 Is Solely Associated with Rheumatoid Arthritis and Systemic Lupus Erythematosus

The genetic association of HLA-DRB1 with rheumatoid arthritis (RA) and systemic lupus erythematosus (SLE) is well documented, but association with other HLA-DR beta genes (HLA-DRB3, HLA-DRB4 and HLA-DRB5) has not been thoroughly studied, despite their similar functions and chromosomal positions. We examined variants in all functional HLA-DR beta genes in RA and SLE patients and controls, down to the amino-acid level, to better understand disease association with the HLA-DR locus. To this end, we improved an existing HLA reference panel to impute variants in all protein-coding HLA-DR beta genes. Using the reference panel, HLA variants were inferred from high-density SNP data of 9,271 RA-control subjects and 5,342 SLE-control subjects. Disease association tests were performed by logistic regression and log-likelihood ratio tests. After imputation using the newly constructed HLA reference panel and statistical analysis, we observed that HLA-DRB1 variants better accounted for the association between MHC and susceptibility to RA and SLE than did the other three HLA-DRB variants. Moreover, there were no secondary effects in HLA-DRB3, HLA-DRB4, or HLA-DRB5 in RA or SLE. Of all the HLA-DR beta chain paralogs, those encoded by HLA-DRB1 solely or dominantly influence susceptibility to RA and SLE.

No gender differences in the frequencies of HLA-DRB3/B4/B5 heterozygotes in newborns and adults in Koreans

HLA class II haplotypes often contain a second expressed HLA-DRB locus tightly linked to the classical HLA-DRB1 locus on the haplotype, which can be either HLA-DRB3, -DRB4 or -DRB5. These encode the HLA-DR51, -DR52 or -DR53 supertypic specificities and mark the ancestral lineages. HLA-DRB3/B4/B5 heterozygote excess in Welsh male newborns has been reported, suggesting a possibility of male-specific major histocompatibility complex (MHC)-mediated prenatal selection. However, it has not been confirmed in newborns of other ethnic groups or in adult populations. We analyzed the HLA-DRB1 and HLA-DRB3/B4/B5 genes in Korean newborns and healthy adults to examine whether MHC-mediated prenatal or postnatal selection exists. A total of 1,038 newborns (cord blood registry, 516 males and 522 females) and 2,082 healthy adults (hematopoietic stem cell donor registry, 1,111 males and 971 females) were HLA typed. HLA-DRB1/B3/B4/B5 DNA typing was performed using Dynal RELI HLA-DRB SSO Kit (Dyanl Biotech, Wirral, U.K.). Genotype frequencies and homozygosity and heterozygosity rates for DRB3/B4/B5 supertypic loci were compared between males and females in newborns and adults. There were no significant differences in the HLA-DRB3/B4/B5 homozygosity and heterozygosity rates between males and females in both newborns and adults. In the comparison between newborns and adults, homozygosity rate was significantly higher in newborn females than in adult females (31.0% vs 25.0%, p=0.01). Whether there is an age-related change from newborns toward adults has not been well studied in other populations, and further studies are warranted. In conclusion, male-specific heterozygosity excess reported in Welsh newborns has not been observed in Korean population, and there might be some ethnic differences in the gender-specific prenatal selection events.

Strategies for developing multi-epitope, subunit-based, chemically synthesized anti-malarial vaccines

An anti-malarial vaccine against the extremely lethal Plasmodium falciparum is desperately needed. Peptides from this parasite's proteins involved in invasion and having high red blood cell-binding ability were identified; these conserved peptides were not immun genic or protection-inducing when used for immunizing Aotus monkeys. Modifying some critical binding residues in these high-activi binding peptides' (HABPs') attachment to red blood cells (RBC) allowed them to induce immunogenicity and protection against expermental challenge and acquire the ability to bind to specific HLA-DRp1* alleles. These modified HABPs adopted certain characterist structural configurations as determined by circular dichroism (CD) and ¹H nuclear magnetic resonance (NMR) associated with certain HLA-DRβ1* haplotype binding activities and characteristics, such as a 2-Å-distance difference between amino acids fitting into HLA-DRp1 Pockets 1 to 9, residues participating in binding to HLA-DR pockets and residues making contact with the TCR, suggesting haplotyp and allele-conscious TCR. This has been demonstrated in HLA-DR-like genotyped monkeys and provides the basis for designing high effective, subunit-based, multi-antigen, multi-stage, synthetic vaccines, for immediate human use, malaria being one of them.

Variations in the electrostatic landscape of class II human leukocyte antigen molecule induced by modifications in the myelin basic protein peptide: a theoretical approach

The receptor-ligand interactions involved in the formation of the complex between Class II Major Histocompatibility Complex molecules and antigenic peptides, which are essential for establishing an adaptive immunological response, were analyzed in the Class II Human Leukocyte Antigen (HLA) - Myelin Basic Protein (MBP) peptide complex (HLA-DRβ1*1501-MBP) using a multipolar molecular electrostatic potential approach. The Human Leukocyte Antigen - peptide complex system was divided into four pockets together with their respective peptide fragment and the corresponding occupying amino acid was replaced by each of the remaining 19 amino acids. Partial atomic charges were calculated by a quantum chemistry approach at the Hatree Fock/3-21*G level, to study the behavior of monopole, dipole and quadrupole electrostatic multipolar moments. Two types of electrostatic behavior were distinguished in the pockets' amino acids: “anchoring” located in Pocket 1 and 4, and “recognition” located in Pocket 4 and 7. According to variations in the electrostatic landscape, pockets were ordered as: Pocket 1>Pocket 9≫Pocket 4≈Pocket 7 which is in agreement with the binding ability reported for Class II Major Histocompatibility Complex pockets. In the same way, amino acids occupying the polymorphic positions β13R, β26F, β28D, β9W, β74A, β47F and β57D were shown to be key for this Receptor-Ligand interaction. The results show that the multipolar molecular electrostatic potential approach is appropriate for characterizing receptor-ligand interactions in the MHC–antigenic peptide complex, which could have potential implications for synthetic vaccine design.

Functional, structural, and immunological compartmentalisation of malaria invasive proteins

Conserved Plasmodium falciparum merozoite high activity binding peptides (HABPs) involved in red blood cell (RBC) invasion which are present in merozoite surface proteins (MSPs) involved in attachment, rolling over RBC, those derived from soluble proteins loosely bound to the membrane, and those present in microneme and rhoptry organelles have an alpha-helical structure and bind with high affinity to HLA-DR52 molecules. On the contrary, conserved HABPs belonging to molecules anchored to the membrane by a GPI tail, or a transmembranal region, or those molecules presenting PEXEL motifs have a strand, turn or unordered configuration and bind with high affinity to HLA-DR53 molecules. Such functional, cellular, structural, and immunological compartmentalisation has tremendous implications in subunit-based, multi-epitope, synthetic, anti-malarial vaccine development.

Conserved extended haplotypes of the major histocompatibility complex: further characterization

Since the complete sequencing of a human major histocompatibility complex (MHC) haplotype, interest in non-human leucocyte antigen (HLA) genes encoded in the MHC has been growing. Non-HLA genes, which outnumber the HLA genes, may contribute to or account for HLA and disease associations. Most information on non-HLA genes has been obtained in separate studies of individual loci. To comprehensively address polymorphisms of relevant non-HLA genes in 'conserved extended haplotypes' (CEH), we investigated 101 International Histocompatibility Workshop reference cell lines and nine additional anonymous samples representing all 37 unambiguously characterized CEHs at MICA, NFKBIL1, LTA, NCR3, AIF1, HSPA1A, HSPA1B, BF, NOTCH4 and a single nucleotide polymorphism (SNP) at HLA-DQA1 as well as MICA, NOTCH4, HSPA1B and all five tumour necrosis factor short tandem repeat (STR) polymorphisms. This work (1) provides an extensive catalogue of MHC polymorphisms in all CEHs, (2) unravels interrelationships between HLA and non-HLA haplotypical lineages, (3) resolves reported typing ambiguities and (4) describes haplospecific markers for a number of CEHs. Analysis also identified a DQA1 SNP and segments containing MHC class III polymorphisms that corresponded with class II (DRB3 and DRB4) lineages. These results portray the MHC where lineages containing non-HLA and HLA variants in linkage disequilibrium may operate in concert and can guide more thorough design and interpretation of HLA-disease relationships.

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.

Duplication, balancing selection and trans-species evolution explain the high levels of polymorphism of the DQA MHC class II gene in voles (Arvicolinae)

Major histocompatibility complex (MHC) genes play important role in host-parasite interactions and parasites are crucial factors influencing the population dynamics of hosts. We described the structure and diversity of exon 2 of the MHC class II DQA gene in three species of voles (Arvicolinae) exhibiting regular multi-annual fluctuations of population density and analysed the processes leading to the observed MHC polymorphism. By using cloning-sequencing methodology and capillary electrophoresis-single strand conformation polymorphism, we described seven sequences in the water, eight in the common, and seven in the bank voles coming from an area of 70 km(2) around the Nozeroy canton in the Jura Mountains (Franche Comté, France). All exon 2 sequences translate to give unique amino acid sequences and positive selection was found to act very intensively on antigen binding sites. We documented the presence of recombination at vole DQA region but its importance in generating allelic polymorphism seems to be relatively limited. For the first time within rodents, we documented the duplication of the DQA gene in all three species with both copies being transcriptionally active. Phylogenetic analysis of allelic sequences revealed extensive trans-species polymorphism within the subfamily although no alleles were shared between species in our data set. We discuss possible role of parasites in forming the recent polymorphism pattern of the DQA locus in voles.

Features of MHC and NK gene clusters

When comparing the immune genome to the genome in general, a higher prevalence for association with disease is the only genetic feature significant in immune genes as a group. However, some genetic features, such as marked levels of polymorphism and gene duplication, are present in subsets of immune genes, namely the Major Histocompatibility Complex (MHC) and Natural Killer (NK) cell receptor gene complexes. In this review, we discuss features of MHC and NK receptor gene clusters, their epistatic interactions, and the impact of both on association to disease.

Analysis of eluted peptides from type 1 diabetes-susceptible HLA class II molecules identified novel islet protein, heparin/heparan sulfate-interacting protein

Identification of peptides derived from pancreatic islet and presented by type 1 diabetes-susceptible MHC class II molecules has great significance to elucidate the pathogenesis of type 1 diabetes. A bulk culture of Epstein-Barr virus-transformed B-cells, which were established from a 22-year-old type 1 diabetic woman with HLA-DR4 and -DQw8, was pulsed with the homogenate of a human embryonic pancreas-derived cell line 1B2C6, and another culture was not pulsed with antigen. Peptide fractions were obtained by treatment of affinity-purified HLA-DR and -DQ molecules with 0.1% trifluoroacetic acid, and were subjected to reverse-phase high performance liquid chromatography (RP-HPLC). The RP-HPLC profiles of peptides derived from DR molecules revealed three peaks that specifically appeared after pulsing, but no such peaks were obtained from DQ molecules. From one of these three peaks, a peptide that consisted of 14 amino acids (AKSXNHTXXNQXRK, where X represents the undetermined amino acids) was identified. This peptide was derived from heparin/heparan sulfate-interacting protein (HIP). Immunostaining of pancreatic sections using antiserum for HIP peptide revealed exclusive staining of the islets. Thus, HIP was identified as an islet protein naturally processed and presented by HLA-DR4 molecules.

Increased heterozygosity for MHC class II lineages in newborn males

In plants, fungi and marine invertebrates, there are genetic compatibility systems to ensure diversity in the offspring. The importance of genetic compatibility in gametic union and selective abortion in vertebrate animals has also been appreciated recently. There have been suggestions that the major histocompatibility complex (HLA in humans) may be a compatibility system in vertebrates. HLA class II haplotypes often contain a second expressed DRB locus which can be either DRB3, DRB4 or DRB5. These encode the supertypical specificities and mark the ancestral lineages. The members of each lineage have related DNA sequences at the main class II locus HLA-DRB1. We analysed 415 newborns at all expressed DRB loci by PCR analysis to seek evidence for sex-specific prenatal selection events. While there was no significant change in heterozygosity rates between males and females at DRB1, the proportion of males carrying two DRB1 specificities from different ancestral lineages was significantly increased (53.7% in males vs 39.3% in females, P = 0.003). The genotypes consisting of phylogenetically most distinct ones, namely the DRB3 and DRB4 haplotypes, showed the most striking difference between sexes (P = 0.007). These results suggested a more favourable outcome for male concepti heterozygous for supertypical haplotypes. Heterozygosity for most divergent haplotypical families ensures the highest degree of functional heterozygosity at the main HLA class II locus DRB1 while increasing the likelihood of heterozygosity also at other MHC loci. Our observations agree with the previously reported heterozygote excess in male newborn rats and mice. Correlations between MHC class II heterozygosity and advertised male quality in deer and pheasant as well as increased reproductive success in MHC class II heterozygous male macaques are examples of postnatal benefits of heterozygosity in males that may be behind the development of prenatal selection mechanisms. The MHC-mediated prenatal selection of males may also be one of the selective events suggested by the very high primary (male-to-female) sex ratio at fertilization reaching close to unity at birth in humans. These results provide an appealing working hypothesis for further studies in humans and other vertebrates.

Effects of genotype matching of feline major histocompatibility complex (FLA) class II DRB on skin-allograft transplantation in cats

In order to confirm the effects of matching of expressed feline major histocompatibility complex (FLA) class II DRB genotype on transplantation immunity in cats, skin-allogeneic transplantation was carried out between cats, in which DRB genes expressed were genotyped by the RT-PCR-RFLP method using group-specific primers. Duration until grafts were rejected was 14.63 +/- 1.69 days (mean +/- standard deviation) in the pairs that had the same type of subgroups, 7.25 +/- 0.71 days in the pairs that had one different type of subgroup and 6.88 +/- 0.35 days in the pairs that had two different types of subgroups. The duration of graft survival in the pairs with the same type of subgroups was significantly longer (P<0.01) than those in the pairs with different types. Although FLA components involved in transplantation immunity should not only be DRB genes, it was suggested that the expressed FLA-DRB genotype might associate with feline transplantation immunity, and that typing and matching of expressed FLA-DRB genes might be one of the important factors in the control of feline transplantation immunity.

Mhc-DQ-DRB haplotype analysis in the rhesus macaque: evidence for a number of different haplotypes displaying a low allelic polymorphism

In the HLA-DRB subregion of man, five major groups of haplotypes, often displaying a remarkable polymorphism, are distinguishable. The polymorphism is thought to be generated by point mutation, microgene conversion and gene rearrangement by recombination. In order to gain insight into the organization of the rhesus macaque major histocompatibility complex (MHC) class II region, DRB genes from monkeys of different origins previously typed for their DQ genes were analyzed. At first DRB haplotypes were deduced from DQ-homozygous monkeys. The stability of these haplotypes was then examined in DQ-heterozygous monkeys by sequence-based typing for the presence of members of the DRB1*03 and DRB1*04 lineage, and for seven single alleles detected on the haplotypes. Six DRB haplotypes linked to the five most frequent and three haplotypes linked to less frequent DQ haplotypes were identified. Six novel DRB alleles were detected. The number of DRB genes per haplotype varied between two and four. The results altogether suggest that in rhesus macaques, in comparison to man, the DQ haplotypes are linked to only a small number of DRB haplotypes, the number and diversity of DRB haplotypes is larger, and the allelic polymorphism of a given haplotype is smaller. The diversity of the DRB haplotypes was partly due to the varying number and identity of genes linked to DRB1*03 and DRB1*04. Furthermore, the number of DRB1 genes themselves varied from zero to two.

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.

DR non-B1 mismatches influence allogeneic MLR-induced TH1- or TH2-like cytokine responses in rhesus monkeys

Human and nonhuman primates have multiple DR B1 and non-B1 alleles. However, the role of mismatched DR non-B1 alleles in primary alloimmune responses is not well understood. Macaques, which share close DNA homologies with human MHC genes and have a high number of beta-chain genes in the DR subregion, are preeminent preclinical models for immunologic studies of transplant tolerance and immunosuppression. In this study, we examined the effect of allogeneic MHC Class II DRB mismatches in Th1- and Th2-like cytokine responses elicited in one-way MLR cultures in rhesus macaques. An ELISPOT method was used to estimate cytokine secretion at the single cell level. Molecular typing for DRB1 and DR non-B1 alleles was performed by a moderate-high resolution PCR-SSP method using a panel of 55 primer pairs covering 74 DRB alleles and clusters. Of 35 unrelated combinations, 66% had multiple (> or = 2) allelic MM at DRB1 and DR non-B1 with no significant correlation between numbers of DRB1 and DR non-B1 mismatches. Pairs with 1 or 0 MM were assigned to a mono/null MM group to obtain sufficient numbers for statistical analysis. The pairs differing by multiple vs. mono/null DRB1 MM showed no significant difference in cytokine prevalence (P = 0.69). In contrast, high IFN-gamma/ IL4 SFC ratios were noted in pairs with multiple vs. mono/null DR non-B1 MM (p = 0.0009). IFN-gamma/IL-10 spot forming cell (SFC) ratios were consistent with IFN-gamma/IL-4 SFC ratios (r = 0.98). Multiple DR non-B1 mismatches showed a trend towards higher MLR proliferative responses, although the stimulation index did not reflect the dominant cytokine response. These observations suggest a bias towards Th1-like cytokine production under allostimulation with multiple DR non-B1 gene products. Further study of the primary structure of DR non-B1 determinants may be helpful in understanding the fine molecular mechanisms governing the regulation of cytokine profiles during allostimulation in primates.

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.

MHC-DRB allelic sequences incorporate distinct intragenic trans-specific segments

The second exon of primate MHC-DRB genes encodes discrete areas of allelic hypervariability (HVR), which are used as the basis for lineage assignments to determine genetic and evolutionary relationships. Comparisons of these regions have led to the "trans-species hypothesis", which proposes that certain MHC alleles from one species are more closely related to those from other species than they are to each other; i.e., that allelic lineages are ancestral in origin. We evaluated this paradigm in an analysis of macaque and baboon MHC-DRB genes using oligotyping and sequencing of 87 new nonhuman primate DRB alleles. A remarkable conservation of sequence motifs in the HVRIII region (codon 60-79) was observed, detected both by hybridization and by sequencing; some of these motifs were found in species such as prosimians that have diverged from the human lineage 50 MYA. However, these fixed HVRIII motif sequences nevertheless occur on a background of diverse lineages suggesting that it is the segmental motif, rather than the allele per se which is trans-specific in origin. Sequences within the first hypervariable region (codons 7-14) identified lineage assignments to several DRB loci (DRB1, DRB3, DRB4, DRB5, DRB6 and DRB7), although a large number of DRB nucleotide sequences did not correspond to a defined allelic motif, suggesting that many of the nonhuman sequences lack human HVRI homologs and have accumulated additional intraspecies variation subsequent to speciation. While there are certain allelic lineages in HVRI that show trans-species conservation, other sequence motifs seem purely species-specific. These differences suggest that HVRI and HVRIII regions have distinct mechanisms for maintenance of trans-specific sequence elements, with different evolutionary histories for segmental nucleotide conservation.

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.

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.

Characterization of three separated exons in the HLA class II DR region of the human major histocompatibility complex

The human major histocompatibility complex, HLA, is a highly polymorphic gene region which includes the DRA and DRB genes. The number of DRB genes differs between haplotypes. The DR4 haplotype seems to be one of the most complex with five DRB loci, DRB1, DRB4, DRB7, DRB8, and DRB9, in addition to the single DRA locus. We determined the nucleotide sequences of three separated DRB exons located between the DRB4 locus and the DRA locus in the DR4 haplotype, two DRB signal-peptide exons (S1 and S3) and one DRB first-domain exon (locus designation DRB9). Sequence comparisons suggest the following order of events for the origin of these exons: DRB9 seems to be the oldest exon and has previously been detected in multiple HLA haplotypes. DRB9 is more divergent than the three other known DRB pseudogenes, all of which have been found in apes. This suggests that DRB9 arose prior to the hominoid divergence. An L1 repeat has been inserted 3' to DRB9. Subsequently, a LTR of the ERV9 retrovirus-like family was inserted into the L1 repeat. Such LTRs have recently been observed in some of the other DRB genes. The pseudogenes DRB7 and DRB8 (containing only exons 3-6) arose after DRB9. Finally, the separated signal peptide exons S1 and S3 were formed. The molecular characterization of these separated DRB exons and insertion elements further clarifies the complex evolutionary history of the HLA-DR region. These selectively neutral exons may serve as useful markers for tracing the phylogeny of HLA haplotypes.