Comparative anatomy of the primate major histocompatibility complex DR subregion: evidence for combinations of DRB genes conserved across species

Genetic variation in the human leukocyte antigen region confers susceptibility to Clostridioides difficile infection

Clostridioides difficile (C. diff.) infection (CDI) is a leading cause of hospital acquired diarrhea in North America and Europe and a major cause of morbidity and mortality. Known risk factors do not fully explain CDI susceptibility, and genetic susceptibility is suggested by the fact that some patients with colons that are colonized with C. diff. do not develop any infection while others develop severe or recurrent infections. To identify common genetic variants associated with CDI, we performed a genome-wide association analysis in 19,861 participants (1349 cases; 18,512 controls) from the Electronic Medical Records and Genomics (eMERGE) Network. Using logistic regression, we found strong evidence for genetic variation in the DRB locus of the MHC (HLA) II region that predisposes individuals to CDI (P > 1.0 × 10-14; OR 1.56). Altered transcriptional regulation in the HLA region may play a role in conferring susceptibility to this opportunistic enteric pathogen.

Characterization of MHC class II B polymorphism in multiple populations of wild gorillas using non-invasive samples and next-generation sequencing

Genes encoded by the major histocompatibility complex (MHC) are crucial for the recognition and presentation of antigens to the immune system. In contrast to their closest relatives, chimpanzees and humans, much less is known about variation in gorillas at these loci. This study explored the exon 2 variation of -DPB1, -DQB1, and -DRB genes in 46 gorillas from four populations while simultaneously evaluating the feasibility of using fecal samples for high-throughput MHC genotyping. By applying strict similarity- and frequency-based analysis, we found, despite our modest sample size, a total of 18 alleles that have not been described previously, thereby illustrating the potential for efficient and highly accurate MHC genotyping from non-invasive DNA samples. We emphasize the importance of controlling for multiple potential sources of error when applying this massively parallel short-read sequencing technology to PCR products generated from low concentration DNA extracts. We observed pronounced differences in MHC variation between species, subspecies and populations that are consistent with both the ancient and recent demographic histories experienced by gorillas.

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.

Major histocompatibility complex and microsatellite variation in two populations of wild gorillas

In comparison to their close relatives the chimpanzees and humans, very little is known concerning the amount and structure of genetic variation in gorillas. Two species of gorillas are recognized and while the western gorillas number in the tens of thousands, only several hundred representatives of the mountain gorilla subspecies of eastern gorillas survive. To analyse the possible effects of these different population sizes, this study compares the variation observed at microsatellite and major histocompatibility complex (MHC) loci in samples of wild western and mountain gorillas, collected using a sampling scheme that targeted multiple social groups within defined geographical areas. Noninvasive samples proved a viable source of DNA for sequence analysis of the second exon of the DRB loci of the MHC. Observed levels of variation at the MHC locus were similar between the two gorilla species and were comparable to those in other primates. Comparison of results from analysis of variation at multiple microsatellite loci found only a slight reduction in heterozygosity for the mountain gorillas despite the relatively smaller population size.

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.

Six mamu-A locus alleles defined by a polymerase chain reaction sequence specific primer method

Rhesus monkeys are relevant models for human diseases and transplantation. In each case, a complete understanding of these models requires knowledge of macaque MHC. Due to high polymorphism and multiple genes per haplotype, it has been difficult to develop a rapid typing method for rhesus monkey MHC class I. We developed a simple and rapid PCR-SSP strategy for rhesus monkey Mamu-A locus typing. Fifty-two rhesus monkeys were included in the study. Six rhesus monkey allel-specific primer pairs were designed based on published Mamu-A locus gene sequences. Allele-specific PCR products ranged in size from 346 to 788 bp; 5' and 3' Mamu-A locus allele specific primers were located in the second and third exons, respectively. Specific PCR product gel purification was followed by direct sequencing, without subcloning, in both directions. Our data showed variability in the number of Mamu-A alleles ranging from 1 to 4 per genotype. The highest frequencies were observed for Mamu-A*02, -A*04, and -A*03 alleles. Thus, we report here the first PCR-SSP typing method for Mamu-A*02, -03, -04, -05, -06, and -07 array of class I alleles. This technique appears to be a highly reproducible and discriminatory method for detecting this subset of class I A locus genes in rhesus monkeys.

Identification of DRB alleles in rhesus monkeys using polymerase chain reaction-sequence-specific primers (PCR-SSP) amplification

Major histocompatibility complex (MHC) class In molecules play a vital role in the regulation of T-cell functions in the mammalian immune system. Two key features characterize the polymorphism of MHC haplotypes in humans and non-human primates: the existence of a large number of alleles, and the high degree of genetic diversity between those alleles. Rhesus monkeys and Chimpanzees have been extensively used as relevant models for human diseases and transplantation We have investigated DRB genes in 19 macaques, members of 3 families, using polymerase chain reaction with sequence-specific primers (PCR-SSP) and denaturing gradient gel electrophoresis (DGGE). After amplification PCR products were purified and subjected direct sequencing. Seven animals (Madison #1) were typed by DDGE also. We report that the DRB haplotypes defined by PCR-SSP exhibit a high degree of concordance with the data obtained by DGGE and direct sequening. Our data show prominent variability in the number of DRB1 alleles ranging from 1-4 per genotype within these families. This analysis demonstrated that most of the amplicons were identical to Mamu-DRB alleles that our PCR primers were to amplify. However, 98-99% similarity was noticed in the case of Mamu-DRB1*0303, Mamu-DRB6*0103 and Mamu-DRB*W201 alleles. The observed mismatches were located in non-polymorphic regions. Thus, family studies in rhesus macaques performed by molecular methods confirmed the multiplicity of Mamu-DRB1 alleles per haplotype and the existence of allelic associations published earlier. In addition, we propose 3 more DRB allele associations (haplotypes): Mamu-DRB1*04-DRB5*03; Mamu-DRB1*04-*DRB*W5; Mamu-DRB1*04*W2. The proposed medium-resolution PCR-SSP technique appears to be a highly reproducible and discriminatory typing method for detecting polymorphisms of DRB genes in rhesus monkeys.

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.

The common marmoset: a new world primate species with limited Mhc class II variability

The common marmoset (Callithrix jacchus) is a New World primate species that is highly susceptible to fatal infections caused by various strains of bacteria. We present here a first step in the molecular characterization of the common marmoset's Mhc class II genes by nucleotide sequence analysis of the polymorphic exon 2 segments. For this study, genetic material was obtained from animals bred in captivity as well as in the wild. The results demonstrate that the common marmoset has, like other primates, apparently functional Mhc-DR and -DQ regions, but the Mhc-DP region has been inactivated. At the -DR and -DQ loci, only a limited number of lineages were detected. On the basis of the number of alleles found, the -DQA and -B loci appear to be oligomorphic, whereas only a moderate degree of polymorphism was observed for two of three Mhc-DRB loci. The contact residues in the peptide-binding site of the Caja-DRB1*03 lineage members are highly conserved, whereas the -DRB*W16 lineage members show more divergence in that respect. The latter locus encodes five oligomorphic lineages whose members are not observed in any other primate species studied, suggesting rapid evolution, as illustrated by frequent exchange of polymorphic motifs. All common marmosets tested were found to share one monomorphic type of Caja-DRB*W12 allele probably encoded by a separate locus. Common marmosets apparently lack haplotype polymorphism because the number of Caja-DRB loci present per haplotype appears to be constant. Despite this, however, an unexpectedly high number of allelic combinations are observed at the haplotypic level, suggesting that Caja-DRB alleles are exchanged frequently between chromosomes by recombination, promoting an optimal distribution of limited Mhc polymorphisms among individuals of a given population. This peculiar genetic make up, in combination with the limited variability of the major histocompatability complex class II repertoire, may contribute to the common marmoset's susceptibility to particular bacterial infections.

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.

DQA-DQB linkage in Old World monkeys

The MHC DQ region in nonhuman primates, as in humans, consists of alpha and beta chains that are polymorphic with strong linkage disequilibrium between certain DQA-DQB alleles. Not only are contemporary HLA class II allelic variants present in evolutionarily distant species, but we demonstrate that linkages between loci also bear ancient roots. In unrelated baboons (Papio cynocephalus anubis) and family segregation analysis of pigtailed macaques (Macaca nemestrina) we found cis-linkages between DQA1*01 and DQB1*05 or *06, between DQA1*05 and DQB1*03, and between DQA1*03 and DQB1*03 alleles, all of which are also prominent in modern humans. In contrast, one linkage that has not been seen in humans, between DQA1*05 and DQB1*06 alleles, was also found. These patterns of selective linkage disequilibrium imply evolutionary mechanisms following the divergence of species that constrain the diversity of haplotypes which evolve.

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.

Maternal-fetal interactions and MHC polymorphism

Two models of maternal-fetal interactions are discussed. In the first, offspring are advantaged if they possess an allele absent in their mother. Polymorphism is maintained because rare alleles have an advantage when present in males. In the second, offspring are disadvantaged if they lack an allele present in their mother. Polymorphism is maintained because rare alleles have an advantage when present in females. Both classes of model are associated with a deficiency of homozygous genotypes. If the artificial assumption of symmetrical selection is relaxed, the second class of model (gestational drive) could account for the otherwise inexplicable absence of MHC polymorphism in some species.

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.

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.

The HLA-DRB9 gene and the origin of HLA-DR haplotypes

HLA-DRB9 is a gene fragment consisting of exon 2 and flanking intron sequences. It is located at the extreme end of the DRB subregion, whose other end is demarcated by the DRB1 locus. We sequenced approximately 1400 base pairs of the segment encompassing the DRB9 locus from eight human haplotypes (DR1, DR10, DR2, DR3, DR5, DR6, DR8, and DR9, the DR4 and DR7 having been sequenced by others earlier), as well as two chimpanzee, five gorillas, one orangutan and one macaque haplotype. The analysis of these sequences indicates that the DRB9 locus, which we estimate to be more than 58 million years (my) old, has been coevolving with the DRB1 locus for the last 4.2 my. As a consequence of this coevolution, the human DRB9 alleles fall into groups that correlate with the DRB1 allelic groups and with the gene organization of the human haplotypes. This observation implies that the present-day HLA-DR haplotype groups (DR1, DR51, DR52, DR8, and DR53) were founded more than 4 my ago and have remained intact (barring minor internal rearrangements that did not recombine the DRB1 and DRB9 genes) for this period of time. The haplotypes have been transmitted during speciations from ancestral to emerging species just like allelic lineages at the DRB1 locus. Thus not only allelic but also haplotype polymorphism evolves trans-specifically.

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.

Allelic diversity at the Mhc-DP locus in rhesus macaques (Macaca mulatta)

Allelic diversity at the major histocompatibility complex class II DP locus of rhesus macaques was studied by sequencing exon 2 of Mamu-DPA1 and -DPB1 genes. The Mamu-DPA1 gene is apparently invariant, whereas the Mamu-DPB1 locus displays polymorphism. Here we report the characterization of 1 Mamu-DPA1 and 13 Mamu-DPB1 alleles which were compared with other available primate Mhc-DPA1 and -DPB1 sequences. As compared with Mhc-DRB and -DQB1, most codons for the contact residues in the antigen binding site of the primate Mhc-DPB1 gene have a relatively low degree of variation in encoding various types of amino acids. In contrast to Mhc-DRB and -DQB, the HLA- and Mamu-DPB1 sequences cluster in a species-specific manner in phylogenetic trees. Mhc-DPB1 polymorphisms, however, are inherited in a transspecies mode of evolution, as is demonstrated by the sharing of lineage members between closely related macaque species. The data demonstrate that the transspecies character of Mhc-DPB1 polymorphism was retained over much shorter periods of time as compared with its sister class II loci, Mhc-DQ and -DR.

Alu elements of the primate major histocompatibility complex

The chromosomal region constituting the major histocompatibility complex (MHC) has undergone complex evolution that is often difficult to decipher. An important aid in the elucidation of the MHC evolution is the presence of Alu elements (repeats) which serve as markers for tracing chromosomal rearrangements. As the first step toward the establishment of sets of evolutionary markers for the MHC, Alu elements present in selected MHC haplotypes of the human species, the gorilla, and the chimpanzee were identified. Restriction fragments of cosmid clones from the libraries of the three species were hybridized with Alu-specific probes, Alu elements were amplified by the polymerase chain reaction, and the amplification products were sequenced. In some cases, sequences of the regions flanking the Alu elements were also obtained. Altogether, 31 new Alu elements were identified, representing six Alu subfamilies. The average density of Alu elements in the MHC is one element per four kilobases (kb) of sequence. Alu elements have apparently been inserted steadily into the MHC over the last 65 million years (my). On average, one Alu element is inserted into the primate MHC every 4 my. Analysis of the human DR3 haplotype supports its origin by duplication from an ancestral haplotype consisting of DRB1 and DRB2 genes. The sharing of an old Alu element by the DRB1 and DRB2 genes, in turn, supports their divergence from a common ancestor more than 55 my ago.

The synonymous substitution rate of the major histocompatibility complex loci in primates

Because the divergence of many allelic lineages at the major histocompatibility complex (MHC) loci predates species divergence, standard methods of calculating synonymous substitution rates are not applicable to this system. We used three alternative methods of rate estimation: one based on the minimum number of substitutions (Dm), another on the nucleotide difference (Dxy), and the third on the net nucleotide difference (Dn). We applied these methods to the protein-encoding sequences of primate MHC class I (A, B, and C) and class II (DRB1) genes. To determine the reliability of the different estimates, we carried out computer simulation. The distribution of the estimates based on Dxy or Dn is generally much broader than that based on Dm. More importantly, the Dm-based method nearly always has the highest probability of recovering true rates, provided that Dm is not smaller than 5. Because of its desirable statistical properties, we used the Dm-based method to estimate the rate of synonymous substitutions. The rate is 1.37 +/- 0.61 for A, 1.84 +/- 0.40 for B, 3.87 +/- 1.05 for C, and 1.18 +/- 0.36 for DRB1 loci, always per site per 10(9) years. Hence despite the extraordinary polymorphism, the mutation rate at the primate MHC loci is no higher than that of other loci.

Conservative evolution of the Mbc-DP region in anthropoid primates

To determine the organization of the DP region in the Mbc of anthropoid primates, we constructed contig maps from cosmid clones of the chimpanzee and orangutan, representatives of the infraorder Catarrhini, as well as of the cotton-top tamarin, a representative of the infraorder Platyrrhini. We found the maps to be remarkably similar to each other and to the previously published map of the human DP region. In each of the four species, the DP region consists of four loci arranged in the same order (DPB2 . . . DPA2 . . . DPB1 . . . DPA1) and in the same transcriptional orientation (tail-to-tail). The regions in the four species are of approximately the same length and many of the restriction sites are shared between species. The inserts of most Alu elements, of a ribosomal protein pseudogene, and of an IgC epsilon-like pseudogene are found in corresponding positions in all four species. The data indicate that the human-type organization of the DP region was established before the divergence of the Catarrhini and Platyrrhini lines more than 37 million years ago and that it has remained principally intact since that time. This conservation of the DP region is in striking contrast to the evolutionary instability of certain other Mbc regions, in particular those occupied by the DRB or C4 and CYP21 loci. We interpret the stability of the DP region as an indication that the region is being phased out functionally.

Mhc-DRB genes of platyrrhine primates

The two infraorders of anthropoid primates, Platyrrhini (New World monkeys) and Catarrhini (Old World monkeys and the hominoids) are estimated to have diverged from a common ancestor 37 million years ago. The major histocompatibility complex class II DRB gene and haplotype polymorphism of the Catarrhini has been characterized in several recent studies. The present study was undertaken to obtain information on the DRB polymorphism of the Platyrrhini. Fifty-five complete exon 2 DRB sequences were obtained from six species of Platyrrhini representing both the Callitrichidae and the Cebidae families. Combined with the results of a parallel contig mapping study, our data indicate that at least three loci (DRB1*03, DRB3, and DRB5) are shared by the Catarrhini and the Platyrrhini. However, the three loci are occupied by functional genes in the former infraorder and mostly by pseudogenes in the latter. Instead of the pseudogenes, the Platyrrhini have evolved a new set of apparently functional genes-DRB11 and DRB*W12 through DRB*W19, which have thus far not been found in the Catarrhini. The DRB*W13, *W14, *W15, *W17, *W18, and *W19 genes seem to be restricted to the Cebidae family, whereas the DRB*W16 locus has so far been documented in the Callitrichidae family only. The DRB alleles of the cotton-top tamarin, and perhaps also those of the common marmoset (both members of the family Callitrichidae), are characterized by low nucleotide diversity, possibly indicating that they diverged from a common ancestral gene relatively recently.

Mhc-DRB and -DQA1 nucleotide sequences of three lowland gorillas. Implications for the evolution of primate Mhc class II haplotypes

Mhc-DRB and -DQA1 second-exon and -DRB 3'-untranslated-region nucleotide sequences of three lowland gorillas with no known family relationship with each other and of two HLA homozygous typing cell lines were determined and compared with published primate Mhc-DRB and -DQA1 sequences. Eleven distinct MhcGogo-DRB second-exon sequences were found, which represent the gorilla counterparts of the HLA-DRB1*03, -DRB1*10, -DRB3, -DRB5, and -DRB6 allelic lineages. One Gogo-DRB second-exon sequence does not have an obvious human counterpart and is tentatively designated Gogo-DRBY*01. The gorilla equivalents of the HLA-DRB2 and -DRB8 loci were identified as judged on Mhc-DRB 3'-untranslated-region sequences. In addition, four different Gogo-DQA1 alleles belonging to three different allelic lineages were detected. The Mhc-DRB-DQA1 haplotypes of these gorillas were deduced based on the obtained Mhc-DRB and -DQA1 sequences and the two published Mhc-DRB haplotypes of the lowland gorilla Sylvia. All deduced Gogo-DRB-DQA1 haplotypes show gene constellations different from known HLA-DRB-DQA1 haplotypes, while some of the Gogo-DRB haplotypes presented here contain more DRB genes than the HLA-DRB haplotypes. Based on phylogenetic trees, bootstrap analyses, and the gorilla, chimpanzee, and human Mhc-DRB haplotypes described, we propose that at least two Mhc-DRB loci, here tentatively designated Mhc-DRBI and -DRBII, existed on an ancient primate Mhc-DRB haplotype. The Mhc-DRB1*01, -DRB1*02 (-DRB1*15 and -DRB1*16), -DRB1*03 (-DRB1*03, -DRB1*08, -DRB1*11, -DRB1*12, -DRB1*13, and DRB1*14), and -DRB1*10 allelic lineages and -DRB3 and -DRBY loci probably evolved from the hypothetical primate Mhc-DRBI locus, whereas the present primate Mhc-DRB2, -DRB4, and -DRB6 loci originate from the ancient Mhc-DRBII locus of this core primate Mhc-DRB haplotype.

Multiplication of Mhc-DRB5 loci in the orangutan: implications for the evolution of DRB haplotypes

The beta chain-encoding (B) class II genes of the primate major histocompatibility complex belong to several families. The DRB family of class II genes is distinguished by the occurrence of haplotype polymorphism--the existence of multiple chromosomal forms differing in length, gene number, and gene combinations, each form occurring at an appreciable frequency in the population. Some of the haplotypes, or fragments thereof, are shared by humans, chimpanzees, and gorillas. In an effort to follow the DRB haplotype polymorphism further back in time, we constructed DRB contig maps of the two chromosomes present in the orangutan cell line CP81. Two types of genes were found in the two haplotypes, Popy-DRB5 and Popy-DRB1*03, the former occurring in two copies and one gene fragment in each haplotype, so that the CP81 cell line contains four complete DRB5 genes and two DRB5 fragments altogether. Since the four genes are more closely related to one another than they are to other DRB5 genes, they must have arisen from a single ancestral copy by multiple duplications. At the same time, however, the two CP81 haplotypes differ considerably in their restriction enzyme sites and in the presence of Alu elements at different positions, indicating that they have been separated for a length of time that exceeds the lifespan of a primate species. Moreover, a segment of about 100 kilobase pairs is shared between the orangutan CP81-1 and the human HLA-DR2 haplotype. These findings indicate that part of the haplotype polymorphism may have persisted for more than 13 million years, which is the estimated time of human-orangutan divergence.

Shared polymorphism between gorilla and human major histocompatibility complex DRB loci

A high degree of polymorphism and high nucleotide diversity mark the functional genes of the major histocompatibility complex (Mbc). Alleles at the different Mbc loci can be classified into distinct lineages that are shared between species and, therefore, are presumed to have been founded before speciation. We have sequenced the most polymorphic part of 25 gorilla Mbc-DRB genes from six individuals. (The DRB genes code for the beta-polypeptide chain of the alpha beta heterodimer that constitutes one family of the class II MHC molecules.) Fifteen of the sequences identify new alleles at four DRB loci; each of the six gorillas was heterozygous at one of the loci at least. Thirteen of the alleles could be assigned to lineages identified previously; the remaining two alleles represent new lineages. All the major human DRB allelic lineages are now known to be shared with apes, and all must have originated before the human-gorilla-chimpanzee divergence more than six million years (my) ago. The presence of some of the gorilla and human lineages in Old World monkeys suggests that these lineages emerged before the divergence of apes and cercopithecids. We argue that the major allelic lineages at the DRB1 locus began to diverge shortly after the rounds of duplication that generated the different DRB loci now found in the hominoids and that this event occurred more than 30 my ago. Comparison of closely related gorilla DRB sequences indicates that polymorphism may be generated by several mechanisms: point mutations, slippage during DNA replication, and recombination. Deduced gene linkages provide evidence for transspecies evolution of haplotype polymorphism.