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.

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.

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.

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.

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.

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.

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.

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.

Unique peptide-binding motif for Mamu-B*037:01: an MHC class I allele common to Indian and Chinese rhesus macaques.

Indian and Chinese rhesus macaques are often used in biomedical research. Genetic analyses of the major histocompatibility class I region have revealed that these macaques display a substantial level of polymorphism at Mamu-A and Mamu-B loci, which have been subject to duplication. Only a few Mamu class I allotypes are characterised for their peptide-binding motifs, although more information of this nature would contribute to a better interpretation of T cell-mediated immune responses. Here, we present the results of the characterisation of the functional properties of Mamu-B*037:01, an allotype commonly encountered in rhesus macaques of Indian and Chinese origin. Mamu-B*037:01 is seen to have a strong preference for acidic amino acids at the third residue, and for arginine, lysine, and tyrosine at the carboxyl terminus. This peptide-binding motif is not described in the human population.

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.

AIDS-protective HLA-B*27/B*57 and chimpanzee MHC class I molecules target analogous conserved areas of HIV-1/SIVcpz.

In the absence of treatment, most HIV-1-infected humans develop AIDS. However, a minority are long-term nonprogressors, and resistance is associated with the presence of particular HLA-B*27/B*57 molecules. In contrast, most HIV-1-infected chimpanzees do not contract AIDS. In comparison with humans, chimpanzees experienced an ancient selective sweep affecting the MHC class I repertoire. We have determined the peptide-binding properties of frequent chimpanzee MHC class I molecules, and show that, like HLA-B*27/B*57, they target similar conserved areas of HIV-1/SIV(cpz). In addition, many animals appear to possess multiple molecules targeting various conserved areas of the HIV-1/SIV(cpz) Gag protein, a quantitative aspect of the immune response that may further minimize the chance of viral escape. The functional characteristics of the contemporary chimpanzee MHC repertoire suggest that the selective sweep was caused by a lentiviral pandemic.

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.

DR haplotype diversity of the cynomolgus macaque as defined by its transcriptome.

The DR region of particular primate species may display allelic polymorphism and gene copy number variation (region configuration polymorphism). The sum of these distinct types of polymorphism is defined as complexity. To date, however, the DR region of cynomolgus macaques (Macaca fascicularis) has been poorly defined. Transcriptome analysis of a pedigreed colony, comprising animals from Indonesia and Indochina, revealed a total of 15 Mafa-DRA and 57 DRB alleles, specifying 28 different region configurations. The DRA alleles can be divided into two distinct lineages. One lineage is polymorphic, but the majority of the amino acid replacements map to the leader peptide. The second lineage is at best oligomorphic, and segregates with one specific Mafa-DRB allele. The number of Mafa-DRB genes ranges from two to five per haplotype. Due to the presence of pseudogenes, however, each haplotype encodes only one to three bona fide DRB transcripts. Depending on the region configuration in which the Mafa-DRB gene is embedded, identical alleles may display differential transcription levels. Region configurations appear to have been generated by recombination-like events. When genes or gene segments are relocated, it seems plausible that they may be placed in the context of distinct transcription control elements. As such, DRB region-related transcription level differences may add an extra layer of polymorphism to this section of the adaptive immune system.

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.

Genomic plasticity of the MHC class I A region in rhesus macaques: extensive haplotype diversity at the population level as revealed by microsatellites.

The Mamu-A genes of the rhesus macaque show different degrees of polymorphism, transcription level variation, and differential haplotype distribution. Per haplotype, usually one "major" transcribed gene is present, A1 (A7), in various combinations with "minor" genes, A2 to A6. In silico analysis of the physical map of a heterozygous animal revealed the presence of similar Mamu-A regions consisting of four duplication units, but with dissimilar positions of the A1 genes on both haplotypes, and in combination with different minor genes. Two microsatellites, D6S2854 and D6S2859, have been selected as potential tools to characterize this complex region. Subsequent analysis of a large breeding colony resulted in the description of highly discriminative patterns, displaying copy number variation in concert with microsatellite repeat length differences. Sequencing and segregation analyses revealed that these patterns are unique for each Mamu-A haplotype. In animals of Indian, Burmese, and Chinese origin, 19, 15, or 9 haplotypes, respectively, could be defined, illustrating the occurrence of differential block duplications and subsequent rearrangements by recombination. The haplotypes can be assigned to 12 unique combinations of genes (region configurations). Although most configurations harbor two transcribed A genes, one or three genes per haplotype are also present. Additionally, haplotypes lacking an A1 gene or with an A1 duplication appear to exist. The presence of different transcribed A genes/alleles in monkeys from various origins may have an impact on differential disease susceptibilities. The high-throughput microsatellite technique will be a valuable tool in animal selection for diverse biomedical research projects.

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.