J Cancer 2018; 9(2):269-274. doi:10.7150/jca.21774
A single nucleotide polymorphism in CYP1B1 leads to differential prostate cancer risk and telomere length
1. Department of Urology, Fudan University Shanghai Cancer Center, Shanghai, China;
2. Department of Oncology, Shanghai Medical College, Fudan University, Shanghai, China.
* These authors contributed equally to this work.
Gu CY, Li GX, Zhu Y, Xu H, Zhu Y, Qin XJ, Bo D, Ye DW. A single nucleotide polymorphism in CYP1B1 leads to differential prostate cancer risk and telomere length. J Cancer 2018; 9(2):269-274. doi:10.7150/jca.21774. Available from http://www.jcancer.org/v09p0269.htm
BACKGROUND: Cytochrome P450 1B1 (CYP1B1) is a key enzyme in its oestrogen metabolism pathway, giving rise to hydroxylation and conjugation. Functionally relevant genetic variants within CYP1B1 may affect the telomere length and subsequently lead to prostate carcinogenesis.
METHODS: We evaluated 8 CYP1B1 tag single nucleotide polymorphisms (SNPs) in 1015 men with prostate cancer (PCa) and 1052 cancer-free controls, and calculated odds ratios (ORs) and 95% confidence intervals (CIs) to estimate their association with risk of PCa. The influence of CYP1B1 SNPs on the relative telomere lengths was then appraised in peripheral blood leukocytes using real-time PCR.
RESULTS: CYP1B1 rs1056836 variant was associated with decreased risk of PCa [odds ratio (OR): 0.80; 95% confidence interval (CI): 0.68-0.99, P = 0.041]. Longer telomere length showed a significantly higher proportion of the CYP1B1 rs1056836 CG/GG genotypes, compared with that of the CC genotype (OR: 1.60, 95% CI: 1.04-2.45).
CONCLUSION: Our findings suggest that genetic variants within CYP1B1 may confer genetic susceptibility to PCa by altering telomere length.
Keywords: CYP1B1, prostate cancer, single nucleotide polymorphism, susceptibility.
Telomeres are specialised hexanucleotide repeats complexed to proteins. They are located at ends of linear chromosomes and protect their structural integrity. Telomere length reduces with increasing chronological age and with accumulation of age-related comorbidities. Given that telomere shortening contributes to cellular senescence, telomere length represents a cellular marker of biological aging. A large number of studies suggest that telomere length is a marker of cancer risk and short leukocyte telomere length are believed to be causative events in malignant transformation.
The hypothesis that estrogen influences telomere length came from several preclinical, observational and interventional studies.This phenomenon was explained by the finding that estrogen stimulates telomerase via estrogen response element present in the TERT promoter region. Additional protective effect of estrogen on telomere length is thought to be mediated through its reduction in oxidative stress. Despite the well-known antioxidant role of estrogen, estrogen metabolites, on the other hand, may paradoxically lead to increased oxidative stress through the formation of depurinating DNA adducts and release of surplus amount of superoxide radicals. Several single nucleotide polymorphisms (SNPs) within the estrogen metabolic pathway have been described to have functional effects on the corresponding enzymes, thereby compromising conjugation and redox buffering capacity of cells. Candidates include CYP1B1, which catalyze the conversion of estrogens to reactive catechol estrogens. Furthermore, we described that a CYP1B1 SNP rs1056836 was significantly associated with biochemical recurrence after radical prostatectomy and relative CYP1B1 mRNA expression.
Encouraged by such findings we have now tested CYP1B1 SNPs for associations with risk of PCa in a case-control study. In addition, a real-time PCR based telomere assay was performed to investigate whether leukocyte telomere was in relation to CYP1B1 SNPs.
Materials and Methods
Study Design and Population
The study subjects were mostly from previously published case-control study. Briefly, 1115 eligible patients recruited into this study were newly diagnosed and histopathologically confirmed primary prostate adenocarcinoma from Fudan University Shanghai Cancer Center (FUSCC) between January 2005 and January 2012, of whom 1015 (92%) agreed to participate in this study. All cases had received no prior chemotherapy or radiotherapy upon recruitment. The tumor stage was determined according to criteria established by the American Joint Committee on Cancer (AJCC) tumor-node-metastasis (TNM) classification system [AJCC Staging Manual, sixth edition, 2002]. Histopathological grading of the specimens was performed according to the Gleason score system. The clinical information including Gleason score, serum PSA level at diagnosis and disease stage were abstracted from the archival medical records. In addition, 1143 age-(± 5yr) and geographical regions-matched cancer-free ethnic Han Chinese controls were recruited from the Taizhou longitudinal (TZL) study conducted during a similar time period. TZL study was a large prospective cohort initiated to explore the environmental and genetic risk factors for common non-communicable diseases. Individuals with a known test of serum PSA > 4 ng/mL present with or without abnormal digital rectal examination were excluded from the control group and those without response to the study participation were excluded (n = 91). All of the participants were interviewed with a questionnaire after a written informed consent was obtained. Blood samples were collected and processed as a routine practice by the FUSCC Tissue Bank (for cases) and the TZL study (for controls). This study was approved by the Institutional Review Board of FUSCC.
SNP Selection and Genotyping
Tagging SNPs were selected in the CYP1B1 gene using the Tagger algorithm with r2 ≥ 0.8, and minor-allele frequencies > 0.05 based on the HapMap population data for Han Chinese in Beijing. There were 9 tagging SNPs (rs1056836, rs10916, rs162562, rs2551188, rs9341250, rs9341266, rs9341248, rs162549, rs1056827).
All these nine SNPs were genotyped by the TaqMan real-time PCR method as described previously[8, 9]. Briefly, DNA isolation was performed by using the Qiagen Blood DNA Mini KIT (Qiagen Inc., Valencia, CA) with the buffy-coat fraction of the blood samples donated by the participants. The results with > 99% call rates and 100% concordance for duplicated specimens were acceptable for further genotyping data analysis. SNP rs9341248 was removed because it did not conform to Hardy-Weinberg equilibrium (HWE) (P < 0.05). Thus, a total of 8 SNPs were included for subsequent analysis.
Measurement of Relative Telomere Length
Relative telomere length measurements were available in 426 patients as described previously. Briefly, telomere length as represented by the telomere repeat copy number to single copy gene (i.e., 36B4) copy number (T/S) ratio was measured using real-time quantitative PCR method on an Applied Biosystems 7900HT. The PCR reaction mixture consisted of SYBR Green Mastermix (Applied Biosystems), 100 nmol/L Tel-1, 900 nmol/L Tel-2, 400 nmol/L 36B4d, 400 nmol/L 36B4u, and 7ng of genomic DNA. The thermal cycling profile was 95°C for 10 min followed by 40 cycles of 95°C for 15 sec and at 56°C (for telomere) or 58°C (for 36B4) for 1 min. Following amplification, a dissociation curve confirmed the specificity of the reaction. During each run, negative and positive controls, a calibrator DNA sample, and a standard curve were included. For each standard curve, 2-fold serial dilutions of a reference DNA sample were used to produce a standard curve in each reaction. Two main steps were involved in telomere length quantification: first, the T/S ratio was determined for each sample based on the standard curve. Second, the ratio for each sample was normalized to the calibrator DNA to standardize sample values across all reaction plates. The laboratory personnel were blinded to disease status. R2 for each standard curve was ≥ 0.99.
For all subjects, the χ2 test was used to assess differences in the frequency distributions of the selected demographic variables and genotypes of 8 SNPs between the cases and controls. The HWE for genotype distribution in controls was tested by a goodness-of-fit χ2 test. Odds ratios (ORs) and 95% confidence intervals (CIs) were calculated by univariable and multivariable unconditional logistic regression models to evaluate associations between the genotypes and risk of PCa without and with adjustment for confounding factors, respectively. Further stratification analyses were conducted to calculate the associations of SNP genotypes with PCa risk by demographic and clinicopathologic variables, followed by the homogeneity Q-tests. Spearman rank correlation was used to investigate associations between telomere length and age. Telomere length was categorized into dichotomies, tertile and quartile based on the distribution. To account for chance associations from multiple comparisons, we used the false-positive report probability (FPRP) to assess the false-positive association findings. All statistical analyses were performed with SAS software (version 9.1; SAS Institute, Cary, NC).
Characteristics of the study subjects
The distributions of demographic characteristics of the subjects are shown in Table 1. Briefly, there were no statistical differences in the distributions of age and smoking status between two groups.
Characteristics of PCa Cases and Controls
N = 1015
N = 1052
|Age||69.1 ± 8.2||68.6 ± 8.9||0.828|
|≤ 64||291 (28.7)||308 (29.3)|
|65-75||496 (48.9)||500 (47.5)|
|> 75||228 (22.5)||244 (23.2)|
|Never||406 (40.0)||412 (39.2)||0.697|
|Ever||609 (60.0)||640 (60.8)|
|PSA value (ng/ml)|
|< 10||180 (19.4)|
|> 20||552 (59.6)|
|≤ 7 (3+4)||317 (31.2)|
|≥ 7 (4+3)||606 (59.7)|
|Stage of disease|
PSA: prostate-specifc antigen
Associations of CYP1B1 SNPs with PCa Risk
The genotype frequencies of 8 CYP1B1 SNPs and their associations with PCa risk are summarized in Table 2. The observed genotype frequencies of 8 SNPs in controls agreed with the Hardy-Weinberg equilibrium. Furthermore, the genotype distribution of CYP1B1 rs1056836 was significantly different between the cases and controls. . In multivariate logistic regression analysis, CYP1B1 rs1056836 was associated with PCa risk. Compared with CC genotype, rs1056836 variant CG/GG genotype was associated with a decreased risk of PCa (OR: 0.80, 95% CI = 0.66-0.97). However, no associations between the other 7 SNPs and PCa risk were observed.
Association of CYP1B1 SNPs with PCa Risk
|SNP||Cases, No. (%)|
n = 1015
|Controls, No. (%)|
n = 1052
|OR (95% CI)||Pa|
|GG||774 (76.3)||816 (77.6)||1.00|
|AG||218(21.5)||215 (20.4)||1.07 (0.87-1.32)||0.533|
|AA||23 (2.3)||21 (2.0)||1.15 (0.63-2.10)||0.649|
|AG/AA vs. GG||1.08 (0.88-1.32)||0.477|
|GG/AG vs. AA||1.13 (0.62-2.06)||0.683|
|TT||787 (77.5)||811 (77.1)||1.00|
|AT||209 (20.6)||222 (21.1)||0.95 (0.77-1.18)||0.636|
|AA||19 (1.9)||19 (1.8)||1.01 (0.53-1.93)||0.971|
|AT/AA vs. TT||0.96 (0.78-1.17)||0.660|
|TT/AT vs. AA||1.02 (0.54-1.95)||0.944|
|TT||808 (79.6)||856 (81.4)||1.00|
|GT||192 (18.9)||183 (17.4)||1.12 (0.89-1.40)||0.333|
|GG||15 (1.5)||13 (1.2)||0.75 (0.50-1.11)||0.152|
|GT/GG vs. TT||1.13 (0.91-1.40)||0.289|
|TT/GT vs. GG||1.21 (0.57-2.56)||0.617|
|AA||695 (68.5)||711 (67.6)||1.00|
|AC||287 (28.3)||305 (29.0)||0.97 (0.80-1.17)||0.719|
|CC||33 (3.3)||36 (3.4)||0.96 (0.59-1.55)||0.854|
|AC/CC vs. AA||0.96 (0.80-1.16)||0.700|
|CC/AC vs.CC||0.97 (0.60-1.56)||0.887|
|CC||541 (53.3)||558 (53.0)||1.00|
|CT||423 (41.7)||421 (40.0)||1.05 (0.88-1.26)||0.564|
|TT||51 (5.0)||73 (6.9)||0.72 (0.49-1.05)||0.086|
|CT/TT vs. CC||1.00 (0.84-1.19)||0.964|
|CC/CT vs.TT||0.70 (0.49-1.02)||0.060|
|CC||955 (94.1)||985 (93.6)||1.00|
|AC||58 (5.7)||64 (6.1)||0.96 (0.66-1.38)||0.806|
|AA||2 (0.2)||3 (0.3)||0.66 (0.11-3.94)||0.644|
|AC/AA vs. CC||0.94 (0.66-1.35)||0.741|
|CC/AC vs. AA||0.66 (0.11-3.95)||0.646|
|CC||668 (65.8)||679 (64.5)||1.00|
|AC||303 (29.9)||327 (31.1)||0.94 (0.78-1.14)||0.543|
|AA||44 (4.3)||46 (4.4)||0.98 (0.64-1.50)||0.908|
|AC/AA vs. CC||0.95 (0.79-1.14)||0.555|
|CC/AC vs. AA||0.99 (0.65-1.52)||0.976|
|CC||728 (71.7)||705 (67.0)||1.00|
|CG||274 (27.0)||323 (30.7)||0.82 (0.68-0.99)||0.041|
|GG||13 (1.3)||24 (2.3)||0.54 (0.27-1.07)||0.079|
|CG/GG vs. CC||0.80 (0.66-0.97)||0.020|
|CC/CG vs. GG||0.57 (0.29-1.14)||0.111|
HWE: Hardy-Weinberg equilibrium. aAdjusted for age and smoking status.
The results were in bold if P < 0.05.
Association between CYP1B1 SNPs and RTL
|SNP||Genotype||Short RTL, n (%)||Long RTL, n (%)||OR (95% CI)a||P|
|rs9341266||GG||58 (27.9)||64 (29.4)||1.00|
|AG/AA||150 (72.1)||154 (70.6)||1.08 (0.71-1.64)||0.737|
|rs162549||TT||44 (21.2)||48 (22.0)||1.00|
|AT/AA||164 (78.8)||170 (78.0)||1.05 (0.66-1.67)||0.828|
|rs10916||TT||35 (16.8)||54 (24.8)||1.00|
|GT/GG||173 (83.2)||164 (75.2)||1.63 (1.01-2.62)||0.045|
|rs162562||AA||60 (28.8)||66 (30.3)||1.00|
|AC/CC||148 (71.2)||152 (69.7)||1.07 (0.71-1.62)||0.747|
|rs2551188||CC||82 (39.4)||103 (47.2)||1.00|
|CT/TT||126 (60.6)||115 (52.8)||1.38 (0.94-2.02)||0.104|
|rs9341250||CC||13 (6.2)||17 (7.8)||1.00|
|AC/AA||195 (93.8)||201 (92.2)||1.27 (0.60-2.68)||0.533|
|rs1056827||CC||67 (32.2)||69 (31.7)||1.00|
|AC/AA||141 (67.8)||149 (68.3)||0.98 (0.65-1.47)||0.901|
|rs1056836||CC||67 (32.2)||50 (22.9)|
|CG/GG||141 (67.8)||168 (77.1)||1.60 (1.04-2.45)||0.033|
RTL: Relative telomere length aAdjusted for age and smoking status.
In the stratified analysis, the decreased risk was also observed in men who were ever smokers, Gleason score ≥7 (4 + 3), and stage III/IV (Supplementary Table 1). Further homogeneity tests suggested, however, that there were no differences in the risk estimates between these strata.
Genetic variations and telomere length
To investigate whether the telomere length had prognostic significance, 426 patients were categorized into dichotomies, tertile and quartile based on their telomere length distribution (range, 0.06-2.06) (Supplementary Table 2). Telomere length was inversely associated with age (r = -0.40, P < 0.001). When participants were dichotomized according to the median telomere length value, significant differences in telomere length by genotype of CYP1B1 rs1056836 and rs10916 were observed (Table 3). Longer telomere length showed a significantly higher proportion of the CG/GG genotypes, compared with that of the CC genotype (OR: 1.60, 95% CI: 1.04-2.45).
FPRP Values for Associations between PCa Risk and CYP1B1 SNPs
The FPRP values at different prior probability levels for all significant findings are summarized in the Supplementary Table 3. When the assumption of prior probability was 0.01, no association was noteworthy in all subgroups for rs1056836 as well as rs10916. In contrast, some greater FPRP values for the other significant associations between CYP1B1 variants and PCa risk suggested some possible false-positive results in the findings, which need further validation in larger studies.
In the current study we found a CYP1B1 inherited variation rs1056836 associated with PCa risk and further showed its biologic association with leukocyte telomere length. Findings on CYP1B1 polymorphisms and susceptibility to PCa have proved inconsistent across different ethnic groups[8-14]. However, previous studies were conducted in small-sized populations, and require confirmation in larger studies.
The SNP rs1056836 is a non-synonymous polymorphism in exon 3, which encodes the heme binding domain, and results in a G-to-C and subsequent amino acid substitutions of Valine-to-Leucine.The SNP has documented functional consequences on transcription, protein structure, and enzyme action[16, 17]. A number of studies have reported a decreased risk of PCa for carriers of the rs1056836 G allele, but the results were inconsistent in different ethnic groups. Meta-analysis results showed that rs1056836 is significantly associated with PCa susceptibility among Asians[18, 19] and our findings support this association. These inconclusive results may have several causes. First, it may result from the distribution of allele frequencies among different ethnicities or between different case sources. We found that the frequency of rs1056836 G allele was significantly lower in Asian controls than in Caucasian or mixed population controls. Second, the sample sizes of the cases and controls in the Asian populations were both noticeably smaller than those in Caucasians or mixed population subjects[18, 19]. The current study has sufficient power confirm the positive associations in Asians. Moreover, the mechanism of PCa development is similar to other cancers and is dependent on the interactions of genetic factors and environmental agents. Environments and lifestyles are very different among individuals of different races. Thus, the aetiology of PCa is unlikely to be explained solely by genetic polymorphisms because hereditary variation alone cannot affect the risk of any disease. CYP1B1 rs1056836 was associated with CYP1B1 mRNA expression and overlapped with PCa risk in the expected opposite direction suggesting potential biological consequences. Furthermore, a recent study showed that overexpression of CYP1B1 stimulated proliferative, migratory and invasive potential of non-tumorigenic PCa cells.
However, exact mechanisms of how CYP1B1 polymorphism contributes to PCa susceptibility requires further illustration. The present dataset offers an attractive possibility to address interactions between telomere length and genotypes. We hypothesized that genetic factors implicated in estrogen metabolism pathways contribute to increased oxidative stress and thereby influence telomere length. Oxidative stress causes telomere attrition by several potential pathways. Reactive oxygen species can directly damage the sensitive triple-G containing telomeres, or indirectly, act by inhibiting telomerase activity and increasing the load of pro-inflammatory cytokines that accelerate immune cell turnover and hence telomere loss. Since estrogen metabolism is a well-known source of reactive oxygen species formation, aberrations in this pathway could compromise the redox buffering capacity of cells, which may in turn lead to increased oxidative burden and telomere attrition. Accordingly, in our study, we show for the first time that functionally relevant SNP in CYP1B1 gene impacts telomere length in PCa patients. CYP1B1 participates in the first step of estrogen metabolism, i.e., the conversion of estrogens to 2- or 4-hydroxyestrogens, and specifically catalyzes the 4-hydroxylation of estrogens. 4-hydroxyestradiol is, in sharp contrast to 2-hydroxyestradiol, inactivated by COMT at a slower rate, and undergoes metabolic redox cycling to generate free radicals and reactive semiquinone and quionone intermediates. The CYP1B1 1056836 variant allele codes for the enzyme with 2-fold higher mRNA compared to wild-type allele, leading to increased 4-hydroxyestradiol production. Consistent with this notion, telomere length was found significantly shorter in individuals carrying the CYP1B1 rs1056836 wild type allele.
A large number of intrinsic and extrinsic factors , such as heredity, epigenetics, aging, stress, immune components, and hormones are all states have been previously associated with shorter telomere length. As such, limitations of the study are mainly related to the scarcity of the relevant exposure data and nature of retrospective study design. The current study was based on the leukocytes telomere length only and not on the prostate tissue telomere length. Because there is an intra-individual synchrony in telomere length across the somatic tissues of humans as evidenced by the strong correlations between the telomere lengths in all tissue types. Moreover, the rates of telomere shortening are similar in the somatic tissues. However, future work aims to directly describe the relationships between prostate tissue, leukocytes telomere lengths and PCa cases. Strengths of the study include a significant number of patients, a candidate gene approach, the high plausibility of the association based on the biologic function of selected candidate genes and an impact on telomere length associated with positive markers.
In summary, significant association with PCa risk for CYP1B1 rs1056836 with telomere length may provide a biological basis for the results obtained in various studies investigating the impact of CYP1B1 SNPs on PCa risk. With these implications, further investigation of SNP‑cancer risk associations in specific subsets of patients with particular environmental or genetic backgrounds, as well as more detailed in vitro and in vivo biological functional studies will be helpful to elucidate how exactly the CYP1B1 genetic variations influenced PCa development.
PSA: prostate-specific antigen; SNP: single nucleotide polymorphism; OR: odds ratio; CI: confidence interval.
The authors have declared that no competing interest exists.
1. Blackburn E H, Greider C W, Szostak J W. Telomeres and telomerase: the path from maize, Tetrahymena and yeast to human cancer and aging. Nat Med. 2006;12(10):1133-1138
2. Benetos A, Dalgard C, Labat C. et al. Sex difference in leukocyte telomere length is ablated in opposite-sex co-twins. Int J Epidemiol. 2014;43(6):1799-1805
3. Song J Y, Kim M J, Jo H H. et al. Antioxidant effect of estrogen on bovine aortic endothelial cells. J Steroid Biochem Mol Biol. 2009;117(1-3):74-80
4. Yue W, Santen R J, Wang J P. et al. Genotoxic metabolites of estradiol in breast: potential mechanism of estradiol induced carcinogenesis. J Steroid Biochem Mol Biol. 2003;86(3-5):477-486
5. Gajjar K, Martin-Hirsch P L, Martin F L. CYP1B1 and hormone-induced cancer. Cancer Lett. 2012;324(1):13-30
6. Gu C Y, Qin X J, Qu Y Y. et al. Genetic variants of the CYP1B1 gene as predictors of biochemical recurrence after radical prostatectomy in localized prostate cancer patients. Medicine (Baltimore). 2016;95(27):e4066
7. Gu C, Qu Y, Zhang G. et al. A single nucleotide polymorphism in ADIPOQ predicts biochemical recurrence after radical prostatectomy in localized prostate cancer. Oncotarget. 2015;6(31):32205-32211
8. Cussenot O, Azzouzi A R, Nicolaiew N. et al. Combination of polymorphisms from genes related to estrogen metabolism and risk of prostate cancers: the hidden face of estrogens. J Clin Oncol. 2007;25(24):3596-3602
9. Beuten J, Gelfond J A, Byrne J J. et al. CYP1B1 variants are associated with prostate cancer in non-Hispanic and Hispanic Caucasians. Carcinogenesis. 2008;29(9):1751-1757
10. Sobti R C, Onsory K, Al-Badran A I. et al. CYP17, SRD5A2, CYP1B1, and CYP2D6 gene polymorphisms with prostate cancer risk in North Indian population. DNA Cell Biol. 2006;25(5):287-294
11. Fukatsu T, Hirokawa Y, Araki T. et al. Genetic polymorphisms of hormone-related genes and prostate cancer risk in the Japanese population. Anticancer Res. 2004;24(4):2431-2437
12. Chang B L, Zheng S L, Isaacs S D. et al. Polymorphisms in the CYP1B1 gene are associated with increased risk of prostate cancer. Br J Cancer. 2003;89(8):1524-1529
13. Tanaka Y, Sasaki M, Kaneuchi M. et al. Polymorphisms of the CYP1B1 gene have higher risk for prostate cancer. Biochem Biophys Res Commun. 2002;296(4):820-826
14. Cicek M S, Liu X, Casey G. et al. Role of androgen metabolism genes CYP1B1, PSA/KLK3, and CYP11alpha in prostate cancer risk and aggressiveness. Cancer Epidemiol Biomarkers Prev. 2005;14(9):2173-2177
15. Holt S K, Kwon E M, Fu R. et al. Association of variants in estrogen-related pathway genes with prostate cancer risk. Prostate. 2013;73(1):1-10
16. Bieche I, Narjoz C, Asselah T. et al. Reverse transcriptase-PCR quantification of mRNA levels from cytochrome (CYP)1, CYP2 and CYP3 families in 22 different human tissues. Pharmacogenet Genomics. 2007;17(9):731-742
17. Beuten J, Gelfond J A, Franke J L. et al. Single and multigenic analysis of the association between variants in 12 steroid hormone metabolism genes and risk of prostate cancer. Cancer Epidemiol Biomarkers Prev. 2009;18(6):1869-1880
18. Yang J, Xu D L, Lu Q. et al. Prostate cancer risk and aggressiveness associated with the CYP1B1 4326C/G (Leu432Val) polymorphism: a meta-analysis of 2788 cases and 2968 controls. Asian J Androl. 2012;14(4):560-565
19. Zhang H, Li L, Xu Y. CYP1B1 polymorphisms and susceptibility to prostate cancer: a meta-analysis. PLoS One. 2013;8(7):e68634
20. Pfeifer G P. Environmental exposures and mutational patterns of cancer genomes. Genome Medicine. 2010;2(8):54
21. Chang I, Mitsui Y, Kim S K. et al. Cytochrome P450 1B1 inhibition suppresses tumorigenicity of prostate cancer via caspase-1 activation. Oncotarget. 2017;8(24):39087-39100
22. Kawanishi S, Oikawa S. Mechanism of telomere shortening by oxidative stress. Ann N Y Acad Sci. 2004;1019:278-284
23. Haendeler J, Hoffmann J, Brandes R P. et al. Hydrogen peroxide triggers nuclear export of telomerase reverse transcriptase via Src kinase family-dependent phosphorylation of tyrosine 707. Mol Cell Biol. 2003;23(13):4598-4610
24. Tsuchiya Y, Nakajima M, Yokoi T. Cytochrome P450-mediated metabolism of estrogens and its regulation in human. Cancer Lett. 2005;227(2):115-124
25. Kaszubowska L. Telomere shortening and ageing of the immune system. J Physiol Pharmacol. 2008;59(Suppl 9):169-186
26. Daniali L, Benetos A, Susser E. et al. Telomeres shorten at equivalent rates in somatic tissues of adults. Nat Commun. 2013;4:1597
Corresponding authors: Bo Dai (bodai1978com) and Ding-Wei Ye (dwyelicom), Department of Urology, Fudan University Shanghai Cancer Center, No. 270 Dong'an Road, Shanghai 200032, People's Republic of China. Tel: 86-21-64175590-2807; Fax: 86-21-64434556.