Single haplotype analysis demonstrates rapid evolution of the killer immunoglobulin-like receptor (KIR) loci in primates

Complete Genomic Assembly of Mauritian Cynomolgus Macaque Killer Ig-like Receptor and Natural Killer Group 2 Haplotypes

Mauritian-origin cynomolgus macaques (MCMs) serve as a powerful nonhuman primate model in biomedical research due to their unique genetic homogeneity, which simplifies experimental designs. Despite their extensive use, a comprehensive understanding of crucial immune-regulating gene families, particularly killer Ig-like receptors (KIR) and NK group 2 (NKG2), has been hindered by the lack of detailed genomic reference assemblies. In this study, we employ advanced long-read sequencing techniques to completely assemble eight KIR and seven NKG2 genomic haplotypes, providing an extensive insight into the structural and allelic diversity of these immunoregulatory gene clusters. Leveraging these genomic resources, we prototype a strategy for genotyping KIR and NKG2 using short-read, whole-exome capture data, illustrating the potential for cost-effective multilocus genotyping at colony scale. These results mark a significant enhancement for biomedical research in MCMs and underscore the feasibility of broad-scale genetic investigations.

The KIR repertoire of a West African chimpanzee population is characterized by limited gene, allele, and haplotype variation

Introduction

The killer cell immunoglobulin-like receptors (KIR) play a pivotal role in modulating the NK cell responses, for instance, through interaction with major histocompatibility complex (MHC) class I molecules. Both gene systems map to different chromosomes but co-evolved during evolution. The human KIR gene family is characterized by abundant allelic polymorphism and copy number variation. In contrast, our knowledge of the KIR repertoire in chimpanzees is limited to 39 reported alleles, with no available population data. Only three genomic KIR region configurations have been mapped, and seventeen additional ones were deduced by genotyping.

Methods

Previously, we documented that the chimpanzee MHC class I repertoire has been skewed due to an ancient selective sweep. To understand the depth of the sweep, we set out to determine the full-length KIR transcriptome - in our MHC characterized pedigreed West African chimpanzee cohort - using SMRT sequencing (PacBio). In addition, the genomic organization of 14 KIR haplotypes was characterized by applying a Cas9-mediated enrichment approach in concert with long-read sequencing by Oxford Nanopore Technologies.

Results

In the cohort, we discovered 35 undescribed and 15 already recorded Patr-KIR alleles, and a novel hybrid KIR gene. Some KIR transcripts are subject to evolutionary conserved alternative splicing events. A detailed insight on the KIR region dynamics (location and order of genes) was obtained, however, only five new KIR region configurations were detected. The population data allowed to investigate the distribution of the MHC-C1 and C2-epitope specificity of the inhibitory lineage III KIR repertoire, and appears to be skewed towards C2.

Discussion

Although the KIR region is known to evolve fast, as observed in other primate species, our overall conclusion is that the genomic architecture and repertoire in West African chimpanzees exhibit only limited to moderate levels of variation. Hence, the ancient selective sweep that affected the chimpanzee MHC class I region may also have impacted the KIR system.

MHC Class I Ligands of Rhesus Macaque Killer Cell Ig-like Receptors

Definition of MHC class I ligands of rhesus macaque killer cell Ig-like receptors (KIRs) is fundamental to NK cell biology in this species as an animal model for infectious diseases, reproductive biology, and transplantation. To provide a more complete foundation for studying NK cell responses, rhesus macaque KIRs representing common allotypes of lineage II KIR genes were tested for interactions with MHC class I molecules representing diverse Macaca mulatta (Mamu)-A, -B, -E, -F, -I, and -AG alleles. KIR-MHC class I interactions were identified by coincubating reporter cell lines bearing chimeric KIR-CD3ζ receptors with target cells expressing individual MHC class I molecules and were corroborated by staining with KIR IgG-Fc fusion proteins. Ligands for 12 KIRs of previously unknown specificity were identified that fell into three general categories: interactions with multiple Mamu-Bw4 molecules, interactions with Mamu-A-related molecules, including allotypes of Mamu-AG and the hybrid Mamu-B*045:03 molecule, or interactions with Mamu-A1*012:01. Whereas most KIRs found to interact with Mamu-Bw4 are inhibitory, most of the KIRs that interact with Mamu-AG are activating. The KIRs that recognize Mamu-A1*012:01 belong to a phylogenetically distinct group of macaque KIRs with a 3-aa deletion in the D0 domain that is also present in human KIR3DL1/S1 and KIR3DL2. This study more than doubles the number of rhesus macaque KIRs with defined MHC class I ligands and identifies interactions with Mamu-AG, -B*045, and -A1*012. These findings support overlapping, but nonredundant, patterns of ligand recognition that reflect extensive functional diversification of these receptors.

Long-read assembly of major histocompatibility complex and killer cell immunoglobulin-like receptor genome regions in cynomolgus macaque

BACKGROUND

The major histocompatibility complex (MHC) and the killer cell immunoglobulin-like receptors (KIR) are key regulators of immune responses. The cynomolgus macaque, an Old World monkey species, can be applied as an important preclinical model for studying human diseases, including coronavirus disease 2019 (COVID-19). Several MHC-KIR combinations have been associated with either a poor or good prognosis. Therefore, macaques with a well-characterized immunogenetic profile may improve drug evaluation and speed up vaccine development. At present, a complete overview of the MHC and KIR haplotype organizations in cynomolgus macaques is lacking, and characterization by conventional techniques is hampered by the extensive expansion of the macaque MHC-B region that complicates the discrimination between genes and alleles.

METHODS

We assembled complete MHC and KIR genomic regions of cynomolgus macaque using third-generation long-read sequencing approach. We identified functional Mafa-B loci at the transcriptome level using locus-specific amplification in a cohort of 33 Vietnamese cynomolgus macaques.

RESULTS

This is the first physical mapping of complete MHC and KIR gene regions in a Vietnamese cynomolgus macaque. Furthermore, we identified four functional Mafa-B loci (B2, B3, B5, and B6) and showed that alleles of the Mafa-I*01, -B*056, -B*034, and -B*001 functional lineages, respectively, are highly frequent in the Vietnamese cynomolgus macaque population.

CONCLUSION

The insights into the MHC and KIR haplotype organizations and the level of diversity may refine the selection of animals with specific genetic markers for future medical research.

KIR3DL05 and KIR3DS02 Recognition of a Nonclassical MHC Class I Molecule in the Rhesus Macaque Implicated in Pregnancy Success

Knowledge of the MHC class I ligands of rhesus macaque killer-cell Ig-like receptors (KIRs) is fundamental to understanding the role of natural killer (NK) cells in this species as a nonhuman primate model for infectious diseases, transplantation and reproductive biology. We previously identified Mamu-AG as a ligand for KIR3DL05. Mamu-AG is a nonclassical MHC class I molecule that is expressed at the maternal-fetal interface of the placenta in rhesus macaques similar to HLA-G in humans. Although Mamu-AG and HLA-G share similar molecular features, including limited polymorphism and a short cytoplasmic tail, Mamu-AG is considerably more polymorphic. To determine which allotypes of Mamu-AG serve as ligands for KIR3DL05, we tested reporter cell lines expressing five different alleles of KIR3DL05 (KIR3DL05*001, KIR3DL05*004, KIR3DL05*005, KIR3DL05*008 and KIR3DL05*X) for responses to target cells expressing eight different alleles of Mamu-AG. All five allotypes of KIR3DL05 responded to Mamu-AG2*01:01, two exhibited dominant responses to Mamu-AG1*05:01, and three had low but detectable responses to Mamu-AG3*03:01, -AG3*03:02, -AG3*03:03 and -AG3*03:04. Since KIR3DL05*X is the product of recombination between KIR3DL05 and KIR3DS02, we also tested an allotype of KIR3DS02 (KIR3DS02*004) and found that this activating KIR also recognizes Mamu-AG2*01:01. Additional analysis of Mamu-AG variants with single amino acid substitutions identified residues in the α1-domain essential for recognition by KIR3DL05. These results reveal variation in KIR3DL05 and KIR3DS02 responses to Mamu-AG and define Mamu-AG polymorphisms that differentially affect KIR recognition.

Remarkably low KIR and HLA diversity in Amerindians reveals signatures of strong purifying selection shaping the centromeric KIR region

The killer-cell immunoglobulin-like receptors (KIR) recognize human leukocyte antigen (HLA) molecules to regulate the cytotoxic and inflammatory responses of natural killer cells. KIR genes are encoded by a rapidly evolving gene family on chromosome 19 and present an unusual variation of presence and absence of genes and high allelic diversity. Although many studies have associated KIR polymorphism with susceptibility to several diseases over the last decades, the high-resolution allele-level haplotypes have only recently started to be described in populations. Here, we use a highly innovative custom next-generation sequencing method that provides a state-of-art characterization of KIR and HLA diversity in 706 individuals from eight unique South American populations: five Amerindian populations from Brazil (three Guarani and two Kaingang); one Amerindian population from Paraguay (Aché); and two urban populations from Southern Brazil (European and Japanese descendants from Curitiba). For the first time, we describe complete high-resolution KIR haplotypes in South American populations, exploring copy number, linkage disequilibrium, and KIR-HLA interactions. We show that all Amerindians analyzed to date exhibit the lowest numbers of KIR-HLA interactions among all described worldwide populations, and that 83-97% of their KIR-HLA interactions rely on a few HLA-C molecules. Using multiple approaches, we found signatures of strong purifying selection on the KIR centromeric region, which codes for the strongest NK cell educator receptors, possibly driven by the limited HLA diversity in these populations. Our study expands the current knowledge of KIR genetic diversity in populations to understand KIR-HLA coevolution and its impact on human health and survival.

Rapid Characterization of Complex Killer Cell Immunoglobulin-Like Receptor (KIR) Regions Using Cas9 Enrichment and Nanopore Sequencing

Long-read sequencing approaches have considerably improved the quality and contiguity of genome assemblies. Such platforms bear the potential to resolve even extremely complex regions, such as multigenic immune families and repetitive stretches of DNA. Deep sequencing coverage, however, is required to overcome low nucleotide accuracy, especially in regions with high homopolymer density, copy number variation, and sequence similarity, such as the MHC and KIR gene clusters of the immune system. Therefore, we have adapted a targeted enrichment protocol in combination with long-read sequencing to efficiently annotate complex KIR gene regions. Using Cas9 endonuclease activity, segments of the KIR gene cluster were enriched and sequenced on an Oxford Nanopore Technologies platform. This provided sufficient coverage to accurately resolve and phase highly complex KIR haplotypes. Our strategy eliminates PCR-induced amplification errors, facilitates rapid characterization of large and complex multigenic regions, including its epigenetic footprint, and is applicable in multiple species, even in the absence of a reference genome.

Rhesus Macaque Activating Killer Immunoglobulin-Like Receptors Associate With Fc Receptor Gamma (FCER1G) and Not With DAP12 Adaptor Proteins Resulting in Stabilized Expression and Enabling Signal Transduction

Activating killer cell immunoglobulin-like receptors (KIR) in macaques are thought to be derived by genetic recombination of the region encoding the transmembrane and intracellular part of KIR2DL4 and a KIR3D gene. As a result, all macaque activating KIR possess a positively charged arginine residue in the transmembrane region. As human KIR2DL4 associates with the FCER1G (also called Fc receptor-gamma, FcRγ) adaptor, we hypothesized that in contrast to human and great ape the activating KIRs of macaques associate with FcRγ instead of DAP12. By applying co-immunoprecipitation of transfected as well as primary cells, we demonstrate that rhesus macaque KIR3DS05 indeed associates with FcRγ and not with DAP12. This association with FcRγ results in increased and substantially stabilized surface expression of KIR3DS05. In addition, we demonstrate that binding of specific ligands of KIR3DS05, Mamu-A1*001 and A1*011, resulted in signal transduction in the presence of FcRγ in contrast to DAP12.

Sequence diversity analyses of an improved rhesus macaque genome enhance its biomedical utility

The rhesus macaque (Macaca mulatta) is the most widely studied nonhuman primate (NHP) in biomedical research. We present an updated reference genome assembly (Mmul_10, contig N50 = 46 Mbp) that increases the sequence contiguity 120-fold and annotate it using 6.5 million full-length transcripts, thus improving our understanding of gene content, isoform diversity, and repeat organization. With the improved assembly of segmental duplications, we discovered new lineage-specific genes and expanded gene families that are potentially informative in studies of evolution and disease susceptibility. Whole-genome sequencing (WGS) data from 853 rhesus macaques identified 85.7 million single-nucleotide variants (SNVs) and 10.5 million indel variants, including potentially damaging variants in genes associated with human autism and developmental delay, providing a framework for developing noninvasive NHP models of human disease.

Efficient Sequencing, Assembly, and Annotation of Human KIR Haplotypes

The homology, recombination, variation, and repetitive elements in the natural killer-cell immunoglobulin-like receptor (KIR) region has made full haplotype DNA interpretation impossible in a high-throughput workflow. Here, we present a new approach using long-read sequencing to efficiently capture, sequence, and assemble diploid human KIR haplotypes. Probes were designed to capture KIR fragments efficiently by leveraging the repeating homology of the region. IDT xGen® Lockdown probes were used to capture 2-8 kb of sheared DNA fragments followed by sequencing on a PacBio Sequel. The sequences were error corrected, binned, and then assembled using the Canu assembler. The location of genes and their exon/intron boundaries are included in the workflow. The assembly and annotation was evaluated on 16 individuals (8 African American and 8 Europeans) from whom ground truth was known via long-range sequencing with fosmid library preparation. Using only 18 capture probes, the results show that the assemblies cover 97% of the GenBank reference, are 99.97% concordant, and it takes only 1.8 haplotigs to cover 75% of the reference. We also report the first assembly of diploid KIR haplotypes from long-read WGS. Our targeted hybridization probe capture and sequencing approach is the first of its kind to fully sequence and phase all diploid human KIR haplotypes, and it is efficient enough for population-scale studies and clinical use. The open and free software is available at https://github.com/droeatumn/kass and supported by a environment at https://hub.docker.com/repository/docker/droeatumn/kass.

The Genetic Mechanisms Driving Diversification of the KIR Gene Cluster in Primates

The activity and function of natural killer (NK) cells are modulated through the interactions of multiple receptor families, of which some recognize MHC class I molecules. The high level of MHC class I polymorphism requires their ligands either to interact with conserved epitopes, as is utilized by the NKG2A receptor family, or to co-evolve with the MHC class I allelic variation, which task is taken up by the killer cell immunoglobulin-like receptor (KIR) family. Multiple molecular mechanisms are responsible for the diversification of the KIR gene system, and include abundant chromosomal recombination, high mutation rates, alternative splicing, and variegated expression. The combination of these genetic mechanisms generates a compound array of diversity as is reflected by the contraction and expansion of KIR haplotypes, frequent birth of fusion genes, allelic polymorphism, structurally distinct isoforms, and variegated expression, which is in contrast to the mainly allelic nature of MHC class I polymorphism in humans. A comparison of the thoroughly studied human and macaque KIR gene repertoires demonstrates a similar evolutionarily conserved toolbox, through which selective forces drove and maintained the diversified nature of the KIR gene cluster. This hypothesis is further supported by the comparative genetics of KIR haplotypes and genes in other primate species. The complex nature of the KIR gene system has an impact upon the education, activity, and function of NK cells in coherence with an individual’s MHC class I repertoire and pathogenic encounters. Although selection operates on an individual, the continuous diversification of the KIR gene system in primates might protect populations against evolving pathogens.

Unparalleled Rapid Evolution of KIR Genes in Rhesus and Cynomolgus Macaque Populations

The killer cell Ig-like receptors (KIR) modulate immune responses through interactions with MHC class I molecules. The KIR region in large cohorts of rhesus and cynomolgus macaque populations were characterized, and the experimental design enabled the definition of a considerable number of alleles (n = 576) and haplotypes, which are highly variable with regard to architecture. Although high levels of polymorphism were recorded, only a few alleles are shared between species and populations. The rapid evolution of allelic polymorphism, accumulated by point mutations, was further confirmed by the emergence of a novel KIR allele in a rhesus macaque family. In addition to allelic variation, abundant orthologous and species-specific KIR genes were identified, the latter of which are frequently generated by fusion events. The concerted action of both genetic mechanisms, in combination with differential selective pressures at the population level, resulted in the unparalleled rapid evolution of the KIR gene region in two closely related macaque species. The variation of the KIR gene repertoire at the species and population level might have an impact on the outcome of preclinical studies with macaque models.

Nomenclature report for killer-cell immunoglobulin-like receptors (KIR) in macaque species: new genes/alleles, renaming recombinant entities and IPD-NHKIR updates

The Killer-cell Immunoglobulin-like Receptors (KIR) are encoded by a diverse group of genes, which are characterized by allelic polymorphism, gene duplications, and recombinations, which may generate recombinant entities. The number of reported macaque KIR sequences is steadily increasing, and these data illustrate a gene system that may match or exceed the complexity of the human KIR cluster. This report lists the names of quality controlled and annotated KIR genes/alleles with all the relevant references for two different macaque species: rhesus and cynomolgus macaques. Numerous recombinant KIR genes in these species necessitate a revision of some of the earlier-published nomenclature guidelines. In addition, this report summarizes the latest information on the Immuno Polymorphism Database (IPD)-NHKIR Database, which contains annotated KIR sequences from four non-human primate species.

Monkeying Around: Using Non-human Primate Models to Study NK Cell Biology in HIV Infections

Natural killer (NK) cells are the major innate effectors primed to eliminate virus-infected and tumor or neoplastic cells. Recent studies also suggest nuances in phenotypic and functional characteristics among NK cell subsets may further permit execution of regulatory and adaptive roles. Animal models, particularly non-human primate (NHP) models, are critical for characterizing NK cell biology in disease and under homeostatic conditions. In HIV infection, NK cells mediate multiple antiviral functions via upregulation of activating receptors, inflammatory cytokine secretion, and antibody dependent cell cytotoxicity through antibody Fc-FcR interaction and others. However, HIV infection can also reciprocally modulate NK cells directly or indirectly, leading to impaired/ineffective NK cell responses. In this review, we will describe multiple aspects of NK cell biology in HIV/SIV infections and their association with viral control and disease progression, and how NHP models were critical in detailing each finding. Further, we will discuss the effect of NK cell depletion in SIV-infected NHP and the characteristics of newly described memory NK cells in NHP models and different mouse strains. Overall, we propose that the role of NK cells in controlling viral infections remains incompletely understood and that NHP models are indispensable in order to efficiently address these deficits.

Two to Tango: Co-evolution of Hominid Natural Killer Cell Receptors and MHC

Natural killer (NK) cells have diverse roles in hominid immunity and reproduction. Modulating these functions are the interactions between major histocompatibility complex (MHC) class I molecules that are ligands for two NK cell surface receptor types. Diverse killer cell immunoglobulin-like receptors (KIR) bind specific motifs encoded within the polymorphic MHC class I cell surface glycoproteins, while, in more conserved interactions, CD94:NKG2A receptors recognize MHC-E with bound peptides derived from MHC class I leader sequences. The hominid lineage presents a choreographed co-evolution of KIR with their MHC class I ligands. MHC-A, -B, and -C are present in all great apes with species-specific haplotypic variation in gene content. The Bw4 epitope recognized by lineage II KIR is restricted to MHC-B but also present on some gorilla and human MHC-A. Common to great apes, but rare in humans, are MHC-B possessing a C1 epitope recognized by lineage III KIR. MHC-C arose from duplication of MHC-B and is fixed in all great apes except orangutan, where it exists on approximately 50% of haplotypes and all allotypes are C1-bearing. Recent study showed that gorillas possess yet another intermediate MHC organization compared to humans. Like orangutans, but unlike the Pan-Homo species, duplication of MHC-B occurred. However, MHC-C is fixed, and the MHC-C C2 epitope (absent in orangutans) emerges. The evolution of MHC-C drove expansion of its cognate lineage III KIR. Recently, position −21 of the MHC-B leader sequence has been shown to be critical in determining NK cell educational outcome. In humans, methionine (−21M) results in CD94:NKG2A-focused education whereas threonine (−21T) produces KIR-focused education. This is another dynamic position among hominids. Orangutans have exclusively −21M, consistent with their intermediate stage in lineage III KIR-focused evolution. Gorillas have both −21M and −21T, like humans, but they are unequally encoded by their duplicated B genes. Chimpanzees have near-fixed −21T, indicative of KIR-focused NK education. Harmonious with this observation, chimpanzee KIR exhibit strong binding and, compared to humans, smaller differences between binding levels of activating and inhibitory KIR. Consistent between these MHC-NK cell receptor systems over the course of hominid evolution is the evolution of polymorphism favoring the more novel and dynamic KIR system.

Extensive Alternative Splicing of KIR Transcripts

The killer-cell Ig-like receptors (KIR) form a multigene entity involved in modulating immune responses through interactions with MHC class I molecules. The complexity of the KIR cluster is reflected by, for instance, abundant levels of allelic polymorphism, gene copy number variation, and stochastic expression profiles. The current transcriptome study involving human and macaque families demonstrates that KIR family members are also subjected to differential levels of alternative splicing, and this seems to be gene dependent. Alternative splicing may result in the partial or complete skipping of exons, or the partial inclusion of introns, as documented at the transcription level. This post-transcriptional process can generate multiple isoforms from a single KIR gene, which diversifies the characteristics of the encoded proteins. For example, alternative splicing could modify ligand interactions, cellular localization, signaling properties, and the number of extracellular domains of the receptor. In humans, we observed abundant splicing for KIR2DL4, and to a lesser extent in the lineage III KIR genes. All experimentally documented splice events are substantiated by in silico splicing strength predictions. To a similar extent, alternative splicing is observed in rhesus macaques, a species that shares a close evolutionary relationship with humans. Splicing profiles of Mamu-KIR1D and Mamu-KIR2DL04 displayed a great diversity, whereas Mamu-KIR3DL20 (lineage V) is consistently spliced to generate a homolog of human KIR2DL5 (lineage I). The latter case represents an example of convergent evolution. Although just a single KIR splice event is shared between humans and macaques, the splicing mechanisms are similar, and the predicted consequences are comparable. In conclusion, alternative splicing adds an additional layer of complexity to the KIR gene system in primates, and results in a wide structural and functional variety of KIR receptors and its isoforms, which may play a role in health and disease.

Diversification of Bw4 Specificity and Recognition of a Nonclassical MHC Class I Molecule Implicated in Maternal-Fetal Tolerance by Killer Cell Ig-like Receptors of the Rhesus Macaque

The rhesus macaque is an important animal model for AIDS and other infectious diseases; however, studies to address NK cell function in this species have been limited by the lack of defined ligands for killer cell Ig-like receptors (KIRs). To identify ligands for rhesus macaque KIRs, we adopted a novel approach based on a pair of stable cell lines. NFAT-responsive luciferase reporter cell lines expressing the extracellular domains of macaque KIRs fused to the transmembrane and cytoplasmic domains of CD28 and CD3ζ were incubated with target cells expressing individual MHC class I molecules, and ligand recognition was detected by the MHC class I-dependent upregulation of luciferase. Using this approach, we found that Mamu-KIR3DL01, -KIR3DL06, -KIR3DL08, and -KIR3DSw08 all recognize Mamu-Bw4 molecules but with differing allotype specificity. In contrast, Mamu-KIR3DL05 recognizes Mamu-A and Mamu-A-related molecules, including Mamu-A1*002 and -A3*13, Mamu-B*036, the product of a recombinant Mamu-B allele with α1 and α2 domain sequences derived from a MHC-A gene, and Mamu-AG*01, a nonclassical molecule expressed on placental trophoblasts that originated from an ancestral duplication of a MHC-A gene. These results reveal an expansion of the lineage II KIRs in macaques that recognize Bw4 ligands and identify a nonclassical molecule implicated in placental development and pregnancy as a ligand for Mamu-KIR3DL05. In addition to offering new insights into KIR-MHC class I coevolution, these findings provide an important foundation for investigating the role of NK cells in the rhesus macaque as an animal model for infectious diseases and reproductive biology.

Nomenclature for the KIR of non-human species

The increasing number of Killer Immunoglobulin-like Receptor (KIR) sequences available for non-human primate species and cattle has prompted development of a centralized database, guidelines for a standardized nomenclature, and minimum requirements for database submission. The guidelines and nomenclature are based on those used for human KIR and incorporate modifications made for inclusion of non-human species in the companion IPD-NHKIR database. Included in this first release are the rhesus macaque (Macaca mulatta), chimpanzee (Pan troglodytes), orangutan (Pongo abelii and Pongo pygmaeus), and cattle (Bos taurus).

Genetics of Natural Killer Cells in Human Health, Disease, and Survival

Natural killer (NK) cells have vital functions in human immunity and reproduction. In the innate and adaptive immune responses to infection, particularly by viruses, NK cells respond by secreting inflammatory cytokines and killing infected cells. In reproduction, NK cells are critical for genesis of the placenta, the organ that controls the supply of oxygen and nutrients to the growing fetus. Controlling NK cell functions are interactions of HLA class I with inhibitory NK cell receptors. First evolved was the conserved interaction of HLA-E with CD94:NKG2A; later established were diverse interactions of HLA-A, -B, and -C with killer cell immunoglobulin-like receptors. Characterizing the latter interactions is rapid evolution, which distinguishes human populations and all species of higher primate. Driving this evolution are the different and competing selections imposed by pathogens on NK cell-mediated immunity and by the constraints of human reproduction on NK cell-mediated placentation. Promoting rapid evolution is independent segregation of polymorphic receptors and ligands throughout human populations.

Structure and Diversity of the Rhesus Macaque Immunoglobulin Loci through Multiple De Novo Genome Assemblies

The rhesus macaque is a critically important animal model in biomedical research, most recently playing a key role in the development of vaccines against human immunodeficiency virus-1. Nevertheless, the immunoglobulin (Ig) loci of macaques are as yet incompletely determined and our understanding of differences between human and macaque humoral immunity remains deficient. We completed a high-coverage, high-quality whole genome sequencing and assembly project with a single rhesus macaque of Indian origin, and partial genome assemblies using genomic molecular targeting of the Ig loci in nine other rhesus macaques of Indian origin. These data indicate that the macaque Ig loci are substantially more diverse than those in humans, including greater sequence diversity and copy-number variation between individuals. It appears likely that such copy-number variation even occurs between allelic loci within individuals. Different Ig gene families in the macaque show distinct relationships to the corresponding human gene families and appear to evolve under different mechanisms. These results raise intriguing questions about the evolution of antigen receptors in primates but also have important practical implications for the design and interpretation of biomedical studies.

AIDS in chimpanzees: the role of MHC genes

The ancestral progenitor of common chimpanzees and bonobos experienced a selective sweep that ravaged its major histocompatibility complex (MHC) class I repertoire. The causative agent was probably an ancestral retrovirus, highly related to the contemporary HIV-1 strain, which initiated the acquired immunodeficiency syndrome pandemic in the human population. As a direct result, MHC class I allotypes with the capability of targeting conserved retroviral elements were enriched in the ancestral progenitor. Even today, the impact can be traced back by studying the functional capacities of the contemporary MHC class I allotypes of common chimpanzees. Viruses, however, have developed several strategies to manipulate the cell-surface expression of MHC class I genes. Monitoring the presence and absence of the MHC class I allotypes on the cell surface is conducted, for instance, by the hosts' gene products of the killer cell immunoglobulin-like receptor (KIR) complex. Hence, one may wonder whether-in the future-any clues with regard to the signature of the MHC class I selective sweep might be unearthed for the KIR genes as well.

KIR3DL01 upregulation on gut natural killer cells in response to SIV infection of KIR- and MHC class I-defined rhesus macaques

Natural killer cells provide an important early defense against viral pathogens and are regulated in part by interactions between highly polymorphic killer-cell immunoglobulin-like receptors (KIRs) on NK cells and their MHC class I ligands on target cells. We previously identified MHC class I ligands for two rhesus macaque KIRs: KIR3DL01 recognizes Mamu-Bw4 molecules and KIR3DL05 recognizes Mamu-A1*002. To determine how these interactions influence NK cell responses, we infected KIR3DL01+ and KIR3DL05+ macaques with and without defined ligands for these receptors with SIVmac239, and monitored NK cell responses in peripheral blood and lymphoid tissues. NK cell responses in blood were broadly stimulated, as indicated by rapid increases in the CD16+ population during acute infection and sustained increases in the CD16+ and CD16-CD56- populations during chronic infection. Markers of proliferation (Ki-67), activation (CD69 & HLA-DR) and antiviral activity (CD107a & TNFα) were also widely expressed, but began to diverge during chronic infection, as reflected by sustained CD107a and TNFα upregulation by KIR3DL01+, but not by KIR3DL05+ NK cells. Significant increases in the frequency of KIR3DL01+ (but not KIR3DL05+) NK cells were also observed in tissues, particularly in the gut-associated lymphoid tissues, where this receptor was preferentially upregulated on CD56+ and CD16-CD56- subsets. These results reveal broad NK cell activation and dynamic changes in the phenotypic properties of NK cells in response to SIV infection, including the enrichment of KIR3DL01+ NK cells in tissues that support high levels of virus replication.

Improved full-length killer cell immunoglobulin-like receptor transcript discovery in Mauritian cynomolgus macaques

Killer cell immunoglobulin-like receptors (KIRs) modulate disease progression of pathogens including HIV, malaria, and hepatitis C. Cynomolgus and rhesus macaques are widely used as nonhuman primate models to study human pathogens, and so, considerable effort has been put into characterizing their KIR genetics. However, previous studies have relied on cDNA cloning and Sanger sequencing that lack the throughput of current sequencing platforms. In this study, we present a high throughput, full-length allele discovery method utilizing Pacific Biosciences circular consensus sequencing (CCS). We also describe a new approach to Macaque Exome Sequencing (MES) and the development of the Rhexome1.0, an adapted target capture reagent that includes macaque-specific capture probe sets. By using sequence reads generated by whole genome sequencing (WGS) and MES to inform primer design, we were able to increase the sensitivity of KIR allele discovery. We demonstrate this increased sensitivity by defining nine novel alleles within a cohort of Mauritian cynomolgus macaques (MCM), a geographically isolated population with restricted KIR genetics that was thought to be completely characterized. Finally, we describe an approach to genotyping KIRs directly from sequence reads generated using WGS/MES reads. The findings presented here expand our understanding of KIR genetics in MCM by associating new genes with all eight KIR haplotypes and demonstrating the existence of at least one KIR3DS gene associated with every haplotype.

Two Orangutan Species Have Evolved Different KIR Alleles and Haplotypes

The immune and reproductive functions of human NK cells are regulated by interactions of the C1 and C2 epitopes of HLA-C with C1-specific and C2-specific lineage III killer cell Ig-like receptors (KIR). This rapidly evolving and diverse system of ligands and receptors is restricted to humans and great apes. In this context, the orangutan has particular relevance because it represents an evolutionary intermediate, one having the C1 epitope and corresponding KIR but lacking the C2 epitope. Through a combination of direct sequencing, KIR genotyping, and data mining from the Great Ape Genome Project, we characterized the KIR alleles and haplotypes for panels of 10 Bornean orangutans and 19 Sumatran orangutans. The orangutan KIR haplotypes have between 5 and 10 KIR genes. The seven orangutan lineage III KIR genes all locate to the centromeric region of the KIR locus, whereas their human counterparts also populate the telomeric region. One lineage III KIR gene is Bornean specific, one is Sumatran specific, and five are shared. Of 12 KIR gene-content haplotypes, 5 are Bornean specific, 5 are Sumatran specific, and 2 are shared. The haplotypes have different combinations of genes encoding activating and inhibitory C1 receptors that can be of higher or lower affinity. All haplotypes encode an inhibitory C1 receptor, but only some haplotypes encode an activating C1 receptor. Of 130 KIR alleles, 55 are Bornean specific, 65 are Sumatran specific, and 10 are shared.

Resurrecting KIR2DP1: A Key Intermediate in the Evolution of Human Inhibitory NK Cell Receptors That Recognize HLA-C

KIR2DP1 is an inactive member of the human lineage III KIR family, which includes all HLA-C-specific receptor genes. The lethal, and only, defect in KIR2DP1 is a nucleotide deletion in codon 88. Fixed in modern humans, the deletion is also in archaic human genomes. KIR2DP1 is polymorphic, with dimorphism at specificity-determining position 44. By repairing the deletion, we resurrected 11 alleles of KIR2DP1^F , the functional antecedent of KIR2DP1 We demonstrate how K44-KIR2DP1^F with lysine 44 recognized C1⁺HLA-C, whereas T44-KIR2DP1^F recognized C2⁺HLA-C. Dimorphisms at 12 other KIR2DP1^F residues modulate receptor avidity or signaling. KIR2DP1 and KIR2DL1 are neighbors in the centromeric KIR region and are in tight linkage disequilibrium. Like KIR2DL1, KIR2DP1 contributed to CenA and CenB KIR haplotype differences. Encoded on CenA, C1-specific K44-KIR2DP1^F were stronger receptors than the attenuated C2-specific T44-KIR2DP1^F encoded on CenB The last common ancestor of humans and chimpanzees had diverse lineage III KIR that passed on to chimpanzees but not to humans. Early humans inherited activating KIR2DS4 and an inhibitory lineage III KIR, likely encoding a C1-specific receptor. The latter spawned the modern family of HLA-C receptors. KIR2DP1^F has properties consistent with KIR2DP1^F having been the founder gene. The first KIR2DP1^F alleles encoded K44-C1 receptors; subsequently KIR2DP1^F alleles encoding T44-C2 receptors evolved. The emergence of dedicated KIR2DL2/3 and KIR2DL1 genes encoding C1 and C2 receptors, respectively, could have led to obsolescence of KIR2DP1^F Alternatively, pathogen subversion caused its demise. Preservation of KIR2DP1^F functional polymorphism was a side effect of fixation of the deletion in KIR2DP1^F by micro gene conversion.

Defining KIR and HLA Class I Genotypes at Highest Resolution via High-Throughput Sequencing

The physiological functions of natural killer (NK) cells in human immunity and reproduction depend upon diverse interactions between killer cell immunoglobulin-like receptors (KIRs) and their HLA class I ligands: HLA-A, HLA-B, and HLA-C. The genomic regions containing the KIR and HLA class I genes are unlinked, structurally complex, and highly polymorphic. They are also strongly associated with a wide spectrum of diseases, including infections, autoimmune disorders, cancers, and pregnancy disorders, as well as the efficacy of transplantation and other immunotherapies. To facilitate study of these extraordinary genes, we developed a method that captures, sequences, and analyzes the 13 KIR genes and HLA-A, HLA-B, and HLA-C from genomic DNA. We also devised a bioinformatics pipeline that attributes sequencing reads to specific KIR genes, determines copy number by read depth, and calls high-resolution genotypes for each KIR gene. We validated this method by using DNA from well-characterized cell lines, comparing it to established methods of HLA and KIR genotyping, and determining KIR genotypes from 1000 Genomes sequence data. This identified 116 previously uncharacterized KIR alleles, which were all demonstrated to be authentic by sequencing from source DNA via standard methods. Analysis of just two KIR genes showed that 22% of the 1000 Genomes individuals have a previously uncharacterized allele or a structural variant. The method we describe is suited to the large-scale analyses that are needed for characterizing human populations and defining the precise HLA and KIR factors associated with disease. The methods are applicable to other highly polymorphic genes.

Co-evolution of MHC class I and variable NK cell receptors in placental mammals

Shaping natural killer (NK) cell functions in human immunity and reproduction are diverse killer cell immunoglobulin-like receptors (KIRs) that recognize polymorphic MHC class I determinants. A survey of placental mammals suggests that KIRs serve as variable NK cell receptors only in certain primates and artiodactyls. Divergence of the functional and variable KIRs in primates and artiodactyls predates placental reproduction. Among artiodactyls, cattle but not pigs have diverse KIRs. Catarrhine (humans, apes, and Old World monkeys) and platyrrhine (New World monkeys) primates, but not prosimians, have diverse KIRs. Platyrrhine and catarrhine systems of KIR and MHC class I are highly diverged, but within the catarrhines, a stepwise co-evolution of MHC class I and KIR is discerned. In Old World monkeys, diversification focuses on MHC-A and MHC-B and their cognate lineage II KIR. With evolution of C1-bearing MHC-C from MHC-B, as informed by orangutan, the focus changes to MHC-C and its cognate lineage III KIR. Evolution of C2 from C1 and fixation of MHC-C drove further elaboration of MHC-C-specific KIR, as exemplified by chimpanzee. In humans, the evolutionary trajectory changes again. Emerging from reorganization of the KIR locus and selective attenuation of KIR avidity for MHC class I are the functionally distinctive KIR A and KIR B haplotypes.

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 and KIR Polymorphisms in Rhesus Macaque SIV Infection

Natural killer lymphocytes are essentially involved as the first line of defense against agents such as viruses and malignant cells. The activity of these cells is regulated via interaction of specific and diverse killer cell immunoglobulin-like receptors (KIR) with the highly polymorphic cognate MHC class I proteins on target cells. Genetic variability of both KIR and MHC-I ligands has been shown to be associated with resistance to many diseases, including infection with the immunodeficiency virus. Disease course and progression to AIDS after infection with human immunodeficiency virus-1 (HIV-1) is essentially influenced by the presence of the stimulatory KIR3DS1 receptor in combination with HLA-Bw4. Knowledge of such genetic interactions that contribute to not only disease resistance but also susceptibility are just as important. Such combined genetic factors were recently reported in the rhesus macaque AIDS model. Here, we review the rhesus macaque MHC class I and KIR gene systems and the role of their polymorphisms in the SIV infection model.

Suppression of a Natural Killer Cell Response by Simian Immunodeficiency Virus Peptides

Natural killer (NK) cell responses in primates are regulated in part through interactions between two highly polymorphic molecules, the killer-cell immunoglobulin-like receptors (KIRs) on NK cells and their major histocompatibility complex (MHC) class I ligands on target cells. We previously reported that the binding of a common MHC class I molecule in the rhesus macaque, Mamu-A1*002, to the inhibitory receptor Mamu-KIR3DL05 is stabilized by certain simian immunodeficiency virus (SIV) peptides, but not by others. Here we investigated the functional implications of these interactions by testing SIV peptides bound by Mamu-A1*002 for the ability to modulate Mamu-KIR3DL05⁺ NK cell responses. Twenty-eight of 75 SIV peptides bound by Mamu-A1*002 suppressed the cytolytic activity of primary Mamu-KIR3DL05⁺ NK cells, including three immunodominant CD8⁺ T cell epitopes previously shown to stabilize Mamu-A1*002 tetramer binding to Mamu-KIR3DL05. Substitutions at C-terminal positions changed inhibitory peptides into disinhibitory peptides, and vice versa, without altering binding to Mamu-A1*002. The functional effects of these peptide variants on NK cell responses also corresponded to their effects on Mamu-A1*002 tetramer binding to Mamu-KIR3DL05. In assays with mixtures of inhibitory and disinhibitory peptides, low concentrations of inhibitory peptides dominated to suppress NK cell responses. Consistent with the inhibition of Mamu-KIR3DL05⁺ NK cells by viral epitopes presented by Mamu-A1*002, SIV replication was significantly higher in Mamu-A1*002⁺ CD4⁺ lymphocytes co-cultured with Mamu-KIR3DL05⁺ NK cells than with Mamu-KIR3DL05^- NK cells. These results demonstrate that viral peptides can differentially affect NK cell responses by modulating MHC class I interactions with inhibitory KIRs, and provide a mechanism by which immunodeficiency viruses may evade NK cell responses.

Natural killer (NK) cells recognize and kill infected cells without prior antigenic stimulation, and thus provide an important early defense against virus infection. NK cell responses in primates are regulated in part through interactions between two highly polymorphic molecules, the killer-cell immunoglobulin-like receptors (KIRs) on NK cells and their major histocompatibility complex (MHC) class I ligands on target cells. Inhibitory KIRs normally suppress NK cell responses through interactions with their MHC class I ligands on the surface of healthy cells. However, when these interactions are perturbed, this inhibition is lost resulting in NK cell activation and killing of the target cell. We investigated the functional implications of simian immunodeficiency virus (SIV) peptides bound by a common MHC class I molecule in the rhesus macaque that stabilize or disrupt binding to an inhibitory KIR. Whereas SIV peptides that stabilized KIR-MHC class I binding suppressed NK cell activation, peptides that disrupted this interaction did not and resulted in NK cell lysis. These findings demonstrate that viral peptides can modulate NK cell responses through KIR-MHC class I interactions, and are consistent with the possibility that human and simian immunodeficiency viruses may acquire changes in epitopes that increase the binding of MHC class I ligands to inhibitory KIRs as a mechanism to suppress NK cell responses.

The killer-cell immunoglobulin-like receptors of macaques

Natural killer (NK) cells play a central role in immune responses through direct cytotoxicity and the release of cytokines that prime adaptive immunity. In simian primates, NK cell responses are regulated by interactions between two highly polymorphic sets of molecules: the killer-cell immunoglobulin-like receptors (KIRs) and their major histocompatibility complex (MHC) class I ligands. KIR-MHC class I interactions in humans have been implicated in the outcome of a number viral diseases and cancers. However, studies to address the role of KIRs in animal models have been limited by the complex immunogenetics and lack of defined ligands for KIRs in non-human primates. Due to the rapid evolution of KIRs, there is little conservation among the KIR genes of different primate species and it is not possible to predict the specificity of KIRs from known KIR-MHC class I interactions in humans. Hence, the MHC class I ligands for KIRs in species other than humans are poorly defined. Here, we review the KIR genes of the rhesus macaque, an important animal model for human immunodeficiency virus infection and other infectious diseases, and the MHC class I ligands that have been identified for KIRs in this species.

The production of KIR-Fc fusion proteins and their use in a multiplex HLA class I binding assay

Soluble recombinant proteins that comprise the extracellular part of a surface expressed receptor attached to the Fc region of an IgG antibody have facilitated the determination of ligand specificity for an array of immune system receptors. Among such receptors is the family of killer cell immunoglobulin-like receptors (KIR) that recognize HLA class I ligands. These receptors, expressed on natural killer (NK) cells and T cells, play important roles in both immune defense and placental development in early pregnancy. Here we describe a method for the production of two domain KIR-Fc fusion proteins using baculovirus infected insect cells. This method is more scalable than traditional mammalian cell expression systems and produces efficiently folded proteins that carry posttranslational modifications found in native KIR. We also describe a multiplex binding assay using the Luminex platform that determines the avidity and specificity of two domain KIR-Fc for a panel of microbeads, each coated with one of 97 HLA class I allotypes. This assay is simple to perform, and represents a major improvement over the assays used previously, which were limited in the number of KIR and HLA class I combinations that could be assayed at any one time. The results obtained from this assay can be used to predict the response of NK cell and T cells when their KIR recognize HLA class I.

Can Selective MHC Downregulation Explain the Specificity and Genetic Diversity of NK Cell Receptors?

Natural killer (NK) cells express inhibiting receptors (iNKRs), which specifically bind MHC-I molecules on the surface of healthy cells. When the expression of MHC-I on the cell surface decreases, which might occur during certain viral infections and cancer, iNKRs lose inhibiting signals and the infected cells become target for NK cell activation (missing-self detection). Although the detection of MHC-I deficient cells can be achieved by conserved receptor-ligand interactions, several iNKRs are encoded by gene families with a remarkable genetic diversity, containing many haplotypes varying in gene content and allelic polymorphism. So far, the biological function of this expansion within the NKR cluster has remained poorly understood. Here, we investigate whether the evolution of diverse iNKRs genes can be driven by a specific viral immunoevasive mechanism: selective MHC downregulation. Several viruses, including EBV, CMV, and HIV, decrease the expression of MHC-I to escape from T cell responses. This downregulation does not always affect all MHC loci in the same way, as viruses target particular MHC molecules. To study the selection pressure of selective MHC downregulation on iNKRs, we have developed an agent-based model simulating an evolutionary scenario of hosts infected with herpes-like viruses, which are able to selectively downregulate the expression of MHC-I molecules on the cell surface. We show that iNKRs evolve specificity and, depending on the similarity of MHC alleles within each locus and the differences between the loci, they can specialize to a particular MHC-I locus. The easier it is to classify an MHC allele to its locus, the lower the required diversity of the NKRs. Thus, the diversification of the iNKR cluster depends on the locus specific MHC structure.

A Coevolutionary Arms Race between Hosts and Viruses Drives Polymorphism and Polygenicity of NK Cell Receptors

Natural killer cell receptors (NKRs) monitor the expression of major histocompatibility class I (MHC-I) and stress molecules to detect unhealthy tissue, such as infected or tumor cells. The NKR gene family shows a remarkable genetic diversity, containing several genes encoding receptors with activating and inhibiting signaling, and varying in gene content and allelic polymorphism. The expansion of the NKR genes is species-specific, with different species evolving alternative expanded NKR genes, which encode structurally different proteins, yet perform comparable functions. So far, the biological function of this expansion within the NKR cluster has remained poorly understood. To study the evolution of NKRs, we have developed an agent-based model implementing a coevolutionary scenario between hosts and herpes-like viruses that are able to evade the immune response by downregulating the expression of MHC-I on the cell surface. We show that hosts evolve specific inhibitory NKRs, specialized to particular MHC-I alleles in the population. Viruses in our simulations readily evolve proteins mimicking the MHC molecules of their host, even in the absence of MHC-I downregulation. As a result, the NKR locus becomes polygenic and polymorphic, encoding both specific inhibiting and activating receptors to optimally protect the hosts from coevolving viruses.

Progression to AIDS in SIV-Infected Rhesus Macaques is Associated with Distinct KIR and MHC class I Polymorphisms and NK Cell Dysfunction

Killer cell immunoglobulin-like receptors (KIR) regulate the activity of natural killer (NK) cells and have been shown to be associated with susceptibility to a number of human infectious diseases. Here, we analyzed NK cell function and genetic associations in a cohort of 52 rhesus macaques experimentally infected with SIVmac and subsequently stratified into high viral load (HVL) and low viral load (LVL) plasma viral loads at set point. This stratification coincided with fast (HVL) and slow (LVL) disease progression indicated by the disease course and critical clinical parameters including CD4+ T cell counts. HVL animals revealed sustained proliferation of NK cells but distinct loss of peripheral blood NK cell numbers and lytic function. Genetic analyses revealed that KIR genes 3DL05, 3DS05, and 3DL10 as well as 3DSW08, 3DLW03, and 3DSW09 are correlated, most likely due to underlying haplotypes. SIV-infection outcome associated with presence of transcripts for two inhibitory KIR genes (KIR3DL02, KIR3DL10) and three activating KIR genes (KIR3DSW08, KIR3DS02, KIR3DS05). Presence of KIR3DL02 and KIR3DSW08 was associated with LVL outcome, whereas presence of KIR3DS02 was associated with HVL outcome. Furthermore, we identified epistasis between KIR and MHC class I alleles as the transcript presence of the correlated genes KIR3DL05, KIR3DS05, and KIR3DL10 increased HVL risk when Mamu-B*012 transcripts were also present or when Mamu-A1*001 transcripts were absent. These genetic associations were mirrored by changes in the numbers, the level of proliferation, and lytic capabilities of NK cells as well as overall survival time and gastro-intestinal tissue viral load.

Old world monkeys and new age science: the evolution of nonhuman primate systems virology

Nonhuman primate (NHP) biomedical models are critical to our understanding of human health and disease, yet we are still in the early stages of developing sufficient tools to support primate genomic research that allow us to better understand the basis of phenotypic traits in NHP models of disease. A mere 7 years ago, the limited NHP transcriptome profiling that was being performed was done using complementary DNA arrays based on human genome sequences, and the lack of NHP genomic information and immunologic reagents precluded the use of NHPs in functional genomic studies. Since then, significant strides have been made in developing genomics capabilities for NHP research, from the rhesus macaque genome sequencing project to the construction of the first macaque-specific high-density oligonucleotide microarray, paving the way for further resource development and additional primate sequencing projects. Complete published draft genome sequences are now available for the chimpanzee ( Chimpanzee Sequencing Analysis Consortium 2005), bonobo ( Prufer et al. 2012), gorilla ( Scally et al. 2012), and baboon ( Ensembl.org 2013), along with the recently completed draft genomes for the cynomolgus macaque and Chinese rhesus macaque. Against this backdrop of both expanding sequence data and the early application of sequence-derived DNA microarrays tools, we will contextualize the development of these community resources and their application to infectious disease research through a literature review of NHP models of acquired immune deficiency syndrome and models of respiratory virus infection. In particular, we will review the use of -omics approaches in studies of simian immunodeficiency virus and respiratory virus pathogenesis and vaccine development, emphasizing the acute and innate responses and the relationship of these to the course of disease and to the evolution of adaptive immunity.

KIR3DL01 recognition of Bw4 ligands in the rhesus macaque: maintenance of Bw4 specificity since the divergence of apes and Old World monkeys

The identification of MHC class I ligands for rhesus macaque killer cell Ig-like receptors (KIRs) is fundamental to our basic understanding of KIR and MHC class I coevolution and to the study of NK cell responses in this nonhuman primate model for AIDS and other viral diseases. In this study, we show that Mamu-KIR3DL01, which is expressed by ∼90% of rhesus macaques, recognizes MHC class I molecules with a Bw4 motif. Primary NK cells expressing Mamu-KIR3DL01 were identified by staining with a mAb which, in this study, was shown to bind Mamu-KIR3DL01 allotypes with an aspartic acid at position 233. The cytolytic activity of Mamu-KIR3DL01(+) NK cells was suppressed by cell lines expressing the Bw4 molecules Mamu-B*007:01, -B*041:01, -B*058:02, and -B*065:01. The Bw4 motif was necessary for Mamu-KIR3DL01 recognition because substitutions in this region abrogated Mamu-KIR3DL01(+) NK cell inhibition. However, the presence of a Bw4 motif was not sufficient for recognition because another Bw4 molecule, Mamu-B*017:01, failed to suppress the cytolytic activity of these NK cells. Replacement of three residues in Mamu-B*017:01, predicted to be KIR contacts based on the three-dimensional structure of the human KIR3DL1-HLA-Bw4 complex, with the corresponding residues at these positions for the other Mamu-Bw4 ligands restored Mamu-KIR3DL01(+) NK cell inhibition. These results define the ligand specificity of one of the most polymorphic and commonly expressed KIRs in the rhesus macaque and reveal similarities in Bw4 recognition by Mamu-KIR3DL01 and human KIR3DL1, despite the absence of an orthologous relationship between these two KIRs or conservation of surface residues predicted to interact with MHC class I ligands.

Quantifying the Protection of Activating and Inhibiting NK Cell Receptors during Infection with a CMV-Like Virus

The responsiveness of natural killer (NK) cells is controlled by balancing signals from activating and inhibitory receptors. The most important ligands of inhibitory NK cell receptors are the highly polymorphic major histocompatibility complex (MHC) class I molecules, which allow NK cells to screen the cellular health of target cells. Although these inhibitory receptor–ligand interactions have been well characterized, the ligands for most activating receptors are still unknown. The mouse cytomegalovirus (MCMV) represents a helpful model to study NK cell-driven immune responses. Many studies have demonstrated that CMV infection can be controlled by NK cells via their activating receptors, but the exact contribution of the different signaling potential (i.e., activating vs. inhibiting) remains puzzling. In this study, we have developed a probabilistic model, which predicts the optimal specificity of inhibitory and activating NK cell receptors needed to offer the best protection against a CMV-like virus. We confirm our analytical predictions with an agent-based model of an evolving host population. Our analysis quantifies the degree of protection of each receptor type, revealing that mixed haplotypes (i.e., haplotypes composed of activating and inhibiting receptors) are most protective against CMV-like viruses, and that the protective effect depends on the number of MHC loci per individual.

Identification and diversity of killer cell Ig-like receptors in Aotus vociferans, a New World monkey

Previous BAC clone analysis of the Platyrrhini owl monkey KIRs have shown an unusual genetic structure in some loci. Therefore, cDNAs encoding KIR molecules from eleven Aotus vociferans monkeys were characterized here; ten putative KIR loci were found, some of which encoded atypical proteins such as KIR4DL and transcripts predicted to encode a D0+D1 configuration (AOTVOKIR2DL1*01v1) which appear to be unique in the Aotus genus. Furthermore, alternative splicing was found as a likely mechanism for producing activator receptors in A. vociferans species. KIR proteins from New World monkeys may be split into three new lineages according to domain by domain phylogenetic analysis. Although the A. vociferans KIR family displayed a high divergence among paralogous genes, individual loci were limited in their genetic polymorphism. Selection analysis showed that both constrained and rapid evolution may operate within the AvKIR family. The frequent alternative splicing (as a likely mechanism generating activator receptors), the presence of KIR4DL and KIR2DL1 (D0+D1) molecules and other data reported here suggest that the KIR family in Aotus has had a rapid evolution, independent from its Catarrhini counterparts.

Bimodal evolution of the killer cell Ig-like receptor (KIR) family in New World primates

The immunoglobulin-like receptor (KIR) gene family in New World primates (Platyrrhini) has been characterized only in the owl monkey (Aotus sp.). To gain a better understanding of the KIR system in Platyrrhini, we analyzed a KIR haplotype in Ateles geoffroyi, and sequenced KIR complementary DNAs (cDNAs) from other three Atelidae species, Ateles hybridus, Ateles belzebuth, and Lagothrix lagotricha. Atelidae expressed a variable set of activating and inhibitory KIRs that diversified independently from their Catarrhini counterparts. They had a unique mechanism to generate activating receptors from inhibitory ones, involving a single nucleotide deletion in exon 7 and a change in the donor splice site of intron 7. The A. geoffroyi haplotype contained at least six gene models including a pseudogene, two coding inhibitory receptors, and three coding activating receptors. The centromeric region was in a tail-to-tail orientation with respect to the telomeric region. The owl monkey KIR haplotype shared this organization, and in phylogenetic trees, the centromeric genes clustered together with those of A. geoffroyi, whereas their telomeric genes clustered independently. KIR cDNAs from the other Atelidae species conformed to this pattern. Signatures of positive selection were found in residues predicted to interact with the major histocompatibility complex. Such signatures, however, primarily explained variability between paralogous genes but not between alleles in a locus. Atelidae, therefore, has expanded the KIR family in a bimodal fashion, where an inverted centromeric region has remained relatively conserved and the telomeric region has diversified by a rapid process of gene duplication and divergence, likely favored by positive selection for ligand binding.

Extensive variation in gene copy number at the killer immunoglobulin-like receptor locus in humans

Killer immunoglobulin-like receptors (KIRs) are involved in the regulation of natural killer cell cytotoxicity. Within the human genome seventeen KIR genes are present, which all contain a large number of allelic variants. The high level of homology among KIR genes has hampered KIR genotyping in larger cohorts, and determination of gene copy number variation (CNV) has been difficult. We have designed a multiplex ligation-dependent probe amplification (MLPA) technique for genotyping and CNV determination in one single assay and validated the results by next-generation sequencing and with a KIR gene-specific short tandem repeat assay. In this way, we demonstrate in a cohort of 120 individuals a high level of CNV for all KIR genes except for the framework genes KIR3DL3 and KIR3DL2. Application of our MLPA assay in segregation analyses of families from the Centre d'Etude du Polymorphisme Humaine, previously KIR-genotyped by classical techniques, confirmed an earlier reported duplication and resulted in the identification of a novel duplication event in one of these families. In summary, our KIR MLPA assay allows rapid and accurate KIR genotyping and CNV detection, thus rendering improved transplantation programs and oncology treatment feasible, and enables more detailed studies on the role of KIRs in human (auto)immunity and infectious disease.

Expression patterns of killer cell immunoglobulin-like receptors (KIR) of NK-cell and T-cell subsets in Old World monkeys

The expression of killer cell immunoglobulin-like receptors (KIR) on lymphocytes of rhesus macaques and other Old World monkeys was unknown so far. We used our recently established monoclonal anti-rhesus macaque KIR antibodies in multicolour flow cytometry for phenotypic characterization of KIR protein expression on natural killer (NK) cells and T cell subsets of rhesus macaques, cynomolgus macaques, hamadryas baboons, and African green monkeys. Similar to human KIR, we found clonal expression patterns of KIR on NK and T cell subsets in rhesus macaques and differences between individuals using pan-KIR3D antibody 1C7 and antibodies specific for single KIR. Similar results were obtained with lymphocytes from the other studied species. Notably, African green monkeys show only a low frequency of KIR3D expressed on CD8⁺ αβT cells. Contrasting human NK cells are KIR-positive CD56bright NK cells and frequencies of KIR-expressing NK cells that are independent of the presence of their cognate MHC class I ligands in rhesus macaques. Interestingly, the frequency of KIR-expressing cells and the expression strength of KIR3D are correlated in γδ T cells of rhesus macaques and CD8⁺ αβT cells of baboons.

The Turkey Ig-like receptor family: identification, expression and function

The chicken leukocyte receptor complex located on microchromosome 31 encodes the chicken Ig-like receptors (CHIR), a vastly expanded gene family which can be further divided into three subgroups: activating CHIR-A, bifunctional CHIR-AB and inhibitory CHIR-B. Here, we investigated the presence of CHIR homologues in other bird species. The available genome databases of turkey, duck and zebra finch were screened with different strategies including BLAST searches employing various CHIR sequences, and keyword searches. We could not identify CHIR homologues in the distantly related zebra finch and duck, however, several partial and complete sequences of CHIR homologues were identified on chromosome 3 of the turkey genome. They were designated as turkey Ig-like receptors (TILR). Using cDNA derived from turkey blood and spleen RNA, six full length TILR could be amplified and further divided according to the typical sequence features into one activating TILR-A, one inhibitory TILR-B and four bifunctional TILR-AB. Since the TILR-AB sequences all displayed the critical residues shown to be involved in binding to IgY, we next confirmed the IgY binding using a soluble TILR-AB1-huIg fusion protein. This fusion protein reacted with IgY derived from various gallinaceous birds, but not with IgY from other bird species. Finally, we tested various mab directed against CHIR for their crossreactivity with either turkey or duck leukocytes. Whereas no staining was detectable with duck cells, the CHIR-AB1 specific mab 8D12 and the CHIR-A2 specific mab 13E2 both reacted with a leukocyte subpopulation that was further identified as thrombocytes by double immunofluorescence employing B-cell, T-cell and thrombocyte specific reagents. In summary, although the turkey harbors similar LRC genes as the chicken, their distribution seems to be distinct with predominance on thrombocytes rather than lymphocytes.

KIR2DL5: An Orphan Inhibitory Receptor Displaying Complex Patterns of Polymorphism and Expression

A recently developed anti-KIR2DL5 (CD158f) antibody has demonstrated KIR2DL5 expression on the surface of NK and T lymphocytes, making it the last functional KIR identified in the human genome. KIR2DL5 belongs to an ancestral lineage of KIR with Ig-like domains of the D0-D2 type, of which KIR2DL4, an HLA-G receptor, is the only other human member. Despite KIR2DL4 and KIR2DL5 being encoded by genes with similar domain usage, several KIR2DL5 functions resemble more closely those of KIR recognizing classical HLA class I molecules - surface-expressed KIR2DL5 inhibits NK cells through the SHP-2 phosphatase and displays a clonal distribution on NK and T lymphocytes. No activating homolog of KIR2DL5 has been described in any species. The genetics of KIR2DL5 is complicated by duplication of its gene in an ancestor of modern humans living ∼1.7 million years ago. Both KIR2DL5 paralogs have undergone allelic diversification; the centromeric gene is most often represented by alleles whose expression is silenced epigenetically through DNA methylation, thus providing a natural system to investigate the regulation of KIR transcription. The role of KIR2DL5 in immunity is not completely understood, in spite of different attempts to define its ligand. Here we revisit the most relevant characteristics of KIR2DL5, an NK-cell receptor possessing a unique combination of genetic, structural, and functional features.

Mutation at positively selected positions in the binding site for HLA-C shows that KIR2DL1 is a more refined but less adaptable NK cell receptor than KIR2DL3

Through recognition of HLA class I, killer cell Ig-like receptors (KIR) modulate NK cell functions in human immunity and reproduction. Although a minority of HLA-A and -B allotypes are KIR ligands, HLA-C allotypes dominate this regulation, because they all carry either the C1 epitope recognized by KIR2DL2/3 or the C2 epitope recognized by KIR2DL1. The C1 epitope and C1-specific KIR evolved first, followed several million years later by the C2 epitope and C2-specific KIR. Strong, varying selection pressure on NK cell functions drove the diversification and divergence of hominid KIR, with six positions in the HLA class I binding site of KIR being targets for positive diversifying selection. Introducing each naturally occurring residue at these positions into KIR2DL1 and KIR2DL3 produced 38 point mutants that were tested for binding to 95 HLA- A, -B, and -C allotypes. Modulating specificity for HLA-C is position 44, whereas positions 71 and 131 control cross-reactivity with HLA-A*11:02. Dominating avidity modulation is position 70, with lesser contributions from positions 68 and 182. KIR2DL3 has lower avidity and broader specificity than KIR2DL1. Mutation could increase the avidity and change the specificity of KIR2DL3, whereas KIR2DL1 specificity was resistant to mutation, and its avidity could only be lowered. The contrasting inflexibility of KIR2DL1 and adaptability of KIR2DL3 fit with C2-specific KIR having evolved from C1-specific KIR, and not vice versa. Substitutions restricted to activating KIR all reduced the avidity of KIR2DL1 and KIR2DL3, further evidence that activating KIR function often becomes subject to selective attenuation.

Human-specific evolution of killer cell immunoglobulin-like receptor recognition of major histocompatibility complex class I molecules

In placental mammals, natural killer (NK) cells are a population of lymphocytes that make unique contributions to immune defence and reproduction, functions essential for survival of individuals, populations and species. Modulating these functions are conserved and variable NK-cell receptors that recognize epitopes of major histocompatibility complex (MHC) class I molecules. In humans, for example, recognition of human leucocyte antigen (HLA)-E by the CD94:NKG2A receptor is conserved, whereas recognition of HLA-A, B and C by the killer cell immunoglobulin-like receptors (KIRs) is diversified. Competing demands of the immune and reproductive systems, and of T-cell and NK-cell immunity-combined with the segregation on different chromosomes of variable NK-cell receptors and their MHC class I ligands-drive an unusually rapid evolution that has resulted in unprecedented levels of species specificity, as first appreciated from comparison of mice and humans. Counterparts to human KIR are present only in simian primates. Observed in these species is the coevolution of KIR and the four MHC class I epitopes to which human KIR recognition is restricted. Unique to hominids is the emergence of the MHC-C locus as a supplier of specialized and superior ligands for KIR. This evolutionary trend is most highly elaborated in the chimpanzee. Unique to the human KIR locus are two groups of KIR haplotypes that are present in all human populations and subject to balancing selection. Group A KIR haplotypes resemble chimpanzee KIR haplotypes and are enriched for genes encoding KIR that bind HLA class I, whereas group B KIR haplotypes are enriched for genes encoding receptors with diminished capacity to bind HLA class I. Correlating with their balance in human populations, B haplotypes favour reproductive success, whereas A haplotypes favour successful immune defence. Evolution of the B KIR haplotypes is thus unique to the human species.

Characterization of killer immunoglobulin-like receptor genetics and comprehensive genotyping by pyrosequencing in rhesus macaques

Background

Human killer immunoglobulin-like receptors (KIRs) play a critical role in governing the immune response to neoplastic and infectious disease. Rhesus macaques serve as important animal models for many human diseases in which KIRs are implicated; however, the study of KIR activity in this model is hindered by incomplete characterization of KIR genetics.

Results

Here we present a characterization of KIR genetics in rhesus macaques (Macaca mulatta). We conducted a survey of KIRs in this species, identifying 47 novel full-length KIR sequences. Using this expanded sequence library to build upon previous work, we present evidence supporting the existence of 22 Mamu-KIR genes, providing a framework within which to describe macaque KIRs. We also developed a novel pyrosequencing-based technique for KIR genotyping. This method provides both comprehensive KIR genotype and frequency estimates of transcript level, with implications for the study of KIRs in all species.

Conclusions

The results of this study significantly improve our understanding of macaque KIR genetic organization and diversity, with implications for the study of many human diseases that use macaques as a model. The ability to obtain comprehensive KIR genotypes is of basic importance for the study of KIRs, and can easily be adapted to other species. Together these findings both advance the field of macaque KIRs and facilitate future research into the role of KIRs in human disease.

Rhesus macaque KIR bind human MHC class I with broad specificity and recognize HLA-C more effectively than HLA-A and HLA-B

Human killer cell immunoglobulin-like receptors (KIR) recognize A3/11, Bw4, C1, and C2 epitopes carried by mutually exclusive subsets of human leukocyte antigen (HLA)-A, -B, and -C allotypes. Chimpanzee and orangutan have counterparts to HLA-A, -B, and -C, and KIR that recognize the A3/11, Bw4, C1, and C2 epitopes, either individually or in combination. Because rhesus macaque has counterparts of HLA-A and -B, but not HLA-C, we expected that rhesus KIR would better recognize HLA-A and -B, than HLA-C. Comparison of the interactions of nine rhesus KIR3D with 95 HLA isoforms, showed the KIR have broad specificity for HLA-A, -B, and -C, but vary in avidity. Considering both the strength and breadth of reaction, HLA-C was the major target for rhesus KIR, followed by HLA-B, then HLA-A. Strong reactions with HLA-A were restricted to the minority of allotypes carrying the Bw4 epitope, whereas strong reactions with HLA-B partitioned between allotypes having and lacking Bw4. Contrasting to HLA-A and -B, every HLA-C allotype bound to the nine rhesus KIR. Sequence comparison of high- and low-binding HLA allotypes revealed the importance of polymorphism in the helix of the α(1) domain and the peptide-binding pockets. At peptide position 9, nonpolar residues favor binding to rhesus KIR, whereas charged residues do not. Contrary to expectation, rhesus KIR bind more effectively to HLA-C, than to HLA-A and -B. This property is consistent with major histocompatibility complex (MHC)-C having evolved in hominids to be a generally superior ligand for KIR than MHC-A and MHC-B.