Natural Killer (NK) cells play critical roles in immune defense and reproduction, yet remain the most poorly understood major lymphocyte population. Because their activation is controlled by a variety of combinatorially expressed activating and inhibitory receptors, NK cell diversity and function are closely linked. To provide an unprecedented understanding of NK cell repertoire diversity, we used mass cytometry to simultaneously analyze 35 parameters, including 28 NK cell receptors, on peripheral blood NK cells from five sets of monozygotic twins and twelve unrelated donors of defined HLA and killer cell immunoglobulin-like receptor (KIR) genotype. This analysis revealed a remarkable degree of NK cell diversity, with an estimated 6,000-30,000 phenotypic populations within an individual and >100,000 phenotypes in this population. Genetics largely determined inhibitory receptor expression, whereas activation receptor expression was heavily environmentally influenced. Therefore, NK cells may maintain self-tolerance through strictly regulated expression of inhibitory receptors, while using adaptable expression patterns of activating and costimulatory receptors to respond to pathogens and tumors. These findings further suggest the possibility that discrete NK cell subpopulations could be harnessed for immunotherapeutic strategies in the settings of infection, reproduction, and transplantation.

Summary Natural killer (NK) cells are lymphocytes having vital functions in innate and adaptive immunity, as well as placental reproduction. Controlling education and functional activity of human NK cells are various receptors that recognize HLA class I on the surface of tissue cells. Epitopes of polymorphic HLA-A,-B and –C are recognized by equally diverse killer cell immunoglobulin-like receptors (KIR). In addition, a peptide cleaved from the leader sequence of HLA-A,-B or –C must bind to HLA-E for it to become a ligand for the conserved CD94:NKG2A receptor. Methionine/threonine dimorphism at position -21 of the leader sequence divides HLA-B allotypes into a majority having -21T that do not supply HLA-E binding peptides and a minority having -21M, that do. Genetic analysis of human populations worldwide shows how haplotypes with -21M HLA-B rarely encode the KIR ligands: Bw4+HLA-B and C2+HLA-C KIR. Thus there are two fundamental forms of HLA haplotype: one preferentially supplying CD94:NKG2A ligands, the other preferentially supplying KIR ligands. -21 HLA-B dimorphism divides the human population into three groups: M/M, M/T and T/T. Mass cytometry and assays of immune function, shows how M/M and M/T individuals have CD94:NKG2A+ NK cells which are better educated, phenotypically more diverse and functionally more potent than those in T/T individuals. Fundamental new insights are given to genetic control of NK cell immunity and the evolution that has limited the number of NK cell receptor ligands encoded by an HLA haplotype. These finding suggest new ways to dissect the numerous clinical associations with HLA class I.

The physiological functions of natural killer (NK) cells in human immunity and reproduction depend upon diverse interactions between killer cell immunoglobulin-like receptors (KIRs) and their HLA class I ligands: HLA-A, HLA-B, and HLA-C. The genomic regions containing the KIR and HLA class I genes are unlinked, structurally complex, and highly polymorphic. They are also strongly associated with a wide spectrum of diseases, including infections, autoimmune disorders, cancers, and pregnancy disorders, as well as the efficacy of transplantation and other immunotherapies. To facilitate study of these extraordinary genes, we developed a method that captures, sequences, and analyzes the 13 KIR genes and HLA-A, HLA-B, and HLA-C from genomic DNA. We also devised a bioinformatics pipeline that attributes sequencing reads to specific KIR genes, determines copy number by read depth, and calls high-resolution genotypes for each KIR gene. We validated this method by using DNA from well-characterized cell lines, comparing it to established methods of HLA and KIR genotyping, and determining KIR genotypes from 1000 Genomes sequence data. This identified 116 previously uncharacterized KIR alleles, which were all demonstrated to be authentic by sequencing from source DNA via standard methods. Analysis of just two KIR genes showed that 22% of the 1000 Genomes individuals have a previously uncharacterized allele or a structural variant. The method we describe is suited to the large-scale analyses that are needed for characterizing human populations and defining the precise HLA and KIR factors associated with disease. The methods are applicable to other highly polymorphic genes.

Interactions between HLA class I molecules and killer-cell immunoglobulin-like receptors (KIR) control natural killer cell (NK) functions in immunity and reproduction. Encoded by genes on different chromosomes, these polymorphic ligands and receptors correlate highly with disease resistance and susceptibility. Although studied at low-resolution in many populations, high-resolution analysis of combinatorial diversity of HLA class I and KIR is limited to Asian and Amerindian populations with low genetic diversity. At the other end of the spectrum is the West African population investigated here: we studied 235 individuals, including 104 mother-child pairs, from the Ga-Adangbe of Ghana. This population has a rich diversity of 175 KIR variants forming 208 KIR haplotypes, and 81 HLA-A, -B and -C variants forming 190 HLA class I haplotypes. Each individual we studied has a unique compound genotype of HLA class I and KIR, forming 1–14 functional ligand-receptor interactions. Maintaining this exceptionally high polymorphism is balancing selection. The centromeric region of the KIR locus, encoding HLA-C receptors, is highly diverse whereas the telomeric region encoding Bw4-specific KIR3DL1, lacks diversity in Africans. Present in the Ga-Adangbe are high frequencies of Bw4-bearing HLA-B*53:01 and Bw4-lacking HLA-B*35:01, which otherwise are identical. Balancing selection at key residues maintains numerous HLA-B allotypes having and lacking Bw4, and also those of stronger and weaker interaction with LILRB1, a KIR-related receptor. Correspondingly, there is a balance at key residues of KIR3DL1 that modulate its level of cell-surface expression. Thus, capacity to interact with NK cells synergizes with peptide binding diversity to drive HLA-B allele frequency distribution. These features of KIR and HLA are consistent with ongoing co-evolution and selection imposed by a pathogen endemic to West Africa. Because of the prevalence of malaria in the Ga-Adangbe and previous associations of cerebral malaria with HLA-B*53:01 and KIR, Plasmodium falciparum is a candidate pathogen.

The human killer cell immunoglobulin-like receptor (KIR) locus comprises two groups of KIR haplotypes, termed A and B. These are present in all human populations but with different relative frequencies, suggesting they have different functional properties that underlie their balancing selection. We studied the genomic organization and functional properties of the alleles of the inhibitory and activating HLA-C receptors encoded by KIR haplotypes. Because every HLA-C allotype functions as a ligand for KIR, the interactions between KIR and HLA-C dominate the HLA class I mediated regulation of human NK cells. The C2 epitope is recognized by inhibitory KIR2DL1 and activating KIR2DS1, whereas the C1 epitope is recognized by inhibitory KIR2DL2 and KIR2DL3. This study shows that the KIR2DL1 and 2DS1 and KIR2DL2/3 alleles form distinctive phylogenetic clades that associate with specific KIR haplotypes. KIR A haplotypes are characterized by KIR2DL1 alleles that encode strong inhibitory C2 receptors and KIR2DL3 alleles encoding weak inhibitory C1 receptors. In striking contrast, KIR B haplotypes are characterized by KIR2DL1 alleles that encode weak inhibitory C2 receptors and KIR2DL2 alleles encoding strong inhibitory C1 receptors. The wide-ranging properties of KIR allotypes arise from substitutions throughout the KIR molecule. Such substitutions can influence cell-surface expression, as well as the avidity and specificity for HLA-C ligands. Consistent with the crucial role of inhibitory HLA-C receptors in self-recognition, and natural killer cell education and response, most KIR haplotypes have both a functional C1 and C2 receptor, despite the considerable variation that occurs in ligand recognition and surface expression.

In sub-Saharan Africans, maternal mortality is unacceptably high, with >400 deaths per 100,000 births compared with <10 deaths per 100,000 births in Europeans. One-third of the deaths are caused by pre-eclampsia, a syndrome arising from defective placentation. Controlling placentation are maternal natural killer (NK) cells that use killer-cell immunoglobulin-like receptor (KIR) to recognize the fetal HLA-C molecules on invading trophoblast. We analyzed genetic polymorphisms of maternal KIR and fetal HLA-C in 484 normal and 254 pre-eclamptic pregnancies at Mulago Hospital, Kampala, Uganda. The combination of maternal KIR AA genotypes and fetal HLA-C alleles encoding the C2 epitope associates with pre-eclampsia [P = 0.0318, odds ratio (OR) = 1.49]. The KIR genes associated with protection are located in centromeric KIR B regions that are unique to sub-Saharan African populations and contain the KIR2DS5 and KIR2DL1 genes (P = 0.0095, OR = 0.59). By contrast, telomeric KIR B genes protect Europeans against pre-eclampsia. Thus, different KIR B regions protect sub-Saharan Africans and Europeans from pre-eclampsia, whereas in both populations, the KIR AA genotype is a risk factor for the syndrome. These results emphasize the importance of undertaking genetic studies of pregnancy disorders in African populations with the potential to provide biological insights not available from studies restricted to European populations.Uganda | pre-eclampsia | NK cells | maternal mortality | KIR

The most polymorphic part of the human genome, the encodes over 160 proteins of diverse function. Half of them, including the and genes, are directly involved in immune responses. Consequently, the region strongly associates with numerous diseases and clinical therapies. Notoriously, the region has been intractable to high-throughput analysis at complete sequence resolution, and current reference haplotypes are inadequate for large-scale studies. To address these challenges, we developed a method that specifically captures and sequences the 4.8-Mbp region from genomic DNA. For 95 homozygous cell lines we assembled, de novo, a set of high-fidelity contigs and a sequence scaffold, representing a mean 98% of the target region. Included are six alternative reference sequences of the human genome that we completed and refined. Characterization of the sequence and structural diversity of the region shows the approach accurately determines the sequences of the highly polymorphic and genes and the complex structural diversity of complement factor It has also uncovered extensive and unexpected diversity in other genes; an example is, which encodes a lung mucin and exhibits more coding sequence alleles than any or gene studied here. More than 60% of the coding sequence alleles analyzed were previously uncharacterized. We have created a substantial database of robust reference haplotype sequences that will enable future population scale studies of this complicated and clinically important region of the human genome.

Djaoud et al. show that Epstein–Barr virus infection triggers two types of human innate immune response, one mediated by the combination of NK cells and γδ T cells and the other committed to a strong NK cell response with little involvement of γδ T cells.