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.

MIC gene polymorphism and haplotype diversity in rhesus macaques

Rhesus macaques (Macaca mulatta) mainly originating from India were analysed for their major histocompatibility complex class I-related (MIC) gene repertoire. Thus far, three distinct genes, designated MIC1, MIC2 and MIC3, have been identified in the rhesus macaque. In addition, an MICD pseudogene has been described mapping apart from the other loci in a telomeric direction. Genomic comparisons and the presence of a characteristic microsatellite in exon 5 suggest that the MIC1 gene is the equivalent of the human MICA gene. Hence, the MIC2 gene, lacking the microsatellite - as do humans -, is considered to be the equivalent of human MICB. The MIC3 gene, a hybrid of MICA and MICB, seems to be generated by a crossing-over event with one breakpoint in intron 3 and accordingly is named MICA/B. Apart from their human counterparts, MICA, MICB and MICA/B cluster in separate branches in the phylogenetic tree, confirming the hybrid character of the MICA/B gene. Population analyses have shown that the various genes display polymorphism, and six MICA, five MICB and three MICA/B alleles have been identified. In the panel of homozygous typing cells, two distinct haplotype configurations have been defined by segregation analyses. Each haplotype comprises an MICB gene in conjunction with either an MICA or an MICA/B gene. Furthermore, the presence of a polymorphic microsatellite in the MICA and MICA/B alleles facilitates speedy and accurate haplotyping.