Interactions between killer cell immunoglobulin-like receptors (KIRs) and human leukocyte antigen (HLA) class I ligands regulate the development and response of human natural killer (NK) cells. Natural selection drove an allele-level group A KIR haplotype and the HLA-C1 ligand to unusually high frequency in the Japanese, who provide a particularly informative population for investigating the mechanisms by which KIR and HLA polymorphism influence NK cell repertoire and function. HLA class I ligands increase the frequencies of NK cells expressing cognate KIR, an effect modified by gene dose, KIR polymorphism, and the presence of other cognate ligand–receptor pairs. The five common Japanese KIR3DLI allotypes have distinguishable inhibitory capacity, frequency of cellular expression, and level of cell surface expression as measured by antibody binding. Although KIR haplotypes encoding 3DL1*001 or 3DL1*005, the strongest inhibitors, have no activating KIR, the dominant haplotype encodes a moderate inhibitor, 3DL1*01502, plus functional forms of the activating receptors 2DL4 and 2DS4. In the population, certain combinations of KIR and HLA class I ligand are overrepresented or underrepresented in women, but not men, and thus influence female fitness and survival. These findings show how KIR–HLA interactions shape the genetic and phenotypic KIR repertoires for both individual humans and the population.
Human killer cell immunoglobulin-like receptors (KIRs) and murine lectin-like Ly49 receptors perform orthologous functions despite their striking difference in protein structure and evolutionary origin (1, 2). As receptors for polymorphic MHC class I, KIR and Ly49 regulate NK cell tolerance and response (3). Independently, the human KIR and murine Ly49 gene families acquired some remarkably similar features, including variegated expression (4–6), signaling pathways (7), and haplotypes varying in gene content, allelic polymorphism, and signaling potential (8). Two groups of human KIR haplotypes are defined (9, 10). Group A haplotypes have a fixed organization of seven genes, most with inhibitory potential; the group B haplotypes have a variable number of KIR genes, many with activating function. Clinical studies have correlated KIR gene content with infection, cancer, autoimmunity, pregnancy syndromes, and transplant outcome (11–16). Most of these associations are with activating KIR genes and group B haplotypes and favor success in reproduction and fighting infection, but also in increased autoimmunity (3).
Although the inhibitory receptors for HLA-A (KIR3DL2), -B (3DL1), -C (2DL1, 2 and 3), and the activating receptor for HLA-G (2DL4) are well studied (17–20), the functional effects of their polymorphism remain largely unexplored. Ligands for the inhibitory 2DL5 and the activating 2DS1, 2, 3, 4, 5, and 3DS1 are unknown or uncertain. The NK cell repertoire depends on both KIR and HLA polymorphisms. Thus, KIR and HLA identical siblings have similar NK cell repertoires, whereas siblings differing at KIR or HLA exhibit a range of phenotype difference (21). Uncertain from that analysis was whether HLA serves to make positive or negative selection. Complicating further investigation of KIR–HLA interactions has been the extent of their diversity. Population comparisons revealed ethnic differences in the ratio of group A and B KIR haplotypes. In East Asian populations, the simpler group A haplotypes dominate, suggesting their KIR diversity might be more simple and experimentally more tractable than for Caucasians (22–24). Here, in-depth analysis of KIR genotype, phenotype, and function in the Japanese shows this hypothesis to be true.
KIR gene content was determined for 132 unrelated Japanese donors. 11 KIR genotypes were distinguished (Fig. 1). With a frequency of 80%, the group A haplotypes dominated. These results extend the preliminary study of 41 donors (22), giving confidence that the results obtained here well represent the KIR system of the Japanese population. Full-length cDNA from KIR expressed in peripheral blood was cloned and sequenced from 53 panel members. In selecting donors and the KIR to be sequenced from each donor, we sampled the full range of the panel's genotypes (Fig. 1) and phenotypes (Fig. 2); 18–25% of donors positive for each KIR were analyzed by cDNA cloning and sequencing, predicting detection of all alleles with >6% frequency.
Sequencing >1,100 clones identified 41 distinct KIR, of which 14 were novel (Fig. 1). Polymorphism concentrated on the group A haplotype genes, particularly those in the telomeric part of the locus. Genes that distinguish group B haplotypes (2DS1, 2, 3, 5, 3DS1, 2DL2, 2DL5) were represented by only one or two variants. Methods that distinguish the KIR alleles were devised and used to type 116 panel members. Each polymorphic KIR gene has a dominant allele, with a frequency of 45–89% compared with 4–17% for the next most common allele (Fig. 1). For 2DL4 and 2DS4, nonfunctional alleles caused by frame-shifting deletions accounted for 19 and 38% of the total alleles, respectively. In contrast, the common nonfunctional 3DL1*004 allele of Caucasians (17%) is absent from the Japanese (Fig. S1).
Polymorphism modulates the level and frequency of KIR3D expression
The variation in KIR phenotype of NK cells from 104 donors was assessed by flow cytometry using four anti-KIR antibodies and correlated with genotype (Fig. 2). The five KIR3DL1 allotypes had different levels of DX9 binding. The levels ranged over an order of magnitude, with the hierarchy 3DL1*005<*007<*001<*020<*01502 (Fig. 3 A, left). Comparison of different anti-KIR3DL1 antibodies and of DX9 binding kinetics showed the allotypic differences are the result of the abundance of cell surface 3DL1, not variable affinity for DX9 (Fig. S2, and references 25–27). Also correlating with binding level was the frequency of cellular expression so that high-binding allotypes (*001, *020, and *01502) were expressed by a larger proportion of NK cells than low-binding allotypes (*005 and *007) (Fig. 3 A, right). The frequency of NK cell expression also increased with gene dose: individuals having two low or high expressing 3DL1 allotypes had mean frequencies of expression that were 1.7- and 1.9-fold greater, respectively, than individuals having one.
Polymorphism also affected the frequency and level of 3DL2 expression. Significant differences were observed between 3DL2*008 (high-binding) and 3DL2*007 and 3DL2*00902 (low-binding). Compared with other allotypes, 3DL2*008 and 3DL2*002 had a higher frequency of NK cell expression (Fig. 3 B). Like 3DL1, a gene dose effect was seen from comparison of individuals with one or two high-expressing allotypes (Fig. 3 C). Lower binding levels and smaller differences between allotypes limited 3DL2 analysis compared with 3DL1. KIR3DL2 heterozygotes did not give the readily interpreted bimodal distributions seen for 3DL1. Furthermore, there were nine 3DL2 allotypes, compared with five for 3DL1, which increased the number of heterozygous combinations and reduced the frequency of the more readily analyzed homozygotes.
For 2DL1/S1, the frequency of NK cells binding EB6 and their binding level both increased with the number of genes encoding KIR that bind the antibody (Fig. 3 C). Although no effects correlating with 2DL1 polymorphism were observed, the panel's high 2DL1*00302 frequency suggests that its high binding could have obscured the effects of low-binding allotypes in heterozygotes. Regarding the anti-2DL2/3+2DS2 antibody (DX27), panel members were either homozygous for 2DL3 or had one copy each of 2DL2, 2DL3, and 2DS2. Donors with the latter genotype had a similar frequency of antibody-binding NK cells, but a substantially increased level of binding (Fig. 3 C).
Cognate HLA class I increases the frequency of NK cells expressing inhibitory KIR
HLA-A, -B, -C types for the panel were determined at high resolution (Fig. S3 A). Both the alleles and their frequencies were consistent with previous studies (28, 29). Hardy-Weinberg equilibrium was observed for KIR and HLA (Fig. S3 B). Distinguishing Japanese from Caucasians are the frequencies of HLA-C ligands for KIR2DL. In Caucasians, 66% of HLA-C allotypes have C1 (the 2DL2/3 ligand) and 34% have C2 (the KIR2DL1 ligand); whereas in Japanese, 92% have C1 and 8% have C2. Thus, almost all Japanese use C1 as a self-ligand for inhibitory KIR, whereas a small minority use C2. The frequency of HLA-A3/11 ligands for 3DL2 is also lower in Japanese (8%) than in Caucasians (22%), whereas they have identical frequency for HLA-Bw4, the 3DL1 ligand (36%).
KIR phenotypes were compared for donors who have or lack a cognate HLA class I ligand for particular KIR. In several combinations, presence of cognate ligand increased the frequency of NK cells expressing the KIR (Fig. 4 A), while decreasing the amount of KIR-specific antibody that is bound (Fig. 4 B). In 2DL1 homozygotes, the presence of cognate C2 almost doubled the frequency of 2DL1+ NK cells while reducing the level of anti-2DL1 binding by 26%. Although the panel's high C1 frequency precluded comparable analysis of 2DL2/3, the two donors who lacked C1 showed decreased frequency of 2DL2/3-expressing cells and reduced levels of antibody binding (Fig. S4). C1 homozygotes had decreased binding of anti-2DL2/3 to NK cells compared with C1/C2 heterozygotes. Presence of Bw4 significantly increased the frequency of cells expressing 3DL1 for donors having two high-binding 3DL1 allotypes and decreased the level of anti-KIR3DL1 bound for donors having one high-binding 3DL1 allotype (Fig. 4, A and B). For donors having one high-binding 3DL1 allotype, paired with a low binding allotype or 3DS1, there was little effect on expression frequency, as also seen for donors having two low-binding allotypes (Fig. 4 C, top).
Cognate HLA class I decreases the frequency of NK cells expressing other inhibitory KIR
For donors with two high-binding 3DL1 allotypes, the increase in frequency of 3DL1+ NK cells is greater in Bw4 heterozygotes than in homozygotes (Fig. 4 C, top). The effect of cognate Bw4 to increase the frequency of 3DL1+ cells was also depressed by the presence of additional cognate HLA–KIR interactions, and the effects were additive. Thus, for individuals who had all four inhibitory ligands (C1, C2, Bw4, and A3/11), the frequencies of 3DL1+ cells are comparable to those in individuals lacking Bw4 (Fig. 4 C, top). Although similar trends were seen for individuals having one high-binding 3DL1 allotype, they did not reach statistical significance. A similar effect of other cognate HLA–KIR interactions was also seen from analysis of the frequency of NK cells expressing 2DL2/3 (Fig. 4 C, bottom).
Allelic polymorphism modulates the inhibitory function of KIR3DL1
Approximately one third of the Japanese population has C1-2DL3 as the only functional inhibitory KIR–HLA class I interaction, while one half of the population can also use the Bw4–3DL1 interaction. Although 2DL3 has little polymorphism in the Japanese, 3DL1 variation is substantial (Fig. 1). Consequently, the Bw4–3DL1 interaction is a potentially important source of functional NK cell diversification in the Japanese population. To address this question, the functions of the five 3DL1 allotypes present in the Japanese panel were compared.
PBMCs from 13 donors, each having a high- and low-binding 3DL1 allotype, were challenged with class I–deficient 721.221 cells and with 221 cells transfected with HLA-B*5801, a 3DL1 ligand. In response to 221 cells, 3DL1+ NK cells produced IFN-γ as detected by intracellular staining. A diminished response was made to 221- B*5801, the inhibition varying 40–90% between the five 3DL1 allotypes (Fig. 5, top). The hierarchy for increasing inhibition was 3DL1*007<*020<*01502<*005<*001. That 3DL1 polymorphism modulates a major NK cell effector function demonstrates its functional importance. Of note, the two allotypes giving strongest inhibition are those expressed at intermediate (3DL1*001) and lowest levels (3DL1*005).
Further modulating each 3DL1 allotype's capacity to influence the NK cell response is its expression frequency in the NK cell population. Consequently, the impact of an allotype such as 3DL1*005 could be restricted to a relatively small number of cells. By taking both factors (the frequencies of expression and inhibition) into account, we estimated the frequency of NK cells in peripheral blood that are inhibitable by B*5801 (Fig. 5, bottom). The difference between this hierarchy (3DL1*007<*005<*020<*01502<*001) and that in Fig. 5′s top half is the position of 3DL1*005, which shifts from the second strongest to second weakest allotype.
Positive and balancing selection of KIR
70 allele-level KIR genotypes were distinguished in the panel. From these genotypes and the sequences of prototypical group A (GenBank accession no. AC011501) and group B (GenBank accession no. AY320039) haplotypes, the structures of Japanese KIR haplotypes were inferred by two independent methods. The 20 most frequent haplotypes were identically ranked (Fig. 6 A). Each was identified in two or more donors and, collectively, they account for 83.4% of the panel's haplotypes. In both gene content and allelic polymorphism, the telomeric part of Japanese KIR haplotypes is more variable than the centromeric part. Together, the 12 deduced A haplotypes comprise 67.7% of total haplotypes compared with 15.7% for the 8 deduced B haplotypes. One A haplotype is markedly more common (35.7%) than any other haplotype (<9%) and, for each of its component KIR genes, this haplotype has the dominant allele.
Application of the Ewens-Watterson test (30, 31) to the polymorphic KIR genes of the group A haplotypes shows that centromeric genes specifying HLA-C receptors (2DL1 and 3) have been subject to positive directional selection (Fig. 6 B), consistent with the high frequencies of the dominant alleles at these genes (Fig. 1). In contrast, the telomeric 2DL4, 3DL1, 2DS4, and 3DL2 genes, including those encoding HLA-A, -B, and -G receptors, appear subject to balancing selection, consistent with their wider range of allele frequencies. Tajima's D statistic (32) was also used to assess the mode and strength of selective pressure on polymorphic sites in group A haplotypes, using composite 7-kb haplotypes constructed from the coding sequences of KIR alleles in their order on the haplotype (Fig. 6 C, within the panel). Aside from confirming directional selection on the 2DL1 and 3 genes, this analysis revealed a strong signature for balancing selection in the telomeric region, one focusing on the ligand-binding domains of 2DL4, 3DL1, and 2DS4 and the signaling domain of 2DL4.
Linkage disequilibrium (LD) can indicate positive selection (33). Haplotype homozygosity revealed relatively high LD throughout the Japanese group A haplotypes; in particular, from 2DL4 through 3DL1 to 2DS4 (Fig. S5). Within this core region, the five 3DL1 alleles are linked with different combinations of alleles for the other two genes. These various cores are further diversified by combination with different 3DL2 alleles, the number depending on the core (Fig. 6 D). For example, 3DL1*005 pairs with four 3DL2 alleles, whereas the haplotype containing the dominant 3DL1*01502 allele is almost homogeneous, which is strong evidence for its positive selection, a selective sweep. As an independent assessment, the Ewens-Watterson homozygosity statistic was calculated for the centromeric and telomeric parts of the group A haplotypes and tested for departure from neutrality. In both parts of the KIR locus, signatures of positive selection were detected (Fig. 6 C, below the panel).
These analyses show that positive selection drove one group A haplotype to high frequency, such that its component alleles are now the dominant alleles of the polymorphic KIR genes. Positive selection also homogenized the centromeric part of the KIR locus. In contrast, balancing selection on telomeric KIR has maintained a variety of other haplotypes and their component alleles at intermediate frequencies. Because of these contrasting selections, a majority of Japanese (∼75%) have the same centromeric KIR genes, including the HLA-C–specific KIR2DL, but differ in the telomeric genes, including the HLA-A– and -B–specific KIR3D.
In addition to the 20 haplotypes shown in Fig. 6 A, 24 rarer haplotypes were deduced by assuming patterns of LD defined previously in other populations (34). Of the 44 haplotypes, 24 were group A and 20 group B. 30 of the 44 haplotypes share some alleles with the most common haplotype (no. 1) (Fig. S6). As well as a high frequency group A haplotype, the Japanese population has an abundance of low frequency A and B haplotypes with centromeric and telomeric motifs distinct from those present in the common haplotypes. These provide a pool of diverse functional potential for natural selection to operate on in the future.
Non-random associations of KIR–HLA in females
Systematic analysis revealed combinations of KIR and HLA genes in the donor panel at non-random frequencies. Stratification by sex showed that the effect was confined to females (Fig. 7). Deleted, nonfunctional forms of 2DL4 (35) preferentially associated with C2, whereas homozygosity for functional, full-length 2DL4 correlated negatively with C2. 3DL1*005 preferentially associated with Bw4, a combination also favoring inhibitory interaction. Although C2 and Bw4 are in strong LD in Caucasians, that is not so for Japanese, for whom only 40% of C2-carrying haplotypes also have HLA-Bw4 (36). Also, only 50% of donors carrying C2 have Bw4 (8 out of 16), a neutral distribution confirming the independence of C2 and Bw4 in the panel. Association with Bw4 is thus an independently associated factor and not a consequence of LD with C2. Weaker trends were HLA-C2 with 3DL1*005 or *001, the most potent inhibitors (Fig. 5), whereas 3DL1*01502, the dominant allotype and less potent inhibitor, correlated negatively with C2.
Six core haplotypes are defined by the five 3DL1 alleles and 3DS1 (Fig. 8). The cores defined by 3DL1*005 and *001 both encode deleted, nonfunctional forms of 2DL4 and 2DS4. Group A haplotypes with these cores lack activating KIR and have a strong inhibitory Bw4 receptor as well as inhibitory receptors for C1 and C2. Consequently, the overrepresented combination of 3DL1*005, nonfunctional 2DL4, Bw4 homozygosity, and C2 (usually accompanied by C1) is one that suggests some selection for strong inhibitory KIR–HLA interactions. Conversely, the underrepresented combinations of functional forms of 2DL4 with C2 suggest selection for less inhibition and more activation.
Selection of NK cell KIR repertoire by cognate HLA class I ligands
In mouse models, acquisition of inhibitory Ly49 receptors by NK cells is informed by interactions with cognate MHC class I ligands expressed on bone marrow stromal cells (37). The common observation is that cognate ligand, usually a transgene product, reduces the frequency of receptor-expressing cells. For humans, we find both increase and decrease in the frequency of receptor expression, which depends on the number of cognate KIR–HLA pairs contributing to NK cell self-tolerance. When acting alone, the cognate ligand increases the frequency of the KIR's expression, but the effect varies with receptor and ligand gene dose. Such increase is consistent with previous observation that NK cell self-tolerance is mediated by inhibitory HLA class I receptors (6). In the presence of other ligand–receptor pairs, the expression frequency for each receptor is reduced from that achieved when it acts alone, and to an extent inversely correlated with the number of functional ligand–receptor pairs.
For mice, it has been proposed that sequential or stochastic expression of inhibitory self-reactive receptors occurs until cumulative inhibition reaches some threshold, whereupon further receptor acquisition stops (38–40). Our results fit with this type of model. That the cognate ligand has not been observed to elevate receptor expression frequency in mice could reflect a difference in the mechanism between the two species or in the mode of analysis. For example, the presence in mice of multiple Ly49 products with overlapping ligand specificity could have obscured the positive effect, as we observed for humans with multiple KIR–HLA pairs.
KIR polymorphism modulates NK cell effector function
For ∼33% of the Japanese population, C1 is the only self-ligand for inhibitory KIR. Because the Japanese are practically monomorphic for KIR2D and polymorphic for KIR3DL1, they provided a good context in which to investigate the function of 3DL1 polymorphism. The five 3DL1 allotypes present in the Japanese differ in their capacity to inhibit NK cell function; they also differ in repertoire, which is the frequency and level of cell surface expression. Much evidence shows that DX9 binding measures the amount of cell surface KIR3DL1 (Fig. S10 and references 25–27), so the allotypic differences observed are not the result of varying affinity for DX9.
The allotypic differences in the frequency and level of KIR3DL1 expression are distinguished from those caused by cognate ligand and represent a ligand-independent influence of 3DL1 polymorphism on the primary NK cell repertoire. The effects of cognate 3DL1–Bw4 interactions are also modulated by KIR3DL1 polymorphism. Notably, the increase in frequency induced by the cognate ligand is greater for high than for low-expressing KIR3DL1 allotypes. The overall trend is for inhibitory capacity to increase with the level and frequency of 3DL1 expression, an exception being 3DL1*005, which combines low frequency and level of expression with strong inhibitory function.
KIR diversity is subject to positive and balancing selection
An outstanding question is whether the haplotypic diversity and allelic polymorphism of KIR genes arose from functional selection or neutral change. From functional assay and statistical analysis, we provide evidence that balancing selection maintains a pool of alleles and haplotypes with distinguishable functions upon which positive selection can act. In the Japanese, one group A KIR haplotype has risen to unusually high frequency, a selection that may have been targeted to the centromeric KIR genes, because the 2DL1*00302 and 2DL3*001 alleles of the selected haplotype are at frequencies approaching 90%. In contrast, in the telomeric parts of the locus, the balancing selection retained diversity as exemplified by the even frequency of other group A and group B haplotype motifs and the distribution of the six core group A haplotype structures (Fig. 8). In these regions, the selected haplotype encodes functional KIR2DL4 and KIR2DS4 forms and a KIR3DL1 allotype giving moderate inhibition. Within the set of Japanese A haplotypes, it is one that combines a full set of functional genes and good potential for both inhibition and activation. Because of this haplotype, the frequency of functional 2DL4 and 2DS4 alleles in the Japanese is much higher than in Caucasians (35, 41). One possible cause of the selective sweep is an epidemic of infectious disease where the selected haplotype enhanced the chance of surviving infection. Alternatively, women carrying this less inhibitory, more activating form of A haplotype could have above average reproductive success, as might be suggested from the results of Hiby et al. (15).
Evidence for selection of KIR–HLA combinations in female genotypes
An unexpected and thought-provoking discovery is the over- and underrepresentation of certain KIR–HLA class I combinations in Japanese females. Because KIR and HLA are on different chromosomes, the overrepresented combinations must confer some selective advantage, whereas the underrepresented combinations are disadvantageous. As only NK cells and T cells express KIR, the selective mechanism likely involves one or both of these cell types. Candidate mechanisms can be distinguished according to when they might act: before conception, during fetal life, or after birth.
Acting before conception and biasing progeny genotype is mate choice. Extensive literature points to MHC polymorphism influencing the choice of sexual partners (42–44), and it is plausible that this biology could extend to combinations of KIR and HLA polymorphisms. For example, preferential choice of Bw4+ partners by Bw4+, 3DL1*005+ individuals would lead to enrichment in the next generations of genotypes combining Bw4 homozygosity with 3DL1*005. For mate choice, a cause for the female bias is hard to envisage.
As only a fraction of human embryos survive implantation and fetal life (45), there is potential for improving survival. There are also sex differences in the mechanisms by which embryos fail to survive and the timing of these events (46). For example, for embryos with normal karyotype, the rate of spontaneous abortion during early pregnancy is greater for females than males, so any combination of genetic factors that reduced the effect would be specifically overrepresented in the female population (46). Candidate mechanisms involve KIR–HLA interactions within the growing fetus, which could contribute to tissue development or defense against infection. NK cells are detectable in fetuses aborted at 6–24 wk of gestation, pointing to their possible contribution (47). Also of possible relevance, KIR3DL3 transcripts are more abundant in the NK cells of females than males and, although the protein appears undetectable in adult cells, it may function during fetal development (48).
After birth, the favored KIR–HLA combinations could improve survival, relative to other combinations, through NK or T cell responses that eliminate infection or cancer, or prevent autoimmune disease. Female specificity could arise in a general way from the hormonal differences between men and women, and their influences on the immune system (49, 50), or more specifically to immune functions under circumstances experienced only by women, such as menstruation, pregnancy, nursing, cancer of breast, or ovary.
Precedent for KIR–HLA combinations influencing human reproduction and survival is the work of Hiby et al. (15). They found that inhibitory combinations of fetal HLA-C2 and maternal homozygosity for A KIR haplotypes predispose Caucasians to preeclampsia. These correlations fit well with the observed interactions of fetal trophoblast with maternal NK cells during placentation. Selection by preeclampsia can also explain the inverse correlation between C2 and A haplotype frequencies in human populations, as exemplified by the low C2 and high group A frequencies in the Japanese. In this context, it is important to note that the correlations reported here do not address the fitness of allogeneic combinations of maternal and fetal factors, but of the KIR and HLA genes brought together in the child's genotype. Consequently, the biology underlying the overrepresented and underrepresented KIR–HLA combinations we find may well be distinct from those described by Hiby et al. (15).
A signature of natural selection on KIR diversity
KIR-expressing NK cells participate in defense against infection and in reproduction, two functions that are essential for the survival of individuals, populations, and species. Both the KIR and their cognate HLA class I ligands are highly variable. In comparison with Caucasians, the Japanese population has a skewed distribution of KIR and HLA variants, reflecting a distinct history of directional and balancing selection. Functional and genetic study of a Japanese panel provides evidence here that the combination of KIR and HLA variation acts to diversify NK cell repertoire and effector function within the population and to facilitate human survival.
Materials And Methods
Human subjects and blood samples.
Blood samples were from 132 healthy, unrelated Japanese individuals (54 males and 78 females, ages 22–73, mean 40 ± 12 yr [SD]) and from related donors in three pedigrees; informed consent was as approved by the Stanford University Institutional Review Board on Human Subjects. PBMCs were prepared with Ficoll-Paque PLUS (GE Healthcare). Genomic DNA and RNA were extracted with DNAzol (Molecular Research Center) and TRIzol (Invitrogen). Complementary DNA was synthesized from RNA, using oligo(dT) and SMART II oligonucleotide primers (SMART-RACE kit; CLONTECH Laboratories, Inc.).
KIR and HLA typing.
KIR genes were typed as described previously (22). 5′RACE used the Advantage 2 PCR kit (CLONTECH Laboratories, Inc.), with the universal primer (UPM) from the SMART-RACE kit or KIR-specific primers as the forward primer combined with KIR-specific reverse primers (Fig. S7). KIR transcripts were cloned into the pGEM T-easy vector (Promega) and sequenced using BigDye terminator 3.1 chemistry and a PRISM 377 DNA sequencer (Applied Biosystems). Sequence assembly and analysis used AlignIR version 2.0 (LICOR). Three or more clones from two or more independent PCR reactions were sequenced to define novel polymorphisms. KIR allele typing by sequence-specific polymorphism PCR used novel primers and those described previously (25, 34, 41, 51) (Figs. S8–S10). PCR used hot-start procedures using AmpliTaq FS on a PE9600 thermalcycler (Applied Biosystems). The juxtaposition of 2DL5, 2DS3, and 2DS5 to either 2DL2 or 3DS1 was determined by long-range PCR (52). Exon 1 was sequenced directly from the PCR products to distinguish 2DL5*002 from *005. In exon 6 of 2DL4, the cluster of 9 or 10 adenines (9A/10A) were typed by sequencing (35). HLA-A, B, C typing was by sequencing locus-specific PCR products (minimally exons 2 and 3, but also exon 4 for some samples) with adaptations from Dunn et al. (53) and PCR-SSOP (Sequence-Specific Oligonucleotide Probing), using LABType SSO Typing Tests (One Lambda).
Flow cytometry and NK cell secretion of IFN-γ.
KIR on CD56+CD3−NK cell surfaces was measured with monoclonal antibodies and flow cytometry (22). Phycoerythrin conjugates of anti-2DL1/2DS1 (EB6; Beckman Coulter), anti-2DL2/2DL3/2DS2 (DX27), and anti-3DL1(DX9) (BD Biosciences), anti-3DL2 (DX31, a gift from Drs. J. Phillips and L. Lanier, DNAX Research Institute, Palo Alto, CA) were used in combination with anti-CD3–PerCP (SK7) and anti-CD56–fluorescein (NCAM16.2) (BD Biosciences) to stain PBMCs. Three-color flow cytometry used a FACScan instrument and CELLQuest software (BD Biosciences). Parameters measured were the percent NK cells binding an antibody and the mean fluorescence intensity (mfi) for the positive cells.
PBMCs from donors heterozygous for high- and low-binding 3DL1 allotypes were cocultured with the class I–deficient 721.221 cell line or 221 transfected with HLA-B*5801. After gating on 3DL1+ NK cells, IFN-γ production was assessed by intracellular staining (54): PBMCs at 6 × 105 cells/well were cocultured with target cells at an optimized effector–target ratio of 1:1 for 9 h; fluorescent anti-CD3 and anti-CD85–PECy5 (GHI/75) (BD Biosciences) were added so that 3DL1+ T cells and the inhibitory effects of LIR1 could be excluded from analysis. Knowing the proportion of donor NK cells expressing each 3DL1 allotype, the inhibitory effect on total NK cells was calculated for each allotype: (% IFN-γ inhibition) × (% allotype-expressing NK cells).
Reconstruction of allele-level KIR haplotypes used the expectation maximization algorithm (Arlequin version 2.001; reference 55) and a Bayesian statistical method using hidden Monte Carlo Markov Chain modeling (PHASE version 2.1; reference 56). LD was calculated for pairs of KIR alleles (Arlequin version 2.001), D′ values >0.5 with P < 0.002 are shown in Fig. S5. The homozygosity statistic F of Ewens (30) and Watterson (31), calculated from KIR allele frequencies was compared with expected homozygosity statistics; significance levels were tested using Slatkin's exact test (57) (Arlequin version 2.001). Tajima's D (32) was calculated for 7-kb haplotypes constructed by combining the coding sequences of KIR alleles on group A haplotypes in order. 10 potential recombination points between haplotypes were first detected by the four-gamete test (58). D values were then calculated in a 300-bp sliding window along the composite haplotypes without extending over the recombination points, and visualized (DnaSP version 4.0; reference 59). For the 9A KIR2DL4 alleles and the 22-bp deletion 2DS4 alleles, each deleted nucleotide was considered an independent position. Haplotype homozygosity was calculated at each gene in group A haplotypes based on Sabeti et al.'s method (33).
Nonparametric tests of significance were used in the analyses described in Figs. 3, 4, 5 and Fig. S4 (two-tailed Mann-Whitney's U-test). Chi-squared tests were used in 2 × 2 contingency tables for pairwise comparisons of genetic associations in Fig. 7. P < 0.05 was considered statistically significant. Fisher's exact test was applied when expected numbers were small (n < 5). Tests used SPSS version 8 software (SPSS Inc). Significance values are as follows: *, P < 0.05; **, P < 0.01; ***, P < 0.005; ****, P < 0.001; and *****, P < 0.0005.
Online supplemental material.
Fig. S1 depicts comparison of KIR allele frequencies between Japanese and Caucasians. Fig. S2 shows kinetics of DX9 binding and phenotype comparison with Z27. Fig. S3 A contains comparison of HLA class I allele distributions between Japanese and Caucasians. Fig. S3 B shows heterozygosity of HLA class I and KIR. Fig. S4 depicts the influence of cognate HLA class I on KIR expression. Fig. S5 displays LD between KIR alleles. Fig. S6 shows the frequency of centromeric and telomeric parts in 44 haplotypes. Fig. S7 contains primers and PCR conditions for KIR cDNA cloning. Fig. S8 depicts primers for KIR allele typing. Fig. S9 shows PCR conditions for KIR allele typing. Fig. S10 depicts primer combinations for KIR allele typing. Online supplemental figures are available.
We thank all blood donors who participated in this study; we also thank M. Minami and his staff for their help in blood sample collection. DX31 was a gift from J. Phillips and L. Lanier. We thank T. Hastie for his advice on statistics.
This work was supported by National Institutes of Health (NIH) grant nos. AI-017892 and AI-022039 (to P. Parham) and a grant from the Leukemia Research Foundation (to M. Yawata).
The authors have no conflicting financial interests.
Abbreviations used: KIR, killer cell immunoglobulin-like receptor; LD, linkage disequilibrium; mfi, mean fluorescence intensity.
M. Yawata and N. Yawata contributed equally to this work.