Abstract
This meta-analysis aims at evaluating the relationships between CYP1A1 genetic polymorphisms and bladder cancer risk. The PubMed, CISCOM, CINAHL, Web of Science, Google Scholar, EBSCO, Cochrane Library, and CBM databases were searched from inception through November 1st, 2013 without language restrictions. Meta-analysis was conducted with the use of the STATA 12.0 software. The relationships were evaluated by calculating the pooled odds ratios (ORs) and their 95% confidence intervals (CIs). Eight case-control studies with a total of 2120 bladder cancer patients and 2061 healthy subjects met the inclusion criteria. Ten common polymorphisms in the CYP1A1 gene were assessed. The results of our meta-analysis suggested that CYP1A1 genetic polymorphisms might be strongly correlated with an increased risk of bladder cancer (allele model: OR=1.23, 95%CI=1.08–1.39, p=0.001; dominant model: OR=1.25, 95%CI=1.07–1.46, p=0.005; respectively), especially for 11599G>C, 2455A>G, 3810T>C, and 113T>C polymorphisms. A subgroup analysis was done to investigate the effect of ethnicity on an individual's risk of bladder cancer. Our results revealed positive significant correlations between CYP1A1 genetic polymorphisms and an increased risk of bladder cancer among Asians (allele model: OR=1.33, 95%CI=1.08–1.65, p=0.009; dominant model: OR=1.37, 95%CI=1.02–1.85, p=0.034; respectively), but not among Caucasians (all p<0.05). Our findings provide convincing evidence that CYP1A1 genetic polymorphisms may contribute to susceptibility to bladder cancer, especially for 11599G>C, 2455A>G, 3810T>C, and 113T>C polymorphisms among Asians.
Introduction
B
Cytochrome-P450 1A1 (CYP1A1) belongs to the large and important superfamily of cytochrome P450 (CYPs) that is involved in the metabolism of many xenobiotic compounds and endogenous lipophilic substances, which can activate procarcinogens into DNA reactive metabolites (Sridhar et al., 2010; van Waterschoot and Schinkel, 2011). Human CYP1A1 gene is located on chromosome 15q22–24 and composed of seven exons (Durocher et al., 1998). It has been commonly known that CYP1A1 is responsible for the activation of procarcinogens to reactive metabolites and is implicated in mediating oxidation, including polycyclic aromatic hydrocarbons and aromatic amines to mutagenic and carcinogenic metabolites (Beedanagari et al., 2009; Jarukamjorn et al., 2010). However, CYP1A1 genetic polymorphisms may lead to alteration or a complete loss in activity of the corresponding metabolizing enzymes, and result in impaired cellular detoxification, which may cause absence of mechanisms for cancer prevention, thereby contributing to susceptibility to bladder carcinogenesis (Ozturk et al., 2011; Cash et al., 2012).
Many previous studies have demonstrated that CYP1A1 genetic polymorphisms may be implicated in the pathogenesis of bladder cancer (Yang et al., 2007; Grando et al., 2009). Nevertheless, contradictory results were also reported in many of the other studies (Brockmoller et al., 1996; Ozturk et al., 2011; Fu and Chen, 2013). Consequently, we performed the present meta-analysis to evaluate the relationships between CYP1A1 genetic polymorphisms and bladder cancer risk.
Materials and Methods
Search strategy
The PubMed, CISCOM, CINAHL, Web of Science, Google Scholar, EBSCO, Cochrane Library, and CBM databases were searched for relevant articles published before November 1st, 2013 without language restrictions. The following keywords and MeSH terms were used: (“SNP” or “mutation” or “genetic polymorphism” or “variation” or “polymorphism” or “single nucleotide polymorphism” or “variant”) and (“urinary bladder neoplasm” or “urinary bladder cancer” or “bladder neoplasm” or “bladder cancer” or “bladder carcinogenesis”) and (“cytochrome P450 1A1” or “CYP1A1” or “cytochrome P-450 CYP1A1”). We also performed a manual search to find other potential articles.
Selection criteria
The included studies should meet all four of the following criteria: (1) the study design should be a clinical cohort or a case-control study that focused on the relationships between CYP1A1 genetic polymorphisms and bladder cancer risk; (2) all patients diagnosed with bladder cancer should be confirmed through histopathologic examinations; (3) the study should provide sufficient information about genotype frequencies of CYP1A1 genetic polymorphisms; and (4) the genotype frequencies should follow the Hardy–Weinberg equilibrium (HWE). If the study could not meet the inclusion criteria, it would be excluded. The most recent or the largest sample size publication was included when the authors published several studies using the same subjects.
Data extraction
Relevant data were systematically extracted from all included studies by two observers by using a standardized form. The researchers collected the following data: language of publication, publication year of article, the first author's surname, geographical location, design of study, sample size, the source of the subjects, genotype frequencies, source of samples, genotyping method, evidence of HWE, and so on.
Quality assessment
Methodological quality was independently assessed by two researchers according to the Newcastle-Ottawa Scale (NOS) criteria (Stang, 2010). The NOS criteria included three aspects: (1) subject selection: 0–4; (2) comparability of subject: 0–2; (3) clinical outcome: 0–3. NOS scores ranged from 0 to 9 with and a score≥7 indicating a good quality.
Statistical analysis
The STATA version 12.0 (Stata Corp, College Station, TX) software was used for meta-analysis. We calculated crude odds ratios (ORs) with their 95% confidence intervals (95%CIs) to evaluate the relationships between common polymorphisms in the CYP1A1 gene and bladder cancer risk. Genotype frequencies of healthy controls were tested for the HWE using the χ2 test. The statistical significance of pooled ORs was assessed by the Z test. The Cochran's Q-statistic and I2 test were used to evaluate potential heterogeneity between studies (Zintzaras and Ioannidis, 2005). If Q test shows a p<0.05 or I2 test exhibits >50%, which indicates significant heterogeneity, the random-effects model was conducted, or else, the fixed-effects model was used. We also performed subgroup and meta-regression analyses to investigate potential sources of heterogeneity. We conducted a sensitivity analysis to assess the influence of single studies on the overall ORs. Begger's funnel plots and Egger's linear regression test were used to investigate publication bias (Peters et al., 2006).
Results
Characteristics of included studies
Initially, the searched keywords identified 53 articles. We reviewed the titles and abstracts of all articles and excluded 24 articles; full texts and data integrity were also reviewed, and 21 articles were further excluded. Finally, eight case-control studies with a total of 2120 bladder cancer patients and 2061 healthy subjects met our inclusion criteria for qualitative data analysis (Brockmoller et al., 1996; Yang et al., 2007; Srivastava et al., 2008; Grando et al., 2009; Villanueva et al., 2009; Ozturk et al., 2011; Berber et al., 2013; Fu and Chen, 2013). Publication years of the eligible studies ranged from 1996 to 2013. Figure 1 shows the selection process of eligible articles. Distribution of the number of topic-related literatures in the electronic database during the last decade is shown in Figure 2. Overall, five studies were conducted among Asians and three studies were conducted among Caucasians. Polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) method was performed in four studies, two studies used allele specific-PCR (AS-PCR) method, while only one study used the TaqMan assay method, and one study used PCR-confronting two-pair primers (PCR-CTPP). Ten common polymorphisms (4889A>G, 462A>G, 728C>T, 11599G>C, 17961C>T, 2455A>G, 3810T>C, 113T>C, 139A>G, and 6235C>T) in CYP1A1 gene were assessed. Genotype frequencies of controls were all in HWE (all p>0.05). None of the studies deviated from the HWE (all p>0.05). NOS scores of all included studies were ≥5. We summarized the study characteristics and methodological quality in Table 1.

Flowchart of literature search and study selection. Eight case-control studies were included in this meta-analysis.

Distribution of the number of topic-related literature in the electronic database during the last decade.
M, male; F, female; SNP, single nucleotide polymorphism; PCR-RFLP, polymerase chain reaction-restriction fragment length polymorphism; PCR-CTPP, polymerase chain reaction-restriction-confronting two-pair primers; AS-PCR, allele specific-polymerase chain reaction-restriction; HWE, Hardy–Weinberg equilibrium; NOS, the Newcastle-Ottawa Scale.
Quantitative data synthesis
Meta-analysis findings on the relationships of CYP1A1 genetic polymorphisms with the risk of bladder cancer were shown in Table 2. The random-effects model was conducted due to obvious heterogeneity existing between studies. The results of our meta-analysis suggested that CYP1A1 genetic polymorphisms might be strongly correlated with an increased risk of bladder cancer (allele model: OR=1.23, 95%CI=1.08–1.39, p=0.001; dominant model: OR=1.25, 95%CI=1.07–1.46, p=0.005; respectively), especially for 11599G>C, 2455A>G, 3810T>C, and 113T>C polymorphisms (Fig. 3).

Forest plots for the relationships between CYP1A1 genetic polymorphisms and bladder cancer risk under the allele and dominant models.
W, wild allele; M, mutant allele; WW, wild homozygote; WM, heterozygote; MM, mutant homozygote; OR, odds ratio; 95%CI, 95% confidence interval.
A subgroup analysis was done to investigate the effect of ethnicity on an individual's risk of bladder cancer (Fig. 4). Our results revealed positive significant correlations between CYP1A1 genetic polymorphisms and an increased risk of bladder cancer among Asians (allele model: OR=1.33, 95%CI=1.08–1.65, p=0.009; dominant model: OR=1.37, 95%CI=1.02–1.85, p=0.034; respectively), but not among Caucasians (all p<0.05). We also performed a subgroup analysis by genotyping method; the results indicated significant associations between CYP1A1 genetic polymorphisms and an increased risk of bladder cancer in both PCR-RFLP and non-PCR-RFLP subgroups (as shown in Table 2).

Subgroup analyses by ethnicity and genotyping method of the relationships between CYP1A1 genetic polymorphisms and bladder cancer risk under the allele and dominant models.
Meta-regression analyses results confirmed that ethnicity might be a main source of heterogeneity (as shown in Table 3). The results of sensitivity analysis indicated that the overall pooled ORs could not be affected by single study (Fig. 5). No evidence for asymmetry was observed in the Begger's funnel plots (Fig. 6). Egger's test also failed to reveal any evidence of publication bias (allele model: t=1.41, p=0.189; dominant model: t=1.33, p=0.212).

Sensitivity analysis of the summary odds ratio coefficients on the relationships between CYP1A1 genetic polymorphisms and bladder cancer risk under the allele and dominant models.

Begger's funnel plot of publication biases on the relationships between CYP1A1 genetic polymorphisms and bladder cancer risk under the allele and dominant models.
SE, standard error; UL, upper limit; LL, lower limit.
Discussion
CYP1A1, a member of CYP1 family, is a hepatic enzyme that is widely distributed in extra-hepatic tissues such as the lungs, kidneys, gastrointestinal tract, throat, placenta, and lymphocytes (Delescluse et al., 2000; Stejskalova and Pavek, 2011). It should be noted that CYP1A1 serves as a phase I detoxification enzyme which is implicated in the activation of pro-carcinogens and protoxin, and the subsequent generation of strong electrophilic compounds is capable of interacting with intracellular macromolecules such as DNA, RNA, proteins, and other nucleophilic group (Suchocki et al., 2010; Masek et al., 2011). Consequently, the metabolic activation may destroy the cell structure and produce inactivated or abnormal enzymes, and inevitably induce CYP1A1 mutations or inhibit its expression, thereby affecting the metabolic and transport pathways of estrogens and evoking programmed cell death or tumorigenesis (Young et al., 2010; Chen et al., 2012). Therefore, it has been postulated that polymorphic variants conferring elevated CYP1A1 enzyme activity may lead to increased production of catechol estrogens and be conducive to the development of bladder cancer (Fontana et al., 2009; Yousef et al., 2012).
In the present meta-analysis, our findings revealed that bladder cancer patients had a higher frequency of CYP1A1 genetic polymorphisms than healthy controls, suggesting that common polymorphisms in the CYP1A1 gene may be causative factors for the incidence of bladder cancer. Although the exact function of CYP1A1 in the pathogenesis of bladder cancer is not yet poorly understood, a potential explanation might be that CYP1A1 genetic polymorphisms may increase the level of its expression through the aryl hydrocarbon receptor (AhR), which is responsible for the xenobiotic activating pathway of several phase I and phase II metabolizing enzymes and involved in various cell signaling pathways critical to cell cycle regulation and normal homeostasis (Nukaya et al., 2010; Mohebati et al., 2012). Furthermore, dysregulation of these pathways caused by CYP1A1 genetic polymorphisms may be closely linked to the development and progression of bladder cancer (Androutsopoulos et al., 2009). A previous study also demonstrated that the A/G and G/G variants of the CYP1A1 gene may contribute to the development of bladder cancer (Grando et al., 2009). Since heterogeneity obviously existed, we performed stratified analyses based on ethnicity and genotyping method. Our findings suggested that CYP1A1 genetic polymorphisms were associated with an increased risk of bladder cancer among Asians, but not among Caucasians, revealing that ethnic differences may cause heterogeneity between studies concerning the role of CYP1A1 genetic polymorphisms in the pathogenesis of bladder cancer. This may be explained that natural selection and random genetic drift may affect an individual's susceptibility to bladder cancer. Meta-regression analysis also confirmed that ethnicity might be a main source of heterogeneity. In short, the results of our meta-analysis were consistent with previous studies that CYP1A1 genetic polymorphisms may contribute to susceptibility to bladder cancer.
The current meta-analysis also had many limitations that should be acknowledged. First, due to relatively small sample size, our results lacked sufficient statistical power to assess the effects of CYP1A1 genetic polymorphisms on susceptibility to bladder cancer. Second, meta-analysis is a retrospective study that may lead to subject selection bias, thereby affecting the reliability of our results. Third, our meta-analysis failed to obtain original data from the included studies, which may limit a further evaluation of potential roles of CYP1A1 genetic polymorphisms in bladder carcinogenesis. Although our study has many limitations, this is the first meta-analysis that focused on the relationships of CYP1A1 genetic polymorphisms with the risk of bladder cancer. Furthermore, we performed a highly sensitive literature search strategy for electronic databases. A manual search of the reference lists from the relevant articles was also conducted to find other potential articles. The selection process of eligible articles was based on strict inclusion and exclusion criteria. Importantly, a rigorous statistical analysis of SNP data provided a basis for pooling of information from individual studies.
In conclusion, our findings provide convincing evidence that CYP1A1 genetic polymorphisms may contribute to susceptibility to bladder cancer, especially for 11599G>C, 2455A>G, 3810T>C, and 113T>C polymorphisms among Asians. Thus, CYP1A1 genetic polymorphisms may be utilized as potential biomarkers for an early diagnosis of bladder cancer. However, due to the limitations mentioned earlier, more research with a larger sample size is needed to provide a more representative statistical analysis.
Footnotes
Acknowledgments
This study was supported by the Science and Technology Research Project of the Higher Education Department of Liaoning Province (No. L2010695). The authors would like to acknowledge the helpful comments on this article received from their reviewers.
Disclosure Statement
The authors declare that they have no conflicts of interest.
