- Research article
- Open Access
Polymorphisms in the cytochrome P450 genes CYP1A2, CYP1B1, CYP3A4, CYP3A5, CYP11A1, CYP17A1, CYP19A1and colorectal cancer risk
BMC Cancer volume 7, Article number: 123 (2007)
Cytochrome P450 (CYP) enzymes have the potential to affect colorectal cancer (CRC) risk by determining the genotoxic impact of exogenous carcinogens and levels of sex hormones.
To investigate if common variants of CYP1A2, CYP1B1, CYP3A4, CYP3A5, CYP11A1, CYP17A1 and CYP19A1 influence CRC risk we genotyped 2,575 CRC cases and 2,707 controls for 20 single nucleotide polymorphisms (SNPs) that have not previously been shown to have functional consequence within these genes.
There was a suggestion of increased risk, albeit insignificant after correction for multiple testing, of CRC for individuals homozygous for CYP1B1 rs162558 and heterozygous for CYP1A2 rs2069522 (odds ratio [OR] = 1.36, 95% confidence interval [CI]: 1.03–1.80 and OR = 1.34, 95% CI: 1.00–1.79 respectively).
This study provides some support for polymorphic variation in CYP1A2 and CYP1B1 playing a role in CRC susceptibility.
Recent studies have provided strong evidence that exposure to carcinogens, such as polycyclic aromatic hydrocarbons, heterocyclic amines and others in the diet, influences the risk of developing colorectal cancer (CRC) [1–3]. In addition to diet, cigarette smoking is another source of carcinogenic exposure relevant to the large bowel, and increased risk of CRC associated with smoking has been documented in several studies (reviewed in ). The genotoxic impact of carcinogen exposure is heavily influenced by a complex array of metabolic pathways, which includes the cytochrome P450 (CYP) enzyme system. The CYP enzymes also participate in the metabolism of a number of endogenous compounds, such as sex hormones and fatty acids, which are increasingly recognized to be relevant to CRC development . Because wide inter-individual variations in activity of many of the CYP enzymes have been related to the existence of genetic polymorphisms, there is an opportunity to look for inherited metabolic susceptibilities to CRC.
In the present study, we have investigated the relationship between single nucleotide polymorphisms (SNP) in CYP1A2, CYP1B1, CYP3A4, CYP3A5, CYP11A1, CYP17A1 and CYP19A1 and CRC risk in a large study of 2,575 cases and 2,707 controls. Our rationale for analysing SNPs mapping to these CYP genes is that CYP3A4 and CYP3A5 are highly expressed in colonic tissue (reviewed in ) and variants of CYP1A2 and CYP1B1 have previously been associated with CRC risk . Additionally, CYP11A1, CYP17A1 and CYP19A1, which encode sex hormone metabolising CYP enzymes, have been linked to risk of other types of cancer [8, 9] and hence represent credible candidates as CRC predisposition genes.
Patients and control subjects
Patients with histologically confirmed adenocarcinomas, ascertained through an ongoing initiative at the Institute of Cancer Research/Royal Marsden Hospital NHS Trust (RMHNHST) (1,474 males, 1,101 females; mean age at diagnosis 59 years; SD ± 10.1) were included in the study. Healthy individuals were recruited at centers throughout the UK as part of the National Cancer Research Network Trial (1999–2002), the Royal Marsden Hospital Trust/Institute of Cancer Research Family History and DNA Registry (1999–2004), or the National Study of Colorectal Cancer Genetics Trial (2004), all established within the United Kingdom. Controls (836 males, 1,871 females; mean age 59 years; SD ± 10.9) were the spouses or unrelated friends of patients with malignancies. None of the controls had a personal history of malignancy at time of ascertainment. All patients and controls were British Caucasians, and there were no obvious differences in the demography of the two groups in terms of place of residence within the UK. Blood samples were obtained with informed consent and ethical review board approval in accordance with the tenets of the Declaration of Helsinki. DNA was extracted from samples using conventional methodologies and quantified using PicoGreen (Invitrogen).
Genotyping of samples was performed using customized Illumina (San Diego, CA) Sentrix Bead Arrays according to the manufacturer's protocols. Annotated flanking sequence information for each SNP was derived from the University of California Santa Cruz (UCSC) Human Genome Browser (Assembly hg17). SNP selection was in part restricted to those that are amenable to genotyping using the array-based platform used. DNA samples with GenCall scores < 0.25 at any locus were considered "no calls." It has been proposed that alleles linked to complex disease are likely to reflect modest changes in gene activity, therefore non-coding SNPs are as likely as coding SNPs to be associated with complex disease risk . Hence our analysis was largely based on the genotyping non-coding polymorphisms in introns and within 2 kb of the mRNA of CYP1A2, CYP1B1, CYP3A4, CYP3A5, CYP11A1, CYP17A1 and CYP19A1. We did however make use of genotypes generated from a genome-wide association study of coding SNPs we have recently performed (two SNPs in CYP1B1 and one in CYP3A4 and CYP3A5A) to assist in haplotype construction. Table 1 provides details of all 20 SNPs analysed in this current report.
The relationships between homozygote and heterozygote carrier status and risk of CRC and between number of putative risk alleles carried and risk of CRC were assessed by means of odds ratios (ORs) with 95% confidence intervals (CIs), calculated using logistic regression (adjusting for age and sex). Where it was not possible to calculate ORs and their CIs by asymptotic methods, an exact approach was implemented. A likelihood ratio test was performed to evaluate the impact of each SNP on case-control status. To test for population stratification, the distribution of genotypes in controls was tested for departures from Hardy-Weinberg equilibrium (HWE). To investigate epistatic interactions, each pair of SNPs was evaluated by fitting a saturated logistic regression model and the log likelihood ratio statistic for comparison with the main effects model computed. To assess the level of linkage disequilibrium (LD) between SNPs, we calculated the pair-wise LD measure D' between consecutive pairs of markers in the genes CYP1B1, CYP3A4, and CYP3A5 using the expectation-maximization algorithm to estimate two-locus haplotype frequencies. The program PHASE [11, 12] which implements a Markov chain Monte Carlo method was used to generate haplotypes. Haplotype frequencies in cases and controls were compared by the chi-squared statistic. We report p-values both before and after adjustment for multiple testing, which was carried out by means of the Bonferroni correction. All computations were undertaken using the statistical software packages Stata (Stata Corp, TX, USA) or LogXact (Cytel Inc., Cambridge, MA, USA).
Genotypes were obtained for 2,561 of the 2,575 cases (99.5%) and 2,695 of the 2,707 controls (99.6%). SNP call rates per sample for each of the 5,256 DNA samples were > 99.8% in cases and > 99.5% in controls. CYP genotypes and allele frequencies in CRC cases and controls are detailed in Table 2. Genotypic frequencies of SNPs in controls were similar to those previously documented in Caucasians and none were found to violate HWE in controls. Table 2 also details the risk associated with each genotype as quantified by the raw OR and the OR adjusted for age and gender, and their associated 95% CIs. Two SNPs showed some evidence of an association with CRC risk: rs162558, which maps 5' to CYP1B1, and rs2069522, which maps 5' to CYP1A2. Homozygotes for rs162558 in CYP1B1 and heterozygotes for rs2069522 in CYP1A2 were associated with a mildly increased risk of CRC (OR = 1.36, 95% CI: 1.03–1.80 and OR = 1.34, 95% CI: 1.00–1.79 respectively), although neither of these findings were significant after correction for multiple testing (adjusted p = 1.0). In the case of rs162558 in CYP1B1, heterozygotes showed no evidence of increased risk. The logistic regression model based on alleles did not yield any further significant associations.
Data on site of tumour and Duke's stage of tumour were available for 2,544 and 587 cases respectively. Of the patients for whom site of tumour was known,1,585 had colon cancer and 959 had rectal tumours. Of the patients for whom data regarding Duke's stage were available, 52 were stage A, 227 were stage B, 305 were stage C and 3 were stage D. Stratification of the analysis by each of these categories did not significantly affect findings.
Estimation of haplotypes in the genes with multiple SNPs resulted in 5, 4 and 3 common haplotypes (i.e. > 1%) for CYP1B1, CYP3A4 and CYP3A5 respectively. There was strong LD between the five common SNPs in CYP1B1 (D' > 0.97 for each pair of SNPs), the four common SNPs in CYP3A4 (D' > 0.95 for each pair of SNPs) and the three SNPs in CYP3A5 (D' > 0.84 for each pair of SNPs). When haplotype frequencies were compared in patients and controls no increased evidence of a relationship with CRC was apparent (Table 3). The joint effects of the polymorphisms on the risk of CRC were evaluated. Twelve pairs of the SNPs displayed nominal evidence of interaction at the 5% level, but none were significant after applying a Bonferroni correction for the 105 likelihood ratio tests undertaken.
Two of the variants we evaluated, CYP1A2 rs2069522 and CYP1B1 rs162558, showed mild evidence of an impact on CRC risk. For CYP1B1 rs162558, homozygosity was associated with increased CRC risk, and for CYP1A2 rs2069522 heterozygosity was associated with risk, although these associations were not significant after correction for multiple testing. The Bonferroni correction is based on the assumption that tests are independent. This is a conservative adjustment for these data due to the high levels of LD in CYP1B1, CYP3A4 and CYP3A5. However, the smallest observed unadjusted p-value of 0.05 would not be significant after correction for multiple testing even if a less conservative adjustment (such as adjusting for the ~7 independent tests once LD is taken into account) were used. Conflicting results from other studies have found these genes to be either associated or not associated with CRC risk [7, 13]. Discrepancies between these previous studies may be due to small sample sizes, confounding interactions with environmental or genetic risk factors, and the choice of SNPs to be genotyped. We have tested 2575 cases and 2707 controls, over 5-fold more than any previous study of CRC involving these genes, and hence had far more power to detect moderate increases in risk.
We acknowledge that the choice of particular SNPs and/or interactions with other risk factors could moderate the observed results in any association study. Previous studies have implicated polymorphisms in the cytochrome P450 genes CYP1A1, CYP2A6, CYP2C9, and CYP2C19 in risk of CRC (reviewed in ). However the small sample sizes of many early association studies may have led to misleading results. While the precise nature of the role of the CYP genes in CRC risk is not yet clear, our findings expand the body of knowledge about these genes, and we hope they will contribute to the development of further research in this area.
As the genotypic frequencies of SNPs in our controls were similar to those previously documented in Caucasians and none were found to violate Hardy-Weinberg equilibrium there is no evidence that cryptic relatedness or population stratification impacted on study findings. Although we do not have participant rate data it seems unlikely that any selection from this would bias study findings (i.e. specific genotype more likely to have been ascertained).
SNP rs2069522 is located -2847 bp relative to CYP1A2, within a putative region of bi-directional suppressive transcriptional control (position -1329 to -4412) for both CYP1A1 and CYP1A2 genes , hence it is possible that this sequence change may influence transcription of CYP1A1 and CYP1A2 proteins. Previously, SNPs in both CYP1A1 and CYP1A2 have been associated with colorectal adenomas and carcinomas [7, 16], therefore the observation of some relationship between CYP1A2 rs2069522 and risk in our study is not without precedent. Moreover, CYP1A2 is critical for the metabolic activation of dietary heterocyclic aromatic amines (HCA) which are mutagenic and have been implicated in the development of CRC [17–23].
SNP rs162558 is located -1548 bp relative to CYP1B1. CYP1B1 activates several human procarcinogens by an aryl hydrocarbon receptor-aryl hydrocarbon receptor nuclear translocator pathway [24, 25]. CYP1B1 is overexpressed in colorectal adenocarcinomas relative to normal colon tissue , and a variant with increased activity towards several substrates including sex hormones has been associated with increased risk of CRC . Our data support the evidence that genetic variants in CYP1B1 may be associated with CRC risk.
CYP19A1 rs10046 has previously been reported to influence levels of estradiol  and possibility impact on breast cancer risk. While estrogens undoubtedly affect CRC risk, in our study we found no evidence that variation in CYP19A1 defined by rs10046 had an impact on CRC risk.
Our study is large compared to contemporaneous studies, and although for the majority of SNPs assayed our analysis has 80% power to detect a relative risk of 1.5, inevitably we cannot entirely exclude a small effect in risk of CRC associated with the variants analysed. Moreover, as there is increasing evidence of a gene-environment effect with respect to CRC risk some of the polymorphisms may mediate CRC risk in the context of specific dietary risk factors. The primary aim of our study design was to acquire enough samples to achieve statistically significant results for the identification of common, low penetrance alleles in CRC. The addition of gene-environment interaction to our aims would require a vastly larger number of samples in order achieve significance. We hope that as technological advances enable studies of this nature to include larger numbers of patient samples, reliable information regarding carcinogen exposure and diet will be informative, but the collection and analysis of environmental factors was unfortunately not possible on the scale of our study.
This study provides some support for polymorphic variation in CYP1A2 and CYP1B1 playing a role in CRC susceptibility.
Cross AJ, Sinha R: Meat-related mutagens/carcinogens in the etiology of colorectal cancer. Environ Mol Mutagen. 2004, 44 (1): 44-55. 10.1002/em.20030.
Lowenfels AB, Anderson ME: Diet and cancer. Cancer. 1977, 39 (4 Suppl): 1809-1814. 10.1002/1097-0142(197704)39:4+<1809::AID-CNCR2820390811>3.0.CO;2-O.
Palmer S: Diet, nutrition, and cancer. Prog Food Nutr Sci. 1985, 9 (3-4): 283-341.
Giovannucci E: An updated review of the epidemiological evidence that cigarette smoking increases risk of colorectal cancer. Cancer Epidemiol Biomarkers Prev. 2001, 10 (7): 725-731.
Rodriguez-Antona C, Ingelman-Sundberg M: Cytochrome P450 pharmacogenetics and cancer. Oncogene. 2006, 25 (11): 1679-1691. 10.1038/sj.onc.1209377.
Ding X, Kaminsky LS: Human extrahepatic cytochromes P450: function in xenobiotic metabolism and tissue-selective chemical toxicity in the respiratory and gastrointestinal tracts. Annu Rev Pharmacol Toxicol. 2003, 43: 149-173. 10.1146/annurev.pharmtox.43.100901.140251.
Landi S, Gemignani F, Moreno V, Gioia-Patricola L, Chabrier A, Guino E, Navarro M, de Oca J, Capella G, Canzian F: A comprehensive analysis of phase I and phase II metabolism gene polymorphisms and risk of colorectal cancer. Pharmacogenet Genomics. 2005, 15 (8): 535-546.
Dunning AM, Dowsett M, Healey CS, Tee L, Luben RN, Folkerd E, Novik KL, Kelemen L, Ogata S, Pharoah PD, Easton DF, Day NE, Ponder BA: Polymorphisms associated with circulating sex hormone levels in postmenopausal women. J Natl Cancer Inst. 2004, 96 (12): 936-945.
Zheng W, Gao YT, Shu XO, Wen W, Cai Q, Dai Q, Smith JR: Population-based case-control study of CYP11A gene polymorphism and breast cancer risk. Cancer Epidemiol Biomarkers Prev. 2004, 13 (5): 709-714.
Carlson CS, Eberle MA, Kruglyak L, Nickerson DA: Mapping complex disease loci in whole-genome association studies. Nature. 2004, 429 (6990): 446-452. 10.1038/nature02623.
Stephens M, Donnelly P: A comparison of bayesian methods for haplotype reconstruction from population genotype data. Am J Hum Genet. 2003, 73 (5): 1162-1169. 10.1086/379378.
Stephens M, Smith NJ, Donnelly P: A new statistical method for haplotype reconstruction from population data. Am J Hum Genet. 2001, 68 (4): 978-989. 10.1086/319501.
Sachse C, Smith G, Wilkie MJ, Barrett JH, Waxman R, Sullivan F, Forman D, Bishop DT, Wolf CR: A pharmacogenetic study to investigate the role of dietary carcinogens in the etiology of colorectal cancer. Carcinogenesis. 2002, 23 (11): 1839-1849. 10.1093/carcin/23.11.1839.
Agundez JA: Cytochrome P450 gene polymorphism and cancer. Curr Drug Metab. 2004, 5 (3): 211-224. 10.2174/1389200043335621.
Ueda R, Iketaki H, Nagata K, Kimura S, Gonzalez FJ, Kusano K, Yoshimura T, Yamazoe Y: A common regulatory region functions bidirectionally in transcriptional activation of the human CYP1A1 and CYP1A2 genes. Mol Pharmacol. 2006, 69 (6): 1924-1930. 10.1124/mol.105.021220.
Moonen H, Engels L, Kleinjans J, Kok T: The CYP1A2-164A-->C polymorphism (CYP1A2*1F) is associated with the risk for colorectal adenomas in humans. Cancer Lett. 2005, 229 (1): 25-31. 10.1016/j.canlet.2004.12.010.
Adamson RH, Thorgeirsson UP, Snyderwine EG, Thorgeirsson SS, Reeves J, Dalgard DW, Takayama S, Sugimura T: Carcinogenicity of 2-amino-3-methylimidazo[4,5-f]quinoline in nonhuman primates: induction of tumors in three macaques. Jpn J Cancer Res. 1990, 81 (1): 10-14.
Giovannucci E, Rimm EB, Stampfer MJ, Colditz GA, Ascherio A, Willett WC: Intake of fat, meat, and fiber in relation to risk of colon cancer in men. Cancer Res. 1994, 54 (9): 2390-2397.
Giovannucci E, Stampfer MJ, Colditz G, Rimm EB, Willett WC: Relationship of diet to risk of colorectal adenoma in men. J Natl Cancer Inst. 1992, 84 (2): 91-98. 10.1093/jnci/84.2.91.
Goldbohm RA, van den Brandt PA, van 't Veer P, Brants HA, Dorant E, Sturmans F, Hermus RJ: A prospective cohort study on the relation between meat consumption and the risk of colon cancer. Cancer Res. 1994, 54 (3): 718-723.
Kono S, Imanishi K, Shinchi K, Yanai F: Relationship of diet to small and large adenomas of the sigmoid colon. Jpn J Cancer Res. 1993, 84 (1): 13-19.
Nagao M, Ushijima T, Watanabe N, Okochi E, Ochiai M, Nakagama H, Sugimura T: Studies on mammary carcinogenesis induced by a heterocyclic amine, 2-amino-1-methyl-6-phenylimidazo[4,5-b]pyridine, in mice and rats. Environ Mol Mutagen. 2002, 39 (2-3): 158-164. 10.1002/em.10047.
Ochiai M, Ogawa K, Wakabayashi K, Sugimura T, Nagase S, Esumi H, Nagao M: Induction of intestinal adenocarcinomas by 2-amino-1-methyl-6-phenylimidazo[4,5-b]pyridine in Nagase analbuminemic rats. Jpn J Cancer Res. 1991, 82 (4): 363-366.
Shehin SE, Stephenson RO, Greenlee WF: Transcriptional regulation of the human CYP1B1 gene. Evidence for involvement of an aryl hydrocarbon receptor response element in constitutive expression. J Biol Chem. 2000, 275 (10): 6770-6776. 10.1074/jbc.275.10.6770.
Shimada T, Hayes CL, Yamazaki H, Amin S, Hecht SS, Guengerich FP, Sutter TR: Activation of chemically diverse procarcinogens by human cytochrome P-450 1B1. Cancer Res. 1996, 56 (13): 2979-2984.
Gibson P, Gill JH, Khan PA, Seargent JM, Martin SW, Batman PA, Griffith J, Bradley C, Double JA, Bibby MC, Loadman PM: Cytochrome P450 1B1 (CYP1B1) is overexpressed in human colon adenocarcinomas relative to normal colon: implications for drug development. Mol Cancer Ther. 2003, 2 (6): 527-534.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2407/7/123/prepub
Grant support: Funding for this work was undertaken with support from Cancer Research UK, CORE, the National Cancer Research Network and the European Union (CCPRB). We gratefully acknowledge the participation of all patients and control individuals.
The author(s) declare that they have no competing interests.
LB analysed the data and drafted the manuscript. EW performed the statistical analysis and helped to draft the manuscript. GS and MR participated in the design and coordination of the study. SP, LW, and MQ carried out the study. RH conceived of the study, participated in its design and coordination, and helped to draft the manuscript.
Lara Bethke, Emily Webb contributed equally to this work.
Rights and permissions
Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
About this article
Cite this article
Bethke, L., Webb, E., Sellick, G. et al. Polymorphisms in the cytochrome P450 genes CYP1A2, CYP1B1, CYP3A4, CYP3A5, CYP11A1, CYP17A1, CYP19A1and colorectal cancer risk. BMC Cancer 7, 123 (2007). https://doi.org/10.1186/1471-2407-7-123
- Single Nucleotide Polymorphism
- Cytochrome P450 Gene
- Single Nucleotide Polymorphism Genotyping
- Single Nucleotide Polymorphism Rs162558
- CYP1A2 Rs2069522