Background: The predictive value of codon 72 arginine homozygosity at the p53 gene for human papilloma virus associated cervical cancer risk remains inconclusive. It has also been proposed that the inheritance of specific germline haplotypes based on three biallelic polymorphisms of p53 (intron 3 16 bp duplication, codon 72 Bst UI (Arg/Pro), and intron 6 Nci I restriction fragment length polymorphism at nucleotide 13494) is a better predictor of various cancer risks.
Aims: To determine the genotype and haplotype frequency of these three p53 polymorphisms in 61 patients with cervical squamous cell carcinoma and 94 ethnically matched controls from the eastern region of India and estimate the risk, if any, of specific genotypes and haplotypes.
Methods: Samples were genotyped by polymerase chain reaction followed by variant specific restriction enzyme digestion. Haplotypes were estimated by the maximum likelihood method using the expectation maximisation algorithm.
Results: Genotype distributions of the three polymorphisms in patients and controls showed a good fit to the Hardy-Weinberg equilibrium. The p53 codon 72 arginine homozygous genotype was significantly over represented in patients compared with controls. Those with the homozygous arginine genotype exhibited a 2.59 fold higher risk of developing squamous cell carcinoma of the uterine cervix. A significant risk was also seen with a combination of two haplotypes, 1–2–1 and 1–2–2.
Conclusion: p53 codon 72 arginine homozygotes appear to be at greater risk of developing squamous cell carcinoma of the uterine cervix. The high risk haplotypes 1–2–1 and 1–2–2 also contain the arginine allele, further strengthening this conclusion.
- CI, confidence interval
- df, degrees of freedom
- HPV, human papillomavirus
- nt, nucleotides
- OR, odds ratio
- PCR, polymerase chain reaction
- RFLP, restriction fragment length polymorphism
- cervical cancer
Statistics from Altmetric.com
- CI, confidence interval
- df, degrees of freedom
- HPV, human papillomavirus
- nt, nucleotides
- OR, odds ratio
- PCR, polymerase chain reaction
- RFLP, restriction fragment length polymorphism
Cervical carcinoma in women is the second most common type of cancer worldwide. The incidence is greater in developing countries, such as India, where this is the most common female malignancy, with almost 100 000 new cases each year.1 This disease is one of the major causes of death from cancer in women, with an overall survival rate of 40%.2 Human papillomavirus (HPV) is considered a necessary cause of all cervical cancers.3 HPV positivity in cervical cancer is estimated to be between 90% and 95%.4 Fifteen HPV types have been classified as high risk types (16, 18, 31, 33, 35, 39, 45, 51, 52, 56, 58, 59, 68, 73, and 82); three have been classified as probable high risk types (26, 53, and 66); and 12 have been classified as low risk types (6, 11, 40, 42, 43, 44, 54, 61, 70, 72, 81, and CP6108).4 Among the high risk HPV genotypes, HPV-16 is the most common in squamous cell carcinoma of the cervix (50–60%), followed by HPV-18, which is present in about 11–15% of cervical cancer cases.4,5 The viral oncoproteins E6 and E7 of HPV mainly contribute to tumorigenesis.5,6 The E6 protein interacts with p53 and inhibits its activity, followed by proteolytic degradation through the ubiquitin pathway.7,8 In addition to inducing the rapid degradation of p53, E6 also binds to and degrades FADD, thus preventing the transmission of apoptotic signals via the Fas pathway.6 Recently it was shown that sustained inactivation of the p53 pathway by the E6 protein is required for maintenance of the proliferative phenotype of HeLa cervical carcinoma cells.9 Because mutation of the p53 gene is a relatively rare event,10 abrogation of p53 function by the E6 protein of HPV is thought to be one of the major events in cervical carcinogenesis.
The possible contribution of inherited polymorphisms in the p53 gene to the development of cancer risk has been the subject of interest after a report that women with the p53 codon 72 arginine homozygous genotype were at higher risk than those with the proline homozygous genotype of developing cervical cancer.11 Since then, there have been several reports worldwide, but no definite conclusion has been reached. Some reports revealed a strong association, whereas others did not.12,13 Two other p53 polymorphisms—a 16 bp duplication at intron 3 and intron 6 Nci I restriction fragment length polymorphism (RFLP) at nucleotide (nt) 13494—have also been studied for their association with different cancer types. Genotypes homozygous for the absence of the Nci I restriction site in intron 6 and the presence of the 16 bp duplication at intron 3 were reported to have an increased risk of lung cancer14 and ovarian cancer.15 The 16 bp duplication at intron 3 was also found to increase the risk of breast cancer by the age of 50 years for both homozygous and carrier genotypes,16 although another study reported no increase in lung cancer risk associated with this particular genotype.17
“There are no reports on the analysis of cervical cancer risk using multiple p53 gene polymorphisms and haplotypes constructed from them”
It has also been proposed that the inheritance of specific germline haplotypes based on these three biallelic polymorphisms of p53 (16 bp duplication in intron 3, codon 72 Arg/Pro, and intron 6 Nci I RFLP at nt 13494) is a better predictor of breast cancer,18–21 colorectal cancer,22 and lung cancer.14 The specific haplotype “p53 1–2–1” exhibited a positive association with breast cancer risk.18–20
To date, to the best of our knowledge, there are no reports on the analysis of cervical cancer risk using multiple p53 gene polymorphisms and haplotypes constructed from them. Therefore, we sought to examine the allele and genotype frequencies for all three p53 polymorphisms (intron 3 16 bp duplication, codon 72 Bst UI (Arg/Pro), and intron 6 Nci I RFLP at nt 13494) in 61 patients with cervical cancer and 94 ethnically matched healthy female controls from eastern India. In our study, haplotypes constructed from these three polymorphisms exhibited distinct differences between the case and control subjects. The association of specific p53 genotypes and/or haplotypes with cervical cancers was examined.
MATERIALS AND METHODS
Sample collection and clinical data
Sixty one freshly operated uterine cervical tissues and peripheral blood leucocytes from the patients were obtained from the hospital section of Chittaranjan National Cancer Institute, Kolkata, India. Three of the samples were dysplastic and the rest were squamous cell carcinoma type, of which 39 were stages I/II and 19 were stages III/IV. Peripheral blood leucocytes were also obtained from 94 ethnically matched, unrelated, normal, healthy female subjects with no history of cervical lesions recruited through a population screening programme for use as controls. Both patients and controls were from the Bengali speaking caste population of the same city, with a similar low income background. The mean (SD) age of the patients was 46.89 (10.29) years and that of the controls was 47.44 (12.53) years. The tissue samples were frozen immediately after collection and stored at −80°C until use. Written consent was obtained from the patients and controls before sample collection. The internal review committee on research using human subjects cleared the project after careful examination of the protocols. High molecular weight genomic DNA was extracted from the peripheral blood leucocytes and from tissues after microdissection according to the standard procedure.23
HPV detection and typing
The presence of HPV in the cervical lesions was detected by performing the polymerase chain reaction (PCR) using primers (MY09 and MY11) from the consensus L1 region.24 Typing of HPV-16/18 in the L1 positive samples was done by means of PCR using specific primers from the E6 region of HPV-1625 and the E7 region of HPV-18.26 The PCR products were electrophoresed on a 2% agarose gel, stained with ethidium bromide, visualised under ultraviolet light, and photographed. For final confirmation of the HPV types, after gel electrophoresis the PCR products were transferred on to a nylon membrane for Southern hybridisation with [32P] labelled HPV type specific probes.27 DNA from the SiHa (for HPV-16) and HeLa (for HPV-18) cell lines and the HPV type specific plasmids were used as positive controls.
p53 gene polymorphism analysis
Three biallelic p53 gene polymorphisms (intron 3 16 bp duplication, codon 72 Bst UI (Arg/Pro), and intron 6 Nci I at nt 13494) for both patients with cervical cancer and controls were analysed by PCR, followed by variant specific restriction enzyme digestion.28 PCR was performed with primers that flanked the two polymorphic sites of the intron 3 16 bp duplication and the exon 4 Bst UI RFLP, as described elsewhere.28 The resulting PCR products were either 432 or 448 bp DNA fragments depending on the absence or presence of the intron 3 16 bp duplication in the template genomic DNA. Length polymorphism was directly evident from 6% polyacrylamide gel analysis of the PCR product. An aliquot of the same PCR product was subjected to restriction digestion with Bst UI (New England Biolabs Inc, Beverly, Massachusetts, USA), which generated DNA fragments of four different sizes: 448 bp, Bst UI digestion resistant (Pro at codon 72) with the intron 3 duplication; 432 bp, Bst UI digestion resistant (Pro at codon 72) without the intron 3 duplication; 246 bp, Bst UI digested (Arg at codon 72) with the intron 3 duplication; 230 bp Bst UI digested (Arg at codon 72) without the intron 3 duplication.
The third polymorphic site (Nci I in intron 6) was PCR amplified separately with flanking primers, as described elsewhere.28 Aliquots (5 μl) of the reaction mixture containing the 913 bp PCR fragment were subjected to restriction digestion with Nci I. The digest contained either a 563 bp DNA fragment (in the absence of the polymorphic Nci I site) or a combination of 286 bp and 277 bp DNA fragments (in the presence of the polymorphic Nci I site), together with an invariant 350 bp DNA fragment as a result of the presence of a non-polymorphic Nci I site in the amplified region of the DNA.
Some of the PCR products (5%) were sequenced (ABI prism 377; Applied Biosystems, Foster City, California, USA) to confirm the genotypes at all loci, which were determined by PCR and PCR–RFLP methods.
Genotypic distributions were examined for a significant departure from the Hardy-Weinberg equilibrium by a goodness of fit χ2 test. A comparison of the haplotype and genotype distributions was performed by χ2 statistics with Yate’s correction whenever required. The frequencies of extended haplotypes were estimated by the maximum likelihood method, using the expectation maximisation algorithm as implemented in the ARLEQUIN software program. Odds ratios (ORs) with 95% confidence intervals (CI) were computed at http://www.hutchon.freeserve.co.uk.
Table 1 depicts the genotype frequency distribution of three different biallelic DNA polymorphisms in the p53 gene for the patients with cervical cancer and the ethnically matched normal female control subjects. The codon 72 Bst UI A1 allele (proline allele) and the intron 6 Nci I A1 alleles are defined by absence of the restriction sites, and the intron 3 16 bp duplication A1 allele is defined by absence of the duplication.17,29 Both the study groups were found to be in Hardy-Weinberg equilibrium at all three p53 loci (table 1). We found a significant difference in the overall genotype frequency distribution only at the codon 72 locus between the case and control subjects (p = 0.038; table 1). Individuals with the arginine homozygous genotype exhibited a 2.59 fold higher risk of developing squamous cell carcinoma of the cervix compared with those with the other two genotypes (OR, 2.59; 95% CI, 1.18 to 5.67). The difference in the numbers of samples genotyped, particularly for the Nci I locus, results from the fact that PCR amplification failed in many of the control subjects, even after repeated attempts.
Table 2 shows the p53 codon 72 genotype distribution in normal controls and patients with HPV positive cervical carcinoma at different histological stages. Overall, we found that 52 of the 61 tumours were HPV positive; 37 and 15 were HPV-16 and HPV-18 positive, respectively. In contrast to a previous report,30 we found no coinfection with HPV-16 and HPV-18 in our samples. One possible explanation for this discrepancy is that we isolated the tumour DNA after microdissection of the tumour tissues, thus enriching the DNA specifically from tumour tissue and minimising the possibility of contamination with DNA from adjacent tissues. It has been reported that tumours are clonal with respect to HPV types, but surrounding tissues may be infected by different subtypes, thus leading to cases of coinfection.31 It is apparent from table 2 that the frequency of HPV-16 positive tumours increases with advancing histological stages of the tumours. We can make no conclusions regarding early HPV infection because of the small number of dysplastic samples. However, a similar increase in HPV infection with progression of histological stages of dysplasia has been reported.30 The frequency of the arginine homozygous genotype was significantly higher in HPV associated cervical cancer cases (30.0%) than in normal controls (15.0%), giving an overall 2.86 fold risk for this genotype (OR, 2.86; 95% CI, 1.27 to 6.44). This risk was increased further when samples were subclassified with respect to tumour stage. Tumours at stage III+IV showed a higher risk (OR, 4.15; 95% CI, 1.42 to 12.2) than those at stage I+II, which showed a moderate risk (OR, 2.24; 95% CI, 0.91 to 5.5). However, the difference in genotypic distribution between HPV positive and HPV negative cases was not significant (χ2 = 0.69; degrees of freedom (df) = 1; p = 0.405). This is probably because very few of the cervical tumours were HPV negative.
Table 3 shows the estimated haplotype frequencies among the three p53 polymorphisms in the samples from patients and controls. The p53 haplotypes are expressed in the order: intron 3–codon 72–intron 6.32 Among the eight possible haplotypes, five were found in normal samples and six in the patients. The 1–2–2 haplotype was most common in samples from both the patients and the controls, whereas the 1–2–1 haplotype was completely missing in the normal controls. There was a significant difference in overall haplotype frequency distribution between the two groups (χ2 = 11.52; df = 5; p = 0.042), but no particular haplotype was found to be associated with the patients. However, when haplotypes 1–2–1 and 1–2–2 were combined, they were found to be associated with cervical carcinogenesis (OR = 1.77; 95% CI, 1.08 to 2.89).
Since it was first reported in 1998,11 the hypothesis that p53 codon 72 polymorphism is a predictor of increased HPV associated cervical cancer risk has generated a considerable amount of controversy.11 Extensive studies have been carried out throughout the world to reveal the role of the homozygous Arg/Arg genotype as a risk factor for HPV associated cervical carcinogenesis. In populations from southeast Asia—such as Taiwan, Thailand, Korea, Japan, Hong Kong, and China—the homozygous Arg/Arg genotype was not found to be a risk for the development of cervical cancer.33–38 Among the Hispanic population from Brazil, Chili, Peru, and in the African population, the Arg/Arg genotype was found to be a risk factor for HPV associated cervical carcinogenesis,39–42 although one study from Argentina failed to find an association.43 Reports from white populations, particularly in Italy and Sweden, are contradictory.44–46 Studies from France, the Czech Republic, Germany, and the UK were unable to find an association,47–50 whereas investigations in Greece, Holland, and Hungary supported the role of the Arg/Arg genotype as a risk factor for HPV associated cervical carcinogenesis.51–53 To date, there are five reports from India examining the association between p53 codon 72 polymorphism and the risk of cervical cancer. Three of these reports found a positive association between the p53 arginine homozygous genotype at codon 72 and HPV associated cervical cancer,30,54,55 whereas one study from southern India failed to find an association.56 In addition, there is one report from eastern India stating that proline homozygosity at codon 72 of p53 is the at risk genotype.57 Because the reports from India were contradictory, we investigated the p53 codon 72 polymorphism and two other p53 polymorphisms in patients with cervical cancer from eastern India. There are two recent meta-analyses of the association between the p53 codon 72 polymorphism and cervical cancer. One of these analyses found a slightly increased risk for adenocarcinoma, but not for squamous cell cancer,13 suggesting that improvement of the predictive usefulness of the p53 codon 72 polymorphism requires more studies that are designed to address the genetic architecture of the p53 codon 72 polymorphism and other risk factors for cervical cancer. Another study found no association between this polymorphism and preinvasive cervical lesions, or invasive cervical cancer of undefined histology, but found a slightly increased risk for squamous cell cancer and adenocarcinoma.12 The authors suggested that two reasons for the contradictory results between these studies are methodological issues in the study designs and confounding by ethnicity. In our present study, both the controls and the patients were taken from the same caste ethnic group, living in the same geographical location.
“The biochemical differences between the two p53 genotypes at codon 72 were elegantly demonstrated by a study showing that the arginine form of the protein was much more susceptible to HPV E6 mediated degradation than the proline form”
Previously, we reported that various subpopulations within the caste group from eastern India do not differ with respect to allele frequency at several other polymorphic loci.58 Thus, the possibility of false associations as a result of population stratification has been ruled out. In addition, in our present study we have included only patients with squamous cell carcinoma of the uterine cervix and have analysed multiple polymorphisms simultaneously. The frequencies of the p53 codon 72 genotypes vary among populations from different countries. In a previous study comprising different ethnic groups, the frequency of arginine homozygotes varied from 12.4% to 53.4% across the world.32 Within India, the frequency of arginine homozygous individuals varied from 14% to 25% in different geographical regions of the country.28,55,56,59
We found a significant over representation of arginine homozygotes at codon 72 of the p53 gene in the patients with cervical cancer. Genotype differences between the case and control individuals at the two other loci were not significant. In our study, 52 of 61 of the patients were HPV-16/18 positive, with HPV-16 sequences being detected in 37 and HPV-18 in 15 of the HPV positive patients. This is in agreement with earlier studies from India.60–62 Arginine homozygotes showed a 2.59 fold higher risk of cervical carcinogenesis, which increased to 2.86 fold when the analysis was restricted to patients with HPV-16/18 positive cervical cancer and normal controls. It was interesting to note that HPV infection increased with the developmental stages of cancer and the presence of the homozygous arginine genotype (tumours at stage I+II: OR, 2.24; 95% CI, 0.91 to 5.5 and stage III+IV: OR, 4.15; 95% CI, 1.42 to 12.2). This supports the hypothesis that individuals with the homozygous arginine genotype are more susceptible to cervical carcinogenesis.
The biochemical differences between the two p53 genotypes at codon 72 were elegantly demonstrated by a study showing that the arginine form of the protein was much more susceptible to HPV E6 mediated degradation than the proline form.11 The regulatory role of intronic sequences is being increasingly recognised.63 Several workers have reported functional differences between the different p53 intron 3 and 6 polymorphisms.14,15,64,65 One study reported lower apoptotic indices and a reduced repair capacity in lymphoblastoid cell lines that harbour the less frequent alleles with the intron 3 16 bp duplication and intron 6 Nci I polymorphisms.14 There were no splicing errors linked to these polymorphisms.15 However, in cell culture analyses these two p53 intronic polymorphisms were not sufficient to impair p53 function during the selection process in neoplastic transformation, and an additional coding region mutation was needed.15 The intronic change of G → C at nt 13964 in intron 6 either caused or was linked to a variation leading to stabilisation and possible inactivation of the p53 protein, thereby contributing to tumorigenesis.65 Intron 4 of p53 has been predicted to contain several transcription factor binding sites,66 and single base pair substitutions in this intron have been shown to disturb the binding of unidentified proteins, resulting in decreased expression of p53.67 Inclusion of intron 3 in murine p53 cDNA has been shown to increase the expression of the gene.5 Overall, the functional role of the two p53 intronic polymorphisms analysed and their modulatory role in cancer risk remain uncertain. Further studies are needed to reveal any regulatory role of these intronic sequences. Alternatively, these intronic p53 polymorphisms may confer an increase in cancer risk through linkage disequilibrium with a functional variant of the p53 gene.
It has been proposed that the inheritance of specific germline haplotypes based on three biallelic polymorphisms of p53 (16 bp duplication in intron 3, codon 72 Arg/Pro, and intron 6 Nci I RFLP at nt 13494) is a better predictor of cancer risk.14,18–22,28 To date, there is no report on specific haplotypes associated with cervical cancer. In our present study, a significant difference in haplotype distribution between the patients with cervical cancer and the normal controls raises the possibility of an association of specific haplotypes with cancer risk (χ2 = 11.52; df = 5; p = 0.042). Although we were unable to ascertain a specific haplotype with increased cervical cancer risk in strict statistical terms, it is evident that the rare haplotype 1–2–1 is completely absent in the control subjects, but not in patients, and another haplotype 1–2–2 showed a modest risk (OR = 1.512; 95% CI, 0.92 to 2.47). This risk was more significant when both the haplotypes, 1–2–1 and 1–2–2, were combined (OR, 1.77; 95% CI, 1.08 to 2.89). However, more samples need to be analysed to support this observation. It is interesting to note that the two haplotypes are different only at the intron 6 Nci I polymorphic site, and both of them contain the arginine allele. This haplotype showed no association in other cancer types,14,18–21 which could be because of differences in the ethnicity of the populations and in cancer sites between these studies.
We suggest that individuals with the p53 codon 72 arginine homozygous genotype are at greater risk of developing squamous cell carcinoma of the uterine cervix. Evidence in favour of this observation is further strengthened by the fact that the specific haplotypes 1–2–2 and 1–2–1 showing a greater risk also contain the arginine allele. The association of these haplotypes could result either from strong linkage disequilibrium among these three sites, which we reported previously,28 or from the fact that the other two intronic polymorphisms contribute functionally to risk by altering the functional property of nascent RNA expression and/or stabilisation, as discussed above.
Take home messages
p53 codon 72 arginine homozygotes appear to be at greater risk of developing squamous cell carcinoma of the uterine cervix
The high risk haplotypes 1–2–1 and 1–2–2 also contain the arginine allele, further strengthening this conclusion
More samples need to be analysed to support this observation
We are thankful to the Director, Chittaranjan National Cancer Institute, Kolkata, India for providing access to patients. We are grateful to all members of the Human Genetics and Genomics Division of IICB and the Department of Oncogene Regulation of CNCI for their kind cooperation and encouragement during the study. SM, CM, and RS are grateful to the Council of Scientific and Industrial Research, New Delhi, India for predoctoral fellowships. This work was supported by research grants from the Department of Biotechnology (grant number BT/MB/05/002/94) and the Department of Science and Technology (grant number SP/SO/ D-75/96) of The Government of India.