Mhc-DRB genes of the pigtail macaque (Macaca nemestrina): implications for the evolution of human DRB genes

Comprehensive characterization of MHC class II haplotypes in Mauritian cynomolgus macaques

There are currently no nonhuman primate models with fully defined major histocompatibility complex (MHC) class II genetics. We recently showed that six common MHC haplotypes account for essentially all MHC diversity in cynomolgus macaques (Macaca fascicularis) from the island of Mauritius. In this study, we employ complementary DNA cloning and sequencing to comprehensively characterize full length MHC class II alleles expressed at the Mafa-DPA, -DPB, -DQA, -DQB, -DRA, and -DRB loci on the six common haplotypes. We describe 34 full-length MHC class II alleles, 12 of which are completely novel. Polymorphism was evident at all six loci including DPA, a locus thought to be monomorphic in rhesus macaques. Similar to other Old World monkeys, Mauritian cynomolgus macaques (MCM) share MHC class II allelic lineages with humans at the DQ and DR loci, but not at the DP loci. Additionally, we identified extensive sharing of MHC class II alleles between MCM and other nonhuman primates. The characterization of these full-length-expressed MHC class II alleles will enable researchers to generate MHC class II transferent cell lines, tetramers, and other molecular reagents that can be used to explore CD4+ T lymphocyte responses in MCM.

High levels of diversity characterize mandrill (Mandrillus sphinx) Mhc-DRB sequences

The major histocompatibility complex (MHC) is highly polymorphic in most primate species studied thus far. The rhesus macaque (Macaca mulatta) has been studied extensively and the Mhc-DRB region demonstrates variability similar to humans. The extent of MHC diversity is relatively unknown for other Old World monkeys (OWM), especially among genera other than Macaca. A molecular survey of the Mhc-DRB region in mandrills (Mandrillus sphinx) revealed extensive variability, suggesting that other OWMs may also possess high levels of Mhc-DRB polymorphism. In the present study, 33 Mhc-DRB loci were identified from only 13 animals. Eleven were wild-born and presumed to be unrelated and two were captive-born twins. Two to seven different sequences were identified for each individual, suggesting that some mandrills may have as many as four Mhc-DRB loci on a single haplotype. From these sequences, representatives of at least six Mhc-DRB loci or lineages were identified. As observed in other primates, some new lineages may have arisen through the process of gene conversion. These findings indicate that mandrills have Mhc-DRB diversity not unlike rhesus macaques and humans.

Study of Cynomolgus monkey (Macaca fascicularis) MhcDRB (Mafa-DRB) polymorphism in two populations

Cynomolgus monkey is one of the macaque species currently used as an animal model for experimental surgery and medicine, in particular, to experiment new drugs or therapy protocols designed for the prevention of allograft rejection. In this field, it is of utmost importance to select histoincompatible recipient-donor pairs. One way to ensure incompatibility between donor and recipient is to check their major histocompatibility complex (MHC) genotypes at the loci playing a determinant role in histocompatibility. We report in this paper on the cynomolgus monkey DRB polymorphism evidenced by sequencing of amplified exon 2 separated either by denaturing gradient gel electrophoresis (DGGE), or by cloning. By the study of 253 unrelated animals from two populations (Mauritius and The Philippines), we characterized 50 exon 2 sequences among which 28 were identical to sequences already reported in Macaca fascicularis or other macaque species (Macaca mulatta, Macaca nemestrina). By cloning and sequencing DRB cDNA, we revealed two additional DRB alleles. Out of the 20 haplotypes that we defined here, only two were found in both populations. The functional impact of DR incompatibility was studied in vitro by mixed lymphocyte culture.

Evolution of murine alpha 1-proteinase inhibitors: gene amplification and reactive center divergence

The organization and sequence of genes encoding the alpha 1-proteinase inhibitor (alpha 1PI), a major serine proteinase inhibitor of the mammalian bloodstream, have been compared in several species, including murine rodents (genus Mus). Analysis of gene copy number indicates that amplification of alpha 1PI genes occurred at some time during evolution of the Mus genus, leading to fixation of a family of about three to five genes in several existing species (e.g., M. domesticus and M. saxicola), and only a single gene in others (e.g., M. caroli). A phylogeny for the various mammalian alpha 1PI mRNAs was constructed based upon synonymous substitutions within coding regions. The mRNAs in different murine species diverged from a common ancestor before the formation of the first species lineages of the Mus genus, i.e., about 10-13 million years ago. Thus, alpha 1PI gene amplification must have occurred prior to Mus speciation; gene families were retained in some, but not all, murine species. The reactive center region of the alpha 1PI polypeptide, which determines target protease specificity, has diverged rapidly during evolution of the Mus species, but not during evolution of other mammalian species included in the analysis. It is likely that this accelerated evolution of the reactive center, which has been noted previously for serine proteinase inhibitors, was driven by some sort of a positive Darwinian selection that was exerted in a taxon-specific manner. We suggest that evolution of alpha 1PI genes of murine rodents has been characterized by both modification of gene copy number and rapid reactive center divergence. These processes may have resulted in a broadened repertoire of proteinase inhibitors that was evolutionarily advantageous during Mus speciation.

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.

Exon-intron organization of fish major histocompatibility complex class II B genes

Major histocompatibility complex (Mhc) molecules bind self and foreign peptides and present them to lymphocytes for recognition. Activation of lymphocytes by Mhc-bound foreign peptides leads to specific immune response against parasites. The Mhc genes have been studied extensively in mammals and birds but much less in other vertebrate classes. In this communication we provide the first description of the exon-intron organization of class II beta-chain-encoding genes from the teleost fish Aulonocara hansbaenschi, family Cichlidae. Each of the genes consists of six exons, E1 through E6, encoding the leader peptide (E1), beta 1 domain (E1+E2), beta 2 domain (E3+E4), connecting peptide (E5), transmembrane region (E5), cytoplasmic domain (E5+E6), and the 3' untranslated region (E6). The exons are separated by relatively short introns, the length of the longest intron being 1.3 kilobase pairs. An important difference between these and all other known class II B genes is that the beta 2 domain-encoding exon is split by an intron 97 base pairs in length. The intron is absent in other teleost fishes such as Brachydanio rerio. A change in the 3' splice site of intron 4 in some of the genes of A. hansbaenschi and of another cichlid fish, Cyphotilapia frontosa, has produced two extra codons at the 5' end of exon 5. Comparison of the A. hansbaenschi coding sequences with those of C. frontosa has revealed a concentration of variability in exon 2 and part of exon 3. Taken together, these observations provide evidence for the existence in cichlid fishes of at least two class II B loci which are functionally equivalent to the corresponding loci in mammals. The exon-intron organization and sequence similarities indicate that the two loci arose by duplication from a common ancestral gene.

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.

Cloning of the beta 2-microglobulin gene in the zebrafish

The beta 2-microglobulin (beta 2m) is a protein found in the serum in a free form and on the cell surface in a form noncovalently associated with the alpha chain of the class I major histocompatibility complex (Mhc) molecules. In mammals, the beta 2m-encoding gene (B2m) is found on a chromosome different from the Mhc proper. We have isolated and characterized the B2m gene of the zebrafish, Brachydanio rerio, family Cyprinidae. We obtained both cDNA and genomic clones of the Brre-B2m gene. The cDNA clones contained the entire coding sequence, the entire 3' untranslated (UT) region, and at least part of the 5'UT region. The genomic clone contained the entire Brre-B2m gene. The coding sequence specifies 97 amino acid residues of the mature protein so that the zebrafish beta 2m is two residues shorter than human and one residue shorter than cattle, fowl, or turkey beta 2m (codons at positions 85 and 86 have been deleted in the Brre-B2m gene). The amino acid and nucleotide sequence similarities between zebrafish and human beta 2m (B2m) are 45% and 59%, respectively. Approximately 24% of the positions are invariant and an additional 9% show only conservative substitutions in comparisons which include all known beta 2m sequences (fish, avian, and mammalian). Most of the conserved positions are in the beta strands (some 47% of the beta-strand positions are conserved in the three vertebrate classes). The Brre-B2m gene consists of four exons separated by three introns. All of the introns are considerably shorter than the corresponding introns in the mammalian B2m genes. The coding sequences of the cDNA and the genomic clones are almost identical but the sequences of the 3'UT regions differ at 1.7% of the sites, suggesting that the genes borne by these clones might have diverged at least 0.7 million years (my) ago. In contrast to the human B2m gene, the Brre-B2m gene shows no bias in the distribution of the CpG dinucleotides: the dinucleotides are distributed evenly along the entire available sequence. The haploid genome of the zebrafish contains only one copy of the B2m gene.

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.