Evaluating the association between p53 codon 72 Arg>pro polymorphism and risk of ovary cancer: a meta-analysis.

AIM: Allelic polymorphism in codon 72 of the p53 tumor suppressor gene causes imbalance of p53 protein expression. Earlier studies have shown association between allelic polymorphism in codon 72 of the p53 gene with risk of ovary cancer (OC); however the results are inconclusive and conflicting. Therefore, we performed this meta-analysis to investigate the relation between p53 codon 72 Arg>Pro polymorphism and overall OC susceptibility.
METHODS: We searched all eligible published studies based on the association between codon 72 of the p53 Arg>Pro polymorphism and risk of OC. Data were pooled together from individual studies and meta-analysis was performed. Pooled odds ratios (ORs) and 95% CI were calculated for allele contrast, homozygous, heterozygous, dominant and recessive genetic models.
RESULTS: A total of twelve studies comprising of 993 OC cases and 1264 healthy controls were included in this meta-analysis. Overall, no significant association was detected for Pro allele carrier (Pro vs. Arg: p = 0.916; OR = 0.980, 95% CI = 0.677 to 1.419), homozygous (Pro/Pro vs. Arg/Arg: p = 0.419; OR = 0.731, 95% CI = 0.341 to 1.564), heterozygous (Arg/Pro vs. Arg/Arg: p = 0.248; OR = 1.237, 95% CI = 0.862 to 1.773), dominant (Pro/Pro+Arg/Pro vsArg/Arg: p = 0.699; OR = 1.089, 95% CI = 0.706 to 1.681), and recessive (Pro/Pro vs Arg/Arg+Arg/Pro: p = 0.329; OR = 0.754, 95% CI = 0.428 to 1.329) genetic models, respectively. Also, in the stratified analysis by ethnicity, no significant association of this polymorphism with risk of OC was found in the Caucasian population.
CONCLUSIONS: This meta-analysis suggested that codon 72 of the p53 Arg>Pro polymorphism may not significantly contribute in ovary cancer susceptibility. However, future large studies with gene-gene and gene-environment interactions are needed to validate these findings.

Introduction

Ovary cancer (OC) is the most common carcinoma among females with poor prognosis. It is the sixth leading cause of death among gynecological malignancies in females worldwide [1], [2]. The etiology of OC is still unclear and epidemiological studies have suggested that susceptibility to OC of an individual is influenced by several genetic factors [3]. However, there is no thorough screening technique for this malignancy, indicating that the identification of a gene related to the risk of OC may improve the early diagnosis and prevention of this deadly disease.

The p53 tumor suppressor gene (TP53 at 17p13), recognized as “the guardian of the genome”, plays a significant role in the cell cycle arrest, senescence, DNA damage repair, regulates the cell cycle and requires loss of function mutations for tumor formation [4].

The p53 protein functions by interfering with central regulators of hypoxia which mediate angiogenesis, and eventually inhibit production of pro-angiogenic factors and endogenous angiogenesis inhibitors [5], [6], [7]. The ability of p53 to eliminate excess, damaged or infected cells by apoptosis is essential for the proper regulation of cell proliferation in multi-cellular organisms. Differential expression of p53 in various cancers and association of serum p53 levels with malignant tumors highlights the significance role of p53 in malignancy [8], [9]. The major modes of TP53 inactivation are single-base substitutions and loss of alleles, with inactivation by viral or cellular proteins [10].

Several polymorphisms have been detected in both coding and non-coding region of this gene [11]. An important single nucleotide germ line polymorphism in the proline- rich domain of exon 4 of p53 gene induces an arginine to proline residue change at amino acid position 72 [12]. The two polymorphic forms (Pro72 and Arg72) of p53 gene have different primary structures, electrophoretic migration and functional properties [13]. The arginine (Arg72) allele increases the ability of p53 to locate to mitochondria and induce cellular death, whereas proline allele (Pro72) impart a lower apoptotic potential and an increased cellular arrest in G1 phase of the cell cycle [14].

Considering the functional significance of p53 gene in carcinogenesis, it is speculated that codon 72 Arg>Pro polymorphism may be a potential susceptibility factor for OC. Lately, several epidemiological case-control studies have evaluated the association between p53codon 72 Arg>Pro polymorphisms and OC risk [15]–[26]. Despite several studies globally, the putative association between p53 codon 72 Arg>Pro genetic polymorphism and OC risk remains uncertain and lacks consensus. Therefore, to derive a more precise conclusion of the possible association between p53 codon 72 Arg>Pro polymorphism and OC risk, a meta-analysis was performed based on eligible published studies.

Materials and Methods

Publication search strategy

We carried out a PubMed (Medline), EMBASE and Google Scholar web database search covering all research articles published with a combination of the following key words: ‘p53 gene (polymorphism OR mutation OR variant) AND ovarian carcinoma or ovary cancer, tumor susceptibility (last updated on November 2013). All the searched studies were retrieved and their reference lists were checked as well for other relevant studies. When, more than one of the same population was included in several publications, only the most recent or complete study was included in this meta-analysis. Since, this is a meta-analysis of published articles based on the association of p53 codon 72 Arg>Pro polymorphism and OC risk, so ethical approval was not required for this study.

Inclusion and exclusion criteria

In order to minimize heterogeneity and ease the appropriate interpretation of this study, published articles included in the current meta-analysis had to meet all the following inclusion criteria: a) studies should have a cross-sectional, case-control or cohort design, b) must evaluated the association between p53 codon 72 Arg>Pro polymorphism and OC risk, c) recruited pathologically or histologically confirmed OC patients and healthy controls, d) have available genotype frequency in case and control, e) and published in English language. Also, when the case-control study was included by more than one research article using the same case series, we selected the study that included the largest number of individuals. On the other side, the major reasons for study exclusion were, overlapping of the data, case-only studies, review articles, and genotype frequencies or numbers are not reported. The study selection procedure has been shown in the form of flow-diagram as Figure S1 (PRISMA Flow Diagram).

Data extraction and quality assessment

For each retrieved research publication, the methodological quality assessment and data extraction were independently abstracted in duplicate using a standard protocol by two independent investigators. Data-collection form was used to guarantee the accuracy of the collected data by stringently following the inclusion-exclusion criteria mentioned above. The main characteristics abstracted from the retrieved studies included the name of the first author, publication year, the country of origin, the number of cases and controls, study type, and genotype frequencies. Cases associated with disagreement on any item of the data from the collected research studies were fully debated with investigators to achieve a final consensus.

Statistical analysis

In order to examine the relationship between p53 codon 72 Arg>Pro polymorphism and OC risk, pooled ORs and their corresponding 95% CIs were estimated. Heterogeneity assumption was examined by the chi-square-based Q-test [27]. Heterogeneity was considered significant when p-value<0.05. The data from single comparison was pooled using fixed effects model [28] when no heterogeneity presented. Otherwise, the random-effects model [29] was used for pooling purpose. Additionally, I2 statistics was employed to quantify inter-study variability and larger values suggested an increasing degree of heterogeneity [30]. Hardy-Weinberg equilibrium (HWE) in the controls was calculated via chi-square test. Funnel plot asymmetry was estimated by Egger's linear regression test which is a type of linear regression methodology to measure the funnel plot asymmetry on the natural logarithm scale of the OR. The significance of the intercept was determined by the t-test considering p-value<0.05 as representation of statistically significant publication bias [31]. Also, the subgroup analysis was carried out by the ethnicity, and the ethnicity was defined mainly as Caucasians. A comparative examination of ‘meta-analysis’ softwares was performed by using url address http://www.meta-analysis.com/pages/comparisons.html. The Comprehensive Meta-Analysis (CMA) Version 2 software program (Biostat, USA) was chosen and utilized to perform all statistical analysis involved in this study.

Results

Characteristics of included studies

According to our selection (inclusion-exclusion) criteria, a total of twelve research articles were finally included through literature search from the PubMed (Medline), EMBASE and Google Scholar web databases in this meta-analysis. All retrieved research publications were examined carefully by reading the titles and abstracts, and the full texts for the potentially relevant research articles were further checked for their aptness for the current meta-analysis. Studies either showing p53 codon 72 Arg>Pro polymorphism to predict survival in OC patients or considering p53 variants as an indicators for response to therapy were excluded straightaway. Similarly, studies investigating the levels of p53 mRNA or protein expression or relevant review articles were also excluded. In the present meta-analysis, only case-control or cohort design studies having frequency of all three genotypes were included. Besides the database search, the reference lists present in the retrieved articles were also checked for other potential research publications (Table 1). Distribution of genotypes, minor allele frequency (MAF) and HWE in the controls and cases have been presented in Table 2.

Table 1

Characteristics of studies included in this meta-analysis.

First Authors and year	Country of origin	Study Design	Cases	Controls	Source of genotyping
Dholariya et al. 2013 [15]	India	HB	100	100	Blood
Malisic et al. 2013 [16]	Serbia	HB	47	70	Tissue
Matei et al. 2012 [17]	Romania	HB	21	21	Blood
Ueda et al. 2006 [18]	Japan	HB	68	95	Blood
Morari et al. 2006 [19]	Brazil	HB	69	222	Blood
Santos et al. 2006 [20]	Portugal	HB	99	188	Blood
Agorastos et al. 2004 [21]	Greece	HB	51	30	Cytobrush
Pegoraro et al. 2003 [22]	South Africa	HB	85	340	Blood
Hogdall et al. 2002 [23]	Denmark	HB	211	83	Blood
Li et al. 2002 [24]	China	HB	39	50	Tissue
Buller et al. 1997 [25]	America	HB	190	52	Blood
Peller et al. 1999 [26]	Israel	HB	13	13	Blood

Table 2

Distribution of p53 polymorphism of twelve studies included in the meta-analysis.

Authors and year	Control				Case				HWE
	Genotype			Minor allele	Genotype			Minor allele
	Arg/Arg	Arg/Pro	Pro/Pro	MAF	Arg/Arg	Arg/Pro	Pro/Pro	MAF	HWEF
Dholariya et al. 2013	62	32	6	0.22	33	50	17	0.42	0.49
Malisic et al. 2013	45	22	3	0.2	22	22	3	0.29	0.88
Matei et al. 2012	7	7	7	0.5	9	6	6	0.42	0.12
Ueda et al. 2006	34	54	7	0.35	21	41	6	0.38	0.02
Morari et al. 2006	117	91	14	0.26	23	46	0	0.33	0.51
Santos et al. 2006	117	58	13	0.22	49	40	10	0.30	0.12
Agorastos et al. 2004	6	19	25	0.69	26	22	3	0.27	0.42
Pegoraro et al. 2003	32	147	161	0.68	14	41	30	0.59	0.85
Hogdall et al. 2002	48	27	8	0.25	118	73	20	0.26	0.16
Li et al. 2002	10	26	14	0.54	14	20	5	0.38	0.74
Buller et al. 1997	30	18	4	0.25	98	79	13	0.27	0.57
Peller et al. 1999	8	5	0	0.19	7	6	0	0.23	0.39

Publication bias

Begg's funnel plot and Egger's test were carried out to evaluate the publication bias among the selected studies for the meta-analysis. The appearance of the shape of funnel plots was seemed to be symmetrical in all the genetic models. The Egger's test was performed to provide the statistical evidence of funnel plot. The results showed lack of publication bias among all comparison models (Table 3).

Table 3

Statistics to test publication bias and heterogeneity in the present meta-analysis.

Comparisons	Egger's regression analysis			Heterogeneity analysis			Model used for the meta-analysis
	Intercept	95% Confidence Interval	p-value	Q-value	Pheterogeneity	I2 (%)
Pro vs. Arg	−1.56	−7.40 to 4.26	0.56	70.75	<0.0001	84.45	Random
Pro/Pro vs. Arg/Arg	−1.61	−6.48 to 3.26	0.47	45.35	<0.0001	77.95	Random
Arg/Pro vs. Arg/Arg	−2.72	−5.95 to 0.51	0.09	30.12	0.002	63.48	Random
Pro/Pro+Arg/Pro vs. Arg/Arg	−3.38	−7.58 to 0.82	0.10	48.78	<0.0001	77.45	Random
Pro/Pro vs. Arg/Arg+Arg/Pro	−0.49	−3.70 to 2.72	0.73	31.08	0.001	67.83	Random

Evaluation of heterogeneity

In order to analyze heterogeneity among the selected studies, Q-test and I2 statistics were employed and heterogeneity was noticed in all the five genetic models. Therefore, random effects model was applied to synthesize the data (Table 3).

Association of p53 codon 72 Arg>Pro polymorphism and OC susceptibility

We pooled all the twelve studies together and it resulted into 1264 controls and 993 OC cases, to review the overall association between p53 codon 72 Arg>Pro polymorphism and OC risk. Overall pooled analysis did not suggest any correlation between p53 codon 72 Arg>Pro polymorphism and OC risk in all the five genetic comparison models, i.e., allele (Pro vs. Arg: p = 0.916; OR = 0.980, 95% CI = 0.677 to 1.419), homozygous (Pro.Pro vs. Arg.Arg: p = 0.419; OR = 0.731, 95% CI = 0.341 to 1.564), heterozygous (Arg.Pro vs. Arg.Arg: p = 0.248; OR = 1.237, 95% CI = 0.862 to 1.773), recessive (Pro.Pro vs. Arg.Arg+Arg.Pro: p = 0.329; OR = 0.754, 95% CI = 0.428 to 1.329) and dominant model (Pro.Pro+Arg.Pro vs. Arg.Arg: p = 0.699; OR = 1.089, 95% CI = 0.706 to 1.681) (Figure 1).

Figure 1

Forest plot of OR with 95% CI of ovary cancer associated with the p53 codon 72 Arg>Pro gene polymorphism.

Black square represent the value of OR and the size of the square indicates the inverse proportion relative to its variance. Horizontal line is the 95% CI of OR. The studies are listed by year of publication.

Subgroup analysis of racial descent

We have analyzed only Caucasian population by study design and participants. This meta-analysis included eight studies (717 cases and 797 controls), Heterogeneity was observed in all genetic models; thus, random effect model was applied to analyze the data. In subgroup analysis publication bias did not exist (Table 4). We did not observe any association of p53 codon 72 Arg>Pro polymorphism with OC risk in Caucasian population in all genetic comparison models, i.e., allele (Pro vs. Arg: p = 0.542; OR = 0.862, 95% CI = 0.535 to 1.388), homozygous (Pro.Pro vs. Arg.Arg: p = 0.307; OR = 0.630, 95% CI = 0.260 to 1.527), heterozygous (Arg.Pro vs. Arg.Arg: p = 0.781; OR = 1.059, 95% CI = 0.707 to 1.586), recessive (Pro.Pro vs. Arg.Arg+Arg.Pro: p = 0.259; OR = 0.686, 95% CI = 0.356 to 1.320) and dominant model (Pro.Pro+Arg.Pro vs. Arg.Arg: p = 0.710; OR = 0.905, 95% CI = 0.534 to 1.532) (Figure 2)

Table 4

Statistics to test publication bias and heterogeneity in the present meta-analysis.

Comparisons	Egger's regression analysis			Heterogeneity analysis			Model used for the meta-analysis
	Intercept	95% Confidence Interval	p-value	Q value	Pheterogeneity	I2 (%)
Pro vs. Arg	−0.65	−7.59 to 6.28	0.82	44.95	<0.0001	84.43	Random
Pro/Pro vs. Arg/Arg	−1.29	−9.14 to 6.54	0.68	26.87	<0.0001	77.67	Random
Arg/Pro vs. Arg/Arg	−1.69	−5.29 to 1.91	0.29	14.34	0.04	51.19	Random
Pro/Pro+Arg/Pro vs. Arg/Arg	−2.34	−7.37 to 2.68	0.29	27.67	<0.0001	74.70	Random
Pro/Pro vs. Arg/Arg+Arg/Pro	−0.17	−4.97 to 4.63	0.93	18.46	0.005	67.51	Random

Figure 2

Forest plot of OR with 95% CI of ovary cancer associated with the p53 codon 72 Arg>Pro gene polymorphism according to the ethnicity group by the random effect model.

Black square represent the value of OR and the size of the square indicates the inverse proportion relative to its variance. Horizontal line is the 95% CI of OR. The studies are listed by year of publication.

Discussion

The p53 tumor suppressor gene, mutated in many human cancer types, suggests its key role in the host's defense against malignancy [6]. The cellular level of the p53 protein is regulated in a complex manner by a negative feedback loop involving ubiquitin medicated degradation [32]. Based upon the nature of the genetic insult, wild-type p53 induces either growth arrest or apoptosis. Any alterations within p53 gene restrict these activities and allow the continuous proliferation of cells, ultimately resulting in the progression to malignancy [33].

Despite remarkable progress in mechanistic understanding of p53 structure and function, the contribution of specific p53 gene polymorphisms to OC risk remains equivocal and proven to be extremely complex biomarkers. Interest in the genetic susceptibility to OC has led to an emerging trend to the study of polymorphisms of genes involved in OC risk. Due to the different roles of p53 gene in human genome, it has been hypothesized that codon 72 Arg>Pro polymorphism is associated with risk of OC. As a result, a large number of studies have been performed to evaluate the association between p53 codon 72 Arg>Pro polymorphism and risk of OC, but the results from different published studies lacks consensus. Inconsistency in results from these studies can be attributed to low statistical power to evaluate the overall effect of the p53 codon 72 Arg>Pro polymorphism with OC risk. The answer to this limitation is a meta-analysis, a powerful tool for investigating the risk factors associated with genetic diseases, which employs quantitative technique to pool the data from individual studies where individual sample sizes are small with lower statistical power, and provides reliable conclusion [34]. Hence, we have done the present meta-analysis from twelve eligible published case-control studies to evaluate the said relation of p53 codon 72 Arg>Pro polymorphism and risk of OC. This study might help to explore a more robust estimate about the role of this polymorphism with OC risk, as combining data from many studies has the advantage of reduced random errors [35].

The overall pooled results of this meta-analysis revealed that p53 codon 72 Arg>Pro polymorphism did not influence an increased or decreased risk of OC in all the five genetic models as per the eligible studies when compared with wild type allele. Even in the stratified analysis by ethnicity, no statistically significant relationship between p53 codon 72 Arg>Pro genotype and OC risk was detected in Caucasian population. One possible explanation is that several other SNPs have been reported in the p53 gene and previous research showed that various other SNPs are related to the susceptibility to OC. It is possible that the analyzed variant does not act as primary susceptibility polymorphism and may be inhibiting p53 function by linking with other functional polymorphism alleles found in linkage disequilibrium (LD). Our results are in agreement with Shen et al. [36], whereas Zhang et al. 2008 has reported a decreased risk [37]. Moreover, pooled and meta-analyses for breast [38], lung [39] and endometrial [40] cancers do not support a significant role for this polymorphism in susceptibility. Susceptibility to OC is a multistep process in which environmental and genetic factors interact closely and a single genetic variant is usually insufficient to predict the risk of this deadly disease.

Heterogeneity between studies is very common in the genetic association studies of meta-analysis. In the present meta-analysis we found inter-study heterogeneity in overall analysis. There are several factors responsible for such heterogeneity, i.e., the genetic backgrounds for cases and controls, diverse genotype distribution of codon 72 Arg>Pro in different ethnic groups and uneven selection criteria for the cases and controls in different studies.

Despite the important findings from our current analysis, we still have to acknowledge some limitations of this study. First, we only included studies published in English language, abstracted and indexed by the selected electronic databases for the data analysis; it is possible that some pertinent reports published in other languages and indexed in other electronic databases may have missed. Second, the result of this meta-analysis was based on unadjusted ORs because not all eligible studies stated adjusted ORs. Third, the role of gene-environment interactions were not considered which may affect the risk of OC. Also, it is worthwhile to mention several strengths of our study. First, we have included significantly more number of cases and controls comparison to the previous meta-analysis study by using effective and efficient search strategy to increase the statistical power of the analysis. Second, the quality of the case-control studies included in the present pooling analysis was satisfactory and met with the pre-set inclusion criteria.

Conclusion

In conclusion, a meta-analysis is a rational approach of data-analysis which pools both statistically significant and non-significant findings from individual studies to improve the statistical performance by increasing the sample size. Our meta-analysis demonstrates that p53 codon 72 Arg>Pro polymorphism might not significantly modulate the OC risk. However, future well designed large studies, particularly stratified by gene-gene and gene-environment interactions might be necessary to clarify the possible role of the p53 codon 72 Arg>Pro polymorphism in the susceptibility to OC.

PRISMA 2009 Flow Diagram. Showing identification and selection of studies for the meta-analysis.

(TIF)

Click here for additional data file.

PRISMA 2009 Checklist.

(DOC)

Click here for additional data file.