Abstract

The genetic variability of the host contributes to the risk of human papillomavirus (HPV)-related cervical disease. Immune response genes to HPV must be investigated to define patients with the highest risk of developing malignant disease. The aim of this study was to investigate the association of polymorphic immune response genes, namely KIR, HLA class I and II, and single-nucleotide polymorphisms (SNPs) of cytokines with HPV-related cervical disease. We selected 79 non-related, admixed Brazilian women from the state of Paraná, southern region of Brazil, who were infected with high carcinogenic risk HPV and present cervical intraepithelial neoplasia grade 3 (CIN3), and 150 HPV-negative women from the same region matched for ethnicity. KIR genes were genotyped using an in-house PCR-SSP. HLA alleles were typed using a reverse sequence-specific oligonucleotide technique. SNPs of TNF −308G>A, IL6 −174G>C, IFNG +874T>A, TGFB1 +869T>C +915G>C, and IL10 −592C>A −819C>T −1082G>A were evaluated using PCR-SSP. The KIR genes were not associated with HPV, although some pairs of i(inhibitory)KIR-ligands occurred more frequently in patients, supporting a role for NK in detrimental chronic inflammatory and carcinogenesis. Some HLA haplotypes were associated with HPV. The associations of INFG and IL10 SNPs potentially reflect impaired or invalid responses in advanced lesions.

1. Introduction

Human papillomavirus (HPV) infections occur frequently in healthy individuals, and high carcinogenic risk (HR) HPV types are a major causal factor for cervical cancer. Persistent infection with one among approximately 15 genotypes of carcinogenic HPV causes almost all cases of cervical cancer; type 16 and HPV-18 account for more than 70% of the cervical cancers detected worldwide. Despite being considered a preventable disease, cervical cancer remains the second most common malignancy among women worldwide, with a higher incidence in underdeveloped countries [1, 2].

The major mechanisms by which HPV contribute to neoplastic initiation and progression involve the activity of two viral oncoproteins, E6 and E7, which interfere with the critical cell cycle tumor suppressive proteins p53 and retinoblastoma (Rb). However, HPV infection alone is not sufficient to induce malignant transformation. The multistep process of tumor formation requires the contribution of other significant cofactors, such as individual genetic variations, intratypic HPV variability, and environmental factors [1, 2]. The genetic variability of the host also plays a role in the risk of developing cervical cancer, especially variability of genes that control the immune response. These highly polymorphic genes are important risk determinants of HPV persistence and disease progression.

The innate immune system comprises the first line of defense following HPV infection. It provides nonspecific protection and enhances the adaptive immune response [3]. Inflammatory cell infiltration occurs in response to HPV tissue damage, with infiltrates consisting initially of neutrophils followed by macrophages and T lymphocytes cells. NK and NKT cells contribute to antiviral innate immune responses. NK cell activation depends on type 1 interferon and proinflammatory cytokines such as IL-12 and IL-18; these cells are able to detect decreased expression of HLA class I in infected and transformed cells [4]. Most cervical HPV infections are cleared or suppressed via cell-mediated immunity: CD4+ and CD8+ T cells are the major effector cells [4], and the Th1 response is associated with clearance of the HPV infection and regression of the cervical cancer [5]. Th2 responses are associated with cervical carcinogenesis [6].

To define patients with the highest risk of developing malignant diseases, the interaction between the host immune response and HPV infection must be investigated. The goal of the present study was to investigate the association of the polymorphic immune response genes, namely, the KIR genes, HLA classes I and II, and SNPs of cytokines, with HPV infection in Brazilian women.

The KIR locus comprises an approximately 150 kb region located on chromosome 19q13.4, which encodes a group of inhibitory and activating KIR molecules. KIRs are key receptors of human natural killer (NK) cells, a subset of lymphocytes that trigger early innate immune response against infection and tumors [7].

The major histocompatibility complex (MHC), also known as the human leukocyte antigen (HLA) complex, located on chromosome 6p21.3, is the most polymorphic genetic system in mammalians and has been studied with regard to a wide variety of diseases of distinct etiology. The fundamental role of the different molecules within the MHC is antigen processing and presentation to the T-cell receptor (TCR), which is crucial for the cell interactions in cell-mediated immunity [8].

Polymorphisms of regulating regions of cytokine genes have been correlated with its production and can confer flexibility in the immune response to the viral infections and cancer biology. Five independent regions were investigated: chromosome 1: IL10 region [9], chromosome 6: tumor necrosis factor (TNF) [10], chromosome 7: interleukin-6 (IL6) [11], chromosome 12:   interferon-gamma (IFNG) [12], and chromosome 19: transforming growth factor-beta (TGFB1) region [13].

2. Materials and Methods

2.1. Patients and Samples

Patients comprised 79 nonrelated, admixed Brazilian women from the state of Paraná in the southern region of Brazil, who were infected with HR-HPV and present cervical intraepithelial neoplasia (CIN) grade 3 (CIN3) and 150 women from the same region matched for ethnicity who were HPV-negative/normal cytology. The study protocol was approved by the ethics committee, and all selected patients signed the free and informed consent form.

In the Paraná state, the degree of the European ancestry is high (80.6%), with a small but significant contribution of African (12.5%) and Amerindian (7.0%) genes according to Probst et al. [14], and the studied populations were considered admixed. The risk of population stratification bias, due to differences in ethnic background between patients and controls, and variations of allele frequencies, according to ethnic background, were minimized by matching patients with control individuals of the same ethnic background, mean age, gender rates, and residence in the same geographical areas.

The patients were diagnosed with high-grade squamous intraepithelial lesion (HSIL) by cytologic smears, CIN3 by histopathology, and also with HR-HPV.

2.1.1. Cytology and Histopathology

The cervical and endocervical material was collected with the aid of an Ayre spatula and a cytobrush for cervical smears and for PCR amplification (suspended in 1 mL of 0.9% NaCl solution and stored at −20°C until analysis). The cytological smears were evaluated and reported according to the Bethesda system as atypical squamous cells of undetermined significance (ASC-US); atypical squamous cells of undetermined significance, which cannot exclude a high-grade squamous intraepithelial lesion (ASC-H); low-grade squamous intraepithelial lesion (LSIL); high-grade squamous intraepithelial lesion (HSIL); in situ or invasive adenocarcinoma (ISCC); or invasive squamous cell carcinoma (SCC).

The cytological criteria for HSIL diagnosis adopted were squamous cells, either isolated or present in small fragments with fewer than ten cells. The cells were the length of the metaplastic cells, showing an increase in the proportion in the nuclear area. The nuclear irregularities, including hyperchromasia, chromatic clustering, irregularity, thickening, or multinucleation, were also used as important cytological criteria [15]. The histopathology findings of biopsy samples were classified as CIN grades I, II, or III, microinvasive or invasive squamous cell carcinoma, or in situ or invasive adenocarcinoma. The histological criteria for CIN were the failure of maturation of the squamous epithelium, nuclear hyperchromasia, and an increased nucleus/cytoplasm ratio. The intensity of these criteria was used to stage the degree of CIN or carcinomas [16]. The HSIL cytological cases included in the present study were confirmed as CIN III by histopatnology.

2.1.2. HPV Molecular Detection

For HPV molecular detection, genomic DNA was extracted using DNAzol (Invitrogen, Carlsbad, CA, USA). HPV polymerase chain reaction (PCR) amplification for HPV was carried out using MY09 (5′-CGTCCMAARGGAWACTGATC-3′)/MY11(5′-GCMCAGGGWCATAAYAATGG-3′), and the PCR product was electrophoresed on a 1.5% agarose gel, stained with 1 μg/mL ethidium bromide, and photodocumented under UV light (approximately 450 bp). Coamplification of the human beta-globin gene (approximately 268 bp) was performed as an internal control, using primers GH20 (5′-GAAGAGCCAAGGACAGGTAC-3′) and PC04 (5′-CAACTTCATCCACG TTCACC-3′) under the same conditions as the HPV PCR. Two types of controls were also included in each reaction series: “no DNA” (negative control) and “HPV-positive DNA” (positive control) [17].

Genotyping was carried out using PCR-based restriction fragment length polymorphism analysis using HpyCH4V (New England Biolabs, Inc., Ipswich, MA, USA). The following HPV genotypes were determined for this genotyping method: HR (16, 18, 31, 33, 35, 39, 45, 51, 52, 56, 58, 59, 68, 66, 73, and 82), UR—undetermined risk (26 and 53), and LR—low-risk (6, 11, 30, 34, 40, 42, 43, 44, 54, 55, 61, 62, 64, 67, 69, 70, 72, 74, 81, 83, 84, and 91). The genotypes were grouped according to the International Agency for Research on Cancer (IARC) based on the carcinogenic potential and evolutionary branch [18].

2.2. Genotyping of KIR, HLA, and Cytokine Genes

Genomic DNA samples were extracted from 150 μL of the buffy coat obtained from 5 mL of EDTA anticoagulant peripheral blood using the EZ-DNA Kit (Biological Industries, Beit Haemek, Israel). The DNA concentration was then determined using a Qubit fluorometer (Life Technologies Corporation, Eugene, Oregon, USA).

All genotyping methods were validated using previously typed and tested reference samples. Positive and negative controls were included in all genotyping method.

2.2.1. KIR Genes Genotyping

Fourteen KIR genes and one pseudogene (KIR2DL1, KIR2DL2, KIR2DL3, KIR2DL4, KIR2DL5, KIR2DS1, KIR2DS2, KIR2DS3, KIR2DS4, KIR2DS5, KIR3DL1, KIR3DL2, KIR3DL3, KIR3DS1, and KIR2DP1) were studied using an in-house polymerase chain reaction using the sequence-specific primer method (PCR-SSP) according to Martin et al. [19] and adapted by Rudnick et al. [20]. Primers were synthesized by Invitrogen (Life Technologies Corporation, Grand Island, NY, USA), and the amplified products were visualized by 2% agarose gel electrophoresis.

2.2.2. HLA Classes I and II Typing

HLA classes I and II allele typing was conducted using the reverse sequence-specific oligonucleotide technique (rSSO; One Lambda Inc., Canoga Park, CA, USA) with Luminex xMap technology (Luminex Corporation, Austin, USA). HLA groups 1 (C1) and 2 (C2) of HLA-C and group Bw4 of HLA-B were defined according to Carrington and Norman [21] and Petersdorf [22].

2.2.3. Genotyping of SNPs in Cytokine Genes

Sequence-specific primer PCR (PCR-SSP; One Lambda Cytokine Genotyping Primer Pack, One Lambda, CA, USA) was performed to genotype the following SNPs: TNF −308G>A (rs1800629), IFNG +874T>A (rs2430561), IL6 −174G>C (rs1800795), IL10 −1082G>A (rs1800896), IL10 −819C>T (rs1800871), IL10 −592C>A (rs1800872), TGFB −509T>C (rs1800469), and TGFB1 +915G>C (rs1800471) (Table 1), according to the manufacturer’s instructions.

2.3. Statistical Analyses

Allele, genotype, and haplotype frequencies of KIR, HLA, and cytokines were calculated by direct counting. Fisher’s exact test and the chi-square test with Yates’ correction were used for statistical comparisons. were considered significant, and values were adjusted by means of the Bonferroni correction to enable multiple comparisons. The odds ratio (OR) was calculated based on the cross product ratio and the exact 95% confidence intervals (CI) using the SISA statistical package (http://www.quantitativeskills.com/sisa/index.htm). Hardy-Weinberg equilibrium [23] was determined by calculating the expected genotype frequencies and comparing them to the observed values using Arlequin software version 3.1 (http://cmpg.unibe.ch/software/arlequin3/).

3. Results

The distributions of allele frequency ratios for all analyzed genes and for KIR haplotype frequencies were in Hardy-Weinberg equilibrium.

There were no significant differences between KIR genes frequencies in patients and controls (Table 2), and the frequency distribution was similar to that reported in another study of the same region [24].

There was no relationship in the frequencies of ligands (C1, C2, Bw4, and HLA-A3/11) and in the combination of KIR-HLA ligands with the HPV disease (Table 3).

The number and type of inhibitory KIR-HLA pairs were evaluated (Table 4), and there was a greater frequency in the patients of three pairs (38.0%) and two pairs (36.6%), followed by one pair (14.1%) and four pairs (11.2%). In the controls, two pairs (46.6%) were the most frequent, followed by three (32.2%), one (11.0%), and four (10.2%) pairs. The pairs KIR3DL1-Bw4 and KIR3DL2-HLA-A3/11 were not detected in either group. Significant difference was observed between patients and controls with respect to the three pairs KIR2DL2/3-C1, KIR3DL1-Bw4, and KIR3DL2-A3/11, which were more frequent in the patients.

Twenty-six KIR haplotypes were observed in HPV patients, 10 of them were not found in controls, and, otherwise, 13 others haplotypes were present only in controls (Figure 1). There were no differences between patients and controls.

No differences were observed in the distribution of HLA-A, B, and DRB1 allele groups between patients and controls (Table 5).

Only the HLA-A*02-HLA-B*51 haplotype showed a reduced frequency in HPV patients (0.006 versus 0.052, , OR = 0.1186, and 95% CI = 0.015–0.8717) than in controls; the other six haplotypes were more frequent in the patients, but a large CI was obtained due to the small number of patients and controls (Table 6).

The cytokine allele frequencies did not differ between HPV patients and controls (Table 7), and the frequencies distribution were consistent with the results of a previous study of the same region [25].

Based on the genotypes, the phenotypes of cytokines production level were inferred (low, intermediate, or high). There were significant differences for the low producer phenotypes of INF-γ, defined by genotype AA [12], which had an increased frequency in patients (35.90 versus 29.59; ; OR = 1.81; 95% CI = 1.18–4.60), and for an intermediate producer of IL-10 (32.05% versus 48.00%; ; OR = 0.1607; 95% CI = 0.28–0.95) defined by the haplotype GCC/ACC and GCC/ATA (Table 8). The GCC/GCC genotype (high producer phenotype of IL-10) [9] was more frequent among patients (20.5%) compared with controls (12%), although this difference was not significant.

4. Discussion

It is widely accepted that cofactors, including endogenous hormones and genetic factors, such as HLA and other genes related to the host immune response, may have important roles in the development of HPV-cervical lesions [26]. Inherited genetic polymorphisms within immune response genes have been shown to be associated with an increased risk of invasive cervical cancer and its immediate precursor, cervical intraepithelial neoplasia grade 3 [27]. An inappropriate innate and specific immune response may increase the risk of lesions and disease progression.

4.1. KIR and Their HLA Ligands

NK cells play an important role in innate immunity against infected and transformed cells as part of the immune surveillance process. KIR genes encode molecules that convey either inhibitory or activating signals (iKIR and aKIR) to NK cells and to a subset of CD8+ T cells. Binding of iKIR (designated 2DL and 3DL) to specific HLA allotypes has been clearly demonstrated and correlated to the ability to inhibit NK cytolysis of target cells bearing these HLA molecules. These interactions are remarkably complex, and synergistic relationship between these polymorphic loci may regulate NK cell-mediated immunity against viral infections [28].

No relationship was found between KIR genes and HPV-related cervical disease in Brazilian patients, consistent with the findings of Song et al. in Korean patients [29]. Although not significant, KIR2DS3, KIR2DS5, and KIR2DL5 were more frequent in the patient group. KIR2DL5, an inhibitory receptor, possesses a combination of genetic, structural, and functional features that make it unique among the KIR [30] and potentially contribute to HPV pathogenesis. KIR3DS1 can induce a persistent, weak inflammatory reaction to HPV that results in continuous tissue injury, similar to HBV-susceptible genes [31]. Carrington et al. [32] found that the presence of the activating KIR3DS1 is related to an increased risk of neoplasia, particularly in the absence of protective inhibitory KIR-HLA. In contrast, Arnheim et al. [33] indicated that the inhibitory allele KIR3DL1 is associated with increased risk of CIN.

The frequencies of the ligand groups (C1, C2 group, Bw4, and HLA-A3/11) did not differ between patients and controls (data not shown). However, Madeleine et al. [34] demonstrated an association between HLA-C subtypes and squamous cell cervical cancer, and Martin et al. [27] showed that C1 (asparagine at position 80) is over represented in women with cervical cancer. In Korean women, HLA-C is associated with HPV-cervical disease: HLA-C*03:03 confers susceptibility whereas HLA-C*01 has a protective effect [29].

In the present work, KIR2DL1-C2 was more frequent in patients (70%) than in controls (38.5%). The strength of NK inhibition varies according to the receptor and the ligand: KIR2DL1-C2 provides a stronger inhibition than other iKIR-HLA [32, 35]. The reduced resistance to viral infections among KIR2DL1-C2-positive individuals may result from the increased inhibition of NK cells. There were significant differences for the three pairs KIR2DL2/3-C1, KIR3DL1-Bw4, and KIR3DL2-A3/11, which displayed an increased frequency in patients. This combination of iKIR and ligands could be associated with persistent inflammatory reactions that play a role in carcinogenesis.

4.2. HLA and Its Association with HPV and CIN

HLA class I and class II proteins are central to host immune responses to viral infections and other pathogens. They are the most polymorphic genes in the human genome, and variations in the peptide binding groove of these proteins influence antigenic specificity. Numerous studies have evaluated the association of HLA with HPV infection and the importance of HLA in the pathogenesis of cervical neoplasia [36]. Similar to other studies [3739], there was no association between HLA specificities and HPV infection in admixed Brazilian women from the state of Paraná. However, the HLA-A*02-B*51 haplotype was associated with resistance to disease. Susceptibility to HPV infection or cervical cancer and precancerous lesion development was associated with the HLA class II: HLA-DRB1 alleles [34, 3951]; HLA-DQB1 alleles [34, 39, 46, 4951]; HLA-DPB1 alleles [51]; and classes I and II haplotypes [30, 34, 4042, 48, 52, 53]. Some alleles and haplotypes had a protective effect against the progression to infection and cancer [34, 3840, 43, 46, 47, 51, 54, 55]. In general, HLA-DQB1*03 increases and DRB1*13 decreases the risk of cervical cancer. In other Brazilian populations, Maciag et al. [44] found that HLA class II polymorphism was involved in genetic susceptibility to HPV infection and cervical cancer: DRB1*15:03, DRB1*04:05, and DQB1*06:02 alleles.

A genome-wide association study of 731.422 SNPs was performed in cervical cancer patients and controls [56]. Three independent loci in MHC region were associated with cervical cancer: the first is adjacent to the MHC class I polypeptide-related sequence A gene (MICA) (rs2516448; OR = 1.42; 95% CI = 1.31 to 1.54; ); the second is between HLA-DRB1 and HLA-DQA1 (rs9272143; OR = 0.67; 95% CI = 0.62 to 0.72; ); and the third is at HLA-DPB2 (rs3117027; OR = 1.25; 95% CI = 1.15 to 1.35; ). Previously reported associations of B*07:02 and DRB1*15:01-DQB1*06:02 with susceptibility to DRB1*13:01-DQA1*01:03-DQB1*06:03 with protection against cervical cancer were confirmed.

The variable results for the association between HLA and disease could be related to the differences in the distribution of HLA in the population; the disease phases (persistence or transitory HPV infection, intraepithelial neoplasia, and cancer); and HPV types. The effects of HLA polymorphisms on cervical carcinogenesis and their biological mechanisms remain unknown. Previous findings suggest a strong link between an inefficient immune response, particularly inefficient cell-mediated and innate immunity, both of which involve classes I and II HLA alleles, and susceptibility to HPV infection. HPV infections are more prevalent and more likely to persist in immunosuppressed individuals.

4.3. Cytokines and HPV

Accumulating epidemiological evidence suggests that polymorphisms in cytokine genes may be involved in the etiology of cervical carcinoma [6]. Th1 cytokines such IFN-γ and TNF-α can induce a cell-mediated immune response, whereas Th2 cytokines such as IL-6 and IL-10 induce predominantly a humoral immune response and immunomodulation of the cellular response. The Th2 cytokine profile is associated with progression to cervical cancer [57].

In the present study, we genotyped SNPs of TNF, IFNG, IL6, IL10, and TGFB1 which are multifunctional cytokine that have been implicated in inflammation, immunity, and cellular organization and have been proposed to play important roles in infection and cancer biology.

Susceptibility to infection was observed in patients with the IFNG +874A/A genotype, which characterized the low producer phenotype of IFN-γ and was more frequent among patients compared with controls. According to Telesheva et al. [58], the outcome of HPV infection is controlled by the interferon component of the immune response: a transitory course of HPV infection is characterized by increased levels of IFN-alpha and IFN-gamma, and persistent infection is related to decreased levels of IFN-alpha.

In the current study, the IL10 GCC/ACC and GCC/ATA genotype, which characterized the intermediate producer phenotype of IL-10, were less frequent in patients, suggesting protection against disease. The IL-10 high producer phenotypes was more frequent in patients, although this increased frequency was not significant and might be related to an immunosuppressive response and development of HPV-positive cervical cancer. Serum levels of IL-10 and its expression in tumor cells are elevated in patients with cervical cancer [56]. IL-10 produced by tumor macrophages induces a regulatory phenotype in T cells and an escape mechanism of the immune response that facilitates tumor growth [59]. The SNP IL10 −1082G>A was not associated with susceptibility to the development of cervical cancer or HPV infection [60].

TGF-β is well known for its antiproliferative effects; however, neoplastic cells often lose their sensitivity to TGF-β. Iancu et al. [61] showed that in human cervical cancer, disruption of the TGF-β signaling pathway might contribute to the malignant progression of cervical dysplasia. In the present study, SNPs of TGFB1 +869, +915 were not associated with HPV infection.

TNF-α, which is secreted mainly by activated macrophages, is an extraordinarily pleiotropic cytokine that has a central role in immune homeostasis, inflammation, and host defense and could be involved in protection against HPV infection by modulating viral replication. Dysregulated TNF expression within the tumor microenvironment appears to favor malignant cell tissue invasion, migration, and ultimately metastasis [62]. Our findings are similar to those reported by Wang et al. [63], who demonstrated that there is no significant association between the TNF −308G>A and HPV infection or cervical cancer. However, our findings differ from those reported for the Argentina population, among whom the high producer allele TNFA −307A was associated with an increased risk for the development of cervical cancer [64].

IL6 encodes a cytokine that plays important roles in the risk for cervical carcinogenesis. In the present study, there was no significant association between IL6 and HPV-related cervical disease. However, a previous report has shown a significant association between the IL6-rs2069837 SNP and an increased risk of cervical cancer [65].

5. Conclusion

The genetic variability of the host contributes to the risk of HPV-related cervical disease. KIR genes were not associated with HPV, although some pairs of iKIR-ligands were more frequent in patients, suggesting that NK cells play a role in detrimental chronic inflammatory conditions and in carcinogenesis. HLA was associated with HPV and participated in the immune response, although its function in carcinogenesis remains unclear. The polymorphic INFG and IL10 genes were associated with the outcome of HPV infection and might be indicative of impaired or invalid immune responses in patients with advanced stage lesions. Additional studies of the immune response to HPV are needed to better define the risk of developing malignant diseases associated with HPV infection.

Acknowledgments

The authors are grateful to all of the volunteers who participated in this study and to both Fabiano Cavalcante de Melo and Marco Antonio Braga for performing the technical procedures.