A Meta-Analysis for Association of XRCC3 rs861539, MTHFR rs1801133, IL-6 rs1800795, IL-12B rs3212227, TNF-α rs1800629, and TLR9 rs352140 Polymorphisms with Susceptibility to Cervical Carcinoma.

Cancer is one of the main public health problems globally with about 18.1 million new cancer cases and 9.6 million cancer deaths in 2018 (Ghelmani et al., 2021; Jarahzadeh et al., 2021; Antikchi et al., 2021). The World Health Organization (WHO) states cervical cancer is the second most frequent cancer among women worldwide (Hamadani et al., 2017; Da Silva et al., 2021) with an estimated 570,000 new cases and 311,000 deaths in 2018 worldwide (Yi et al., 2020). It is reported that 85% of cervical cancer cases and 87% of the cervical cancer deaths occur in less developed countries (WHO/ICO, 2010). With over 500,000 cases of cervical cancer reported each year, nearly 80 percent of those are in developing countries, including Africa with 68,000, an estimated 77,000 in Latin America and the Caribbean, and 245,000 in Asia (Arbyn et al., 2020). Cervical cancer is driven by persistent infection with one of 15 carcinogenic human papillomavirus (HPV) types (Ghaemmaghamiet al., 2008; Chen et al., 2011). However, the impact of HPV type and intratypic variants on patient outcomes is far less understood (Karimi-Zarchi et al., 2013; Baghestani et al., 2018; Rader et al., 2019). To date, the occurrence and development of cervical cancer is suggested to be associated with persistent HPV infection (Ghaemmaghami et al., 2008; Zarchi et al., 2010; Pal and Kundu, 2020). The DNA of HPV integrates into the host cell genome and disrupts the open reading frame and causes overexpression of E6 and E7 genes (Binesh et al., 2012; Soheili et al., 2016). However, the specific molecular mechanisms and potential single gene require further studies. Familial based studies and evaluation of inherited genetic variations revealed that host genetic factors have a role in cervical cancer pathogenesis (Yang et al., 2020; Ahmadi et al., 2021). Recently, a comprehensive review study showed that CDK1, CCNB1, ITGB1, FN1, MMP9 and STAT1 played different roles in the progression of cervical cancer through different signaling pathways (Sayad, Ahmadi, Nekouian, et al., 2020; Yi et al., 2020). Moreover, studies have shown that cervical cancer has a heritable genetic component. However, as other complex disease, the identified loci only explain a minority of the risk of cervical cancer. This gap of the heritability explained by genetic markers genome-wide association studies and the heritability identified in familial studies has been termed as ‘missing heritability’ (Chen et al., 2016; Leo et al., 2017; Sayad, Ahmadi, Moradi, et al., 2020). Thus, our understanding of the genetic basis of cervical cancer is still limited. In this meta-analysis, we explore the role of XRCC3 rs861539, MTHFR rs1801133, IL-6 rs1800795, IL-12B rs3212227, TNF-α rs1800629, and TLR9 rs352140 polymorphisms in susceptibility to cervical cancer.

Seyedeh Fatemeh Parsaeian, Fatemeh Asadian, Mojgan Karimi-Zarchi: conceptualization, investigation. Seyed Alireza Dastgheib, Sepideh Setayesh: Software, original draft preparation. Atiyeh Javaheri, Razieh Sadat Tabatabaie: Investigation. Fatemeh Asadian: Investigation, writing. Fatemeh Asadian, Hossein Neamatzadeh: Methodology, software. Hossein Golestanpour, Hossein Neamatzadeh: Formal analysis, investigation. Hossein Golestanpour: Project administration. Atiyeh Javaheri, Razieh Sadat Tabatabaie: Writing, reviewing, editing.

This article does not contain any studies with human participants or animals performed by any of the authors.

Not applicable for this manuscript.

Characteristics of the Enrolled Studies The flowchart of selection of studies and reasons for exclusion is presented in Figure 1. There were 1,013 articles relevant to our search words and 395 duplicated studies were excluded. Then, 207 studies removed after reading titles and abstracts. Following the inclusion exclusion criteria 352 studies were excluded. Finally, a total of 59 studies were included in this meta-analysis. There are seven studies with 1112 cases and 1233 controls on XRCC3 rs861539 (He et al., 2008; Xiao, 2009; Settheetham-Ishida et al., 2011; Djansugurova et al., 2013; Pérez et al., 2013; Al-Harbi et al., 2017), 14 studies with 2694 cases and 3349 controls MTHFR rs1801133 (Lambropoulos et al., 2003; Sull et al., 2004; Kang et al., 2005; Zoodsma et al., 2005; Delgado-Enciso et al., 2006; Ma et al., 2006; Nandan et al., 2008; Shekari et al., 2008; Kohaar et al., 2010; Yang et al., 2011; Prasad and Wilkhoo, 2011; Tong et al., 2011; Mostowska et al., 2011; von Keyserling et al., 2011), four studies with 1121 cases and 1109 controls on IL-12B rs3212227 (Chen et al., 2009; Do Carmo Vasconcelos De Carvalho et al., 2012; Roszak, Mostowska, et al., 2012), seven studies with 1452 cases and 2186 controls on IL-6 rs1800795 (Nogueira De Souza et al., 2006; Gangwar et al., 2009; Grimm et al., 2011; Shi et al., 2013; de Lima Júnior et al., 2016; Pu et al., 2016; Zidi et al., 2017), 20 studies with 4,781 cases and 4909 controls on TNF-α rs1800629 (Stanczuk et al., 2019; Jang et al., 2001; Calhoun et al., 2002; Gostout et al., 2003; Deshpande et al., 2005; Duarte et al., 2005; Govan et al., 2006; Kohaar et al., 2007; Singh et al., 2009; Zu et al., 2010; Ivansson et al., 2010; Zuo et al., 2011; Wang et al., 2011, 2012; Badano et al., 2012; Barbisan et al., 2012; Sousa et al., 2014; Roszak et al., 2015; Zidi et al., 2015) and seven studies with 1743 cases and 2292 controls on TLR9 rs352140 (Pandey et al., 2011; Roszak, Lianeri, et al., 2012; Lai et al., 2013; Bi, 2014; Zidi et al., 2016; Jin et al., 2017). The details of included studies were shown in Tables 1 and 2 . All those 22 studies were reported in English and Chinese. Among those 59 studies, 31 studies were from Asian populations, 18 studies from Caucasian populations, five studies from African populations, and five studies were from mixed populations. The studies have been carried out in China, Thailand, Argentina, Kazakhstan, Brazil, Saudi Arabia, Greece, Korea, Netherlands, Mexico, India, Poland, Germany, Austria, Tunisia, Zimbabwe, USA, South Africa, Sweden, and Portugal. The sample size of cases ranged from 21 to 636, while the sample size of controls ranged from 73 to 800 in the controls. Eight different methods including AS-PCR, PCR-RFLP, Direct Sequencing, SnapShot, TaqMan, LDR-PCR, ARMS-PCR, qRT-PCR and HMR were used to genotyping those polymorphisms. Hardy-Weinberg equilibrium (HWE) test was calculated for all publications and P<0.05 was considered as a departure from HWE ( Table 1 ). Quantitative Data Synthesis XRCC3 RS861539 The summary for the association of XRCC3 RS861539 polymorphism with cervical cancer risk are shown in Table 3 . Pooled data revealed that XRCC3 RS861539 polymorphism was significantly associated with susceptibility to cervical cancer under the heterozygote genetic model (TC vs. CC: OR = 1.00, 95% CI 1.066-1.585, p = 0.009). Moreover, stratified analysis by ethnicity revealed that the polymorphisms was significantly associated with cervical cancer in Asian women under three genetic models, i.e., allele (T vs. C: OR= 1.302, 95% CI 1.076-1.576, P= 0.007), heterozygote (TC vs. CC: OR = 1.441, 95% CI 1.113-1.867, p = 0.006) and dominant (TT+TC vs. CC: OR = 1.469, 95% CI 1.148-1.880, p = 0.002), but not in Caucasian women. TNF-α rs1800629 The summary for the association of TNF-α rs1800629 polymorphism with cervical cancer risk are shown in Table 3 . Pooled data from all eligible studies indicated that the TNF-α rs1800629 polymorphism was associated with cervical cancer risk in overall population under four genetic models i.e., allele (A vs. G: OR = 1.277, 95% CI = 1.104-1.477, P = 0.001), homozygote (AA vs. GG: OR = 1.333, 95% CI = 1.062-1.674, P = 0.013), heterozygote (AG vs. GG: OR = 1.307, 95% CI = 1.064-1.605, P = 0.011), and dominant (AA+AG vs. GG: OR = 1.324, 95% CI = 1.104-1.587, P = 0.002). The subgroup analysis by ethnicity also showed that this polymorphism was associated with cervical cancer in Caucasian women under three genetic models i.e., allele (A vs. G, OR = 1.242, 95% CI = 1.043-1.478, P = 0.015); homozygote (AA vs. GG, OR = 1.586, 95% CI = 1.147-2.193, P = 0.005), and recessive (AA vs. AG+GG, OR = 1.569, 95% CI = 1.137-2.165, P = 0.006) and in African women under two genetic models i.e., heterozygote (AG vs. GG, OR = 1.670, 95% CI = 1.228-2.270, P = 0.001) and dominant (AA+AG vs. GG, OR = 1.453, 95% CI = 1.111-1.902, P = 0.006), but not in Asian women. MTHFR rs1801133 The summary for the association of MTHFR rs1801133 polymorphism with cervical cancer risk are shown in Table 3 . Pooled results showed that the MTHFR rs1801133 polymorphism was not associated with cervical cancer in overall population and by ethnicity. IL-12B rs3212227 The summary for the association of IL-12B rs3212227 polymorphism with cervical cancer risk are shown in Table 4 . Pooled data showed that the IL-12B rs3212227 polymorphism was not associated with cervical cancer globally. Stratified analysis by ethnicity revealed that this polymorphism was associated with cervical cancer in Asian women under two genetic models i.e., heterozygote (CA vs. AA: OR = 1.349, 95% CI 1.032-1.762, p=0.028) and dominant (CC+CA vs. AA: OR = 1.340, 95% CI 1.041-1.725, p =0.023). IL-6 rs1800795 The summary for the association of IL-6 rs1800795 polymorphism with cervical cancer risk are shown in Table 4 . Pooled data from all eligible studies showed that the IL-6 rs1800795 polymorphism was significantly associated with cervical cancer risk under four genetic models, i.e., allele (C vs. G: OR = 1.294, 95% CI 1.071-1.564, p= 0.007), homozygote (CC vs. GG: OR = 1.633, 95% CI 1.059-2.520, p= 0.027), dominant (CC+CG vs. GG: OR = 1.312, 95% CI 1.048-1.643, p= 0.018) and recessive (CC vs. CG+GG: OR = 1.592, 95% CI 1.268-1.999, p≤0.001). Moreover, subgroup analysis by ethnicity revealed an increased risk of cervical cancer in Asian women. TLR9 rs352140 The summary for the association of TLR9 rs352140 polymorphism with cervical cancer risk are shown in Table 4 . Pooled results showed that the TLR9 rs352140 polymorphism was not associated with cervical cancer in overall population and by ethnicity. Test of heterogeneity In the current meta-analysis the χ 2 test and I 2 statistics were used for assessing the heterogeneity of the included studies. Results indicate that there was statistical heterogeneity for XRCC3 rs861539, MTHFR rs1801133, IL-6 rs1800795, IL-12B rs3212227, TNF-α rs1800629, and TLR9 rs352140 polymorphisms under most genetic models. Thus, the random effect model was used for evaluating the pooled OR and 95% CI for those models. Subgroup analyses showed that ethnicity of participants might contribute to part of heterogeneity. Sensitivity Analysis The process of performing a meta-analysis involves a sequence of decisions and it is important to perform a sensitivity analysis in order to assess the impact of different decisions on pooled data. Thus, a sensitivity analysis was carried out by excluding a single study in turn on pooled ORs. The results showed that no individual study had an influence on the pooled OR all involved polymorphisms at XRCC3 rs861539, MTHFR rs1801133, IL-6 rs1800795, IL-12B rs3212227, TNF-α rs1800629, and TLR9 rs352140 polymorphisms, suggesting the stability of our findings. Moreover, excluding HWE deviated studies suggested that there were no independent studies that significantly influenced our pooled data. Publication Bias Begg’s funnel plot and Egger’s test were used for evaluating publication bias. As shown in Tables 3 and 4 , the Egger’s test results showed that there was no publication bias for the XRCC3 rs861539, MTHFR rs1801133, IL-6 rs1800795, IL-12B rs3212227, TNF-α rs1800629, and TLR9 rs352140 polymorphisms under all five genetic models. Moreover, Begg’s funnel did not statistically revealed a significant publication bias in any of the models for all involved polymorphisms. Thus, the publication bias tests revealed that our pooled ORs were reliable. Characteristics of Studies Included in the Meta-Analysis Characteristics of the Case-Control Studies Included in the Meta-Analyses Results of the Association of XRCC3 RS861539, MTHFR rs1801133 and IL-12B rs3212227 Polymorphisms with Cervical Cancer Risk NA, not applicable Summary Risk Estimates for Association of IL-6 rs1800795, TNF-α rs1800629 and TLR9 rs352140 Polymorphisms with Cervical Cancer Risk

The authors declare that they have no conflict of interest.

Genetic factors have been shown to influence the susceptibility of patients to various diseases and have attracted increasing attention (Motamedi et al., 2012; Mazaheri et al., 2014; Kabiri Rad et al., 2018). Several efforts have been made to identify the genetic susceptibility factors underlying development of cervical cancer. However, only a few polymorphisms have shown consistency among studies. In this meta-analysis, the association of XRCC3 rs861539, MTHFR rs1801133, IL-6 rs1800795, IL-12B rs3212227, TNF-α rs1800629 and TLR9 rs352140 polymorphisms with susceptibility to cervical carcinoma was assessed by including all relevant studies. Yuan et al., (2021) in a meta-analysis based on 15 studies with 5,740 cases and 9,931 controls revealed that there was no significant association between the XRCC3 Thr241Met and the risk of gynecological malignancies. However, their subgroup analysis by ethnicity showed that XRCC3 Thr241Met was associated with an increased risk of gynecological malignancies in Asians. Moreover, their stratified analysis by cancer type indicated that XRCC3 Thr241Met associated with cervical cancer in Asians (CT vs. CC: OR=1.50, 95%CI=1.04-2.14; TT vs. CC: OR=3.14, 95%CI=1.38-7.14; CT+TT vs. CC: OR=1.64, 95% CI=1.17-2.31). Abbas et al., (2010) in a case-control study with 260 cervical cancer cases and 265 controls evaluated the association of XRCC1+399A/G, XRCC2+31467G/A and XRCC3+18067C/T polymorphisms with cervical cancer in Indian women. Their results showed that XRCC2+31479G/A and XRCC3+18067C/T polymorphisms were not associated with cervical cancer. However, they showed that he XRCC1+399A/G is linked with cervical cancer in the Indian population. Al-Harbi et al., (2017) in a study among 232 cervical cancer cases and 313 control subjects evaluated the association of CDKN1A C31A, ATM G1853A, HDM2 T309G, TGFB1 T10C, XRCC1 G399A, and XRCC3 C241T with cervical cancer among Saudi Arabian women. They showed that the TGFB1 T10C and XRCC1 G399A polymorphisms were associated with cervical cancer risk. Our pooled results showed that the MTHFR rs1801133 polymorphism was not associated with cervical cancer in overall population and by ethnicity. Silva et al., (2019)reported that there were no differences in the genotypic and allelic distribution of MTHFR C677T polymorphism between remission (with the presence of pre-neoplastic lesions) and Persistence (with the presence of pre-neoplastic lesions). The same authors, in another study showed that MTHFR C677T polymorphism was not associated with cervical cancer and HPV infection (Silva et al., 2019). However, Sohrabi et al., in a case-control study evaluated the association of MTHFR A1298C and C677T variants among in 50 cervical intraepithelial neoplasia cases, 98 HPV-positive subjects and 47 non-cancerous/non-HPV patients as healthy controls. Their results showed that MTHFR 1298 CC is more likely to be a potential risk factor for HPV-cervical cancer progression (Sohrabi et al., 2020). Zhou et al., (2020) indicated that the MTHFR rs4846048 enhanced the risk of cervical cancer through association with miR-522. Gong et al., (2018) reported that MTHFR C677T polymorphism was not associated with the risk of cervical cancer or cervical intra-epithelial neoplasia, while, the MTHFR A1298C polymorphism could increase the risk of both cervical cancer and cervical intra-epithelial neoplasia. In the current study pooled data indicated that the TNF-α rs1800629 polymorphism was associated with cervical cancer risk. Wang et al., (2012) in meta-analysis based 27 studies showed that both TNF-α -238 and -308 G/A polymorphisms could be used to identity individual with elevated susceptibility to cervical cancer in by ethnicity (Bi, 2014). Behboodi et al., (2021) in a study among 153 Iranian cervical cancer cases and 292 free cancer subjects demonstrated that TNF-α rs1800629 was associated with increased level and risk of developing cervical cancer. However, Duvlis et al., (2020) reported that TNF-a-238G/A and TNF-a-308 G/T polymorphisms were not associated with the risk of HPV associated cervical intraepithelial lesions or cervical cancer cases in Macedonian women compared to controls. Moreover, Traore et al., (2020) revealed that TNF-308 G/A or IL18-607C/A polymorphisms were not associated with HPV infection among Burkina Faso women. Pooled data from all eligible studies showed that the IL-6 rs1800795 polymorphism was significantly associated with cervical cancer risk in overall population and among Asian women. Similarly, Duan et al., (2018) in a meta-analysis based on 7 studies showed that the IL-6 rs1800795 polymorphism is associated with risk of cervical cancer in overall. Pu et al., (2016) in a study with 360 cervical cancer cases and 728 healthy subjects showed that this polymorphism is risk factor for cervical cancer development in Chinese women. de Souza et al., (2006) in a case-control study with 56 cases and 253 controls evaluated the association of IL-6 rs1800795 polymorphism with cervical cancer in a Brazilian population. In 2017, Zidi et al., (2017) in a case-control study evaluated the effects of six different genetic variants atIL-6 with risk of cervical cancer. The study revealed that the IL-6 rs1800795 polymorphism has a has a protective role in cervical cancer development in Tunisian women.de Lima Júnior et al., (2016) revealed that the IL-6 rs1800795 polymorphism was not associated with HPV infection and healthy controls in in Brazilian women. The current meta-analysis data showed that the IL-12B rs3212227 genetic variant was not associated with susceptibility to cervical cancer, while; stratified analysis showed that this polymorphism was associated with cervical cancer risk among Asian women. Zheng et al., (2017) in a meta-analysis based on 33 articles with 10,587 cancer cases and 12,040 healthy subjects assessed the genetic association of IL-12B rs3212227 with cancer risk. Their pooled data indicated that IL-12B rs3212227 polymorphism was associated with cancer risk in overall. However, subgroup analysis by cancer type showed that IL-12B rs3212227 was not associated with cervical cancer. In 2016, Chang et al., (2015) in a meta-analysis based on 5 studies with 2552 cervical cases and 2232 healthy subjects showed that this polymorphism did not associate with risk of cervical cancer. In 2012, de Carvalho et al., found that the IL-12B rs3212227 variant has a protective role in development of cervical cancer in Brazilian women (Do Carmo Vasconcelos De Carvalho et al., 2012). Similarly, Han et al., (2008) found that IL-12B rs3212227 did not associate with cervical cancer risk in Korean population. However, Roszak et al., (2012) showed that this polymorphism was associated with increased risk of cervical cancer among polish women. The current meta-analysis data showed that TLR9 rs352140 was not associated with susceptibility to cervical cancer. Nath et al., (2020) in a study showed that TLR4/9 polymorphisms are associated with increased HPV16/18 infection susceptibility and cervical squamous cell carcinoma risk among the women of Jharkhand state. In another study in India showed the TLR4 and TLR9 polymorphisms and haplotypes with hrHPV infection and cervical cancer risk (Pandey et al., 2019). Martínez-Campos et al., (2017) in a case-control study revealed that TLR9 rs187084 is a risk factor for HPV infection, squamous intraepithelial cervical lesion, and uterine cervical neoplasm in Mexican female population. Jin et al., (2017) also evaluated the association of some variants at Toll-like receptors gene with cervical cancer among 420 Chinese cervical cancer patients and 842 controls. Their results showed that mutant alleles of TLR2 rs3775290, TLR4 rs7873784, and TLR9 rs352140 were associated with increased cervical cancer risk. Similarly, Yang et al., (2020) in meta-analysis based on eleven studies indicated that the TLR9 rs187084 and rs352140 polymorphisms may contribute to development of cervical cancer, but not TLR2-196 to -174 del/ins polymorphism. In summary, our pooled data indicated that the XRCC3 RS861539, TNF-α rs1800629, and IL-6 rs1800795 genetic variants were associated with susceptibility to cervical cancer globally. However, the MTHFR rs1801133, IL-12B rs3212227 and TLR9 rs352140 variants were not associated. Larger and more rigorous studies among different ethnicities are needed to further evaluate these associations with cervical cancer.

The datasets generated during and/or analyzed during this study are available from the corresponding author on reasonable request.

Literate Search Strategy We performed a comprehensive literature search on PubMed, Scopus, China National Knowledge Infrastructure, Wanfang, VIP Information Chinese Journal Service Platform, and China Biology Medicine disc, PubMed, EMBASE, Web of Science and the Cochrane Library, Google Scholar, Cochrane Library, EMBASE, Scientific Information Database (SID), WanFang, VIP, Chinese Biomedical Database (CBD), Scientific Electronic Library Online (SciELO), China National Knowledge Infrastructure (CNKI), IranDoc and Egyptian Knowledge Bank (EKB) Journals database to identify all relevant studies on the association of XRCC3 rs861539, MTHFR rs1801133, IL-6 rs1800795, IL-12B rs3212227, TNF-α rs1800629 and TLR9 rs352140 with cervical cancer up to 30th June 2021. The following terms were used in various combinations in this search: (‘’Uterine Cervical Neoplasm‘’ OR ‘’Cervix Cancer’’ OR ‘’Cervical Cancer’’ OR ‘’Cervical Neoplasm’’ OR ‘’Cervical Carcinoma’’ OR ‘’Squamous Cell Carcinoma’’ OR ‘’Adenocarcinoma’’) AND (‘’X-Ray Repair Cross Complementing 3’’ OR ‘’XRCC3’’ OR ‘’rs861539’’ OR ‘’g.4880C>G’’ OR ‘’c.-237C>G’’ OR ‘’n.54-321G>C’’) AND (‘’Methylene tetrahydrofolate reductase’’ OR ‘’MTHFR’’ OR ‘’677C>T’’ OR ‘’rs1801133’’ OR ‘’g.14783C>T’’ OR ‘’c.788C>T’’ OR ‘’p.Ala263Val’’ OR ‘’p.Ala222Val’’) AND (‘’Interleukin 12’’ OR ‘’IL-12’’ OR ‘’1188A>C’’ OR ‘’rs3212227’’ OR ‘’ g.158742950T>G’’ OR ‘’ g.19532A>C’’ OR ‘’ c.*159A>C’’) AND (‘Interleukin 6’’ OR ‘’IL-6’’ OR ‘’-174G>C’’ OR ‘’rs1800795’’ OR ‘’ g.4880C>G’’ OR ‘’c.-274C>G’’ OR ‘’ n.54-321G>C’’) AND (‘’Toll like receptors’’ OR ‘’TLR’’ ‘’ OR ‘’rs352140’’ OR ‘’2848G/A’’) AND (‘’Gene’’ OR ‘’Single-Nucleotide Polymorphism” or “SNP” OR ‘’Polymorphism’’ OR ‘’Genotype’ OR ‘’Allele’’ OR ‘’Variant’’ OR ‘’Variation’’ OR ‘’Mutation’’ OR ‘’Mutant’’). The search was carried out in English, Chinese and Persian. When overlapping data on the same cases were included in more than one publication, only the one with the larger sample size was selected. Moreover, the reference list of the retrieved studies and reviews manually checked to identify more potential eligible studies. Eligibility criteria To obtain the papers, the following criteria should be met to include the papers in our study: 1) studies with case-control or cohort design; 2) studies reported original data; 3) studies appraised association of XRCC3 rs861539, MTHFR rs1801133, IL-6 rs1800795, IL-12B rs3212227, TNF-α rs1800629 and TLR9 rs352140 polymorphisms with cervical cancer; and 4) studies with available and sufficient data for calculating an odds ratio (OR) with 95% confidence interval (CI). The following exclusion criteria were also used: 1) studies on other polymorphism at XRCC3, MTHFR, IL-12B, IL-6, and TNF-α genes; 2) animal studies or in vitro studies; 3) studies with in sufficient data on genotype frequencies or which the number of genotypes and alleles could not be ascertained; 4) linkage studies; 5) family based studies including sibling, twins and trios-parents studies; 6) abstracts, case reports, commentaries, editorials, conference articles, reviews, proceedings and meta-analyses; and 7) duplicates studies or overlapping data. Data Extraction Two authors reviewed and extracted necessary information independently in accordance with our inclusion criteria. For conflicting data, the authors carried out discussions until a consensus was reached. If they could not reach a consensus, disagreement was resolved by the third author who participated in the discussion. For each eligible study the following data was collected: first author name, year of publication, country of origin, ethnicity (Caucasian, Asian, African, Mixed populations), genotyping methods, sample size, allele and genotype frequency of XRCC3 rs861539, MTHFR rs1801133, IL-12B rs3212227, IL-6 rs1800795, TNF-α rs1800629 and TLR9 rs352140 polymorphisms in cervical cancer cases and controls, Minor Allele Frequency (MAFs) and Hardy-Weinberg equilibrium (HWE) in healthy controls. In this meta-analysis different case-control groups or cohorts in one publication were considered as independent studies. We did not define any minimum sample size to include in this meta-analysis. The ‘‘mixed’’ group means mixed or unknown populations. If more than one study was published by the same author(s) using repeated or overlapped data, the studies with the largest sample size or the most recently published study was included to the meta-analysis. If selected articles did not reported necessary data the corresponding authors was contacted by email to request the missing data. Statistical Analysis The comparison between the XRCC3 rs861539, MTHFR rs1801133, IL-6 rs1800795, IL-12B rs3212227, TNF-α rs1800629 and TLR9 rs352140 polymorphism and cervical cancer was assessed using pooled odds ratio (OR) and 95% confidence interval (CI). The pooled data were calculated under five genetic models, i.e., allele (B vs. A), homozygote (BB vs. AA), heterozygote (BA vs. BB), dominant (BB+BA vs. AA) and recessive (BB vs. BA+AA), in which a ‘‘A’’ denotes a major allele; ‘‘B’’ denotes a minor allele. The χ 2 test and I 2 statistics were used to assess whether there was between-study heterogeneity, in which P >.10 and I 2  < 50% could be considered that there was no statistical heterogeneity between the research results. The fixed effect model was selected for data consolidation; P ≤.10 and I 2 ≥50% could be considered that there was statistical heterogeneity between the research results, and a random effect model was used for data consolidation. If there was significant heterogeneity among studies, the random effects model (DerSimonian and Laird) was used; otherwise, the fixed-effects model (Mantel and Haenszel) was acceptable. The Hardy-Weinberg equilibrium of the control group was evaluated using the χ 2 test, and the expected and actual genotype frequencies of the control group were compared. In this meta-analysis, P-values of ≤.05 were considered statistically significant (Azadi-Yazdi et al., 2017; Mirjalili et al., 2019; Dastgheib et al., 2020). In addition, the sensitivity analysis was performed by excluding HWE-violating studies. Potential publication bias was evaluated using the Egger’s linear regression test and Begg’s quantitative test. The asymmetric plot of Egger’s test and the P-value of Begg’s test less than 0.05 were considered a significant publication bias. If there was evidence of publication bias (P<0.05), trim and fill method was applied to adjust for the effect of publication bias (Jafari et al., 2020). All statistical analyses were performed using Comprehensive Meta-Analysis (CMA) Software version 2.0 (Biostat, Englewood, USA). All tests were two-sided, and the P< 0.05 was considered statistically significant.