Mediators of Inflammation
Volume 2018 (2018), Article ID 9127146, 11 pages
https://doi.org/10.1155/2018/9127146

Research Article

The Associations between Toll-Like Receptor 9 Gene Polymorphisms and Cervical Cancer Susceptibility

Sijuan Tian,¹ Liping Zhang,² Ting Yang,¹ Xing Wei,¹ Li Zhang,¹ Yang Yu,¹ Yang Li,¹ Di Cao,¹ and Xiaofeng Yang¹

¹Department of Obstetrics and Gynecology, The First Affiliated Hospital of Xi’an Jiaotong University, 277 West Yanta Road, Xi’an, Shaanxi 710061, China
²Department of Obstetrics and Gynecology, Taihe Hospital, Hubei University of Medicine, Shiyan, Hubei 442000, China

Correspondence should be addressed to Xiaofeng Yang; yxf73@163.com

Received 28 January 2018; Accepted 27 June 2018; Published 25 July 2018

Academic Editor: Giuseppe Valacchi

Copyright © 2018 Sijuan Tian 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.

This meta-analysis systematically reviews the association between Toll-like receptor 9 polymorphisms and the risk of cervical cancer. Case-control studies focused on the association were collected from the PubMed, Web of Science, Cochrane Library, Embase, MEDLINE, CNKI, VIP, and Wanfang databases from inception to July 2017. We screened the studies and assessed the methodological quality of the included studies and extracted data. A meta-analysis was performed using RevMan 5.3 and Stata 12.0 software. Pooled odds ratios and 95% confidence intervals were employed to evaluate the strength of the associations between Toll-like receptor 9 polymorphisms and cervical cancer risk. A total of 9 studies comprising 3331 cervical cancer patients and 4109 healthy controls met the inclusion criteria. Of these, 8 studies contained information about G2848A (rs352140) and 4 studies contained information about −1486T/C (rs187084). Our results revealed that the associations between rs187084 and cervical cancer risk in the dominant model () and heterozygous model () were significant, with 1.30- and 1.32-fold increases in susceptibility, respectively, compared to that in the wild-type model. However, rs352140 was not related to cervical cancer regardless of whether the subgroup analysis was conducted (). In conclusion, there is a significant correlation between rs187084 and cervical cancer risk with the minor C allele increasing the risk of occurrence of cervical cancer. However, rs352140 is not associated with the occurrence of cervical cancer.

1. Introduction

Cervical cancer is the fourth most common cancer in women in terms of both incidence and mortality worldwide [1]. According to global cancer statistics, there were approximately 527,600 new cervical cancer cases and 265,700 deaths in 2012. In low-income countries, cervical cancer ranks as second in incidence and is the third leading cause of cancer-associated death among women. In China, the estimated new cancer cases and deaths were 98,900 and 30,500, respectively, in 2015, presenting an upward trend [2, 3]. Human papillomavirus (HPV) was identified as a principal cause of cervical cancer [4, 5]. However, 90% of HPV genital infections can be spontaneously cleared, and few of these progress to cervical cancer [6], suggesting that other pathogeneses and aetiologies might contribute to cervical carcinogenesis.

Toll-like receptors (TLRs) are transmembrane proteins that recognize pathogen-associated molecular patterns (PAMPs), conserved structural motifs in bacteria, fungi, and viruses [7]. TLRs initiate the innate immune response and further modulate acquired immune response, playing an important role in inflammation and carcinogenesis. TLR9 recognizes nonmethylated CpG islands in viral DNA and activates the immune system. In the past few years, evidence suggests that TLR9 expression increases according to the histopathological grade of the cervical pathological process and the HPV E6 and E7 oncoproteins deregulate the expression and function of TLR9 [8, 9].

It has been noted that the TLR9 gene is polymorphic and associated with various cancers, including cervical cancer, prostate cancer, oesophageal cancer, gastric cancer, breast cancer, colorectal carcinoma, and lymphoma [10]. Many single-nucleotide polymorphisms (SNPs) of TLR9 have been studied, including rs352140, rs187084, rs352149, rs445676, and rs5743836, according to the NCBI database (https://www.ncbi.nlm.nih.gov/snp). Among these SNPs, two prominent variants, namely, rs352140 (G2848A) and rs187084 (−1486T/C), are frequently observed to be related to cervical cancer susceptibility. However, the results are contentious. Therefore, we conducted this meta-analysis to estimate the association between cervical cancer risk and the most concerning two SNPs of TLR9, rs352140 and rs187084.

2. Materials and Methods

2.1. Literature Search Strategy

We reviewed the PubMed, Web of Science, the Cochrane Library, Embase, MEDLINE, CNKI, VIP, and Wanfang databases systematically and comprehensively. The search terms are the following: Toll-Like Receptors [Mesh] or TLR, polymorphism/variant/mutation/SNP, Uterine Cervical Neoplasm [Mesh]/cervix cancer/cervical cancer, and the combinations of these. Additionally, we searched the reference lists of all identified articles manually for more studies.

2.2. Inclusion and Exclusion Criteria

Studies included needed to meet the following criteria: (1) a focus on the association between the TLR9 gene polymorphisms (rs352140 and rs187084) and the risk of cervical cancer, (2) human study subjects, (3) case-control studies, (4) available and sufficient genotype distribution of data to calculate odds ratios (ORs) and corresponding 95% confidence intervals (CIs), and (5) diagnoses based on cervical biopsy pathology. Additionally, if there were duplicate studies, we made sure that the most recent or the most complete one was included. If it did not satisfy the criteria above, the article was excluded.

2.3. Data Extraction and Synthesis

Two investigators extracted relevant data from all the eligible studies independently. A third reviewer was invited to participate in the work when some disagreement occurred; consensus was ultimately reached by discussion. These characteristics were collected from each study: the first author, publication year, race, total numbers of cases and controls, study design, source of controls, genotyping method, and evidence of HWE in controls.

2.4. Quality Assessment

The quality of included studies was assessed using the Newcastle-Ottawa scale, including three categories: selection, comparability, and outcome. Additional eight items were used to assess the methodology of each qualified study. The highest score was 9. Studies with a score of more than 7 were considered as high quality. A study awarded a score of 0–3, 4–6, or 7–9 was considered as a low-, moderate-, or high-quality study, respectively.

2.5. Statistical Analysis

The odds ratios (ORs) and 95% confidence interval (CIs) were applied to assess the strength of the correlation between SNPs and cervical cancer susceptibility. A -test revealed statistical significance when . ² and statistics were employed to detect heterogeneity among different studies. There was no heterogeneity if ² < 50% and and a fixed effects model was used; otherwise, we thought that heterogeneity existed in the incorporated populations and a random effects model was used instead. Subsequently, we conducted a subgroup analysis according to race. Hardy-Weinberg equilibrium (HWE) was evaluated by ² test in control groups with indicating a deviation from HWE. Sensitivity analysis was utilized to estimate the robustness and stability of the meta-analysis results by deleting all the studies one by one. Next, Begg’s funnel plot and Egger’s test were used to evaluate publication bias. For each SNP, five genetic models were evaluated to assess the correlation with cervical cancer susceptibility: the allele model, dominant model, recessive model, heterozygote model, and homozygous model. The statistical analyses were performed using RevMan 5.3 and Stata 12.0 software. All values were two sided, and was considered to be statistically significant.

3. Results

3.1. Characteristics of Included Studies

By searching the electronic databases systematically, we initially retrieved 72 articles (Figure 1). After excluding duplicate studies, 32 articles remained. Further reviewing of the titles and abstracts of the identified studies allowed the removal of 22 articles. Of those removed, 16 were clearly irrelevant to TLR9 polymorphisms, 4 were review papers, and 2 were meta-analyses. We downloaded the remaining 10 articles as full-text reports and reviewed them carefully. One record was excluded for containing duplicate samples. Finally, 9 case-control studies containing 3331 cases and 4109 controls were included in the meta-analysis, among which 8 studies were about rs352140 (G2848A) and 4 articles were about rs187084 (−1486T/C). All studies were based on Caucasian or Chinese Han populations.

Figure 1: Flow diagram of searching procedure.

These 9 case-control studies were published between 2011 and 2017. Four studies were performed with Caucasians, while five papers were based on Chinese Han populations. All the studies were performed using PCR-RFLP to identify the polymorphism sites except two studies, which used the TaqMan and Illumina GoldenGate methods. Two of the nine control groups were hospital-based groups, while the rest were population-based groups (Table 1).

Table 1: Characteristics of the included studies.

First author
Year
Race
Number (case/control)
Study design
Source of controls
Genotyping method
Study quality (NOS)

Bi [11]
2014
Chinese Han
102/100
CC
Population
PCR-RFLP
7

Bodelon [12]
2014
Caucasian
876/1100
CC
Population
Illumina GoldenGate
8

Chen [13]
2012
Chinese Han
712/717
CC
Population
PCR-RFLP
7

Jin [14]
2017
Chinese Han
420/842
CC
Hospital
PCR-RFLP
7

Lai [15]
2013
Chinese Han
120/100
CC
Hospital
PCR-RFLP
8

Pandey [16]
2011
Caucasian
200/200
CC
Population
PCR-RFLP
7

Roszak [17]
2012
Caucasian
426/460
CC
Population
PCR-RFLP
7

Xu [18]
2017
Chinese Han
345/330
CC
Population
TaqMan
7

Zidi [19]
2016
Caucasian
130/260
CC
Population
PCR-RFLP
8

CC: case-control; PCR: polymerase chain reaction; RFLP: restriction fragment length polymorphism.

We applied the Newcastle-Ottawa scale (NOS) to estimate the quality of the nine included studies [20]. The results showed that all studies were of high quality (Table 1). The distributions of the genotypes and allele frequencies of rs352140 (G2848A) and rs187084 (−1486T/C) are shown in Table 2. The distributions of the genotypes in the nine control groups were in accordance with Hardy-Weinberg equilibrium (HWE) except for two studies [15, 16] (Table 2).

Table 2: TLR9 polymorphism genotype distribution and allele frequency in cases and controls.

First author
Genotype (N)
Allele frequency (N)
HWE (p value)

Case
Control
Case
Control

2848G>A
Total
GG
GA
AA
Total
GG
GA
AA
G
A
G
A

Bi
102
33
58
11
100
31
47
22
124
80
109
91
0.601

Bodelon
876
NA
NA
NA
1100
NA
NA
NA
787
965
1036
1164
0.81

Jin
420
208
160
52
842
543
257
42
576
264
1343
341
0.111

Lai
120
98
14
8
100
97
2
1
210
30
196
4
<0.005

Pandey
200
59
115
26
200
59
112
29
233
167
230
170
0.039

Roszak
426
87
230
109
460
122
235
103
404
448
479
441
0.614

XU
345
135
163
47
330
131
152
47
433
257
414
246
0.786

Zidi
130
42
48
40
260
83
117
60
132
128
283
237
0.134

−1486T>C
Total
TT
TC
CC
Total
TT
TC
CC
T
C
T
C
HWE

Bi
102
25
68
9
100
26
54
20
118
86
106
94
0.401

Chen
694
246
346
102
715
289
319
107
838
550
897
533
0.220

Lai
120
118
1
1
100
99
0
1
237
3
198
2
<0.005

Roszak
426
141
206
79
460
193
203
64
488
364
589
331
0.367

HWE: Hardy-Weinberg equilibrium.

3.2. Meta-Analysis Results

There were 8 studies on rs352140 that included 2619 cases and 3392 healthy controls. The meta-analysis results did not show a statistical relationship between rs352140 and the risk of cervical cancer in any of the five genetic models: allele model (OR = 1.20, ) (Figure 2), dominant model (OR = 1.30, ), recessive model (OR = 1.23, ), heterozygote model (OR = 1.24, ), or homozygous genetic model (OR = 1.34, ) (Table 3). Considering the heterogeneity among studies, a subgroup analysis was performed and stratified by race. As presented in Table 3, the results were still stable (), but the A allele was correlated to increasing cervical cancer susceptibility in Caucasians based on the allele genetic model (OR = 1.11, ) (Figure 3). After excluding the two studies that were not in accordance with HWE, the pooled results did not change (data not shown).

Figure 2: Forest plots of the association between TLR9 rs352140 polymorphism and cervical cancer risk in the allele genetic model.

Table 3: Meta-analysis results of rs352140 based on five genetic models.

Genetic models
OR (95% CI)
p value
Heterogeneity
Effects model

I ² (%)
p value

Allele model (A versus G)

Overall
1.20 (0.97, 1.47)
0.09
83
0.00001
R

Race

    Chinese Han
1.46 (0.86, 2.48)
0.16
91
<0.0001
R

    Caucasian
1.11 (1.01, 1.22)
0.03
0
0.62
R

Dominant model (AA + GA versus GG)

Overall
1.30 (0.97, 1.74)
0.08
73
0.001
R

Race

    Chinese Han
1.54 (0.91, 2.61)
0.11
83
0.0006
R

    Caucasian
1.17 (0.91, 1.49)
0.22
17
0.30
R

Recessive model (AA versus GA + GG)

Overall
1.23 (0.81, 1.86)
0.34
77
0.0002
R

Race

    Chinese Han
1.35 (0.54, 3.39)
0.52
87
<0.0001
R

    Caucasian
1.20 (0.94, 1.51)
0.14
0
0.39
R

Heterozygous genetic model (GA versus GG)

Overall
1.24 (0.96, 1.59)
0.10
58
0.03
R

Race

    Chinese Han
1.42 (0.94, 2.15)
0.09
67
0.03
R

    Caucasian
1.09 (0.81, 1.49)
0.56
38
0.20
R

Homozygous genetic model (AA versus GG)

Overall
1.34 (0.83, 2.15)
0.23
78
0.0001
R

Race

    Chinese Han
1.53 (0.57, 4.15)
0.40
88
<0.0001
R

    Caucasian
1.30 (0.98, 1.73)
0.07
0
0.42
R

F: fixed-effect model; R: random-effect model; OR: odds ratio; 95% CI: 95% confidence interval.

Figure 3: Subgroup analysis of the association between TLR9 rs352140 polymorphism and the risk of cervical cancer stratified by race in the allele genetic model.

With regard to rs187084, there were 4 records involving 1342 cases and 1375 controls. The pooled ORs suggested a significant association between the SNP and cervical cancer risk: allele model (OR = 1.15, ), dominant model (OR = 1.30, ) (Figure 4), and heterozygote model (OR = 1.32, ); nevertheless, it was not related to cancer susceptibility on the recessive model or homozygous genetic model () (Table 4). We then performed subgroup analysis based on ethnicity to reassess the relationship between rs187084 and cervical cancer risk. The results were stable. Data revealed that TLR9 rs187084 was related to cervical cancer risk of Chinese Han population in the dominant model (OR = 1.22, ) (Figure 5) and heterozygous genetic model (OR = 1.28, ). There was no statistical correlation in the recessive model (OR = 0.70, ) or homozygous genetic model (OR = 0.88, ) (Figure 6, Table 4). Furthermore, the same results were present after eliminating a study that deviated from HWE.

Figure 4: Forest plots of the association between TLR9 rs187084 polymorphism and cervical cancer risk in the dominant genetic model.

Table 4: Meta-analysis results of rs187084 based on five genetic models.

Genetic models
OR (95% CI)
p value
Heterogeneity
Effects model

I ² (%)
p value

Allele model (C versus T)

Overall
1.15 (1.03, 1.29)
0.02
43
0.15
F

Race

    Chinese Han
1.06 (0.92, 1.23)
0.39
0
0.38
F

    Caucasian
1.33 (1.10, 1.61)
0.004
—
—
—

Dominant model (CC + CT versus TT)

Overall
1.30 (1.11, 1.53)
0.002
0
0.73
F

Race

    Chinese Han
1.22 (1.00, 1.50)
0.05
0
0.90
F

    Caucasian
1.46 (1.11, 1.92)
0.007
—
—
—

Recessive model (CC versus CT + TT)

Overall
0.94 (0.59, 1.50)
0.79
63
0.04
R

Race

    Chinese Han
0.70 (0.34, 1.44)
0.34
52
0.12
R

    Caucasian
1.41 (0.98, 2.02)
0.06
—
—
—

Heterozygous genetic model (CT versus TT)

Overall
1.32 (1.11, 1.57)
0.002
0
0.95
F

Race

    Chinese Han
1.28 (1.03, 1.59)
0.02
0
0.92
F

    Caucasian
1.39 (1.04, 1.86)
0.03
—
—
—

Homozygous genetic model (CC versus TT)

Overall
1.13 (0.71, 1.80)
0.60
56
0.08
R

Race

    Chinese Han
0.88 (0.49, 1.59)
0.67
31
0.24
R

    Caucasian
1.69 (1.14, 2.51)
0.009
—
—
—

F: fixed-effect model; R: random-effect model; OR: odds ratio; 95% CI: 95% confidence interval.

Figure 5: Subgroup analysis of the association between TLR9 rs187084 polymorphism and cervical cancer risk stratified by race in the dominant genetic model.

Figure 6: Subgroup analysis of the association between TLR9 rs187084 polymorphism and cervical cancer risk stratified by race in the homozygous genetic model.

With regard to rs187084, there were 4 records involving 1342 cases and 1375 controls. The pooled ORs suggested a significant association between the SNP and cervical cancer risk: allele model (OR = 1.15, ), dominant model (OR = 1.30, ) (Figure 4), and heterozygote model (OR = 1.32, ). However, the SNP was not related to cancer susceptibility in the recessive model or homozygous genetic model () (Table 4). We performed a subgroup analysis based on race to reassess the relationship between rs187084 and cervical cancer risk; the results were stable. Data revealed that TLR9 rs187084 was related to cervical cancer risk in the Chinese Han population in the dominant model (OR = 1.22, ) (Figure 5) and heterozygous genetic model (OR = 1.28, ). There was no statistical correlation in the recessive model (OR = 0.70, ) or homozygous genetic model (OR = 0.88, ) (Figure 6, Table 4). Furthermore, the same results were present after eliminating a study that deviated from HWE.

3.3. Detection for Heterogeneity

As presented in Table 3, there was great heterogeneity among studies relating to rs352140 in all genetic models (² > 50%, ). In consideration of this, we employed a random effects model for the meta-analysis. Additionally, the subgroup analysis was stratified by race to eliminate heterogeneity (Figure 3). It was clearly decreased in Caucasians as shown: allele model (² = 0%, ), dominant model (² = 17%, ), recessive model (² = 0%, ), heterozygote model (² = 38%, ), and homozygous genetic model (² = 0%, ). This implied that race might be a confounding factor and source of heterogeneity, whereas the pooled ORs were robust substantially.

As for heterogeneity among studies on rs187084, there was no significant heterogeneity in the allele model (² = 43%, ), dominant model (² = 0%, ), or heterozygous genetic model (² = 0%, ). However, heterogeneity was present in the recessive model (² = 63%, ) and homozygous genetic model (² = 56%, ). Similarly, a random effects model and subgroup analysis were conducted to evaluate heterogeneity. We found that the heterogeneity decreased in the Chinese Han population based on the homozygous genetic model (² = 31%, ) (Figure 6), but no significant differences were evident in the pooled results (Table 4).

3.4. Sensitivity Analysis

As mentioned before, two studies focused on rs352140 and one focused on rs187084 were not consistent with the balance of HWE in the control groups (). However, the results did not change substantially after removing them. Furthermore, a sensitivity analysis was utilized to evaluate the stability of the meta-analysis by deleting all the studies one by one. However, the overall ORs did not change significantly in any of the genetic models, indicating that the meta-analysis was robust and stable. The sensitivity analysis of the associations between the two SNPs and cervical cancer susceptibility based on the dominant genetic model is shown in Figures 7(a) and 7(b).

(a)

(b)

Figure 7: Sensitivity analysis of the association between TLR9 SNPs and risk of cervical cancer in the dominant genetic model. (a) rs352140. (b) rs187084.

3.5. Publication Bias

To detect publication bias, Begg’s test and Egger’s test were conducted. The results showed that there was no significant evidence of publication bias for either rs352140 or rs187084 based on the five genetic models () (Figures 8(a) and 8(b)).

(a)

(b)

Figure 8: Publication bias of TLR9 polymorphisms in the allele model. (a) rs352140 (Begg’s test: , Egger’s test: ). (b) rs187084 (Begg’s test: , Egger’s test: ).

4. Discussion

Persistent infection with HR-HPV is a causative factor in the progression of cervical cancer [21–23]. Cigarette smoking, age, and sexual activity may influence susceptibility to cervical cancer [3, 24].

TLR9 recognizes unmethylated CpG DNA and promotes the production of inflammatory cytokines that contribute to tumour immunity [25]. In addition, elevated TLR9 levels with persistent HPV infection can drive inflammation, thereby contributing to cervical cancer risk [26]. In recent years, a number of studies have focused on TLR9 gene polymorphisms and cervical cancer susceptibility, but the conclusions are questionable. For rs352140, several studies revealed that rs352140 was not related to cervical cancer; however, the study of Lai et al. supported that rs352140 was associated with an increased risk of cervical cancer, and the study of Pandy et al. showed that A allele had borderline significance conferring a marginally increased risk for advanced cervical cancer [11, 15, 16, 18]. With regard to rs187084, some studies indicated that rs352140 contributed to cervical cancer susceptibility, while the study of Bi et al. showed that the SNP was not associated with cervical carcinogenesis [11, 13]. The findings were even controversial when focused on the same ethnic group. Xu et al. and Bi et al. agreed that there was no relationship between rs352140 and the risk of cervical cancer, while Lai et al. and Jin et al. held the opposite point of view [11, 14, 15, 18]. Similarly, the studies of Lai et al. and Bi et al. indicated that rs187084 was not related to cervical cancer, while the study of Chen et al. supported a different conclusion [11, 13, 15].

Our results indicated that rs352140 is not associated with cervical cancer risk, whereas rs187084 increases susceptibility to cervical cancer. As shown in Table 4, rs187084 is associated with cervical cancer susceptibility in the allele model, dominant genetic model, and heterozygous genetic model but not in the recessive genetic model or homozygous genetic model. We postulated that the minor C allele at the −1486 locus might be a dominant mutation and that C allele carriers had an increased risk for cervical cancer. This is consistent with our conclusion that rs187084 is a risk factor for cervical cancer. However, too many reports of associations between genetic variants and common cancer sites are false positives [27]. Therefore, assessing the false positive report probability (FPRP) and using it to decide whether our finding is deserving of attention are important. Considering that our meta-analysis is a large sample and multicentre study, we applied 0.2 as the cut-off value for FPRP. According to Wacholder’s method, we evaluated the FPRPs for the dominant model (CC + CT versus TT) and heterozygous genetic model (CT versus TT). The FPRPs were 0.142 and 0.154, respectively [27]. Obviously, both FPRP results were below the cut-off value, suggesting that rs187084 is related to cervical cancer susceptibility and that our conclusion is noteworthy.

Subgroup analysis of rs352140 showed a significant decline in heterogeneity in the Caucasian population, implying that the heterogeneity may result from racial differences. Moreover, it was statistically significant that the A allele elevated the cervical cancer risk in the allele genetic model (OR = 1.11, ) while the pooled results in other genetic models remained steady. The evaluated FPRP was 0.751, which was greater than 0.2, implying a false positive association between rs352140 and cervical cancer risk among the Caucasian population [27]. Subgroup analysis and sensitivity analysis of rs187084 were also conducted. The pooled results remained identical, presenting an obvious association between rs187084 and increased cervical cancer risk and suggesting that our meta-analysis results were reliable and stable.

We searched for TLR9 SNPs in the NCBI database and found that rs352140 was a synonymous codon. This, to some extent, explained the reason for the lack of relevance regarding cervical cancer risk. It has been identified that regulatory polymorphisms controlling gene expression can be localized hundreds of kb away from the genes they influence [28, 29]. In this way, rs187084 is potentially a functional variant located in the TLR9 promoter region; so, we propose that this SNP contributes to the regulation of TLR9 at a basal transcript level [30, 31]. Moreover, rs187084 is located near the region that interacts with the HPV16 E6 and E7 oncoproteins [32]. Therefore, it may regulate transcription during HR-HPV infection and increase the risk of cervical cancer. TLR9 plays a fundamental role in pathogen recognition and the activation of innate immunity and was found to be overexpressed in cervical cancer [33, 34]. IFN-γ and TNF-α are two important cytokines that participate in the innate immune response. In addition to this, they stimulate antigen-presenting cells to express MHC II, enhancing antigen presentation and improving the adaptive immune response. The minor C allele of rs187084 led to significantly lower expression of protective cytokines IFN-γ and TNF-α than the T allele [35]. In addition, rs187084 C allele increased the transcriptional activity of TLR9 and provoked a higher level of gene expression [35, 36].

A previous meta-analysis on the relationship between TLR9 polymorphisms rs352140 and rs187084 and cervical cancer risk conducted by Mu et al. is consistent with our results [37]. There were five studies of rs352140 and four case-control studies of rs187084 in the study of Mu et al. that assessed the correlation between TLR9 SNPs and cervical cancer risk. Compared to the study of Mu et al., our meta-analysis included three additional qualified studies [14, 18, 19] to assess the association and therefore obtained more reliable and persuasive conclusions. Moreover, in order to reduce the heterogeneity caused by racial differences, we performed a subgroup analysis stratified by race. Thus, our findings are more convincing and rigorous.

Several limitations are present in our meta-analysis. First, even though we retrieved records from electronic databases systematically, there were likely still some unpublished studies we could not acquire because of their negative results, leading to potential bias. Second, the incidence of cervical cancer is highest in sub-Saharan Africa, Latin America, the Caribbean, and Melanesia, where people of African origin account for the majority [1]. However, there were no statistics and studies of interest focused on men of African descent even though several databases were searched. This led to our study focusing on Caucasian and Chinese Han populations, resulting in a bias in the relationship that we are concerned. Additionally, although we considered the effect of HPV on our conclusions and attempted to perform a subgroup analysis, few of the included studies contained HPV data. Therefore, we were unable to conduct a subgroup analysis on HPV infection. Future studies containing more comprehensive information are needed to obtain more reliable conclusions.

In conclusion, this meta-analysis demonstrated that there is no association between TLR9 rs352140 (G2848A) and cervical cancer susceptibility, while rs187084 (−1486T/C) has a statistical correlation with the risk of cervical cancer. The rs187084 C allele is considered a deleterious allele and increases cervical cancer risk. Further large, well-designed case-control studies are needed to authenticate these results.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there is no conflict of interest regarding the publication of this paper.

Acknowledgments

The research is supported by the National Natural Science Foundation of China (no. 81472428) and New Century Excellent Talents in University (NCET-12-0441).

References

L. A. Torre, F. Bray, R. L. Siegel, J. Ferlay, J. Lortet-Tieulent, and A. Jemal, “Global cancer statistics, 2012,” CA: A Cancer Journal for Clinicians, vol. 65, no. 2, pp. 87–108, 2015.

W. Chen, R. Zheng, P. D. Baade et al., “Cancer statistics in China, 2015,” CA: A Cancer Journal for Clinicians, vol. 66, no. 2, pp. 115–132, 2016.

Q. Zhang, W. Xie, F. Wang et al., “Epidemiological investigation and risk factors for cervical lesions: cervical cancer screening among women in rural areas of Henan province China,” Medical Science Monitor, vol. 22, pp. 1858–1865, 2016.

N. Muñoz, F. X. Bosch, S. de Sanjosé et al., “Epidemiologic classification of human papillomavirus types associated with cervical cancer,” The New England Journal of Medicine, vol. 348, no. 6, pp. 518–527, 2003.

S. de Sanjose, M. Diaz, X. Castellsague et al., “Worldwide prevalence and genotype distribution of cervical human papillomavirus DNA in women with normal cytology: a meta-analysis,” The Lancet Infectious Diseases, vol. 7, no. 7, pp. 453–459, 2007.

M. Zhao, Y. Li, X. Wei et al., “Negative immune factors might predominate local tumor immune status and promote carcinogenesis in cervical carcinoma,” Virology Journal, vol. 14, no. 1, p. 5, 2017.

A. M. Mehta, M. Mooij, I. Branković, S. Ouburg, S. A. Morré, and E. S. Jordanova, “Cervical carcinogenesis and immune response gene polymorphisms: a review,” Journal of Immunology Research, vol. 2017, 12 pages, 2017.

Y. Hao, J. L. Yuan, A. Abudula, A. Hasimu, N. Kadeer, and X. Guo, “TLR9 expression in uterine cervical lesions of Uyghur women correlate with cervical cancer progression and selective silencing of human papillomavirus 16 E6 and E7 oncoproteins in vitro,” Asian Pacific Journal of Cancer Prevention, vol. 15, no. 14, pp. 5867–5872, 2014.

U. A. Hasan, E. Bates, F. Takeshita et al., “TLR9 expression and function is abolished by the cervical cancer-associated human papillomavirus type 16,” Journal of Immunology, vol. 178, no. 5, pp. 3186–3197, 2007.

G.-X. Wan, Y. W. Cao, W. Q. Li, Y. C. Li, W. J. Zhang, and F. Li, “Associations between TLR9 polymorphisms and cancer risk: evidence from an updated meta-analysis of 25,685 subjects,” Asian Pacific Journal of Cancer Prevention, vol. 15, no. 19, pp. 8279–8285, 2014.

X. Bi, X. H. Xu, and X. Yu, “Study on TLRs polymorphisms and cervical cancer susceptibility,” Progress in Obstetrics & Gynecology, vol. 23, no. 7, pp. 520–523, 2014.

C. Bodelon, M. M. Madeleine, L. G. Johnson et al., “Genetic variation in the TLR and NF-κB pathways and cervical and vulvar cancer risk: a population-based case-control study,” International Journal of Cancer, vol. 134, no. 2, pp. 437–444, 2014.

X. Chen, S. Wang, L. Liu et al., “A genetic variant in the promoter region of Toll-like receptor 9 and cervical cancer susceptibility,” DNA and Cell Biology, vol. 31, no. 5, pp. 766–771, 2012.

Y. Jin, S. Qiu, N. Shao, and J. Zheng, “Association of toll-like receptor gene polymorphisms and its interaction with HPV infection in determining the susceptibility of cervical cancer in Chinese Han population,” Mammalian Genome, vol. 28, no. 5-6, pp. 213–219, 2017.

Z.-Z. Lai, Ni-Zhang, X. L. Pan, and L. Song, “Toll-like receptor 9 (TLR9) gene polymorphisms associated with increased susceptibility of human papillomavirus-16 infection in patients with cervical cancer,” The Journal of International Medical Research, vol. 41, no. 4, pp. 1027–1036, 2013.

S. Pandey, B. Mittal, M. Srivastava et al., “Evaluation of Toll-like receptors 3 (c.1377C/T) and 9 (G2848A) gene polymorphisms in cervical cancer susceptibility,” Molecular Biology Reports, vol. 38, no. 7, pp. 4715–4721, 2011.

A. Roszak, M. Lianeri, A. Sowińska, and P. P. Jagodziński, “Involvement of toll-like receptor 9 polymorphism in cervical cancer development,” Molecular Biology Reports, vol. 39, no. 8, pp. 8425–8430, 2012.

X. Xu, C. Yin, Z. Yan, Y. Hu, and Z. Yang, “The correlation of TLR9 gene polymorphism with occurrence and development of cervical cancer in Yunnan Han population,” Journal of Guizhou Medical University, vol. 1, pp. 21–25, 2017.

S. Zidi, I. Sghaier, E. Gazouani, A. Mezlini, and B. Yacoubi-Loueslati, “Evaluation of Toll-like receptors 2/3/4/9 gene polymorphisms in cervical cancer evolution,” Pathology Oncology Research, vol. 22, no. 2, pp. 323–330, 2016.

G. A. Wells, B. J. Shea, D. O'Connell et al., “The Newcastle–Ottawa Scale (NOS) for assessing the quality of non-randomized studies in meta-analysis,” Available at : http://www.ohri.ca/programs/clinical_epidemiology/oxford.asp.

S.-J. Lee, A. Yang, T.-C. Wu, and C.-F. Hung, “Immunotherapy for human papillomavirus-associated disease and cervical cancer: review of clinical and translational research,” Journal of Gynecologic Oncology, vol. 27, no. 5, 2016.

X. Castellsagué, “Natural history and epidemiology of HPV infection and cervical cancer,” Gynecologic Oncology, vol. 110, no. 3, pp. S4–S7, 2008.

S. S. Wang, P. Gonzalez, K. Yu et al., “Common genetic variants and risk for HPV persistence and progression to cervical cancer,” PLoS One, vol. 5, no. 1, article e8667, 2010.

S. Campino, J. Forton, S. Auburn et al., “TLR9 polymorphisms in African populations: no association with severe malaria, but evidence of cis-variants acting on gene expression,” Malaria Journal, vol. 8, no. 1, p. 44, 2009.

S. Cornélie, J. Hoebeke, A. M. Schacht et al., “Direct evidence that Toll-like receptor 9 (TLR9) functionally binds plasmid DNA by specific cytosine-phosphate-guanine motif recognition,” The Journal of Biological Chemistry, vol. 279, no. 15, pp. 15124–15129, 2004.

F. Cannella, A. Pierangeli, C. Scagnolari et al., “TLR9 is expressed in human papillomavirus-positive cervical cells and is overexpressed in persistent infections,” Immunobiology, vol. 220, no. 3, pp. 363–368, 2015.

S. Wacholder, S. Chanock, M. Garcia-Closas, L. El Ghormli, and N. Rothman, “Assessing the probability that a positive report is false: an approach for molecular epidemiology studies,” Journal of the National Cancer Institute, vol. 96, no. 6, pp. 434–442, 2004.

B. E. Stranger, A. C. Nica, M. S. Forrest et al., “Population genomics of human gene expression,” Nature Genetics, vol. 39, no. 10, pp. 1217–1224, 2007.

J. T. Forton, I. A. Udalova, S. Campino, K. A. Rockett, J. Hull, and D. P. Kwiatkowski, “Localization of a long-range cis-regulatory element of IL13 by allelic transcript ratio mapping,” Genome Research, vol. 17, no. 1, pp. 82–87, 2006.

K. Tao, M. Fujii, S. I. Tsukumo et al., “Genetic variations of Toll-like receptor 9 predispose to systemic lupus erythematosus in Japanese population,” Annals of the Rheumatic Diseases, vol. 66, no. 7, pp. 905–909, 2007.

M. W. Ng, C. S. Lau, T. M. Chan, W. H. S. Wong, and Y. L. Lau, “Polymorphisms of the toll-like receptor 9 (TLR9) gene with systemic lupus erythematosus in Chinese,” Rheumatology, vol. 44, no. 11, pp. 1456-1457, 2005.

L. Hamann, C. Glaeser, A. Hamprecht, M. Gross, A. Gomma, and R. R. Schumann, “Toll-like receptor (TLR)-9 promotor polymorphisms and atherosclerosis,” Clinica Chimica Acta, vol. 364, no. 1-2, pp. 303–307, 2006.

C. Martínez-Campos, M. Bahena-Román, K. Torres-Poveda, A. I. Burguete-García, and V. Madrid-Marina, “TLR9 gene polymorphism −1486T/C (rs187084) is associated with uterine cervical neoplasm in Mexican female population,” Journal of Cancer Research and Clinical Oncology, vol. 143, no. 12, pp. 2437–2445, 2017.

E. Fehri, E. Ennaifer, M. Ardhaoui et al., “Expression of Toll-like receptor 9 increases with progression of cervical neoplasia in Tunisian women - a comparative analysis of condyloma, cervical intraepithelial neoplasia and invasive carcinoma,” Asian Pacific Journal of Cancer Prevention, vol. 15, no. 15, pp. 6145–6150, 2014.

D. Bharti, A. Kumar, R. S. Mahla et al., “The role of TLR9 polymorphism in susceptibility to pulmonary tuberculosis,” Immunogenetics, vol. 66, no. 12, pp. 675–681, 2014.

J. Fischer, A. N. R. Weber, S. Böhm et al., “Sex-specific effects of TLR9 promoter variants on spontaneous clearance of HCV infection,” Gut, vol. 66, no. 10, pp. 1829–1837, 2017.

X. Mu, J. Zhao, X. Yuan et al., “Gene polymorphisms of Toll-like receptor 9 −1486T/C and 2848G/A in cervical cancer risk,” International Journal of Gynecological Cancer, vol. 25, no. 7, pp. 1173–1178, 2015.