Trans-species origin of Mhc-DRB polymorphism in the chimpanzee

Major histocompatibility complex complement (MHC) Bf alleles show trans species evolution between man and chimpanzee

HLA and disease studies by using single allele statistics have been fruitless during the last 40 years for explaining association pathogenesis of the associated diseases.Other approaches are necessary to untangle this puzzle. We aim to revisit complement alleleism in humans and primates for both studying MHC and disease association to complotypes and extended MHC haplotypes in order to also explain the positive directional selection of maintaining immune response genes (complement, MHC adaptive and MHC non-specific genes) that keeps these three type of genes together in a short chromosome stretch (MHC) for million years. These genes may be linked to conjointly avoid microbes attack and autoimmunity. In the present paper, it is obtained a new Bf chimpanzee allele, provisionaly named Patr-Bf*A:01,that differs from other Bf alleles by having CTG at eleventh codon of exon 2 in order to start the newly suggested methodology and explain functional and evolutionary MHC obscure aspects. Exons 1 to 6 of Ba fragment of Bf gene were obtained from chimpanzee. This new chimpanzee Factor B allele (Patr-Bf*A:01) is to be identical to a infrequent human Bf allele (SNP rs641153); it stresses the strong evolutive pressure upon certain alleles that are trans specific. It also may apply to MHC extended haplotipes which may conjointly act to start an adequate immune response. It is the first time that a complement MHC class III allele is described to undergo trans species evolution,in contrast to class I and class II alleles which had already been reported . Allelism of complement factors are again proposed for studying MHC complement genes, complotypes, and extended MHC haplotypes which may be more informative that single MHC marker studies.

Identification of major histocompatibility complex genotypes associated with resistance to an amphibian emerging infectious disease

Amphibian chytridiomycosis, caused by Batrachochytrium dendrobatidis (Bd), emerged from Asia and spread globally. By comparing functional MHC IIß1 alleles from an Asian Bd-resistant anuran species (Bufo gargarizans) with those of an Australasian Bd-susceptible species (Litoria caerulea), we identified MHC genotypes associated with Bd resistance. These alleles encode a glycine deletion (G90β1) and adjacent motifs in the deepest pathogen-derived peptide-binding groove. Every Bd-resistant individual, but no susceptible individuals, possessed at least one allele encoding the variant. We detected trans-species polymorphism at the end of the MHC IIβ1 sequences. The G90β1 deletion was encoded by different alleles in the two species, suggesting it may have evolved independently in each species rather than having been derived from a common ancestor. These results are consistent with a scenario by which MHC adaptations that confer resistance to the pathogen have evolved by convergent evolution. Immunogenetic studies such as this are critical to ongoing conservation efforts.

Trans-specific polymorphism and the convergent evolution of supertypes in major histocompatibility complex class II genes in darters (Etheostoma)

Major Histocompatibility Complex (MHC) genes are one of the most polymorphic gene groups known in vertebrates. MHC genes also exhibit allelic variants that are shared among taxa, referred to as trans-specific polymorphism (TSP). The role that selection plays in maintaining such high diversity within species, as well as TSP, is an ongoing discussion in biology. In this study, we used deep-sequencing techniques to characterize MHC class IIb gene diversity in three sympatric species of darters. We found at least 5 copies of the MHC gene in darters, with 126 genetic variants encoding 122 unique amino acid sequences. We identified four supertypes based on the binding properties of proteins encoded by the sequences. Although each species had a unique pool of variants, many variants were shared between species pairs and across all three species. Phylogenetic analysis showed that the variants did not group together monophyletically based on species identity or on supertype. An expanded phylogenetic analysis showed that some darter alleles grouped together with alleles from other percid fishes. Our findings show that TSP occurs in darters, which suggests that balancing selection is acting at the genotype level. Supertypes, however, are most likely evolving convergently, as evidenced by the fact that alleles do not form monophyletic groups based on supertype. Our research demonstrates that selection may be acting differently on MHC genes at the genotype and supertype levels, selecting for the maintenance of high genotypic diversity while driving the convergent evolution of similar MHC phenotypes across different species.

A human-specific allelic group of the MHC DRB1 gene in primates

Background

Diversity among human leukocyte antigen (HLA) molecules has been maintained by host-pathogen coevolution over a long period of time. Reflecting this diversity, the HLA loci are the most polymorphic in the human genome. One characteristic of HLA diversity is long-term persistence of allelic lineages, which causes trans-species polymorphisms to be shared among closely related species. Modern humans have disseminated across the world after their exodus from Africa, while chimpanzees have remained in Africa since the speciation event between humans and chimpanzees. It is thought that modern humans have recently acquired resistance to novel pathogens outside Africa. In the present study, we investigated HLA alleles that could contribute to this local adaptation in humans and also studied the contribution of natural selection to human evolution by using molecular data.

Results

Phylogenetic analysis of HLA-DRB1 genes identified two major groups, HLA Groups A and B. Group A formed a monophyletic clade distinct from DRB1 alleles in other Catarrhini, suggesting that Group A is a human-specific allelic group. Our estimates of divergence time suggested that seven HLA-DRB1 Group A allelic lineages in humans have been maintained since before the speciation event between humans and chimpanzees, while chimpanzees possess only one DRB1 allelic lineage (Patr-DRB1*03), which is a sister group to Group A. Experimental data showed that some Group A alleles bound to peptides derived from human-specific pathogens. Of the Group A alleles, three exist at high frequencies in several local populations outside Africa.

Conclusions

HLA Group A alleles are likely to have been retained in human lineages for a long period of time and have not expanded since the divergence of humans and chimpanzees. On the other hand, most orthologs of HLA Group A alleles may have been lost in the chimpanzee due to differences in selective pressures. The presence of alleles with high frequency outside of Africa suggests these HLA molecules result from the local adaptations of humans. Our study helps elucidate the mechanism by which the human adaptive immune system has coevolved with pathogens over a long period of time.

The research of W.E. Mayer (1953-2012): a spectrum of immune systems

Over a period of some 20 years, Werner Eugen Mayer played a significant role in establishing a framework for molecular studies of Mhc genes in multiple vertebrates. His work largely concerned gene isolation, sequencing, and related bioinformatic analyses both for the Mhc and for immune system genes of about 200 species, ranging from apes, monkeys, rodents, and marsupials, through to birds, bony fishes, and lampreys. In addition to his exploration of diverse Mhc genes, Werner is remembered for playing a critical role in the development of two important insights into the evolution of immune systems. His was among the first published DNA sequence-based descriptions of trans-species evolution of Mhc alleles, including the first description of the long-lived polymorphisms shared by humans and chimpanzees. This research opened the way for using Mhc polymorphisms in demographic analyses. The second important insight in which he played a prominent role involved the characterization of immune cells and their expressed genes in the lamprey, a jawless vertebrate. His findings helped to indicate the considerable degree to which extant immune mechanisms were co-opted in the creation of the adaptive immune system of jawed vertebrates.

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.

The functional gene diversity in natural populations over postglacial areas: the shaping mechanisms behind genetic composition of longnose dace (Rhinichthys cataractae) in northeastern North America

The diversity of functional genes and the related processes are important issues for conservation biology. This is especially relevant for populations that have suffered from demographic reduction as a consequence of the processes of postglacial colonization. In this perspective, the aims of the present study are (1) to quantify the genetic diversity of functional genes and (2) to disentangle the long- and short-term effects of natural selection that shapes genetic diversity from those of drift, mutation, and allopatric fragmentation. This research was conducted using an extensive genetic polymorphism analysis of populations of longnose dace (Rhinichthys cataractae) living over an area once covered by Pleistocene glaciations. The sequence and diversity of one exon of three genes (MHC IIβ, growth hormone, and trypsin) were jointly analyzed with non-coding nuclear loci from 27 populations; these populations were sampled over four major basins of northeastern North America. The survey revealed a surprisingly low allelic richness, especially for the MHC gene, considering the number of individuals and populations sampled. The results suggest that there is a complex mixture of different evolutionary processes shaping the level of polymorphism among longnose dace. While our study underlines the importance of the short-term effects of neutral processes and the major impact of post-glacial colonization on gene diversity, locally dependent balancing selection was detected on MHC. From this perspective, our results support an understanding of the importance of drift on functional gene diversity but also highlight the transient effects of natural selection on allelic composition, even in populations that show drastic reduction of genetic diversity.

Extensive DRB region diversity in cynomolgus macaques: recombination as a driving force

The DR region of primate species is generally complex and displays diversity concerning the number and combination of distinct types of DRB genes present per region configuration. A highly variable short tandem repeat (STR) present in intron 2 of nearly all primate DRB genes can be utilized as a quick and accurate high through-put typing procedure. This approach resulted previously in the description of unique and haplotype-specific DRB-STR length patterns in humans, chimpanzees, and rhesus macaques. For the present study, a cohort of 230 cynomolgus monkeys, including self-sustaining breeding groups, has been examined. MtDNA analysis showed that most animals originated from the Indonesian islands, but some are derived from the mainland, south and north of the Isthmus of Kra. Haplotyping and subsequent sequencing resulted in the detection of 118 alleles, including 28 unreported ones. A total of 49 Mafa-DRB region configurations were detected, of which 28 have not yet been described. Humans and chimpanzees possess a low number of different DRB region configurations in concert with a high degree of allelic variation. In contrast, however, allelic heterogeneity within a given Mafa-DRB configuration is even less frequently observed than in rhesus macaques. Several of these region configurations appear to have been generated by recombination-like events, most probably propagated by a retroviral element mapping within DRB6 pseudogenes, which are present on the majority of haplotypes. This undocumented high level of DRB region configuration-associated diversity most likely represents a species-specific strategy to cope with various pathogens.

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.

Estimation of the species-specific mutation rates at the DRB1 locus in humans and chimpanzee

To estimate the species-specific mutation rates at the DRB1 locus in humans and chimpanzee, we analyzed the nucleotide sequence of a 37.6-kb chimpanzee chromosomal segment containing the entire Patr-DRB1*0701 allele and the flanking nongenic region and we compared it with two corresponding human sequences containing the HLA-DRB1*070101 allele using the sequence of HLA-DRB1*04011 as an outgroup. Because the allelic pair of HLA-DRB1*070101 and Patr-DRB1*0701 shows the lowest number of substitutions between the two species, it appears that these sequences diverged close to the time of the humans-chimpanzee divergence (6 million years ago). Alignment of the nucleotide sequences for HLA-DRB1*070101 and Patr-DRB1*0701 alleles showed that they share a high degree of similarity, suggesting that the studied chromosomal segments with these sequences have not been subjected to recombination since the humans-chimpanzee divergence. Comparison of the flanking 10.6 kb of nongenic sequences revealed an average of 41.5 and 83 single nucleotide substitutions in humans and chimpanzee, respectively. Thus, the species-specific nucleotide substitution rates in the flanking nongenic region were estimated to be 6.53 x 10(-10) and 1.31 x 10(-9) per site per year in humans and chimpanzee, respectively. Unexpectedly, the estimated rate in humans was twofold lower than in chimpanzee (P < 10(-3), Tajima's relative rate test) and lower than the average substitution rate in the human genome. Because the nucleotide substitution rate in nongenic regions free from selection is expected to be equal to the mutation rate, the estimated substitution rate should correspond to the species-specific mutation rate at the DRB1 locus. Our results strongly suggest that the mutation rate at DRB1 locus differs among species.

Owl monkey MHC-DRB exon 2 reveals high similarity with several HLA-DRB lineages

One hundred and ten novel MHC-DRB gene exon 2 nucleotide sequences were sequenced in 96 monkeys from three owl monkey species (67 from Aotus nancymaae, 30 from Aotus nigriceps and 13 from Aotus vociferans). Owl monkeys, like humans, have high MHC-DRB allele polymorphism, revealing a striking similarity with several human allele lineages in the peptide binding region and presenting major convergence with DRB lineages from several Catarrhini (humans, apes and Old World monkeys) rather than with others New World monkeys (Platyrrhini). The parallelism between human and Aotus MHC-DRB reveals additional similarities regarding variability pattern, selection pressure and physicochemical constraints in amino acid replacements. These observations concerning previous findings of similarity between the Aotus immune system molecules and their human counterparts affirm this specie's usefulness as an excellent animal model in biomedical research.

Identification and analysis of MHC class II DRB1 (Patr-DRB1) alleles in chimpanzees

The MHC-DRB1 gene is known to display the most extensive allelic polymorphisms among MHC class II genes. We attempted the selective identification of chimpanzee (Pan troglodytes) DRB1 (Patr-DRB1) alleles using the polymerase chain reaction (PCR) technique in three steps: first, we performed Patr-DRB1*02 lineage-specific 8-kb PCR for *02 lineage detection in each chimpanzee; second, we performed 620-bp PCR for amplification of full-length exon 2; and finally, we carried out an insert check using the pattern of microsatellite repeat length variability. In the genomic DNA of 23 chimpanzees, nine Patr-DRB1 alleles containing two new alleles were detected. Our approach provides a relatively effective method of identifying Patr-DRB1 alleles in individual chimpanzees and should also contribute to our understanding of the features of MHC molecules in non-human primates.

The evolution of the immunoglobulin heavy chain variable region (IgV H ) in Leporids: an unusual case of transspecies polymorphism

In domestic rabbit (Oryctolagus cuniculus), three serological types have been distinguished at the variable domain of the antibody H chain, the so-called V(H) a allotypes a1, a2, and a3. They correspond to highly divergent allelic lineages of the V(H) 1 gene, which is the gene rabbit utilizes in more than 80% of VDJ rearrangements. The sharing of serological V(H) a markers between rabbit and snowshoe hare (Lepus americanus) has suggested that the large genetic distances between rabbit V(H) 1 alleles (9-14% nucleotide differences) can be explained by unusually long lineage persistence times (transspecies polymorphism). Because this interpretation of the serological data is uncertain, we have determined the nucleotide sequences of V(H) genes expressed in specimens of Lepus species. Two sequence groups were distinguished, one of which occurred only in hare specimen displaying serological motifs of the rabbit V(H) a-a2 allotype. Sequences of this group are part of a monophyletic cluster containing the V(H) 1 sequences of the rabbit a2 allotype. The fact that this "transspecies a2 cluster" did not include genes of other rabbit V(H) a allotypes (a1, a3, and a4) is incompatible with the existence of a common V(H) a ancestor gene within the species, and suggests that the divergence of the V(H) a lineages preceded the Lepus vs Oryctolagus split. The sequence data are furthermore compatible with the hypothesis that the V(H)a polymorphism can be two times older than the divergence time between the Lepus and Oryctolagus lineages, which was estimated at 16-24 million years.

Independent origin of functional MHC class II genes in humans and New World monkeys

In previous studies, major histocompatibility complex (MHC) class II DP, DQ, and DR families of genes were characterized in different primate species mostly on the basis of their second exon sequences. Resemblances were found between Old World monkey (OWM) and New World monkey (NWM) genes and were interpreted as being the result of transspecies evolution. Subsequent analysis of intron sequences of catarrhine and platyrrhine DRB genes, however, revealed that the amplifiable genes were not, in fact, orthologous. To test other DRB genes and other families of the class II region Southern blot hybridizations were carried out with tamarin genomic DNA using probes specific for the third exons of the tamarin DQA, DQB, DPB, and DRB genes. The hybridizing bands were extracted from the gel and the third exons of the genes were amplified by PCR, cloned, and sequenced. With two exceptions, all NWM class II genes were found to group separately from the human sequences. Only the sequences of one nonfunctional DQB locus appeared to be more closely related to human genes than to other platyrrhine DQB genes. In the DRB family one gene was found that grouped with sheep and strepsirhine DRB sequences and might represent an old gene lineage. To extend the sequences to the second exon, long PCRs were performed on tamarin genomic DNA. This approach was successful for five of the ten third exon sequences. From these data, we conclude that at least the functional MHC class II genes have expanded independently in catarrhines and platyrrhines.

Trans-speciation maintenance in the MHC region of a polymorphism which includes a polymorphic dinucleotide locus, and the de novo arisal of a polymorphic tetranucleotide microsatellite

Alleles and the surrounding regions of DQCAR, a dinucleotide repeat tightly linked to HLA-DQB1, were sequenced in a range of primate species including man. Three polymorphic regions can usefully be defined in the description of these sequences: the dinucleotide GT repeat itself, the anonymous region 5' of this repeat, and a variable CTGT repeat in the 3' region. The 5' sequence displayed six alleles in the individuals studied. One of these alleles was invariably associated with substitutions in the GT repeat and absence of the CTGT repeat, the others with pure, polymorphic GT repeats and variation in the numbers of CTGT repeats. Haplotypes can be classified by the allele in the 5' region. Those carrying allele 1 were only found in man, those with allele 2 in man, chimpanzee and gorilla. The third haplotype (indicated by the presence of allele 3) was found in chimpanzee, gorilla and orang-utan, the fourth in chimpanzee and gibbon, the fifth in baboon, guenon and mangabey and the sixth in guenon and macaque. The alleles in the 5' region, but from different species, are thus often more similar than alleles from the same species, a phenomenon already shown for some HLA genes. This suggests that major histocompatibility sequences and surrounding sequences shared a correlated evolutionary history. The new polymorphic tetranucleotide microsatellite (CTGT, 3rd region) has possibly arisen de novo from the pre-existing dinucleotide GT. This study provides information not only on the molecular evolution of this particular microsatellite but also of the trans-speciation maintenance of polymorphism of its surrounding sequences.

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.

Tracing the origin of HLA-DRB1 alleles by microsatellite polymorphism

We analyzed the origin of allelic diversity at the class II HLA-DRB1 locus, using a complex microsatellite located in intron 2, close to the polymorphic second exon. A phylogenetic analysis of human, gorilla, and chimpanzee DRB1 sequences indicated that the structure of the microsatellite has evolved, primarily by point mutations, from a putative ancestral (GT)x(GA)y-complex-dinucleotide repeat. In all contemporary DRB1 allelic lineages, with the exception of the human *04 and the gorilla *08 lineages, the (GA)y repeat is interrupted, often by a G-->C substitution. In general, the length of the 3' (GA)y repeat correlates with the allelic lineage and thus evolves more slowly than a middle (GA)z repeat, whose length correlates with specific alleles within the lineage. Comparison of the microsatellite sequence from 30 human DRB1 alleles showed the longer 5' (GT)x to be more variable than the shorter middle (GA)z and 3' (GA)y repeats. Analysis of multiple samples with the same exon sequence, derived from different continents, showed that the 5' (GT)x repeat evolves more rapidly than the middle (GA)z and the 3' (GA)y repeats, which is consistent with findings of a higher mutation rate for longer tracts. The microsatellite-repeat-length variation was used to trace the origin of new DRB1 alleles, such as the new *08 alleles found in the Cayapa people of Ecuador and the Ticuna people of Brazil.

Major histocompatibility complex class II polymorphisms in primates

In the past decade, the major histocompatibility complex (MHC) class II region of several primate species has been investigated extensively. Here we will discuss the similarities and differences found in the MHC class II repertoires of primate species including humans, chimpanzees, rhesus macaques, cotton-top tamarins and common marmosets. Such types of comparisons shed light on the evolutionary stability of MHC class II alleles, lineages and loci as well as on the evolutionary origin and biological significance of haplotype configurations.

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.

Comparative genome organization of the major histocompatibility complex: lessons from the Felidae

The mammalian major histocompatibility complex (MHC) has taught both immunologists and evolutionary biologists a great deal about the patterns and processes that have led to immune defenses. Driven principally by human and mouse studies, comparative MHC projects among other mammalian species offer certain advantages in connecting MHC genome characters to natural situations. We have studied the MHC in the domestic cat and in several wild species of Felidae. Our observations affirm class I and class II homology with other mammalian orders, derivative gene duplications during the Felidae radiation, abundant persistent trans-species allele polymorphism, recombination-derived amino acid motifs, and inverted ratios of non-synonymous to silent substitutions in the MHC peptide-binding regions, consistent with overdominant selection in class I and II genes. MHC diversity as quantified in population studies is a powerful barometer of historic demographic reduction for several endangered species including cheetahs, Asiatic lions, Florida panthers and tigers. In two cases (Florida panther and cheetah), reduced MHC variation may be contributing to uniform population sensitivity to emerging infectious pathogens. The Felidae species, nearly all endangered and monitored for conservation concerns, have allowed a glimpse of species adaptation, mediated by MHC divergence, using comparative inferences drawn from human and mouse models.

Analysis of intron sequences at the class II HLA-DRB1 locus: implications for the age of allelic diversity

Analyses of the coding sequences of HLA class II alleles have revealed high similarity between species, indicating that much of the polymorphism predates the separation of human (Homo) and chimpanzee (Pan), 4-7.4 million years ago. Recent studies of the intron sequences of alleles provide support for a much more recent origin and rapid generation of HLA alleles. At the DRB1 locus, intron analysis indicates that most of the allelic lineages have diverged from each other before the separation of Homo and Pan, consistent with the exon analysis. However, the intron sequences of alleles within lineages are almost identical, indicating that many of the alleles have been generated after the divergence of the Homo and Pan lineages.

Mhc diversity in Pacific salmon: population structure and trans-species allelism

Geographic variation at an Mhc class I A1 exon was surveyed in 14 populations of coho salmon (Oncorhynchus kisutch) and 15 populations of chinook salmon (O. tshawytscha) inhabiting rivers of British Columbia, Canada. A total of 2,504 fish were sampled using denaturing gradient gel electrophoresis (DGGE), which distinguished 17 alleles in coho salmon and 20 alleles in chinook salmon. Heterozygosity at the A1 locus was moderately high for both coho (0.7) and chinook (0.6) salmon, but sequence divergence was low, with mean inter- and intraspecific nucleotide similarities of approximately 0.96. In a maximum parsimony tree, all of the observed alleles clustered into two trans-specific lineages. Within each lineage, coho and chinook alleles tended to fall into species-specific subclusters. Much of the intraspecific allelic variation within each lineage could be accounted for by nonsynonymous point mutation, indicative of balancing selection. The FST values for both coho (0.11) and chinook (0.13) salmon indicated that much of the allelic diversity was partitioned among populations. Neighbor-joining analyses of A1 allelic frequencies among coho and chinook salmon populations showed strong patterns of geographic differentiation similar to those based on neutral genetic markers such as microsatellite loci. Both natural selection and the salmonid zoogeographic history of frequent population bottlenecks have shaped the patterns of diversity observed at this and other Mhc exons in Pacific salmonids.

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.

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.

Trans-species polymorphism of class II Mhc loci in danio fishes

A characteristic feature of the major histocompatibility complex (Mhc) polymorphism in mammals is the existence of allelic lineages shared by related species. This trans-species polymorphism has thus far been documented only in primates, rodents, and artiodactyls. In this communication we provide evidence that it also exists in cyprinid (bony) fishes at the class II A and B loci coding for the alpha and beta polypeptide chains of the class II alpha:beta heterodimers. The study has focused on three species of the family Cyprinidae, subfamily Rasborinae: the zebrafish (Danio rerio), the giant danio (D. malabaricus), and the pearl danio (D. albolineatus). The polymerase chain reaction was used to amplify and then sequence intron 1 and exon 2 of the class II B loci and exon 2 of the class II A loci in these species. Phylogenetic analysis of the sequences revealed the existence of allelic lineages whose divergence predates the divergence of the three species at both the A and B loci. The lineages at the B locus in particular are separated by large genetic distances. The polymorphism is concentrated in the peptide-binding region sites and is apparently maintained by balancing selection. Sharing of this unique Mhc feature by both bony fishes and mammals suggests that the main function of the Mhc (presentation of peptides to T lymphocytes) has not changed during the last 400 million years of its evolution.

Evolution of Mhc class II polymorphism: the rise and fall of class II gene function in primates

The substitution rate at the codons implicated at ARS of Mhc class II genes has previously been shown to be heavily biased towards nonsynonymous substitutions, indicative of positive selection for polymorphism. Based on our analysis of the number of synonymous changes at codons outside putative ARS in primates, the average age of the polymorphism at class II loci was found to increase in the following order: DPB1, DRB3, DRB5, DRB1, DRB4, DQB1, DQA1. For DRB loci, nonsynonymous changes were found to exceed synonymous changes at HLA-DRB1, DRB3 and DRB5, while no evidence of deviations from equal rates of synonymous and nonsynonymous substitutions were found for DRB6. The pattern of substitutions at the DRB loci of most Catarrhini species indicates constant positive selection at ARS codons over the evolutionary period examined. An exception to the relatively stable selection pattern between species exhibited by most loci is the appearance of polymorphism under positive selection at DRB4 only in the regular chimpanzee. The ds/dn ratios for DQA1 and DQB1 alleles are lower than for the most polymorphic DRB genes. Since the dn/ds ratio of ARS codons may be positively correlated to the ds for non-ARS codons, at least for DQB1, caution must be exercised in interpreting the low ratio for the DQ genes as an indication of weaker selection. The DQA1 allelic lineages show different dn/ds ratios, consistent with the hypothesis that the lineages are constrained from evolving in relation to the diversity of the interacting DQB1 alleles. In contrast to all other class II loci, DPB1 appears to have been subjected to strong positive selection only in the human lineage, and may represent the most conspicuous example of an Mhc locus acquiring an altered function in antigen presentation.

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.

A uniquely high level of recombination at the HLA-B locus

Major histocompatibility complex (MHC) loci are some of the most polymorphic genes in the animal kingdom. Recently, it has been suggested that although most of the human MHC loci are relatively stable, the HLA-B locus can undergo rapid changes, especially in isolated populations. To investigate the mechanisms of HLA-B evolution we have compared the sequences of 19 HLA-B homologues from chimpanzees and bonobos to 65 HLA-B sequences. Analysis of the chimpanzee and bonobo HLA-B homologues revealed that despite obvious similarities between chimpanzee and human alleles in exon 2, there was little conservation of exon 3 between humans and the two chimpanzee species. This finding suggests that, unlike all other HLA loci, recombination has characterized the HLA-B locus and its homologues for over 5 million years.

Gel electrophoretic analysis of rhesus macaque major histocompatibility complex class II DR molecules

Rhesus macaque MHC class II DR molecules were isolated from radiolabeled B-cell line extracts by immunoprecipitation with the mAbs 7.3.19.1 and B8.11.2 and subsequently analyzed by 2D-gel electrophoresis. The B-cell lines used for this study were obtained from monkeys that are homozygous for the Mamu-DR region as defined by serologic techniques. Some of these animals have been selectively bred and originate from consanguineous matings. These analyses show that monkeys with the same allotyping may express different types of DR molecules. As in humans, the number of DR molecules expressed per haplotype is not constant and varies from 1 to 3, depending on the serologically defined Mamu-DR specificity, whereas it has been shown that the number of Mamu-DRB genes present per haplotype varies from 2 to 6. Therefore the present study also demonstrates that some of the rhesus macaque DR regions contain one or more pseudogenes.

Resolution of the HLA-DRB6 puzzle: a case of grafting a de novo-generated exon on an existing gene

HLA-DRB6, one of the human major histocompatibility complex genes, lacks exon 1, which normally codes for the leader and the first four amino acid residues of the mature protein. Because it also lacks the HLA promoter, it was surprising to find that the gene is transcribed at a low level in a chimpanzee B-lymphoblastoma cell line, in which the DRB6 homolog is truncated as in humans. The study designed to resolve the paradox has revealed that a retrovirus related to the mouse mammary tumor viruses was inserted into intron 1 of DRB6 > 23 million years ago. The insertion was either accompanied or followed by the deletion of exon 1 and the promoter region of DRB6. In the 3' long terminal repeat of the retrovirus, however, an open reading frame for a new exon arose, which codes for a sequence of mostly hydrophobic amino acid residues; the sequence could function as a leader for the truncated DRB6 gene. This new exon has a functional donor splice site at its 3' end, which enables it to be spliced in register with DRB6 exon 2. Upstream from the new exon is a promoter enabling transcription of the DRB6 gene. Besides providing an example of a de novo generation of an exon, the study suggests a potential mechanism for generating new genes through the replacement of old exons with newly generated ones.

The biologic importance of conserved major histocompatibility complex class II motifs in primates

Phylogenetic comparisons of polymorphic second-exon sequences of MHC class II DRB genes showed that equivalents of the HLA-DRB1*03 alleles are present in various nonhuman primate species such as chimpanzees, gorillas, and rhesus macaques. These alleles must root from ancestral structure(s) that were once present in a progenitor species that lived about 35 million years ago. Due to accumulation of genetic variation, however, sequences that cluster into a lineage are generally unique to a species. To investigate the biologic importance of such conservation and variation, the peptide-binding capacity of various Mhc-DRB1*03 lineage members was studied. Primate Mhc-DRB1*03 lineage members successfully binding the p3-13 peptide of the 65-kD heat-shock protein of Mycobacterium tuberculosis/leprae share a motif that maps to the floor of the peptide-binding site. Apart from that, some rhesus macaque MHC class-II-positive cells were able to present the p3-13 peptide to HLA-DR17-restricted T cells whereas cells obtained from great ape species failed to do so. Therefore, these studies open ways to understand which MHC polymorphisms have been maintained in evolution and which MHC residues are essential for peptide binding and T-cell recognition.

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.