Role of KIR and CD16A genotypes in colorectal carcinoma genetic risk and clinical stage

Background NK cell cytotoxicity is regulated by the types of the interaction between killer immunoglobulin-like receptors (KIRs) and human leukocyte antigen (HLA) class I ligands on target cells and the different binding affinity of the Fcγ receptor IIIA (CD16A) for IgG-coated tumor cells. Thus, it is conceivable that KIR and CD16A gene contents may contribute to the function of NK cells by modulating an immune response in the colorectal carcinoma (CRC) microenvironment. This hypothesis is supported by recent evidence suggesting that NK cells improve the clinical course of CRC patients by enhancing the anti-CRC effect of CD8 + T cells. This information provides the rationale to test the hypothesis whether the independent KIR segregation and specificity, as well as CD16A gene polymorphisms, have an impact on CRC. Methods Using polymerase chain reaction-sequence-specific primers (PCR-SSP) and sequence-based typing (SBT), we investigated KIR/HLA-C complex and CD16A (48H/R/L,158V/F) gene polymorphisms in 52 CRC patients and 61 local healthy controls (LCTRs). Results The allele frequency (AF) of at least five activating KIR (aKIRs) of the B haplotype (p = 0.036, OR 0.204), KIR2DL2 (p = 0.047, OR 0.2616), and KIR2DS2 genes (5.8 vs LCTR 13.8 % and vs. Fasano’s CTR 16.3 %, p = 0.05, OR 0.3145), in the absence of their cognate HLA-C1 ligands, were significantly associated with a reduced genetic risk of CRC. In contrast, CD16A-48H polymorphism was positively associated with an increased genetic risk of CRC (p = 0.05, OR 2.761). The latter was also found to be correlated with advanced stages of disease [III and IV (p = 0.03, OR 3.625)]. Conclusions Our data suggest that the analysis of aKIRs and KIR2DL2 gene and CD16A-48H may be of interest for the identification of individuals at reduced and increased genetic risk of CRC, respectively. Electronic supplementary material The online version of this article (doi:10.1186/s12967-016-1001-y) contains supplementary material, which is available to authorized users.


Background
Compelling evidence suggests that a pro-inflammatory solid tumor microenvironment composed of CD8 + T lymphocytes is associated with a favorable course of disease [1][2][3]. In contrast, a solid tumor microenvironment infiltrated by anti-inflammatory cells including tumorassociated macrophages (TAMs) [4] and regulatory T cells (Tregs) [5] is an unfavorable prognostic factor, and it is associated with poor prognosis. However, CRC seems to be an exception since inflammation, regardless the type of cellular infiltration including TAMs [6][7][8][9] and Tregs [10], is associated with a favorable prognosis. In this context, however, the role of NK cells is still debated [11,12]. The argument involving the role of NK cells in CRC is in part due to the narrowed degree of NK

Open Access
Journal of Translational Medicine *Correspondence: giuseppe.sconocchia@cnr.it 1 Laboratory of Tumor Immunology and Immunotherapy, CNR Institute of Translational Pharmacology (IFT), Via Fosso del Cavaliere 100, 00133 L' Aquila, Rome, Italy Full list of author information is available at the end of the article cell infiltration of solid tumors [12]. The reduced NK cell infiltration in the tumor bed can be, at least in part, explained by the detrimental effect of malignant cells on NK cells. Besides, there is also conflicting information about the clinical relevance of NK cells in CRC [13][14][15]. Recent findings, however, suggest that intra-tumor NK cell infiltration improves the clinical outcome of CRC by enhancing the protective role of tumor-infiltrating CD8 + T cells [16]. In addition, loss or down-regulation of HLA class I molecules on autologous CRC cells promotes NK cell cytotoxicity [17]. Because KIRs are independently segregated, an autologous KIR-HLA mismatched may exist. As a consequence, NK cells can trigger an enhanced cytotoxicity, as compared with NK cells expressing self-recognizing inhibitory KIRs [18].
KIRs consist of immunoglobulin-like receptors and heterodimers [19]. Immunoglobulin-like KIRs, with a long intracellular tail (KIRL), mediate inhibitory activity while those with a short intracellular tail (KIRS) can function as activating molecules. The HLA-C antigens are the best-known KIR ligands. Based on polymorphisms at positions 77 and 80, in the α1-domain of the β chain, KIRs can be identified into two additional groups: the HLA-C1 that carries an asparagine residue at position 80, and the HLA-C2 that has a lysine residue at position 80 [20,21]. To date, KIR genotyping is utilized for the analysis of gene content and categorization of A/B haplotypes as well as for prediction of NK cell reactivity in autologous and allo NK cell-based immunotherapy [22,23]. The A haplotype contains several inhibitory KIR genes and only one activating KIR. Conversely, B haplotype displays more activating KIR genes. Based on the linkage disequilibrium between particular alleles of different KIR loci, within the B haplotype, two gene clusters are identified. The first group consists of four centromeric KIR genes (C4 group), while the second group contains telomeric genes (T4 group) [24]. Recently, several studies have indicated a relationship between KIR genotype and cancer, but there is scant information about a link between KIR genes and CRC [25][26][27].
A primary mechanism by which NK cells kill CRC cells is the antibody-dependent cellular cytotoxicity (ADCC). This function is a result of the engagement of the Fc fragment of an antibody to the CD16A antigen on NK cell surface. The extent of ADCC relies on the binding affinity of the Fc fragment of the mAb to CD16A that is in turn regulated by CD16A polymorphisms [28]. In fact, the presence of valine/valine (V/V), but not valine/phenylalanine (V/F) at position 158, enhances ADCC and human IgG1 binding to CD16A [29]. Among monoclonal antibodies (mAbs), the anti-epidermal-growth factor receptor (EGFR), namely cetuximab has been successfully utilized for targeting therapy of EGFR+ malignancies including CRC [30]. Furthermore, cetuximab -activated NK cells have been shown to exert an indirect anti-tumor effect by promoting tumor antigen-anti-EGFR-specific T cells [31]. Thus, the genetic distribution of KIRs and CD16A is likely to have a substantial impact on NK cell phenotype and function that may affect the CRC pathogenesis. This information provides the rationale for testing the effect of KIR and CD16A genotypes on the genetic susceptibility to develop CRC and some clinical features.

Patient and control samples
The case-control study was performed using peripheral blood cells obtained from Italian CRC patients, recruited at the Policlinic of Tor Vergata between August 2010 and December 2012. This study (protocol number: 133/10) was approved by an independent, Institutional Review Board (IRB) Committee at the Teaching Hospital of Tor Vergata, Rome, Italy. Fifty-two CRC patients (31 males and 21 females), aged between 32 and 88 years, without a previous history of inflammatory bowel disease (IBD) and autoimmune diseases, were included in this protocol. Staging of the disease was assessed according to the pathological classification of tumor-node-metastasis (pTNM), as shown in Table 1. Control population comprised 61 healthy subjects (26 males and 35 females) originating from CRC patients geographical area, hereafter referred to as local controls (LCTRs) and historical controls [32,33].

CD16A genotyping by sequence-based typing
The genotyping of CD16A-158V/F and CD16A-48 leucine (L)/arginine (R)/histidine (H) was performed on genomic DNA by polymerase chain reaction (PCR) using primers previously described [36]. Briefly, PCR reactions were set up with 50 ng of genomic DNA per 50 µl reaction and amplified products were purified and sequenced using Big Dye Terminator Chemistry (Applied Biosystems, Foster City, CA) on an ABI Prism 3130 Genetic Analyzer (PE Applied Biosystems). Sequencing reactions were purified, and the analysis of results was obtained by alignment of the processed sequences and the reference human CD16A gene.

Statistical analysis
Pearson X 2 test was used for comparison of the gene frequency of the patient group with controls. Fisher's exact test was used when the expected difference between the two groups was small. The odds ratio (OR) and its 95 % confidence intervals (CI) were calculated. A p value of 0.05 or less was considered to be significant. The influence of KIR genes with or without association with their HLA-C ligands was assessed using the binary logistic regression analysis (IBM SPSS Statistics v. 19). The SPSS statistical package, Version 19 for Windows, was used for data and statistical analysis. CD16A haplotype frequencies were assessed by an expectation-maximum (EM) algorithm. Thus, samples containing one allele were considered homozygous, and that allele was counted twice in the analysis. For multilocus genotypic data, the maximum likelihood was estimated using an EM algorithm when the gametic phase was not known. To determine whether the observed diploid genotypes were the product of a random gamete union, the Hardy-Weinberg equilibrium was calculated by the Guo and Thomson exact test [37]. The analysis of molecular polymorphisms at CD16A-48 and 158, within the population, was performed using the ARLEQUIN V.3.0 POPULATION GENETICS software.

Association of five aKIRs of the B haplotype with a reduced susceptibility to CRC
The extent of KIR gene differences in 52 CRC patients and controls was assessed by SSP genotyping of 16 different NK receptor genes (data provided in Additional file 1: Table S1). To determine whether a distinct KIR haplotype frequency is associated with a reduced or increased risk for a genetic susceptibility to CRC, we compared the A and B haplotype frequencies of CRC patients with those of LCTRs. In both sample groups, the A/B haplotypes were the most represented (LCTR 70.5 % and CRC 73.1 %). Within the B haplotype, we considered different centromeric and telomeric genes, and we noted a reduced percentage of the C4T4 subtypes in CRC patients as compared to LCTR (5.8 vs. 16.4 %, Additional file 2: Table S2). More importantly, we found a significantly increased frequency of aKIR genes (≥5) in the B haplotype of LCTRs than that of CRC patients (16.4 vs. 3.9 %, p = 0.036, OR 0.204). These results were further validated when the frequency of aKIRs in the B haplotype of Italian historical controls [32] was compared to that of CRC patients (p = 0.0005), Table 2. These data suggest that the presence of ≥5 aKIR genes may protect healthy subjects from CRC.

Enhanced KIR2DL2 and KIR2DS2 allele frequencies, in the absence of their HLA-C ligand, is associated with a reduced susceptibility to CRC
Since the NK receptor/HLA ligand interaction controls NK-cell function influencing both the licensing and intensity of the NK cell immune response, we examined such a combination in CRC patients and controls by their cognate HLA-C ligand. We found that both frequencies of single HLA-C alleles and groups of KIRs, categorized according to the allelic residues at position 80 (Asn/Lys), were not significantly different between CRC patients and controls. Nevertheless, there was a statistically  [33]. These results suggest that the presence of KIR2DL2 and KIR2DS2, in the absence of their respective HLA C1 ligands, may protect one from the disease ( Table 3). The KIR2DS2 results were further supported by a binary logistic regression analysis (p = 0.008, OR 0.186).

The impact of CD16A-48 allele variants on CRC genetic susceptibility and CRC pathological staging
To determine the role of CD16A genotypes in the genetic susceptibility to CRC and pathological stage, we compared the 559G/T (158, V/F) and 230A/G/T (48, H/R/L) allele polymorphisms in CRC patients and LCTRs (Table 4). CD16A-48H was significantly associated with CRC (p = 0.05, OR 2.761) and stages III-IV (p = 0.03, OR 3.625). Considering CD16A allele polymorphisms in haplotype (

Discussion
Our study strengthens the protective role of NK cells in CRC that is supported by our KIR results in which a prevalence of the B haplotype hosting ≥5 aKIR genes, KIR2DS2 and KIR2DL2, in the absence of their HLA-C ligands, protects one from developing CRC. On the contrary, the contribution of NK cells to the CD16A-48H genotype-increased susceptibility to CRC and its association with advanced stages of disease is not readily apparent.
It is noteworthy that most of the studies addressing the role of NK cells in CRC have been performed utilizing immunohistochemistry. This methodology provides  In this study, we have compared the KIR genotype allele frequency of CRC patients to that of healthy local and historical controls, expecting to find a reduced allele frequency of activating KIRs and inhibitory KIRs, in the absence of their HLA-C ligands in CRC patients.
In support of our hypothesis, we found a significantly reduced frequency of ≥5 aKIRs, belonging to the B haplotype, in CRC patients when compared to LTCRs. In agreement with our study, other additional scientific reports, performed by KIR genotyping of CRC patients, have demonstrated an association of multiple aKIRs with long-term disease-free survival [38] and favorable responses to chemotherapy in patients with metastatic CRC [27]. An obvious interpretation of these data is that the presence of multiple KIRs may overcome the inhibitory effects of the KIRs on NK cells increasing the level of NK cell activation. On the other hand, it has been shown that NK cells are capable of eliminating autologous CRC   [39]. According to the missing self-hypothesis, an enhanced level of NK cell activation may also be applied to the interpretation of the protective effect of KIR2DL2 in healthy subjects, lacking their HLA-C1 ligand genes.
In contrast, the analysis of the significance of KIR2DS2 is harder to comprehend since it is not clear whether KIR2DS2 can activate NK cells in the absence of their HLA-C1 ligands.
The unique protective role of KIR2DL2 over the other KIRs may rely on the fact that KIR2DL2 is one of the most common KIR, nearly present in all individuals. Moreover, KIR2DL2 is also expressed by NK-T lymphocytes leading to an increased number of effector cells with anti-CRC potential [40]. The protective role of KIR2DL2, in the absence of the HLA-C1 ligands, is not restricted to CRC but is also demonstrated in inflammatory bowel disease, including Crohn's disease [41], underlying the importance of this KIR in NK cell immunity. In contrast, it is quite difficult to distinguish the role of KIR2DS2 from that of KIR2DL2. Considering the protective effect is achieved in the absence of HLA-C1 ligands, while an active KIR2DS2 needs to bind the HLA-C1 ligands, there is a chance that this effect may not be related to KIR2DS2 but is rather restricted to the influence of KIR2DL2 gene [42].
The complexity of this topic is underlined by the evidence that the presence of ≥5 aKIRs is associated with human papilloma virus-induced cervical neoplasia [43]. Furthermore, a second study has shown that the presence of an increased number of aKIR genes is also found to be associated with the Epstein-Barr virus-linked nasopharyngeal carcinoma (NPC) [44].
Thus, the presence of ≥5 activating KIRs has a different effect in solid tumors being protective in CRC and predisposing in cervical and nasopharyngeal carcinomas. One can speculate that this difference may reflect the possibility that NK cells are involved in the generation of two types of inflammation with opposite effects on tumor growth and development. In this context, NK cells may enhance the anti-CRC effect of CD8 + T cells by releasing interferon γ which up-regulates the expression of the HLA-class I and class II on the surface of CRC cells [45]. In contrast, NK cell activation, following virusinduced malignancies, such as EBV or HPV, may trigger mechanism(s) of tumor immunoevasion from the NK cell immune surveillance. On the other hand, it has been clearly shown that NK cells can damage dendritic cells and macrophages, both involved in the presentation of tumor-associated antigen(s) to T cells [46].
The extent of ADCC mediated by CD16A, in the presence of mAbs, depends on the cell types including NK cells, monocytes, and granulocytes and CD16A polymorphisms. This is because the latter possesses the ability to bind the Fc fragment of mAbs with different affinity, affecting the extent of ADCC against tumor cells [47,48]. Then, guided by recent findings indicating that some CD16A polymorphisms predict the therapeutic effectiveness of cetuximab in EGFR+ cancer patients [49], we found that only CD16A-48H gene variant is involved, being more associated with CRC patients and advanced stages of CRC. This association may reflect a status of immunodeficiency in CRC patients since the presence of CD16A-48H gene is linked to patients with episodes of recurrent viral infection of the respiratory tract supporting an in vivo NK cell abnormality [50]. By the available data, it is not possible to identify a cause/effect relationship between CD16A-48H phenotype and NK cell immunodeficiency. However, the presence of CD16A-48H results in a loss of the Leu11c/B73.1 epitope on NK cells that leads to the inhibition of the IgG binding to the CD16A [51,52]. In this light, we speculate that CD16A-CD48H may favor CRC development by an impaired ADCC against CRC cells.

Conclusions
Our study sustains the possibility that aKIR genotyping of the B haplotype, KIR2DL2 and HLA-C1 may be utilized to identify healthy individuals at reduced risk of developing CRC, while the genotyping of the CD16A-158V/48H haplotype is of interest for the identification of those at increased risk of CRC and tumor staging.