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

The 2024 IPD-MHC database update: a comprehensive resource for major histocompatibility complex studies

The IPD-MHC Database project (http://www.ebi.ac.uk/ipd/mhc/) serves as a comprehensive and expertly curated repository for major histocompatibility complex (MHC) sequences from non-human species, providing the necessary infrastructure and tools to study the function and evolution of this highly polymorphic genomic region. In its latest version, the IPD-MHC database has expanded both in content and in the tools for data visualization and comparison. The database now hosts over 18 000 MHC alleles from 125 species, organized into eleven taxonomic groups, all manually curated and named by the Comparative MHC Nomenclature Committee. A cetacean section has recently been included, offering researchers valuable data to study the immune system of whales, dolphins, and porpoises, as well establishing the official nomenclature platform for the Cetacea Leukocyte Antigens (CeLA). In response to user demand and reflecting broader trends in bioinformatics and immunogenetics, IPD-MHC now includes the predicted tertiary structure of over 8000 alleles and allows comparison and visualisation of allele variation within and between species at single residue resolution. These latest developments maintain the critically important link between official nomenclature of curated alleles and the ability to analyse this complex polymorphism using the most up to date methods within a single repository.

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.

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.

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.

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.

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.