Human and Rhesus Macaque KIR Haplotypes Defined by Their Transcriptomes

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.

HLA/MHC and KIR characterization in humans and non-human primates using Oxford Nanopore Technologies and Pacific Biosciences sequencing platforms

The gene products of the HLA/MHC and KIR multigene families are important modulators of the immune system and are associated with health and disease. Characterization of the genes encoding these receptors has been integrated into different biomedical applications, including transplantation and reproduction biology, immune therapies and in fundamental research into disease susceptibility or resistance. Conventional short-read sequencing strategies have shown their value in high throughput typing, but are insufficient to uncover the entire complexity of the highly polymorphic HLA/MHC and KIR gene systems. The implementation of single-molecule and real-time sequencing platforms, offered by Pacific Biosciences (PacBio) and Oxford Nanopore Technologies (ONT), revolutionized the fields of genomics and transcriptomics. Using fundamentally distinct principles, these platforms generate long-read data that can unwire the plasticity of the HLA/MHC and KIR genes, including high-resolution characterization of genes, alleles, phased haplotypes, transcription levels and epigenetics modification patterns. These insights might have profound clinical relevance, such as improved matching of donors and patients in clinical transplantation, but could also lift disease association studies to a higher level. Even more, a comprehensive characterization may refine animal models in preclinical studies. In this review, the different HLA/MHC and KIR characterization approaches using PacBio and ONT platforms are described and discussed.

Comparative genetics of KIR haplotype diversity in humans and rhesus macaques: the balancing act

The role of natural killer (NK) cells is tightly modulated by interactions of killer cell immunoglobulin-like receptors (KIR) with their ligands of the MHC class I family. Several characteristics of the KIR gene products are conserved in primate evolution, like the receptor structures and the variegated expression pattern. At the genomic level, however, the clusters encoding the KIR family display species-specific diversity, reflected by differential gene expansions and haplotype architecture. The human KIR cluster is extensively studied in large cohorts from various populations, which revealed two KIR haplotype groups, A and B, that represent more inhibitory and more activating functional profiles, respectively. So far, genomic KIR analyses in large outbred populations of non-human primate species are lacking. In this study, we roughly quadrupled the number of rhesus macaques studied for their KIR transcriptome (n = 298). Using segregation analysis, we defined 112 unique KIR region configurations, half of which display a more inhibitory profile, whereas the other half has a more activating potential. The frequencies and functional potential of these profiles might mirror the human KIR haplotype groups. However, whereas the human group A and B KIR haplotypes are confined to largely fixed organizations, the haplotypes in macaques feature highly variable gene content. Moreover, KIR homozygosity was hardly encountered in this panel of macaques. This study exhibits highly diverse haplotype architectures in humans and macaques, which nevertheless might have an equivalent effect on the modulation of NK cell activity.

The Genomic Organization of the LILR Region Remained Largely Conserved Throughout Primate Evolution: Implications for Health And Disease

The genes of the leukocyte immunoglobulin-like receptor (LILR) family map to the leukocyte receptor complex (LRC) on chromosome 19, and consist of both activating and inhibiting entities. These receptors are often involved in regulating immune responses, and are considered to play a role in health and disease. The human LILR region and evolutionary equivalents in some rodent and bird species have been thoroughly characterized. In non-human primates, the LILR region is annotated, but a thorough comparison between humans and non-human primates has not yet been documented. Therefore, it was decided to undertake a comprehensive comparison of the human and non-human primate LILR region at the genomic level. During primate evolution the organization of the LILR region remained largely conserved. One major exception, however, is provided by the common marmoset, a New World monkey species, which seems to feature a substantial contraction of the number of LILR genes in both the centromeric and the telomeric region. Furthermore, genomic analysis revealed that the killer-cell immunoglobulin-like receptor gene KIR3DX1, which maps in the LILR region, features one copy in humans and great ape species. A second copy, which might have been introduced by a duplication event, was observed in the lesser apes, and in Old and New World monkey species. The highly conserved gene organization allowed us to standardize the LILR gene nomenclature for non-human primate species, and implies that most of the receptors encoded by these genes likely fulfill highly preserved functions.

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.

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.

Characterization of 100 extended major histocompatibility complex haplotypes in Indonesian cynomolgus macaques

Many medical advancements—including improvements to anti-rejection therapies in transplantation and vaccine development—rely on preclinical studies conducted in cynomolgus macaques (Macaca fascicularis). Major histocompatibility complex (MHC) class I and class II genes of cynomolgus macaques are orthologous to human leukocyte antigen complex (HLA) class I and class II genes, respectively. Both encode cell-surface proteins involved in cell recognition and rejection of non-host tissues. MHC class I and class II genes are highly polymorphic, so comprehensive genotyping requires the development of complete databases of allelic variants. Our group used PacBio circular consensus sequencing of full-length cDNA amplicons to characterize MHC class I and class II transcript sequences for a cohort of 293 Indonesian cynomolgus macaques (ICM) in a large, pedigreed breeding colony. These studies allowed us to expand the existing database of Macaca fascicularis (Mafa) alleles by identifying an additional 141 MHC class I and 61 class II transcript sequences. In addition, we defined co-segregating combinations of allelic variants as regional haplotypes for 70 Mafa-A, 78 Mafa-B, and 45 Mafa-DRB gene clusters. Finally, we defined class I and class II transcripts that are associated with 100 extended MHC haplotypes in this breeding colony by combining our genotyping analyses with short tandem repeat (STR) patterns across the MHC region. Our sequencing analyses and haplotype definitions improve the utility of these ICM for transplantation studies as well as infectious disease and vaccine research.

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.

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.