• KIR2DL5B is associated with poor molecular response and transformation-free survival in CML patients enrolled to the TIDEL-II study.

  • KIR genotyping would select out high risk CML patients at baseline and allow better targeting of novel interventions.

Killer immunoglobulin-like receptors (KIRs) on natural killer (NK) cells have been shown to predict for response in chronic phase-chronic myeloid leukemia (CP-CML) patients treated with tyrosine kinase inhibitors. We performed KIR genotyping in 148 newly diagnosed CP-CML patients treated with a novel sequential imatinib/nilotinib strategy aimed at achievement of optimal molecular responses at defined time points. We found the presence of KIR2DL5B to be associated with inferior transformation-free survival and event-free survival and an independent predictor of inferior major molecular response (BCR-ABL1 ≤0.1%) and molecular response 4.5 (BCR-ABL1 ≤0.0032%). This suggests a critical early role for NK cells in facilitating response to imatinib that cannot be overcome by subsequent intensification of therapy. KIR genotyping may add valuable prognostic information to future baseline predictive scoring systems in CP-CML patients and facilitate optimal frontline treatment selection.

The majority of chronic phase-chronic myeloid leukemia (CP-CML) patients treated with tyrosine kinase inhibitors (TKIs) have excellent transformation-free survival (TFS), with only 2% to 7% progressing to accelerated phase or blast crisis.1  Prognostic biomarkers that reliably identify these high-risk patients remain elusive, although individual immune response profiles may contribute to the differences in treatment outcomes.

Killer immunoglobulin-like receptors (KIRs), expressed by natural killer (NK) cells, are an integral element of the innate immune response.2  The KIR genetic locus has 16 highly polymorphic genes, coding for either activating or inhibitory KIRs. KIRs participate in NK cell-mediated cell killing of virally infected or tumor cells, in which the ligand may be altered or missing. Normal host cells with appropriate ligand expression are protected by interaction through the inhibitory KIRs.3 

KIR alleles are highly individualized, and each of the 16 genes may either be present or absent. These genes are organized into 2 broad haplotypes: A and B.3  KIR genotype variations are associated with cancer treatment outcomes and predisposition to immune disorders.4,5  Importantly, donor B haplotypes are associated with improved survival in the context of unrelated allogeneic hematopoietic cell transplants for acute myeloid leukemia.6 

In this study, we investigated the link between KIR genotypes and treatment outcomes in CP-CML patients treated in the Therapeutic Intensification in De Novo Leukaemia (TIDEL-II) study, using dose adapted frontline imatinib followed by rapid switching to nilotinib for failing to achieve consensus optimal responses.7,8 

In TIDEL-II, all 210 patients started treatment with imatinib 600 mg/day.9  The imatinib dose was escalated to 800 mg/day if serum imatinib trough levels at day 22 were <1000 ng/mL. Patients who failed to achieve any of the predetermined molecular targets (BCR-ABL1 ≤10%, ≤1%, and ≤0.1% at 3, 6, and 12 months respectively, on the international scale) were switched to nilotinib 400 mg twice daily with or without an antecedent trial of imatinib 800 mg/day. KIR genotyping was performed in 148 patients (samples unavailable in 62 patients) with the Genotyping SSP Kit (Invitrogen, Carlsbad, CA). Genotypes were correlated with molecular response and survival outcomes. Survival probabilities were calculated using the Kaplan-Meier method.10  Cumulative incidence probabilities were calculated using Fine and Gray’s method, with study withdrawal for any reason as competing risks.11  The Akaike information criterion was used for optimum model selection in competing risk multivariate analyses. All statistical analyses were done using R.12-14  This research was approved by the institutional Human Research Ethics Committee and conducted in accordance with the Declaration of Helsinki.

Characteristics and outcomes of the 148 patients in the KIR substudy were similar to the TIDEL-II cohort in entirety (supplemental Table 1 available on the Blood Web site). Univariate analyses correlated achievement of major molecular response (MMR; BCR-ABL1 ≤0.1%) with early molecular response (EMR; BCR-ABL1 ≤10% at 3 months), Sokal score, and the absence of KIR2DL5B, KIR2DS3, and KIR2DL2 (Table 1). EMR is a particularly strong predictor of subsequent molecular outcomes in our cohort, as demonstrated in other studies.8  Among different KIR genotypes, patients who had the KIR2DL5B gene (KIR2DL5BPOS) had inferior achievement of MMR (Figure 1A) and molecular response 4.5 (MR4.5; BCR-ABL1 ≤ 0.0032%; Figure 1B). KIR2DL5BPOS was also associated with inferior TFS and event-free survival (EFS) (Figure 1C-D; supplemental Tables 2 and 3) but not overall survival (data not shown).

Table 1

Variables associated with the achievement of MMR: univariate and multivariate analyses

n (%)HR(95% CI)P value
Univariate analysis     
 KIR     
  KIR2DL5B 31 (21%) 0.423 (0.262-0.682) <.001 
  KIR2DS3 44 (30%) 0.547 (0.368-0.811) .003 
  KIR2DL2 83 (56%) 0.607 (0.428-0.859) .005 
  KIR2DS2 82 (55%) 0.71 (0.501-1.01) .055 
  KIR3DS1 41 (28%) 0.717 (0.481-1.07) .1 
  KIR2DS5 51 (34%) 0.758 (0.519-1.11) .15 
  KIR2DS1 56 (38%) 0.833 (0.575-1.21) .33 
  KIR2DL5A 44 (30%) 0.836 (0.562-1.25) .38 
  KIR2DL3 138 (93%) 1.29 (0.544-3.05) .57 
  KIR2DS4DEL* 129 (87%) 1.1 (0.668-1.82) .7 
  KIR2DS4 60 (41%) 1.02 (0.716-1.44) .93 
 Sokal     
  Low vs High  0.648 (0.446-0.942) .023 
  Low vs Intermediate  0.842 (0.591-1.198) .34 
  Intermediate vs High  0.77 (0.512-1.16) .21 
 BCR-ABL1 at 3 months IS     
  1-10% vs <1%  0.459 (0.3019-0.699) .003 
  >10% vs <1%  0.164 (0.0749-0.361) <.001 
  >10% vs 1-10%  0.358 (0.164-0.782) <.001 
Multivariate analysis     
 KIR     
  KIR2DL5B  0.52 (0.284-0.951) .034 
  KIR2DL2  0.62 (0.3804-1.01) .055 
BCR-ABL1 at 3 months, IS     
  1-10% vs <1%  0.164 (0.2395-0.585) <.001 
  >10% vs <1%  0.13 (0.0605-0.279) <.001 
  >10% vs 1-10%  0.347 (0.164-0.735) <.001 
n (%)HR(95% CI)P value
Univariate analysis     
 KIR     
  KIR2DL5B 31 (21%) 0.423 (0.262-0.682) <.001 
  KIR2DS3 44 (30%) 0.547 (0.368-0.811) .003 
  KIR2DL2 83 (56%) 0.607 (0.428-0.859) .005 
  KIR2DS2 82 (55%) 0.71 (0.501-1.01) .055 
  KIR3DS1 41 (28%) 0.717 (0.481-1.07) .1 
  KIR2DS5 51 (34%) 0.758 (0.519-1.11) .15 
  KIR2DS1 56 (38%) 0.833 (0.575-1.21) .33 
  KIR2DL5A 44 (30%) 0.836 (0.562-1.25) .38 
  KIR2DL3 138 (93%) 1.29 (0.544-3.05) .57 
  KIR2DS4DEL* 129 (87%) 1.1 (0.668-1.82) .7 
  KIR2DS4 60 (41%) 1.02 (0.716-1.44) .93 
 Sokal     
  Low vs High  0.648 (0.446-0.942) .023 
  Low vs Intermediate  0.842 (0.591-1.198) .34 
  Intermediate vs High  0.77 (0.512-1.16) .21 
 BCR-ABL1 at 3 months IS     
  1-10% vs <1%  0.459 (0.3019-0.699) .003 
  >10% vs <1%  0.164 (0.0749-0.361) <.001 
  >10% vs 1-10%  0.358 (0.164-0.782) <.001 
Multivariate analysis     
 KIR     
  KIR2DL5B  0.52 (0.284-0.951) .034 
  KIR2DL2  0.62 (0.3804-1.01) .055 
BCR-ABL1 at 3 months, IS     
  1-10% vs <1%  0.164 (0.2395-0.585) <.001 
  >10% vs <1%  0.13 (0.0605-0.279) <.001 
  >10% vs 1-10%  0.347 (0.164-0.735) <.001 

In all proportional hazard analyses involving KIRs, patients negative for the genotype are assigned a hazard ratio (HR) of 1.0. The genes KIR3DL3, KIR3DP1, KIR2DL4, and KIR3DL2 have prevalence of close to 100% in our cohort (similar to their known prevalence in the general population) and were excluded from further analysis. Univariate analyses revealed KIR2DL5B, KIR2DS3, KIR2DL2, and early molecular response (EMR; BCR-ABL1% at 3 months) to be significantly correlated with the achievement of MMR. Models including all combination of these variables were tested. Results from the optimum model selected using the Akaike information criteria is listed here and include only the input variables of KIR2DL2, KIR2DL5B, and EMR. Only the latter two were shown to be independently associated with molecular response. The inferior prognostic risk conferred by having BCR-ABL1 >10% at 3 months is particularly strong, even though 15 of 18 of the EMR failure patients had exposure to the more potent TKI nilotinib (2 transformed to BC and 1 patient withdrew prior to switching). KIR2DL5B and EMR were also associated with the achievement of MR4.5 in a multivariate analysis (KIR2DL5BPOS: HR, 0.50; P = .05; BCR-ABL1 >10% at 3 months: HR, 0.07; P = .005; 1-10% BCR-ABL1 at 3 months: HR, 0.36; P = .0002). Age and female sex did not have a statistically significant correlation with the achievement of MMR; nilotinib exposure was associated with inferior achievement of MMR, as is expected from the treatment schema of the TIDEL-II study (data not listed). (The multivariate model that included all 3 KIR genes with EMR is not listed here, as it does not optimally describe the data, and collinearity interferes with accurate reporting of hazard ratios).

*

KIR2DS4DEL refers to the 2DS4 allelic variant bearing a 22-bp deletion, also referred to as KIR1D.

Figure 1

Correlation between KIR2DL5B and treatment outcomes in patients treated with imatinib and nilotinib in the TIDEL-II study.KIR2DL5BPOS patients have inferior achievement of molecular responses: (A) cumulative incidence of MMR and (B) MR4.5. KIR2DL5BPOS patients also have inferior TFS and EFS compared with KIR2DL5BNEG patients. (C) TFS events include transformation to accelerated phase and blast crisis, as well as death from any cause. (D) EFS events include TFS events and loss of MMR or BCR-ABL1 increasing to a level >1% from a nadir ≤1%, kinase domain mutations, and discontinuation of TIDEL-II treatment (imatinib and/or nilotinib) for any cause. The difference in overall survival as segregated by KIR2DL5B status is not statistically significant (data not shown).

Figure 1

Correlation between KIR2DL5B and treatment outcomes in patients treated with imatinib and nilotinib in the TIDEL-II study.KIR2DL5BPOS patients have inferior achievement of molecular responses: (A) cumulative incidence of MMR and (B) MR4.5. KIR2DL5BPOS patients also have inferior TFS and EFS compared with KIR2DL5BNEG patients. (C) TFS events include transformation to accelerated phase and blast crisis, as well as death from any cause. (D) EFS events include TFS events and loss of MMR or BCR-ABL1 increasing to a level >1% from a nadir ≤1%, kinase domain mutations, and discontinuation of TIDEL-II treatment (imatinib and/or nilotinib) for any cause. The difference in overall survival as segregated by KIR2DL5B status is not statistically significant (data not shown).

Close modal

KIR2DL2POS and KIR2DS3POS were also associated with inferior achievement of MMR and are commonly coinherited with KIR2DL5B as an allele.15  In our cohort, all 31 patients positive for KIR2DL5B also had the KIR2DL2 gene, and 27 of 31 (87%) also had the KIR2DS3 gene (supplemental Figure 1). These frequencies are similar to data from the Allele Frequency Database compared with predominantly Caucasian cohorts.16,17  Due to linkage disequilibrium between KIR genes, we suspected the prognostic power of KIR2DL2 and KIR2DS3 were secondary to their association with KIR2DL5B. Indeed, a multivariate analysis including the 3 KIR genotypes, Sokal score, and EMR as variables demonstrated the KIR2DL5B status and EMR to be the only independent predictors for achievement of MMR (Table 1). Consistent with an independent relationship, there was no correlation between KIR2DL5B genotype and EMR: 5 of 31 KIR2DL5BPOS patients failed to achieve EMR compared with 13 of 116 KIR2DL5BNEG patients (16% vs 11%, P = .54, EMR data missing for 1 patient).

Ours is the first study showing the prognostic significance of KIR genotypes in patients treated (sequentially) with nilotinib. In TIDEL-II, patients at risk of inferior outcomes were rapidly switched from imatinib to nilotinib. KIR2DL5BPOS and KIR2DL5BNEG patients who received nilotinib treatment were similar in proportion (26% vs 21%). However, their outcomes are different: 1 of 8 (13%) and 11 of 25 (56%) subsequently achieved MMR in the KIR2DL5BPOS and KIR2DL5BNEG groups, respectively (supplemental Table 4). Our findings suggest that even with the potent second-generation TKI nilotinib, KIR genotypes, a predetermined genetic host factor, may still be one of the most discriminatory prognostic markers available at baseline. Although the biological mechanism that underpins this observed association remains to be elucidated, KIR2DL5B encodes an inhibitory “orphan” KIR receptor (its ligand is unknown),18  but its absence may increase efficiency of NK-mediated killing of leukemic stem cells. The importance of immune responses in CML disease control with TKI therapy is consistent with effects seen with allogeneic stem cell transplant, donor lymphocyte infusions, and interferon-alfa treatment, in which the therapeutic effect is wholly or partially secondary to T- and NK-cell activity.19,20  Interestingly, increased numbers of oligoclonal T and NK cells occur with dasatinib therapy and are associated with a better prognosis.21  Additionally, certain KIR genotypes have previously been reported to be overrepresented in CML patients.22  However, KIR genotype frequencies observed in our cohort are similar to those observed in other Caucasian populations reported in the Allele Frequency Database.

Several groups have claimed a link between KIR genotypes and treatment response in CP-CML. In a British cohort treated with imatinib 400 mg/day, KIR2DS1NEG was correlated with superior achievement of complete cytogenetic response; KIR2DL5B was correlated with response by univariate analysis only in this study.23  However, in a follow-up report, there was no correlation between KIR genotypes and outcomes in their dasatinib-treated patients.24  In contrast, Kreutzman et al25  showed an association between superior molecular response and KIR2DL5A/BNEG in a dasatinib-treated cohort. Observations of the British imatinib and Scandinavian cohorts, together with the current study, demonstrate the strongest link between KIR2DL5B of the KIR cluster and treatment outcomes.

EMR is also known to be strongly associated with treatment outcomes and may also be used to guide treatment decisions. The independent prognostic significance of EMR and KIR suggests their prognostic information may be additive. However, EMR information is only available 3 months after treatment commencement. In contrast, KIR2DL5B can identify, at baseline, the 20% of patients with a transformation risk of ∼10% over 2 years vs the 80% of patients with a transformation risk of <3%. Incorporation of KIR2DL5B in future baseline prognostic scores, together with other predictive markers, may enable targeted interventions to improve outcomes at the earliest available opportunity, if the presence and magnitude of KIR2DL5B’s prognostic influence can be further validated and confirmed in prospective studies. Additionally, further understanding the role of KIR and cytotoxic cells in eradicating leukemic clones will translate to specific salvage treatment strategies aimed at boosting immune responses in poor-risk patients thus identified, such as the addition of interferon-alfa.

Presented at the American Society of Hematology Meeting, San Francisco, CA, December 6-9, 2014.

The online version of this article contains a data supplement.

The publication costs of this article were defrayed in part by page charge payment. Therefore, and solely to indicate this fact, this article is hereby marked “advertisement” in accordance with 18 USC section 1734.

This paper presents results of a correlative study of TIDEL-II, a clinical trial sponsored by the Australasian Leukaemia and Lymphoma Group. A list of contributing investigators of this clinical study appears in supplemental Materials. The authors acknowledge the ongoing assistance of Stephanie Arbon, Bronwyn Cambareri, Jodi Braley, and all the site coordinators of the TIDEL-II study.

This study was supported financially by Novartis, Australian National Health and Medical Research Council grant APP1059165, and the Royal Adelaide Hospital Contributing Haematologists’ Committee.

Novartis had no role in the analysis or presentation of these data.

Contribution: D.T.Y. contributed to research design, gathered and analyzed the data, and drafted the manuscript; S.B. and D.L.W. contributed to research design and critically reviewed the manuscript; T.P.H. and A.S.Y. contributed to research design, supervised the research, and critically reviewed the manuscript; and C.T. and L.V. performed the research and critically reviewed the manuscript.

Conflict-of-interest disclosure: D.L.W., D.T.Y., and T.P.H. acted as a consultant or advisor to and received research funding and honoraria from Novartis, BMS, and Ariad; T.P.H. acted as a consultant or advisor to and received honoraria from Pfizer; S.B. acted as a consultant or advisor to and received research funding and honoraria from Novartis and BMS and received research funding from Ariad; A.S.Y. received research funding and honoraria from Novartis. D.T.Y. received scholarship funding from the Leukaemia Foundation of Australia and the Royal Adelaide Research Foundation AR Clarkson scholarship. The remaining authors declare no competing financial interests.

Correspondence: Agnes Yong, SA Pathology, PO Box 14, Rundle Mall, Adelaide, SA 5000, Australia; e-mail: agnes.yong@sa.gov.au.

1
Hughes
 
T
White
 
D
Which TKI? An embarrassment of riches for chronic myeloid leukemia patients.
Hematology Am Soc Hematol Educ Program
2013
, vol. 
2013
 (pg. 
168
-
175
)
2
Valiante
 
NM
Uhrberg
 
M
Shilling
 
HG
et al. 
Functionally and structurally distinct NK cell receptor repertoires in the peripheral blood of two human donors.
Immunity
1997
, vol. 
7
 
6
(pg. 
739
-
751
)
3
Parham
 
P
Guethlein
 
LA
Pregnancy immunogenetics: NK cell education in the womb?
J Clin Invest
2010
, vol. 
120
 
11
(pg. 
3801
-
3804
)
4
Purdy
 
AK
Campbell
 
KS
Natural killer cells and cancer: regulation by the killer cell Ig-like receptors (KIR).
Cancer Biol Ther
2009
, vol. 
8
 
23
(pg. 
2211
-
2220
)
5
Zitvogel
 
L
Apetoh
 
L
Ghiringhelli
 
F
Kroemer
 
G
Immunological aspects of cancer chemotherapy.
Nat Rev Immunol
2008
, vol. 
8
 
1
(pg. 
59
-
73
)
6
Cooley
 
S
Trachtenberg
 
E
Bergemann
 
TL
et al. 
Donors with group B KIR haplotypes improve relapse-free survival after unrelated hematopoietic cell transplantation for acute myelogenous leukemia.
Blood
2009
, vol. 
113
 
3
(pg. 
726
-
732
)
7
NCCN. Clinical Practice Guidelines in Oncology Chronic Myelogenous Leukemia. Version 2. Available at: http://www.nccn.org/professionals/physician_gls/f_guidelines.asp#site. Accessed January 2, 2014
8
Baccarani
 
M
Deininger
 
MW
Rosti
 
G
et al. 
European LeukemiaNet recommendations for the management of chronic myeloid leukemia: 2013.
Blood
2013
, vol. 
122
 
6
(pg. 
872
-
884
)
9
Yeung
 
DT
Osborn
 
MP
White
 
DL
et al. 
Australasian Leukaemia and Lymphoma Group
TIDEL-II: first-line use of imatinib in CML with early switch to nilotinib for failure to achieve time-dependent molecular targets.
Blood
2015
, vol. 
125
 
6
(pg. 
915
-
923
)
10
Kaplan
 
EL
Meier
 
P
Nonparametric estimation from incomplete observations.
J Am Stat Assoc
1958
, vol. 
53
 
282
(pg. 
457
-
481
)
11
Fine
 
JP
Gray
 
RJ
A proportional hazards model for the subdistribution of a competing risk.
J Am Stat Assoc
1999
, vol. 
94
 
446
(pg. 
496
-
509
)
12
R Core Team
 
R: a language and environment for statistical computing. Available at: http://www.R-project.org/. Accessed April 9, 2012
13
Gray
 
R
 
cmprsk: subdistribution analysis of competing risks. R package version 2.2-4. Available at: http://CRAN.R-project.org/package=cmprsk. Accessed March, 21, 2013
14
Therneau
 
T
 
A package for survival analysis in S. R package version 2.37-7. Available at: http://CRAN.R-project.org/package=survival. Accessed January 2, 2014
15
Gourraud
 
PA
Meenagh
 
A
Cambon-Thomsen
 
A
Middleton
 
D
Linkage disequilibrium organization of the human KIR superlocus: implications for KIR data analyses.
Immunogenetics
2010
, vol. 
62
 
11-12
(pg. 
729
-
740
)
16
Gonzalez-Galarza
 
FF
Christmas
 
S
Middleton
 
D
Jones
 
AR
Allele frequency net: a database and online repository for immune gene frequencies in worldwide populations.
Nucleic Acids Res
2011
, vol. 
39
 
Database issue
(pg. 
D913
-
D919
)
17
Gonzalez
 
A
Meenagh
 
A
Sleator
 
C
Middleton
 
D
Investigation of killer cell immunoglobulin-like receptor (KIR) gene diversity: KIR2DL2, KIR2DL5 and KIR2DS5.
Tissue Antigens
2008
, vol. 
72
 
1
(pg. 
11
-
20
)
18
Cisneros
 
E
Moraru
 
M
Gómez-Lozano
 
N
López-Botet
 
M
Vilches
 
C
KIR2DL5: an orphan inhibitory receptor displaying complex patterns of polymorphism and expression.
Front Immunol
2012
, vol. 
3
 pg. 
289
 
19
Chalandon
 
Y
Passweg
 
JR
Guglielmi
 
C
et al. 
Chronic Malignancies Working Party of the European Group for Blood and Marrow Transplantation (EBMT)
Early administration of donor lymphocyte infusions upon molecular relapse after allogeneic hematopoietic stem cell transplantation for chronic myeloid leukemia: a study by the Chronic Malignancies Working Party of the EBMT.
Haematologica
2014
, vol. 
99
 
9
(pg. 
1492
-
1498
)
20
Kreutzman
 
A
Rohon
 
P
Faber
 
E
et al. 
Chronic myeloid leukemia patients in prolonged remission following interferon-α monotherapy have distinct cytokine and oligoclonal lymphocyte profile.
PLoS One
2011
, vol. 
6
 
8
pg. 
e23022
 
21
Kreutzman
 
A
Juvonen
 
V
Kairisto
 
V
et al. 
Mono/oligoclonal T and NK cells are common in chronic myeloid leukemia patients at diagnosis and expand during dasatinib therapy.
Blood
2010
, vol. 
116
 
5
(pg. 
772
-
782
)
22
La Nasa
 
G
Caocci
 
G
Littera
 
R
et al. 
Homozygosity for killer immunoglobin-like receptor haplotype A predicts complete molecular response to treatment with tyrosine kinase inhibitors in chronic myeloid leukemia patients.
Exp Hematol
2013
, vol. 
41
 
5
(pg. 
424
-
431
)
23
Marin
 
D
Gabriel
 
IH
Ahmad
 
S
et al. 
KIR2DS1 genotype predicts for complete cytogenetic response and survival in newly diagnosed chronic myeloid leukemia patients treated with imatinib.
Leukemia
2012
, vol. 
26
 
2
(pg. 
296
-
302
)
24
Ali
 
S
Sergeant
 
R
O’Brien
 
SG
et al. 
Dasatinib may overcome the negative prognostic impact of KIR2DS1 in newly diagnosed patients with chronic myeloid leukemia.
Blood
2012
, vol. 
120
 
3
(pg. 
697
-
698
)
25
Kreutzman
 
A
Jaatinen
 
T
Greco
 
D
et al. 
 
Killer-cell immunoglobulin-like receptor gene profile predicts good molecular response to dasatinib therapy in chronic myeloid leukemia. Exp Hematol. 2012;40(11):906-913
Sign in via your Institution