Comprehensive characterization of MHC class II haplotypes in Mauritian cynomolgus macaques

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

Rhesus macaque inhibitory and activating KIR3D interact with Mamu-A-encoded ligands

Specific interactions between killer cell Ig-like receptors (KIRs) and MHC class I ligands have not been described in rhesus macaques despite their importance in biomedical research. Using KIR-Fc fusion proteins, we detected specific interactions for three inhibitory KIRs (3DLW03, 3DL05, 3DL11) and one activating KIR (3DS05). As ligands we identified Macaca mulatta MHC (Mamu)-A1- and Mamu-A3-encoded allotypes, among them Mamu-A1*001:01, which is well known for association with slow progression to AIDS in the rhesus macaque experimental SIV infection model. Interactions with Mamu-B or Mamu-I molecules were not found. KIR3DLW03 and KIR3DL05 differ in their binding sites to their shared ligand Mamu-A1*001:01, with 3DLW03 depending on presence of the α1 domain, whereas 3DL05 depends on both the α1 and α2 domains. Fine-mapping studies revealed that binding of KIR3DLW03 is influenced by presence of the complete Bw4 epitope (positions 77, 80-83), whereas that of KIR3DL05 is mainly influenced by amino acid position 77 of Bw4 and positions 80-83 of Bw6. Our findings allowed the successful prediction of a further ligand of KIR3DL05, Mamu-A1*002:01. These functional differences of rhesus macaque KIR3DL molecules are in line with the known genetic diversification of lineage II KIRs in macaques.     

Multiple KIR gene polymorphisms are associated with plasma viral loads in SIV-infected rhesus macaques

Innate immune mechanisms play a deterministic role in the rate of disease progression during acute infection in HIV infected humans and SIV infection of non-human primates. The role NK cells play in mediating such an effect has thus gained importance. One of the major sets of molecules that regulate NK cell function are the killer cell immunoglobulin-like molecules (KIR’s). Our laboratory has previously shown an association of KIR3DL alleles 13 and 14 with high plasma viral loads in a cohort of SIV-infected rhesus macaques. To gain a more detailed understanding of the role of KIR polymorphisms, our laboratory herein conducted studies of three additional KIR loci and show that select KIR3DH alleles appear to be more strongly associated with high plasma viral loads than KIR3DL alleles 13 and 14. In addition, we herein document the existence of additional new alleles for the KIR1D, KIR2DL4, and the KIR3DH loci.     

Characterisation of mouse monoclonal antibodies against rhesus macaque killer immunoglobulin-like receptors KIR3D

Killer immunoglobulin-like receptors (KIRs) represent a highly polymorphic and diverse gene family in rhesus macaques. Analyses of the respective gene products have been hampered until now due to non-availability of specific monoclonal antibodies and failure of cross-reactivity of anti-human KIR antibodies. We utilised one activating (KIR3DSW08) and two inhibitory (KIR3DLW03 and KIR3DL05) rhesus macaque KIR-Fc fusion proteins for generation of monoclonal antibodies in mice. Besides broadly reacting ones, we obtained anti-rhesus macaque KIR antibodies with intermediate and with single specificity. These monoclonal antibodies were tested for binding to a panel of rhesus macaque KIR proteins after heterologous expression on transiently transfected cells. Epitope mapping identified two polymorphic regions that are located next to each other in the mature KIR proteins. The availability of monoclonal antibodies against rhesus macaque KIR proteins will enable future studies on KIR at the protein level in rhesus macaques as important animal models of human infectious diseases.     

KIR polymorphisms modulate peptide-dependent binding to an MHC class I ligand with a Bw6 motif

Molecular interactions between killer immunoglobulin-like receptors (KIRs) and their MHC class I ligands play a central role in the regulation of natural killer (NK) cell responses to viral pathogens and tumors. Here we identify Mamu-A1*00201 (Mamu-A*02), a common MHC class I molecule in the rhesus macaque with a canonical Bw6 motif, as a ligand for Mamu-KIR3DL05. Mamu-A1*00201 tetramers folded with certain SIV peptides, but not others, directly stained primary NK cells and Jurkat cells expressing multiple allotypes of Mamu-KIR3DL05. Differences in binding avidity were associated with polymorphisms in the D0 and D1 domains of Mamu-KIR3DL05, whereas differences in peptide-selectivity mapped to the D1 domain. The reciprocal exchange of the third predicted MHC class I-contact loop of the D1 domain switched the specificity of two Mamu-KIR3DL05 allotypes for different Mamu-A1*00201-peptide complexes. Consistent with the function of an inhibitory KIR, incubation of lymphocytes from Mamu-KIR3DL05(+) macaques with target cells expressing Mamu-A1*00201 suppressed the degranulation of tetramer-positive NK cells. These observations reveal a previously unappreciated role for D1 polymorphisms in determining the selectivity of KIRs for MHC class I-bound peptides, and identify the first functional KIR-MHC class I interaction in the rhesus macaque. The modulation of KIR-MHC class I interactions by viral peptides has important implications to pathogenesis, since it suggests that the immunodeficiency viruses, and potentially other types of viruses and tumors, may acquire changes in epitopes that increase the affinity of certain MHC class I ligands for inhibitory KIRs to prevent the activation of specific NK cell subsets.     

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. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Simian immunodeficiency virus SIVmac239 infection of major histocompatibility complex-identical cynomolgus macaques from Mauritius

Nonhuman primates are widely used to study correlates of protective immunity in AIDS research. Successful cellular immune responses have been difficult to identify because heterogeneity within macaque major histocompatibility complex (MHC) genes results in quantitative and qualitative differences in immune responses. Here we use microsatellite analysis to show that simian immunodeficiency virus (SIV)-susceptible cynomolgus macaques (Macaca fascicularis) from the Indian Ocean island of Mauritius have extremely simple MHC genetics, with six common haplotypes accounting for two-thirds of the MHC haplotypes in feral animals. Remarkably, 39% of Mauritian cynomolgus macaques carry at least one complete copy of the most frequent MHC haplotype, and 8% of these animals are homozygous. In stark contrast, entire MHC haplotypes are rarely conserved in unrelated Indian rhesus macaques. After intrarectal infection with highly pathogenic SIVmac239 virus, a pair of MHC-identical Mauritian cynomolgus macaques mounted concordant cellular immune responses comparable to those previously reported for a pair of monozygotic twins infected with the same strain of human immunodeficiency virus. Our identification of relatively abundant SIV-susceptible, MHC-identical macaques will facilitate research into protective cellular immunity.     

Characterization of 47 MHC class I sequences in Filipino cynomolgus macaques

Cynomolgus macaques (Macaca fascicularis) provide increasingly common models for infectious disease research. Several geographically distinct populations of these macaques from Southeast Asia and the Indian Ocean island of Mauritius are available for pathogenesis studies. Though host genetics may profoundly impact results of such studies, similarities and differences between populations are often overlooked. In this study we identified 47 full-length MHC class I nucleotide sequences in 16 cynomolgus macaques of Filipino origin. The majority of MHC class I sequences characterized (39 of 47) were unique to this regional population. However, we discovered eight sequences with perfect identity and six sequences with close similarity to previously defined MHC class I sequences from other macaque populations. We identified two ancestral MHC haplotypes that appear to be shared between Filipino and Mauritian cynomolgus macaques, notably a Mafa-B haplotype that has previously been shown to protect Mauritian cynomolgus macaques against challenge with a simian/human immunodeficiency virus, SHIV_89.6P. We also identified a Filipino cynomolgus macaque MHC class I sequence for which the predicted protein sequence differs from Mamu-B*17 by a single amino acid. This is important because Mamu-B*17 is strongly associated with protection against simian immunodeficiency virus (SIV) challenge in Indian rhesus macaques. These findings have implications for the evolutionary history of Filipino cynomolgus macaques as well as for the use of this model in SIV/SHIV research protocols. Kevin J. Campbell and Ann M. Detmer contributed equally to this work.     

Influence of interleukin IL-2 and IL-12 + IL-18 on surface expression of immunoglobulin-like receptors KIR2DL1, KIR2DL2, and KIR3DL2 in natural killer cells

Natural killer (NK) cells express killer cell inhibitory receptors (KIRs) that recognize polymorphic class I MHC molecules. In the present study, we analyze the modulatory effect of IL-2 alone or a combination of IL-12 with IL-18 on surface expression of killer cell immunoglobulin-like receptors KIR2DL1, KIR2DL2, and KIR3DL2 in NK cells. Thus, it was found that IL-2 causes a significant increase in the proportion of cells with given studied receptors. Stimulation by a mixture of IL-12 and IL-18 caused significant increase in the fraction of cells with the KIR2DL1 and KIR2DL2, however no significant change in the percentage of cells with KIR3DL2 receptor on their surface was observed. The results of the study show the presence of KIRs on both resting and activated NK cells, this may suggest that KIRs have also an important role in the regulatory processes after activation of this subpopulation of cells.     

Pathogenicity of simian-human immunodeficiency virus SHIV-89.6P and SIVmac is attenuated in cynomolgus macaques and associated with early T-lymphocyte responses

Because most studies of AIDS pathogenesis in nonhuman primates have been performed in Indian-origin rhesus macaques (Macaca mulatta), little is known about lentiviral pathogenicity and control of virus replication following infection of alternative macaque species. Here, we report the consequences of simian-human immunodeficiency virus SHIV-89.6P and SIVmac251 infection in cynomolgus (Macaca fascicularis) and rhesus macaques of Chinese origin. Compared to the pathogenicity of the same viruses in Indian rhesus macaques, both cynomolgus and Chinese rhesus macaques showed lower levels of plasma virus. By 9 to 10 months after infection, both viruses became undetectable in plasma more frequently in cynomolgus than in either Chinese or Indian rhesus macaques. Furthermore, after SHIV-89.6P infection, CD4+ T-cell numbers declined less and survival was longer in cynomolgus and Chinese rhesus macaques than in Indian rhesus macaques. This attenuated pathogenicity was associated with gamma interferon ELISPOT responses to Gag and Env that were generated earlier and of higher frequency in cynomolgus than in Indian rhesus macaques. Cynomolgus macaques also developed higher titer neutralizing antibodies against SHIV-89.6 at 10 and 20 weeks postinoculation than Indian rhesus macaques. These studies demonstrate that the pathogenicity of nonhuman primate lentiviruses varies markedly based on the species or geographic origin of the macaques infected and suggest that the cellular immune responses may contribute to the control of pathogenicity in cynomolgus macaques. While cynomolgus and Chinese rhesus macaques provide alternative animal models of lentiviral infection, the lower levels of viremia in cynomolgus macaques limit the usefulness of infection of this species for vaccine trials that utilize viral load as an experimental endpoint.     

Evolutionary patterns of killer cell Ig-like receptor genes in Old World monkeys

Killer cell Ig-like receptors (KIRs) modulate the cytotoxic effects of Natural Killer cells. KIR genes are encoded in the Leucocyte Receptor Complex and are characterized by their high haplotypic diversity and polymorphism. The KIR system has been studied in only three species of Old World monkeys, the rhesus macaque, the cynomolgus macaque, and the sabaeus monkey, displaying a complexity rivaling that of hominids (human and apes). Here we analyzed bacterial artificial chromosome draft sequences spanning the KIR haplotype of three other Old World monkeys, the vervet monkey (Chlorocebus aethiops), the olive baboon (Papio anubis) and the colobus monkey (Colobus guereza). A total of 25 KIR gene models were identified in these species, predicted to encode receptors with 1, 2, and 3 extracellular Ig domains, all of them with long cytoplasmic domains having two putative ITIMs, although three had a positively charged residue in the transmembrane domain. Sequence and phylogenetic analyses showed that most Old World monkeys shared five classes of KIR loci: i) KIR2DL5/3DL20 in the most centromeric region, followed by ii) the single Ig domain-encoding locus KIR1D, iii) the pseudogene KIR2DP, iv) the conserved KIR2DL4, and v) the highly diversified KIR3DL/H loci in the telomeric half of the cluster. An exception to this pattern was the KIR haplotype of the colobus monkey that lacked the KIR1D, KIR2DP, and KIR2DL4 loci of the central region of the cluster. Thus, Old World monkeys display a broad spectrum of KIR haplotype variation that has been generated upon an ancestral haplotype architecture by gene duplication, gene deletion, and non-homologous recombination.     

Diversity of the killer cell Ig-like receptors of rhesus monkeys

Because the killer cell Ig-like receptors (KIRs) have only been characterized in humans and chimpanzees, we do not have a full understanding of their evolutionary history. Therefore, cDNAs encoding the KIR molecules of five rhesus monkeys were characterized, and were found to differ from the KIR molecules identified in humans and chimpanzees. Whereas only one KIR2DL4 molecule is detected in humans and chimpanzees, two distinct KIR2DL4 homologues were identified in the monkeys. Although the two human KIR3DL molecules are limited in their polymorphism, the KIR3DL homologues in the monkeys were highly polymorphic. Up to five KIR3DL homologues were identified in each monkey that was studied, and eleven distinct KIR3DL molecules were detected in the five rhesus monkeys. Two novel families of KIR molecules were identified in the rhesus monkeys, KIR3DH and KIR1D. The KIR3DH molecules have three Ig domains, transmembrane domains homologous to KIR2DL4 molecules that contain an arginine, and short cytoplasmic domains. With these features, the KIR3DH molecules resemble the activating forms of the human KIR molecules. The KIR1D molecule encodes only one complete Ig domain before a frame-shift in the second Ig domain occurs, leading to early termination of the molecule. Multiple splice variants of KIR1D exist that encode at least one Ig domain, as well as transmembrane and cytoplasmic domains. The extensive diversity of the rhesus monkey KIR3DL homologues and the novel KIR3DH and KIR1D molecules suggests that the KIR family of molecules has evolved rapidly during the evolution of primates.     

Unusually high frequency MHC class I alleles in Mauritian origin cynomolgus macaques

Acute shortages of Indian origin Rhesus macaques significantly hinder HIV/AIDS research. Cellular immune responses are particularly difficult to study because only a subset of animals possess MHC class I (MHC I) alleles with defined peptide-binding specificities. To expand the pool of nonhuman primates suitable for studies of cellular immunity, we defined 66 MHC I alleles in Cynomolgus macaques (Macaca fascicularis) of Chinese, Vietnamese, and Mauritian origin. Most MHC I alleles were found only in animals from a single geographic origin, suggesting that Cynomolgus macaques from different origins are not interchangeable in studies of cellular immunity. Animals from Mauritius may be particularly valuable because >50% of these Cynomolgus macaques share the MHC class I allele combination Mafa-B*430101, Mafa-B*440101, and Mafa-B*460101. The increased MHC I allele sharing of Mauritian origin Cynomolgus macaques may dramatically reduce the overall number of animals needed to study cellular immune responses in nonhuman primates while simultaneously reducing the confounding effects of genetic heterogeneity in HIV/AIDS research.     

Evolution of killer cell Ig-like receptor (KIR) genes: definition of an orangutan KIR haplotype reveals expansion of lineage III KIR associated with the emergence of MHC-C

Orangutan (Pongo pygmaeus) MHC-C appears less evolved than human HLA-C: Popy-C is not fixed and its alleles encode only one (C1) of the two motifs for killer cell Ig-like receptor (KIR) ligands. To assess the structure and complexity of the orangutan KIR locus, the complete nucleotide sequence of an orangutan KIR haplotype was determined. The PopyKIR locus is flanked by LILR and FCAR and consists of seven genes and pseudogenes, two novel and five corresponding to known cDNA. Distinguishing all KIRs in this rapidly evolving KIR locus from the KIR3DX1 gene is an LTR33A/MLT1D element in intron 3. These two forms of KIR represent lineages that originated by duplication of a common ancestor. The conserved, framework regions of primate KIR loci comprise the 5' part of a lineage V KIR, the 3' part of a pseudogene, the complete 2DL4 gene, and the 3' part of a lineage II KIR. Although previously defined PopyKIR2DL4 alleles contain premature termination codons, the sequenced haplotype's PopyKIR2DL4 allele encodes a full-length protein. A model for KIR evolution is proposed. Distinguishing the orangutan KIR haplotype from the proposed common ancestor of primate KIR haplotypes is an increased number to give three lineage III KIR genes in the centromeric part of the locus, the site for most human lineage III genes encoding HLA-C specific KIR. Thus, expansion of lineage III KIR is associated with emergence of MHC-C.     

Genetic substructure in cynomolgus macaques (Macaca fascicularis) on the island of Mauritius

Background

Nonhuman primates are commonly used in biomedical research as animal models of human disease and behavior. Compared to common rodent models, nonhuman primates are genetically, physiologically, behaviorally and neurologically more similar to humans owing to more recent shared ancestry and therefore provide the advantage of greater translational validity in preclinical studies. The cynomolgus macaque (Macaca fascicularis) is one of the most commonly used nonhuman primates in academic and industry settings, yet population genetic research has revealed significant substructure throughout the species distribution that may confound studies. Cynomolgus monkeys introduced to Mauritius specifically have previously been thought to maintain the least genetic heterogeneity of all cynomolgus monkeys, although recent work, including work from our lab, suggests macaques from Mauritius too may harbor cryptic substructure.

Results

To evaluate putative substructure in Mauritian cynomolgus macaques, we designed a panel of 96 single nucleotide polymorphisms based on preliminary findings from previous work to screen 246 of cynomolgus monkeys from two primary suppliers. Results from this study support substructure in Mauritian macaques and suggest a minimum of two populations and maybe three on Mauritius, with moderate admixture.

Conclusion

These findings inform the natural history of these monkeys suggesting either a previously unrecognized physical or ecological barrier to gene flow on Mauritius and/or the breakdown of historic substructure resulting from the history of macaque introduction to the island. These findings are relevant to ongoing research using these models in part because of increased appreciation of segregating common variation with functional effects and may be used to better inform animal selection in preclinical research.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-748) contains supplementary material, which is available to authorized users.     

Distribution of natural killer cell immunoglobulin-like receptor sequences in three ethnic groups

Killer cell immunoglobulin-like receptors (KIRs) are members of a group of molecules that specifically recognize HLA class I ligands and are found on subsets of human lymphopoetic cells. The number of KIR loci can vary between individuals, resulting in a heterogeneous array of possible KIR genes. The range of observed profiles has been explained by the occurrence of two haplotype families termed A and B which can be distinguished on the basis of certain KIR sequences. Here we attempted to determine whether the frequencies of putative KIR loci and the two haplotype groups vary in three ethnically defined, healthy, and unrelated control populations, namely UK Caucasoid (n=136), Palestinian (n=105) and Thai (n=119). We molecularly typed genomic DNA for the presence of 12 putative KIR loci, KIR2DL1, KIR2DL2, KIR2DL3, KIR2DL4, KIR3DL1, KIR3DL2, KIR2DS1, KIR2DS2, KIR2DS3, KIR2DS4, KIR2DS5, and KIR3DS1, using modified PCR sequence-specific primers. The patterns of KIR locus frequencies combined with the similar linkage disequilibrium values suggest that there was a distinction in the distribution of the two broad haplotype groups between the populations studied. The A haplotype was always the most prevalent, but the ratio of A to B varied between populations. The frequency of B haplotype was highest in the Palestinians and lowest in the Thais (Pc<0.0001).     

The Structure of the Atypical Killer Cell Immunoglobulin-like Receptor,KIR2DL4

The engagement of natural killer cell immunoglobulin-like receptors (KIRs) with their target ligands, human leukocyte antigen (HLA) molecules, is a critical component of innate immunity. Structurally, KIRs typically have either two (D1-D2) or three (D0-D1-D2) extracellular immunoglobulin domains, with the D1 and D2 domain recognizing the α1 and α2 helices of HLA, respectively, whereas the D0 domain of the KIR3DLs binds a loop region flanking the α1 helix of the HLA molecule. KIR2DL4 is distinct from other KIRs (except KIR2DL5) in that it does not contain a D1 domain and instead has a D0-D2 arrangement. Functionally, KIR2DL4 is also atypical in that, unlike all other KIRs, KIR2DL4 has both activating and inhibitory signaling domains. Here, we determined the 2.8 Å crystal structure of the extracellular domains of KIR2DL4. Structurally, KIR2DL4 is reminiscent of other KIR2DL receptors, with the D0 and D2 adopting the C2-type immunoglobulin fold arranged with an acute elbow angle. However, KIR2DL4 self-associated via the D0 domain in a concentration-dependent manner and was observed as a tetramer in the crystal lattice by size exclusion chromatography, dynamic light scattering, analytical ultracentrifugation, and small angle x-ray scattering experiments. The assignment of residues in the D0 domain to forming the KIR2DL4 tetramer precludes an interaction with HLA akin to that observed for KIR3DL1. Accordingly, no interaction was observed to HLA by direct binding studies. Our data suggest that the unique functional properties of KIR2DL4 may be mediated by self-association of the receptor.     

Divergence population genetic analysis of hybridization between rhesus and cynomolgus macaques

The geographic ranges of rhesus ( Macaca mulatta ) and cynomolgus ( M. fascicularis ) macaques adjoin in Indochina where they appear to hybridize. We used published and newly generated DNA sequences from 19 loci spanning ~20 kb to test whether introgression has occurred between these macaque species. We studied introgression at the level of nuclear DNA and distinguished between incomplete lineage sorting of ancestral polymorphisms or interspecific gene flow. We implemented a divergence population genetics approach by fitting our data to an isolation model implemented in the software IMa. The model that posits no gene flow from the rhesus into the cynomolgus macaque was rejected ( P  = 1.99 × 10⁻⁸). Gene flow in this direction was estimated as 2 Nm ~1.2, while gene flow in the reverse direction was nonsignificantly different from zero ( P  = 0.16). The divergence time between species was estimated as ~1.3 million years. Balancing selection, a special case of incomplete sorting, was taken into consideration, as well as potential crossbreeding in captivity. Parameter estimates varied between analyses of subsets of data, although we still rejected isolation models. Geographic sampling of the data, where samples of cynomolgus macaques derived from Indochina were excluded, revealed a lost signature of gene flow, indicating that interspecific gene flow is restricted to mainland Indochina. Our results, in conjunction with those by others, justify future detailed analyses into the genetics of reproductive barriers and reticulate evolution in these two genome-enabled primates. Future studies of the natural hybridization between rhesus and cynomolgus macaques would expand the repertoire of systems available for speciation studies in primates.