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

High KIR diversity in Uganda and Botswana children living with HIV

Killer-cell immunoglobulin-like receptors (KIRs) are essential components of the innate immune system found on the surfaces of natural killer (NK) cells. The KIRs encoding genes are located on chromosome 19q13.4 and are genetically diverse across populations. KIRs are associated with various disease states including HIV progression, and are linked to transplantation rejection and reproductive success. However, there is limited knowledge on the diversity of KIRs from Uganda and Botswana HIV-infected paediatric cohorts, with high endemic HIV rates. We used next-generation sequencing technologies on 312 (246 Uganda, 66 Botswana) samples to generate KIR allele data and employed customised bioinformatics techniques for allelic, allotype and disease association analysis. We show that these sample sets from Botswana and Uganda have different KIRs of different diversities. In Uganda, we observed 147 vs 111 alleles in the Botswana cohort, which had a more than 1 % frequency. We also found significant deviation towards homozygosity for the KIR3DL2 gene for both rapid (RPs) and long-term non-progressors (LTNPs)in the Ugandan cohort. The frequency of the bw4-80I ligand was also significantly higher among the LTNPs than RPs (8.9 % Vs 2.0%, P-value: 0.032). In the Ugandan cohort, KIR2DS4*001 (OR: 0.671, 95 % CI: 0.481-0.937, FDR adjusted Pc=0.142) and KIR2DS4*006 (OR: 2.519, 95 % CI: 1.085-5.851, FDR adjusted Pc=0.142) were not associated with HIV disease progression after adjustment for multiple testing. Our study results provide additional knowledge of the genetic diversity of KIRs in African populations and provide evidence that will inform future immunogenetics studies concerning human disease susceptibility, evolution and host immune responses.

An archaic HLA class I receptor allele diversifies natural killer cell-driven immunity in First Nations peoples of Oceania

Genetic variation in host immunity impacts the disproportionate burden of infectious diseases that can be experienced by First Nations peoples. Polymorphic human leukocyte antigen (HLA) class I and killer cell immunoglobulin-like receptors (KIRs) are key regulators of natural killer (NK) cells, which mediate early infection control. How this variation impacts their responses across populations is unclear. We show that HLA-A∗24:02 became the dominant ligand for inhibitory KIR3DL1 in First Nations peoples across Oceania, through positive natural selection. We identify KIR3DL1∗114, widespread across and unique to Oceania, as an allele lineage derived from archaic humans. KIR3DL1∗114+NK cells from First Nations Australian donors are inhibited through binding HLA-A∗24:02. The KIR3DL1∗114 lineage is defined by phenylalanine at residue 166. Structural and binding studies show phenylalanine 166 forms multiple unique contacts with HLA-peptide complexes, increasing both affinity and specificity. Accordingly, assessing immunogenetic variation and the functional implications for immunity are fundamental toward understanding population-based disease associations.

Full-resolution HLA and KIR gene annotations for human genome assemblies

The human leukocyte antigen (HLA) genes and the killer cell immunoglobulin-like receptor (KIR) genes are critical to immune responses and are associated with many immune-related diseases. Located in highly polymorphic regions, it is difficult to study them with traditional short-read alignment-based methods. Although modern long-read assemblers can often assemble these genes, using existing tools to annotate HLA and KIR genes in these assemblies remains a nontrivial task. Here, we describe Immuannot, a new computation tool to annotate the gene structures of HLA and KIR genes and to type the allele of each gene. Applying Immuannot to 56 regional and 212 whole-genome assemblies from previous studies, we annotate 9931 HLA and KIR genes and found that almost half of these genes, 4068, have novel sequences compared with the current Immuno Polymorphism Database (IPD). These novel gene sequences are represented by 2664 distinct alleles, some of which contained nonsynonymous variations, resulting in 92 novel protein sequences. We demonstrate the complex haplotype structures at the two loci and report the linkage between HLA/KIR haplotypes and gene alleles. We anticipate that Immuannot will speed up the discovery of new HLA/KIR alleles and enable the association of HLA/KIR haplotype structures with clinical outcomes in the future.

Genetic complexity of killer-cell immunoglobulin-like receptor genes in human pangenome assemblies

The killer-cell immunoglobulin-like receptor (KIR) gene complex, a highly polymorphic region of the human genome that encodes proteins involved in immune responses, poses strong challenges in genotyping owing to its remarkable genetic diversity and structural intricacy. Accurate analysis of KIR alleles, including their structural variations, is crucial for understanding their roles in various immune responses. Leveraging the high-quality genome assemblies from the Human Pangenome Reference Consortium (HPRC), we present a novel bioinformatic tool, the structural KIR annoTator (SKIRT), to investigate gene diversity and facilitate precise KIR allele analysis. In 47 HPRC-phased assemblies, SKIRT identifies a recurrent novel KIR2DS4/3DL1 fusion gene in the paternal haplotype of HG02630 and maternal haplotype of NA19240. Additionally, SKIRT accurately identifies eight structural variants and 15 novel nonsynonymous alleles, all of which are independently validated using short-read data or quantitative polymerase chain reaction. Our study has discovered a total of 570 novel alleles, among which eight haplotypes harbor at least one KIR gene duplication, six haplotypes have lost at least one framework gene, and 75 out of 94 haplotypes (79.8%) carry at least five novel alleles, thus confirming KIR genetic diversity. These findings are pivotal in providing insights into KIR gene diversity and serve as a solid foundation for understanding the functional consequences of KIR structural variations. High-resolution genome assemblies offer unprecedented opportunities to explore polymorphic regions that are challenging to investigate using short-read sequencing methods. The SKIRT pipeline emerges as a highly efficient tool, enabling the comprehensive detection of the complete spectrum of KIR alleles within human genome assemblies.

HLA and KIR genetic association and NK cells in anti-NMDAR encephalitis

Introduction

Genetic predisposition to autoimmune encephalitis with antibodies against N-methyl-D-aspartate receptor (NMDAR) is poorly understood. Given the diversity of associated environmental factors (tumors, infections), we hypothesized that human leukocyte antigen (HLA) and killer-cell immunoglobulin-like receptors (KIR), two extremely polymorphic gene complexes key to the immune system, might be relevant for the genetic predisposition to anti-NMDAR encephalitis. Notably, KIR are chiefly expressed by Natural Killer (NK) cells, recognize distinct HLA class I allotypes and play a major role in anti-tumor and anti-infection responses.

Methods

We conducted a Genome Wide Association Study (GWAS) with subsequent control-matching using Principal Component Analysis (PCA) and HLA imputation, in a multi-ethnic cohort of anti-NMDAR encephalitis (n=479); KIR and HLA were further sequenced in a large subsample (n=323). PCA-controlled logistic regression was then conducted for carrier frequencies (HLA and KIR) and copy number variation (KIR). HLA-KIR interaction associations were also modeled. Additionally, single cell sequencing was conducted in peripheral blood mononuclear cells from 16 cases and 16 controls, NK cells were sorted and phenotyped.

Results

Anti-NMDAR encephalitis showed a weak HLA association with DRB1*01:01~DQA1*01:01~DQB1*05:01 (OR=1.57, 1.51, 1.45; respectively), and DRB1*11:01 (OR=1.60); these effects were stronger in European descendants and in patients without an underlying ovarian teratoma. More interestingly, we found increased copy number variation of KIR2DL5B (OR=1.72), principally due to an overrepresentation of KIR2DL5B*00201. Further, we identified two allele associations in framework genes, KIR2DL4*00103 (25.4% vs. 12.5% in controls, OR=1.98) and KIR3DL3*00302 (5.3% vs. 1.3%, OR=4.44). Notably, the ligands of these KIR2DL4 and KIR3DL3, respectively, HLA-G and HHLA2, are known to act as immune checkpoint with immunosuppressive functions. However, we did not find differences in specific KIR-HLA ligand interactions or HLA-G polymorphisms between cases and controls. Similarly, gene expression of CD56dim or CD56bright NK cells did not differ between cases and controls.

Discussion

Our observations for the first time suggest that the HLA-KIR axis might be involved in anti-NMDAR encephalitis. While the genetic risk conferred by the identified polymorphisms appears small, a role of this axis in the pathophysiology of this disease appears highly plausible and should be analyzed in future studies.

IGSF8 is an innate immune checkpoint and cancer immunotherapy target

Antigen presentation defects in tumors are prevalent mechanisms of adaptive immune evasion and resistance to cancer immunotherapy, whereas how tumors evade innate immunity is less clear. Using CRISPR screens, we discovered that IGSF8 expressed on tumors suppresses NK cell function by interacting with human KIR3DL2 and mouse Klra9 receptors on NK cells. IGSF8 is normally expressed in neuronal tissues and is not required for cell survival in vitro or in vivo. It is overexpressed and associated with low antigen presentation, low immune infiltration, and worse clinical outcomes in many tumors. An antibody that blocks IGSF8-NK receptor interaction enhances NK cell killing of malignant cells in vitro and upregulates antigen presentation, NK cell-mediated cytotoxicity, and T cell signaling in vivo. In syngeneic tumor models, anti-IGSF8 alone, or in combination with anti-PD1, inhibits tumor growth. Our results indicate that IGSF8 is an innate immune checkpoint that could be exploited as a therapeutic target.

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.

Full resolution HLA and KIR genes annotation for human genome assemblies

The HLA (Human Leukocyte Antigen) genes and the KIR (Killer cell Immunoglobulin-like Receptor) genes are critical to immune responses and are associated with many immune-related diseases. Located in highly polymorphic regions, they are hard to be studied with traditional short-read alignment-based methods. Although modern long-read assemblers can often assemble these genes, using existing tools to annotate HLA and KIR genes in these assemblies remains a non-trivial task. Here, we describe Immuannot, a new computation tool to annotate the gene structures of HLA and KIR genes and to type the allele of each gene. Applying Immuannot to 56 regional and 212 whole-genome assemblies from previous studies, we annotated 9,931 HLA and KIR genes and found that almost half of these genes, 4,068, had novel sequences compared to the current Immuno Polymorphism Database (IPD). These novel gene sequences were represented by 2,664 distinct alleles, some of which contained non-synonymous variations resulting in 92 novel protein sequences. We demonstrated the complex haplotype structures at the two loci and reported the linkage between HLA/KIR haplotypes and gene alleles. We anticipate that Immuannot will speed up the discovery of new HLA/KIR alleles and enable the association of HLA/KIR haplotype structures with clinical outcomes in the future.

High-resolution characterization of KIR genes polymorphism in healthy subjects from the Bulgarian population-A pilot study

Although killer-cell immunoglobulin-like receptor (KIR) gene content has been widely studied in health and disease, with the advancement of next-generation sequencing (NGS) technology the high-resolution characterization of this complex gene region has become achievable. KIR allele-level diversity has lately been described across human populations. The present study aimed to analyze for the first time the allele-level polymorphism of nine KIR genes in 155 healthy, unrelated individuals from the Bulgarian population by applying NGS. The highest degree of polymorphism was detected for the KIR3DL3 gene with 40 observed alleles at five-digit resolution in total, 22 of which were common. On the other hand, the KIR3DS1 gene was found to have the lowest degree of polymorphism among the studied KIR genes with one common allele: KIR3DS1*01301 (31.6%). To better understand KIR allelic associations and patterns in Bulgarians, we have estimated the pairwise linkage disequilibrium (LD) for the 10 KIR loci, where KIR2DL3*00501 allele was found in strong LD with KIR2DL1*00101 (D' = 1.00, R2  = 0.742). This is the first study investigating KIR polymorphism at the allele level in a population from the South-East European region. Considering the effect of the populationally shaped KIR allelic polymorphism on NK cell function, this data could lead to a better understanding of the genetic heterogeneity of this region and can be carried into clinical practice by improvement of the strategies taken for NK-mediated diseases.

HLA-DPB1*13:01 associates with enhanced, and KIR2DS4*001 with diminished protection from developing severe COVID-19

Extreme polymorphism of HLA and killer-cell immunoglobulin-like receptors (KIR) differentiates immune responses across individuals. Additional to T cell receptor interactions, subsets of HLA class I act as ligands for inhibitory and activating KIR, allowing natural killer (NK) cells to detect and kill infected cells. We investigated the impact of HLA and KIR polymorphism on the severity of COVID-19. High resolution HLA class I and II and KIR genotypes were determined from 403 non-hospitalized and 1575 hospitalized SARS-CoV-2 infected patients from Italy collected in 2020. We observed that possession of the activating KIR2DS4*001 allotype is associated with severe disease, requiring hospitalization (OR = 1.48, 95% CI 1.20-1.85, pc  = 0.017), and this effect is greater in individuals homozygous for KIR2DS4*001 (OR = 3.74, 95% CI 1.75-9.29, pc  = 0.003). We also observed the HLA class II allotype, HLA-DPB1*13:01 protects SARS-CoV-2 infected patients from severe disease (OR = 0.49, 95% CI 0.33-0.74, pc  = 0.019). These association analyses were replicated using logistic regression with sex and age as covariates. Autoantibodies against IFN-α associated with COVID-19 severity were detected in 26% of 156 hospitalized patients tested. HLA-C*08:02 was more frequent in patients with IFN-α autoantibodies than those without, and KIR3DL1*01502 was only present in patients lacking IFN-α antibodies. These findings suggest that KIR and HLA polymorphism is integral in determining the clinical outcome following SARS-CoV-2 infection, by influencing the course both of innate and adaptive immunity.

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.

Breaking down barriers: The potential of smarter CAR-engineered NK cells against solid tumors

Natural killer (NK) cells are considered to be the foremost fighters of our innate immune system against foreign invaders and thus tend to promptly latch onto the virus-infected and tumor/cancerous cells, killing them through phagocytosis. At present, the application of genetically engineered Chimeric antigen receptor (CAR) receptors ensures a guaranteed optimistic response with NK cells and would not allow the affected cells to dodge or escape unchecked. Hence the specificity and uniqueness of CAR-NK cells over CAR-T therapy make them a better immunotherapeutic choice to reduce the load of trafficking of numerous tumor cells near the healthy cell populations in a more intact way than offered by CAR-T immunotherapy. Our review mainly focuses on the preclinical, clinical, and recent advances in clinical research trials and further strategies to achieve an augmented and efficient cure against solid tumors.

Non-Expressed Donor KIR3DL1 Alleles May Represent a Risk Factor for Relapse after T-Replete Haploidentical Hematopoietic Stem Cell Transplantation

KIR3DL1 alleles are expressed at different levels on the natural killer (NK) cell surface. In particular, the non-expressed KIR3DL1*004 allele appears to be common in Caucasian populations. However, the overall distribution of non-expressed KIR3DL1 alleles and their clinical relevance after T-replete haploidentical hematopoietic stem cell transplantation (hHSCT) with post-transplant cyclophosphamide remain poorly documented in European populations. In a cohort of French blood donors (N = 278), we compared the distribution of expressed and non-expressed KIR3DL1 alleles using next-generation sequencing (NGS) technology combined with multi-color flow cytometry. We confirmed the predominance of the non-expressed KIR3DL1*004 allele. Using allele-specific constructs, the phenotype and function of the uncommon KIR3DL1*019 allotype were characterized using the Jurkat T cell line and NKL transfectants. Although poorly expressed on the NK cell surface, KIR3DL1*019 is retained within NK cells, where it induces missing self-recognition of the Bw4 epitope. Transposing our in vitro observations to a cohort of hHSCT patients (N = 186) led us to observe that non-expressed KIR3DL1 HSC grafts increased the incidence of relapse in patients with myeloid diseases. Non-expressed KIR3DL1 alleles could, therefore, influence the outcome of hHSCT.

Software update: Interpreting killer-cell immunoglobulin-like receptors from whole genome sequence data with PING

Here, we demonstrate improvements to our bioinformatic pipeline, PING, which provides high-resolution genotyping of killer-cell immunoglobulin-like receptor (KIR) sequencing data, that expand the method to provide KIR interpretation from whole genome sequencing (WGS) data. We evaluated performance using synthetic sequence datasets and real-world data from the 1000 Genomes Project (1KGP). PING demonstrated high exonic genotyping performance on the synthetic sequence dataset meant to approximate real-world data at 95% accuracy (N = 1366). This result was mirrored in the analysis of 1KGP European data (N = 215) with most genes showing near or below 5% frequency of unresolved exonic genotypes, which is an important indicator for genotyping errors in real-world data. An analysis into the distributions of genotyping errors for the synthetic sequence datasets gave insights into how to further improve genotype accuracy. Similarly, an analysis into ambiguous exonic genotype frequencies for the 1KGP European data, which showed high rates of unresolved genotypes, highlighted that an effective phasing method will be an impactful future additional to the PING workflow. Together, these results demonstrate that PING can effectively provide high-resolution KIR genotyping on WGS data.

Donor selection for adoptive immunotherapy with NK cells in AML patients: Comparison between analysis of lytic NK cell clones and phenotypical identification of alloreactive NK cell repertoire

Natural killer (NK) cell-based adoptive immunotherapy in leukemia patients is an emerging field of interest based on clinical evidence of efficacy and safety. Elderly acute myeloid leukemia (AML) patients have been successfully treated with NK cells from HLA-haploidentical donors, especially when high amounts of alloreactive NK cells were infused. The aim of this study was comparing two approaches to define the size of alloreactive NK cells in haploidentical donors for AML patients recruited in two clinical trials with the acronym "NK-AML" (NCT03955848), and "MRD-NK". The standard methodology was based on the frequency of NK cell clones capable of lysing the related patient-derived cells. The alternative approach consisted of the phenotypic identification of freshly derived NK cells expressing, as inhibitory receptors, only the inhibitory KIR(s) specific for the mismatched KIR-Ligand(s) (HLA-C1, HLA-C2, HLA-Bw4). However, in KIR2DS2⁺ donors and HLA-C1⁺ patients, the unavailability of reagents staining only the inhibitory counterpart (KIR2DL2/L3) may lead to an underestimated identification of the alloreactive NK cell subset. Conversely, in the case of HLA-C1 mismatch, the alloreactive NK cell subset could be overestimated due to the ability of KIR2DL2/L3 to recognize with low-affinity also HLA-C2. Especially in this context, the additional exclusion of LIR1-expressing cells might be relevant to refine the size of the alloreactive NK cell subset. We could also associate degranulation assays, using as effector cells IL-2 activated donor peripheral blood mononuclear cells (PBMC) or NK cells upon co-culture with the related patient target cells. The donor alloreactive NK cell subset always displayed the highest functional activity, confirming its identification accuracy by flow cytometry. Despite the phenotypic limitations and considering the proposed corrective actions, a good correlation was shown by the comparison of the two investigated approaches. In addition, the characterization of receptor expression on a fraction of NK cell clones revealed expected but also few unexpected patterns. Thus, in most instances, the quantification of phenotypically defined alloreactive NK cells from PBMC can provide data similar to the analysis of lytic clones, with several advantages, such as a shorter time to achieve the results and, perhaps, higher reproducibility/feasibility in many laboratories.

Myelodysplastic Syndrome associated TET2 mutations affect NK cell function and genome methylation

Myelodysplastic syndromes (MDS) are clonal hematopoietic disorders, representing high risk of progression to acute myeloid leukaemia, and frequently associated to somatic mutations, notably in the epigenetic regulator TET2. Natural Killer (NK) cells play a role in the anti-leukemic immune response via their cytolytic activity. Here we show that patients with MDS clones harbouring mutations in the TET2 gene are characterised by phenotypic defects in their circulating NK cells. Remarkably, NK cells and MDS clones from the same patient share the TET2 genotype, and the NK cells are characterised by increased methylation of genomic DNA and reduced expression of Killer Immunoglobulin-like receptors (KIR), perforin, and TNF-α. In vitro inhibition of TET2 in NK cells of healthy donors reduces their cytotoxicity, supporting its critical role in NK cell function. Conversely, NK cells from patients treated with azacytidine (#NCT02985190; https://clinicaltrials.gov/ ) show increased KIR and cytolytic protein expression, and IFN-γ production. Altogether, our findings show that, in addition to their oncogenic consequences in the myeloid cell subsets, TET2 mutations contribute to repressing NK-cell function in MDS patients.

Association of KIR2DL5, KIR2DS5, and KIR2DS1 allelic variation and atopic dermatitis

Natural killer cells (NK) have been associated with the pathophysiology of atopic dermatitis (AD). NK function is regulated by killer cell Ig-like receptor family (KIR) receptors that interact with HLA ligands. The study goal was to focus on allelic variation in genes KIR2DL5, KIR2DS5, and KIR2DS1 with respect to AD. This was a case-control study of individuals with (n = 313) and without (n = 176) AD. Associations were estimated using logistic regression. The prevalence of KIR2DL5 was 52.5% (95% CI 48.0,57.0), KIR2DS5 was 33.0% (28.8,37.3), and KIR2DS1 was 33.6% (29.4,38.0). The presence of the KIR2DL5*001:01 increased the odds of having AD by about 86% (odds ratio (OR): 1.86(1.23,2.82) p = 0.003). The risk for individuals homozygous for KIR2DL5*001:01 was even greater (OR: 2.16 (95% CI 1.31,3.53) p = 0.0023). The odds of having AD with KIR2DL5*001:01 was similar in Whites and Blacks. Allelic variation in KIR2DS5 and KIR2DS1 was not associated with AD. There is no known HLA binding ligand for KIR2DL5. The effect of KIR2DL5*001:01 increased in the presence of HLA-B*-21TT leader sequence (2.46(1.37,4.41) p = 0.0025) and the HLA-C2 ligand (2.07 (1.37,4.41, p = 0.000002). Our study shows an independent association of the KIR2DL5*001:01 with AD and is the first study to associate AD with KIR allelic variation.

KIR Allelic Variation and the Remission of Atopic Dermatitis Over Time

Atopic dermatitis (AD) is a common chronic skin disease. Although generally thought to be a disease of T-cell dysregulation, recent studies have suggested that immune dysregulation of NK cells is also important. Killer cell Ig-like receptors (KIRs) are involved with NK cell regulation. The Pediatric Eczema Elective Registry is a U.S. nationwide longitudinal cohort with up to 10 y of follow-up in which 655 children had DNA available for full allelic KIR sequencing. Every 6 mo, AD activity was reported by Pediatric Eczema Elective Registry children. Using generalized estimating equations, we evaluated the association of KIR allelic variation in concert with known HLA binding ligands and whether the child reported AD in "remission" (no skin lesions and not using AD medication). KIR2DS4*001:01 (odds ratio 0.53, 95% CI [0.32, 0.88]) and KIR2DL4*001:02 (0.54, [0.33, 0.89]) in the presence of C*04:01 had the largest effect on decreasing the likelihood of AD remission. The haplotype KIR 2DL4*001:02 ∼ 2DS4*001:01 ∼ 3DL2*002:01 (0.77, [0.60, 0.99]) was also associated with a decreased likelihood of AD remission. Our findings add to the general body of evidence of a growing literature on the importance of NK cells with respect to the immunopathogenesis and natural history of AD.

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.

Allelic Polymorphisms of Killer Immunoglobulin-Like Receptor Genes in Malay and Orang Asli Populations of Peninsular Malaysia

Next-generation DNA sequencing (NGS) technology advancements provide new insight into the level of variation in killer immunoglobulin-like receptor (KIR) genes. High resolution allele genotyping of seven KIR genes was conducted among 94 unrelated Malay and Orang Asli (OA) individuals of Peninsular Malaysia. A manual bioinformatics analysis is performed and optimised by Sanger sequencing method. The Malays expressed a total of 22 alleles, as compared to only 15 alleles in the OA population. In total, 12 centromeric and 9 telomeric allelic haplotypes were identified in the Malays, whereas 8 centromeric and 5 telomeric allelic haplotypes were identified in the OA. The KIR2DL1, KIR2DL3, and KIR2DS4 genes exhibited a high degree of variation and balanced distribution in the Malay and OA populations. On the other hand, KIR2DL4, KIR3DL1, KIR3DL2 and KIR3DL3 genes exhibited a high degree of conservation, with less number of alleles identified and the dominance of a single allele at high frequency. High-resolution KIR allele genotyping has revealed unique sequence variations and allelic haplotypes between individuals and populations. The distributions of KIR alleles and haplotypes are useful for genetic population studies and serve as a baseline for future transplantation matching and disease association research.

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.

High-resolution human KIR genotyping

Killer immunoglobulin-like receptors (KIR) regulate the function of natural killer cells through interactions with various ligands on the surface of cells, thereby determining whether natural killer (NK) cells are to be activated or inhibited from killing the cell being interrogated. The genes encoding these proteins display extensive variation through variable gene content, copy number and allele polymorphism. The combination of KIR genes and their ligands is implicated in various clinical settings including haematopoietic stem cell and solid organ transplant and infectious disease progression. The determination of KIR genes has been used as a factor in the selection of optimal stem cell donors with haplotype variations in recipient and donor giving differential clinical outcomes. Methods to determine KIR genes have primarily involved ascertaining the presence or absence of genes in an individual. With the more recent introduction of massively parallel clonal next-generation sequencing and single molecule very long read length third-generation sequencing, high-resolution determination of KIR alleles has become feasible. Determining the extent and functional impact of allele variation has the potential to lead to further optimisation of clinical outcomes as well as a deeper understanding of the functional properties of the receptors and their interactions with ligands. This review summarizes recently published high-resolution KIR genotyping methods and considers the various advantages and disadvantages of the approaches taken. In addition the application of allele level genotyping in the setting of transplantation and infectious disease control is discussed.

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-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.