BioMed Research International

Volume 2013 (2013), Article ID 146079, 11 pages

http://dx.doi.org/10.1155/2013/146079

Research Article

The Association of the Immune Response Genes to Human Papillomavirus-Related Cervical Disease in a Brazilian Population

Amanda Vansan Marangon,¹ Gláucia Andreia Soares Guelsin,¹ Jeane Eliete Laguila Visentainer,² Sueli Donizete Borelli,² Maria Angélica Ehara Watanabe,³ Márcia Edilaine Lopes Consolaro,⁴ Katiany Rizzieri Caleffi-Ferracioli,⁴ Cristiane Conceição Chagas Rudnick,⁵ and Ana Maria Sell²

¹Campinas State University, Rua Carlos Chagas 480, Barão Geraldo, 13083878 Campinas, SP, Brazil

²Basic Health Sciences Department, Maringa State University, Avenue Colombo 5790, 87020900 Maringá, PR, Brazil

³Pathology Sciences Department, Londrina State University, Rodovia Celso Garcia Cid (PR 445), 86051-970 Londrina, PR, Brazil

⁴Department of Clinical Analysis and Biomedicine, Maringa State University, Avenue Colombo 5790, 87020900 Maringá, PR, Brazil

⁵Program of Biosciences Applied to Pharmacy, Department of Clinical Analysis and Biomedicine, Maringa State University, Avenue Colombo 5790, 87020900 Maringá, PR, Brazil

Received 29 March 2013; Revised 5 June 2013; Accepted 15 June 2013

Academic Editor: Enrique Medina-Acosta

Copyright © 2013 Amanda Vansan Marangon et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

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.

Table 1: Cytokine gene SNPs interrogated in this study.

Cytokine gene
Gene chromosome location
SNP designation in the kit
dbSNP-ID
Location

TNF
6p21.3
−308 G/A
rs1800629
Promoter

IFNG
12q14
+874 T/A
rs2430561
Intron

IL6
7p21
−174 G/C
rs1800795
Promoter

IL10
1q31-q32
−1082 A/G
rs1800896
Promoter

−819 C/T
rs1800871
Promoter

−592 C/A
rs1800872
Promoter

TGFB1
19q13.1
−509 T/C (or 869 T/C)
rs1800469
Promoter

+915 G/C
rs1800471
Exon 1

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].

Table 2: Frequencies of KIR genes in HPV patients and controls.

KIR genes
HPV patients
Controls

%

%

2DL1
68
95.8
116
98.3

2DL2
35
49.3
55
46.6

2DL3
63
88.7
104
88.1

2DL5
43
60.6
62
52.5

2DP1
68
95.8
115
97.4

2DS1
30
42.2
50
42.4

2DS2
37
52.1
57
48.3

2DS3
26
36.6
36
30.5

2DS4
65
91.5
108
91.5

2DS5
28
39.4
40
33.9

3DL1
65
91.5
109
92.4

3DS1
31
43.6
48
40.7

2DL4, 3DL2,  3DL3, and  3DP1
71
100
118
100

KIR gene frequencies were similar in both the groups ().

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).

Table 3: Distribution of KIR and HLA ligands in HPV patients and controls.

KIR and HLA ligands
HPV patients
Controls

%

%

2DL1-C2
46
67.65
78
67.24

2DL1 without C2
22
32.35
38
32.75

2DL2-C1
28
80.0
43
78.18

2DL2 without C1
7
20.0
12
21.81

2DL3-C1
54
85.71
87
83.65

2DL3 without C1
9
14.29
17
16.34

3DL1-Bw4
44
67.69
73
66.97

3DL1 without Bw4
21
32.31
36
33.02

2DS1-C2
22
73.33
30
60.00

2DS1 without C2
8
26.67
20
40.00

2DS2-C1
30
81.08
45
78.94

2DS2 without C1
7
18.92
12
21.05

3DS1-Bw4
21
67.74
30
62.50

3DS1 without Bw4
10
32.26
18
37.50

Bw4: HLA-A*23, 24, and 32; HLA-B*08, 13, 27, 37, 44, 51, 52, 53, 57, and 58.

Group C1: HLA-C*01, 03, 07, 08, 12, 14, and 16.

Group C2: HLA-C*02, 04, 05, 06, 07, 15, 17, and 18.

Difference was not observed ().

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.

Table 4: Combinations of inhibitory KIR-HLA pairs and their frequencies in HPV and control Brazilian women from Paraná, Southern Brazil.

Number of pairs
KIR-HLA
HPV patients

(%)
Control

(%)

1 pair
2DL2/3-C1
3 (30.0)
7 (61.5)

2DL1-C2
7 (70.0)
4 (38.5)

2 pairs
2DL2/3-C1, 3DL1-Bw4
12 (46.2)
16 (34.5)

2DL2/3-C1, 2DL1-C2
7 (26.9)
14 (29.0)

2DL1-C2, 3DL1-Bw4
4 (15.4)
11 (23.6)

2DL1-C2, 3DL2-A3/11
2 (7.7)
4 (9.1)

2DL2/3-C1, 3DL2-A3/11
1 (3.9)
2 (3.7)

3 pairs
2DL1-C2, 2DL2/3-C1, 3DL1-Bw4
14 (51.9)
25 (76.3)

2DL2/3-C1, 3DL1-Bw4, 3DL2-A3/11a
7 (25.9)
2 (5.3)

2DL1-C2, 2DL2/3-C1, 3DL2-A3/11
6 (23.0)
5 (15.8)

2DL1-C2, 3DL1-Bw4, 3DL2-A3/11
0 (0)
1 (2.6)

4 pairs
2DL1-C2, 2DL2/3-C1, 3DL1-Bw4, 3DL2-A3/11
8 (100)
12 (100)

; OR = 3.42; 95% CI = 2.45–18.22.

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.

Figure 1: KIR genotypes frequencies in HPV patients and controls from Paraná, Southern Brazil.

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

Table 5: HLA allele frequencies in HPV patients and control groups.

HLA-A allele types
HLA-B allele types
HLA-DRB1 allele types

HPV patients
Controls

HPV patients
Controls

HPV patients
Controls

A*01
14 (9.0)
31 (10.3)
B*07
8 (5.1)
17 (5.7)
DRB1*01
20 (12.8)
51 (17.0)

A*02
37 (23.7)
67 (22.3)
B*08
11 (7.0)
15 (5)
DRB1*03
17 (10.9)
23 (7.7)

A*03
21 (13.5)
29 (9.7)
B*13
3 (1.9)
6 (2)
DRB1*04
8 (5.1)
25 (8.3)

A*11
11 (7.0)
21 (7.0)
B*14
9 (5.8)
14 (4.7)
DRB1*07
21 (13.5)
33 (11.0)

A*23
7 (4.5)
13 (4.3)
B*15
22 (14.1)
26 (8.7)
DRB1*08
10 (6.4)
17 (5.7)

A*24
20 (12.8)
31 (10.3)
B*18
6 (3.8)
21 (7)
DRB1*09
3 (19.0)
4 (1.3)

A*25
2 (1.2)
7 (2.3)
B*27
6 (3.8)
10 (3.3)
DRB1*10
6 (3.8)
7 (2.3)

A*26
2 (1.2)
13 (4.3)
B*35
18 (11.5)
39 (13)
DRB1*11
27 (17.3)
34 (11.3)

A*29
8 (5.1)
14 (4.7)
B*37
2 (1.3)
4 (1.3)
DRB1*12
4 (2.6)
6 (2.0)

A*30
12 (7.7)
15 (5.0)
B*38
2 (1.3)
3 (1)
DRB1*13
17 (10.9)
47 (15.7)

A*31
6 (3.8)
13 (4.3)
B*39
5 (3.2)
10 (3.3)
DRB1*14
8 (5.1)
14 (4.7)

A*32
2 (1.3)
9 (3.0)
B*40
5 (3.2)
15 (5)
DRB1*15
8 (5.1)
25 (8.3)

A*33
3 (1.9)
14 (4.7)
B*41
1 (0.6)
6 (2)
DRB1*16
7 (4.5)
14 (4.7)

A*34
1 (0.6)
1 (0.3)
B*42
1 (0.6)
10 (0.3)

A*66
1 (0.6)
3 (1.0)
B*44
21 (13.5)
34 (11.3)

A*68
9 (5.8)
16 (5.3)
B*45
4 (2.6)
5 (1.7)

B*48
1 (0.6)
1 (0.3)

B*49
4 (2.6)
9 (3.0)

B*50
2 (1.3)
8 (2.7)

B*51
12 (7.7)
25 (8.3)

B*52
4 (2.6)
5 (1.7)

B*53
1 (0.6)
10 (3.3)

B*55
1 (0.6)
3 (1.0)

B*57
4 (2.6)
9 (3.0)

B*58
3 (1.9)
4 (1.3)

: number of alleles; : number of individuals; %: alleles frequencies.

Difference was not observed between both groups ().

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).

Table 6: HLA haplotype frequencies with significant differences between HPV patients and controls.

Haplotypes
Patients (hf)
Controls (hf)

OR
CI

HLA-A*02–HLA-B*51
1 (0.006)
16 (0.052)
0.006
0.1186
0.015–0.8717

HLA-A*03–HLA-DRB1*11
7 (0.043)
1 (0.003)
0.002
150.828
1.712–115.230

HLA-B*14–HLA-DRB1*13
5 (0.032)
1 (0.003)
0.017
108.906
1.146–85.504

HLA-B*15–HLA-DRB1*07
5 (0.029)
1 (0.003)
0.017
108.906
1.146–85.505

HLA-B*15–HLA-DRB1*11
7 (0.048)
1 (0.003)
0.002
150.828
1.712–115.230

HLA-B*44–HLA-DRB1*01
6 (0.036)
2 (0.005)
0.018
49.028
1.188–29.885

HLA-B*44–HLA-DRB1*11
7 (0.048)
1 (0.003)
0.002
150.828
1.712–115.230

: haplotype numbers; hf: haplotype frequencies (%); : value; OR: odds ratio; CI (95%): 95% confidence interval.

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].

Table 7: Cytokines alleles and genotypes frequencies in HPV patients and controls.

Cytokine alleles and genotypes
HPV patients

Controls

Cytokine alleles and genotypes
HPV patients

Controls

TNF  −308

IL10  −1082

    G
139 (88.0)
174 (86.1)
    G
50 (32.1)
72 (36)

    A
19 (12.0)
28 (13.9)
    A
106 (68)
128 (64)

    G/G
62 (78.5)
73 (72.3)
    G/G
9 (11.5)
12 (12)

    G/A
15 (19.0)
28 (27.7)
    G/A
32 (41.0)
48 (48)

    A/A
2 (2.5)
0 (0)
    A/A
37 (47.4)
40 (40)

INFG  +874

IL10  −819

    T
69 (44.2)
88 (44.9)
    C
98 (62.8)
129 (64.5)

    A
87 (55.8)
108 (55.1)
    T
58 (37.2)
71 (35.5)

    T/T
19 (24.4)
19 (19.4)
    C/C
31 (39.8)
42 (42)

    T/A
31 (39.7)
50 (51.0)
    C/T
36 (46.2)
45 (45)

    A/A
28 (35.9)
29 (29.6)
    T/T
11 (14.1)
13 (13)

IL6  −174

IL10  −592

    G
116 (73.4)
132 (65.4)
    C
98 (62.8)
129 (64.5)

    C
42 (26.6)
70 (34.7)
    A
58 (37.2)
71 (35.5)

    G/G
44 (55.7)
45 (44.6)
    C/C
31 (39.8)
42 (42)

    G/C
28 (35.4)
42 (41.6)
    C/A
36 (46.2)
45 (45)

    C/C
7 (8.9)
14 (13.9)
    A/A
11 (14.1)
13 (13)

TGFB1  +869

TGFB1  +915

    T
89 (56.3)
108 (54)
    G
141 (89.2)
190 (95)

    C
69 (43.7)
92 (46)
    C
17 (10.7)
10 (5)

    T/T
26 (32.9)
25 (25)
    G/G
63 (79.8)
90 (90)

    T/C
37 (46.8)
58 (58)
    G/C
15 (19.0)
10 (10)

    C/C
16 (20.3)
17 (17)
    C/C
1 (1.3)
0 (0)

: number of observed alleles and genotypes; %: allele and genotype frequencies.

Difference was not observed between both groups ().

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.

Table 8: Expected phenotype frequencies according to genotypes for the cytokines TNF-, IFN-, IL-6, IL-10, and TGF-β1.

Phenotypes
Genotypes
Patients ()  

(%)
Controls ()  

(%)

TNF

Low
    G/G
62 (78.48)
73 (72.28)

High
    G/A
17 (21.52)
28 (27.72)

    A/A

IL6

High
    G/G
72 (91.14)
87 (86.14)

    G/C

Low
    C/C
7 (8.86)
14 (13.86)

INFG

High
    T/T
19 (24.36)
19 (19.39)

Intermediate
    T/A
31 (39.74)
50 (51.02)

Lowb
    A/A
28 (35.90)
29 (29.59)

IL10

High
    GCC/GCC
16 (20.51)
12 (12.00)

Intermediatea
    GCC/ACC
25 (32.05)
48 (48.00)

    GCC/ATA

Low
    ACC/ACC
37 (47.44)
40 (40.00)

    ACC/ATA

    ATA/ATA

TGFB1

High
    T/T  G/G
52 (65.82)
78 (78.00)

    T/C  G/G

Intermediate
    T/C  G/C
22 (27.85)
17 (17.00)

    C/C  G/G

    T/T  G/C

Low
    C/C  G/C
5 (6.33)
5 (5.00)

    C/C  C/C

    T/T  C/C

    T/C  C/C

: number of excepted phenotypes according to genotypes.

%: frequencies.

a; OR = 0.1607; 95% CI = 0.276–0.947.

b; OR = 1.81; 95% CI = 1.178–4.604.

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 [37–39], 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, 39–51]; HLA-DQB1 alleles [34, 39, 46, 49–51]; HLA-DPB1 alleles [51]; and classes I and II haplotypes [30, 34, 40–42, 48, 52, 53]. Some alleles and haplotypes had a protective effect against the progression to infection and cancer [34, 38–40, 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.

References

M. Schiffman, P. E. Castle, J. Jeronimo, A. C. Rodriguez, and S. Wacholder, “Human papillomavirus and cervical cancer,” The Lancet, vol. 370, no. 9590, pp. 890–907, 2007.

D. Subramanya and P. D. Grivas, “HPV and cervical cancer: updates on an established relationship,” Postgraduate Medicine, vol. 120, no. 4, pp. 7–13, 2008.

C. D. Woodworth, “HPV innate immunity,” Frontiers in Bioscience, vol. 7, pp. d2058–d2071, 2002.

M. H. Hibma, “The immune response to papillomavirus during infection persistence and regression,” The Open Virology Journal, vol. 6, pp. 241–248, 2012.

T. Sasagawa, H. Takagi, and S. Makinoda, “Immune responses against human papillomavirus (HPV) infection and evasion of host defense in cervical cancer,” Journal of Infection and Chemotherapy, vol. 18, no. 6, pp. 807–815, 2012.

A. Sharma, M. Rajappa, A. Satyam, and M. Sharma, “Cytokines ( and ) in patients with advanced cervical cancer undergoing neoadjuvant chemoradiation correlation with treatment response,” International Journal of Gynecological Cancer, vol. 19, no. 7, pp. 1269–1275, 2009.

S. Rajagopalan and E. O. Long, “Understanding how combinations of HLA and KIR genes influence disease,” Journal of Experimental Medicine, vol. 201, no. 7, pp. 1025–1029, 2005.

R. N. Germain and D. H. Margulies, “The biochemistry and cell biology and antigen processing and presentation,” Annual Review of Immunology, vol. 11, pp. 403–450, 1993.

D. M. Turner, D. M. Williams, D. Sankaran, M. Lazarus, P. J. Sinnott, and I. V. Hutchinson, “An investigation of polymorphism in the interleukin-10 gene promoter,” European Journal of Immunogenetics, vol. 24, no. 1, pp. 1–8, 1997.

K. M. Kroeger, K. S. Carville, and L. J. Abraham, “The -308 tumor necrosis factor-α promoter polymorphism effects transcription,” Molecular Immunology, vol. 34, no. 5, pp. 391–399, 1997.

D. Fishman, G. Faulds, R. Jeffey et al., “The effect of novel polymorphisms in the interleukin-6 (IL-6) gene on IL-6 transcription and plasma IL-6 levels, and an association with systemic-onset juvenile chronic arthritis,” Journal of Clinical Investigation, vol. 102, no. 7, pp. 1369–1376, 1998.

V. Pravica, A. Asderakis, C. Perrey, A. Hajeer, P. J. Sinnott, and I. V. Hutchinson, “In vitro production of IFN-γ correlates with CA repeat polymorphism in the human IFN-γ gene,” European Journal of Immunogenetics, vol. 26, no. 1, pp. 1–3, 1999.

M. R. Awad, A. El-Gamel, P. Hasleton, D. M. Turner, P. J. Sinnott, and I. V. Hutchinson, “Genotypic variation in the transforming growth factor-β1 gene: association with transforming growth factor-β1 production, fibrotic lung disease, and graft fibrosis after lung transplantation,” Transplantation, vol. 66, no. 8, pp. 1014–1020, 1998.

C. M. Probst, E. P. Bompeixe, N. F. Pereira et al., “HLA polymorphism and evaluation of European, African, and Amerindian contribution to the white and mulatto populations from Paraná, Brazil,” Human Biology, vol. 72, no. 4, pp. 597–617, 2000.

D. Solomon and R. Nayar, Bethesda System for Cervical-Vaginal Cytology, Revinter, Rio de Janeiro, Brazil, 2005.

D. C. Szurkus and T. A. Harrison, “Loop excision for high-grade squamous intraepithelial lesion on cytology: correlation with colposcopic and histologic findings,” American Journal of Obstetrics and Gynecology, vol. 188, no. 5, pp. 1180–1182, 2003.

W. Qu, G. Jiang, Y. Cruz et al., “PCR detection of human papillomavirus: comparison between MY09/MY11 and GP5+/GP6+ primer systems,” Journal of Clinical Microbiology, vol. 35, no. 6, pp. 1304–1310, 1997.

E. Santiago, L. Camacho, M. L. Junquera, and F. Vázquez, “Full HPV typing by a single restriction enzyme,” Journal of Clinical Virology, vol. 37, no. 1, pp. 38–46, 2006.

M. P. Martin, G. Nelson, J.-H. Lee et al., “Cutting edge: susceptibility to psoriatic arthritis: influence of activating killer Ig-like receptor genes in the absence of specific HLA-C alleles,” Journal of Immunology, vol. 169, no. 6, pp. 2818–2822, 2002.

C. C. C. Rudnick, G. A. S. Guelsin, A. V. Marangon, D. S. A. Franceschi, A. M. Sell, and J. E. L. Visentainer, “Methodology optimization for KIR genotyping,” Jornal Brasileiro de Patologia e Medicina Laboratorial, vol. 46, no. 3, pp. 215–224, 2010.

M. Carrington and P. Norman, The KIR Gene Cluster, National Library of Medicine (US), NCBI, Bethesda, Md, USA, 2003.

E. W. Petersdorf, “Risk assessment in haematopoietic stem cell transplantation: histocompatibility,” Best Practice and Research: Clinical Haematology, vol. 20, no. 2, pp. 155–170, 2007.

S. W. Guo and E. A. Thompson, “Performing the exact test of Hardy-Weinberg proportion for multiple alleles,” Biometrics, vol. 48, no. 2, pp. 361–372, 1992.

C. C. C. Rudnick, D. S. A. Franceschi, A. V. Marangon, G. A. S. Guelsin, A. M. Sell, and J. E. L. Visentainer, “Killer cell immunoglobulin-like receptor gene diversity in a Southern Brazilian population from the state of Paraná,” Human Immunology, vol. 69, no. 12, pp. 872–876, 2008.

J. E. L. Visentainer, A. M. Sell, G. C. Da Silva et al., “TNF, IFNG, IL6, IL10 and TGFB1 gene polymorphisms in South and Southeast Brazil,” International Journal of Immunogenetics, vol. 35, no. 4-5, pp. 287–293, 2008.

E. M. Burd, “Human papillomavirus and cervical cancer,” Clinical Microbiology Reviews, vol. 16, no. 1, pp. 1–17, 2003.

M. P. Martin, I. B. Borecki, Z. Zhang et al., “HLA-Cw group 1 ligands for KIR increase susceptibility to invasive cervical cancer,” Immunogenetics, vol. 62, no. 11-12, pp. 761–765, 2010.

M. Han, M. Fallena, Y. Guo, and P. Stastny, “Natural killer cell crossmatch: functional analysis of inhibitory killer immunoglobulin-like receptors and their HLA ligands,” Human Immunology, vol. 68, no. 6, pp. 507–513, 2007.

M. J. Song, C. W. Lee, J. H. Kim, et al., “Association of KIR genes and HLA-C alleles with HPV-related uterine cervical disease in Korean women,” Tissue Antigens, vol. 81, no. 03, pp. 164–170, 2013.

S.-I. Yusa, T. L. Catina, and K. S. Campbell, “KIR2DL5 can inhibit human NK cell activation via recruitment of Src homology region 2-containing protein tyrosine phosphatase-2 (SHP-2),” Journal of Immunology, vol. 172, no. 12, pp. 7385–7392, 2004.

Z.-M. Lu, Y.-L. Jiao, Z.-L. Feng et al., “Polymorphisms of killer cell immunoglobulin-like receptor gene: possible association with susceptibility to or clearance of hepatitis B virus infection in Chinese Han population,” Croatian Medical Journal, vol. 48, no. 6, pp. 800–806, 2007.

M. Carrington, S. Wang, M. P. Martin et al., “Hierarchy of resistance to cervical neoplasia mediated by combinations of killer immunoglobulin-like receptor and human leukocyte antigen loci,” Journal of Experimental Medicine, vol. 201, no. 7, pp. 1069–1075, 2005.

L. Arnheim, J. Dillner, and C. B. Sanjeevi, “A population-based cohort study of KIR genes and genotypes in relation to cervical intraepithelial neoplasia,” Tissue Antigens, vol. 65, no. 3, pp. 252–259, 2005.

M. M. Madeleine, B. Brumback, K. L. Cushing-Haugen et al., “Human leukocyte antigen class II and cervical cancer risk: a population-based study,” Journal of Infectious Diseases, vol. 186, no. 11, pp. 1565–1574, 2002.

P. Parham, “MHC class I molecules and KIRS in human history, health and survival,” Nature Reviews Immunology, vol. 5, no. 3, pp. 201–214, 2005.

K. Chattopadhyay, “A comprehensive review on host genetic susceptibility to human papillomavirus infection and progression to cervical cancer,” Indian Journal of Human Genetics, vol. 17, no. 3, pp. 132–144, 2011.

Y.-C. Yang, T.-Y. Chang, Y.-J. Lee et al., “HLA-DRB1 alleles and cervical squamous cell carcinoma: experimental study and meta-analysis,” Human Immunology, vol. 67, no. 4-5, pp. 331–340, 2006.

S. Ades, A. Koushik, E. Duarte-Franco et al., “Selected class I and class II HLA alleles and haplotypes and risk of high-grade cervical intraepithelial neoplasia,” International Journal of Cancer, vol. 122, no. 12, pp. 2820–2826, 2008.

A. Guzalinuer, A. Mihrinsa, S.-Q. Zhang, H. Li, N. Gulishare, and G.-Q. Zhang, “Association between HPV infection and HLA-DQB1 alleles polymorphism in the cervical carcinogenesis in Uyghur women in southern Xinjiang,” Zhonghua Zhong Liu Za Zhi, vol. 32, no. 7, pp. 492–496, 2010.

I. Kohaar, S. Hussain, N. Thakur et al., “Association between human leukocyte antigen class II alleles and human papillomavirus-mediated cervical cancer in Indian women,” Human Immunology, vol. 70, no. 4, pp. 222–229, 2009.

D. M. Hernández-Hernández, R. M. Cerda-Flores, T. Juárez-Cedillo et al., “Human leukocyte antigens I and II haplotypes associated with human papillomavirus 16-positive invasive cervical cancer in Mexican women,” International Journal of Gynecological Cancer, vol. 19, no. 6, pp. 1099–1106, 2009.

Y. Ben Othmane, E. Ghazouani, A. Mezlini, et al., “HLA class II susceptibility to cervical cancer among Tunisian women,” Bull Cancer, vol. 99, no. 9, pp. 81–86, 2012.

F. A. Castro, K. Haimila, I. Sareneva et al., “Association of HLA-DRB1, interleukin-6 and cyclin D1 polymorphisms with cervical cancer in the Swedish population—a candidate gene approach,” International Journal of Cancer, vol. 125, no. 8, pp. 1851–1858, 2009.

P. C. Maciag, N. F. Schlecht, P. S. A. Souza, E. L. Franco, L. L. Villa, and M. L. Petzl-Erler, “Major histocompatibility complex class II polymorphisms and risk of cervical cancer and human papillomavirus infection in Brazilian women,” Cancer Epidemiology Biomarkers & Prevention, vol. 9, no. 11, pp. 1183–1191, 2000.

J. Cervantes, C. Lema, L. V. Hurtado et al., “HLA-DRB1*1602 allele is positively associated with HPV cervical infection in Bolivian Andean women,” Human Immunology, vol. 64, no. 9, pp. 890–895, 2003.

K. Eiguchi, S. Tatti, L. V. Alonio et al., “Association of DRB1 and DQB1 HLA class II polymorphisms in high-grade and neoplastic cervical lesions of women from Argentina,” Journal of Lower Genital Tract Disease, vol. 12, no. 4, pp. 262–268, 2008.

C. Alaez-Verson, J. Berumen-Campos, A. Munguía-Saldaña et al., “HPV-16 and HLA-DRB1 alleles are associated with cervical carcinoma in Mexican Mestizo women,” Archives of Medical Research, vol. 42, no. 5, pp. 421–425, 2011.

L.-C. Chuang, C.-Y. Hu, H.-C. Chen et al., “Associations of human leukocyte antigen class II genotypes with human papillomavirus 18 infection and cervical intraepithelial neoplasia risk,” Cancer, vol. 118, no. 1, pp. 223–231, 2012.

D. D. Dao, C. H. Sierra-Torres, S. C. Robazetti et al., “HLA-DQB1 and cervical cancer in Venezuelan women,” Gynecologic Oncology, vol. 96, no. 2, pp. 349–354, 2005.

C. Lema, A. L. Fuessel-Haws, L. R. Lewis et al., “Association between HLA-DQB1 and cervical dysplasia in Vietnamese women,” International Journal of Gynecological Cancer, vol. 16, no. 3, pp. 1269–1277, 2006.

J. Liang, A. Xu, Y. Xie, A. O. Awonuga, and Z. Lin, “Some but not all of HLA-II alleles are associated with cervical cancer in Chinese women,” Cancer Genetics and Cytogenetics, vol. 187, no. 2, pp. 95–100, 2008.

Y. Wu, B. Liu, W. Lin et al., “Human leukocyte antigen class II alleles and risk of cervical cancer in China,” Human Immunology, vol. 68, no. 3, pp. 192–200, 2007.

P. S. De Araujo Souza and L. L. Villa, “Genetic susceptibility to infection with human papillomavirus and development of cervical cancer in women in Brazil,” Mutation Research, vol. 544, no. 2-3, pp. 375–383, 2003.

K. Matsumoto, H. Maeda, A. Oki et al., “HLA class II DRB1*1302 allele protects against progression to cervical intraepithelial neoplasia grade 3: a multicenter prospective cohort study,” International Journal of Gynecological Cancer, vol. 22, no. 3, pp. 471–478, 2012.

P. K. S. Chan, J. L. K. Cheung, T.-H. Cheung et al., “HLA-DQB1 polymorphisms and risk for cervical cancer: a case-control study in a southern Chinese population,” Gynecologic Oncology, vol. 105, no. 3, pp. 736–741, 2007.

D. Chen, I. Juko-Pecirep, J. Hammer, et al., “Genome-wide association study of susceptibility loci for cervical cancer,” Journal of the National Cancer Institute, vol. 105, no. 9, pp. 624–633, 2013.

Q. Feng, H. Wei, J. Morihara, et al., “Th2 type inflammation promotes the gradual progression of HPV-infected cervical cells to cervical carcinoma,” Gynecologic Oncology, vol. 127, no. 2, pp. 412–419, 2012.

L. F. Telesheva, V. F. Dolgushina, O. S. Abramovskikh, et al., “Cytokine status of cervical mucus in women with transitory and persistent course of papillomavirus infection,” Zhurnal Mikrobiologii, Epidemiologii, i Immunobiologii, vol. 4, pp. 118–121, 2012.

A. Bolpetti, J. S. Silva, L. L. Villa, and A. P. Lepique, “Interleukin-10 production by tumor infiltrating macrophages plays a role in Human Papillomavirus 16 tumor growth,” BMC Immunology, vol. 11, article 27, 2010.

G. Barbisan, L. O. Pérez, A. Contreras, et al., “TNF- and IL-10 promoter polymorphisms, HPV infection, and cervical cancer risk,” Tumor Biology, vol. 33, no. 5, pp. 1549–1556, 2012.

I. V. Iancu, A. Botezatu, C. D. Goia-Ruşanu et al., “TGF-beta signalling pathway factors in HPV-induced cervical lesions,” Roumanian Archives of Microbiology and Immunology, vol. 69, no. 3, pp. 113–118, 2010.

S. Mocellin and D. Nitti, “TNF and cancer: the two sides of the coin,” Frontiers in Bioscience, vol. 13, no. 7, pp. 2774–2783, 2008.

N. Wang, D. Yin, S. Zhang et al., “TNF-alpha rs1800629 polymorphism is not associated with HPV infection or cervical cancer in the Chinese population,” PLoS ONE, vol. 7, no. 9, Article ID e44952, 2012.

I. Badano, S. M. Stietz, T. G. Schurr et al., “Analysis of TNFα promoter SNPs and the risk of cervical cancer in urban populations of Posadas (Misiones, Argentina),” Journal of Clinical Virology, vol. 53, no. 1, pp. 54–59, 2012.

T. Y. Shi, M. L. Zhu, J. He, et al., “Polymorphisms of the Interleukin 6 gene contribute to cervical cancer susceptibility in Eastern Chinese women,” Human Genetics, vol. 132, no. 3, pp. 301–312, 2013.