Author Correction: Association of TLR4 and TLR9 polymorphisms and haplotypes with cervical cancer susceptibility

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

www.nature.com/scientificreports www.nature.com/scientificreports/ TLRs are a part of innate immune system and significantly contribute in battling bacteria, viruses and other pathogens, and provide anti-tumor immunity 11 . TLRs serve as the initiator of inflammatory response generated by various factors including infection and tissue injury. Briefly, TLRs after binding to exogenous microbial or endogenous-tissue injury generated ligands activate transcription factors via adaptor protein myeloid differentiation factor 88 (MyD88) or MyD88 adaptor-like/Toll-interleukin 1 receptor domain-containing adaptor protein (Mal/TIRAP) leading to cytokines production and activation of adaptive immune response 12 .
To date, ten functional TLRs designated as TLR1 to TLR10 are expressed in humans by immune and certain non-immune cells. Of these TLRs, TLR1, 2, 4, 5, 6 and 10 are found on the cell surfaces whereas TLR3, 7, 8 and 9 are located in the endosomes or endoplasmic reticulum 13 . TLRs have also been implicated in the initiation, progression and metastasis of tumors 14,15 . Aberrant expression of different TLRs including TLR4 and TLR9 have been detected in gastric 16,17 , ovarian 18 , colorectal 19 , lung 20 , breast 21 , prostate 22 as well as cervical cancers 23 . Furthermore, Hasan et al. 24 , reported the involvement of HPV16 E6 and E7 oncoproteins in the inhibition of TLR9 transcription, leading to decreased immune response and escape for HPV16.
Moreover, as inflammation is now considered as one of the crucial carcinogenic factors 12,25 , genetic variability in inflammation-associated TLR genes has revealed their potential role in influencing the susceptibility to pathogenic infections and development of cancer 10 . Of the various TLRs, TLR4 is known to recognise exogenous ligands such as lipopolysaccharide (LPS), fusion (F) protein of respiratory syncytial virus as well as endogenous ligand like heat shock proteins (HSP60, HSP70) and high mobility group box 1 (HMGB1) [26][27][28] , whereas TLR9 recognizes unmethylated CpG-rich bacterial and viral DNA 29 .
Reports on the influence of TLR4 and TLR9 single nucleotide polymorphisms (SNPs) in cervical cancer susceptibility are limited as well as conflicting [30][31][32] . In the case of TLR4 polymorphisms, Asp299Gly (rs4986790) and Thr399Ile (rs4986791) were shown to be associated with tumor progression, however, no direct association of these SNPs was found in case-control set up 33,34 . Among the common TLR9 polymorphisms -1486 T/C (rs187084) and C2848T (rs352140) polymorphisms were found to be the risk factors for cervical cancer [35][36][37] . Conversely, a study by Pandey et al. 38 reported no association of TLR9 C2848T polymorphism with cervical cancer, however, the same SNP was marginally associated with advanced cancer stages. Jin et al. 39 reported a significant difference in the distribution of minor alleles of TLR4 3′ UTR SNP rs7873784 C/G and TLR9 SNP G2848A in cervical cancer and HPV positive cases. However, in the same study group, the other TLR4 SNPs (rs4986791, rs11536889) were not associated with cervical cancer.
Considering the importance of chronic inflammation in carcinogenesis as well as the influence of TLR genes' polymorphisms in inflammation and cancer susceptibility, the present study was designed to investigate the role of four TLR4 (rs4986790, rs10759931, rs11536889 and rs1927911) and equal number of TLR9 (rs187084, rs5743836, rs352140 and rs352139) SNPs in HPV infection and cervical cancer susceptibility.

Results
Clinico-demographic characteristics. Mean age of cervical cancer patients (52.4 ± 11.6 years) and controls (51.8 ± 11.8 years) was comparable without any statistically significant difference (p = 0.625). However, features such as age at marriage (p < 0.001), age at first childbirth (p < 0.001) and parity (p < 0.0001) showed statistically significant difference between the cases and controls. All the cervical cancer cases were histopathologically diagnosed as squamous cell carcinoma type. Clinical staging of cervical cancer biopsies was performed as per the FIGO guidelines that revealed 9 (8.2%), 39 (35.5%), 55 (50%) and 7 (6.3%) patients in Stage I, II, III and IV respectively. The detailed demographic and clinicopathologic features of patients are presented in Table 1 www.nature.com/scientificreports www.nature.com/scientificreports/ combined frequency of HPV 16 and 18 was found to be 68% (75/110). Moreover, two out of 141 control subjects (1.4%) were also detected positive for HPV consensus sequences, of which one (0.7%) carried HPV16 DNA.
Genotype distributions. All the TLR4 and TLR9 SNPs within the control population were in agreement with the Hardy-Weinberg equilibrium except for TLR4 SNP rs11536889. However, the polymorphism was retained as its homozygous genotype GG was not detected in any of the study subjects which could be a probable reason for its deviation from the Hardy-Weinberg equilibrium.
A significant difference in the distribution of genotype frequencies between the cases and the controls were observed for TLR4 SNPs rs11536889 (p = 0.013) and rs1927911 (p = 0.04) as well as TLR9 SNPs rs187084 (p = 0.01) and rs352139 (p = 0.04). The distribution of genotypes for TLR4 and TLR9 are shown in Supplementary Table S1  )] and the major allele of TLR9 SNP rs187084 were also varied significantly between patients and controls, conferring their association with the cervical cancer risk. Genotypic and allelic association between TLR4 and TLR9 variants and cervical cancer risk is presented in Table 3. A comparative analysis between early (stage I + II) and late (stage III + IV) stages revealed heterozygous genotypes of TLR9 rs187084 [p = 0.011, age-adjusted OR = 0.283 (0.107-0.749)] and rs352140 [p = 0.015, age-adjusted OR = 0.304 (0.117-0.790)] to be associated with early stage cervical cancer. However, none of the TLR4 SNPs shown significant association with early or late stages of cancer (Table 4).

Haplotype analysis. Linkage disequilibrium (LD) analysis revealed two SNPs of each TLR4
(rs10759931 aka rs11536858, rs1927911) and TLR9 (rs352139, rs187084) genes in strong LD (Fig. 1). The haplotypes were generated using the four SNPs of each TLR4 and TLR9 genes among the cases and controls (Table 5). Six common haplotype of TLR4 (frequency > 5%) and TLR9 (frequency > 2.5%) showed an accumulated frequency of 86.1% and 79.6% respectively in controls. Among HPV 16 and 18 positive patients TLR4 and TLR9 haplotypes revealed an accumulated frequency of 85.5% and 84.7% respectively. Distribution of TLR4 haplotypes differed significantly in HPV 16 and 18 infected cases (Pglobal = 0.045) as compared to control, whereas no such difference was detected while evaluating the TLR9 haplotypes (Pglobal = 0.493) ( Among cervical cancer cases, TLR4 and TLR9 haplotypes revealed an accumulated frequency of 85% and 83.1% respectively. Results of the global test score showed a significant difference in haplotype distribution between patients and controls in the case of TLR4 variants (Pglobal = 0.0033), while no significant difference was obtained for TLR9 variants (Pglobal = 0.227) ( Table 5). Furthermore, the TLR4 haplotype GTAC [p = 0.047, OR = 1.77 (1.00-3.13)] and TLR9 haplotype GATC [p = 0.019, OR = 3.95 (1.15-13.50)] were found to be associated with the increased risk of cervical cancer whereas the TLR4 haplotype GCAG [p = 0.0076, OR = 0.39 (0.19-0.79)] was significantly associated with decreased risk of cervical cancer. Furthermore, within cases, haplotypes analysis did not reveal an association of either TLR4 (Pglobal = 0.733) or TLR9 (Pglobal = 0.546) haplotypes with the early or late stages of cervical cancer (Supplementary Table S2).

Discussion
The influence of TLR polymorphisms is gradually increasing in the field of biomarkers study in various diseases including cancer 10 . In the present study, we investigated the role of the common TLR4 and TLR9 SNPs in susceptibility to HPV infection and cervical cancer among the study subjects from Gujarat, India. Considering the influence of hrHPVs in cervical carcinogenesis, we first analyzed the prevalence of two major hrHPVs HPV 16 and 18 that revealed a frequency of 68% as compared to nearly 71% and 78% prevalence globally as well as in India respectively 40 . However, a previous report 41 , from the same geographic region as of ours found 60% of the patients to be infected with HPV 16 and 18. The difference in the percentage of hrHPV detection, though not very high, can be attributed to the variation in the sample size as the number of patients in the present study were more than double as reported by Patel et al. 41 . A higher prevalence of approximately 21% HPV infection other than HPV 16 and 18 in our study subjects highlights the necessity of genotyping other hrHPVs to identify additional prevailing HPVs.
We further analyzed polymorphisms present in UTRs, exons, and introns of TLR4 and TLR9 genes. The variations in UTRs are known to influence ribosome recognition, termination and post-transcriptional modification which may alter the expression and functionality of a particular protein 42 . We found a mixed association of different 3′ UTR and 5′ UTR SNPs of TLR4 and TLR9 genes in our study subjects, suggesting a probable role of these SNPs in disease susceptibility.
TLR9 promoter SNP rs187084 (-1486T/C) conferred a increased risk to HPV 16 and 18 infection and cervical cancer. A similar association of TLR9 rs187084 polymorphism with an increased risk of cervical cancer was reported among Polish and Chinese women 35,36 . Our results on TLR9 rs187084 polymorphism are in good agreement with the recent meta-analyses 30,31 that supported a significant role of rs187084 in cervical cancer risk. Within cases, TLR9 rs187084 showed over presentation in early-stage cancer compared with late stages. Interestingly, we did not find an association of another TLR9 promoter SNP rs5743836 (−1237T/C) with HPV infection and/ or cervical cancer risk. Our result supports the observation of Oliveira et al. 43 who reported no association of TLR9 promoter SNP rs5743836 with HPV infection or clearance in healthy Brazilian women. Even though no direct role of TLR9 promoter SNPs has been reported in cervical cancer, the T allele of TLR9 promoter SNP rs187084 (−1486T/C) together with G allele of intronic rs352139 A/G SNP have been suggested to down regulate TLR9 expression in systemic lupus erythematosus 44 . The T allele of rs5743836 (−1237T/C) has been suggested to be associated with high basal promoter activity 45 and C allele with higher affinity to NF-κB binding, causing increased production of proinflammatory cytokines 46 .
With regard to TLR4 promoter SNP rs10759931, no association was observed either with HPV infection or cervical cancer risk. However, the same SNP has been reported to be associated with prostate and gastric www.nature.com/scientificreports www.nature.com/scientificreports/ cancers risk 47,48 . The homozygous AA genotype of TLR4 rs10759931 has been reported to be associated with high TLR4 expression in symptomatic atherosclerotic patients compared to non-symptomatic and healthy individuals carrying GG or GA genotypes 49 . They found that the two alleles of rs10759931 differ in their binding affinity to GATA-2 transcriptional factor. Furthermore, we observed the 3′ UTR heterozygous genotype GC of TLR4 rs11536889 to be associated with increased risk of cervical cancer in our study subjects. A similar observation was found in bladder cancer 50 , however, the association status of this SNP with other cancers was inconsistent 32 . Moreover, the G allele of TLR4 rs11536889 3′ UTR SNP has been suggested to play a key role in inhibiting TLR4 translation in monocytes 51 . However, expression analysis of TLR4 and TLR9 genes may provide more insights into the functional role of these UTR SNPs in cervical cancer risk.
Additionally, we analyzed a synonymous and a non-synonymous SNP of TLR9 and TLR4 genes respectively. Even though a synonymous change does not alter incorporation of amino acid, it has been observed that such SNPs can alter mRNA splicing, stability, and structure as well as protein folding thereby affecting the function of the subsequent protein 52 . We did not find a significant effect of TLR9 synonymous SNP rs352140 (G2848A; Pro545Pro) with cervical cancer risk which is in good agreement with a recent meta-analysis by Tian et al. 30 . An association of G2848A SNP with early stages of cervical cancer was detected in our study subjects which is in contrast to the report of Pandey et al. 38 who observed an association of the same SNP with the late stage  www.nature.com/scientificreports www.nature.com/scientificreports/ cervical cancer in North Indian women. However, Roszak et al. 35 reported an association of C2848T SNP along with -1486T/C SNP with cervical cancer risk in the Polish population. Similarly, the Han Chinese women carrying TLR9 rs352140 (G2848A) GA/AA genotype along with HPV16 infection showed an increased risk of cervical cancer compared to women with GG genotype 35,53 .
With regard to non-synonymous SNP rs4986790 (A896G; Asp299Gly) of TLR4, intriguingly, we found the heterozygous AG genotype (Asp/Gly) to be strongly linked to HPV 16/18 infection suggesting a queering effect of the amino acid change as no interaction of HPV capsid proteins with TLR4 is known yet. The amino acid change is reported to affect van der Waals interaction and hydrogen bonding in the leucine-rich repeats of TLR4, thereby modulating its surface properties that may affect the binding of TLR4 ligand such as LPS 54 . Although HPV is not a known TLR4 ligand, our paradoxical observation warrants a meticulous investigation. Furthermore, we observed a significant association of minor allele G (Gly) of Asp299Gly polymorphism with cervical cancer risk, however, no genotypic association was found. Similarly, in North Indian women, no association of TLR4 Asp299Gly polymorphism, in addition to another common TLR4 Thr399Ile polymorphism with cervical cancer risk was observed by Pandey et al. 33 . Moreover, Asp299Gly polymorphism has been found to be contradictorily associated with different cancer types including cervical cancer 32 .
A growing body of evidence suggests a potential role of intronic SNPs located either in exon/ intron boundaries, intron splice enhancer, branchpoint site or outside the exon-intron splice junctions in regulating gene expression 55 . It has also been observed that intronic SNPs in one gene can affect the expression of a far located gene 55 . Congruously, we observed a significant difference in the distribution of genotypes of TLR9 intronic rs352139 A/G SNP between cases and controls, however, none of its genotypes or allele was associated with cervical cancer risk. On the other hand, the heterozygous genotype of TLR4 intronic rs1927911 SNP was significantly associated with cervical cancer risk which is in agreement with the observation of Song et al. 47 in prostate cancer. However, in hepatocellular carcinoma, the same SNP showed a protective effect 56 .
As haplotypes are considered more informative than SNPs 57 , we generated haplotypes from different combinations of TLR4 and TLR9 SNPs. The TLR4 haplotype GTAC was linked with a significant increase in cervical cancer risk in addition to the TLR9 haplotype GATC that also showed association with increased HPV 16 and 18 infections. Intriguingly, another TLR4 haplotype GCAG showed a significant association with decreased cervical cancer risk as well as acquiring the hrHPV infection, suggesting its protective role. Moreover, to understand the influence of TLR4 and TLR9 haplotypes on tumor progression, we correlated the haplotypes with early (I and II) and late (III and IV) tumor stages. However, none of the haplotypes showed association with clinical aggressiveness. Since these haplotypes included both risk as well as protective alleles, a crucial role of TLR4 and TLR9 polymorphisms may be envisaged towards HPV infection and cervical cancer susceptibility.
To identify the strong coinheritance of the SNPs we calculated linkage disequilibrium of TLR4 and TLR9 SNPs, wherein TLR4 rs10759931 and rs1927911, and TLR9 rs187084 and rs352139 were in strong LD, evincing strong influence of these inherited variations in cervical cancer. Intriguingly, we observed that in both the genes strong LD was detected between SNPs of 5′ UTR and the first intron only. Conceptually there should be a decrease in linkage disequilibrium with a decrease in distance between two loci. However, our study revealed SNP pairs in both TLR4 and TLR9 genes that did not follow the standard notion. For example, in TLR4, SNP pair rs10759931:rs4986790 with a distance of 11.1 Kb showed strong LD (D′ = 0.54) as compared to another SNP pair rs4986790:rs11536889 that had a shorter distance of 2.8Kb (D′ = 0.12). Similarly, TLR9 SNP pair rs352140:rs187084 (distance = 4.3 kb) was in strong LD (D′ = 0.5) compared to SNP pair rs5743836:rs187084 (D′ = 0.04) with shorter distance of 0.24 kb among them. Our LD analysis is in agreement with the observations of Stephens et al. 57 who suggested that distance between the SNPs does not have a significant impact on the level of LD. Various SNP pairs of TLR4 and TLR9 genes, their genetic distance and D′ values are shown in Supplementary Table S3.
Although our results suggest a significant role of TLR4 and TLR9 polymorphisms in cervical cancer, the study has some vital limitations too. Firstly, the selection bias cannot be excluded as it was a hospital-based case-control study, Moreover, the size of the study population needed augmentation to increase the statistical power, which is one of the major limiting factors among the numerous cancer case-controls studies worldwide. Additionally, in vivo expression analysis would have reflected the effect of SNPs on the expression pattern of TLR4 and TLR9.
To our knowledge, this is the first comprehensive analysis of TLR4 and TLR9 SNPs and haplotypes to understand their role in cervical cancer. Our results suggest moderate to strong impact of TLR4 and TLR9 polymorphisms in susceptibility to hrHPV infection and cervical cancer. Additional research on large and varied ethnic populations is warranted to precisely understand the impact of both the genes in HPV infection and cervical cancer risk.
Methods study subjects. The study comprised of 110 untreated cervical cancer patients and 141 healthy controls recruited from 2012 to 2017; from Shree Krishna Hospital, Karamsad, Anand; and Sir Sayajirao General Hospital and Medical College, Vadodara, India. The sample types included primary histopathologically diagnosed cervical cancer biopsies and cytologically confirmed normal cervical smears from healthy controls. The clinical staging of cervical cancer samples was done as per The International Federation of Gynecology and Obstetrics (FIGO) recommendations. The study subjects belonging to Gujarati ethnicity were comparable in age and non-relatives of each other. The patients manifesting multiple cancers and those who underwent radiation or chemotherapy were excluded from the study. The inclusion criteria of healthy controls included the absence of cancer history in family and cervix related disorders such as cervicitis, warts, pre-cancerous and cancerous lesions. Additionally, sample collection was avoided from the women undergoing menstruation. All experiments were performed in accordance with the relevant guidelines and regulations. The study was approved by the Institutional Review Board, Ashok and Rita Patel Institute of Physiotherapy, CHARUSAT, Changa, Anand; Institutional Ethics Committee, HP Patel Centre for Medical Care and Education, Karamsad and Institutional Ethics Committee for Human www.nature.com/scientificreports www.nature.com/scientificreports/ Research (IECHR) Medical College and SSG Hospital, Vadodara, India. Written informed consent was obtained from all the study subjects. DNA extraction. The samples were collected in chilled phosphate buffered saline and were either processed immediately or stored at −20 °C till further processing. DNA was isolated using standard Proteinase-K phenol-chloroform extraction method. In the case of a low number of cervical cells, spin-column based DNA isolation kit (NucleoSpin Tissue, Macherey-Nagel, Germany) was utilized. The quality and quantity of extracted DNA were determined using ethidium bromide-stained 1% agarose gel on a GelDoc system (BioRad, USA) and NanoDrop 2000 (Thermofisher, USA).
HpV detection. HPV detection was first carried out using consensus Gp5+/Gp6+ primers followed by type-specific primers for the detection of hrHPV 16   www.nature.com/scientificreports www.nature.com/scientificreports/ primer and reverse primer, 1X ROX reference Dye II and 25 ng of template DNA. The positive controls for HPV 16 and 18 were obtained as a part of participation in the Global HPV Proficiency Study, Equalis, Uppsala, Sweden. β-globin gene served as an internal control while in the negative control DNA was replaced with PCR grade nuclease-free water. All the reactions were performed in duplicates. Touchdown thermal profile for HPV detection by consensus primers and thermal cycling conditions for HPV 16 and 18 detections along with the details of primer sequence and amplicon size is mentioned in Supplementary Table S4.

Genotype analyses.
A total of eight SNPs, four each of TLR4 (rs4986790, rs10759931, rs11536889, rs1927911) and TLR9 (rs187084, rs5743836, rs352140, rs352139) genes were analyzed either using Polymerase Chain Reaction and Restriction Fragment Length Polymorphism (PCR-RFLP) or Allele-Specific PCR (AS-PCR). The selection of SNPs was carried out using SNP database of NCBI (https://www.ncbi.nlm.nih.gov/snp/). The SNPs were selected on the basis of (1) Genetic region: In this criteria the SNPs were selected to cover different regions of gene, for example, exon, intron and UTRs, (2) Global minor allele frequency: The SNPs with minor allele frequency > 5% were evaluated for association analysis (3) Frequent association of SNPs with different inflammation associated cancers: To fulfil the above criteria literature survey was conducted using PubMed and random web search. The characteristics of TLR4 and TLR9 SNPs included in this study are shown in Supplementary Table S5. Sequences of primers specific for each SNP, amplicon size and thermal profile is mentioned in Supplementary Table S6. A typical PCR of 25 µl contained 50 to 100 ng genomic DNA, 0.1 mM dNTP mix, 0.1 µM of each oligonucleotide primer and 0.8U Taq DNA polymerase (Kapabiosystems, USA). All the reactions were performed on an MJ Mini Thermal Cycler (BioRad, USA). Except for TLR9 rs352139 polymorphism that was genotyped using AS-PCR, the rest of the SNPs were subjected to restriction digestion using 5U of respective restriction enzymes procured from New England Biolabs, USA. For the identification of SNPs by RFLP, the associated restriction enzymes, incubation temperature and time, digested products, genotypes and mode of www.nature.com/scientificreports www.nature.com/scientificreports/ visualization is detailed in Supplementary Table S7. The amplified, as well as restriction digested products, were visualized on a GelDoc system (BioRad, USA). statistical analysis. Alterations in demographic features among cases and controls were compared using student t-test and chi-square test for continuous and categorical variables respectively. Age of study subjects was expressed as mean ± standard deviation. Deviation from Hardy-Weinberg equilibrium was determined by the χ 2 goodness-of-fit test. Pearson's χ 2 test was used to evaluate the difference of the SNP distribution among cases and controls. Genotypic and allelic association of SNPs with the disease were estimated using χ 2 and Fisher's exact test. Unconditional logistic regression analysis was performed to compute age-adjusted odds ratio (OR). All the statistical analysis was performed on the Statistical Package for Social Sciences version 24.0 (SPSS, USA). Tests of statistical significance were two-sided and taken as significant when the p-value was less than 0.05. Haplotype block structure and linkage disequilibrium (LD) structure were determined by Haploview (v4.2) and Locusview (v2.0). The D′ values were computed using the default algorithm created by Gabriel et al. 58 at 95% confidence interval. Haplotypes were estimated using an accelerated EM algorithm similar to the partition/ ligation method as described by Qin et al. 59 . Sum of the fractional likelihoods of each individual for each haplotype was used to obtain a count for case-control association tests. Global score test was performed using FAMHAP software v19 to evaluate the differences in haplotype frequency distribution among cases and controls. Association of the individual haplotype with cervical cancer as well as HPV infection was measured by the χ 2 test.

Data Availability
All data generated or analysed during this study are included in this published article (and its Supplementary Information files). www.nature.com/scientificreports www.nature.com/scientificreports/ www.nature.com/scientificreports www.nature.com/scientificreports/