Correction
15 May 2015: De Re V, Caggiari L, De Zorzi M, Repetto O, Zignego AL, et al. (2015) Correction: Genetic Diversity of the KIR/HLA System and Susceptibility to Hepatitis C Virus-Related Diseases. PLOS ONE 10(5): e0128849. https://doi.org/10.1371/journal.pone.0128849 View correction
Figures
Abstract
Background
The variability in the association of host innate immune response to Hepatitis C virus (HCV) infection requires ruling out the possible role of host KIR and HLA genotypes in HCV-related disorders: therefore, we therefore explored the relationships between KIR/HLA genotypes and chronic HCV infection (CHC) as they relate to the risk of HCV-related hepatocarcinoma (HCC) or lymphoproliferative disease progression.
Methods and Findings
We analyzed data from 396 HCV-positive patients with CHC (n = 125), HCC (118), and lymphoproliferative diseases (153), and 501 HCV-negative patients. All were HIV and HBV negative. KIR-SSO was used to determine the KIR typing. KIR2DL5 and KIR2DS4 variants were performed using PCR and GeneScan analysis. HLA/class-I genotyping was performed using PCR-sequence-based typing. The interaction between the KIR gene and ligand HLA molecules was investigated. Differences in frequencies were estimated using Fisher’s exact test, and Cochran-Armitage trend test. The non-random association of KIR alleles was estimated using the linkage disequilibrium test. We found an association of KIR2DS2/KIR2DL2 genes, with the HCV-related lymphoproliferative disorders. Furthermore, individuals with a HLA-Bw6 KIR3DL1+ combination of genes showed higher risk of developing lymphoma than cryoglobulinemia. KIR2DS3 gene was found to be the principal gene associated with chronic HCV infection, while a reduction of HLA-Bw4 + KIR3DS1+ was associated with an increased risk of developing HCC.
Conclusions
Our data highlight a role of the innate-system in developing HCV-related disorders and specifically KIR2DS3 and KIR2D genes demonstrated an ability to direct HCV disease progression, and mainly towards lymphoproliferative disorders. Moreover the determination of KIR3D/HLA combination of genes direct the HCV progression towards a lymphoma rather than an hepatic disease. In this contest IFN-α therapy, a standard therapy for HCV-infection and lymphoproliferative diseases, known to be able to transiently enhance the cytotoxicity of NK-cells support the role of NK cells to counterstain HCV-related and lymphoproliferative diseases.
Citation: De Re V, Caggiari L, De Zorzi M, Repetto O, Zignego AL, Izzo F, et al. (2015) Genetic Diversity of the KIR/HLA System and Susceptibility to Hepatitis C Virus-Related Diseases. PLoS ONE 10(2): e0117420. https://doi.org/10.1371/journal.pone.0117420
Academic Editor: Golo Ahlenstiel, University of Sydney, AUSTRALIA
Received: July 28, 2014; Accepted: December 21, 2014; Published: February 20, 2015
Copyright: © 2015 De Re et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited
Data Availability: All relevant data are within the paper.
Funding: Italian Association for Cancer Research (AIRC n.10266) grant funding to part of the research reagents. Laura caggiari, Mariangela De Zorzi and Ombretta Repetto fellowships are funded by 5‰ Dir Sc. CRO. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
The World Health Organization (WHO) estimates that about 3% of the world population is infected with the Hepatitis C virus (HCV), and three to four million individuals are newly infected each year. Although new antiviral treatments are very promising [1], today only a minority of patients successfully clear up HCV infections, and the remaining patients (60–85%) develop chronic infection.
The patients with persistent infection are at risk of developing chronic liver lesions, ranging from minimal inflammation to cirrhosis and hepatocellular carcinoma (HCC). Moreover, chronic infections could cause several extrahepatic diseases, including the type II mixed cryoglobulinemia (MC), a systemic autoimmune disease characterized by a monoclonal/oligoclonal proliferation of B-cells. This produces immune complex-mediated disorders and may evolve in a small fraction of patients to frank malignant B-cell non-Hodgkin’s lymphoma (NHL) [2].
The mechanisms whereby HCV establishes chronic infection, autoimmune diseases, HCC or NHL are still poorly understood. Viral factors and determinants of the host are both involved in HCV clearance and disease pathogenesis. Spontaneous viral clearance has been associated with a low viral diversity and strong CD8+ and CD4+ T-cell responses targeting a broad range of viral epitopes [3]. At the same time, infectious agents are known to develop a high variety of mechanisms, which finally lead to an ineffectual T-cell response and a persistent infection. Currently, the exact mechanism responsible for the T-cell failure of HCV infection is not well-defined, although several evidences suggest that HCV mutations within human leukocyte antigen (HLA)-restricted epitopes, which are no longer recognized by T-cells and neutralizing antibodies, are one of the most potent strategies utilized by HCV. Moreover, other mechanisms potentially contributing to virus persistence include the HCV interference with host cellular components and signaling pathways [4]. Additionally, alcohol consumption, age, sex and genetics cofactors are also important predisposing conditions.
Natural killer (NK) cells are part of the innate immune system that plays an important role in HCV infection by killing infected/altered cells through a direct and antibody-dependent cell-mediated (ADCC) cytotoxic action and through the release of cytokines that activate and go back other effector cells of the immune system [5;6].
NK cell activity is mainly regulated by NK receptors, which include the killer cell Ig like receptors (KIR), encoded by a family of activating and inhibitory genes located on the human chromosome 19. At present, a total of 16 KIR human genes have been characterized, of which 8 are NK-cell inhibitors (KIR2DL1–5, 3DL1–3), 6 are activators (KIR2DS1–5, 3DS1) and 2 are pseudo genes (KIR2DP1, 3DP1). Contrary to T-cell, which recognize tumor antigens presented by HLA molecules, NK cells can also recognize the “loss” of HLA induced by transformation, cellular stress or infection. Different levels of NK-cell activation are modulated by various KIR interactions with their cognate HLA ligand, resulting in an overall balance of signals from activating versus inhibitory receptors. Depending on both the inheritance of stochastic KIR and HLA genes, as well as levels of KIR/HLA surface expression levels, a cell tolerance or rejection occurs [7].
NK cells have also been reported to play a role in HCV clearance [5]. In particular, NK cells were demonstrated to mediate the inhibition of HCV-replication and to exert a targeted cytotoxic action against targeted cells given that NK cells isolated from healthy donors kill HCV-replicating cells and secrete IFN-γ [8;9]. However, there are evidences that HCV-mediated interferences with the action of NK cells [10;11].
HLA-C1 KIR2DL3+ in homozygosis has been associated with HCV clearance in several studies [12;13] but the occurrence of this association was not always observed. At the same time, HLA-C1 KIR2DL3+ has also been associated with sustained virologic response to anti-HCV therapy [10;11]. Moreover, a protective role for HLA-Bw4 KIR3DS1+ against liver disease progression has been proposed [14;15].
Liver NK cells are abundant in the liver and play beneficial roles in inhibiting viral infection, tumor cell growth and liver fibrosis, but they can also play detrimental roles in stimulating liver injury and attenuating liver regeneration by directly killing target cells and the production of IFN-γ. Conversely, inhibition of NK cells may play a defensive role by protecting liver cells from inflammation and injury, as supported by some experimental models [16–18]. Both the frequency of liver and peripheral NK cells and their function were found to decrease in HCC patients, suggesting a potential protective role, although many aspects of their exact involvement in HCC still needs to be deciphered [5;19–22]. The purpose of this study is to examine KIR genotypes and HLA KIR compounds in chronic HCV-infected patients with HCC and malignant extrahepatic lymphoproliferative diseases in order to investigate the potential effect of the host HLA/KIR profile on HCV-related disease progression.
Methods
Ethics Statement
CRO Aviano National Cancer Institute Review Board (IRB) have approved the study. Ethical principles espresse in the Declaration of Helsinki has been respected. Informed written or oral consent have been obtained from the participants as reported in the Methods section of the manuscript
Patient characteristics
A total of 396 HCV-infected (HCV-positive) and 501 non HCV-infected (HCV-negative) Italian patients were selected for the present study. Patients were recruited from 9 Italian centers (CRO Aviano National Cancer Institute, Pordenone, “S. Maria degli Angeli” General Hospital, Pordenone; University of Bari Medical School, Bari; “S. Maria della Misericordia” General hospital, Udine; University of Florence, Florence; “IRCCS Casa Sollievo della Sofferenza,” San Giovanni Rotondo, Bari; National Cancer Institute “Fondazione Pascale,” Naples; Ospedale “S. Gerardo di Monza,” Milan; University of Catania, Catania).
The diagnosis for chronic HCV infection is based essentially on the positivity of antibodies against HCV and serum transaminase (ALT) levels. Test immunoenzimatico (III-generation EIA) against HCV-core and HCV-non structural antigens were used. There are several commercial tests for the determination of the genotypes: the mostly used was the test Inno-Lipa. As regard the diagnosis of HCC, clinical centers used the standard criteria listed in the European Association for the Study of the Liver (EASL) that incorporate both invasive and noninvasive measures [23] Noninvasive criteria include two imaging techniques, both demonstrating a focal lesion >2 cm in diameter with features of arterial hypervascularization. Most coauthors are members of the A.L.CRI (Associazione Italiana per la Lotta alle Crioglobulinemie) an association that proposed and validate the consensus protocol for detection and typing of cryoglobulins [24]. NHL in the course of HCV infection has been confirmed by histopathologists based on WHO classification. [25].
Both groups of patients were tested negative for HBV and HIV infections. Among HCV-positive patients, 125 had a chronic HCV infection but without HCC or any sign/symptom of definite MC or NHL (CHC n = 125), patients with HCV-related hepatocellular carcinoma (HCC n = 118), and patients with a lymphoproliferative disease (n = 153) including a group affected by either a definite cryoglobulinemic syndrome according to previously described criteria [23] (MC n = 75) or a definite malignant B-cell non-Hodgkin’s lymphoma (NHL n = 78). Our reference series of HCV negative cohort members included 501 individuals with negative virological tests for HCV, HBV or HIV with no clinical evidence of neoplastic or autoimmune disorders (n = 79), with a celiac disorder (n = 76) or with a gastric /intestinal neoplasia (n = 346). Informations concerning the demographic and clinical characteristics of the patient’s groups were reported on Table 1.
Genomic DNA was extracted from peripheral blood sample using the EZ1 DNA blood kit and the BioRobot EZ1 Workstation (Qiagen Inc., Valencia, CA).
Samples were collected after obtaining informed written consent.
.
PCR-SSP KIR typing
All the KIR/HLA typing has been performed in a unique center (CRO). Genomic DNA was used to determine the genotype of KIR genes using Lifecodes KIR-SSO typing kit for use with Luminex (Gen-Probe Transplant Diagnostic, GTI for Italy) according to the manufacturer’s instructions. The presence or absence of the following 16 KIR genes was identified: KIR2DL1, KIR2DL2, KIR2DL3, KIR2DL4, KIR2DL5, KIR2DS1, KIR2DS2, KIR2DS3, KIR2DS4, KIR2DS5, KIR3DL1, KIR3DL2, KIR3DL3, KIR3DS1, KIR2DP1 and KIR3DP1.
KIR gene profiles were determined by the presence or absence of each KIR gene in a given individual. This method of KIR typing does not allow the direct determination of KIR2DL2 copy number. Instead, we used the allelic nature of KIR2DL2 and KIR2DL3 at the 2DL2/2DL3 locus to infer the number of copies of KIR2DL2.
KIR2DL5 A and B subtyping was performed using polymerase chain reaction (PCR), as we previously reported [24]. KIR2DS4 was typed for encoded cell-surface receptor (full) or a truncated protein variant with loss of the transmembrane and cytoplasmatic domains (del) by GeneScan analysis with a 6-FAM labeled primer [24]. The deleted variant of KIR2DS4 was not anchored to the cell membrane but was encoded for a soluble form of the protein that is potentially secreted and likely lacks function(s).
KIR genotype assessment
KIR genotype was assigned according to the database (http://www.allelefrequencies.net/). All genotype contained KIR2DL4, KIR3DL2, KIR3DL3, KIR2DP1 and KIR3DP1 framework genes. In addition, genotypes were assumed to contain either 2DL5A or 2DL5B, and 2DS4full or 2DS4del.
In the assessment of the KIR genotype, group B genotypes were defined by the presence of one or more of the following genes: KIR2DL5, KIR2DS1, KIR2DS2, KIR2DS3, KIR2DS5 and KIR3DS1. Conversely, stable group A genotype was defined by the absence of all these genes and by the presence of KIR3DL1, KIR2DL1, KIR2DL3 and KIR2DS4 genes. A comparison was done with the frequencies reported in the database for the Caucasian population [25] and linkage disequilibrium organization of the complex human KIR superlocus [26].
Moreover, centromeric (Cent) and telomeric (Tel) regions splits in half the KIR genotype; KIR3DL3 and KIR3DP1 delimited the centromeric part of the KIR locus, whereas KIR2DL4 and KIR3DL2 delimited the telomeric part (Fig. 1). KIR2DL5, KIR2DS5 and KIR2DS3 genes can be present both in centromeric and telomeric locations.
A total of 42 KIR gene profiles were identified (panel A). Genotype ID number reported are those from the reference database [28]. The presence of KIR genes is indicated by grey shanding. Genotypes AA and BX according to criteria reported in material and method section are indicated in the first column.°Only ID Genotypes present in at least two patients were reported in tables. Distribution of genotype ID1 and ID4 for HCV-negative patients and for patients with HCV infection and different outcomes (panel B). Distribution of KIR3DL1, KIR2DL3, KIR2DS4, KIR2DS3 genes ID4 for HCV-negative patients and for patients with HCV infection and different outcomes (panel C). * Fisher’s Exact test, p<0.05.
SBT HLA typing
HLA genotyping was performed using PCR-sequence-based typing (PCR-SBT) with primers specific for each class I locus (A, B and C) [24].
Two novel HLA class I alleles (A*02:374 and B*51:141) were identified in this study and submitted in the GenBank nucleotide sequence database.
KIR/HLA interactions
Several interactions between KIR genes and HLA ligands were contemplated. HLA-C molecules, the KIR-ligands, were classified as either HLA-C1 or HLA-C2 on the basis of dimorphisms at position 80. So far, known KIR-HLA pairs are HLA-C2 epitope (Asparagine at position 77, Lysine at position 80) 2DL1/2DS1+ and HLA-C1 epitope (Serine at position 77, Asparagine at position 80) 2DL2/3/2DS2+, HLA-Bw4 3DL1/3DS1+, HLA-C*04 2DS4+ and HLA-A3/A11 3DL2+.
HLA-Bw4–80I (Isoleucine at position 80) was recognized to bind KIR3DL1 with higher affinity than Bw4–80T (Treonine at position 80), so we also analyzed the frequencies of Bw4-containing allotypes based on their dimorphism at position 80 [30].
Ligands for the lasting KIRs remain elusive.
Statistical analysis
The genotypes of KIR gene and KIR-HLA pair frequencies were individually determined by direct counting of the individual who tested positive for a specific (pair of) gene. Linkage disequilibrium (LD) between pairwise KIR genes was assessed with Cramer’s V correlation coefficient. Differences in frequencies were estimated by using Fisher’s exact test and trend in frequencies were analyzed with Cochran-Armitage trend test (z2). A p-value of <0.05 was considered significant.
Results
KIR genotypes and single KIR genes
In this study, we found 65 different KIR genotype profiles: the 38 that were present in at least two patients are reported in Fig. 1A. The framework of KIR2DL4, KIR3DL2, KIR3DL3 and KIR3DP1 genes were present in all individuals; thus, they are not reported in the Figure. Two new genotypes, indicated as “new” in the identification number column (ID), were found and reported in the S1 Table.
The genotype ID 4 frequency in patients with CHC was found to be lower than in HCV-positive patients with lymphoproliferative disorder (8% versus 20.9%, Fisher’s exact test, p<0.01, Fig. 1B). Furthermore, the genotype ID 4 frequency in patients with lymphoproliferative disorder was higher than in HCV-negative patients (20.9% versus 12.4%, p = 0.01, Fig. 1B) and higher than that found in the reference Caucasoid ethnicity (ranging between 18.0% and 6.0% in Portugal and in Iran East Azezebaijan, respectively) [25] suggesting that the genotype ID 4 is in relation with HCV-related lymphoproliferative disorders, particularly NHL.
Classic linkage disequilibrium (LD) studies identified two distinct regions in the KIR cluster around the KIR2DL4 gene: a centromeric (cen) region, which seems to be driven by the KIR2DL5 and KIR2DL2/3 loci, and a telomeric (tel) region driven by the KIR3DL1/S1 locus. The overall patterns of LD of KIR receptors were similar in all groups of patients examined in this study (Fig. 2). The genotype ID4 found to be associated with HCV-lymphoproliferative disorders is distinguished from the most prevalent ID 1 for the presence of the KIR2DS2 in strong LD with the KIR2DL2 gene.
A stretch of 14 kb DNA that interconnects KIR3DP1 and KIR2DL4 divides the KIR genotype into two halves. The centromeric half is delimited by 3DL3 and 3DP1, while the telomeric half is delimited by 2DL4 and 3DL2. There is different KIR gene content, due to a recombination of these genes, in KIR genotypes across individuals and populations. The framework genes, present in all genotypes are shown in grey boxes; genes encoding activating KIR are in red color; and those for inhibitory receptors are in blue color. KIR2DL4 encodes a receptor that has both inhibitory and activating functions The KIR2DP1 and 3DP1 (green) are pseudogenes that do not express a receptor. Pairwise D’ LD based on Cramer’s V correlation coefficient between the presence and absence of different KIR genes in four groups of patients (panel B-E) B: HCV negative; C: Chronic HCV; D: Hepatocellular carcinoma; E: Lymphoproliferative disease. The KIR cluster genetic polymorphism is considered as the presence or absence of KIR genes.
When the frequency of single KIR genes was compared a number of statistically significant differences in HCV-positive groups (Fig. 1C) was identified. Among these genes, the frequency of KIR2DS3 gene was lower in HCC malignancies than in CHC (30.5% in HCC versus 43.2% in CHC, p = 0.05, Fig. 1C). Moreover, although the differences in HCV-negative groups did not reach a statistical significance, the frequency of this gene was found to be lower in HCV-negative individuals (36.1%) than in CHC patients (43.2%). The possibility of inheriting both an absence of the KIR2DS3 gene and KIR2DS2/KIR2DL2 genes is about 47–48% in the worldwide population, and this data underlines the specific KIR2DS3 association with CHC [28].
To contrast the findings for the KIR2DS3 that indicated the frequency was higher in patients with CHC, both KIR3DL1 and KIR2DS4 gene frequencies were lower in CHC than in HCC cases (Fig. 1C). Both of these genes are in strong LD (Fig. 2), and according to LD they are rarely inherited (≤25%) along with the KIR2DS3 gene.
Centromeric and Telomeric KIR locus
Subdividing the KIR cluster into the cent and tel regions identified by the LD study, simplifies the description of the haplotype structure of the KIR genes. This case series (Table 2) found 9 different centromeric gene motifs (Cent 1–9), 8 telomeric gene motifs (Tel 1–8) and 6 centromeric/telomeric gene (Cent/Tel 1–6) motifs.
With regards to the centromeric KIR regions, the Cent 2 gene-content motif was more frequent in patients with lymphoproliferative disorder than in patients with CHC (49.0% versus 34.4%, p = 0.02, Fig. 3A). In contrast, in lymphoproliferative disorder patients, Cent 6 motif was lower than in CHC patients (9.1% versus 18.4%, p = 0.03, Fig. 3A). The Cent2 and Cent6 motifs were present in about 44.3% and 7.6% of the HCV-negative patients. The KIR2DL3 gene distinguishes Cent 2 from Cent 6 motifs and was more frequent in lymphoproliferative disorder compared to CHC: however, this difference did not reach statistical significance (87.6% versus 80.8%, Fig. 1C).
Distribution of Cent2 and Cent6 locus (panel A, Table 2). Distribution of Tel3 and Tel6 loci (panel B, Table 2). Distribution of Cent/Tel1 loci (panel C, Table 2). Distribution of Cent/Tel1 with A/B KIR2DL5 variant (panel D, Table 3). Distribution of Cent-2DS3/5 (1) loci (panel E, Table 3). Distribution of Tel-2DS4 (3) AND Tel-2DS4 (8) locus (panel F, Table 3). * Fisherx’s Exact test, p<0.05.
Regarding the telomeric KIR motifs, the Tel 3 motif was completely absent in CHC patients, but present in HCV-negative patients (0% versus 3.6%, Fisher’s exact test, p = 0.03 Table 2). Tel 3 motif was found to be characterized by the absence of the KIR3DS1 gene.
Tel 6 motif characterized by the absence of both the KIR3DL1 and 2DS4 genes and the presence of the KIR3DS1 and 2DS1 genes was more frequent in CHC patients than in HCC (8.8% versus 1.7%, p = 0.02, Fig. 3B). The Tel 6 motif was found in 6.3% of HCV-negative patients.
Centromeric/telomeric KIR locus
The Cent/Tel1 motif was characterized by the presence of the KIR2DL5, KIR2DS3 and KIR2DS5 genes, and its frequency was lower in lymphoproliferative disorder cases than in CHC cases (6.5% versus 14.4%, p = 0.04, Fig. 3C) In HCV-negative patients, it was found to be 13.8%. Among these genes, the KIR2DS3 compared with CHC patients was less expressed in HCC patients (30.5% versus 43.2%, p = 0.05, Fig. 1C) and in lymphoproliferative disorder patients (32.0% versus 43.2%) but the difference was not statistical significant.
Genomic blocks sharing specific sets of alleles in strong linkage disequilibrium
In European-American individuals, linkage disequilibrium was found to be particularly strong between the KIR2DL5-KIR2DS3/KIR2DS5, KIR2DS3/KIR2DS5-KIR2DL1 and KIR3DL1-KIR2DS4 gene pairs [28]. Moreover, the KIR2DL5 gene A and B variants are known to discriminate the centromeric (B) from the telomeric (A) component of the KIR gene locus [29], and two functionally related variants of the KIR2DS4 gene (full and del) are reported. The study combined all of this information in order to estimate the weight of single gene/variant discrimination. The presence of both 2DL5 A/B variants and the 2DS3 with 2DS5 genes (Cent/Tel1 A/B, Table 3) was less frequent in lymphoproliferative disorders compared to patients with chronic HCV infection (2.6% versus 11.2%, p < 0.01, Fig. 3D). Rather, among HCV-negative patients the frequency of this motif was 10.1%. Furthermore, the presence of both the KIR2DS3 and the KIR2DS5 genes with KIR2DL1 gene (Cent-2DS3/5 (1), Table 3) was lower in lymphoproliferative disorders cases compared to CHC cases (6.5 versus 14.4%, p = 0.04, Fig. 3E), and its frequency was about 13.8% in HCV-negative patients.
The KIR3DL1-KIR2DS4 block is part of a conserved genotype termed “ancestral genotype.” The presence of both the KIR3DL1 gene and the deleted form of the KIR2DS4 gene (KIR2DS4 del variant, Tel-2DS4 (3), Table 3), which generates a premature stop codon, revealed a frequency higher in HCC cases (70.3%) compared to CHC cases (52.0%, p < 0.01, Fig. 3F). In contrast, the absence of both the KIR3DL1 and KIR2DS4 genes (Tel-2DS4 (8), Table 3) was less frequent in HCC patients than in CHC patients (1.7% versus 8.8%, p = 0.02, Fig. 3F). The frequencies of these two motifs in HCV-negative cases were 63.3% and 6.2% respectively (Table 3).
HLA frequencies
In a previous study, we found a different association between HLA with HCV-related MC and HLA with NHL [33]; to this we analyzed HLA allele frequencies by dividing the group of patients with lymphoproliferative disorders in the benign MC disease and the more malignant NHL disease.
The frequencies of HLA class I ligands (HLA-A, HLA-B and HLA-C) are reported in Table 4 for both groups of HCV-positive patients and HCV-negative individuals.
Compared to all the other groups, the MC group showed a significantly different frequency of HLA genes. The proportion of patients with HLA-A*03 or HLA-A*11 alleles was found to be lower in lymphoproliferative disorders cases than in CHC (19.6% versus 32.8%, p < 0.01). The proportion of HCV-negative individuals with theses alleles was also found to be lower (25.4%).
In addition, differences in frequencies were found for the HLA-Bw4 ligand: the proportion of patients with at least one HLA-Bw4 alleles (Bw4–80I or Bw4–80T) was higher in MC cases than in NHL cases (76.0% versus 56.4%, p = 0.02). This difference was related to the presence of Bw4–80I variant. In fact, the frequency of patients with at least one Bw4–80I allele was higher in MC patients than in NHL (60.0% versus 32.0%, p < 0.01).
KIR / HLA pair frequencies
Known KIR gene-HLA ligand combinations were examined. Statistical difference between KIR-HLA pair frequencies are reported in Table 5. The frequency of KIR3DL1 gene was found to be higher in HCC cases than in CHC (98.3% versus 90.4%, p = 0.01). Furthermore, the proportion of patients having KIR3DL1 with at least one HLA-Bw4 allele (Bw4–80I or Bw4–80T) was lower in patients with NHL than in HCC and MC patients (55.4% versus, 69.8% and 75.7%, p = 0.05 and p = 0.01, respectively). Moreover, a significant negative trend was detected for the 3DL1/HLA-B genotype frequencies in NHL patients (from 42.3% for 3DL1/Bw6-Bw6 to 2.6% for both 3DL1/Bw4-I80-Bw4-I80 and 3DL1/Bw4-T80-Bw4-T80 genotypes, Cochran-Armitage trend test z2 = 1.98, p = 0.047) (Fig. 4).
It is showed a negative trend association between the KIR3DL1/HLA-Bw4 with the presence of HCV-related lymphoma.
The frequency of KIR2DS4 was found to be higher in HCC patients compared to CHC patients (98.3% versus 91.2%, p = 0.02) but independent from the presence of the HLA-Cw04 ligand (Table 5).
The frequency of the KIR3DL2 gene combining with either HLA-A03+ or HLA-A11+ ligands was higher in CHC patients than in MC patients (32.8% versus, 14.7% p < 0.01) (Table 5).
Activator/Inhibitor KIR genes and HLA matched ligands. Two different patterns of KIR gene and related combinations of HLA ligands werec examinated (Table 6): the inhibitor motif with all of the 2DL3/2DL2/2DL1/3DL1 KIR genes and the activator motif with the 2DS2/2DS1/3DS1KIR genes.
The frequency of patients having a full set of KIR-inhibitor genes was found to be higher among lymphoproliferative disorders cases than CHC cases (47.7% versus 32.8%, Fisher’s exact test, p = 0.01). However, in patients with the full set of KIR-inhibitors the additional presence of HLA ligands (HLA-C1, HLA-C2 and HLA-Bw4) did not show significant differences in frequency. This suggests that the function of the full set of KIR-inhibitors was not related to the HLA ligands.
The frequency of patients with a full activator motif was found to be lower, but not statistically significant, in patients with lymphoproliferative disorders (19.0%) than in others patients (about 24–26%). However, the proportion of patients showing both activator motif and at least one HLA-Bw4 ligand was found to be lower in HCC patients than in CHC patients (50% versus 78.1%, p = 0.03).
Discussion
The main result of our study is to have shown a significant association between KIR/HLA genotypes to direct toward HCV-related lymphoproliferative diseases. By means of a well represented group of HCV-infected patients, this study compared patients developing a chronic infection with those developing either a malignant hepatic disorder (HCC), or a lymphoproliferative disease. Focus was placed on the identification of KIR gene receptors involved in these different HCV-related diseases as their identification may lead to better understanding of HCV-pathogenesis.
Previous studies identified KIR2DS3 [31] and HLA-C1+ KIR2DL3+, especially when in homozygous [12;13;30] as NK cell-associated KIR genes at a higher frequency in patients who resolved infection compared with CHC patients. Even though a statistically significant difference was not found, we found a trend of reduction of KIR2DL3 gene frequency and an increase in KIR2DS3 gene frequency in CHC patients compared to HCV-negative patients (Fig. 1C). However, this trend was not found in HCV-related malignancies (Fig. 1C). Thus, these data highlight that, in patients who are unable to eliminate the virus, KIR genes that are usually beneficial for HCV eradication are conversely associated to malignancies. Our results are in agreement with several reports describing the functional impairment of NK cells in chronically HCV-infected patients and an influence of HCV viral load towards an increased risk of HCC [32–35] and lymphoproliferative disorders [40].
The decreased association of KIR2DL3 with CHC was restricted to the single KIR2DL3 gene since the centromeric region Cent 2 and Cent 6 frequencies (different only for the presence of the KIR2DL3 gene) among HCV-related groups was in contrast (Table 2). Furthermore, results suggest that KIR2DL3 was more related with lymphoproliferative disorders when the KIR2DL2 gene was also present (Cent 2 compared to Cent 1) (Table 2 and Fig. 3A). In a study on ligand-instructed models of NK-cell education, the recognition of HLA by an inhibitory KIR2DL2 receptor was demonstrated to suppress the subsequent expression of a second KIR2DL1 receptor [37]. Based on this, we proposed that NK cell inhibition derived from the presence of both the KIR2DL2 and KIR2DL3 genes may be weaker that those derived from NK having the presence of the KIR2DL1 gene, and suggested that this reduction in NK cell inhibition was associated with lymphoproliferative disorder.
As regard the KIR2DS3 gene found associated with CHC (Fig. 1C), an analysis of the telomeric/centromeric regions showed that the region including both the 2DS3 and the KIR2DS5 genes (Cent/Tel1) was mostly reduced in lymphoproliferative disorders patients, while the motif excluding the KIR2DS5 gene (Cent/Tel3) was more reduced in HCC patients (Table 2, Fig. 3). These data suggest that the presence of the KIR2DS5 gene might have a protective effect against lymphoproliferative progression in KIR2DS3+ individuals. Although it is tempting to speculate a direct role for KIR2DL2, KIR2DL3, KIR2DS3 and KIR2DS5 in HCV-related lymphoproliferations, it is likely that these genes are all surrogate markers of the same KIR genotype found at a higher frequency in HCV-related lymphoproliferations compared to CHC cases (genotype ID number 4, Fig. 1A and B). Furthermore, a pool of specific KIR/HLA interactions is thought to play a key role in determining whether the lymphoproliferative disorders result in MC cases rather than in the more malignant NHL. Indeed, HLA-B gene analysis showed a negative trend between HLA-Bw4 variants in terms of ligand provision for KIR3DL1 receptors in malignant NHL (Table 4, 5). Data also showed that interaction of HLA-Bw4–80I with the inhibitory KIR3DL1 receptor may influence MC development; thus, individuals with an HLA-Bw4–80I+ KIR3DL1+ genotype would have a low risk of developing NHL (Fig. 4). This data is intriguing; it has been demonstrated that a reduction in NK activation results from a decrease in the number of strongest inhibitor KIR/HLA combinations (i.e., HLA-C2 > HLA-C1> HLA-Bw4), which is in turn coupled with a decrease in the number of activator KIR/HLA pairs [38]. Overall, the number of inhibitor KIR genes was found to be higher in lymphoproliferative disorder patients than in other groups (Table 6), and simultaneously lymphoproliferative disorders showed a reduction in NK activation due to a lower frequency of the activator KIR genes KIR2DS3 and KIR2DS5 (Fig. 3). It is thus hypothesized that an overall inhibition of NK cells may facilitate lymphoproliferative development, and that the involvement of KIR3DL1 is mainly associated with NHL. The higher frequency of KIR3DL1/Bw6 combination found in NHL (Fig. 4) suggests that a HLA-Bw6/antigen-derived complex originating during HCV infection might constitute a signal for the inhibitory KIR3DL1 receptor, as previously demonstrated in other viral situations [39–41]. Although further studies are necessary to confirm these hypotheses, it is now known that IFN-α therapy, a standard therapy for HCV-infection and HCV-related lymphoproliferative diseases, transiently enhances cytotoxicity of NK cells and triggers their activation, rendering these cells more effective against both infection and lymphoproliferation [42].
The role of the KIR genotype in HCV-related HCC is not quite so clear. NK cells are highly enriched in the liver and have been demonstrated to be able to eliminate HCV-infected and transformed cells, especially after Interferon-α stimulation [43]. A protective role of NK cells to resolve HCV infection and liver disease progression has also been demonstrated [44]. Moreover, among patients with liver transplantation and the HLA-C1+ KIR2DL3+ genotype, those having a higher number of activating KIR genes showed a minor HCV recurrence and progression of HCV-related disease [45;46]. Thus, overall data indicates that NK cells and the linked KIR genotype might influence not only HCV viral load, but also the risk of progression to HCC. The most reliable model to explain this finding is a weaker inhibition of NK cells, thus a greater NK activator function may lead to a more efficient elimination of the HCV-infected cells, ultimately protecting the liver from disease progression [47]. Our study results support this association: we found: (i) a higher frequency of KIR3DL1 and KIR2DS4 genes, both representative of the AA-genotype characterized by a reduction of activator KIR genes [48] in HCC cases than in CHC cases (Fig. 1C) and (ii) an association of HCC with the presence of KIR3DL1 and deleted nonfunctional variant of the KIR2DS4D gene, the unique activator gene presents in the AA-genotype (Tel-2DS4 (3), Table 3 and Fig. 3F). Furthermore, in our series, a main characteristic of HCC was the frequent absence of HLA-Bw4+ KIR3DS1+ interaction (Table 6). A protective effect of this interaction on both cirrhosis and HCC has been previously reported in several studies [12;14;32]. Our data confirm this association and highlight its role in counterbalancing the inhibitory effect of HLA-Bw4+ KIR3DL1 interaction.
This study, to the authors’ knowledge, is the first report highlighting the potential impact of the immunogenetic background to address the development of HCV-related lymphoproliferative disorders. HCV-related disorders are multifactorial and complex diseases, but our study highlighted that the persistence of HCV infection after antiviral-therapy remains a significant factor associated with the development of HCC and lymphorpoliferative diseases. In this context, the suggestion of a role where the KIR/HLA genetic background is affecting tumor development and clinical response in the set of HCV-infected patients is intriguing and requires further investigations disclosing perspectives for the development of combining NK-based therapies with conventional antitumoral and anti-HCV treatments.
Supporting Information
S1 Table. Rare KIR genotypes.
27 rare genotype were identified and the ID number reported are those from the reference database [28]. The presence of KIR genes is indicated by the presence of X symbol. Genotypes AA and BX according to criteria reported in material and method section are indicated in the first column.
https://doi.org/10.1371/journal.pone.0117420.s001
(DOCX)
Author Contributions
Conceived and designed the experiments: VD LC. Performed the experiments: LC MD. Analyzed the data: VD LC MD AG. Wrote the paper: VD LC OR AM MT FB AZ PD. Agree with manuscript results and conclusions: VD LC MD OR AZ FI MT FB AM DS VR SD PP EV MB CM ML AG AZ RC PD. Enrolled patients: AZ FI MT FB AM DS VR SD PP EV MB CM ML RC.
References
- 1. Afdhal N, Zeuzem S, Kwo P, Chojkier M, Gitlin N, et al. (2014) Ledipasvir and Sofosbuvir for Untreated HCV Genotype 1 Infection. N Engl J Med. pmid:25548814
- 2. De Re V, De Vita S, Sansonno D, Toffoli G (2008) Mixed cryoglobulinemia syndrome as an additional autoimmune disorder associated with risk for lymphoma development. Blood 111: 5760. pmid:18544704
- 3. Park SH, Shin EC, Capone S, Caggiari L, De Re V, et al. (2012) Successful vaccination induces multifunctional memory t-cell precursors associated with early control of hepatitis C virus. Gastroenterology 143: 1048–1060. pmid:22705008
- 4. Thimme R, Binder M, Bartenschlager R (2012) Failure of innate and adaptive immune responses in controlling hepatitis C virus infection. FEMS Microbiology Reviews 36: 663–683. pmid:22142141
- 5. Cheent K, Khakoo SI (2011) Natural killer cells and hepatitis C: Action and reaction. Gut 60: 268–278. pmid:20876777
- 6. Orange JS (2002) Human natural killer cell deficiencies and susceptibility to infection. Microbes and Infection 4: 1545–1558. pmid:12505527
- 7. Moesta AK, Parham P (2012) Diverse functionality among human NK cell receptors for the C1 epitope of HLA-C: KIR2DS2, KIR2DL2, and KIR2DL3. Front Immunol 3:336. eCollection;%2012.: 336. . pmid:23189078
- 8. Larkin J, Bost A, Glass JI, Tan SL (2006) Cytokine-activated natural killer cells exert direct killing of hepatoma cells harboring hepatitis C virus replicons. J Interferon Cytokine Res 26: 854–865. pmid:17238828
- 9. Stegmann KA, Bjorkstrom NK, Veber H, Ciesek S, Riese P, et al. (2010) Interferon-alpha-induced TRAIL on natural killer cells is associated with control of hepatitis C virus infection. Gastroenterology 138: 1885–1897. pmid:20334827
- 10. Crotta S, Brazzoli M, Piccioli D, Valiante NM, Wack A (2010) Hepatitis C virions subvert natural killer cell activation to generate a cytokine environment permissive for infection. J Hepatol 52: 183–190. pmid:20015567
- 11. Li Y, Zhang T, Ho C, Orange JS, Douglas SD, et al. (2004) Natural killer cells inhibit hepatitis C virus expression. J Leukoc Biol 76: 1171–1179. pmid:15339939
- 12. Khakoo SI, Thio CL, Martin MP, Brooks CR, Gao X, et al. (2004) HLA and NK cell inhibitory receptor genes in resolving hepatitis C virus infection. Science 305: 872–874. pmid:15297676
- 13. Knapp S, Usama W, Hegazy D, Brackenbury L, Guha IN, et al. (2010) Consistent beneficial effects of killer cell immunoglobulin-like receptor 2dl3 and group 1 human leukocyte antigen-c following exposure to hepatitis c virus. Hepatology 51: 1168–1175. pmid:20077564
- 14. Lopez-Vazquez A, Rodrigo L, Martinez-Borra J, Perez R, Rodriguez M, et al. (2005) Protective effect of the HLA-Bw4I80 epitope and the killer cell immunoglobulin-like receptor 3DS1 gene against the development of hepatocellular carcinoma in patients with hepatitis C virus infection. J Infect Dis 192: 162–165. pmid:15942906
- 15. Littera R, Zamboni F, Tondolo V, Fantola G, Chessa L, et al. (2013) Absence of activating killer immunoglobulin-like receptor genes combined with hepatitis C viral genotype is predictive of hepatocellular carcinoma. Hum Immunol 74: 1288–1294. pmid:23756163
- 16. Miyagi T, Takehara T, Tatsumi T, Kanto T, Suzuki T, et al. (2003) CD1d-mediated stimulation of natural killer T cells selectively activates hepatic natural killer cells to eliminate experimentally disseminated hepatoma cells in murine liver. Int J Cancer 106: 81–89. pmid:12794761
- 17. Tatsumi T, Takehara T, Yamaguchi S, Sasakawa A, Sakamori R, et al. (2007) Intrahepatic delivery of alpha-galactosylceramide-pulsed dendritic cells suppresses liver tumor. Hepatology 45: 22–30. pmid:17187419
- 18. Miyagi T, Takehara T, Tatsumi T, Suzuki T, Jinushi M, et al. (2004) Concanavalin a injection activates intrahepatic innate immune cells to provoke an antitumor effect in murine liver. Hepatology 40: 1190–1196. pmid:15486936
- 19. Bray SM, Vujanovic L, Butterfield LH (2011) Dendritic cell-based vaccines positively impact natural killer and regulatory T cells in hepatocellular carcinoma patients. Clinical and Developmental Immunology 2011. pmid:22347323
- 20. Cai L, Zhang Z, Zhou L, Wang H, Fu J, et al. (2008) Functional impairment in circulating and intrahepatic NK cells and relative mechanism in hepatocellular carcinoma patients. Clinical Immunology 129: 428–437. pmid:18824414
- 21. Gao B, Radaeva S, Park O (2009) Liver natural killer and natural killer T cells: Immunobiology and emerging roles in liver diseases. Journal of Leukocyte Biology 86: 513–528. pmid:19542050
- 22. Zhao F, Korangy F, Greten TF (2012) Cellular immune suppressor mechanisms in patients with hepatocellular carcinoma. Digestive Diseases 30: 477–482. pmid:23108303
- 23. Bruix J, Sherman M, Lloret JM, Beaugrand M, Lencioni R, et al. (2001) Clinical management of hepatocellular carcinoma. Conclusions of the Barcelona-2000 EASL conference. European Association for the Study of the Liver. J Hepatol.35: 421–430. pmid:11592607
- 24. Quartuccio L, Isola M, Corazza L, Ramos-Casals M, Retamozo S, et al. 2014 Validation of the classification criteria for cryoglobulinaemic vasculitis. Rheumatology (Oxford). 53: 2209–2213. pmid:24994905
- 25.
X-Jaffe ES, Harris NL, Stein H, Vardiman JW (2001). Pathology and Genetics of Tumours of Haematopoietic and Lymphoid Tissues. World Health Organization Classification of Tumours. Lyon, France: IARC Press; 2001.
- 26. De Re V, Gragnani L, Fognani E, Piluso A, Izzo F, et al. (2014) Impact of Immunogenetic IL28B Polymorphism on Natural Outcome of HCV Infection. Biomed Res Int 2014:710642. Epub;%2014 Feb 26.: 710642. . pmid:24707497
- 27. Caggiari L, Toffoli G, De Re V, Orzes N, Spina M, et al. (2011) KIR/HLA combination associated with the risk of complications in celiac disease. Int J Biol Markers 26: 221–228. pmid:22180175
- 28. Gonzalez-Galarza FF, Christmas S, Middleton D, Jones AR (2011) Allele frequency net: a database and online repository for immune gene frequencies in worldwide populations. Nucleic Acids Res 39: D913–D919 pmid:21062830
- 29. .Gourraud PA, Meenagh A, Cambon-Thomsen A, Middleton D (2010) Linkage disequilibrium organization of the humanKIR superlocus: implications for KIR data analyses. Immunogenetics 62: 729–740 pmid:20878401
- 30. De Re V, Caggiari L, de Zorzi M, Talamini R, Racanelli V, et al. (2014) Genetic diversity of the KIR/HLA system and outcome of patients with metastatic colorectal cancer treated with chemotherapy. PLoS One 9: e84940. pmid:24497922
- 31. Pyo CW, Guethlein LA, Vu Q, Wang R, Abi-Rached L, et al. (2010) Different patterns of evolution in the centromeric and telomeric regions of group A and B genotypes of the human killer cell Ig-like receptor locus. PLoS ONE 5. pmid:21283510
- 32. Gomez-Lozano N, Gardiner CM, Parham P, Vilches C (2002) Some human KIR genotypes contain two KIR2DL5 genes: KIR2DL5A and KIR2DL5B. Immunogenetics 54: 314–319. pmid:12185535
- 33. De Re V, Caggiari L, Monti G, Libra M, Spina M, Dolcetti R, et al (2010). HLA DR-DQ combination associated with the increased risk of developing human HCV positive non-Hodgkin’s lymphoma is related to the type II mixed cryoglobulinemia. Tissue Antigens.75:127–35 pmid:20002609
- 34. Zuniga J, Romero V, Azocar J, Terreros D, Vargas-Rojas MI, et al. (2009) Protective KIR-HLA interactions for HCV infection in intravenous drug users. Mol Immunol 46: 2723–2727. pmid:19552960
- 35. Dring MM, Morrison MH, McSharry BP, Guinan KJ, Hagan R, et al. (2011) Innate immune genes synergize to predict increased risk of chronic disease in hepatitis C virus infection. Proc Natl Acad Sci U S A 108: 5736–5741. pmid:21402922
- 36. Cariani E, Pilli M, Zerbini A, Rota C, Olivani A, et al. (2013) HLA and killer immunoglobulin-like receptor genes as outcome predictors of hepatitis C virus-related hepatocellular carcinoma. Clin Cancer Res 19: 5465–5473. pmid:23938290
- 37. Chang KC, Wu YY, Hung CH, Lu SN, Lee CM, et al. (2013) Clinical-guide risk prediction of hepatocellular carcinoma development in chronic hepatitis C patients after interferon-based therapy. Br J Cancer 109: 2481–2488. pmid:24084770
- 38. Lee MH, Yang HI, Lu SN, Jen CL, Yeh SH, et al. (2010) Hepatitis C virus seromarkers and subsequent risk of hepatocellular carcinoma: long-term predictors from a community-based cohort study. J Clin Oncol %20;28: 4587–4593. pmid:20855826
- 39. Shindoh J, Hasegawa K, Matsuyama Y, Inoue Y, Ishizawa T, et al. (2013) Low hepatitis C viral load predicts better long-term outcomes in patients undergoing resection of hepatocellular carcinoma irrespective of serologic eradication of hepatitis C virus. J Clin Oncol %20;31: 766–773. pmid:23129744
- 40. Kawamura Y, Ikeda K, Arase Y, Yatsuji H, Sezaki H, et al. (2007) Viral elimination reduces incidence of malignant lymphoma in patients with hepatitis C. Am J Med 120: 1034–1041. pmid:18060923
- 41. Schonberg K, Sribar M, Enczmann J, Fischer JC, Uhrberg M (2011) Analyses of HLA-C-specific KIR repertoires in donors with group A and B genotypes suggest a ligand-instructed model of NK cell receptor acquisition. Blood 117: 98–107. pmid:20935255
- 42. Ahlenstiel G, Martin MP, Gao X, Carrington M, Rehermann B (2008) Distinct KIR/HLA compound genotypes affect the kinetics of human antiviral natural killer cell responses. J Clin Invest 118: 1017–1026. pmid:18246204
- 43. Sharma D, Bastard K, Guethlein LA, Norman PJ, Yawata N, et al. (2009) Dimorphic motifs in D0 and D1+D2 domains of killer cell Ig-like receptor 3DL1 combine to form receptors with high, moderate, and no avidity for the complex of a peptide derived from HIV and HLA-A*2402. J Immunol 183: 4569–4582. pmid:19752231
- 44. Parsons MS, Wren L, Isitman G, Navis M, Stratov I, et al. (2012) HIV infection abrogates the functional advantage of natural killer cells educated through KIR3DL1/HLA-Bw4 interactions to mediate anti-HIV antibody-dependent cellular cytotoxicity. J Virol 86: 4488–4495. pmid:22345455
- 45. Fadda L, ’Connor GM, Kumar S, Piechocka-Trocha A, Gardiner CM, et al. (2011) Common HIV-1 peptide variants mediate differential binding of KIR3DL1 to HLA-Bw4 molecules. J Virol 85: 5970–5974. pmid:21471246
- 46. Hu X, Jiang Y, Li X, Gao Y, Guo X, et al. (2013) Long-Term Effect on Natural Killer Cells by Interferon-alpha Therapy on the Outcomes of HCV Infection. J Interferon Cytokine Res. pmid:24367931
- 47. Stegmann KA, Bjorkstrom NK, Ciesek S, Lunemann S, Jaroszewicz J, et al. (2012) Interferon alpha-stimulated natural killer cells from patients with acute hepatitis C virus (HCV) infection recognize HCV-infected and uninfected hepatoma cells via DNAX accessory molecule-1. J Infect Dis 205: 1351–1362. pmid:22457290
- 48. Yamagiwa S, Matsuda Y, Ichida T, Honda Y, Takamura M, et al. (2008) Sustained response to interferon-alpha plus ribavirin therapy for chronic hepatitis C is closely associated with increased dynamism of intrahepatic natural killer and natural killer T cells. Hepatol Res 38: 664–672. pmid:18328072
- 49. de Arias AE, Haworth SE, Belli LS, Burra P, Pinzello G, et al. (2009) Killer cell immunoglobulin-like receptor genotype and killer cell immunoglobulin-like receptor-human leukocyte antigen C ligand compatibility affect the severity of hepatitis C virus recurrence after liver transplantation. Liver Transpl 15: 390–399. pmid:19326408
- 50. Legaz I, Lopez-Alvarez MR, Campillo JA, Moya-Quiles MR, Bolarin JM, et al. (2013) KIR gene mismatching and KIR/C ligands in liver transplantation: consequences for short-term liver allograft injury. Transplantation 95: 1037–1044. pmid:23478359
- 51. Winter CC, Gumperz JE, Parham P, Long EO, Wagtmann N (1998) Direct binding and functional transfer of NK cell inhibitory receptors reveal novel patterns of HLA-C allotype recognition. J Immunol 161: 571–577. pmid:9670929
- 52. Martin AM, Kulski JK, Gaudieri S, Witt CS, Freitas EM, et al. (2004) Comparative genomic analysis, diversity and evolution of two KIR genotypes A and B. Gene 335:121–31. pmid:15194195