Acute lymphoblastic leukemia of pre-B cells (pre-B ALL) is the most frequent form of leukemia affecting children in Western countries. Evidence is accumulating that genetic factors play an important role in conferring susceptibility/resistance to leukemia in children. In this regard, activating killer-cell immunoglobulin-like receptor (KIR) genes are of particular interest. Humans may inherit different numbers of the 6 distinct activating KIR genes. Little is known about the impact of this genetic variation on the innate susceptibility or resistance of humans to the development of B-ALL. We addressed this issue by performing a case-control study in Canadian children of white origin. Our results show that harboring activating KIR genes is associated with reduced risk for developing B-ALL in these children. Of the 6 activating KIR genes, KIR2DS2 was maximally associated with decreased risk for the disease (P = 1.14 × 10−7). Furthermore, our results showed that inheritance of a higher number of activating KIR genes was associated with significant reductions in risk for ALL in children. These results were also consistent across different ALL phenotypes, which included children with pre-T cell ALL. Our study provides novel insights concerning the pathogenesis of childhood leukemia in white children and has implications for the development of new immunotherapies for this cancer.

Acute lymphoblastic leukemia (ALL) is the predominant form of leukemia occurring in children and accounts for approximately 75% of all leukemia cases. It is a heterogeneous group of white blood cell cancers that usually involve B- or T-cell precursors (blasts) and are accordingly called pre-B ALL (B-ALL) or pre-T ALL (T-ALL). In Western countries, B-ALL is the predominant form of ALL in children and comprises 70%-80% of the pediatric ALL cases, whereas T-ALL is a rare tumor in these populations and accounts for approximately 15% of the cases.1,2  The exact causes of leukemia are not known. Several environmental factors, including infections at an early age, exposure to ionizing radiation and certain chemicals (eg, benzene), and parental use of alcohol and tobacco, have been recognized as risk factors for ALL in children. Epidemiologic studies suggest that the genetic makeup of an individual plays an important role in determining his or her innate resistance or susceptibility to the development of this cancer. In this respect, we and others3-6  have shown that mutations in genes that compromise a person's ability to metabolize carcinogens, maintain DNA integrity, adequately respond to oxidative stress, and maintain cell-cycle check points increase the risk for developing ALL in children (reviewed in Stiller1 ). Recent genome-wide association studies have identified single-nucleotide polymorphisms in several new genes that could enhance susceptibility for childhood ALL. The genes include the AT-rich interactive domain (ARID)-5B, CCAAT/enhancer binding protein-ϵ (CEBPE), hydroxy acid oxidase (HAO)-1, erythrocyte membrane protein band 4.1-like (EPB4.1L)-2, and IKAROS family zinc finger protein (IKZF)-1.7-9 

Genome-wide association studies are very powerful and have the ability to identify key genetic determinants implicated in conferring susceptibility or resistance to human diseases, including B-ALL. An inherent limitation of these studies is that currently available chips cover only areas of the genome that have been well characterized for single-nucleotide polymorphisms and small insertion-deletions. Unfortunately, there are numerous regions of the genome that have not yet been well characterized, and variations therein have not been well represented in these chips. For such regions, a candidate gene approach would be necessary to identify potential susceptibility/resistance genes. One such region in the genome is the leukocyte receptor complex, which is located on human chromosome 19q13.4 and houses the killer-cell immunoglobulin-like receptor (KIR) genes (reviewed in Yawata et al,10  Parham,11  and Ianello et al12 ). Investigation of this region for potential genetic determinants for B-ALL is of great interest because some KIR genes have been shown to be associated with development of certain forms of leukemia in humans.13 

The KIR gene family comprises 14 functional genes. The genes encode synonymous receptors that are expressed mainly on natural killer (NK) cells and regulate their functional activities.12  NK cells have been well documented for their ability to kill cancer cells in the body, including leukemic ones.14,15  Of the KIR genes, 6 encode receptors that, on binding with their cognate ligands on target cells, activate NK cells and trigger NK cell–mediated killing and secretion of soluble mediators. Hence, these genes have been named as activating KIR genes and include KIR2DS1-3, KIR2DS4, and KIR3DS1.11,12  Although all KIR genes are highly polymorphic, an important variation with respect to the activating KIR genes is that humans may carry 1 or more of these genes. This variation in the inheritance of activating KIR genes results from differences in KIR haplotypes, which have been divided into 2 groups, A and B. The group A haplotypes are less variable and carry only 1 activating KIR gene, KIR2DS4, which often has a 22-bp deletion in exon 5 and encodes a soluble nonfunctional receptor.16  Thus, humans with the AA genotype may not carry any functional activating KIR gene. The group B haplotypes are more variable and may carry a full complement or a subset of the activating KIR genes. Some of these haplotypes lack KIR2DS4 but contain 1 or more of the other 5 activating KIR genes. Several studies have shown that KIR genotypes and the number of activating KIR genes inherited by a person may affect resistance/susceptibility to infectious agents, development of autoimmune diseases, reproductive efficiency, GVHD, and cancer17-22  (reviewed in Kulkarni et al23  and Farag et al24 ). As mentioned above, they have also been shown to protect humans from certain types of leukemia.13  However, little is known concerning the impact of activating KIR gene content of humans on their innate resistance/susceptibility to develop B-ALL. In the present study, we show that activating KIR genes are associated with decreased risk for development of B-ALL in children. Furthermore, we also show that the inheritance of a larger number of these genes in humans is associated with an enhanced reduction in the risk for developing this malignancy.

Patient populations

To investigate associations between activating KIR genes and B-ALL, we performed a case-control study at the Centre Hospitalier Universitaire Sainte-Justine Research Center in Montreal. Cases of B-ALL (n = 102) were recruited from the hematology/oncology clinic of the hospital during 1989-2004 and included children diagnosed before 18 years of age. Diagnosis of the leukemia was based on established criteria that included histologic examination of the blood and bone marrow smears, as well as immunophenotyping of the lymphoblasts. The relevant demographic, clinical, and histopathologic data were extracted from the medical records of these patients or acquired from the children's parents and are shown in Table 1.

Table 1

Demographic and clinical parameters of the patients

ParameterLeukemia type
B-ALLB-ALLT-ALL
Ethnicity* French Non-French Mixed 
Number 100 45 30 
Sex    
    Male 56 (56.0) 24 (53.3) 21 (70.0) 
    Female 54 (54.0) 21 (46.7) 9 (30.0) 
Immunophenotype    
    Pre-B 96 (96.0) 45 (100.0) 0 (0.0) 
    Pre-pre-B 4 (4.0) 0 (0.0) 0 (0.0) 
    Pre-T 0 (0.0) 0 (0.0) 30 (100.0) 
Age    
    ≤ 1 y 2 (2.0) 0 (0.0) 
    1-10 y 72 (72.0) 41 15 (50.0) 
    > 10 y 26 (26.0) 15 (50.0) 
Chromosomal translocations    
    None 67 (67.0) 10 (22.2) 10 (33.3) 
    t (12:21) 19 (19.0) 1 (2.2) 0 (0.0) 
    Others 8 (8.0) 1 (2.2) 2 (6.7) 
    ND 6 (6.2) 33 (73.4) 18 (60.0) 
Relapses    
    0 89 (89.0) ND ND 
    1 11 (11.0)   
DNA Index    
    1 63 (63.0) 12 (26.6) 2 (6.7) 
    > 1 35 (35.0) 12 (26.6) 23 (76.7) 
    ND 2 (2.0) 21 (46.7) 5 (16.7) 
ParameterLeukemia type
B-ALLB-ALLT-ALL
Ethnicity* French Non-French Mixed 
Number 100 45 30 
Sex    
    Male 56 (56.0) 24 (53.3) 21 (70.0) 
    Female 54 (54.0) 21 (46.7) 9 (30.0) 
Immunophenotype    
    Pre-B 96 (96.0) 45 (100.0) 0 (0.0) 
    Pre-pre-B 4 (4.0) 0 (0.0) 0 (0.0) 
    Pre-T 0 (0.0) 0 (0.0) 30 (100.0) 
Age    
    ≤ 1 y 2 (2.0) 0 (0.0) 
    1-10 y 72 (72.0) 41 15 (50.0) 
    > 10 y 26 (26.0) 15 (50.0) 
Chromosomal translocations    
    None 67 (67.0) 10 (22.2) 10 (33.3) 
    t (12:21) 19 (19.0) 1 (2.2) 0 (0.0) 
    Others 8 (8.0) 1 (2.2) 2 (6.7) 
    ND 6 (6.2) 33 (73.4) 18 (60.0) 
Relapses    
    0 89 (89.0) ND ND 
    1 11 (11.0)   
DNA Index    
    1 63 (63.0) 12 (26.6) 2 (6.7) 
    > 1 35 (35.0) 12 (26.6) 23 (76.7) 
    ND 2 (2.0) 21 (46.7) 5 (16.7) 

Values are n (%).

ND indicates not determined or unknown.

*

For ethnicity, French indicates whites of French-Canadian ancestry; Non-French, whites of non–French Canadian ancestry; and Mixed, white individuals of French or non–French Canadian ancestry.

Given the potential heterogeneity in the distribution of KIR gene frequencies within seemingly ethnically homogeneous populations, control subjects (n = 245) were selected from different sources to enhance population representation. These included children visiting the orthopedic department of the study hospital for minor fractures, their siblings, and randomly selected children and adults from the general population. The control subjects were restricted to people who were without any type of leukemia or any other cancer or autoimmune disorder. These control subjects have been used previously for replication or confirmation of associations between several susceptibility genes and complex diseases such as inflammatory bowel diseases.25-27  It is noteworthy that the KIR gene frequencies within the control population are reflective of the distribution within other populations of similar ancestry (Belgian and French; http://www.allelefrequencies.net; Table 2), which further demonstrates that bias, if any, in the selection of control subjects was limited. Both case subjects and control subjects were restricted to those with self-reported French Canadian ancestry and who resided in Quebec at the time of recruitment.

Table 2

Association between activating KIR genes and risk for B-ALL in white patients of French-Canadian ancestry

GenePatients (n = 100), n (%)Control subjects (n = 245), n (%)OR (95% CI)P
2DS1 28 (28.3)** 102 (41.6) 0.55 (0.33-0.92) .02 
2DS2 20 (20.0) 138 (56.3) 0.19 (0.11-0.34) 5.40 × 10−9 
2DS3 14 (14.1)** 83 (33.9) 0.32 (0.17-0.60) 3.70 × 10−4 
2DS4* 38 (39.2)** 139 (56.7) 0.49 (0.30-0.79) .004 
2DS5 19 (19.2)** 104 (42.4) 0.32 (0.18-0.56) 7.50 × 10−5 
3DS1 16 (16.0) 101 (41.2) 0.27 (0.15-0.49) 1.57 × 10−5 
GenePatients (n = 100), n (%)Control subjects (n = 245), n (%)OR (95% CI)P
2DS1 28 (28.3)** 102 (41.6) 0.55 (0.33-0.92) .02 
2DS2 20 (20.0) 138 (56.3) 0.19 (0.11-0.34) 5.40 × 10−9 
2DS3 14 (14.1)** 83 (33.9) 0.32 (0.17-0.60) 3.70 × 10−4 
2DS4* 38 (39.2)** 139 (56.7) 0.49 (0.30-0.79) .004 
2DS5 19 (19.2)** 104 (42.4) 0.32 (0.18-0.56) 7.50 × 10−5 
3DS1 16 (16.0) 101 (41.2) 0.27 (0.15-0.49) 1.57 × 10−5 

The prefix KIR is not shown with the gene names.

OR indicates odds ratio; and CI, confidence interval.

*

KIR2DS4 was considered absent when the PCR amplified only the 22-bp deleted amplicon. It was considered present when the nonmutant amplicon was observed on the gel with or without the mutant band.

**

Percentages are based on 99 patients, as these genes could not be genotyped in 1 patient.

To assess whether associations between the KIR genes and ALL differed according to the ALL phenotype and ethnicity, we further investigated these associations in a sample of 30 T-ALL pediatric patients of mixed white ancestry and in another sample of 40 pre-B ALL pediatric patients of non-French ancestry. The characteristics of these 2 cohorts are also given in Table 1.

Blood or saliva samples were collected from the participants as a source of DNA. Genomic DNA from these samples was extracted by a salting out technique or with commercial kits (Oragene, DNA Genotek). The quantity and quality of the extracted DNA samples were determined by UV spectrophotometry. The samples were coded and kept at the BioBank maintained at the Centre Hospitalier Universitaire Sainte-Justine Research Center.

The study was approved by the institutional ethics committee of the Centre Hospitalier (CHU) Ste-Justine, and written informed consent was acquired from all participants in accordance with the Declaration of Helsinki.

KIR genotyping

Genomic DNA samples from patients and control subjects were analyzed for the presence or absence of individual activating KIR genes by use of PCR and published sequence-specific primers as described previously.16,19,28-33  Standard PCR protocols were used. The general procedure included initial denaturation at 95°C for 5 minutes followed by 30-35 cycles, each comprising 30 seconds of denaturation at 95°C, 30 seconds of annealing at 60-68°C, and 2 minutes of extension at 72°C. The final cycle included extension for 10 minutes. The number of PCR cycles and annealing temperatures were adjusted for each gene in preliminary experiments. The presence or absence of the genes was determined by running the reaction products on 2%-2.5% agarose gels. The gels were stained with ethidium bromide, scanned, and imaged. The appearance of the amplicon band of the expected molecular weight was considered to indicate the presence of the gene in the genomic DNA sample. DNA samples known to be positive or negative for the gene were used as positive and negative controls, respectively. The PCR reactions were also run without the DNA template as a safeguard against false-positive cases, as well as for a housekeeping gene, GAPDH, as a control for false-negative cases. Despite all the precautions, the PCR reactions for 2 patients did not succeed, and these patients were excluded from the analyses.

Statistical analysis

On the basis of the genotyping results, patients were classified as those in whom the gene was present or absent. Because we could not distinguish between heterozygotes and homozygotes for a particular gene, it was not possible to examine deviation from Hardy-Weinberg equilibrium for the studied genes. Initial analysis included examination of the distribution of the genes in case and control subjects with χ2 tests. Subsequently, logistic regression analysis was performed. A single model for each gene as the independent variable was first fit. Because many of the KIR genes studied are in high linkage disequilibrium with each other, a multivariate model that included all 6 genes was then fit to assess the independent effect of each gene after controlling for the effects of other genes. In addition to the effect of individual genes, we also assessed the combined effects of harboring 1 or more activating genes. For this purpose, we stratified the samples according to the number of genes harbored (≤ 1, 2-3, and ≥ 4). Odds ratios and corresponding 95% confidence intervals were estimated. All analyses were performed with STATA software (Version 10; STATA Corp).

A total of 100 B-ALL cases of French Canadian ancestry and 245 control subjects were successfully genotyped for the 6 activating KIR genes. Table 2 shows results from univariate analysis for the association between the KIR genes and risk for ALL. The frequencies of all the activating KIR genes were lower in patients than in control subjects, and all of them significantly (P < .05) decreased the risk for developing B-ALL. Similar results were observed in the unrelated white cohort of B-ALL patients of non–French Canadian ancestry. As shown in Table 3, all of the genes showed significant associations in this sample except for the KIR2DS4 gene. A similar pattern of associations was evident for the unrelated sample of white T-ALL patients. As shown in Table 4, significant associations were observed with 3 genes (KIR2DS1, KIR2DS2, and KIR2DS3). In this smaller cohort, associations did not reach statistical significance (P > .05) in the case of the other 3 genes (KIR2DS4, KIR2DS5, and KIR3DS1), potentially because of low power. On the basis of the similarity of associations across the ALL phenotypes, we performed an analysis that combined all of the ALL patients. A multivariate logistic regression analysis of the combined cohort that accounted for the effects of other genes showed that each gene independently conferred protection for ALL (Table 5). The strongest associations were noted with the KIR2DS2 gene (odds ratio 0.27, 95% confidence interval 0.16-0.43, P = 1.14 × 10−7), with a similar pattern of association across each of the different ALL phenotypes.

Table 3

Association between activating KIR genes and risk for B-ALL in white patients of Canadian non-French ancestry

GenePatients (n = 45), n (%)Control subjects (n = 245), n (%)OR (95% CI)zP
2DS1 5 (11.1) 102 (41.6) 0.18 (0.07-0.46) −3.54 4.0 × 10−4 
2DS2 10 (22.2) 138 (56.3) 0.22 (0.10-0.47) −3.96 7.5 × 10−5 
2DS3 3 (6.8) 83 (33.9) 0.14 (0.04-0.46) −3.22 .001 
2DS4* 21 (46.7) 139 (56.7) 0.67 (0.35-1.26) −1.24 .21 
2DS5 7 (15.6) 104 (42.4) 0.25 (0.11-0.58) −3.22 .001 
3DS1 9 (20.0) 101 (41.2) 0.36 (0.16-0.77) −2.61 .009 
GenePatients (n = 45), n (%)Control subjects (n = 245), n (%)OR (95% CI)zP
2DS1 5 (11.1) 102 (41.6) 0.18 (0.07-0.46) −3.54 4.0 × 10−4 
2DS2 10 (22.2) 138 (56.3) 0.22 (0.10-0.47) −3.96 7.5 × 10−5 
2DS3 3 (6.8) 83 (33.9) 0.14 (0.04-0.46) −3.22 .001 
2DS4* 21 (46.7) 139 (56.7) 0.67 (0.35-1.26) −1.24 .21 
2DS5 7 (15.6) 104 (42.4) 0.25 (0.11-0.58) −3.22 .001 
3DS1 9 (20.0) 101 (41.2) 0.36 (0.16-0.77) −2.61 .009 

The prefix KIR is not shown with the gene names.

OR indicates odds ratio; and CI, confidence interval.

*

KIR2DS4 was considered absent when the PCR amplified only the 22-bp deleted amplicon. It was considered present when the nonmutant amplicon was observed on the gel with or without the mutant band.

Table 4

Association between activating KIR genes and risk for T-ALL in white

GenePatients (n = 30), n (%)Control subjects (n = 245), n (%)OR (95% CI)P
2DS1 6 (20.0) 102 (41.6) 0.35 (0.14-0.89) .027 
2DS2 9 (30.0) 138 (56.3) 0.32 (0.15-0.75) .009 
2DS3 2 (6.8) 83 (33.9) 0.14 (0.03-0.60) .008 
2DS4* 12 (40.0) 139 (56.7) 0.51 (0.23-1.10) .086 
2DS5 8 (26.7) 104 (42.4) 0.49 (0.21-1.15) .100 
3DS1 8 (26.7) 101 (41.2) 0.52 (0.22-1.21) .130 
GenePatients (n = 30), n (%)Control subjects (n = 245), n (%)OR (95% CI)P
2DS1 6 (20.0) 102 (41.6) 0.35 (0.14-0.89) .027 
2DS2 9 (30.0) 138 (56.3) 0.32 (0.15-0.75) .009 
2DS3 2 (6.8) 83 (33.9) 0.14 (0.03-0.60) .008 
2DS4* 12 (40.0) 139 (56.7) 0.51 (0.23-1.10) .086 
2DS5 8 (26.7) 104 (42.4) 0.49 (0.21-1.15) .100 
3DS1 8 (26.7) 101 (41.2) 0.52 (0.22-1.21) .130 

The prefix KIR is not shown with the gene names.

OR indicates odds ratio; and CI, confidence interval.

*

KIR2DS4 was considered absent when the PCR amplified only the 22-bp deleted amplicon. It was considered present when the nonmutant amplicon was observed on the gel with or without the mutant band.

Table 5

Association between activating KIR genes and ALL, including patients with B-ALL and T-ALL

GenePatients (n = 175), n (%)Control subjects (n = 245), n (%)Univariate analysis
Adjusted analysis*
OR (95% CI)POR (95% CI)P
2DS1 39 (22.4) 102 (41.6) 0.41 (0.26-0.63) 5.1 × 10−5 0.48 (0.29-0.80) .005 
2DS2 39 (22.3) 138 (56.3) 0.22 (0.14-0.34) 1.5 × 10−11 0.27 (0.16-0.44) 1.14 × 10−7 
2DS3 19 (10.9) 83 (33.9) 0.24 (0.14-0.41) 2.8 × 10−7 0.44 (0.24-0.80) .008 
2DS4 71 (41.3) 139 (56.7) 0.54 (0.36-0.79) .002 0.52 (0.33-0.82) .005 
2DS5 34 (19.5) 104 (42.4) 0.33 (0.21-0.52) 1.5 × 10−6 0.43 (0.26-0.72) .001 
3DS1 33 (18.7) 101 (41.2) 0.33 (0.21-0.52) 2.1 × 10−6 0.53 (0.32-0.89) .016 
GenePatients (n = 175), n (%)Control subjects (n = 245), n (%)Univariate analysis
Adjusted analysis*
OR (95% CI)POR (95% CI)P
2DS1 39 (22.4) 102 (41.6) 0.41 (0.26-0.63) 5.1 × 10−5 0.48 (0.29-0.80) .005 
2DS2 39 (22.3) 138 (56.3) 0.22 (0.14-0.34) 1.5 × 10−11 0.27 (0.16-0.44) 1.14 × 10−7 
2DS3 19 (10.9) 83 (33.9) 0.24 (0.14-0.41) 2.8 × 10−7 0.44 (0.24-0.80) .008 
2DS4 71 (41.3) 139 (56.7) 0.54 (0.36-0.79) .002 0.52 (0.33-0.82) .005 
2DS5 34 (19.5) 104 (42.4) 0.33 (0.21-0.52) 1.5 × 10−6 0.43 (0.26-0.72) .001 
3DS1 33 (18.7) 101 (41.2) 0.33 (0.21-0.52) 2.1 × 10−6 0.53 (0.32-0.89) .016 

The prefix KIR is not shown with the gene names.

OR indicates odds ratio; and CI, confidence interval.

*

Logistic regression analysis adjusted for the effect of other genes in the table.

KIR2DS4 was considered absent when the PCR amplified only the 22-bp deleted amplicon. It was considered present when the nonmutant amplicon was observed on the gel with or without the mutant band.

Because of haplotypic variations, humans may vary with respect to the number of inherited activating KIR genes, and therefore, we investigated the impact of the number of inherited activating KIR genes on the risk for developing this leukemia (including all of the ALL phenotypes). We observed that carrying a higher number of KIR genes was strongly associated with increasing protection for ALL (Table 6). Compared with children who harbored 1 or no genes, children who harbored 4 or more KIR genes were least susceptible to ALL (odds ratio 0.06, 95% confidence interval 0.03-0.13, P = 1.1 × 10−15).

Table 6

Association between the numbers of inherited activating KIR genes and risk for ALL (B-ALL and T-ALL)

Number of genesPatients, n (%)Control subjects, n (%)OR (95% CI)P
≤ 1 107 (61.1) 48 (19.6) Reference  
2-3 58 (33.1) 122 (49.8) 0.21 (0.13-0.34) 5.6 × 10−11 
≥ 4 10 (5.7) 75 (30.6) 0.06 (0.03-0.13) 1.1 × 10−15 
Number of genesPatients, n (%)Control subjects, n (%)OR (95% CI)P
≤ 1 107 (61.1) 48 (19.6) Reference  
2-3 58 (33.1) 122 (49.8) 0.21 (0.13-0.34) 5.6 × 10−11 
≥ 4 10 (5.7) 75 (30.6) 0.06 (0.03-0.13) 1.1 × 10−15 

OR indicates odds ratio; and CI, confidence interval.

In this case-control study, we investigated whether activating KIR genes were associated with susceptibility or resistance to B-ALL in Canadian children. We observed that all 6 activating KIR genes tended to reduce risk and confer protection against the acquisition of this type of leukemia after adjustment for the effects of individual genes (P values as low as 10−7). However, the decrease in risk was not the same with different genes. For example, the maximum reduction in the risk for developing B-ALL was conferred by inheritance of the KIR2DS2 gene. Furthermore, the present data strongly suggest that harboring increasing numbers of KIR genes provides additional protection against the acquisition of ALL. Furthermore, observed associations with the KIR genes extended to other phenotypes of ALL (T-ALL) and to white individuals of differing ancestry (whites of non-French ancestry).

To the best of our knowledge, this is first study that has described associations of activating KIR genes with reduction in risk for developing childhood ALL. A few earlier studies focused on adult patients. One group34,35  (reviewed in Verheyden et al15 ) has shown potential associations of certain inhibitory KIR-HLA gene interactions but did not demonstrate independent effects of activating KIR genes with chronic lymphoid and chronic myeloid leukemias in a Belgian population. No associations with ALL were noted in that study, probably because of low power (n = 8). Furthermore, the researchers did not discriminate between functional and mutant variants of the KIR2DS4 gene. As mentioned in the introduction, this gene is frequently mutated and shows a 22-bp deletion,16  which results in the expression of a nonfunctional receptor. In all the analyses performed in the present study, we considered those who showed only this 22-bp mutant form as being negative for the gene. Subjects were considered positive for the gene when a nonmutated amplicon band was seen on the gel with or without the presence of the mutated band. In another study conducted in Polish and German individuals, Giebel et al17  found a protective association between the KIR2DS4 gene and chronic myeloid leukemia, but no associations with ALL were evident. Similar to the previous study, the lack of association with ALL may have been the result of the small number of ALL patients (n = 21) investigated. A more recent study in a Chinese population13  has shown that activating KIR genes were more frequent in chronic myeloid leukemia patients than in healthy control subjects, although differences achieved statistical significance only for the KIR2DS4 gene. In contrast, that group found significant protective associations between the KIR2DS3 and ALL. The present findings related to the KIR2DS3 gene and ALL are consistent with that study. The lack of associations noted for the other KIR genes in the Chinese study may be the result of low power, age-of-onset effects, different types of leukemia, or population-specific differences in the frequencies of the KIR genes.

It is noteworthy that all the previous studies13,17,34,35  examined associations between KIR genes and adult-onset leukemia. Because the present study was focused on childhood-onset leukemia, direct comparisons between it and the earlier studies may not be valid and conclusive. This is particularly the case given that early-onset leukemia may involve different causes and hence differential susceptibilities vis-à-vis adult-onset leukemia. Furthermore, even within childhood leukemia, causes of ALL may differ between different types (B-ALL, T-ALL, and ALL with unspecified cell types). In the pediatric context, the present findings suggest that associations between the KIR genes may be common to at least 2 ALL phenotypes, B-ALL and T-ALL. Because of the lower prevalence of T-ALL, we could only manage to test a limited number (n = 30) of these cases in the present study. Despite their limited number, frequencies of all activating KIR genes were lower in the case subjects and reached significance (P < .05) for 3 of the genes (KIR2DS1, KIR2DS2, and KIR2DS3; Table 4). These results suggest that all of the activating KIR genes could have shown significant associations with this type of childhood leukemia had there been a larger number of case subjects. Certainly, additional studies will be required to confirm whether the findings observed here are unique to early-onset ALL. In the context of case-control studies performed for examining genetic associations, confounding because of population stratification (ethnicity) is a potential concern; however, samples in the present study were homogeneous with regard to ethnicity. Residual stratification cannot be ruled out; however, such stratification, if any, is unlikely to account for the strong associations noted in the study.

Ideally, we would have investigated the expression of ligands for activating KIR genes in the case and control subjects in the present study. Unfortunately, the ligands for these receptors remain largely unknown. Unlike their inhibitory counterparts, activating KIRs bind to HLA class I antigens with very low affinities.36  In fact, it is believed that these activating receptors may bind to the ligands that are expressed de novo in cancer or virus-infected cells, or they may recognize foreign peptides bound with HLA class I antigens. In the case of KIR2DS4, the receptor was described as binding an unidentified ligand expressed by primary melanoma cells.37  Another activating KIR gene, KIR3DS1, which was shown to be protective against B-ALL in the present study, has been associated with protection against infection with HIV type 1, as well as with delayed progression of the infection toward AIDS.18,22  In fact, KIR3DS1 is the collective name given to several allelic variants of KIR3DL1. These variants encode a short-tailed activating version of the receptor. KIR3DL1 is one of the most polymorphic genes in the KIR family. It exists in 59 allelic variants, of which 13 encode the short-tailed activating version of the receptor and hence are collectively designated as KIR3DS1, whereas the remaining 46 allelic variants encode a long-tailed inhibitory version of the receptor and are designated as KIR3DL1.10-12  Researchers have discovered that KIR3DL1 receptors bind HLA-Bw4 allotypes that have an isoleucine (Ile) at position 80 in their α-1 domain. The activating versions of the receptor were shown to bind the HLA-Bw4 allotypes with threonine (Thr) at position 80. However, subsequent studies have not confirmed these observations.10-12  Thus, like other short-cytoplasmic-tail activating KIR, the exact ligands for KIR3DS1 remain unknown. Furthermore, in the absence of antibodies (commercial or otherwise) that bind specifically with activating KIRs, functional studies on these receptors and their ligands must await advances in these areas.

It is noteworthy that in addition to NK cells, activating KIRs are also expressed on the surface of certain subsets of T cells. It is believed that activation of these receptors decreases the activation threshold of the immune cells and helps in controlling cancers and viral infections.23,24  The present findings related to childhood ALL fit with this paradigm and provide important novel insights concerning the immunopathogenesis of this malignancy. Because of their expression on the cell surface, KIRs are amenable to therapeutic interventions with antibodies or peptide mimics, unlike other target proteins that may be expressed inside cells. Furthermore, appropriately selected NK cells may be used as effector cells in future immunotherapeutic approaches against childhood leukemia.

In conclusion, the present study demonstrates novel associations between activating KIR genes and protection against childhood ALL.

An Inside Blood analysis of this article appears at the front of this issue.

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.

The authors thank the Canadian Cancer Society Research Institute for support.

Contribution: Z.A. and S.S. designed and performed the experiments and wrote the first version of the manuscript; A.I. and O.D. helped in data interpretation and in writing the manuscript; M.D., C.-I.R., and D.S. helped in patient selection; D.K.A. provided important insights on the disease, analyzed and interpreted the results, and contributed to the manuscript; A.A. conceived the idea, supervised the study, and contributed to the manuscript; and all authors read and contributed to the final version of the manuscript.

Conflict-of-interest disclosure: The authors declare no competing financial interests.

Correspondence: Ali Ahmad, Laboratory of Innate Immunity, CHU Sainte-Justine Research Center/Department of Microbiology & Immunology, 3175 Cote Ste Catherine, Montreal, QC, H3T 1C5, Canada; e-mail: ali.ahmad@recherche-ste-justine.qc.ca.

1
Stiller
 
CA
Epidemiology and genetics of childhood cancer.
Oncogene
2004
, vol. 
23
 
38
(pg. 
6429
-
6444
)
2
Kaatsch
 
P
Epidemiology of childhood cancer.
Cancer Treat Rev
2010
, vol. 
36
 
4
(pg. 
277
-
285
)
3
Sinnett
 
D
Krajinovic
 
M
Labuda
 
D
Genetic susceptibility to childhood acute lymphoblastic leukemia.
Leuk Lymphoma
2000
, vol. 
38
 
5–6
(pg. 
447
-
462
)
4
Krajinovic
 
M
Sinnett
 
H
Richer
 
C
Labuda
 
D
Sinnett
 
D
Role of NQO1, MPO and CYP2E1 genetic polymorphisms in the susceptibility to childhood acute lymphoblastic leukemia.
Int J Cancer
2002
, vol. 
97
 
2
(pg. 
230
-
236
)
5
Healy
 
J
Belanger
 
H
Beaulieu
 
P
Lariviere
 
M
Labuda
 
D
Sinnett
 
D
Promoter SNPs in G1/S checkpoint regulators and their impact on the susceptibility to childhood leukemia.
Blood
2007
, vol. 
109
 
2
(pg. 
683
-
692
)
6
Krajinovic
 
M
Lamothe
 
S
Labuda
 
D
et al. 
Role of MTHFR genetic polymorphisms in the susceptibility to childhood acute lymphoblastic leukemia.
Blood
2004
, vol. 
103
 
1
(pg. 
252
-
257
)
7
Trevino
 
LR
Yang
 
W
French
 
D
et al. 
Germline genomic variations associated with childhood acute lymphoblastic leukemia.
Nat Genet
2009
, vol. 
41
 
9
(pg. 
1001
-
1005
)
8
Papaemmanuil
 
E
Hosking
 
FJ
Vijayakrishnan
 
J
et al. 
Loci on 7p12.2, 10q21.2 and 14q11.2 are associated with risk of childhood acute lymphoblastic leukemia.
Nat Genet
2009
, vol. 
41
 
9
(pg. 
1006
-
1010
)
9
Han
 
S
Lee
 
K-M
Park
 
SK
et al. 
Genome wide association study of childhood acute lymphoblastic leukemia in Korea.
Leuk Res
2010
, vol. 
34
 
10
(pg. 
1271
-
1274
)
10
Yawata
 
M
Yawata
 
N
Abi-Rached
 
L
Parham
 
P
Variation within the human killer cell immunoglobulin-like receptor (KIR) gene family.
Crit Rev Immunol
2002
, vol. 
22
 
5–6
(pg. 
463
-
482
)
11
Parham
 
P
Immunogenetics of killer cell immunoglobulin-like receptors.
Mol Immunol
2005
, vol. 
42
 
4
(pg. 
459
-
462
)
12
Iannello
 
A
Debbeche
 
O
Samarani
 
S
Ahmad
 
A
Antiviral NK cell responses in HIV infection, I: NK cell receptor genes as determinants of resistance to HIV infection and progression to AIDS.
J Leukoc Biol
2008
, vol. 
84
 
1
(pg. 
1
-
26
)
13
Zhang
 
Y
Wang
 
B
Ye
 
S
et al. 
Killer cell immunoglobulin-like receptor gene polymorphisms in patients with leukemia: possible association with susceptibility to the disease.
Leukemia Res
2010
, vol. 
34
 
1
(pg. 
55
-
58
)
14
Ruggeri
 
L
Capanni
 
M
Urbani
 
E
et al. 
Effectiveness of donor natural killer cell alloreactivity in mismatched hematopoietic transplants.
Science
2002
, vol. 
295
 
5562
(pg. 
2097
-
2100
)
15
Verheyden
 
S
Demanet
 
C
NK cell receptors and their ligands in leukemia.
Leukemia
2008
, vol. 
22
 
2
(pg. 
249
-
257
)
16
Maxwell
 
LD
Williams
 
F
Gilmore
 
P
Meenagh
 
A
Middleton
 
D
Investigation of killer-cell immunoglobulin-like receptor gene diversity, II: KIR2DS4.
Hum Immunol
2004
, vol. 
65
 
6
(pg. 
613
-
621
)
17
Giebel
 
S
Nowak
 
I
Krawczyk-Kulis
 
M
Holowiecki
 
H
Kyrcz-Krzemien Kursnierczyk
 
P
Association of KIR2DS4 and its variant KIR1D with leukemia.
Leukemia
2008
, vol. 
22
 
11
(pg. 
2129
-
2130
)
18
Martin
 
MP
Gao
 
X
Lee
 
JH
et al. 
Epistatic interaction between KIR3DS1 and HLA-B delays progression to AIDS.
Nat Genet
2002
, vol. 
31
 
4
(pg. 
429
-
434
)
19
Hiby
 
SE
Walker
 
JJ
O'Shaughnessy
 
KM
et al. 
Combinations of maternal KIR and foetal HLA-C genes influence the risk of preeclampsia and reproductive success.
J Exp Med
2004
, vol. 
200
 
8
(pg. 
957
-
965
)
20
Ruggeri
 
L
Capanni
 
M
Urbani
 
E
et al. 
Effectiveness of donor natural killer cell alloreactivity in mismatched hematopoietic transplants.
Science
2002
, vol. 
295
 
5562
(pg. 
2097
-
2100
)
21
Ahlenstiel
 
G
Martin
 
MP
Gao
 
X
Carrington
 
M
Rehermann
 
B
Distinct KIR/HLA compound genotypes affect the kinetics of human antiviral natural killer cell responses.
J Clin Invest
2008
, vol. 
118
 
3
(pg. 
1017
-
1026
)
22
Boulet
 
S
Sharafi
 
S
Simic
 
N
et al. 
Increased proportion of KIR3DS1 homozygotes in HIV-exposed uninfected individuals.
AIDS
2008
, vol. 
22
 
5
(pg. 
595
-
599
)
23
Kulkarni
 
S
Martin
 
MP
Carrington
 
M
The Yin and Yang of HLA and KIR in human disease.
Semin Immunol
2008
, vol. 
20
 
20
(pg. 
343
-
352
)
24
Farag
 
SS
Fehniger
 
TA
Ruggeri
 
L
Velardi
 
A
Caligiuri
 
MA
Natural killer cell receptors: new biology and insights into the graft-versus-leukemia effect.
Blood
2002
, vol. 
100
 
6
(pg. 
1935
-
1947
)
25
Amre
 
DK
Mack
 
DR
Morgan
 
K
et al. 
Autophagy gene ATG16L1 but not IRGM is associated with Crohn's disease in Canadian children.
Inflamm Bowel Dis
2009
, vol. 
15
 
4
(pg. 
501
-
507
)
26
Amre
 
DK
Mack
 
D
Israel
 
D
et al. 
Association between genetic variants in the IL-23R gene and early-onset Crohn's disease: results from a case-control and family-based study among Canadian children.
Am J Gastroenterol
2008
, vol. 
103
 
3
(pg. 
615
-
620
)
27
Amre
 
DK
Mack
 
DR
Morgan
 
K
et al. 
Investigation of reported associations between the 20a13 and 21q22 loci and pediatric-onset Crohn's disease in Canadian children.
Am J Gastroenterol
2009
, vol. 
104
 
11
(pg. 
2824
-
2828
)
28
Uhrberg
 
M
Parham
 
P
Wernet
 
P
Definition of gene content for nine common group B haplotypes of the Caucasoid population: KIR haplotypes contain between seven and eleven genes.
Immunogenetics
2002
, vol. 
54
 
4
(pg. 
221
-
229
)
29
Uhrberg
 
M
Valiante
 
NM
Shum
 
BP
et al. 
Human diversity in killer cell inhibitory receptor genes.
Immunity
1997
, vol. 
7
 
6
(pg. 
753
-
763
)
30
Vilches
 
C
Rajalingam
 
R
Uhrberg
 
M
Gardiner
 
C
Young
 
NT
Parham
 
P
KIR2DL5, a novel killer-cell receptor with a D0-D1 configuration of Ig-like domains.
J Immunol
2000
, vol. 
164
 
11
(pg. 
5797
-
5804
)
31
Yawata
 
M
Yawata
 
N
McQueen
 
KL
et al. 
Predominance of group A haplotypes in Japanese associated with diverse NK cell repertoires of KIR expression.
Immunogenetics
2002
, vol. 
54
 
8
(pg. 
543
-
550
)
32
Zhi-ming
 
L
Yu-lian
 
J
Zhao-li
 
F
et al. 
Polymorphisms of killer-cell immunoglobulin-like receptor gene: possible association with susceptibility to or clearance of hepatitis B virus infection in Chinese Han population.
Croat Med J
2007
, vol. 
48
 
6
(pg. 
800
-
806
)
33
Gomez-Lozano
 
N
Vilches
 
C
Genotyping of human killer-cell immunoglobulin-like receptor genes by polymerase chain reaction with sequence-specific primers: an update.
Tissue Antigens
2002
, vol. 
59
 
3
(pg. 
184
-
193
)
34
Verheyden
 
S
Demanet
 
C
Susceptibility to myeloid and lymphoid leukemia is mediated by distinct inhibitory KIR-HLA ligand interactions.
Leukemia
2006
, vol. 
20
 
8
(pg. 
1437
-
1438
)
35
Verheyden
 
S
Bernier
 
M
Demanet
 
C
Identification of natural killer cell receptor phenotypes associated with leukemia.
Leukemia
2004
, vol. 
18
 
12
(pg. 
2002
-
2007
)
36
Katz
 
G
Markel
 
G
Mizrahi
 
S
Arnon
 
TI
Mandelboim
 
O
Recognition of HLA-Cw4 but not HLA-Cw6 by the NK cell receptor killer cell Ig-like receptor two-domain short-tail number 4.
J Immunol
2001
, vol. 
166
 (pg. 
7260
-
7267
)
37
Katz
 
G
Gazit
 
R
Arnon
 
TI
et al. 
MHC class I-independent recognition of NK-activating receptor KIR2DS4.
J Immunol
2004
, vol. 
173
 
3
(pg. 
1819
-
1825
)

Author notes

*

Z.A. and S.S. contributed equally to this study.

Sign in via your Institution