Meiotic recombination generates rich diversity in NK cell receptor genes, alleles, and haplotypes

High-resolution KIR and HLA genotyping in three Chinese ethnic minorities reveals distinct origins

Polymorphism of killer-cell immunoglobulin-like receptors (KIRs) and their HLA class I ligands impacts the effector activity of cytotoxic NK cell and T cell subsets. Therefore, understanding the extent and implications of KIR and HLA class I genetic polymorphism across various populations is important for immunological and medical research. In this study, we conducted a high-resolution investigation of KIR and HLA class I diversity in three distinct Chinese ethnic minority populations. We studied the She, Yugur, and Tajik, and compared them with the Zhejiang Han population (Zhe), which represents the majority Southern Han ethnicity. Our findings revealed that the Tajik population exhibited the most diverse KIR copy number, allele, and haplotype diversity among the four populations. This diversity aligns with their proposed ancestral origin, closely resembling that of Iranian populations, with a relatively higher presence of KIR-B genes, alleles, and haplotypes compared with the other Chinese populations. The Yugur population displayed KIR distributions similar to those of the Tibetans and Southeast Asians, whereas the She population resembled the Zhe and other East Asians, as confirmed by genetic distance analysis of KIR. Additionally, we identified 12.9% of individuals across the three minority populations as having KIR haplotypes characterized by specific gene block insertions or deletions. Genetic analysis based on HLA alleles yielded consistent results, even though there were extensive variations in HLA alleles. The observed variations in KIR interactions, such as higher numbers of 2DL1-C2 interactions in Tajik and Yugur populations and of 2DL3-C1 interactions in the She population, are likely shaped by demographic and evolutionary mechanisms specific to their local environments. Overall, our findings offer valuable insights into the distribution of KIR and HLA diversity among three distinct Chinese ethnic minority populations, which can inform future clinical and population studies.

Determination for KIR genotype and allele copy number via real-time quantitative PCR method

Killer cell immunoglobulin-like receptor (KIR) and human leukocyte antigen (HLA) play crucial roles in regulating NK cell activity. Here, we report a real-time quantitative PCR (qPCR) to genotype all KIR genes and their copy numbers simultaneously. With 18 pairs of locus-specific primers, we identified KIR genes by Ct values and determined KIR copy number using the 2-∆Ct method. Haplotypes were assigned based on KIR gene copy numbers. The real-time qPCR results were consistent with the NGS method, except for one sample with KIR2DL5 discrepancy. qPCR is a multiplex method that can identify KIR copy number, which helps obtain a relatively accurate haplotype structure, facilitating increased KIR research in laboratories where NGS or other high-resolution methods are not available.

Geny: A Genotyping Tool for Allelic Decomposition of Killer Cell Immunoglobulin-Like Receptor Genes

Accurate genotyping of Killer cell Immunoglobulin-like Receptor (KIR) genes plays a pivotal role in enhancing our understanding of innate immune responses, disease correlations, and the advancement of personalized medicine. However, due to the high variability of the KIR region and high level of sequence similarity among different KIR genes, the currently available genotyping methods are unable to accurately infer copy numbers, genotypes and haplotypes of individual KIR genes from next-generation sequencing data. Here we introduce Geny, a new computational tool for precise genotyping of KIR genes. Geny utilizes available KIR haplotype databases and proposes a novel combination of expectation-maximization filtering schemes and integer linear programming-based combinatorial optimization models to resolve ambiguous reads, provide accurate copy number estimation and estimate the haplotype of each copy for the genes within the KIR region. We evaluated Geny on a large set of simulated short-read datasets covering the known validated KIR region assemblies and a set of Illumina short-read samples sequenced from 25 validated samples from the Human Pangenome Reference Consortium collection and showed that it outperforms the existing genotyping tools in terms of accuracy, precision and recall. We envision Geny becoming a valuable resource for understanding immune system response and consequently advancing the field of patient-centric medicine.

Exceptional diversity of KIR and HLA class I in Egypt

Genetically determined variation of killer cell immunoglobulin like receptors (KIR) and their HLA class I ligands affects multiple aspects of human health. Their extreme diversity is generated through complex interplay of natural selection for pathogen resistance and reproductive health, combined with demographic structure and dispersal. Despite significant importance to multiple health conditions of differential effect across populations, the nature and extent of immunogenetic diversity is under-studied for many geographic regions. Here, we describe the first high-resolution analysis of KIR and HLA class I combinatorial diversity in Northern Africa. Analysis of 125 healthy unrelated individuals from Cairo in Egypt yielded 186 KIR alleles arranged in 146 distinct centromeric and 79 distinct telomeric haplotypes. The most frequent haplotypes observed were KIR-A, encoding two inhibitory receptors specific for HLA-C, two that are specific for HLA-A and -B, and no activating receptors. Together with 141 alleles of HLA class I, 75 of which encode a KIR ligand, we identified a mean of six distinct interacting pairs of inhibitory KIR and HLA allotypes per individual. We additionally characterize 16 KIR alleles newly identified in the study population. Our findings place Egyptians as one of the most highly diverse populations worldwide, with important implications for transplant matching and studies of immune-mediated diseases.

KIR genotype and haplotype repertoire in Kuwaiti healthy donors, hematopoietic cell transplant recipients and healthy family members

The gene complex located on chromosome 19q13.4 encodes the Killer-cell Immunoglobulin-like Receptors (KIRs), which exhibit remarkable polymorphism in both gene content and sequences. Further, the repertoire of KIR genes varies within and between populations, creating a diverse pool of KIR genotypes. This study was carried out to characterize KIR genotypes and haplotypes among 379 Arab Kuwaiti individuals including 60 subjects from 20 trio families, 49 hematopoietic cell transplantation (HCT) recipients and 270 healthy Kuwaiti volunteer HCT donors. KIR Genotyping was performed by a combination of reverse sequence specific oligonucleotide probes (rSSO) and/or Real Time PCR. The frequencies of KIR genes in 270 healthy Kuwaiti volunteer donors were compared to previously reported frequencies in other populations. In addition, we compared the differences in KIR repertoire of patients and healthy donors to investigate the reproducibility of previously reported significant differences between patients with hematological malignancies and healthy donors. The observed frequencies in our cohort volunteer HCT donors was comparable to those reported in neighboring Arab populations. The activating genes KIR2DS1, KIR2DS5 and KIR3DS1 and the inhibitory gene KIR2DL5 were significantly more frequent in patients compared to healthy donors, however, none of the previously reported differences were reproducible in our Kuwaiti cohort. This report is the first description of KIR gene carrier frequency and haplotype characterization in a fairly large cohort of the Kuwaiti population, which may have implications in KIR based HCT donor selection strategies.

KIR-HLA gene diversities and susceptibility to lung cancer

Killer-cell immunoglobulin-like receptors (KIR) are essential for acquiring natural killer (NK) cell effector function, which is modulated by a balance between the net input of signals derived from inhibitory and activating receptors through engagement by human leukocyte antigen (HLA) class I ligands. KIR and HLA loci are polygenic and polymorphic and exhibit substantial variation between individuals and populations. We attempted to investigate the contribution of KIR complex and HLA class I ligands to the genetic predisposition to lung cancer in the native population of southern Iran. We genotyped 16 KIR genes for a total of 232 patients with lung cancer and 448 healthy controls (HC), among which 85 patients and 178 HCs were taken into account for evaluating combined KIR-HLA associations. KIR2DL2 and 2DS2 were increased significantly in patients than in controls, individually (OR 1.63, and OR 1.42, respectively) and in combination with HLA-C1 ligands (OR 1.99, and OR 1.93, respectively). KIR3DS1 (OR 0.67) and 2DS1 (OR 0.69) were more likely presented in controls in the absence of their relative ligands. The incidence of CxTx subset was increased in lung cancer patients (OR 1.83), and disease risk strikingly increased by more than fivefold among genotype ID19 carriers (a CxTx genotype that carries 2DL2 in the absence of 2DS2, OR 5.92). We found that genotypes with iKIRs > aKIRs (OR 1.67) were more frequently presented in lung cancer patients. Additionally, patients with lung cancer were more likely to carry the combination of CxTx/2DS2 compared to controls (OR 2.04), and iKIRs > aKIRs genotypes in the presence of 2DL2 (OR 2.05) increased the likelihood of lung cancer development. Here we report new susceptibility factors and the contribution of KIR and HLA-I encoding genes to lung cancer risk, highlighting an array of genetic effects and disease setting which regulates NK cell responsiveness. Our results suggest that inherited KIR genes and HLA-I ligands specifying the educational state of NK cells can modify lung cancer risk.

Immunogenomics of Killer Cell Immunoglobulin-Like Receptor (KIR) and HLA Class I: Coevolution and Consequences for Human Health

Interactions of killer cell immunoglobin-like receptors (KIR) with human leukocyte antigens (HLA) class I regulate effector functions of key cytotoxic cells of innate and adaptive immunity. The extreme diversity of this interaction is genetically determined, having evolved in the ever-changing environment of pathogen exposure. Diversity of KIR and HLA genes is further facilitated by their independent segregation on separate chromosomes. That fetal implantation relies on many of the same types of immune cells as infection control places certain constraints on the evolution of KIR interactions with HLA. Consequently, specific inherited combinations of receptors and ligands may predispose to specific immune-mediated diseases, including autoimmunity. Combinatorial diversity of KIR and HLA class I can also differentiate success rates of immunotherapy directed to these diseases. Progress toward both etiopathology and predicting response to therapy is being achieved through detailed characterization of the extent and consequences of the combinatorial diversity of KIR and HLA. Achieving these goals is more tractable with the development of integrated analyses of molecular evolution, function, and pathology that will establish guidelines for understanding and managing risks. Here, we present what is known about the coevolution of KIR with HLA class I and the impact of their complexity on immune function and homeostasis.

Epitope characterization of a monoclonal antibody that selectively recognizes KIR2DL1 allotypes

Killer immunoglobulin‐like receptor (KIR) genes code for a family of inhibitory and activating receptors, finely tuning NK cell function. Numerous studies reported the relevance of KIR allelic polymorphism on KIR expression, ligand affinity, and strength in signal transduction. Although KIR variability, including gene copy number and allelic polymorphism, in combination with HLA class I polymorphism, impacts both KIR expression and NK cell education, only a precise phenotypic analysis can define the size of the different KIR^pos NK cell subsets. In this context, reagents recognizing a limited number of KIRs is essential. In this study, we have characterized the specificity of an anti‐KIR mAb termed HP‐DM1. Testing its binding to HEK‐293T cells transfected with plasmids coding for different KIRs, we demonstrated that HP‐DM1 mAb exclusively reacts with KIR2DL1. Using site‐directed mutagenesis, we identified the four amino acids relevant for HP‐DM1 recognition: M44, S67, R68, and T70. HP‐DM1 mAb binds to a conformational epitope including M44, the residue crucial for HLA‐C K80 recognition by KIR2DL1. Based on the HP‐DM1 epitope characterization, we could extend its reactivity to all KIR2DL1 allotypes identified except for KIR2DL1*022 and, most likely, KIR2DL1*020, predicting that it does not recognize any other KIR with the only exception of KIR2DS1*013. Moreover, by identifying the residues relevant for HP‐DM1 binding, continuously updating of its reactivity will be facilitated.

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.

Allele imputation for the killer cell immunoglobulin-like receptor KIR3DL1/S1

Highly polymorphic interaction of KIR3DL1 and KIR3DS1 with HLA class I ligands modulates the effector functions of natural killer (NK) cells and some T cells. This genetically determined diversity affects severity of infections, immune-mediated diseases, and some cancers, and impacts the course of immunotherapies, including transplantation. KIR3DL1 is an inhibitory receptor, and KIR3DS1 is an activating receptor encoded by the KIR3DL1/S1 gene that has more than 200 diverse and divergent alleles. Determination of KIR3DL1/S1 genotypes for medical application is hampered by complex sequence and structural variation, requiring targeted approaches to generate and analyze high-resolution allele data. To overcome these obstacles, we developed and optimized a model for imputing KIR3DL1/S1 alleles at high-resolution from whole-genome SNP data. We designed the model to represent a substantial component of human genetic diversity. Our Global imputation model is effective at genotyping KIR3DL1/S1 alleles with an accuracy ranging from 88% in Africans to 97% in East Asians, with mean specificity of 99% and sensitivity of 95% for alleles >1% frequency. We used the established algorithm of the HIBAG program, in a modification named Pulling Out Natural killer cell Genomics (PONG). Because HIBAG was designed to impute HLA alleles also from whole-genome SNP data, PONG allows combinatorial diversity of KIR3DL1/S1 with HLA-A and -B to be analyzed using complementary techniques on a single data source. The use of PONG thus negates the need for targeted sequencing data in very large-scale association studies where such methods might not be tractable.

KIR gene content imputation from single-nucleotide polymorphisms in the Finnish population

The killer cell immunoglobulin-like receptor (KIR) gene cluster on chromosome 19 encodes cell surface glycoproteins that bind class I human leukocyte antigen (HLA) molecules as well as some other ligands. Through regulation of natural killer (NK) cell activity KIRs participate in tumour surveillance and clearing viral infections. KIR gene gene copy number variation associates with the outcome of transplantations and susceptibility to immune-mediated diseases. Inferring KIR gene content from genetic variant data is therefore desirable for immunogenetic analysis, particularly in the context of growing biobank genome data collections that rely on genotyping by microarray. Here we describe a stand-alone and freely available gene content imputation for 12 KIR genes. The models were trained using 807 Finnish biobank samples genotyped for 5900 KIR-region SNPs and analysed for KIR gene content with targeted sequencing. Cross-validation results demonstrate a high mean overall accuracy of 98.5% (95% CI [97.0-99.2]%) which compares favourably with previous methods including short-read sequencing based approaches.

Abundant CpG-sequences in human genomes inhibit KIR3DL2-expressing NK cells

Killer Immunoglobulin-like Receptors (KIR) comprise a diverse, highly polymorphic family of cell-surface glycoproteins that are principally expressed by Natural Killer (NK) cells. These innate immune lymphocytes fulfill vital functions in human reproduction and immune responses to viral infection. KIR3DL2 is an inhibitory NK cell receptor that recognizes a common epitope of the HLA-A3 and HLA-A11 class I glycoproteins of the major histocompatibility complex. KIR3DL2 also binds exogenous DNA containing the CpG motif. This interaction causes internalization of the KIR-DNA. Exogenous CpG-DNA typically activates NK cells, but the specificity of KIR3DL2-DNA binding and internalization is unclear. We hypothesized that KIR3DL2 binds exogenous DNA in a sequence-specific manner that differentiates pathogen DNA from self-DNA. In testing this hypothesis, we surveyed octameric CpG-DNA sequences in the human genome, and in reference genomes of all bacteria, fungi, viruses, and parasites, with focus on medically relevant species. Among all pathogens, the nucleotides flanking CpG motifs in the genomes of parasitic worms that infect humans are most divergent from those in the human genome. We cultured KIR3DL2⁺NKL cells with the commonest CpG-DNA sequences in either human or pathogen genomes. DNA uptake was negatively correlated with the most common CpG-DNA sequences in the human genome. These CpG-DNA sequences induced inhibitory signaling in KIR3DL2⁺NKL cells. In contrast, KIR3DL2⁺NKL cells lysed more malignant targets and produced more IFNγ after culture with CpG-DNA sequences prevalent in parasitic worms. By applying functional immunology to evolutionary genomics, we conclude that KIR3DL2 allows NK cells to differentiate self-DNA from pathogen DNA.

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.

High-throughput Interpretation of Killer-cell Immunoglobulin-like Receptor Short-read Sequencing Data with PING

The killer-cell immunoglobulin-like receptor (KIR) complex on chromosome 19 encodes receptors that modulate the activity of natural killer cells, and variation in these genes has been linked to infectious and autoimmune disease, as well as having bearing on pregnancy and transplant outcomes. The medical relevance and high variability of KIR genes makes short-read sequencing an attractive technology for interrogating the region, providing a high-throughput, high-fidelity sequencing method that is cost-effective. However, because this gene complex is characterized by extensive nucleotide polymorphism, structural variation including gene fusions and deletions, and a high level of homology between genes, its interrogation at high resolution has been thwarted by bioinformatic challenges, with most studies limited to examining presence or absence of specific genes. Here, we present the PING (Pushing Immunogenetics to the Next Generation) pipeline, which incorporates empirical data, novel alignment strategies and a custom alignment processing workflow to enable high-throughput KIR sequence analysis from short-read data. PING provides KIR gene copy number classification functionality for all KIR genes through use of a comprehensive alignment reference. The gene copy number determined per individual enables an innovative genotype determination workflow using genotype-matched references. Together, these methods address the challenges imposed by the structural complexity and overall homology of the KIR complex. To determine copy number and genotype determination accuracy, we applied PING to European and African validation cohorts and a synthetic dataset. PING demonstrated exceptional copy number determination performance across all datasets and robust genotype determination performance. Finally, an investigation into discordant genotypes for the synthetic dataset provides insight into misaligned reads, advancing our understanding in interpretation of short-read sequencing data in complex genomic regions. PING promises to support a new era of studies of KIR polymorphism, delivering high-resolution KIR genotypes that are highly accurate, enabling high-quality, high-throughput KIR genotyping for disease and population studies.

HLA Class I Molecules as Immune Checkpoints for NK Cell Alloreactivity and Anti-Viral Immunity in Kidney Transplantation

Natural killer (NK) cells are innate lymphocytes that can kill diseased- or virally-infected cells, mediate antibody dependent cytotoxicity and produce type I immune-associated cytokines upon activation. NK cells also contribute to the allo-immune response upon kidney transplantation either by promoting allograft rejection through lysis of cells of the transplanted organ or by promoting alloreactive T cells. In addition, they protect against viral infections upon transplantation which may be especially relevant in patients receiving high dose immune suppression. NK cell activation is tightly regulated through the integrated balance of signaling via inhibitory- and activating receptors. HLA class I molecules are critical regulators of NK cell activation through the interaction with inhibitory- as well as activating NK cell receptors, hence, HLA molecules act as critical immune checkpoints for NK cells. In the current review, we evaluate how NK cell alloreactivity and anti-viral immunity are regulated by NK cell receptors belonging to the KIR family and interacting with classical HLA class I molecules, or by NKG2A/C and LILRB1/KIR2DL4 engaging non-classical HLA-E or -G. In addition, we provide an overview of the methods to determine genetic variation in these receptors and their HLA ligands.

Following Transplantation for Acute Myelogenous Leukemia, Donor KIR Cen B02 Better Protects against Relapse than KIR Cen B01

In the treatment of acute myelogenous leukemia with allogeneic hematopoietic cell transplantation, we previously demonstrated that there is a greater protection from relapse of leukemia when the hematopoietic cell transplantation donor has either the Cen B/B KIR genotype or a genotype having two or more KIR B gene segments. In those earlier analyses, KIR genotyping could only be assessed at the low resolution of gene presence or absence. To give the analysis greater depth, we developed high-resolution KIR sequence-based typing that defines all the KIR alleles and distinguishes the expressed alleles from those that are not expressed. We now describe and analyze high-resolution KIR genotypes for 890 donors of this human transplant cohort. Cen B01 and Cen B02 are the common CenB haplotypes, with Cen B02 having evolved from Cen B01 by deletion of the KIR2DL5, 2DS3/5, 2DP1, and 2DL1 genes. We observed a consistent trend for Cen B02 to provide stronger protection against relapse than Cen B01 This correlation indicates that protection depends on the donor having inhibitory KIR2DL2 and/or activating KIR2DS2, and is enhanced by the donor lacking inhibitory KIR2DL1, 2DL3, and 3DL1. High-resolution KIR typing has allowed us to compare the strength of the interactions between the recipient's HLA class I and the KIR expressed by the donor-derived NK cells and T cells, but no clinically significant interactions were observed. The trend observed between donor Cen B02 and reduced relapse of leukemia points to the value of studying ever larger transplant cohorts.

Adaptive Admixture of HLA Class I Allotypes Enhanced Genetically Determined Strength of Natural Killer Cells in East Asians

Human natural killer (NK) cells are essential for controlling infection, cancer, and fetal development. NK cell functions are modulated by interactions between polymorphic inhibitory killer cell immunoglobulin-like receptors (KIR) and polymorphic HLA-A, -B, and -C ligands expressed on tissue cells. All HLA-C alleles encode a KIR ligand and contribute to reproduction and immunity. In contrast, only some HLA-A and -B alleles encode KIR ligands and they focus on immunity. By high-resolution analysis of KIR and HLA-A, -B, and -C genes, we show that the Chinese Southern Han (CHS) are significantly enriched for interactions between inhibitory KIR and HLA-A and -B. This enrichment has had substantial input through population admixture with neighboring populations, who contributed HLA class I haplotypes expressing the KIR ligands B*46:01 and B*58:01, which subsequently rose to high frequency by natural selection. Consequently, over 80% of Southern Han HLA haplotypes encode more than one KIR ligand. Complementing the high number of KIR ligands, the CHS KIR locus combines a high frequency of genes expressing potent inhibitory KIR, with a low frequency of those expressing activating KIR. The Southern Han centromeric KIR region encodes strong, conserved, inhibitory HLA-C-specific receptors, and the telomeric region provides a high number and diversity of inhibitory HLA-A and -B-specific receptors. In all these characteristics, the CHS represent other East Asians, whose NK cell repertoires are thus enhanced in quantity, diversity, and effector strength, likely augmenting resistance to endemic viral infections.

High-Resolution Characterization of KIR Genes in a Large North American Cohort Reveals Novel Details of Structural and Sequence Diversity

The KIR (killer-cell immunoglobulin-like receptor) region is characterized by structural variation and high sequence similarity among genes, imposing technical difficulties for analysis. We undertook the most comprehensive study to date of KIR genetic diversity in a large population sample, applying next-generation sequencing in 2,130 United States European-descendant individuals. Data were analyzed using our custom bioinformatics pipeline specifically designed to address technical obstacles in determining KIR genotypes. Precise gene copy number determination allowed us to identify a set of uncommon gene-content KIR haplotypes accounting for 5.2% of structural variation. In this cohort, KIR2DL4 is the framework gene that most varies in copy number (6.5% of all individuals). We identified phased high-resolution alleles in large multi-locus insertions and also likely founder haplotypes from which they were deleted. Additionally, we observed 250 alleles at 5-digit resolution, of which 90 have frequencies ≥1%. We found sequence patterns that were consistent with the presence of novel alleles in 398 (18.7%) individuals and contextualized multiple orphan dbSNPs within the KIR complex. We also identified a novel KIR2DL1 variant, Pro151Arg, and demonstrated by molecular dynamics that this substitution is predicted to affect interaction with HLA-C. No previous studies have fully explored the full range of structural and sequence variation of KIR as we present here. We demonstrate that pairing high-throughput sequencing with state-of-art computational tools in a large cohort permits exploration of all aspects of KIR variation including determination of population-level haplotype diversity, improving understanding of the KIR system, and providing an important reference for future studies.

High-Resolution Analysis Identifies High Frequency of KIR-A Haplotypes and Inhibitory Interactions of KIR With HLA Class I in Zhejiang Han

Killer cell immunoglobulin-like receptors (KIR) interact with human leukocyte antigen (HLA) class I molecules, modulating critical NK cell functions in the maintenance of human health. Characterizing the distribution and characteristics of KIR and HLA allotype diversity across defined human populations is thus essential for understanding the multiple associations with disease, and for directing therapies. In this study of 176 Zhejiang Han individuals from Southeastern China, we describe diversity of the highly polymorphic KIR and HLA class I genes at high resolution. KIR-A haplotypes, which carry four inhibitory receptors specific for HLA-A, B or C, are known to associate with protection from infection and some cancers. We show the Chinese Southern Han from Zhejiang are characterized by a high frequency of KIR-A haplotypes and a high frequency of C1 KIR ligands. Accordingly, interactions of inhibitory KIR2DL3 with C1+HLA are more frequent in Zhejiang Han than populations outside East Asia. Zhejiang Han exhibit greater diversity of inhibitory than activating KIR, with three-domain inhibitory KIR exhibiting the greatest degree of polymorphism. As distinguished by gene copy number and allele content, 54 centromeric and 37 telomeric haplotypes were observed. We observed 6% of the population to have KIR haplotypes containing large-scale duplications or deletions that include complete genes. A unique truncated haplotype containing only KIR2DL4 in the telomeric region was also identified. An additional feature is the high frequency of HLA-B*46:01, which may have arisen due to selection pressure from infectious disease. This study will provide further insight into the role of KIR and HLA polymorphism in disease susceptibility of Zhejiang Chinese.

The combinatorial diversity of KIR and HLA class I allotypes in Peninsular Malaysia

Killer cell immunoglobulin-like receptors (KIRs) interact with polymorphic human leucocyte antigen (HLA) class I molecules, modulating natural killer (NK) cell functions and affecting both the susceptibility and outcome of immune-mediated diseases. The KIR locus is highly diverse in gene content, copy number and allelic polymorphism within individuals and across geographical populations. To analyse currently under-represented Asian and Pacific populations, we investigated the combinatorial diversity of KIR and HLA class I in 92 unrelated Malay and 75 Malaysian Chinese individuals from the Malay Peninsula. We identified substantial allelic and structural diversity of the KIR locus in both populations and characterized novel variations at each analysis level. The Malay population is more diverse than Malay Chinese, likely representing a unique history including admixture with immigrating populations spanning several thousand years. Characterizing the Malay population are KIR haplotypes with large structural variants present in 10% individuals, and KIR and HLA alleles previously identified in Austronesian populations. Despite the differences in ancestries, the proportion of HLA allotypes that serve as KIR ligands is similar in each population. The exception is a significantly reduced frequency of interactions of KIR2DL1 with C2+ HLA-C in the Malaysian Chinese group, caused by the low frequency of C2+ HLA. One likely implication is a greater protection from preeclampsia, a pregnancy disorder associated with KIR2DL1, which shows higher incidence in the Malay than in the Malaysian Chinese. This first complete, high-resolution, characterization of combinatorial diversity of KIR and HLA in Malaysians will form a valuable reference for future clinical and population studies.

Diversity of KIR genes and their HLA-C ligands in Ugandan populations with historically varied malaria transmission intensity

BACKGROUND

Malaria is one of the most serious infectious diseases in the world. The malaria burden is greatly affected by human immunity, and immune responses vary between populations. Genetic diversity in KIR and HLA-C genes, which are important in immunity to infectious diseases, is likely to play a role in this heterogeneity. Several studies have shown that KIR and HLA-C genes influence the immune response to viral infections, but few studies have examined the role of KIR and HLA-C in malaria infection, and these have used low-resolution genotyping. The aim of this study was to determine whether genetic variation in KIR and their HLA-C ligands differ in Ugandan populations with historically varied malaria transmission intensity using more comprehensive genotyping approaches.

METHODS

High throughput multiplex quantitative real-time PCR method was used to genotype KIR genetic variants and copy number variation and a high-throughput real-time PCR method was developed to genotype HLA-C1 and C2 allotypes for 1344 participants, aged 6 months to 10 years, enrolled from Ugandan populations with historically high (Tororo District), medium (Jinja District) and low (Kanungu District) malaria transmission intensity.

RESULTS

The prevalence of KIR3DS1, KIR2DL5, KIR2DS5, and KIR2DS1 genes was significantly lower in populations from Kanungu compared to Tororo (7.6 vs 13.2%: p = 0.006, 57.2 vs 66.4%: p = 0.005, 33.2 vs 46.6%: p < 0.001, and 19.7 vs 26.7%: p = 0.014, respectively) or Jinja (7.6 vs 18.1%: p < 0.001, 57.2 vs 63.8%: p = 0.048, 33.2 vs 43.5%: p = 0.002, and 19.7 vs 30.4%: p < 0.001, respectively). The prevalence of homozygous HLA-C2 was significantly higher in populations from Kanungu (31.6%) compared to Jinja (21.4%), p = 0.043, with no significant difference between Kanungu and Tororo (26.7%), p = 0.296.

CONCLUSIONS

The KIR3DS1, KIR2DL5, KIR2DS5 and KIR2DS1 genes may partly explain differences in transmission intensity of malaria since these genes have been positively selected for in places with historically high malaria transmission intensity. The high-throughput, multiplex, real-time HLA-C genotyping PCR method developed will be useful in disease-association studies involving large cohorts.

KIR2DL4 genetic diversity in a Brazilian population sample: implications for transcription regulation and protein diversity in samples with different ancestry backgrounds

KIR2DL4 is an important immune modulator expressed in natural killer cells; HLA-G is its main ligand. We have characterized the KIR2DL4 genetic diversity by considering the promoter, all exons, and all introns in a highly admixed Brazilian population sample and by using massively parallel sequencing. We introduce a molecular method to amplify and to sequence the complete KIR2DL4 gene. To avoid the mapping bias and genotype errors commonly observed in gene families, we have developed and validated a bioinformatic pipeline designed to minimize these errors and applied it to survey the variability of 220 individuals from the State of São Paulo, southeastern Brazil. We have also compared the KIR2DL4 genetic diversity in the Brazilian cohort with the diversity previously reported by the 1000Genomes consortium. KIR2DL4 presents high linkage disequilibrium throughout the gene, with coding sequences associated with specific promoters. There are few but divergent promoter haplotypes. We have also detected many new KIR2DL4 sequences, all bearing nucleotide exchanges in introns and encoding previously described proteins. Exons 3 and 4, which encode the external domains, are the most variable. The ancestry background influences the KIR2DL4 allele frequencies and must be considered for association studies regarding KIR2DL4.

The Role of the HLA Class I α2 Helix in Determining Ligand Hierarchy for the Killer Cell Ig-like Receptor 3DL1

HLA class I molecules that represent ligands for the inhibitory killer cell Ig-like receptor (KIR) 3DL1 found on NK cells are categorically defined as those HLA-A and HLA-B allotypes containing the Bw4 motif, yet KIR3DL1 demonstrates hierarchical recognition of these HLA-Bw4 ligands. To better understand the molecular basis underpinning differential KIR3DL1 recognition, the HLA-A^Bw4 family of allotypes were investigated. Transfected human 721.221 cells expressing HLA-A*32:01 strongly inhibited primary human KIR3DL1⁺ NK cells, whereas HLA-A*24:02 and HLA-A*23:01 displayed intermediate potency and HLA-A*25:01 failed to inhibit activation of KIR3DL1⁺ NK cells. Structural studies demonstrated that recognition of HLA-A*24:02 by KIR3DL1 used identical contacts as the potent HLA-B*57:01 ligand. Namely, the D1-D2 domains of KIR3DL1 were placed over the α1 helix and α2 helix of the HLA-A*24:02 binding cleft, respectively, whereas the D0 domain contacted the side of the HLA-A*24:02 molecule. Nevertheless, functional analyses showed KIR3DL1 recognition of HLA-A*24:02 was more sensitive to substitutions within the α2 helix of HLA-A*24:02, including residues Ile¹⁴² and Lys¹⁴⁴ Furthermore, the presence of Thr¹⁴⁹ in the α2 helix of HLA-A*25:01 abrogated KIR3DL1⁺ NK inhibition. Together, these data demonstrate a role for the HLA class I α2 helix in determining the hierarchy of KIR3DL1 ligands. Thus, recognition of HLA class I is dependent on a complex interplay between the peptide repertoire, polymorphisms within and proximal to the Bw4 motif, and the α2 helix. Collectively, the data furthers our understanding of KIR3DL1 ligands and will inform genetic association and immunogenetics studies examining the role of KIR3DL1 in disease settings.

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.

MRLR: unraveling high-resolution meiotic recombination by linked reads

MOTIVATION

Meiotic recombination facilitates the transmission of exchanged genetic material between homologous chromosomes and plays a crucial role in increasing the genetic variations in eukaryotic organisms. In humans, thousands of crossover events have been identified by genotyping related family members. However, most of these crossover regions span tens to hundreds of kb, which is not sufficient resolution to accurately identify the crossover breakpoints in a typical trio family.

RESULTS

We have developed MRLR, a software using 10X linked reads to identify crossover events at a high resolution. By reconstructing the gamete genome, MRLR only requires a trio family dataset and can efficiently discover the crossover events. Using MRLR, we revealed a fine-scale pattern of crossover regions in six human families. From the two closest heterozygous alleles around the crossovers, we determined that MRLR achieved a median resolution 4.5 kb. This method can delineate a genome-wide landscape of crossover events at a precise scale, which is important for both functional and genomic features analysis of meiotic recombination.

AVAILABILITY AND IMPLEMENTATION

MRLR is freely available at https://github.com/ChongLab/MRLR, implemented in Perl.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

KIR Variation in Iranians Combines High Haplotype and Allotype Diversity With an Abundance of Functional Inhibitory Receptors

Natural killer (NK) cells are innate lymphocytes that eliminate infected and transformed cells. They discriminate healthy from diseased tissue through killer cell Ig-like receptor (KIR) recognition of HLA class I ligands. Directly impacting NK cell function, KIR polymorphism associates with infection control and multiple autoimmune and pregnancy syndromes. Here we analyze KIR diversity of 241 individuals from five groups of Iranians. These five populations represent Baloch, Kurd, and Lur, together comprising 15% of the ethnically diverse Iranian population. We identified 159 KIR alleles, including 11 not previously characterized. We also identified 170 centromeric and 94 telomeric haplotypes, and 15 different KIR haplotypes carrying either a deletion or duplication encompassing one or more complete KIR genes. As expected, comparing our data with those representing major worldwide populations revealed the greatest similarity between Iranians and Europeans. Despite this similarity we observed higher frequencies of KIR3DL1^*001 in Iran than any other population, and the highest frequency of HLA-B^*51, a Bw4-containing allotype that acts as a strong educator of KIR3DL1^*001⁺ NK cells. Compared to Europeans, the Iranians we studied also have a reduced frequency of 3DL1^*004, which encodes an allotype that is not expressed at the NK cell surface. Concurrent with the resulting high frequency of strong viable interactions between inhibitory KIR and polymorphic HLA class I, the majority of KIR-A haplotypes characterized do not express a functional activating receptor. By contrast, the most frequent KIR-B haplotype in Iran expresses only one functional inhibitory KIR and the maximum number of activating KIR. This first complete, high-resolution, characterization of the KIR locus of Iranians will form a valuable reference for future clinical and population studies.

Haplotype-Based Analysis of KIR-Gene Profiles in a South European Population-Distribution of Standard and Variant Haplotypes, and Identification of Novel Recombinant Structures

Inhibitory Killer-cell Immunoglobulin-like Receptors (KIR) specific for HLA class I molecules enable human natural killer cells to monitor altered antigen presentation in pathogen-infected and tumor cells. KIR genes display extensive copy-number variation and allelic polymorphism. They organize in a series of variable arrangements, designated KIR haplotypes, which derive from duplications of ancestral genes and sequence diversification through point mutation and unequal crossing-over events. Genomic studies have established the organization of multiple KIR haplotypes—many of them are fixed in most human populations, whereas variants of those have less certain distributions. Whilst KIR-gene diversity of many populations and ethnicities has been explored superficially (frequencies of individual genes and presence/absence profiles), less abundant are in-depth analyses of how such diversity emerges from KIR-haplotype structures. We characterize here the genetic diversity of KIR in a sample of 414 Spanish individuals. Using a parsimonious approach, we manage to explain all 38 observed KIR-gene profiles by homo- or heterozygous combinations of six fixed centromeric and telomeric motifs; of six variant gene arrangements characterized previously by us and others; and of two novel haplotypes never detected before in Caucasoids. Associated to the latter haplotypes, we also identified the novel transcribed KIR2DL5B^*0020202 allele, and a chimeric KIR2DS2/KIR2DL3 gene (designated KIR2DL3^*033) that challenges current criteria for classification and nomenclature of KIR genes and haplotypes.

Estimation of German KIR Allele Group Haplotype Frequencies

The impact of the highly polymorphic Killer-cell immunoglobulin-like receptor (KIR) gene cluster on the outcome of hematopoietic stem cell transplantation (HCST) is subject of current research. To further understand the involvement of this gene family into Natural Killer (NK) cell-mediated graft-versus-leukemia reactions, knowledge of haplotype structures, and allelic linkage is of importance. In this analysis, we estimate population-specific KIR haplotype frequencies at allele group resolution in a cohort of n = 458 German families. We addressed the polymorphism of the KIR gene complex and phasing ambiguities by a combined approach. Haplotype inference within first-degree family relations allowed us to limit the number of possible diplotypes. Structural restriction to a pattern set of 92 previously described KIR copy number haplotypes further reduced ambiguities. KIR haplotype frequency estimation was finally accomplished by means of an expectation-maximization algorithm. Applying a resolution threshold of ½ n, we were able to identify a set of 551 KIR allele group haplotypes, representing 21 KIR copy number haplotypes. The haplotype frequencies allow studying linkage disequilibrium in two-locus as well as in multi-locus analyses. Our study reveals associations between KIR haplotype structures and allele group frequencies, thereby broadening our understanding of the KIR gene complex.

KIR specificity and avidity of standard and unusual C1, C2, Bw4, Bw6 and A3/11 amino acid motifs at entire HLA:KIR interface between NK and target cells, the functional and evolutionary classification of HLA class I molecules

Natural killer (NK) cells make vital contributions to the immune system and the reproductive system. Notably, NK cells of donor origin can recognize and kill residual leukaemic cells and cure malignant patients in hematopoietic stem cell (HSC) transplant setting. NK cell function is regulated by KIRs that recognize cognate HLA class I molecules on target cells, depending on their amino acid residues. In review, we addressed the question of binding capacity and avidity of HLA class I molecules to different killer cell immunoglobulin-like receptors (KIRs) depending on all interacting amino acid residues both on HLA and KIR side. We searched PubMed database and analysed available HLA:KIR crystallographic data for amino acid residues in HLA molecules, those physically involved in binding KIRs (termed here the "entire KIR interface"). Within entire KIR interface, we selected five functional sequence motifs (14-19, 66-76, 77-84, 88-92 and 142-151) and classified them according to the conservation of their amino acid sequences among 8,942 HLA class I molecules. Although some conserved amino acid motifs were shared by different groups of KIR ligands, the HLA motif combinations were exclusive for the ligand groups. In 135 common HLA class I molecules with known HLA:KIR recognition, we found 54 combinations of five motifs in each of the KIR-binding interfaces (C1, C2, Bw4, A3/11) and conserved non-KIR-binding interfaces. Based on the entire KIR interface, this analysis allowed to classify 8,942 HLA class I molecules into KIR specificity groups. This functional and evolutionary classification of entire KIR interfaces provides a tool for unambiguously predicting HLA:KIR interactions for common and those HLA molecules that have not yet been functionally tested. Considering the entire KIR interface in HLA class I molecules, functional interactions of HLA and KIR can be predicted in immune responses, reproduction and allotransplantation. Further functional studies are needed on the HLA:KIR interaction variations caused by the repertoires of peptides presented by HLA molecules and KIR polymorphisms at allelic level.

Killer Ig-Like Receptors (KIRs): Their Role in NK Cell Modulation and Developments Leading to Their Clinical Exploitation

Natural killer (NK) cells contribute to the first line of defense against viruses and to the control of tumor growth and metastasis spread. The discovery of HLA class I specific inhibitory receptors, primarily of killer Ig-like receptors (KIRs), and of activating receptors has been fundamental to unravel NK cell function and the molecular mechanisms of tumor cell killing. Stemmed from the seminal discoveries in early '90s, in which Alessandro Moretta was the major actor, an extraordinary amount of research on KIR specificity, genetics, polymorphism, and repertoire has followed. These basic notions on NK cells and their receptors have been successfully translated to clinical applications, primarily to the haploidentical hematopoietic stem cell transplantation to cure otherwise fatal leukemia in patients with no HLA compatible donors. The finding that NK cells may express the PD-1 inhibitory checkpoint, particularly in cancer patients, may allow understanding how anti-PD-1 therapy could function also in case of HLA class Ineg tumors, usually susceptible to NK-mediated killing. This, together with the synergy of therapeutic anti-checkpoint monoclonal antibodies, including those directed against NKG2A or KIRs, emerging in recent or ongoing studies, opened new solid perspectives in cancer therapy.

Diversity of Killer Cell Immunoglobulin-Like Receptors and Disease

Natural killer (NK) cells are bone marrow-derived large granular lymphocytes defined by CD3^negCD56^pos and represent 5% to 25% of peripheral blood mononuclear cell fraction of the healthy humans. NK cells have a highly specific and sophisticated target cell recognition receptor system arbitrated by the integration of signals triggered by a multitude of inhibitory and activating receptors. Human NK cells express distinct families of receptors, including (1) killer cell immunoglobulin-like receptors, (2) killer cell lectin-like receptors, (3) leukocyte immunoglobulin-like receptors, and (4) natural cytotoxicity receptors.

Conservation, Extensive Heterozygosity, and Convergence of Signaling Potential All Indicate a Critical Role for KIR3DL3 in Higher Primates

Natural killer (NK) cell functions are modulated by polymorphic killer cell immunoglobulin-like receptors (KIR). Among 13 human KIR genes, which vary by presence and copy number, KIR3DL3 is ubiquitously present in every individual across diverse populations. No ligand or function is known for KIR3DL3, but limited knowledge of expression suggests involvement in reproduction, likely during placentation. With 157 human alleles, KIR3DL3 is also highly polymorphic and we show heterozygosity exceeds that of HLA-B in many populations. The external domains of catarrhine primate KIR3DL3 evolved as a conserved lineage distinct from other KIR. Accordingly, and in contrast to other KIR, we show the focus of natural selection does not correspond exclusively to known ligand binding sites. Instead, a strong signal for diversifying selection occurs in the D1 Ig domain at a site involved in receptor aggregation, which we show is polymorphic in humans worldwide, suggesting differential ability for receptor aggregation. Meanwhile in the cytoplasmic tail, the first of two inhibitory tyrosine motifs (ITIM) is conserved, whereas independent genomic events have mutated the second ITIM of KIR3DL3 alleles in all great apes. Together, these findings suggest that KIR3DL3 binds a conserved ligand, and a function requiring both receptor aggregation and inhibitory signal attenuation. In this model KIR3DL3 resembles other NK cell inhibitory receptors having only one ITIM, which interact with bivalent downstream signaling proteins through dimerization. Due to the extensive conservation across species, selection, and other unusual properties, we consider elucidating the ligand and function of KIR3DL3 to be a pressing question.

KIR/HLA Gene Profile Implication in Systemic Sclerosis Patients from Mexico

Introduction

Systemic Sclerosis (SSc) is an autoimmune, inflammatory, and multisystemic disease characterized by the presence of autoantibodies and fibrosis. The pathogenesis involves the interaction between immune system cells such as macrophages, NK cells, T cells, and B cells. Killer-cell Immunoglobulin-like Receptors (KIR) are expressed in NK cells and some T cell subsets that recognize HLA class I molecules as ligands and are involved in regulating the activation and inhibition of these cells. The KIR family consists of 14 genes and two pseudogenes; according to the gene content, the genotype could be AA and Bx. The aim of this study was to evaluate the association between KIR/HLA genes and genotypes with SSc and the clinical characteristics.

Methods

We included 50 SSc patients and 90 Control Subjects (CS). Genotyping of KIR, HLA-C, -Bw4, and -A^∗03/^∗11 was made by SSP-PCR.

Results

In SSc patients, a higher frequency of KIR2DL2 (p = 0.0007, p′ = 0.011), KIR2DS4del (p = 0.001, p′ = 0.021), and HLA-C2 (p = 0.02, p′ = 0.09) was found. This is the first study to evaluate the frequency of HLA-A^∗03/^∗11 in SSc patients, of which a low frequency was found in both groups. Compound genotypes KIR2DL2+/HLA-C1+ or KIR2DL2+/HLA-C2+ have a higher frequency in SSc patients. The Bx genotype was the most frequent and was associated with risk to SSc (p = 0.007, OR = 3.1, 95% CI = 1.4–7.9, p′ = 0.014). The genotypes with a higher iKIR number than aKIR (iKIR > aKIR) were found in all individuals; genotypes with 7-8 iKIR genes were increased in SSc patients. We do not find an association between the KIR genes with the clinical characteristics.

Conclusion

The results suggest that KIR2DL2 and 2DS4del could have a risk role in the development of SSc, but not with clinical manifestations.

Report from the Killer-cell Immunoglobulin-like Receptors (KIR) component of the 17th International HLA and Immunogenetics Workshop

The goals of the KIR component of the 17th International HLA and Immunogenetics Workshop (IHIW) were to encourage and educate researchers to begin analyzing KIR at allelic resolution, and to survey the nature and extent of KIR allelic diversity across human populations. To represent worldwide diversity, we analyzed 1269 individuals from ten populations, focusing on the most polymorphic KIR genes, which express receptors having three immunoglobulin (Ig)-like domains (KIR3DL1/S1, KIR3DL2 and KIR3DL3). We identified 13 novel alleles of KIR3DL1/S1, 13 of KIR3DL2 and 18 of KIR3DL3. Previously identified alleles, corresponding to 33 alleles of KIR3DL1/S1, 38 of KIR3DL2, and 43 of KIR3DL3, represented over 90% of the observed allele frequencies for these genes. In total we observed 37 KIR3DL1/S1 allotypes, 40 for KIR3DL2 and 44 for KIR3DL3. As KIR allotype diversity can affect NK cell function, this demonstrates potential for high functional diversity worldwide. Allelic variation further diversifies KIR haplotypes. We determined KIR3DL3 ∼ KIR3DL1/S1 ∼ KIR3DL2 haplotypes from five of the studied populations, and observed multiple population-specific haplotypes in each. This included 234 distinct haplotypes in European Americans, 191 in Ugandans, 35 in Papuans, 95 in Egyptians and 86 in Spanish populations. For another 35 populations, encompassing 642,105 individuals we focused on KIR3DL2 and identified another 375 novel alleles, with approximately half of them observed in more than one individual. The KIR allelic level data gathered from this project represents the most comprehensive summary of global KIR allelic diversity to date, and continued analysis will improve understanding of KIR allelic polymorphism in global populations. Further, the wealth of new data gathered in the course of this workshop component highlights the value of collaborative, community-based efforts in immunogenetics research, exemplified by the IHIW.

Major histocompatibility complex (MHC) fragment numbers alone - in Atlantic cod and in general - do not represent functional variability

This correspondence concerns a publication by Malmstrøm et al. in Nature Genetics in October 2016. Malmstrøm et al. made an important contribution to fish phylogeny research by using low-coverage genome sequencing for comparison of 66 teleost (modern bony) fish species, with 64 of those 66 belonging to the species-rich clade Neoteleostei, and with 27 of those 64 belonging to the order Gadiformes. For these 66 species, Malmstrøm et al. estimated numbers of genes belonging to the major histocompatibility complex (MHC) class I lineages U and Z and concluded that in teleost fish these combined numbers are positively associated with, and a driving factor of, the rates of establishment of new fish species (speciation rates). They also claimed that functional genes for the MHC class II system molecules MHC IIA, MHC IIB, CD4 and CD74 were lost in early Gadiformes. Our main criticisms are (1) that the authors did not provide sufficient evidence for presence or absence of intact functional MHC class I or MHC class II system genes, (2) that they did not discuss that an MHC subpopulation gene number alone is a very incomplete measure of MHC variance, and (3) that the MHC system is more likely to reduce speciation rates than to enhance them. Furthermore, their use of the Ornstein-Uhlenbeck model is a typical example of overly naïve use of that model system. In short, we conclude that their new model of MHC class I evolution, reflected in their title "Evolution of the immune system influences speciation rates in teleost fish", is unsubstantiated, and that their "pinpointing" of the functional loss of the MHC class II system and all the important MHC class II system genes to the onset of Gadiformes is preliminary, because they did not sufficiently investigate the species at the clade border.

Analysis of KIR3DP1 Polymorphism Provides Relevant Information on Centromeric KIR Gene Content

Four killer cell Ig-like receptor (KIR) genes, collectively referred to as framework genes, characterize almost all KIR haplotypes. In particular, KIR3DL3 and KIR3DL2 mark the ends of the locus, whereas KIR3DP1 and KIR2DL4 are located in the central part. A recombination hot spot, mapped between KIR3DP1 and KIR2DL4, splits the haplotypes into two regions: a centromeric (Cen) region (spanning from KIR3DL3 to KIR3DP1) and a telomeric region (from KIR2DL4 to KIR3DL2), both varying in KIR gene content. In this study, we analyzed KIR3DP1 polymorphism in a cohort of 316 healthy, unrelated individuals. To this aim, we divided KIR3DP1 alleles into two groups by the use of a sequence-specific primer- PCR approach. Our data clearly indicated that KIR3DP1 alleles present on haplotypes carrying Cen-A or Cen-B1 regions differ from those having Cen-B2 motifs. Few donors (∼3%) made exceptions, and they were all, except one, characterized by uncommon haplotypes, including either KIR deletions or KIR duplications. Consequently, as KIR2DL1 is present in Cen-A and Cen-B1 regions but absent in Cen-B2 regions, we demonstrated that KIR3DP1 polymorphism might represent a suitable marker for KIR2DL1 gene copy number analysis. Moreover, because Cen-B1 and Cen-B2 regions are characterized by different KIR3DP1 alleles, we showed that KIR3DP1 polymorphism analysis also provides information to dissect between Cen-B1/Cen-B1 and Cen-B1/Cen-B2 donors. Taken together, our data suggest that the analysis of KIR3DP1 polymorphism should be included in KIR repertoire evaluation.

Major histocompatibility complex (MHC) fragment numbers alone - in Atlantic cod and in general - do not represent functional variability

This correspondence concerns a publication by Malmstrøm et al. in Nature Genetics in October 2016. Malmstrøm et al. made an important contribution to fish phylogeny research by using low-coverage genome sequencing for comparison of 66 teleost (modern bony) fish species, with 64 of those 66 belonging to the species-rich clade Neoteleostei, and with 27 of those 64 belonging to the order Gadiformes. For these 66 species, Malmstrøm et al. estimated numbers of genes belonging to the major histocompatibility complex (MHC) class I lineages U and Z and concluded that in teleost fish these combined numbers are positively associated with, and a driving factor of, the rates of establishment of new fish species (speciation rates). They also claimed that functional genes for the MHC class II system molecules MHC IIA, MHC IIB, CD4 and CD74 were lost in early Gadiformes. Our main criticisms are (1) that the authors did not provide sufficient evidence for presence or absence of intact functional MHC class I or MHC class II system genes, (2) that they did not discuss that an MHC subpopulation gene number alone is a very incomplete measure of MHC variance, and (3) that the MHC system is more likely to reduce speciation rates than to enhance them. Furthermore, their use of the Ornstein-Uhlenbeck model is a typical example of overly naïve use of that model system. In short, we conclude that their new model of MHC class I evolution, reflected in their title “Evolution of the immune system influences speciation rates in teleost fish”, is unsubstantiated, and that their “pinpointing” of the functional loss of the MHC class II system and all the important MHC class II system genes to the onset of Gadiformes is preliminary, because they did not sufficiently investigate the species at the clade border.

Different Selected Mechanisms Attenuated the Inhibitory Interaction of KIR2DL1 with C2+ HLA-C in Two Indigenous Human Populations in Southern Africa

The functions of human NK cells in defense against pathogens and placental development during reproduction are modulated by interactions of killer cell Ig-like receptors (KIRs) with HLA-A, -B and -C class I ligands. Both receptors and ligands are highly polymorphic and exhibit extensive differences between human populations. Indigenous to southern Africa are the KhoeSan, the most ancient group of modern human populations, who have highest genomic diversity worldwide. We studied two KhoeSan populations, the Nama pastoralists and the ≠Khomani San hunter-gatherers. Comprehensive next-generation sequence analysis of HLA-A, -B, and -C and all KIR genes identified 248 different KIR and 137 HLA class I, which assort into ∼200 haplotypes for each gene family. All 74 Nama and 78 ≠Khomani San studied have different genotypes. Numerous novel KIR alleles were identified, including three arising by intergenic recombination. On average, KhoeSan individuals have seven to eight pairs of interacting KIR and HLA class I ligands, the highest diversity and divergence of polymorphic NK cell receptors and ligands observed to date. In this context of high genetic diversity, both the Nama and the ≠Khomani San have an unusually conserved, centromeric KIR haplotype that has arisen to high frequency and is different in the two KhoeSan populations. Distinguishing these haplotypes are independent mutations in KIR2DL1, which both prevent KIR2DL1 from functioning as an inhibitory receptor for C2⁺ HLA-C. The relatively high frequency of C2⁺ HLA-C in the Nama and the ≠Khomani San appears to have led to natural selection against strong inhibitory C2-specific KIR.

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.

Association between KIR genes and dust mite sensitization in a Brazilian population

BACKGROUND

Killer cell immunoglobulin-like receptors (KIRs), found on the surface of natural killer (NK) cells, play a key role in controlling the innate response. Such response depends on a series of cellular interactions between these receptors and HLA activating/inhibiting ligands. Atopic diseases have been associated with genes that regulate cytokine production and HLA genes, which may either protect or predispose to hypersensitivity.

OBJECTIVE

To verify an association study of KIR genes with sensitization to the following mites: Dermatophagoides farinae, Dermatophagoides pteronyssinus, and Blomia tropicalis.

METHODS

A total of 341 children aged up to 14 years, were classified as mite-sensitive or mite-insensitive after undergoing a skin prick test for immediate allergic reactions. The presence/absence of KIR genes and their human leukocyte antigen (HLA) ligands was determined by polymerase chain reaction-sequence specific oligonucleotide (PCR-SSO) with the commercial kit LabType™ using Luminex™.

RESULTS

The frequencies of KIR genes and their respective class I HLA ligands and the frequency of haplotypes were performed in sensitive and insensitive individuals, and no significant differences were found.

CONCLUSION

Our results suggest no influence of KIR genes on resistance/susceptibility to sensitization to dust mites.

Novel KIR genotypes and gene copy number variations in northeastern Thais

KIR (Killer Immunoglobulin-like Receptor) variants influence immune responses and are genetic factors in disease susceptibility. Using sequence-specific priming PCR, we have previously described the diversity of KIR genes in term of presence/absence in northeastern Thais (NETs). To provide additional resolution beyond conventional methods, quantitative PCR was applied to determine KIR copy number profiles. Novel expanded and contracted KIR copy number profiles were identified at cumulatively high frequencies. These all comprise haplotypes with duplication (6·9%) or deletion (2·7%) of KIR3DL1/S1 along with adjacent genes. Five expanded KIR profiles comprised haplotypes with duplications of KIR2DP1, 2DL1, 3DP1, 2DL4, 3DL1/S1 and 2DS1/4, whereas two contracted profiles contained only a single copy of KIR3DP1, 3DL1/S1 and 2DL4. Using a KIR haplotype prediction program (KIR Haplotype Identifier), 14% of NET haplotypes carried atypical haplotypes based on the gene copy number data.

Distribution of KIR genes and KIR2DS4 gene variants in two Mexican Mestizo populations

Killer immunoglobulin-like receptors (KIR) are transmembrane proteins that regulate NK and T cell subsets by recognizing HLA-I molecules as ligands. The KIR gene family consists of 16 genes, located at chromosome 19q13.4. KIR gene frequencies vary among populations. In Mexico, HLA and genetic ancestry studies show that Mestizo populations have different genetic backgrounds based on admixture with European, African, and Asian ancestry. This study aimed to evaluate the frequencies of KIR genes and genotypes in Guerrero and Jalisco, two Mexican Mestizo populations located in the south and the west of the country, respectively, and to compare these frequencies with those of other populations. KIR genotyping was performed by SSP-PCR. We observed that KIR gene frequencies were similar in both populations. There were 24 genotypes observed in Guerrero, 38 genotypes observed in Jalisco, 15 genotypes shared in both populations and 32 genotypes unique to one population or the other. In 10 individuals, nine novel genotypes were identified. KIR2DS4 gene variants showed significant differences: The KIR2DS4full gene was more common in Guerrero (p<0.0001), and the KIR2DS4del variant was more common in Jalisco (p<0.05). Differences in KIR2DS4 gene variants and genotypic profiles could be influenced by the genetic admixture in both regions.

Hematopoietic stem cell transplantation: Improving alloreactive Bw4 donor selection by genotyping codon 86 of KIR3DL1/S1

KIR3DL1 is a natural killer (NK) cell receptor that recognizes the Bw4 epitope of human leukocyte antigen (HLA) class I molecules. Following hematopoietic stem cell transplantation for patients lacking Bw4, KIR3DL1-expressing NK cells from Bw4-positive donors can be alloreactive and eliminate tumor cells. However, KIR3DL1 alleles having T instead of C at nucleotide 320 (encoding leucine 86 instead of serine 86) are not expressed on the cell surface. Thus, not all individuals testing positive for KIR3DL1 are optimal donors for Bw4-negative recipients. Therefore, we developed a method for genotyping codon 86, which was validated by its perfect correlation with NK cell phenotype for 100 donors of diverse KIR3DL1/S1 genotype. We typed 600 donors and found that ∼12.2% had the KIR3DL1 gene, but did not express cell-surface KIR3DL1. By contrast, high-expressing allotypes were identified when haplotypes from four families with duplicated KIR3DL1/S1 genes were characterized at high resolution. Identifying donors who have KIR3DL1 but lack cell-surface KIR3DL1 would refine donor selection. With this technique, the number of individuals identified who may not be optimal donors for Bw4-negative patients increases by threefold, when compared with standard methods. Taken together, we propose that allele typing of killer cell Ig-like receptor (KIR) polymorphisms should become a standard practice when selecting donors.

HLA-Bw4-I-80 Isoform Differentially Influences Clinical Outcome As Compared to HLA-Bw4-T-80 and HLA-A-Bw4 Isoforms in Rituximab or Dinutuximab-Based Cancer Immunotherapy

Killer-cell immunoglobulin-like receptors (KIRs) are a family of glycoproteins expressed primarily on natural killer cells that can regulate their function. Inhibitory KIRs recognize MHC class I molecules (KIR-ligands) as ligands. We have reported associations of KIRs and KIR-ligands for patients in two monoclonal antibody (mAb)-based trials: (1) A Children’s Oncology Group (COG) trial for children with high-risk neuroblastoma randomized to immunotherapy treatment with dinutuximab (anti-GD2 mAb) + GM-CSF + IL-2 + isotretinion or to treatment with isotretinoin alone and (2) An Eastern Cooperative Oncology Group (ECOG) trial for adults with low-tumor burden follicular lymphoma responding to an induction course of rituximab (anti-CD20 mAb) and randomized to treatment with maintenance rituximab or no-maintenance rituximab. In each trial, certain KIR/KIR-ligand genotypes were associated with clinical benefit for patients randomized to immunotherapy treatment (immunotherapy in COG; maintenance rituximab in ECOG) as compared to patients that did not receive the immunotherapy [isotretinoin alone (COG); no-maintenance (ECOG)]. Namely, patients with both KIR3DL1 and its HLA-Bw4 ligand (KIR3DL1+/HLA-Bw4+ genotype) had improved clinical outcomes if randomized to immunotherapy regimens, as compared to patients with the KIR3DL1+/HLA-Bw4+ genotype randomized to the non-immunotherapy regimen. Conversely, patients that did not have the KIR3DL1+/HLA-Bw4+ genotype showed no evidence of a difference in outcome if receiving the immunotherapy vs. no-immunotherapy. For each trial, HLA-Bw4 status was determined by assessing the genotypes of three separate isoforms of HLA-Bw4: (1) HLA-B-Bw4 with threonine at amino acid 80 (B-Bw4-T80); (2) HLA-B-Bw4 with isoleucine at amino acid 80 (HLA-B-Bw4-I80); and (3) HLA-A with a Bw4 epitope (HLA-A-Bw4). Here, we report on associations with clinical outcome for patients with KIR3DL1 and these separate isoforms of HLA-Bw4. Patients randomized to immunotherapy with KIR3DL1+/A-Bw4+ or with KIR3DL1+/B-Bw4-T80+ had better outcome vs. those randomized to no-immunotherapy, whereas for those with KIR3DL1+/B-Bw4-I80+ there was no evidence of a difference based on immunotherapy vs. no-immunotherapy. Additionally, we observed differences within treatment types (either within immunotherapy or no-immunotherapy) that were associated with the genotype status for the different KIR3DL1/HLA-Bw4-isoforms. These studies suggest that specific HLA-Bw4 isoforms may differentially influence response to these mAb-based immunotherapy, further confirming the involvement of KIR-bearing cells in tumor-reactive mAb-based cancer immunotherapy.

NK Cells: Uncertain Allies against Malaria

Until recently, studies of natural killer (NK) cells in infection have focused almost entirely on their role in viral infections. However, there is an increasing awareness of the potential for NK cells to contribute to the control of a wider range of pathogens, including intracellular parasites such as Plasmodium spp. Given the high prevalence of parasitic diseases in the developing world and the devastating effects these pathogens have on large numbers of vulnerable people, investigating interactions between NK cells and parasitized host cells presents the opportunity to reveal novel immunological mechanisms with the potential to aid efforts to eradicate these diseases. The capacity of NK cells to produce inflammatory cytokines early after malaria infection, as well as a possible role in direct cytotoxic killing of malaria-infected cells, suggests a beneficial impact of NK cells in this disease. However, NK cells may also contribute to overproduction of pro-inflammatory cytokines and the consequent immunopathology. As comparatively little is known about the role of NK cells later in the course of infection, and growing evidence suggests that heterogeneity in NK cell responses to malaria may be influenced by KIR/HLA interactions, a better understanding of the mechanisms by which NK cells might directly interact with parasitized cells may reveal a new role for these cells in the course of malaria infection.

Deciphering the killer-cell immunoglobulin-like receptor system at super-resolution for natural killer and T-cell biology

Killer-cell immunoglobulin-like receptors (KIRs) are components of two fundamental biological systems essential for human health and survival. First, they contribute to host immune responses, both innate and adaptive, through their expression by natural killer cells and T cells. Second, KIR play a key role in regulating placentation, and hence reproductive success. Analogous to the diversity of their human leucocyte antigen class I ligands, KIR are extremely polymorphic. In this review, we describe recent developments, fuelled by methodological advances, that are helping to decipher the KIR system in terms of haplotypes, polymorphisms, expression patterns and their ligand interactions. These developments are delivering deeper insight into the relevance of KIR in immune system function, evolution and disease.

qKAT: a high-throughput qPCR method for KIR gene copy number and haplotype determination

Killer cell immunoglobulin-like receptors (KIRs), expressed on natural killer cells and T cells, have considerable biomedical relevance playing significant roles in immunity, pregnancy and transplantation. The KIR locus is one of the most complex and polymorphic regions of the human genome. Extensive sequence homology and copy number variation makes KIRs technically laborious and expensive to type. To aid the investigation of KIRs in human disease we developed a high-throughput, multiplex real-time polymerase chain reaction method to determine gene copy number for each KIR locus. We used reference DNA samples to validate the accuracy and a cohort of 1698 individuals to evaluate capability for precise copy number discrimination. The method provides improved information and identifies KIR haplotype alterations that were not previously visible using other approaches.

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.

Killer cell immunoglobulin-like receptor 3DL1 polymorphism defines distinct hierarchies of HLA class I recognition

Rossjohn, Brooks, Vivian, and colleagues provide the most complete picture to date of the impact of KIR3DL1 polymorphism on HLA class I recognition, which can be used to both reevaluate previous work on the involvement of KIR3DL1 in disease as well as inform future disease association studies.

Natural killer (NK) cells play a key role in immunity, but how HLA class I (HLA-I) and killer cell immunoglobulin-like receptor 3DL1 (KIR3DL1) polymorphism impacts disease outcome remains unclear. KIR3DL1 (*001/*005/*015) tetramers were screened for reactivity against a panel of HLA-I molecules. This revealed different and distinct hierarchies of specificity for each KIR3DL1 allotype, with KIR3DL1*005 recognizing the widest array of HLA-I ligands. These differences were further reflected in functional studies using NK clones expressing these specific KIR3DL1 allotypes. Unexpectedly, the Ile/Thr80 dimorphism in the Bw4-motif did not categorically define strong/weak KIR3DL1 recognition. Although the KIR3DL1*001, *005, and *015 polymorphisms are remote from the KIR3DL1–HLA-I interface, the structures of these three KIR3DL1–HLA-I complexes showed that the broader HLA-I specificity of KIR3DL1*005 correlated with an altered KIR3DL1*005 interdomain positioning and increased mobility within its ligand-binding site. Collectively, we provide a generic framework for understanding the impact of KIR3DL1 polymorphism on the recognition of HLA-I allomorphs.

KIR and HLA haplotype analysis in a family lacking the KIR 2DL1-2DP1 genes

The killer cell immunoglobulin-like receptor (KIR) gene cluster exhibits extensive allelic and haplotypic diversity that is observed as presence/absence of genes, resulting in expansion and contraction of KIR haplotypes and by allelic variation of individual KIR genes. We report a case of KIR pseudogene 2DP1 and 2DL1 gene absence in members of one family with the children suffering from acute myelogenous leukemia (AML). Killer cell immunoglobulin-like receptor low resolution genotyping was performed by the polymerase chain reaction (PCR)-sequence-specific primers (SSP)/sequence-specific oligonucleotide (SSO) method and haplotype assignment was done by gene content analysis. Both parents and the maternal grandfather, shared the same Cen-B2 KIR haplotype, containing KIR 3DL3, -2DS2, -2DL2 and -3DP1 genes. The second haplotype in the KIR genotype of the mother and grandfather was Tel-A1 with KIR 2DL4 (normal and deleted variant), -3DL1, -22 bp deletion variant of the 2DS4 gene and -3DL2, while the second haplotype in the KIR genotype of the father was Tel-B1 with 2DL4 (normal variant), -3DS1, -2DL5, -2DS5, -2DS1 and 3DL2 genes. Haplotype analysis in all three offsprings revealed that the children inherited the Cen-B2 haplotype with the same gene content but two of the children inherited a deleted variant of the 2DL4 gene, while the third child inherited a normal one. The second haplotype of all three offspring contained KIR 2DL4, -2DL5, -2DS1, -2DS4 (del 22bp variant), -2DS5, -3DL1 and -3DL2 genes, which was the basis of the assumption that there is a hybrid haplotype and that the present 3DL1 gene is a variant of the 3DS1 gene. Due to consanguinity among the ancestors, the results of KIR segregation analysis showed the existence of a very rare KIR genotype in the offspring. The family who is the subject of this case is even more interesting because the father was 10/10 human leukocyte antigen (HLA)-matched to his daughter, all members of the family have the “best” donor KIR-B content and the presence of a rare KIR genotype with KIR 2DP1-2DL1 genes absence.