Co-evolution of human leukocyte antigen (HLA) class I ligands with killer-cell immunoglobulin-like receptors (KIR) in a genetically diverse population of sub-Saharan Africans

Unveiling the Significance of HLA and KIR Diversity in Underrepresented Populations

Human leukocyte antigen (HLA) molecules and their relationships with natural killer (NK) cells, specifically through their interaction with killer-cell immunoglobulin-like receptors (KIRs), exhibit robust associations with the outcomes of diverse diseases. Moreover, genetic variations in HLA and KIR immune system genes offer limitless depths of complexity. In recent years, a surge of high-powered genome-wide association studies (GWASs) utilizing single nucleotide polymorphism (SNP) arrays has occurred, significantly advancing our understanding of disease pathogenesis. Additionally, advances in HLA reference panels have enabled higher resolution and more reliable imputation, allowing for finer-grained evaluation of the association between sequence variations and disease risk. However, it is essential to note that the majority of these GWASs have focused primarily on populations of Caucasian and Asian origins, neglecting underrepresented populations in Latin America and Africa. This omission not only leads to disparities in health care access but also restricts our knowledge of novel genetic variants involved in disease pathogenesis within these overlooked populations. Since the KIR and HLA haplotypes prevalent in each population are clearly modelled by the specific environment, the aim of this review is to encourage studies investigating HLA/KIR involvement in infection and autoimmune diseases, reproduction, and transplantation in underrepresented populations.

25 years of the IPD-IMGT/HLA Database

Twenty-five years ago, in 1998, the HLA Informatics Group of the Anthony Nolan Research Institute released the IMGT/HLA Database. Since this time, this online resource has acted as the repository for the numerous variant sequences of HLA alleles named by the WHO Nomenclature Committee for Factors of the HLA System. The IPD-IMGT/HLA Database has provided a stable, highly accessible, user-friendly repository for this work. During this time, the technology underlying HLA typing has undergone significant changes. Next generation sequencing (NGS) has superseded previous methodologies of HLA typing and can generate large amounts of high-resolution sequencing data. This has resulted in a drastic increase in the number and complexity of sequences submitted to the database. The challenge for the IPD-IMGT/HLA Database has been to maintain the highest standards of curation, while supporting the core set of tools and functionality to our users with increased numbers of submissions and sequences. Traditional methods of accessing and presenting data have been challenged and new methods utilising new computing technologies have had to be developed to keep pace and support a shifting user demographic.

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.

Development of cost-effective and fast KIR genotyping by multiplex PCR-SSP

Killer-cell immunoglobulin-like receptors (KIR) control natural killer (NK) cell functions by recognizing HLA molecules and modulating the activity of NK cells. The KIR gene cluster contains polymorphic and highly homologous genes. Diversity of the KIR region is achieved through differences in gene content, allelic polymorphism, and gene copy number, which result in unrelated individuals having different KIR genotypes and individualized immune responses that are relevant to multiple aspects of human health and disease. Therefore, KIR genotyping is increasingly used in epidemiological studies. Here, we developed multiplex polymerase chain reaction with sequence-specific primers (PCR-SSP) to compensate for the shortcomings of the conventional PCR-SSP method, which is most commonly used for KIR analysis. Multiplex PCR-SSP method involves six multiplex reactions that detect 16 KIR genes and distinguish variant types of some KIR genes by adding two reactions. The assay was evaluated in a blind survey using a panel of 40 reference DNA standards from the UCLA KIR Exchange Program. The results are 100% concordant with the genotype determined using Luminex-based reverse sequence-specific oligonucleotide typing systems. Additionally, we investigated the currently known 16 KIR genes and their common variants in 120 unrelated Korean individuals. The results were consistent with the KIR genotype previously reported by Hwang et al. This multiplex PCR-SSP is an efficient method for analyzing KIR genotypes in both small- and large-scale studies with minimal labor, reagents, and DNA. Furthermore, by providing a better definition of KIR polymorphisms it can contribute to developments in immunogenetics.

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.

Malaria-driven adaptation of MHC class I in wild bonobo populations

The malaria parasite Plasmodium falciparum causes substantial human mortality, primarily in equatorial Africa. Enriched in affected African populations, the B*53 variant of HLA-B, a cell surface protein that presents peptide antigens to cytotoxic lymphocytes, confers protection against severe malaria. Gorilla, chimpanzee, and bonobo are humans' closest living relatives. These African apes have HLA-B orthologs and are infected by parasites in the same subgenus (Laverania) as P. falciparum, but the consequences of these infections are unclear. Laverania parasites infect bonobos (Pan paniscus) at only one (TL2) of many sites sampled across their range. TL2 spans the Lomami River and has genetically divergent subpopulations of bonobos on each side. Papa-B, the bonobo ortholog of HLA-B, includes variants having a B*53-like (B07) peptide-binding supertype profile. Here we show that B07 Papa-B occur at high frequency in TL2 bonobos and that malaria appears to have independently selected for different B07 alleles in the two subpopulations.

Characterization of the immunoglobulin lambda chain locus from diverse populations reveals extensive genetic variation

Immunoglobulins (IGs), crucial components of the adaptive immune system, are encoded by three genomic loci. However, the complexity of the IG loci severely limits the effective use of short read sequencing, limiting our knowledge of population diversity in these loci. We leveraged existing long read whole-genome sequencing (WGS) data, fosmid technology, and IG targeted single-molecule, real-time (SMRT) long-read sequencing (IG-Cap) to create haplotype-resolved assemblies of the IG Lambda (IGL) locus from 6 ethnically diverse individuals. In addition, we generated 10 diploid assemblies of IGL from a diverse cohort of individuals utilizing IG-Cap. From these 16 individuals, we identified significant allelic diversity, including 36 novel IGLV alleles. In addition, we observed highly elevated single nucleotide variation (SNV) in IGLV genes relative to IGL intergenic and genomic background SNV density. By comparing SNV calls between our high quality assemblies and existing short read datasets from the same individuals, we show a high propensity for false-positives in the short read datasets. Finally, for the first time, we nucleotide-resolved common 5-10 Kb duplications in the IGLC region that contain functional IGLJ and IGLC genes. Together these data represent a significant advancement in our understanding of genetic variation and population diversity in the IGL locus.

Analysis of KIR gene variants in The Cancer Genome Atlas and UK Biobank using KIRCLE

BACKGROUND

Natural killer (NK) cells represent a critical component of the innate immune system's response against cancer and viral infections, among other diseases. To distinguish healthy host cells from infected or tumor cells, killer immunoglobulin receptors (KIR) on NK cells bind and recognize Human Leukocyte Antigen (HLA) complexes on their target cells. However, NK cells exhibit great diversity in their mechanism of activation, and the outcomes of their activation are not yet understood fully. Just like the HLAs they bind, KIR receptors exhibit high allelic diversity in the human population. Here we provide a method to identify KIR allele variants from whole exome sequencing data and uncover novel associations between these variants and various molecular and clinical correlates.

RESULTS

In order to better understand KIRs, we have developed KIRCLE, a novel method for genotyping individual KIR genes from whole exome sequencing data, and used it to analyze approximately sixty-thousand patient samples in The Cancer Genome Atlas (TCGA) and UK Biobank. We were able to assess population frequencies for different KIR alleles and demonstrate that, similar to HLA alleles, individuals' KIR alleles correlate strongly with their ethnicities. In addition, we observed associations between different KIR alleles and HLA alleles, including HLA-B*53 with KIR3DL2*013 (Fisher's exact FDR = 7.64e-51). Finally, we showcased statistically significant associations between KIR alleles and various clinical correlates, including peptic ulcer disease (Fisher's exact FDR = 0.0429) and age of onset of atopy (Mann-Whitney U FDR = 0.0751).

CONCLUSIONS

We show that KIRCLE is able to infer KIR variants accurately and consistently, and we demonstrate its utility using data from approximately sixty-thousand individuals from TCGA and UK Biobank to discover novel molecular and clinical correlations with KIR germline variants. Peptic ulcer disease and atopy are just two diseases in which NK cells may play a role beyond their "classical" realm of anti-tumor and anti-viral responses. This tool may be used both as a benchmark for future KIR-variant-inference algorithms, and to better understand the immunogenomics of and disease processes involving KIRs.

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.

Donor with HLA-C2 is associated with acute rejection following liver transplantation in southern Chinese

Apart from presenting peptides to T cells, class I HLA molecules serve as ligands for KIRs and regulate the response of NK cells. The role played by HLA and KIR in the acute rejection (AR) following liver transplantation has been controversial. In this retrospective study, we assessed the influence of class I HLA alleles, HLA matching between donor-recipient pairs, recipient KIR and donor HLA ligands on AR following liver transplantation in southern Chinese. 143 recipients and 78 donors obtained from a single transplant center were included in the study cohort. 33 recipients with histologically confirmed AR were observed. We found that the incidence of AR did not correlate with donor or recipient class I HLA alleles and HLA matching. Neither recipient KIR gene nor the KIR genotype was associated with AR, moreover, high-resolution genotyping of 14 functional KIR genes of recipients showed that no KIR allele was independently associated with AR. However, the frequency of HLA-C2⁺ donor significantly increased in AR group compared with NAR group (52.9% vs. 24.6%, P = 0.03). In the presence of HLA-C2 by the donor allograft, AR was more frequently observed in recipients with normal expressed KIR2DS4 (43.8% vs. 15.0%, P = 0.03). Donor with HLA-C2 is therefore a major determinant of AR, which can confer risk effect in liver transplantation. Our findings can provide valuable clues for better understanding pathogenesis of AR and have important clinical implications in liver transplantation for Chinese. This article is protected by copyright. All rights reserved.

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.

High Resolution HLA ∼A, ∼B, ∼C, ∼DRB1, ∼DQA1, and ∼DQB1 Diversity in South African Populations

Background: Lack of HLA data in southern African populations hampers disease association studies and our understanding of genetic diversity in these populations. We aimed to determine HLA diversity in South African populations using high resolution HLA ∼A, ∼B, ∼C, ∼DRB1, ∼DQA1 and ∼DQB1 data, from 3005 previously typed individuals.

Methods: We determined allele and haplotype frequencies, deviations from Hardy-Weinberg equilibrium (HWE), linkage disequilibrium (LD) and neutrality test. South African HLA class I data was additionally compared to other global populations using non-metrical multidimensional scaling (NMDS), genetic distances and principal component analysis (PCA).

Results: All loci strongly (p < 0.0001) deviated from HWE, coupled with excessive heterozygosity in most loci. Two of the three most frequent alleles, HLA ∼DQA1*05:02 (0.2584) and HLA ∼C*17:01 (0.1488) were previously reported in South African populations at lower frequencies. NMDS showed genetic distinctness of South African populations. Phylogenetic analysis and PCA clustered our current dataset with previous South African studies. Additionally, South Africans seem to be related to other sub-Saharan populations using HLA class I allele frequencies.

Discussion and Conclusion: Despite the retrospective nature of the study, data missingness, the imbalance of sample sizes for each locus and haplotype pairs, and induced methodological difficulties, this study provides a unique and large HLA dataset of South Africans, which might be a useful resource to support anthropological studies, disease association studies, population based vaccine development and donor recruitment programs. We additionally provide simulated high resolution HLA class I data to augment the mixed resolution typing results generated from this study.

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.

Heterogeneity of Genetic Admixture Determines SLE Susceptibility in Mexican

Systemic Lupus Erythematosus (SLE) is an autoimmune inflammatory disorder for which Major Histocompatibility Complex (MHC) genes are well identified as risk factors. SLE patients present different clinical phenotypes, which are partly explained by admixture patterns variation among Mexicans. Population genetic has insight into the high genetic variability of Mexicans, mainly described through HLA gene studies with anthropological and biomedical importance. A prospective, case-control study was performed. In this study, we recruited 146 SLE patients, and 234 healthy individuals were included as a control group; both groups were admixed Mexicans from Mexico City. The HLA typing methods were based on Next Generation Sequencing and Sequence-Based Typing (SBT). The data analysis was performed with population genetic programs and statistical packages. The admixture estimations based on HLA-B and -DRB1 revealed that SLE patients have a higher Southwestern European ancestry proportion (48 ± 8%) than healthy individuals (30 ± 7%). In contrast, Mexican Native American components are diminished in SLE patients (44 ± 1%) and augmented in Healthy individuals (63 ± 4%). HLA alleles and haplotypes' frequency analysis found variants previously described in SLE patients from Mexico City. Moreover, a conserved extended haplotype that confers risk to develop SLE was found, the HLA-A^∗29:02∼C^∗16:01∼B^∗44:03∼DRB1^∗07:01∼DQB1^∗02:02, pC = 0.02, OR = 1.41. Consistent with the admixture estimations, the origin of all risk alleles and haplotypes found in this study are European, while the protection alleles are Mexican Native American. The analysis of genetic distances supported that the SLE patient group is closer to the Southwestern European parental populace and farthest from Mexican Native Americans than healthy individuals. Heterogeneity of genetic admixture determines SLE susceptibility and protection in Mexicans. HLA sequencing is helpful to determine susceptibility alleles and haplotypes restricted to some populations.

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.

HLAs, TCRs, and KIRs, a Triumvirate of Human Cell-Mediated Immunity

In all human cells, human leukocyte antigen (HLA) class I glycoproteins assemble with a peptide and take it to the cell surface for surveillance by lymphocytes. These include natural killer (NK) cells and γδ T cells of innate immunity and αβ T cells of adaptive immunity. In healthy cells, the presented peptides derive from human proteins, to which lymphocytes are tolerant. In pathogen-infected cells, HLA class I expression is perturbed. Reduced HLA class I expression is detected by KIR and CD94:NKG2A receptors of NK cells. Almost any change in peptide presentation can be detected by αβ CD8⁺ T cells. In responding to extracellular pathogens, HLA class II glycoproteins, expressed by specialized antigen-presenting cells, present peptides to αβ CD4⁺ T cells. In comparison to the families of major histocompatibility complex (MHC) class I, MHC class II and αβ T cell receptors, the antigenic specificity of the γδ T cell receptors is incompletely understood.

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.

Search for MHC/TCR-Like Systems in Living Organisms

Highly polymorphic loci evolved many times over the history of species. These polymorphic loci are involved in three types of functions: kind recognition, self-incompatibility, and the jawed vertebrate adaptive immune system (AIS). In the first part of this perspective, we reanalyzed and described some cases of polymorphic loci reported in the literature. There is a convergent evolution within each functional category and between functional categories, suggesting that the emergence of these self/non-self recognition loci has occurred multiple times throughout the evolutionary history. Most of the highly polymorphic loci are coding for proteins that have a homophilic interaction or heterophilic interaction between linked loci, leading to self or non-self-recognition. The highly polymorphic MHCs, which are involved in the AIS have a different functional mechanism, as they interact through presented self or non-self-peptides with T cell receptors, whose diversity is generated by somatic recombination. Here we propose a mechanism called “the capacity of recognition competition mechanism” that might contribute to the evolution of MHC polymorphism. We propose that the published cases corresponding to these three biological categories represent a small part of what can be found throughout the tree of life, and that similar mechanisms will be found many times, including the one where polymorphic loci interact with somatically generated loci.

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.

Haplotype Motif-Based Models for KIR-Genotype Informed Selection of Hematopoietic Cell Donors Fail to Predict Outcome of Patients With Myelodysplastic Syndromes or Secondary Acute Myeloid Leukemia

Results from registry studies suggest that harnessing Natural Killer (NK) cell reactivity mediated through Killer cell Immunoglobulin-like Receptors (KIR) could reduce the risk of relapse after allogeneic Hematopoietic Cell Transplantation (HCT). Several competing models have been developed to classify donors as KIR-advantageous or disadvantageous. Basically, these models differ by grouping donors based on distinct KIR–KIR–ligand combinations or by haplotype motif assignment. This study aimed to validate different models for unrelated donor selection for patients with Myelodysplatic Syndromes (MDS) or secondary Acute Myeloid Leukemia (sAML). In a joint retrospective study of the European Society for Blood and Marrow Transplantation (EBMT) and the Center for International Blood and Marrow Transplant Research (CIBMTR) registry data from 1704 patients with secondary AML or MDS were analysed. The cohort consisted mainly of older patients (median age 61 years) with high risk disease who had received chemotherapy-based reduced intensity conditioning and anti-thymocyte globulin prior to allogeneic HCT from well-matched unrelated stem cell donors. The impact of the predictors on Overall Survival (OS) and relapse incidence was tested in Cox regression models adjusted for patient age, a modified disease risk index, performance status, donor age, HLA-match, sex-match, CMV-match, conditioning intensity, type of T-cell depletion and graft type. KIR genes were typed using high-resolution amplicon-based next generation sequencing. In univariable and multivariable analyses none of the models predicted OS and the risk of relapse consistently. Our results do not support the hypothesis that optimizing NK-mediated alloreactivity is possible by KIR-genotype informed selection of HLA-matched unrelated donors. However, in the context of allogeneic transplantation, NK-cell biology is complex and only partly understood. KIR-genes are highly diverse and current assignment of haplotype motifs based on the presence or absence of selected KIR genes is over-simplistic. As a consequence, further research is highly warranted and should integrate cutting edge knowledge on KIR genetics, and NK-cell biology into future studies focused on homogeneous groups of patients and treatment modalities.

Identification of a tumor-specific allo-HLA-restricted γδTCR

γδT cells are key players in cancer immune surveillance because of their ability to recognize malignant transformed cells, which makes them promising therapeutic tools in the treatment of cancer. However, the biological mechanisms of how γδT-cell receptors (TCRs) interact with their ligands are poorly understood. Within this context, we describe the novel allo-HLA-restricted and CD8α-dependent Vγ5Vδ1TCR. In contrast to the previous assumption of the general allo-HLA reactivity of a minor fraction of γδTCRs, we show that classic anti-HLA-directed, γδTCR-mediated reactivity can selectively act on hematological and solid tumor cells, while not harming healthy tissues in vitro and in vivo. We identified the molecular interface with proximity to the peptide-binding groove of HLA-A*24:02 as the essential determinant for recognition and describe the critical role of CD8 as a coreceptor. We conclude that alloreactive γδT-cell repertoires provide therapeutic opportunities, either within the context of haplotransplantation or as individual γδTCRs for genetic engineering of tumor-reactive T cells.

Hla-C genetic diversity and evolutionary insights in two samples from Brazil and Benin

Human leukocyte antigen-C (HLA-C) is a classical HLA class I molecule that binds and presents peptides to cytotoxic T lymphocytes in the cell surface. HLA-C has a dual function because it also interacts with Killer-cell immunoglobulin-like receptors (KIR) receptors expressed in natural killer and T cells, modulating their activity. The structure and diversity of the HLA-C regulatory regions, as well as the relationship among variants along the HLA-C locus, are poorly addressed, and few population-based studies explored the HLA-C variability in the entire gene in different population samples. Here we present a molecular and bioinformatics method to evaluate the entire HLA-C diversity, including regulatory sequences. Then, we applied this method to survey the HLA-C diversity in two population samples with different demographic histories, one highly admixed from Brazil with major European contribution, and one from Benin with major African contribution. The HLA-C promoter and 3'UTR were very polymorphic with the presence of few, but highly divergent haplotypes. These segments also present conserved sequences that are shared among different primate species. Nucleotide diversity was higher in other segments rather than exons 2 and 3, particularly around exon 5 and the second half of the 3'UTR region. We detected evidence of balancing selection on the entire HLA-C locus and positive selection in the HLA-C leader peptide, for both populations. HLA-C motifs previously associated with KIR interaction and expression regulation are similar between both populations. Each allele group is associated with specific regulatory sequences, reflecting the high linkage disequilibrium along the entire HLA-C locus in both populations.

Variations in killer-cell immunoglobulin-like receptor and human leukocyte antigen genes and immunity to malaria

Malaria is one of the deadliest infectious diseases in the world. Immune responses to Plasmodium falciparum malaria vary among individuals and between populations. Human genetic variation in immune system genes is likely to play a role in this heterogeneity. Natural killer (NK) cells produce inflammatory cytokines in response to malaria infection, kill intraerythrocytic Plasmodium falciparum parasites by cytolysis, and participate in the initiation and development of adaptive immune responses to plasmodial infection. These functions are modulated by interactions between killer-cell immunoglobulin-like receptors (KIRs) and human leukocyte antigens (HLAs). Therefore, variations in KIR and HLA genes can have a direct impact on NK cell functions. Understanding the role of KIRs and HLAs in immunity to malaria can help to better characterize antimalarial immune responses. In this review, we summarize the different KIRs and HLAs associated with immunity to malaria thus far.

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.

IPD-IMGT/HLA Database

The IPD-IMGT/HLA Database, http://www.ebi.ac.uk/ipd/imgt/hla/, currently contains over 25 000 allele sequence for 45 genes, which are located within the Major Histocompatibility Complex (MHC) of the human genome. This region is the most polymorphic region of the human genome, and the levels of polymorphism seen exceed most other genes. Some of the genes have several thousand variants and are now termed hyperpolymorphic, rather than just simply polymorphic. The IPD-IMGT/HLA Database has provided a stable, highly accessible, user-friendly repository for this information, providing the scientific and medical community access to the many variant sequences of this gene system, that are critical for the successful outcome of transplantation. The number of currently known variants, and dramatic increase in the number of new variants being identified has necessitated a dedicated resource with custom tools for curation and publication. The challenge for the database is to continue to provide a highly curated database of sequence variants, while supporting the increased number of submissions and complexity of sequences. In order to do this, traditional methods of accessing and presenting data will be challenged, and new methods will need to be utilized to keep pace with new discoveries.

Genetic diversity affects the nanoscale membrane organization and signaling of natural killer cell receptors

Genetic diversity in human natural killer (NK) cell receptors is linked to resistance and susceptibility to many diseases. Here, we tested the effect of this diversity on the nanoscale organization of killer cell immunoglobulin-like receptors (KIRs). Using superresolution microscopy, we found that inhibitory KIRs encoded by different genes and alleles were organized differently at the surface of primary human NK cells. KIRs that were found at low abundance assembled into smaller clusters than those formed by KIRs that were more highly abundant, and at low abundance, there was a greater proportion of KIRs in clusters. Upon receptor triggering, a structured interface called the immune synapse assembles, which facilitates signal integration and controls NK cell responses. Here, triggering of low-abundance receptors resulted in less phosphorylation of the downstream phosphatase SHP-1 but more phosphorylation of the adaptor protein Crk than did triggering of high-abundance receptors. In cells with greater KIR abundance, SHP-1 dephosphorylated Crk, which potentiated NK cell spreading during activation. Thus, genetic variation modulates both the abundance and nanoscale organization of inhibitory KIRs. That is, as well as the number of receptors at the cell surface varying with genotype, the way in which these receptors are organized in the membrane also varies. Essentially, a change in the average surface abundance of a protein at the cell surface is a coarse descriptor entwined with changes in local nanoscale clustering. Together, our data indicate that genetic diversity in inhibitory KIRs affects membrane-proximal signaling and, unexpectedly, the formation of activating immune synapses.

Differentiation and adaptation of natural killer cells for anti-malarial immunity

Natural killer cells employ a diverse arsenal of effector mechanisms to target intracellular pathogens. Differentiation of natural killer (NK) cell activation pathways occurs along a continuum from reliance on innate pro-inflammatory cytokines and stress-induced host ligands through to interaction with signals derived from acquired immune responses. Importantly, the degree of functional differentiation of the NK cell lineage influences the magnitude and specificity of interactions with host cells infected with viruses, bacteria, fungi, and parasites. Individual humans possess a vast diversity of distinct NK cell clones, each with the capacity to vary along this functional differentiation pathway, which - when combined - results in unique individual responses to different infections. Here we summarize these NK cell differentiation events, review evidence for direct interaction of malaria-infected host cells with NK cells and assess how innate inflammatory signals induced by malaria parasite-associated molecular patterns influence the indirect activation and function of NK cells. Finally, we discuss evidence that anti-malarial immunity develops in parallel with advancing NK differentiation, coincident with a loss of reliance on inflammatory signals, and a refined capacity of NK cells to target malaria parasites more precisely, particularly through antibody-dependent mechanisms.

Distribution of Killer-Cell Immunoglobulin-Like Receptor Genes and Combinations of Their Human Leucocyte Antigen Ligands in 11 Ethnic Populations in China

The aim of this study was to analyze the distribution of killer-cell immunoglobulin-like receptor (KIR) genes and their human leucocyte antigen (HLA) ligand combinations in different original ethnic populations in China, and thus, to provide relevant genomic diversity data for the future study of viral infections, autoimmune diseases, and reproductive fitness. A total of 1119 unrelated individuals from 11 ethnic populations—including Hani, Jinuo, Lisu, Nu, Bulang, Wa, Dai, Maonan, Zhuang, Tu, and Yugu—from four original groups, were included. The presence/absence of the 16 KIR loci were detected, and the KIR gene’s phenotype, genotype, and haplotype A and B frequencies, as well as KIR ligand’s HLA allotype and KIR–HLA pairs for each population, were calculated. Principal component analysis and phylogenetic trees were constructed to compare the characteristics of the KIR and KIR–HLA pair distributions of these 11 populations. In total, 92 KIR genotypes were identified, including six new genotypes. The KIR and its HLA ligands had a distributed diversity in 11 ethnic populations in China, and each group had its specific KIR and KIR–HLA pair profile. The difference among the KIR–HLA pairs between northern and southern groups, but not among the four original groups, may reflect strong pressure from previous or ongoing infectious diseases, which have a significant impact on KIR and its HLA combination repertoires.

Human NK cell receptor KIR2DS4 detects a conserved bacterial epitope presented by HLA-C

Natural killer (NK) cells have an important role in immune defense against viruses and cancer. Activation of human NK cell cytotoxicity toward infected or tumor cells is regulated by killer cell immunoglobulin-like receptors (KIRs) that bind to human leukocyte antigen class I (HLA-I). Combinations of KIR with HLA-I are genetically associated with susceptibility to disease. KIR2DS4, an activating member of the KIR family with poorly defined ligands, is a receptor of unknown function. Here, we show that KIR2DS4 has a strong preference for rare peptides carrying a Trp at position 8 (p8) of 9-mer peptides bound to HLA-C*05:01. The complex of a peptide bound to HLA-C*05:01 with a Trp at p8 was sufficient for activation of primary KIR2DS4+ NK cells, independent of activation by other receptors and of prior NK cell licensing. HLA-C*05:01+ cells that expressed the peptide epitope triggered KIR2DS4+ NK cell degranulation. We show an inverse correlation of the worldwide allele frequency of functional KIR2DS4 with that of HLA-C*05:01, indicative of functional interaction and balancing selection. We found a highly conserved peptide sequence motif for HLA-C*05:01-restricted activation of human KIR2DS4+ NK cells in bacterial recombinase A (RecA). KIR2DS4+ NK cells were stimulated by RecA epitopes from multiple human pathogens, including Helicobacter, Chlamydia, Brucella, and Campylobacter. We predict that over 1,000 bacterial species could activate NK cells through KIR2DS4, and propose that human NK cells also contribute to immune defense against bacteria through recognition of a conserved RecA epitope presented by HLA-C*05:01.

Fluctuating and Geographically Specific Selection Characterize Rapid Evolution of the Human KIR Region

The killer-cell immunoglobulin-like receptor (KIR) region comprises a fast-evolving family of genes that encode receptors for natural killer (NK) cells and have crucial role in host defense. Evolution of KIR was examined in the context of the human genome. Gene-content diversity and single nucleotide polymorphisms (SNP) in the KIR genes and flanking regions were compared to >660,000 genome-wide SNPs in over 800 individuals from 52 populations of the human genome diversity panel (HGDP). KIR allelic diversity was further examined using next generation sequencing in a subset of 56 individuals. We identified the SNP rs587560 located in KIR3DL3 as a marker of KIR2DL2 and KIR2DL3 and, consequently, Cen A and Cen B haplotypes. We also show that combinations of two KIR2DL4 SNPs (rs35656676 and rs592645) distinguish KIR3DL1 from KIR3DS1 and also define the major KIR3DL1 high- and low-expressing alleles lineages. Comparing the diversity of the SNPs within the KIR region to remainder of the genome, we observed a high diversity for the centromeric KIR region consistent with balancing selection (p < 0.01); in contrast, centromeric KIR diversity is significantly reduced in East Asian populations (p < 0.01), indicating purifying selection. By analyzing SNP haplotypes in a region spanning ~500 kb that includes the KIR cluster, we observed evidence of strong positive selection in Africa for high-expressing KIR3DL1 alleles, favored over the low-expressing alleles (p < 0.01). In sharp contrast, the strong positive selection (p < 0.01) that we also observed in the telomeric KIR region in Oceanic populations tracked with a high frequency of KIR3DS1. In addition, we demonstrated that worldwide frequency of high-expression KIR3DL1 alleles was correlated with virus with virus (r = 0.64, p < 10⁻⁶) and protozoa (r = 0.69, p < 10⁻⁶) loads, which points to selection globally on KIR3DL1 high-expressing alleles attributable to pathogen exposure.

Diversity of KIR, HLA Class I, and Their Interactions in Seven Populations of Sub-Saharan Africans

HLA class I and KIR sequences were determined for Dogon, Fulani, and Baka populations of western Africa, Mbuti of central Africa, and Datooga, Iraqw, and Hadza of eastern Africa. Study of 162 individuals identified 134 HLA class I alleles (41 HLA-A, 60 HLA-B, and 33 HLA-C). Common to all populations are three HLA-C alleles (C1⁺C*07:01, C1⁺C*07:02, and C2⁺C*06:02) but no HLA-A or -B Unexpectedly, no novel HLA class I was identified in these previously unstudied and anthropologically distinctive populations. In contrast, of 227 KIR detected, 22 are present in all seven populations and 28 are novel. A high diversity of HLA A-C-B haplotypes was observed. In six populations, most haplotypes are represented just once. But in the Hadza, a majority of haplotypes occur more than once, with 2 having high frequencies and 10 having intermediate frequencies. The centromeric (cen) part of the KIR locus exhibits an even balance between cenA and cenB in all seven populations. The telomeric (tel) part has an even balance of telA to telB in East Africa, but this changes across the continent to where telB is vestigial in West Africa. All four KIR ligands (A3/11, Bw4, C1, and C2) are present in six of the populations. HLA haplotypes of the Iraqw and Hadza encode two KIR ligands, whereas the other populations have an even balance between haplotypes encoding one and two KIR ligands. Individuals in these African populations have a mean of 6.8-8.4 different interactions between KIR and HLA class I, compared with 2.9-6.5 for non-Africans.

qKAT: Quantitative Semi-automated Typing of Killer-cell Immunoglobulin-like Receptor Genes

Killer cell immunoglobulin-like receptors (KIRs) are a set of inhibitory and activating immune receptors, on natural killer (NK) and T cells, encoded by a polymorphic cluster of genes on chromosome 19. Their best-characterized ligands are the human leukocyte antigen (HLA) molecules that are encoded within the major histocompatibility complex (MHC) locus on chromosome 6. There is substantial evidence that they play a significant role in immunity, reproduction, and transplantation, making it crucial to have techniques that can accurately genotype them. However, high-sequence homology, as well as allelic and copy number variation, make it difficult to design methods that can accurately and efficiently genotype all KIR genes. Traditional methods are usually limited in the resolution of data obtained, throughput, cost-effectiveness, and the time taken for setting up and running the experiments. We describe a method called quantitative KIR semi-automated typing (qKAT), which is a high-throughput multiplex real-time polymerase chain reaction method that can determine the gene copy numbers for all genes in the KIR locus. qKAT is a simple high-throughput method that can provide high-resolution KIR copy number data, which can be further used to infer the variations in the structurally polymorphic haplotypes that encompass them. This copy number and haplotype data can be beneficial for studies on large-scale disease associations, population genetics, as well as investigations on expression and functional interactions between KIR and HLA.

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.

Conserved and variable natural killer cell receptors: diverse approaches to viral infections

Natural killer (NK) cells are lymphocytes of the innate immune system with essential roles during viral infections. NK cell functions are mediated through a repertoire of non-rearranging inhibitory and activating receptors that interact with major histocompatibility complex (MHC)-peptide complexes on the surface of infected cells. Recent work studying the conserved CD94-NKG2A and variable killer cell immunoglobulin-like receptor-MHC systems suggest that these two receptor families may have subtly different properties in terms of interactions with MHC class I bound peptides, and in recognition of down-regulation of MHC class I. In this review, we discuss how these properties generate diversity in the NK cell response to viruses.

Allele-Level KIR Genotyping of More Than a Million Samples: Workflow, Algorithm, and Observations

The killer-cell immunoglobulin-like receptor (KIR) genes regulate natural killer cell activity, influencing predisposition to immune mediated disease, and affecting hematopoietic stem cell transplantation (HSCT) outcome. Owing to the complexity of the KIR locus, with extensive gene copy number variation (CNV) and allelic diversity, high-resolution characterization of KIR has so far been applied only to relatively small cohorts. Here, we present a comprehensive high-throughput KIR genotyping approach based on next generation sequencing. Through PCR amplification of specific exons, our approach delivers both copy numbers of the individual genes and allelic information for every KIR gene. Ten-fold replicate analysis of a set of 190 samples revealed a precision of 99.9%. Genotyping of an independent set of 360 samples resulted in an accuracy of more than 99% taking into account consistent copy number prediction. We applied the workflow to genotype 1.8 million stem cell donor registry samples. We report on the observed KIR allele diversity and relative abundance of alleles based on a subset of more than 300,000 samples. Furthermore, we identified more than 2,000 previously unreported KIR variants repeatedly in independent samples, underscoring the large diversity of the KIR region that awaits discovery. This cost-efficient high-resolution KIR genotyping approach is now applied to samples of volunteers registering as potential donors for HSCT. This will facilitate the utilization of KIR as additional selection criterion to improve unrelated donor stem cell transplantation outcome. In addition, the approach may serve studies requiring high-resolution KIR genotyping, like population genetics and disease association studies.

Inhibitory killer cell immunoglobulin-like receptors strengthen CD8+ T cell-mediated control of HIV-1, HCV, and HTLV-1

Killer cell immunoglobulin-like receptors (KIRs) are expressed predominantly on natural killer cells, where they play a key role in the regulation of innate immune responses. Recent studies show that inhibitory KIRs can also affect adaptive T cell-mediated immunity. In mice and in human T cells in vitro, inhibitory KIR ligation enhanced CD8+ T cell survival. To investigate the clinical relevance of these observations, we conducted an extensive immunogenetic analysis of multiple independent cohorts of HIV-1-, hepatitis C virus (HCV)-, and human T cell leukemia virus type 1 (HTLV-1)-infected individuals in conjunction with in vitro assays of T cell survival, analysis of ex vivo KIR expression, and mathematical modeling of host-virus dynamics. Our data suggest that functional engagement of inhibitory KIRs enhances the CD8+ T cell response against HIV-1, HCV, and HTLV-1 and is a significant determinant of clinical outcome in all three viral infections.

High-Resolution Genetic and Phenotypic Analysis of KIR2DL1 Alleles and Their Association with Pre-Eclampsia

Killer-cell Ig-like receptor (KIR) genes are inherited as haplotypes. They are expressed by NK cells and linked to outcomes of infectious diseases and pregnancy in humans. Understanding how genotype relates to phenotype is difficult because of the extensive diversity of the KIR family. Indeed, high-resolution KIR genotyping and phenotyping in single NK cells in the context of disease association is lacking. In this article, we describe a new method to separate NK cells expressing allotypes of the KIR2DL1 gene carried by the KIR A haplotype (KIR2DL1A) from those expressing KIR2DL1 alleles carried by the KIR B haplotype (KIR2DL1B). We find that in KIR AB heterozygous individuals, different KIR2DL1 allotypes can be detected in both peripheral blood and uterine NK cells. Using this new method, we demonstrate that both blood and uterine NK cells codominantly express KIR2DL1A and KIR2DL1B allotypes but with a predominance of KIR2DL1A variants, which associate with enhanced NK cell function. In a case-control study of pre-eclampsia, we show that KIR2DL1A, not KIR2DL1B, associates with increased disease risk. This method will facilitate our understanding of how individual KIR2DL1 allelic variants affect NK cell function and contribute to disease risk.

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.

Applications of Immunopharmacogenomics: Predicting, Preventing, and Understanding Immune-Mediated Adverse Drug Reactions

Adverse drug reactions (ADRs) are a significant health care burden. Immune-mediated adverse drug reactions (IM-ADRs) are responsible for one-fifth of ADRs but contribute a disproportionately high amount of that burden due to their severity. Variation in human leukocyte antigen ( HLA) genes has emerged as a potential preprescription screening strategy for the prevention of previously unpredictable IM-ADRs. Immunopharmacogenomics combines the disciplines of immunogenomics and pharmacogenomics and focuses on the effects of immune-specific variation on drug disposition and IM-ADRs. In this review, we present the latest evidence for HLA associations with IM-ADRs, ongoing research into biological mechanisms of IM-ADRs, and the translation of clinical actionable biomarkers for IM-ADRs, with a focus on T cell-mediated ADRs.

Evolutionary history and adaptation of a human pygmy population of Flores Island, Indonesia

Flores Island, Indonesia, was inhabited by the small-bodied hominin species Homo floresiensis, which has an unknown evolutionary relationship to modern humans. This island is also home to an extant human pygmy population. Here we describe genome-scale single-nucleotide polymorphism data and whole-genome sequences from a contemporary human pygmy population living on Flores near the cave where H. floresiensis was found. The genomes of Flores pygmies reveal a complex history of admixture with Denisovans and Neanderthals but no evidence for gene flow with other archaic hominins. Modern individuals bear the signatures of recent positive selection encompassing the FADS (fatty acid desaturase) gene cluster, likely related to diet, and polygenic selection acting on standing variation that contributed to their short-stature phenotype. Thus, multiple independent instances of hominin insular dwarfism occurred on Flores.

KIR3DL1 alleles and their epistatic interactions with human leukocyte antigen class I influence resistance and susceptibility to HIV-1 acquisition in the Pumwani sex worker cohort

OBJECTIVE

To determine the associations of KIR3DL1/S1(3DL1/S1) and its epistatic interactions with human leukocyte antigen class I (HLA-I) alleles with resistance and susceptibility to HIV-1.

DESIGN

Despite repeated exposure to HIV-1, a subset of women enrolled in the Pumwani sex worker cohort remain HIV uninfected. Previous studies have shown that specific HLA class I and II alleles were associated with this natural immunity. In this study, we investigated the association of 3DL1/S1 and its epistatic interactions with HLA-I, with resistance or susceptibility to HIV-1 acquisition.

METHODS

We used a sequence-based typing method to genotype 3DL1/S1 of 641 women in this cohort. The association of 3DL1/S1 and its epistatic interactions with HLA-I were analyzed using SPSS statistics software.

RESULTS

3DL1041 is enriched in the HIV-1-resistant women [P = 0.009, Pc = 0.0468, odds ratio (OR): 3.359, 95% confidence interval (CI): 1.39-8.32], whereas, 3DL1020 was associated with susceptibility to HIV-1 infection before correction for multiple comparisons (P = 0.029, Pc = 0.0858, OR: 0.316, 95%CI: 0.10-1.04). Epistatic interactions between several 3DL1 alleles and specific HLA-I alleles were observed. Among them the cocarriage of 3DL1041 with Bw4 (P = 1E - 05, Pc = 0.0015, OR: 13.33, 95%CI: 3.43-51.9), or Bw6 (P = 0.008, Pc = 0.272, OR: 3.92, 95%CI: 1.51-10.17), increased the odds of remaining HIV-1 uninfected. Further, 3DL1041+/Bw4+ women who entered the cohort HIV negative remained uninfected (P = 0.032, Pc = 0.0858). Cocarriage of 3DL101501 with C02 : 10 (P = 2.73E - 07, Pc = 7.0954E - 06), B15 : 03 (P = 3.21E - 04, Pc = 0.0042), A24 supertype (P = 8.89E - 04, Pc = 0.0077), or A23 : 01 (P = 0.0036, Pc = 0.0236) was associated with increased susceptibility to seroconversion.

CONCLUSION

The effects of interactions between 3DL1 and HLA-I alleles on resistance/susceptibility to HIV-1 infection suggest that innate immunity plays an important role in HIV-1 acquisition and should be studied and explored for HIV prevention.

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.

KIR-like activating natural killer cell receptors and their association with complicated malaria in north India

Killer immunoglobulin-like receptors (KIRs) genomic regions have been suggested to influence malaria pathogenesis and infection susceptibility. KIRs are known as activating natural killer (NK) cell receptors, which upon binding to their corresponding human leukocyte antigen (HLA) ligands cause lysis of any infected cell. We have examined the potential association of KIR genes with complicated malaria (CM) among north Indians in this study and further evaluated the KIR receptor-HLA ligand association on the severity of the disease considering the uncomplicated malaria (UCM) subjects as control. Molecular profiling of KIR and HLA was carried out using the PCR-SSP method. Susceptible association was found for individuals possessing KIR2DS2 (OR=1.76, p-value=0.0390), KIR2DL1 (OR=2.87, p-value=0.0005) and KIR2DL3 (OR=2.74, p-value=0.0011) genes with CM. This was supported by the strong linkage disequilibrium observed for 2DS2-2DL2 (D́=0.87, r²=0.54) with CM. Whereas the receptor-ligand association has revealed risk association against KIR2DS2-HLAC1 (OR=2.08, p-value=0.0229), KIR2DL3-HLAC1 (OR=1.79, p-value=0.0301), and KIR2DL1-HLAC2 (OR=2.10, p-value=0.0175) combinations for complicated malaria. The frequency of different KIR genes are more or less similar to that observed in African population showing not much genetic diversity at KIR level in context to malarial infection. In conclusion, our data indicates KIR gene loci differentially influenced the malarial outcome in north Indians and in particular the KIR2DS2 gene appeared to be associated with disease severity.