Abstract
Interactions of human natural killer (NK) cell inhibitory receptors with polymorphic HLA-A, -B and -C molecules educate NK cells for immune surveillance against tumor cells. The KIR A haplotype encodes a distinctive set of HLA-specific NK cell inhibiting receptors having strong influence on immunity. We observed higher frequency of KIR A homozygosity among 745 healthy Chinese Southern Han than 836 adult patients representing three types of leukemia: ALL (OR = 0.68, 95% CI = 0.52–0.89, p = 0.004), AML (OR = 0.76, 95% CI = 0.59–0.98, p = 0.034), and CML (OR = 0.72 95% CI = 0.51–1.0, ns). We observed the same trend for NHL (OR = 0.47 95% CI = 0.26–0.88 p = 0.017). For ALL, the protective effect of the KIR AA genotype was greater in the presence of KIR ligands C1 (Pc = 0.01) and Bw4 (Pc = 0.001), which are tightly linked in East Asians. By contrast, the C2 ligand strengthened protection from CML (Pc = 0.004). NK cells isolated from KIR AA individuals were significantly more cytotoxic toward leukemic cells than those from other KIR genotypes (p < 0.0001). These data suggest KIR allotypes encoded by East Asian KIR A haplotypes are strongly inhibitory, arming NK cells to respond to leukemogenic cells having altered HLA expression. Thus, the study of populations with distinct KIR and HLA distributions enlightens understanding of immune mechanisms that significantly impact leukemia pathogenesis.
Introduction
Natural killer (NK) cells can detect and eliminate leukemia cells (). To perform this function, NK cells express multiple inhibitory and activating receptors specific for the HLA class I proteins that are expressed on tissue cells (–). Polymorphic inhibitory receptors educate NK cells, enabling them to kill tumors or infected cells having altered or reduced HLA class I expression. In turn, the activating receptors can complement this function by identifying foreign or neo-antigens (–). Study of NK cell receptor polymorphism is thus critical for understanding the anti-tumor effects of NK cells, and to facilitate their subsequent manipulation for use in immunotherapy ().
Killer cell immunoglobulin-like receptors (KIR) are the most polymorphic of the human NK cell receptors, and are characterized by extreme variation in gene number (). The KIR locus comprises up to 13 functional genes, KIR2DL1-5, KIR2DS1-5, KIR3DL1/S1, and KIR3DL2-3, where those designated with an “L” suffix encode inhibitory receptors and those with an “S” suffix encode activating receptors. Previous population, transplantation, and disease association studies have identified two functionally distinct groups of KIR haplotypes that are present in all populations (, ). KIR A haplotypes express four inhibitory receptors specific for polymorphic HLA class I (KIR2DL1, KIR2DL3, KIR3DL1, KIR3DL2), plus zero or one activating receptor (KIR2DS4) also specific for HLA class I. Multiple distinct KIR B haplotypes express fewer inhibitory receptors but more activating receptors than the KIR A haplotypes ().
Although previous examinations of the role of KIR diversity in leukemia control have yielded potentially conflicting results (), a common feature spanning multiple different populations and leukemia subtypes (–) is they are consistent with the KIR A haplotype being more frequent in healthy controls than in the patient groups. Because AA homozygous individuals possess all four inhibitory receptors and few activating receptors specific for HLA class I, their NK cells may therefore be strongly educated to detect any changes in HLA expression that can occur on leukemia cells (). However, both the frequency of the KIR A haplotype, and the abundance of KIR ligands are highly variable across human populations (, ). In this context, examination of specific well-defined human populations is critical for understanding the role of KIR in leukemia control. Here we interrogate such effects in three leukemia subtypes sampled in a large cohort from Southern China, where the KIR A haplotype is very frequent.
Materials and Methods
We studied 836 adult ALL, AML or CML patients, 225 pediatric leukemia patients and 745 healthy adult blood donors. We also studied 46 adult NHL patients, all having diffuse large B-cell lymphoma. The mean age of the adult patients analyzed was 31.4 years, controls 32.7 years, and pediatric cases 11.2 years. The patients were recruited from the hematopoietic stem cell transplantation (HSCT) program in the Shenzhen Blood Center, from August 1999 through June 2015. Leukemia diagnosis followed the French-American-British (FAB) classification, based on morphological observations. CML patients were not included if they had elected to take TKI (Imatinib) therapy instead of transplant. The CML patients included were from a study published in China (). The controls were all unrelated healthy blood donors from Shenzhen Blood Center. All subjects are self-identified as Chinese Southern Han. Written informed consent was obtained from all study participants, and from the parents of all pediatric leukemia patients. The cohorts are described in Supplemental Figure 1.
KIR and HLA Genotyping
The presence or absence of all 13 functional KIR genes and two pseudogenes was determined from genomic DNA using the PCR-SSP based KIR Ready Gene kit (Inno-Train Diagnostik GmbH, Kronberg Im Taunus, Germany). It is difficult to distinguish KIR AB heterozygotes from KIR BB homozygotes, but the KIR AA genotype is readily identifiable (, ). Thus, individuals were designated as homozygous for the KIR A haplotype if they tested positive for KIR2DL1, KIR2DL3, KIR2DL4, KIR2DS4, KIR3DL1, KIR3DL2, KIR3DL3, KIR2DP1, KIR3DP1, and no other KIR. KIR2DL3 and KIR3DL1/S1 alleles were determined by Sanger sequencing from a subset of 306 controls, as described (–). HLA class I genes were sequenced as previously described (), and assigned genotypes at two field (distinct polypeptide sequence) resolution using the ASSIGN4.7 software (Conexio Genomics, Applecross, Australia).
Statistical Analysis
For comparison of KIR AA frequencies, data from other populations were obtained from Allelefrequencies.net (). The database was accessed in August 2018 and all populations having KIR AA genotype frequency data were used (N = 154). The populations were divided into major geographic origins as indicated (Amerindian, East Asian, European, Oceanian, sub-Saharan African) (). The frequencies of KIR AA genotypes were compared across these major population groups using ANOVA, and, when comparing two populations, using unpaired t-tests for comparison of means. The tests were performed using GraphPad software (https://www.graphpad.com/). Odds ratio and confidence intervals were calculated using MedCalc (https://www.medcalc.org/calc/). Where appropriate, the analyses were corrected for multiple comparisons using the Bonferroni method.
NK Cell Cytotoxicity Assays
A total of 15 healthy Chinese Southern Han subjects (six of them with KIR AA genotype) were recruited from Shenzhen Blood Center. PBMC were separated by density gradient from 10 ml of peripheral blood and cultured for 13–14 days using the HANK cell in vitro preparation kit (Hank Bioengineering, Shenzhen, China). This process uses irradiated feeder cells expressing IL15 to promote NK cell expansion as described () for use in immunotherapy (, ), and we use it here to produce NK cells for in vitro assay. The expansion process increases the number and activity of NK cells but does not affect KIR expression (, , ). After this process, the mean purity of CD56+CD3− cells was 86 ± 4.5%. Antibodies used were: UCHT1 mouse anti-human CD3-FITC and B159 mouse anti-human CD56-APC (BD Biosciences, CA, USA). Target cells (1 × 106) were incubated for 10 min at 37°C in 1 mL PBS containing 2.5 μmol/L CFSE and then washed in PBS. Target cells were K562, which is a leukemia derived cell line that lacks HLA class I expression, and seven primary leukemic blasts from ALL and AML patients. NK cells and target cells were combined at an effector-target-ratio of 20:1, and incubated in RPMI-1640 complete medium for 4 h at 37°C with 5% CO2. Dead cells were stained using 7-AAD (BD Biosciences, CA, USA) for 15 min as described () and measured using a BD ACCURI C6 flow cytometer (BD Biosciences, CA, USA). Spontaneous lysis of target cells was measured by including a control without NK cells, and this value was subtracted from the total percentage of specific target cell lysis. All experiments were performed in duplicate.
Detection of CD107a Expression
A total of 18 healthy Chinese Southern Han individuals (11 of them with KIR AA genotype) were recruited from Shenzhen Blood Center. NK cells were isolated as above and incubated with K562 cells at a ratio of 5:1 for 5 h at 37°C. Optimal expression of CD107a on the cell surface of degranulating NK cells occurs between 4~6 h post stimulation (). After the first 1 h, 4 μL monensin (BD Biosciences, CA, USA) was added to each 6 mL culture. The cells were then incubated with 2.5 uL human Fc block (BD Biosciences, CA, USA) at 25°C for 15 min, then stained with CD56, CD3, and CD107a mAbs. Antibodies used were: UCHT1 mouse anti-human CD3-FITC, H4A3 mouse anti-human CD107a-PECy5, and B159 mouse anti-human CD56-APC (all BD Biosciences, CA, USA). Flow cytometric analysis was performed using an ACCURI C6 machine (BD Biosciences, CA, USA) Surface expression of CD107a was assessed in CD56+CD3− cells (). To detect spontaneous degranulation, a control sample without target cells was included in every experiment and this information was used for background subtraction of the results.
Results
KIR A Homozygosity Protects Chinese Han From Leukemia and Non-Hodgkin's Lymphoma
The KIR A haplotype is found in all human populations to varying frequency and is associated with NK cell function in protection from infectious diseases and in susceptibility to preeclampsia (, ). Because NK cells also function in tumor control, we tested whether KIR A affects the probability of developing leukemia. We first demonstrated that the frequency of individuals who are homozygous for the KIR A haplotype is significantly greater in East Asians than in other populations (p < 0.0001: Figure 1). We therefore examined a Chinese Southern Han cohort. In this cohort, the frequency of the KIR AA genotype is significantly higher among healthy individuals (55.3%) than adult leukemia patients (47.1%: OR = 0.72, 95% CI = 0.59–0.88, p = 0.0012) or pediatric leukemia patients (44.9%: OR = 0.6, 95% CI = 0.48–0.89, p = 0.006) (Figure 2A). When analyzed by disease subgroup, the KIR AA genotype was more frequent in healthy controls than in ALL, AML or CML patients (Figure 2B). For adult ALL (OR = 0.68, 95% CI = 0.52–0.89, p = 0.004) and AML (OR = 0.76, 95% CI = 0.59–0.98, p = 0.034), and pediatric ALL (OR = 0.64, 95% CI = 0.45–0.92, p = 0.017) the differences are statistically significant (Figure 2B). Similar trends were observed for CML and pediatric AML but did not reach statistical significance (p = 0.052 and p = 0.1, respectively), likely due to smaller sample sizes. We also analyzed non-Hodgkin's lymphoma (NHL) and again observed a protective effect of the KIR AA genotype (OR = 0.47, 95% CI = 0.26–0.88, p = 0.017), which was statistically significant in spite of the small sample size (N = 46: Figure 2A). The consistent reduction of frequency in all the disease groups analyzed strongly suggests Chinese Southern Han individuals having the KIR AA genotype are protected from developing leukemia or lymphoma.
Figure 1
Figure 2
In KIR A Homozygotes, KIR Ligands Differentially Protect From Leukemia Subtypes
To further explore our findings, we analyzed the HLA class I ligands for KIR. HLA class I allele frequencies can vary significantly between the ethnic groups and geographic regions of China, for example between Northern and Southern Han (, ). Validating the equivalent compositions of our disease and control subjects, there were no significant differences in HLA-A, -B or -C allele or haplotype frequency spectra between controls and the full patient group (Supplemental Figures 2A,B) or from previously published frequencies of HLA class I alleles of Chinese Southern Han (). Neither was there a difference between groups in the frequencies of the −21M variant of HLA-B (Supplemental Figure 2C), which can influence NK cell education ().
KIR bind differentially to specific amino acid motifs present on non-overlapping subsets of HLA class I variants (allotypes). The HLA allotypes are HLA-A*03 and -A*11 (the A3/11 motif: KIR3DL2 ligands), HLA-A and -B allotypes that contain the Bw4 motif (KIR3DL1 ligands), HLA-B and -C allotypes containing the C1 motif (primarily KIR2DL2/3 ligands), and HLA-C allotypes containing the C2 motif (primarily KIR2DL1 ligands). We observed trends for reduction of A3/11 frequency in the adult ALL patients (60.4 vs. 53%, p = 0.02: Figure 2C) and reduction of C2 frequency in the adult AML (29.2 vs. 23.5%, p = 0.047) and CML (29.2 vs. 19.4%, p = 0.01) patients relative to controls. However, these findings did not remain significant after accounting for multiple testing by correction for the four ligands analyzed (Figure 2C). By contrast, combinatorial analysis of KIR and their HLA class I ligands revealed clear distinctions between the three leukemia subsets (Figures 2C,D). The KIR AA plus Bw4 compound genotype was more frequent in healthy individuals (42.4%) than in subjects with ALL (30.2% adult, p = 0.0002, and 31.0% pediatric p = 0.011) or AML (34.8% adult, p = 0.017 and 32.5% pediatric, ns). Similarly, the KIR AA plus C1 compound genotype was significantly more frequent amongst the healthy controls (54.5%) than either the ALL (44.5% adult, p = 0.003 and 43.7% pediatric, p = 0.018) or AML (47.5% adult, p = 0.033 and 44.6% pediatric, ns) patients. For adult ALL, the strength of the effect was consistently greater when considering KIR in combination with HLA class I than when considering KIR genotype alone. These associations remained significant after correction for multiple testing (Pc = 0.001 Bw4+HLA, and Pc = 0.012 C1+HLA: Figure 2C). A reduced frequency of each of the three KIR/ligand combinations was observed in AML and NHL patients, but neither association remained significant after correction, suggesting there is likely limited improvement in protection when considering ligand over that offered by KIR alone.
Analysis of representative world populations revealed distinctively strong linkage disequilibrium between Bw4+HLA-B and C1+HLA-C in East Asians (Figure 3A). HLA-B*46, which is the only frequent HLA-B allele that expresses a C1 motif, is characteristic to East Asia, common in the Chinese Southern Han (~15%), and is linked to C1+HLA-C (, ). Also in the Chinese Southern Han, all KIR A haplotypes express KIR2DL3*001, which is an inhibitory receptor that interacts with C1+HLA-C (). KIR2DL3*001 also predominates in Japanese, another East Asian population (). Similarly, the majority of East Asian KIR A haplotypes (>70%) express either KIR3DL1*001 or KIR3DL1*015, which are strongly-inhibitory KIR3DL1 allotypes specific for Bw4+HLA (). Of note, the non-expressed KIR3DL1*004 allele that is common in Europeans () is absent from East Asians (–). Thus, in addition to the high frequency of KIR A haplotypes, a given Chinese Southern Han KIR A haplotype is more than twice as likely to express KIR3DL1*001 or *015 than a European KIR A haplotype (Figure 3B). Taken together, these observations indicate that Chinese Southern Han having the KIR AA genotype are more likely to have multiple interactions between strongly inhibitory KIR and their corresponding HLA ligands than are European individuals, who are more commonly studied. Thus, these strongly inhibitory KIR/HLA ligand interactions may be important factors in the NK cell control of ALL and to a lesser degree AML, in Chinese Southern Han.
Figure 3
In contrast to ALL and AML, the strength of protection from developing CML offered by the KIR AA genotype was enhanced by the C2 ligand (OR = 0.34, 95% CI = 0.18–0.65, Pc = 0.004) to a much greater degree than by the C1/Bw4 ligands (Figure 2C). Although the sample size of the CML group is relatively small, the strength of the effect provides compelling evidence that Chinese Southern Han individuals having the KIR AA plus C2+HLA-C genotype are less likely to develop CML than those who do not have this compound genotype. KIR2DL1, which is encoded by the KIR A haplotype, has high specificity for C2+HLA-C, but the strength of inhibition resulting from the interaction varies across KIR2DL1 subtypes (, ). In East Asians, only two KIR2DL1 subtypes predominate: strongly inhibitory KIR2DL1*003 encoded by KIR A haplotypes, and weakly inhibitory KIR2DL1*004 encoded by KIR B haplotypes () (and ZD, PJN unpublished). NK cells expressing KIR2DL1*003 are highly effective at killing cells that have reduced C2+HLA-C expression, whereas those expressing KIR2DL1*004 are significantly weaker in this regard (, ). Thus, individuals who have the KIR AA genotype in this population likely have two genomic copies of the strongly inhibiting KIR2DL1 allotype, KIR2DL1*003, and should be best equipped at recognizing and eliminating tumor cells having aberrant C2+HLA-C expression. This could take the form of reduction or loss of HLA expression (, ) or presentation of tumor neoantigen peptides by the C2+HLA-C molecules, abrogating interaction with KIR2DL1 (, ). This result thus identifies a role for strongly functional KIR2DL1 allotypes in protection from CML in Chinese Han. This finding strikes strong parallel with several reproductive diseases including preeclampsia, where the KIR AA genotype in combination with C2 ligand impacts the risk of developing the disease (), and highlights the need for study of populations in addition to Europeans for understanding mechanisms of protection from disease ().
In summary, the KIR AA genotype confers differential protection against ALL, AML, CML and NHL; whereas the C2 ligand enhances KIR AA protection against CML, the C1 and Bw4 ligands enhance KIR AA protection against ALL, and the KIR AA genotype alone protects from AML and NHL.
NK Cells From KIR A Homozygous Individuals Are Strongly Cytotoxic to Leukemia Cells
To investigate whether individuals who are KIR AA+ have stronger responses to leukemia than other individuals do, we measured the cytotoxicity against leukemia cells, of NK cells from healthy individuals with and without KIR AA. We used standard assays for NK cell activity (, , , ) to study six healthy Chinese Southern Han subjects having the KIR AA genotype, and nine without it. NK cells from individuals with the KIR AA genotype demonstrated significantly greater cytotoxicity against the HLA class I negative K562 leukemia cell line, a standard target used to measure NK cell cytotoxicity, than NK cells isolated from individuals without the AA genotype (mean killing 65 vs. 49%, p = 0.002; Figure 4A). These data suggest that NK cells from KIR AA individuals are better educated to respond to loss of HLA class I expression than are those from individuals carrying other KIR genotypes.
Figure 4
Next, the same NK cells were tested for their cytotoxicity against primary leukemic cells, which were isolated from five patients with AML and two with ALL (Supplemental Figure 3). Again, NK cells isolated from individuals who have the KIR AA genotype had a stronger cytotoxic response than those isolated from non-AA individuals (29.9 vs. 20.5%, p = 0.001; Figure 4B). In line with previous observations on European subjects (), we observed variability in the baseline killing of leukemia cells across NK cell donors. This may be due to differential education by donor HLA class I allotypes. Although all the donors and patients studied possess the C1 ligand (Supplemental Figure 3), the very high number of subjects required precludes from completely controlling for HLA variation. Nevertheless, a higher degree of cytotoxicity was consistently observed for NK cells from KIR AA+ donors as compared to KIR AA− donors regardless of the source of leukemia cell targets (p < 0.0001 by paired t-test; Figure 4C). This result was further validated in an independent set of 18 healthy individuals and a complementary method, which measured the cell surface expression of CD107a that accompanies degranulation of NK cells (). Here, a mean of 25.4 % NK cells from KIRAA+ subjects expressed CD107a following contact with K562 cells, vs. 17.8% from non-AA subjects (p = 0.03: Figure 4D). In summary, these findings show that in Chinese Southern Han, NK cells from individuals with the KIR AA genotype exert greater cytotoxicity against leukemic cells than those from individuals having any other KIR genotype.
Discussion
KIR are highly polymorphic and the distribution of KIR haplotypes differs markedly across human populations (). Because multiple previous studies were consistent with, but had not identified, a role for the KIR AA genotype in protection from leukemia, we studied a cohort of Chinese Southern Han from East Asia, where the KIR A haplotype is at highest frequency. We showed that possession of the KIR AA genotype protects individuals from developing either of three types of leukemia and one lymphoma analyzed. These results were further refined by considering the HLA ligands for KIR in the analysis. Importantly, in vitro functional assays demonstrated that NK cells from healthy individuals with KIR AA genotype had strong reactivity against leukemic cells.
Because there have been conflicting reports on the role of KIR in protection from leukemia, we sought to study a large number of patients and controls that were matched by clearly defined age and ethnicity. We studied over 1,000 patients and 745 controls and observed no significant difference in the overall frequency distribution of HLA class I alleles between the two groups or compared with previous studies of Chinese Southern Han (). Because allele frequency spectra of these highly polymorphic genes characterize and distinguish human populations (), this finding validates the self-defined ancestry of the cohorts as well as their suitability for use in these comparisons. We studied subjects from China, where the mean age for leukemia diagnosis is younger than that seen for Europeans (–), and compared them with healthy blood donors having a similar age distribution. The high frequency of KIR AA homozygous individuals observed in our control group is consistent with multiple previous studies of healthy Chinese Southern Han (, –). We also observed a difference in KIR AA frequency between AML patients and controls of similar magnitude as that seen in a previous study of Chinese subjects from Singapore (). One caveat is that we studied patients selected as potential transplantation recipients, but we consider it unlikely that KIR/HLA genotype affects this initial decision. Moreover, our finding that possession of the KIR AA genotype consistently protects from at least three types of leukemia and one lymphoma was further validated by in vitro functional assays of NK cell cytotoxicity and lytic granule release.
Our results suggest that any contradictory findings in the association of KIR/HLA and leukemia that have been reported are likely to be driven by the relative population frequencies of the KIR A haplotypes and the KIR allotypes they encode. Many previous studies are concordant with the results presented here (, , –, ), while others that differ could be explained by cryptic differences in ancestries of patients and controls, resulting in differences in KIR A haplotype frequencies, and the alleles present on those haplotypes (, , –). We hypothesize that the strongly inhibitory KIR allotypes that characterize East Asian KIR haplotypes efficiently educate NK cells to detect changes in HLA expression that can occur in leukemic cells (, ). These changes may include loss or reduction of HLA expression, or alterations to the peptide repertoire that may be detected by NK cells through KIR binding (, ). That the same strongly-inhibiting KIR allotypes are present in other populations, albeit to lower frequencies (), suggests their effect on leukemia control will be similar, but much larger sample sizes will be required to detect them. We studied a population with relatively low genetic diversity, such that possession of the KIR A haplotype is accompanied by a uniform set of alleles of the component genes. High resolution analyses will likely be required to account for the substantial allelic diversity of any further population groups analyzed ().
Although the KIR AA genotype appears to protect against leukemia and NHL in Chinese Southern Hans, there were differences in the HLA ligand required to strengthen this effect. These differences may be due to the distinct maturation stages and the non-overlapping subsets of lineage markers and fusion proteins expressed by the leukemia subtypes (, ). The clearest ligand effect was observed for CML, which is likely due to the relative uniformity of this leukemia class (). The least defined ligand effect was observed for AML, likely due to the highly heterogeneous nature of this disease (). Another ligand for inhibitory NK cell receptors is HLA-E, which can present peptide fragments derived from the expression of HLA-A, -B or -C. CD94/NKG2A educates and modulates NK cells through binding HLA-E, complementing the role of KIR and extending the capacity of NK cells to monitor HLA class I expression (, ). Although CD94/NKG2A is conserved, diversity across individuals occurs because only some HLA-B alleles encode the −21M variant peptide that binds to HLA-E. Accordingly, a recent study showed −21M HLA-B associated with increased survival for AML patients receiving immune therapy (). However, as we observed no significant differences in HLA allele frequencies (Supplemental Figure 2), there were no significant differences in frequency of the HLA-B −21M variant between the patients and controls of our study. This indicates that HLA-B −21M polymorphism is unlikely to affect leukemia susceptibility in Chinese Southern Han.
Although the endpoint of our analysis was disease diagnosis, a limitation is that we were not able to perform functional tests of NK cells extracted from leukemia patients, or phenotypic analyses of HLA expression by leukemic cells. The exception is the K562 target, which is a leukemia cell known to lack HLA class I expression. Another limitation is that the relatively low sample numbers we obtained for the functional tests did not allow for thorough analysis of KIR ligand mismatches between NK cell donor and leukemia target cells. It will be important in future work to examine how changes in HLA expression or neo-antigen peptide presentation may affect the NK cell response to cells that could be leukemogenic, and ultimately how they may influence disease course. These kinds of experiments require pre-genotyping of large numbers of individuals, and will be greatly facilitated by involving populations where the respective genotypes are at high frequency, including the Chinese Southern Han.
In conclusion, this study enhances our understanding of the immune response to leukemia in the world's largest population group, shows the benefit of including specific well-defined populations in study of diseases involving highly polymorphic genes, and identifies a role for NK cell education by inhibitory KIR in control of leukemia. In light of recent studies that promote expansion beyond the usual European populations analyzed (), the benefits of studying well-defined population groups is becoming increasingly apparent.
Statements
Data availability statement
All datasets generated for this study are included in the manuscript and/or the Supplementary Files.
Ethics statement
The study was approved by the ethics review board of Shenzhen Blood Center, Shenzhen, Guangdong, China. Written informed consent was obtained from all study participants.
Author contributions
ZD, SC, JZha, WH, MC, and PN: experiment and analysis design. ZD, YQ, RC, JZhu, JZhe, MZ, GZ, LH, HZ, LL, and WZ: data generation and interpretation. ZD, YQ, QY, MM, XG, and PN: Data analysis. ZD, JZha, XG, MM, MC, and PN: Manuscript preparation. All authors: manuscript review and approval.
Acknowledgments
This study was supported by the National Natural Science Foundation of China (Number: 81373158). PN was supported by start-up funds from the University of Colorado and by NIH U01 AI090905 and U19 NS095774. YQ, MM, XG, and MC were supported by federal funds from the National Cancer Institute, National Institutes of Health (NIH) under contract HHSN261200800001E. The content of this publication does not necessarily reflect the views or policies of the Department of Health and Human Services, nor does mention of trade names, commercial products, or organizations imply endorsement by the US Government. YQ, MM, XG, and MC were also supported by the Intramural Research Program of the NIH, Frederick National Lab, Center for Cancer Research. We thank the patients and marrow registry donors for generously providing blood samples for this study. Thanks to Lisbeth Guethlein and Jill Hollenbach for critical comments.
Conflict of interest
MZ is an employee of Shenzhen Hank Bioengineering Institute. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2019.01646/full#supplementary-material
References
1.
LocatelliFPendeDFalcoMDella ChiesaMMorettaAMorettaL. NK cells mediate a crucial graft-versus-leukemia effect in haploidentical-HSCT to cure high-risk acute leukemia. Trends Immunol. (2018) 39:577–90. 10.1016/j.it.2018.04.009
2.
VivierERauletDHMorettaACaligiuriMAZitvogelLLanierLLet al. Innate or adaptive immunity? The example of natural killer cells. Science. (2011) 331:44–9. 10.1126/science.1198687
3.
LongEOKimHSLiuDPetersonMERajagopalanS. Controlling natural killer cell responses: integration of signals for activation and inhibition. Annu Rev Immunol. (2013) 31:227–58. 10.1146/annurev-immunol-020711-075005
4.
SaundersPMVivianJPO'connorGMSullivanLCPymmPRossjohnJet al. A bird's eye view of NK cell receptor interactions with their MHC class I ligands. Immunol Rev. (2015) 267:148–66. 10.1111/imr.12319
5.
DasJKhakooSI. NK cells: tuned by peptide?Immunol Rev. (2015) 267:214–27. 10.1111/imr.12315
6.
FreudAGMundy-BosseBLYuJCaligiuriMA. The broad spectrum of human natural killer cell diversity. Immunity. (2017) 47:820–33. 10.1016/j.immuni.2017.10.008
7.
HolzemerAGarcia-BeltranWFAltfeldM. Natural killer cell interactions with classical and non-classical human leukocyte antigen class I in HIV-1 infection. Front Immunol. (2017) 8:1496. 10.3389/fimmu.2017.01496
8.
DaherMRezvaniK. Next generation natural killer cells for cancer immunotherapy: the promise of genetic engineering. Curr Opin Immunol. (2018) 51:146–53. 10.1016/j.coi.2018.03.013
9.
UhrbergMValianteNMShumBPShillingHGLienert-WeidenbachKCorlissBet al. Human diversity in killer cell inhibitory receptor genes. Immunity. (1997) 7:753–63. 10.1016/S1074-7613(00)80394-5
10.
BashirovaAAMartinMPMcvicarDWCarringtonM. The killer immunoglobulin-like receptor gene cluster: tuning the genome for defense. Annu Rev Genomics Hum Genet. (2006) 7:277–300. 10.1146/annurev.genom.7.080505.115726
11.
BeziatVHiltonHGNormanPJTraherneJA. Deciphering the killer-cell immunoglobulin-like receptor system at super-resolution for natural killer and T-cell biology. Immunology. (2017) 150:248–64. 10.1111/imm.12684
12.
HsuKCChidaSGeraghtyDEDupontB. The killer cell immunoglobulin-like receptor (KIR) genomic region: gene-order, haplotypes and allelic polymorphism. Immunol Rev. (2002) 190:40–52. 10.1034/j.1600-065X.2002.19004.x
13.
AugustoDG. The impact of KIR polymorphism on the risk of developing cancer: not as strong as imagined?Front Genet. (2016) 7:121. 10.3389/fgene.2016.00121
14.
VerheydenSBernierMDemanetC. Identification of natural killer cell receptor phenotypes associated with leukemia. Leukemia. (2004) 18:2002–7. 10.1038/sj.leu.2403525
15.
VerheydenSDemanetC. Susceptibility to myeloid and lymphoid leukemia is mediated by distinct inhibitory KIR-HLA ligand interactions. Leukemia. (2006) 20:1437–8. 10.1038/sj.leu.2404279
16.
ZhangYWangBYeSLiuSLiuMShenCet al. Killer cell immunoglobulin-like receptor gene polymorphisms in patients with leukemia: possible association with susceptibility to the disease. Leuk Res. (2010) 34:55–8. 10.1016/j.leukres.2009.04.022
17.
VejbaesyaSSae-TamPKhuhapinantASrinakD. Killer cell immunoglobulin-like receptors in Thai patients with leukemia and diffuse large B-cell lymphoma. Hum Immunol. (2014) 75:673–6. 10.1016/j.humimm.2014.04.004
18.
MisraMKPrakashSMoulikNRKumarAAgrawalS. Genetic associations of killer immunoglobulin like receptors and class I human leukocyte antigens on childhood acute lymphoblastic leukemia among north Indians. Hum Immunol. (2016) 77:41–6. 10.1016/j.humimm.2015.10.009
19.
ShenMLinnYCRenEC. KIR-HLA profiling shows presence of higher frequencies of strong inhibitory KIR-ligands among prognostically poor risk AML patients. Immunogenetics. (2016) 68:133–44. 10.1007/s00251-015-0888-4
20.
ReusingSBManserAREnczmannJMulderAClaasFHCarringtonMet al. Selective downregulation of HLA-C and HLA-E in childhood acute lymphoblastic leukaemia. Br J Haematol. (2016) 174:477–80. 10.1111/bjh.13777
21.
Gonzalez-GalarzaFFTakeshitaLYSantosEJKempsonFMaiaMHDa SilvaALet al. Allele frequency net 2015 update: new features for HLA epitopes, KIR and disease and HLA adverse drug reaction associations. Nucleic Acids Res. (2015) 43:D784–8. 10.1093/nar/gku1166
22.
GuethleinLANormanPJHiltonHHParhamP. Co-evolution of MHC class I and variable NK cell receptors in placental mammals. Immunol Rev. (2015) 267:259–82. 10.1111/imr.12326
23.
DengZZhenJWangDHeLZouH. [Association of genetic polymorphisms of KIR-HLA system with chronic myeloid leukemia among ethnic Hans from southern China]. Zhonghua Yi Xue Yi Chuan Xue Za Zhi. (2017) 34:53–7. 10.3760/cma.j.issn.1003-9406.2017.01.012
24.
PyoCWGuethleinLAVuQWangRAbi-RachedLNormanPJet al. Different patterns of evolution in the centromeric and telomeric regions of group A and B haplotypes of the human killer cell Ig-like receptor locus. PLoS ONE. (2010) 5:e15115. 10.1371/journal.pone.0015115
25.
NormanPJAbi-RachedLGendzekhadzeKKorbelDGleimerMRowleyDet al. Unusual selection on the KIR3DL1/S1 natural killer cell receptor in Africans. Nat Genet. (2007) 39:1092–9. 10.1038/ng2111
26.
ZhuFMJiangKLvQFZhangWHeJXuXGet al. [Polymorphism of killer cell immunoglobulin-like receptor 2DL3 gene in Zhejiang Han population]. Zhonghua Yi Xue Yi Chuan Xue Za Zhi. (2007) 24:88–90. 10.3760/j.issn:1003-9406.2007.01.021
27.
ZhenJHeLXuYZhaoJYuQZouHet al. Allelic polymorphism of KIR2DL2/2DL3 in a southern Chinese population. Tissue Antigens. (2015) 86:362–7. 10.1111/tan.12681
28.
ZhenJWangDHeLZouHXuYGaoSet al. Genetic profile of KIR and HLA in southern Chinese Han population. Hum Immunol. (2014) 75:59–64. 10.1016/j.humimm.2013.09.006
29.
ImaiCIwamotoSCampanaD. Genetic modification of primary natural killer cells overcomes inhibitory signals and induces specific killing of leukemic cells. Blood. (2005) 106:376–83. 10.1182/blood-2004-12-4797
30.
ZhangMDanielSHuangYChanceyCHuangQLeiYFet al. Anti-West Nile virus activity of in vitro expanded human primary natural killer cells. BMC Immunol. (2010) 11:3. 10.1186/1471-2172-11-3
31.
QinZChenJZengJNiuLXieSWangXet al. Effect of NK cell immunotherapy on immune function in patients with hepatic carcinoma: a preliminary clinical study. Cancer Biol Ther. (2017) 18:323–30. 10.1080/15384047.2017.1310346
32.
FujisakiHKakudaHShimasakiNImaiCMaJLockeyTet al. Expansion of highly cytotoxic human natural killer cells for cancer cell therapy. Cancer Res. (2009) 69:4010–7. 10.1158/0008-5472.CAN-08-3712
33.
RomeeRRosarioMBerrien-ElliottMMWagnerJAJewellBASchappeTet al. Cytokine-induced memory-like natural killer cells exhibit enhanced responses against myeloid leukemia. Sci Transl Med. (2016) 8:357ra123. 10.1126/scitranslmed.aaf2341
34.
SchmidIKrallWJUittenbogaartCHBraunJGiorgiJV. Dead cell discrimination with 7-amino-actinomycin D in combination with dual color immunofluorescence in single laser flow cytometry. Cytometry. (1992) 13:204–8. 10.1002/cyto.990130216
35.
AlterGMalenfantJMAltfeldM. CD107a as a functional marker for the identification of natural killer cell activity. J Immunol Methods. (2004) 294:15–22. 10.1016/j.jim.2004.08.008
36.
PenackOGentiliniCFischerLAsemissenAMScheibenbogenCThielEet al. CD56dimCD16neg cells are responsible for natural cytotoxicity against tumor targets. Leukemia. (2005) 19:835–40. 10.1038/sj.leu.2403704
37.
MoffettAColucciF. Co-evolution of NK receptors and HLA ligands in humans is driven by reproduction. Immunol Rev. (2015) 267:283–97. 10.1111/imr.12323
38.
XuSYinXLiSJinWLouHYangLet al. Genomic dissection of population substructure of Han Chinese and its implication in association studies. Am J Hum Genet. (2009) 85:762–74. 10.1016/j.ajhg.2009.10.015
39.
ChenCHYangJHChiangCWKHsiungCNWuPEChangLCet al. Population structure of Han Chinese in the modern Taiwanese population based on 10,000 participants in the Taiwan Biobank project. Hum Mol Genet. (2016) 25:5321–31. 10.1093/hmg/ddw346
40.
TrachtenbergEVinsonMHayesEHsuYMHoutchensKErlichHet al. HLA class I (A, B, C) and class II (DRB1, DQA1, DQB1, DPB1) alleles and haplotypes in the Han from southern China. Tissue Antigens. (2007) 70:455–63. 10.1111/j.1399-0039.2007.00932.x
41.
HorowitzADjaoudZNemat-GorganiNBlokhuisJHiltonHGBeziatVet al. Class I HLA haplotypes form two schools that educate NK cells in different ways. Sci Immunol. (2016) 1:eaag1672. 10.1126/sciimmunol.aag1672
42.
SolbergODMackSJLancasterAKSingleRMTsaiYSanchez-MazasAet al. Balancing selection and heterogeneity across the classical human leukocyte antigen loci: a meta-analytic review of 497 population studies. Hum Immunol. (2008) 69:443–64. 10.1016/j.humimm.2008.05.001
43.
ExcoffierLLischerHE. Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Mol Ecol Resour. (2010) 10:564–567. 10.1111/j.1755-0998.2010.02847.x
44.
YawataMYawataNDraghiMLittleAMPartheniouFParhamP. Roles for HLA and KIR polymorphisms in natural killer cell repertoire selection and modulation of effector function. J Exp Med. (2006) 203:633–45. 10.1084/jem.20051884
45.
TaoSDHeYMYingYLHeJZhuFMLvHJ. KIR3DL1 genetic diversity and phenotypic variation in the Chinese Han population. Genes Immun. (2014) 15:8–15. 10.1038/gene.2013.55
46.
DengZZhenJZhuBZhangGYuQWangDet al. Allelic diversity of KIR3DL1/3DS1 in a southern Chinese population. Hum Immunol. (2015) 76:663–6. 10.1016/j.humimm.2015.09.017
47.
BariRBellTLeungWHVongQPChanWKDas GuptaNet al. Significant functional heterogeneity among KIR2DL1 alleles and a pivotal role of arginine 245. Blood. (2009) 114:5182–90. 10.1182/blood-2009-07-231977
48.
HiltonHGGuethleinLAGoyosANemat-GorganiNBushnellDANormanPJet al. Polymorphic HLA-C receptors balance the functional characteristics of KIR haplotypes. J Immunol. (2015) 195:3160–70. 10.4049/jimmunol.1501358
49.
YuJHellerGChewningJKimSYokoyamaWMHsuKC. Hierarchy of the human natural killer cell response is determined by class and quantity of inhibitory receptors for self-HLA-B and HLA-C ligands. J Immunol. (2007) 179:5977–89. 10.4049/jimmunol.179.9.5977
50.
JonckerNTFernandezNCTreinerEVivierERauletDH. NK cell responsiveness is tuned commensurate with the number of inhibitory receptors for self-MHC class I: the rheostat model. J Immunol. (2009) 182:4572–80. 10.4049/jimmunol.0803900
51.
SimMJMalakerSAKhanAStowellJMShabanowitzJPetersonMEet al. Canonical and cross-reactive binding of NK cell inhibitory receptors to HLA-C allotypes is dictated by peptides bound to HLA-C. Front Immunol. (2017) 8:193. 10.3389/fimmu.2017.00193
52.
NakimuliAChazaraOHibySEFarrellLTukwasibweSJayaramanJet al. A KIR B centromeric region present in Africans but not Europeans protects pregnant women from pre-eclampsia. Proc Natl Acad Sci USA. (2015) 112:845–50. 10.1073/pnas.1413453112
53.
FauriatCAnderssonSBjorklundATCarlstenMSchafferMBjorkstromNKet al. Estimation of the size of the alloreactive NK cell repertoire: studies in individuals homozygous for the group A KIR haplotype. J Immunol. (2008) 181:6010–9. 10.4049/jimmunol.181.9.6010
54.
CmaCMA. [Guideline for management of chronic myeloid leukemia (2013)]. Zhonghua Xue Ye Xue Za Zhi. (2013) 34:464–70. 10.3760/cma.j.issn.0253-2727.2013.05.021
55.
HeXFWangQRCenJNQiuHYSunANChenSNet al. [EVI1 expression, clinical and cytogenetical characteristics in 447 patients with acute myeloid leukemia]. Zhonghua Xue Ye Xue Za Zhi. (2016) 37:936–41. 10.3760/cma.j.issn.0253-2727.2016.11.002
56.
YuWJWangYXuLPZhangXHLiuKYHuangXJ. [Influence of donor-recipient sex matching on outcomes of haploidentical hematopoietic stem cell transplantation for acute leukemia]. Zhonghua Xue Ye Xue Za Zhi. (2018) 39:398–403. 10.3760/cma.j.issn.0253-2727.2018.05.010
57.
JiangKZhuFMLvQFYanLX. Distribution of killer cell immunoglobulin-like receptor genes in the Chinese Han population. Tissue Antigens. (2005) 65:556–63. 10.1111/j.1399-0039.2005.00412.x
58.
WuGQZhaoYMLaiXYYangKLZhuFMZhangWet al. Distribution of killer-cell immunoglobulin-like receptor genes in Eastern mainland Chinese Han and Taiwanese Han populations. Tissue Antigens. (2009) 74:499–507. 10.1111/j.1399-0039.2009.01366.x
59.
ShiLShiLTaoYLinKLiuSYuLet al. Distribution of killer cell immunoglobulin-like receptor genes and combinations with HLA-C ligands in an isolated Han population in southwest China. Tissue Antigens. (2011) 78:60–4. 10.1111/j.1399-0039.2011.01683.x
60.
PamukGETozkirHUyanikMSGurkanHDuymazJPamukON. Natural killer cell killer immunoglobulin-like gene receptor polymorphisms in non-Hodgkin lymphoma: possible association with clinical course. Leuk Lymphoma. (2015) 56:2902–7. 10.3109/10428194.2015.1014361
61.
MiddletonDDilerASMeenaghASleatorCGourraudPA. Killer immunoglobulin-like receptors (KIR2DL2 and/or KIR2DS2) in presence of their ligand (HLA-C1 group) protect against chronic myeloid leukaemia. Tissue Antigens. (2009) 73:553–60. 10.1111/j.1399-0039.2009.01235.x
62.
AlmalteZSamaraniSIannelloADebbecheODuvalMInfante-RivardCet al. Novel associations between activating killer-cell immunoglobulin-like receptor genes and childhood leukemia. Blood. (2011) 118:1323–8. 10.1182/blood-2010-10-313791
63.
BaborFManserARFischerJCScherenschlichNEnczmannJChazaraOet al. KIR ligand C2 is associated with increased susceptibility to childhood ALL and confers an elevated risk for late relapse. Blood. (2014) 124:2248–51. 10.1182/blood-2014-05-572065
64.
De SmithAJWalshKMLadnerMBZhangSXiaoCCohenFet al. The role of KIR genes and their cognate HLA class I ligands in childhood acute lymphoblastic leukemia. Blood. (2014) 123:2497–503. 10.1182/blood-2013-11-540625
65.
StringarisKMarinDBarrettAJHillsRSobieskiCCaoKet al. KIR gene haplotype: an independent predictor of clinical outcome in MDS patients. Blood. (2016) 128:2819–23. 10.1182/blood-2016-05-713099
66.
NormanPJHollenbachJANemat-GorganiNMarinWMNorbergSJAshouriEet al. Defining KIR and HLA class I genotypes at highest resolution via high-throughput sequencing. Am J Hum Genet. (2016) 99:375–91. 10.1016/j.ajhg.2016.06.023
67.
HaferlachTKohlmannAWieczorekLBassoGKronnieGTBeneMCet al. Clinical utility of microarray-based gene expression profiling in the diagnosis and subclassification of leukemia: report from the international microarray innovations in leukemia study group. J Clin Oncol. (2010) 28:2529–37. 10.1200/JCO.2009.23.4732
68.
ArberDAOraziAHasserjianRThieleJBorowitzMJLe BeauMMet al. The 2016 revision to the World Health Organization classification of myeloid neoplasms and acute leukemia. Blood. (2016) 127:2391–405. 10.1182/blood-2016-03-643544
69.
CheredaBMeloJV. Natural course and biology of CML. Ann Hematol. (2015) 94(Suppl. 2):S107–21. 10.1007/s00277-015-2325-z
70.
LeyTJMillerCDingLRaphaelBJMungallAJRobertsonAet al. Genomic and epigenomic landscapes of adult de novo acute myeloid leukemia. N Engl J Med. (2013) 368:2059–74. 10.1056/NEJMoa1301689
71.
BraudVMAllanDSO'callaghanCASoderstromKD'andreaAOggGSet al. HLA-E binds to natural killer cell receptors CD94/NKG2A, B and C. Nature. (1998) 391:795–9. 10.1038/35869
72.
HallnerABernsonEHusseinBAEwald SanderFBruneMAureliusJet al. The HLA-B−21 dimorphism impacts on NK cell education and clinical outcome of immunotherapy in acute myeloid leukemia. Blood. (2019) 133:1479–88. 10.1182/blood-2018-09-874990
73.
WojcikGLGraffMNishimuraKKTaoRHaesslerJGignouxCRet al. Genetic analyses of diverse populations improves discovery for complex traits. Nature. (2019) 570:514–518. 10.1038/s41586-019-1310-4
Summary
Keywords
leukemia, AML, ALL, CML, NHL, HLA, KIR, Chinese Southern Han
Citation
Deng Z, Zhao J, Cai S, Qi Y, Yu Q, Martin MP, Gao X, Chen R, Zhuo J, Zhen J, Zhang M, Zhang G, He L, Zou H, Lu L, Zhu W, Hong W, Carrington M and Norman PJ (2019) Natural Killer Cells Offer Differential Protection From Leukemia in Chinese Southern Han. Front. Immunol. 10:1646. doi: 10.3389/fimmu.2019.01646
Received
13 April 2019
Accepted
02 July 2019
Published
16 July 2019
Volume
10 - 2019
Edited by
Alexander David Barrow, The University of Melbourne, Australia
Reviewed by
Amir Horowitz, Icahn School of Medicine at Mount Sinai, United States; Nicole Bernard, McGill University, Canada
Updates
Copyright
© 2019 Deng, Zhao, Cai, Qi, Yu, Martin, Gao, Chen, Zhuo, Zhen, Zhang, Zhang, He, Zou, Lu, Zhu, Hong, Carrington and Norman.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Paul J. Norman paul.norman@ucdenver.edu
This article was submitted to NK and Innate Lymphoid Cell Biology, a section of the journal Frontiers in Immunology
†These authors have contributed equally to this work
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.