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.

Behcet disease is a chronic, relapsing-remitting autoinflammatory syndrome with a strong HLA-B*51 association. In this paper, we describe a human cohort of 267 individuals with Behcet disease and 445 matched controls from a tertiary referral center in the U.K. HLA-B*51 was confirmed as a genetic risk factor in this group (p = 0.0006, Bonferroni-Dunn correction for multiple testing [Pc] = 0.0192, odds ratio [OR] 1.92, 95% confidence interval [CI] 1.33-2.76). KIR3DL1/S1 allele-level analysis indicated that low-expressing KIR3DL1/S1 alleles in combination with KIR3DS1 increased the risk of developing Behcet disease (KIR3DL1(Low)/KIR3DS1: p = 0.0004, Pc = 0.0040, OR 2.47, 95% CI 1.43-4.25), whereas high-expressing KIR3DL1/S1 alleles in combination with a null-expressing KIR3DL1 reduced the risk of disease (KIR3DL1(HIGH)/KIR3DL(NULL): p = 0.0035, Pc = 0.0350, OR 0.53, 95% CI 0.33-0.87). Behcet disease can manifest as a purely mucocutaneous disease or can involve other organ systems such as the eyes. In the U.K. cohort studied in this study, KIR3DL1(Low/)KIR3DS1 increased the risk of ophthalmic disease (p = 1.2 x 10(-5), OR 3.92, 95% CI 2.06-7.47), whereas KIR3DL1(HIGH)/KIR3DL1(NULL) reduced the risk of having purely mucocutaneous disease (p = 0.0048, OR 0.45, 95% CI 0.25-0.81). To our knowledge, this is the first analysis of KIR3DL1/S1 allelic variation in Behcet disease and may provide insight into the pathogenic role of HLA-B*51 and its interaction with KIR3DL1/S1.

Natural killer (NK) cells can recognize virus-infected and stressed cells(1) using activating and inhibitory receptors, many of which interact with HLA class I. Although early studies also suggested a functional impact of HLA class II on NK cell activity(2,3), the NK cell receptors that specifically recognize HLA class II molecules have never been identified. We investigated whether two major families of NK cell receptors, killer-cell immunoglobulin-like receptors (KIRs) and natural cytotoxicity receptors (NCRs), contained receptors that bound to HLA class II, and identified a direct interaction between the NK cell receptor NKp44 and a subset of HLA-DP molecules, including HLA-DP401, one of the most frequent class II allotypes in white populations(4). Using NKp444 zeta(+) reporter cells and primary human NKp44(+) NK cells, we demonstrated that interactions between NKp44 and HLA-DP401 trigger functional NK cell responses. This interaction between a subset of HLA-DP molecules and NKp44 implicates HLA class II as a component of the innate immune response, much like HLA class I. It also provides a potential mechanism for the described associations between HLA-DP subtypes and several disease outcomes, including hepatitis B virus infection(5-7), graft-versus-host diseases and inflammatory bowel disease(9,10).

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(+) T*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.

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.

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/31, 13 of KIR3DL2 and 18 of KIR3DL3. Previously identified alleles, corresponding to 33 alleles of KIR3DL1/31, 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 similar to KIR3DL1/S1 similar to 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.

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 not equal 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 similar to 200 haplotypes for each gene family. All 74 Nama and 78 not equal 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 not equal 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 not equal Khomani San appears to have led to natural selection against strong inhibitory C2-specific KIR.