Analysis of macaque BTN3A genes and transcripts in the extended MHC: conserved orthologs of human γδ T cell modulators

Butyrophilins (BTN), specifically BTN3A, play a central role in the modulation of γδ T cells, which are mainly present in gut and mucosal tissues. BTN3A1 is known, for example, to activate Vγ9Vδ2 T cells by means of a phosphoantigen interaction. In the extended HLA region, three genes are located, designated BTN3A1, BTN3A2 and BTN3A3, which were also defined in rhesus macaques. In contrast to humans, rhesus monkeys have an additional gene, BTN3A3Like, which has the features of a pseudogene. cDNA analysis of 32 Indian rhesus and 16 cynomolgus macaques originating from multiple-generation families revealed that all three genes are oligomorphic, and the deduced amino acids display limited variation. The macaque BTN3A alleles segregated together with MHC alleles, proving their location in the extended (Major Histocompatibility Complex) MHC. BTN3A nearly full-length transcripts of macaques and humans cluster tightly together in the phylogenetic tree, suggesting that the genes represent true orthologs of each other. Despite the limited level of polymorphism, 15 Mamu- and 14 Mafa-BTN3A haplotypes were defined, and, as in humans, all three BTN3A genes are transcribed in PBMCs and colon tissues. In addition to regular full-length transcripts, a high number of various alternative splicing (AS) products were observed for all BTN3A alleles, which may result in different isoforms. The comparable function of certain subsets of γδ T cells in human and non-human primates in concert with high levels of sequence conservation observed for the BTN3A transcripts presents the opportunity to study these not yet well understood molecules in macaques as a model species.

Determining Mhc-DRB profiles in wild populations of three congeneric true lemur species by noninvasive methods

The major histocompatibility complex (MHC) is a highly polymorphic and polygenic genomic region that plays a crucial role in immune-related diseases. Given the need for comparative studies on the variability of immunologically important genes among wild populations and species, we investigated the allelic variation of MHC class II DRB among three congeneric true lemur species: the red-fronted lemur (Eulemur rufifrons), red-bellied lemur (Eulemur rubriventer), and black lemur (Eulemur macaco). We noninvasively collected hair and faecal samples from these species across different regions in Madagascar. We assessed DRB exon 2 polymorphism with a newly developed primer set, amplifying nearly all non-synonymous codons of the antigen-binding sites. We defined 26 DRB alleles from 45 individuals (17 alleles from E. rufifrons (N = 18); 5 from E. rubriventer (N = 7); and 4 from E. macaco (N = 20). All detected alleles are novel and show high levels of nucleotide (26.8%) and non-synonymous codon polymorphism (39.4%). In these lemur species, we found neither evidence of a duplication of DRB genes nor a sharing of alleles among sympatric groups or allopatric populations of the same species. The non-sharing of alleles may be the result of a geographical separation over a long time span and/or different pathogen selection pressures. We found dN/dS rates > 1 in the functionally important antigen recognition sites, providing evidence for balancing selection. Especially for small and isolated populations, quantifying and monitoring DRB variation are recommended to establish successful conservation plans that mitigate the possible loss of immunogenetic diversity in lemurs.

A Specialist Macaque MHC Class I Molecule with HLA-B*27-like Peptide-Binding Characteristics

In different macaque species, the MHC A2*05 gene is present in abundance, and its gene products are characterized by low cell-surface expression and a highly conserved peptide-binding cleft. We have characterized the peptide-binding motif of Mamu-A2*05:01, and elucidated the binding capacity for virus-derived peptides. The macaque A2*05 allotype prefers the basic amino acid arginine at the second position of the peptide, and hydrophobic and polar amino acids at the C-terminal end. These preferences are shared with HLA-B*27 and Mamu-B*008, molecules shown to be involved in elite control in human HIV type 1 and macaque SIV infections, respectively. In contrast, however, Mamu-A2*05 preferentially binds 8-mer peptides. Retention in the endoplasmic reticulum seems to be the cause of the lower cell-surface expression. Subsequent peptide-binding studies have illustrated that Mamu-A2*05:01 is able to bind SIV-epitopes known to evoke a strong CD8⁺ T cell response in the context of the Mamu-B*008 allotype in SIV-infected rhesus macaques. Thus, the macaque A2*05 gene encodes a specialized MHC class I molecule, and is most likely transported to the cell surface only when suitable peptides become available.

No postcopulatory selection against MHC-homozygous offspring: Evidence from a pedigreed captive rhesus macaque colony

The heterozygosity status of polymorphic elements of the immune system, such as the major histocompatibility complex (MHC), is known to increase the potential to cope with a wider variety of pathogens. Pre- and postcopulatory processes may regulate MHC heterozygosity. In a population where mating occurs among individuals that share identical MHC haplotypes, postcopulatory selection may disfavour homozygous offspring or ones with two MHC haplotypes identical to its mother. We tested these ideas by determining the incidence of MHC-heterozygous and MHC-homozygous individuals in a pedigreed, partially consanguineous captive rhesus monkey colony. Bayesian statistics showed that when parents share MHC haplotypes, the distribution of MHC-heterozygous and MHC-homozygous individuals significantly fitted the expected Mendelian distribution, both for the complete MHC haplotypes, and for MHC class I or II genes separately. Altogether, we found in this captive colony no evidence for postcopulatory selection against MHC-homozygous individuals. However, the distribution of paternally and maternally inherited MHC haplotypes tended to differ significantly from expected. Individuals with two MHC haplotypes identical to their mother were underrepresented and offspring with MHC haplotypes identical to their father tended to be overrepresented. This suggests that postcopulatory processes affect MHC haplotype combination in offspring, but do not prevent low MHC heterozygosity.

Major histocompatibility complex haplotyping and long-amplicon allele discovery in cynomolgus macaques from Chinese breeding facilities

Very little is currently known about the major histocompatibility complex (MHC) region of cynomolgus macaques (Macaca fascicularis; Mafa) from Chinese breeding centers. We performed comprehensive MHC class I haplotype analysis of 100 cynomolgus macaques from two different centers, with animals from different reported original geographic origins (Vietnamese, Cambodian, and Cambodian/Indonesian mixed-origin). Many of the samples were of known relation to each other (sire, dam, and progeny sets), making it possible to characterize lineage-level haplotypes in these animals. We identified 52 Mafa-A and 74 Mafa-B haplotypes in this cohort, many of which were restricted to specific sample origins. We also characterized full-length MHC class I transcripts using Pacific Biosciences (PacBio) RS II single-molecule real-time (SMRT) sequencing. This technology allows for complete read-through of unfragmented MHC class I transcripts (~1100 bp in length), so no assembly is required to unambiguously resolve novel full-length sequences. Overall, we identified 311 total full-length transcripts in a subset of 72 cynomolgus macaques from these Chinese breeding facilities; 130 of these sequences were novel and an additional 115 extended existing short database sequences to span the complete open reading frame. This significantly expands the number of Mafa-A, Mafa-B, and Mafa-I full-length alleles in the official cynomolgus macaque MHC class I database. The PacBio technique described here represents a general method for full-length allele discovery and genotyping that can be extended to other complex immune loci such as MHC class II, killer immunoglobulin-like receptors, and Fc gamma receptors.

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

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

The orthologs of HLA-DQ and -DP genes display abundant levels of variability in macaque species

The human major histocompatibility complex (MHC) region encodes three types of class II molecules designated HLA-DR, -DQ, and -DP. Both the HLA-DQ and -DP gene region comprise a duplicated tandem of A and B genes, whereas in macaques, only one set of genes is present per region. A substantial sequencing project on the DQ and DP genes in various macaque populations resulted in the detection of previously 304 unreported full-length alleles. Phylogenetic studies showed that humans and macaques share trans-species lineages for the DQA1 and DQB1 genes, whereas the DPA1 and DPB1 lineages in macaques appear to be species-specific. Amino acid variability plot analyses revealed that each of the four genes displays more allelic variation in macaques than is encountered in humans. Moreover, the numbers of different amino acids at certain positions in the encoded proteins are higher than in humans. This phenomenon is remarkably prominent at the contact positions of the peptide-binding sites of the deduced macaque DPβ-chains. These differences in the MHC class II DP regions of macaques and humans suggest separate evolutionary mechanisms in the generation of diversity.

Fifty-one full-length major histocompatibility complex class II alleles in the olive baboon (Papio anubis)

Here we report 51 novel major histocompatibility complex (MHC) class II alleles in a group of related olive baboons.

Co-evolution of the MHC class I and KIR gene families in rhesus macaques: ancestry and plasticity

Researchers dealing with the human leukocyte antigen (HLA) class I and killer immunoglobulin receptor (KIR) multi‐gene families in humans are often wary of the complex and seemingly different situation that is encountered regarding these gene families in Old World monkeys. For the sake of comparison, the well‐defined and thoroughly studied situation in humans has been taken as a reference. In macaques, both the major histocompatibility complex class I and KIR gene families are plastic entities that have experienced various rounds of expansion, contraction, and subsequent recombination processes. As a consequence, haplotypes in macaques display substantial diversity with regard to gene copy number variation. Additionally, for both multi‐gene families, differential levels of polymorphism (allelic variation), and expression are observed as well. A comparative genetic approach has allowed us to answer questions related to ancestry, to shed light on unique adaptations of the species’ immune system, and to provide insights into the genetic events and selective pressures that have shaped the range of these gene families.

Complex MHC Class I Gene Transcription Profiles and Their Functional Impact in Orangutans

MHC haplotypes of humans and the African great ape species have one copy of the MHC-A, -B, and -C genes. In contrast, MHC haplotypes of orangutans, the Asian great ape species, exhibit variation in the number of gene copies. An in-depth analysis of the MHC class I gene repertoire in the two orangutan species, Pongo abelii and Pongo pygmaeus, is presented in this article. This analysis involved Sanger and next-generation sequencing methodologies, revealing diverse and complicated transcription profiles for orangutan MHC-A, -B, and -C. Thirty-five previously unreported MHC class I alleles are described. The data demonstrate that each orangutan MHC haplotype has one copy of the MHC-A gene, and that the MHC-B region has been subject to duplication, giving rise to at least three MHC-B genes. The MHC-B*03 and -B*08 lineages of alleles each account for a separate MHC-B gene. All MHC-B*08 allotypes have the C1-epitope motif recognized by killer cell Ig-like receptor. At least one other MHC-B gene is present, pointing to MHC-B alleles that are not B*03 or B*08. The MHC-C gene is present only on some haplotypes, and each MHC-C allotype has the C1-epitope. The transcription profiles demonstrate that MHC-A alleles are highly transcribed, whereas MHC-C alleles, when present, are transcribed at very low levels. The MHC-B alleles are transcribed to a variable extent and over a wide range. For those orangutan MHC class I allotypes that are detected by human monoclonal anti-HLA class I Abs, the level of cell-surface expression of proteins correlates with the level of transcription of the allele.

Novel DRA alleles extracted from seven macaque cohorts

In this document, we report the detection of 37 DRA alleles in macaque cohorts.

Differential recombination dynamics within the MHC of macaque species

A panel of 15 carefully selected microsatellites (short tandem repeats, STRs) has allowed us to study segregation and haplotype stability in various macaque species. The STRs span the major histocompatibility complex (MHC) region and map in more detail from the centromeric part of the Mhc-A to the DR region. Two large panels of Indian rhesus and Indonesian/Indochinese cynomolgus macaques have been subjected to pedigree analysis, allowing the definition of 161 and 36 different haplotypes and the physical mapping of 10 and 5 recombination sites, respectively. Although most recombination sites within the studied section of the Indian rhesus monkey MHC are situated between the Mhc-A and Mhc-B regions, the resulting recombination rate for this genomic segment is low and similar to that in humans. In contrast, in Indonesian/Indochinese macaques, two recombination sites, which appear to be absent in rhesus macaques, map between the class III and II regions. As a result, the mean recombination frequency of the core MHC, Mhc-A to class II, is higher in Indonesian/Indochinese cynomolgus than in Indian rhesus macaques, but as such is comparable to that in humans. The present communication demonstrates that the dynamics of recombination 'hot/cold spots' in the MHC, as well as their frequencies, may differ substantially between highly related macaque species.

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.

Multilocus definition of MHC haplotypes in pedigreed cynomolgus macaques (Macaca fascicularis)

Cynomolgus macaques (Macaca fascicularis) are used widely in biomedical research, and the genetics of their MHC (Mhc-Mafa) has become the focus of considerable attention in recent years. The cohort of Indonesian pedigreed macaques that we present here was typed for Mafa-A, -B, and -DR, by sequencing, as described in earlier studies. Additionally, the DRB region of these animals was characterised by microsatellite analyses. In this study, full-length sequencing of Mafa-DPA/B and -DQA/B in these animals was performed. A total of 75 different alleles were observed; 22 of which have not previously been reported, plus 18 extended exon 2 alleles that were already known. Furthermore, two microsatellites, D6S2854 and D6S2859, were used to characterise the complex Mafa-A region. Sequencing and segregation analyses revealed that the length patterns of these microsatellites are unique for each Mafa-A haplotype. In this work, we present a pedigreed colony of approximately 120 cynomolgus macaques; all of which are typed for the most significant polymorphic MHC class I and class II markers. Offspring of these pedigreed animals are easily characterised for their MHC by microsatellite analyses on the Mafa-A and -DRB regions, which makes the cumbersome sequencing analyses redundant.

Male Mating Tactics in Captive Rhesus Macaques (Macaca mulatta): The Influence of Dominance, Markets, and Relationship Quality

Male mating success in a multimale-multifemale group can depend on several variables: body condition, dominance, coalitions, "friendship," or an exchange of services for mating access. Exchange patterns may also be determined by market effects or social relationships. We studied the mating tactics of males in a captive, multimale-multifemale group of rhesus macaques and the resulting patterns of mating and paternity to determine the influence of dominance rank, mating markets, and relationship quality on their mating tactics. Male rank was positively related to the total number of copulations and the number of mating partners, but did not explain male mating distribution completely. Moreover, male fertilization success was not related to male rank. Males did not exchange grooming for mating access on the same day and neither the supply nor the rank (as a proxy for quality) of receptive females affected the amount of male grooming, suggesting that market effects did not explain male mating access. However, there was a positive correlation between long-term grooming patterns of both males and females and mating access, indicating that social relationships were important for male mating access. Paternity data revealed that these social relationships were also important for male reproductive success. We conclude that both male rank and male-female "friendship" determined male mating access in these rhesus macaques, but that "friendship" was more important in determining paternity, emphasizing the importance of intersex social bonds in male mating success in multimale primate societies. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1007/s10764-011-9552-5) contains supplementary material, which is available to authorized users.

The mosaic of KIR haplotypes in rhesus macaques

To further refine and improve biomedical research in rhesus macaques, it is necessary to increase our knowledge concerning both the degree of allelic variation (polymorphism) and diversity (gene copy number variation) in the killer cell immunoglobulin-like receptor (KIR) gene cluster. Pedigreed animals in particular should be studied, as segregation data will provide clues to the linkage of particular KIR genes/alleles segregating on a haplotype and to its gene content as well. A dual strategy allowed us to screen the presence and absence of genes and the corresponding transcripts, as well as to track differences in transcription levels. On the basis of this approach, 14 diverse KIR haplotypes have been described. These haplotypes consist of multiple inhibitory and activating Mamu-KIR genes, and any gene present on one haplotype may be absent on another. This suggests that the cost of accelerated evolution by recombination may be the loss of certain framework genes on a haplotype.

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.

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

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

Evidence for balancing selection acting on KIR2DL4 genotypes in rhesus macaques of Indian origin

The interaction of killer-cell immunoglobulin-like receptors (KIR) and their respective major histocompatibility complex (MHC) ligands can alter the activation state of the natural killer (NK) cell. In both humans and rhesus macaques, particular types of non-classical MHC class I molecules are predominantly expressed on the trophoblast. In humans, human leukocyte antigen G has been demonstrated to act as a ligand for KIR2DL4, present on all NK cells, whereas Mamu-AG may execute a similar function in rhesus macaques. During primate evolution, orthologues of KIR2DL4 appear to have been highly conserved, suggesting strong purifying selection. A cohort of 112 related and unrelated rhesus macaques of mostly Indian origin were selected to study their KIR2DL4 genes for the occurrence of polymorphism. Comparison of the proximal region provided evidence for strong conservative selection acting on the exons encoding the Ig domains. As is found in humans, in the Indian rhesus macaque population, two different KIR2DL4 entities are encountered, which differ for their intra-cellular signalling motifs. One genotype contains a complex mutation in the distal region of exon 9, which negates a serine/threonine kinase site. Furthermore, both allelic entities are present in a distribution, which suggests that balancing selection is operating on these two distinct forms of KIR2DL4.

Compound evolutionary history of the rhesus macaque MHC class I B region revealed by microsatellite analysis and localization of retroviral sequences

In humans, the single polymorphic B locus of the major histocompatibility complex is linked to the microsatellite MIB. In rhesus macaques, however, haplotypes are characterized by the presence of unique combinations of multiple B genes, which may display different levels of polymorphism. The aim of the study was to shed light on the evolutionary history of this highly complex region. First, the robustness of the microsatellite MIB-linked to almost half of the B genes in rhesus macaques (Mamu-B)–for accurate B haplotyping was studied. Based on the physical map of an established haplotype comprising 7 MIB loci, each located next to a certain Mamu-B gene, two MIB loci, MIB1 and MIB6, were investigated in a panel of MHC homozygous monkeys. MIB1 revealed a complex genotyping pattern, whereas MIB6 analysis resulted in the detection of one or no amplicon. Both patterns are specific for a given B haplotype, show Mendelian segregation, and even allow a more precise haplotype definition than do traditional typing methods. Second, a search was performed for retroelements that may have played a role in duplication processes as observed in the macaque B region. This resulted in the description of two types of duplicons. One basic unit comprises an expressed Mamu-B gene, adjacent to an HERV16 copy closely linked to MIB. The second type of duplicon comprises a Mamu-B (pseudo)gene, linked to a truncated HERV16 structure lacking its MIB segment. Such truncation seems to coincide with the loss of B gene transcription. Subsequent to the duplication processes, recombination between MIB and Mamu-B loci appears to have occurred, resulting in a hyperplastic B region. Thus, analysis of MIB in addition to B loci allows deciphering of the compound evolutionary history of the class I B region in Old World monkeys.

A splice site mutation converts an inhibitory killer cell Ig-like receptor into an activating one

The killer cell Ig-like receptor (KIR) 3DH protein in rhesus macaques (Macaca mulatta) is thought to be an activating one because it contains a charged arginine in its transmembrane domain and has a truncated cytoplasmic domain. MmKIR3DH has thus far been characterized by an analysis of cDNA. Its presence and polymorphism has been further investigated by examining mRNA transcripts and genomic sequences in families. Multiple copies of MmKIR3DH are present per animal, suggesting that the gene has been duplicated on some haplotypes. All transcripts are truncated and lack exon 8. Investigation of the gene itself shows that exon 8 is present, intact, and homologous to MmKIR2DL4. However, there is a mutation in the donor splice site of intron 8, which is absent in MmKIR2DL4 genomic sequences. This mutation introduces a frameshift, subsequently resulting in a premature stopcodon. To further verify this mutation, a cohort of unrelated animals from different geographical locations was examined, and both exon 8 and the splice site mutation were seen to be present in their MmKIR3DH genes. The data suggest that the splice site mutation causes the truncation of the MmKIR3DH transcript and the subsequent loss of its inhibitory motifs further downstream. Loss of inhibitory potential through different mutations is observed in other primate species as well, suggesting convergent evolution; however, this is the first report to document that a mutation in an intron produces a similar effect.

Genomic plasticity of the immune-related Mhc class I B region in macaque species

Background

In sharp contrast to humans and great apes, the expanded Mhc-B region of rhesus and cynomolgus macaques is characterized by the presence of differential numbers and unique combinations of polymorphic class I B genes per haplotype. The MIB microsatellite is closely linked to the single class I B gene in human and in some great apes studied. The physical map of the Mhc of a heterozygous rhesus monkey provides unique material to analyze MIB and Mamu-B copy number variation and then allows one to decipher the compound evolutionary history of this region in primate species.

Results

In silico research pinpointed 12 MIB copies (duplicons), most of which are associated with expressed B-genes that cluster in a separate clade in the phylogenetic tree. Generic primers tested on homozygous rhesus and pedigreed cynomolgus macaques allowed the identification of eight to eleven MIB copies per individual. The number of MIB copies present per haplotype varies from a minimum of three to six in cynomolgus macaques and from five to eight copies in rhesus macaques. Phylogenetic analyses highlight a strong transpecific sharing of MIB duplicons. Using the physical map, we observed that, similar to MIB duplicons, highly divergent Mamu-B genes can be present on the same haplotype. Haplotype variation as reflected by the copy number variation of class I B loci is best explained by recombination events, which are found to occur between MIBs and Mamu-B.

Conclusion

The data suggest the existence of highly divergent MIB and Mamu-B lineages on a given haplotype, as well as variable MIB and B copy numbers and configurations, at least in rhesus macaque. Recombination seems to occur between MIB and Mamu-B loci, and the resulting haplotypic plasticity at the individual level may be a strategy to better cope with pathogens. Therefore, evolutionary inferences based on the multiplicated MIB loci but also other markers close to B-genes appear to be promising for the study of B-region organization and evolution in primates.

A snapshot of the Mamu-B genes and their allelic repertoire in rhesus macaques of Chinese origin

The major histocompatibility complex class I gene repertoire was investigated in a large panel of rhesus macaques of Chinese origin. As observed in Indian animals, subjects of Chinese derivation display Mamu-B gene copy number variation, and the sum of expressed genes varies among haplotypes. In addition, these genes display differential transcription levels. The majority of the Mamu-B alleles discovered during this investigation appear to be unique for the population studied. Only one particular Mamu-B haplotype is shared between Indian and Chinese animals, and it must have been present in the progenitor stock. Hence, the data highlight the fact that most allelic polymorphism, and most of the Mamu-B haplotypes themselves, are of relatively recent origin and were most likely generated after the separation of the Indian and Chinese rhesus macaque populations.

Reshuffling of ancient peptide binding motifs between HLA-DRB multigene family members: old wine served in new skins

In most primate species, the class II region of the Major Histocompatibility Complex (MHC) displays diversity with regard to gene copy number and combination of DRB genes present per region configuration. Some of these loci exhibit extremely high levels of allelic variability, whereas others display only moderate levels of polymorphism. To understand the evolutionary history of the various HLA-DR region genes, a large number of full-length sequences of rhesus macaques, chimpanzees and humans were determined. The exon-intron organisation of the DRA gene, displaying only low levels of polymorphism, appears to have been highly conserved during primate evolution. The physical length of various DRB genes/alleles, however, fluctuates significantly in primates due to the presence of indels (insertions/deletions), mainly mapping to intron 1. Phylogenetic evidence supports the notion that the generation of new DRB genes is a dynamic and steadily ongoing process. Indeed, most of the primate DRB alleles investigated represent relatively young entities, possessing species-unique sequences. This seems to contradict the current view that the highly similar peptide binding motifs of many HLA-, Patr- and Mamu-DR molecules, encoded by exon 2 of the DRB gene, represent old entities, which predate primate speciation. As no evidence was found for convergent evolution, the combination of these two observations indicates that ancient peptide binding motifs are frequently reshuffled among duplicated members of the HLA-DRB multigene family.

Pinpointing a selective sweep to the chimpanzee MHC class I region by comparative genomics

Chimpanzees experienced a reduction of the allelic repertoire at the major histocompatibility complex (MHC) class I A and B loci, which may have been caused by a retrovirus belonging to the simian immunodeficiency virus (SIV) family. Extended MHC haplotypes were defined in a pedigreed chimpanzee colony. Comparison of genetic variation at microsatellite markers mapping inside and outside the Mhc region was carried out in humans and chimpanzees to investigate the genomic extent of the repertoire reduction. Multilocus demographic analyses underscored that chimpanzees indeed experienced a selective sweep that mainly targeted the chromosomal segment carrying the Mhc class I region. Probably due to genetic linkage, the sweep also affected other polymorphic loci, mapping in the close vicinity of the Mhc class I region genes. Nevertheless, although the allelic repertoire at particular Mhc class I and II loci appears to be limited, naturally occurring recombination events allowed the establishment of haplotype diversity after the sweep. However, recombination did not have sufficient time to erase the signal of the selective sweep.

A non-human primate model for gluten sensitivity

BACKGROUND AND AIMS

Gluten sensitivity is widespread among humans. For example, in celiac disease patients, an inflammatory response to dietary gluten leads to enteropathy, malabsorption, circulating antibodies against gluten and transglutaminase 2, and clinical symptoms such as diarrhea. There is a growing need in fundamental and translational research for animal models that exhibit aspects of human gluten sensitivity.

METHODS

Using ELISA-based antibody assays, we screened a population of captive rhesus macaques with chronic diarrhea of non-infectious origin to estimate the incidence of gluten sensitivity. A selected animal with elevated anti-gliadin antibodies and a matched control were extensively studied through alternating periods of gluten-free diet and gluten challenge. Blinded clinical and histological evaluations were conducted to seek evidence for gluten sensitivity.

RESULTS

When fed with a gluten-containing diet, gluten-sensitive macaques showed signs and symptoms of celiac disease including chronic diarrhea, malabsorptive steatorrhea, intestinal lesions and anti-gliadin antibodies. A gluten-free diet reversed these clinical, histological and serological features, while reintroduction of dietary gluten caused rapid relapse.

CONCLUSIONS

Gluten-sensitive rhesus macaques may be an attractive resource for investigating both the pathogenesis and the treatment of celiac disease.

MIC gene polymorphism and haplotype diversity in rhesus macaques

Rhesus macaques (Macaca mulatta) mainly originating from India were analysed for their major histocompatibility complex class I-related (MIC) gene repertoire. Thus far, three distinct genes, designated MIC1, MIC2 and MIC3, have been identified in the rhesus macaque. In addition, an MICD pseudogene has been described mapping apart from the other loci in a telomeric direction. Genomic comparisons and the presence of a characteristic microsatellite in exon 5 suggest that the MIC1 gene is the equivalent of the human MICA gene. Hence, the MIC2 gene, lacking the microsatellite - as do humans -, is considered to be the equivalent of human MICB. The MIC3 gene, a hybrid of MICA and MICB, seems to be generated by a crossing-over event with one breakpoint in intron 3 and accordingly is named MICA/B. Apart from their human counterparts, MICA, MICB and MICA/B cluster in separate branches in the phylogenetic tree, confirming the hybrid character of the MICA/B gene. Population analyses have shown that the various genes display polymorphism, and six MICA, five MICB and three MICA/B alleles have been identified. In the panel of homozygous typing cells, two distinct haplotype configurations have been defined by segregation analyses. Each haplotype comprises an MICB gene in conjunction with either an MICA or an MICA/B gene. Furthermore, the presence of a polymorphic microsatellite in the MICA and MICA/B alleles facilitates speedy and accurate haplotyping.

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

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

MHC class I A region diversity and polymorphism in macaque species

The HLA-A locus represents a single copy gene that displays abundant allelic polymorphism in the human population, whereas, in contrast, a nonhuman primate species such as the rhesus macaque (Macaca mulatta) possesses multiple HLA-A-like (Mamu-A) genes, which parade varying degrees of polymorphism. The number and combination of transcribed Mamu-A genes present per chromosome display diversity in a population of Indian animals. At present, it is not clearly understood whether these different A region configurations are evolutionarily stable entities. To shed light on this issue, rhesus macaques from a Chinese population and a panel of cynomolgus monkeys (Macaca fascicularis) were screened for various A region-linked variations. Comparisons demonstrated that most A region configurations are old entities predating macaque speciation, whereas most allelic variation (>95%) is of more recent origin. The latter situation contrasts the observations of the major histocompatibility complex class II genes in rhesus and cynomolgus macaques, which share a high number of identical alleles (>30%) as defined by exon 2 sequencing.

Reactivation by exon shuffling of a conserved HLA-DR3-like pseudogene segment in a New World primate species

The common marmoset (Callithrix jacchus), a New World monkey species with a limited MHC class II repertoire, is highly susceptible to certain bacterial infections. Genomic analysis of exon 2 sequences documented the existence of only one DRB region configuration harboring three loci. Two of these loci display moderate levels of allelic polymorphism, whereas the -DRB*W12 gene appears to be monomorphic. This study shows that only the Caja-DRB*W16 and -DRB*W12 loci produce functional transcripts. The Caja-DRB1*03 locus is occupied by a pseudogene, given that most of the transcripts, if detected at all, show imperfections and are present at low levels. Moreover, two hybrid transcripts were identified that feature the evolutionarily conserved peptide-binding motif characteristic for the Caja-DRB1*03 gene. Thus, the severely reduced MHC class II repertoire in common marmosets has been expanded by reactivation of a pseudogene segment as a result of exon shuffling.

Extensive sharing of MHC class II alleles between rhesus and cynomolgus macaques

In contrast to rhesus monkeys, substantial knowledge on cynomolgus monkey major histocompatibility complex (MHC) class II haplotypes is lacking. Therefore, 17 animals, including one pedigreed family, were thoroughly characterized for polymorphic Mhc class II region genes as well as their mitochondrial DNA (mtDNA) sequences. Different cynomolgus macaque populations appear to exhibit unique mtDNA profiles reflecting their geographic origin. Within the present panel, 10 Mafa-DPB1, 14 Mafa-DQA1, 12 Mafa-DQB1, and 35 Mafa-DRB exon 2 sequences were identified. All of these alleles cluster into lineages that were previously described for rhesus macaques. Moreover, about half of the Mafa-DPB1, Mafa-DQA1, and Mafa-DQB1 alleles and one third of the Mafa-DRB exon 2 sequences are identical to rhesus macaque orthologues. Such a high level of Mhc class II allele sharing has not been reported for primate species. Pedigree analysis allowed the characterization of nine distinct Mafa class II haplotypes, and seven additional ones could be deduced. Two of these haplotypes harbor a duplication of the Mafa-DQB1 locus. Despite extensive allele sharing, rhesus and cynomolgus monkeys do not appear to possess identical Mhc class II haplotypes, thus illustrating that new haplotypes were generated after speciation by recombination-like processes.

Microsatellite typing of the rhesus macaque MHC region

To improve the results gained by serotyping rhesus macaque major histocompatibility complex (MHC) antigens, molecular typing techniques have been established for class I and II genes. Like the rhesus macaque Mamu-DRB loci, the Mamu-A and -B are not only polymorphic but also polygenic. As a consequence, sequence-based typing of these genes is time-consuming. Therefore, eight MHC-linked microsatellites, or short tandem repeats (STRs), were evaluated for their use in haplotype characterization. Polymorphism analyses in rhesus macaques of Indian and Chinese origin showed high STR allelic diversity in both populations but different patterns of allele frequency distribution between the groups. Pedigree data for class I and II loci and the eight STRs allowed us to determine extended MHC haplotypes in rhesus macaque breeding groups. STR sequencing and comparisons with the complete rhesus macaque MHC genomic map allowed the exact positioning of the markers. Strong linkage disequilibria were observed between Mamu-DR and -DQ loci and adjacent STRs. Microsatellite typing provides an efficient, robust, and quick method of genotyping and deriving MHC haplotypes for rhesus macaques regardless of their geographical origin. The incorporation of MHC-linked STRs into routine genetic tests will contribute to efforts to improve the genetic characterization of the rhesus macaque for biomedical research and can provide comparative information about the evolution of the MHC region.

Reduced MIC gene repertoire variation in West African chimpanzees as compared to humans

The human major histocompatibility complex class I chain-related (MIC) genes are members of a multicopy family showing similarity to the classical HLA-A, HLA-B, and HLA-C genes. Only the MICA and MICB genes produce functional transcripts. In chimpanzees, however, only one MIC gene is expressed, showing an intermediate character, resulting from a deletion fusing the MICA and MICB gene segments together. The present population study illustrates that all chimpanzee haplotypes sampled possess the hybrid MICA/B gene. In contrast to the human situation this gene displays reduced allelic variation. The observed repertoire reduction of the chimpanzee MICA/B gene is in conformity with the severe repertoire condensation documented for Patr-B locus lineages, probably due to the close proximity of both genes.

Unparalleled complexity of the MHC class I region in rhesus macaques

The highly polymorphic gene products of the classical MHC class I genes in humans (HLA-A, HLA-B, and HLA-C) play a critical role in the immune defense against intracellular infections. Because non-human primates are important models for AIDS vaccine research, rhesus monkeys from a thoroughly pedigreed and serotyped colony were subjected to full-length cDNA analysis of MHC class I gene transcripts. Rhesus macaques express multiple dominant Mamu-A and Mamu-B transcripts (majors) per chromosome, which are characterized by high expression levels. The presence of additional cDNAs with low levels of expression (minors) suggests evidence for transcriptional control of MHC class I genes. Moreover, phylogenetic analyses illustrate that most of the Mamu-A and Mamu-B loci/lineages identified display no or only limited levels of allelic polymorphism. Thus, MHC class I diversity in rhesus macaques is typified by the existence of an unmatched high number of Mamu-A and Mamu-B region configurations that exhibit polymorphism with regard to the number and combination of transcribed loci present per chromosome.

Evolutionary stability of MHC class II haplotypes in diverse rhesus macaque populations

A thoroughly characterized breeding colony of 172 pedigreed rhesus macaques was used to analyze exon 2 of the polymorphic Mamu- DPB1, -DQA1, -DQB1, and - DRB loci. Most of the monkeys or their ancestors originated in India, though the panel also included animals from Burma and China, as well as some of unknown origin and mixed breeds. In these animals, mtDNA appears to correlate with the aforementioned geographic origin, and a large number of Mamu class II alleles were observed. The different Mamu- DPB1 alleles were largely shared between monkeys of different origin, whereas in humans particular alleles appear to be unique for ethnic populations. In contrast to Mamu-DPB1, the highly polymorphic - DQA1/DQB1 alleles form tightly linked pairs that appear to be about two-thirds population specific. For most of the DQA1/DQB1 pairs, Mamu- DRB region configurations present on the same chromosome have been ascertained, resulting in 41 different -DQ/DRB haplotypes. These distinct DQ/DRB haplotypes seem to be specific for monkeys of a determined origin. Thus, in evolutionary terms, the Mamu-DP, -DQ, and -DR regions show increasing instability with regard to allelic polymorphism, such as for -DP/DQ, or gene content and allelic polymorphism, such as for -DR, resulting in population-specific class II haplotypes. Furthermore, novel haplotypes are generated by recombination-like events. The results imply that mtDNA analysis in combination with Mhc typing is a helpful tool for selecting animals for biomedical experiments.

Evidence for an ancient selective sweep in the MHC class I gene repertoire of chimpanzees

MHC class I molecules play an essential role in the immune defense against intracellular infections. The hallmark of the MHC is its extensive degree of polymorphism at the population level. However, the present comparison of MHC class I gene intron variation revealed that chimpanzees have experienced a severe repertoire reduction at the orthologues of the HLA-A, -B, and -C loci. The loss of variability predates the (sub)speciation of chimpanzees and did not effect other known gene systems. Therefore the selective sweep in the MHC class I gene may have resulted from a widespread viral infection. Based on the present results and the fact that chimpanzees have a natural resistance to the development of AIDS, we hypothesize that the selective sweep was caused by the chimpanzee-derived simian immunodeficiency virus (SIVcpz), the closest relative of HIV-1, or a closely related retrovirus. Hence, the contemporary chimpanzee populations represent the offspring of AIDS-resistant animals, the survivors of a HIV-like pandemic that took place in the distant past.

Extensive Mhc-DQB variation in humans and non-human primate species

Non-human primates are often used in biomedical research, and the application of these animals as a model in immune-related diseases necessitates the characterisation of their MHC system. In particular, the MHC class II regions of the chimpanzee (Pan troglodytes), the rhesus macaque (Macaca mulatta)and the common marmoset (Callithrix jacchus) have been subject of molecular biological studies in recent years. In this study the emphasis was on MHC class II genes of another macaque species, Macaca fascicularis(crab eating macaque or cynomolgous monkey). The exon 2 of the Mhc-DQB gene (Mafa-DQB) was sequenced in each of a random panel of 60 non-pedigreed cynomolgous monkeys. This resulted in the detection of 23 Mafa-DQB1alleles that had not previously been published. In addition, unreported alleles were found in chimpanzees, rhesus macaques, orang-utans (Pongo pygmaeus) and stump-tailed macaques (Macaca arctoides), of which a few individuals were included in this study. Phylogenetic analyses confirm the trans-species model of evolution of the MHC-DQBlineages, in which a group of major alleles is passed on in the phylogeny, and has led to the sharing of allelic lineages by different species of non-human primates. The sharing of alleles is observed only for closely related macaque species. Furthermore, this manuscript provides an overview of all published, and whenever necessary corrected, non-human primate Mhc-DQB exon 2 alleles.