Marit K. van der Wiel

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.

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.

Haplotype diversity generated by ancient recombination-like events in the MHC of Indian rhesus macaques

The Mamu-A, Mamu-B, and Mamu-DRB genes of the rhesus macaque show several levels of complexity such as allelic heterogeneity (polymorphism), copy number variation, differential segregation of genes/alleles present on a haplotype (diversity) and transcription level differences. A combination of techniques was implemented to screen a large panel of pedigreed Indian rhesus macaques (1,384 individuals representing the offspring of 137 founding animals) for haplotype diversity in an efficient and inexpensive manner. This approach allowed the definition of 140 haplotypes that display a relatively low degree of region variation as reflected by the presence of only 17 A, 18 B and 22 DRB types, respectively, exhibiting a global linkage disequilibrium comparable to that in humans. This finding contrasts with the situation observed in rhesus macaques from other geographic origins and in cynomolgus monkeys from Indonesia. In these latter populations, nearly every haplotype appears to be characterised by a unique A, B and DRB region. In the Indian population, however, a reshuffling of existing segments generated “new” haplotypes. Since the recombination frequency within the core MHC of the Indian rhesus macaques is relatively low, the various haplotypes were most probably produced by recombination events that accumulated over a long evolutionary time span. This idea is in accord with the notion that Indian rhesus macaques experienced a severe reduction in population during the Pleistocene due to a bottleneck caused by geographic changes. Thus, recombination-like processes appear to be a way to expand a diminished genetic repertoire in an isolated and relatively small founder population.

The extreme plasticity of killer cell Ig-like receptor (KIR) haplotypes differentiates rhesus macaques from humans

NK cells are essential in shaping immune responses and play an important role during pregnancy and in controlling infections. Killer cell immunoglobulin‐like receptors (KIRs) educate the NK cell and determine its state of activation. Our goal was to determine how the KIR repertoire of the rhesus macaque (Macaca mulatta) has been shaped during evolution. The presence or absence of 22 KIR gene groups was determined in 378 animals. Some unexpected observations were made in an outbred colony comprising animals of different origins. For instance, the KIR region appears to be highly plastic, and an unprecedented number of genotypes and haplotypes was observed. In contrast to humans, there is no distinction between group A and B haplotypes in the rhesus macaque, suggesting that different selective forces may be operative. Moreover, specific genes appear to be either present or absent in animals of different geographic origins. This extreme plasticity may have been propelled by co‐evolution with the rhesus macaque MHC class I region, which shows signatures of expansion. The mosaic‐like complexity of KIR genotypes as observed at the population level may represent an effective strategy for surviving epidemic infections.

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.

The repertoire of MHC class I genes in the common marmoset: evidence for functional plasticity.

In humans, the classical antigen presentation function of major histocompatibility complex (MHC) class I molecules is controlled by the human leukocyte antigen HLA-A, HLA-B and HLA-C loci. A similar observation has been made for great apes and Old World monkey species. In contrast, a New World monkey species such as the cotton-top tamarin (Saguinus oedipus) appears to employ the G locus for its classical antigen presentation function. At present, little is known about the classical MHC class I repertoire of the common marmoset (Callithrix jacchus), another New World monkey that is widely used in biomedical research. In the present population study, no evidence has been found for abundant transcription of classical I class genes. However, in each common marmoset, four to seven different G-like alleles were detected, suggesting that the ancestral locus has been subject to expansion. Segregation studies provided evidence for at least two G-like genes present per haplotype, which are transcribed by a variety of cell types. The alleles of these Caja-G genes cluster in separate lineages, suggesting that the loci diversified considerably after duplication. Phylogenetic analyses of the introns confirm that the Caja-G loci cluster in the vicinity of HLA-G, indicating that both genes shared an ancestor. In contrast to HLA-G, Caja-G shows considerable polymorphism at the peptide-binding sites. This observation, together with the lack of detectable transcripts of A and B-like genes, indicates that Caja-G genes have taken over the function of classical class I genes. These data highlight the extreme plasticity of the MHC class I gene system.

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.

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.

Human and Rhesus Macaque KIR Haplotypes Defined by Their Transcriptomes

The killer-cell Ig-like receptors (KIRs) play a central role in the immune recognition in infection, pregnancy, and transplantation through their interactions with MHC class I molecules. KIR genes display abundant copy number variation as well as high levels of polymorphism. As a result, it is challenging to characterize this structurally dynamic region. KIR haplotypes have been analyzed in different species using conventional characterization methods, such as Sanger sequencing and Roche/454 pyrosequencing. However, these methods are time-consuming and often failed to define complete haplotypes, or do not reach allele-level resolution. In addition, most analyses were performed on genomic DNA, and thus were lacking substantial information about transcription and its corresponding modifications. In this paper, we present a single-molecule real-time sequencing approach, using Pacific Biosciences Sequel platform to characterize the KIR transcriptomes in human and rhesus macaque (Macaca mulatta) families. This high-resolution approach allowed the identification of novel Mamu-KIR alleles, the extension of reported allele sequences, and the determination of human and macaque KIR haplotypes. In addition, multiple recombinant KIR genes were discovered, all located on contracted haplotypes, which were likely the result of chromosomal rearrangements. The relatively high number of contracted haplotypes discovered might be indicative of selection on small KIR repertoires and/or novel fusion gene products. This next-generation method provides an improved high-resolution characterization of the KIR cluster in humans and macaques, which eventually may aid in a better understanding and interpretation of KIR allele–associated diseases, as well as the immune response in transplantation and reproduction.

MHC class I diversity of olive baboons (Papio anubis) unravelled by next-generation sequencing

The olive baboon represents an important model system to study various aspects of human biology and health, including the origin and diversity of the major histocompatibility complex. After screening of a group of related animals for polymorphisms associated with a well-defined microsatellite marker, subsequent MHC class I typing of a selected population of 24 animals was performed on two distinct next-generation sequencing (NGS) platforms. A substantial number of 21 A and 80 B transcripts were discovered, about half of which had not been previously reported. Per animal, from one to four highly transcribed A alleles (majors) were observed, in addition to ones characterised by low transcripion levels (minors), such as members of the A*14 lineage. Furthermore, in one animal, up to 13 B alleles with differential transcription level profiles may be present. Based on segregation profiles, 16 Paan-AB haplotypes were defined. A haplotype encodes in general one or two major A and three to seven B transcripts, respectively. A further peculiarity is the presence of at least one copy of a B*02 lineage on nearly every haplotype, which indicates that B*02 represents a separate locus with probably a specialistic function. Haplotypes appear to be generated by recombination-like events, and the breakpoints map not only between the A and B regions but also within the B region itself. Therefore, the genetic makeup of the olive baboon MHC class I region appears to have been subject to a similar or even more complex expansion process than the one documented for macaque species.