Abstract
Background: The gene encoding the phase I enzyme cytochrome P4502D6 (CYP2D6) has been previously investigated for its potential predictive role in the efficacy of breast cancer treatments such as tamoxifen, but its role in breast cancer susceptibility is unclear. This study aims to evaluate the association between germ line variations in CYP2D6 and breast cancer susceptibility.
Methods: DNA samples from 13,472 cases and controls were genotyped for seven known functional variants [minor allele frequency (MAF) ≥ 0.01] and five single nucleotide polymorphisms (SNP) that tag common genetic variation (MAF > 0.05) in CYP2D6.
Results: One relatively rare functional variant, CYP2D6*6, (MAF = 0.01) showed a modest increased association with breast cancer susceptibility (Ptrend = 0.02; OR = 1.32; 95% CI = 1.04–1.68). All other functional and tagSNPs showed no association with breast cancer susceptibility.
Conclusions: Common variants of CYP2D6 do not play a significant role in breast cancer susceptibility. However, this study raises questions regarding the role of rare variants, such as CYP2D6*6, in breast cancer susceptibility which merit further investigation.
Impact: This large case–control study, involving 13,472 women, found no evidence of any association between common CYP2D6 gene variants and breast cancer susceptibility. However, one relatively rare functional variant CYP2D6*6 showed a modest association with breast cancer susceptibility, indicating that the role of rare CYP2D6 variants in breast cancer risk is unclear and requires further investigation in an adequately powered study. Cancer Epidemiol Biomarkers Prev; 20(6); 1255–8. ©2011 AACR.
Introduction
Genome-wide association studies provide empirical evidence that part of the inherited genetic component of breast cancer risk is due to common alleles (1). Coverage of common genetic variation by the various genome-wide single nucleotide polymorphism (SNP) arrays is incomplete. Coverage of variants with a minor allele frequency (MAF) of less than 0.05 is very poor. Breast cancer development is multifactorial, involving genetic and environmental factors.
Genes encoding enzymes within metabolic pathways may be important because of their role in detoxification of xenobiotic compounds. Conversely, they may potentiate the risk of cancer by activating carcinogens. Environmental exposures can affect cancer susceptibility genes (2). Cytochrome P4502D6 (CYP2D6) is a polymorphic gene involved in phase I metabolism (3). Many CYP2D6 variants have altered function which may influence cancer susceptibility. Variants are categorized into poor (PM); intermediate (IM); extensive (EM); and ultra (UM) metabolizers.
The role of CYP2D6 in breast cancer susceptibility is uncertain. Nine published studies, with heterogeneous methodology, study design, and analysis, have given conflicting results. Four studies found an association between CYP2D6 variants and breast cancer susceptibility (4). However, 5 studies showed either no supportive evidence or contradictory evidence (5).
Methods and Materials
These methods and materials have previously been published (3).
Cases and controls
The SEARCH breast cancer case–control study and the characteristics of the cases used in this research have been described previously (3, 6). The total number of cases and controls available for analysis were 6,640 and 6,832, respectively. Eastern Region Multicentre Research Ethics Committee approved the study. Participants provided written informed consent.
Selection of tagging and functional SNPs
Selected tagSNPs tagged all known common SNPs (MAF ≥ 0.05) in the CYP2D6 gene, with a pairwise correlation of r2 > 0.8. Seven functional SNPs were successfully genotyped (3).
Genotyping
Genotyping was conducted using TaqMan (Applied Biosystems) according to manufacturer's instructions. Primers and probes were supplied by Applied Biosystems as Assays-by-Design. Plates were read on the ABI Prism 7900 using the Sequence Detection Software.
A nested PCR approach was used for CYP2D6*4 and a biplex TaqMan real-time quantification assay was used to identify the CYP2D6*5 (gene deletion) and CYP2D6*UM (gene duplication) variants (3).
Statistics
For each polymorphism, deviation of the genotype frequencies from those expected under Hardy–Weinberg equilibrium was assessed using a χ2 test. Genotype frequencies in cases and controls were compared using a χ2 test with 2 degrees of freedom (df; Pheterogeneity) and the Cochrane–Armitage trend test (χ2 on 1df) for the per allele breast cancer risk (Ptrend). The relative risks of breast cancer for heterozygotes and rare homozygotes, relative to common homozygotes, were estimated as ORs with associated 95% CIs.
The admixture maximum likelihood (AML) test (7) tests the null hypothesis that none of the SNPs are associated with disease compared with the alternative that 1 or more of SNPs are associated.
Comparison of haplotype frequencies by haplotype block, and subject-specific expected haplotype indicators were calculated using an in-house program (6). Haplotypes with greater than 2% frequency were considered “common.” Rare haplotypes were pooled. Logistic regression generated haplotype specific risks with respect to all other common haplotypes present (and “rare” grouping) as ORs with 95% CIs. All analyses were conducted in Intercooled Stata version 10.
Results
Using HapMap data, 5 tagSNPs tagged 46 common SNPs (MAF ≥ 0.05) in a 100-kb region centered on CYP2D6. Seven key functional SNPs, (MAF ≥ 0.01), were identified. All variants were genotyped in 13,472 cases and controls.
The AML test did not indicate any evidence for an association with breast cancer susceptibility [Ptrend = 0.24 linkage disequilibrium (LD) block1; Ptrend = 0.78 for LD block2].
Only 1 SNP had a modest association with increased breast cancer susceptibility (at Ptrend ≤ 0.05). CYP2D6*6 heterozygotes showed an increased breast cancer risk (P = 0.02; OR = 1.32; 95% CI = 1.04–1.68). No rare homozygotes for CYP2D6*6 (MAF = 0.01) were detected. No association was seen with any other tagSNPs or functional SNPs (Table 1).
Breast cancer risk associated with CYP2D6 functional and tagSNPs
SNPs . | Metabolizer status . | MAFf . | h OR (95% Cl; heterozygotes vs. common homozygotes) . | OR (95% Cl; rare homozygotes vs. common homozygotes) . | Heterogeneity P (2df) . | Trend test P (1df) . |
---|---|---|---|---|---|---|
Functional SNPs | ||||||
CYP2D6*41 | (IM) | 0.09 | 0.98 [0.89–1.07] | 0.95 [0.71–1.28] | 0.85 | 0.58 |
CYP2D6*04 | (PM) | 0.2 | 1.00 [0.92–1.07] | 0.86 [0.73–1.02] | 0.23 | 0.26 |
CYP2D6*05a | (PM) | 0.04 | N/Ac | 3.08 [0.83–11.37] | 0.07 | 0.09 |
CYP2D6*06 | (PM) | 0.01 | 1.33 [1.05–1.69] | N/Ag | 0.02 | 0.02 |
CYP2D6*09 | (IM) | 0.03 | 1.04 [0.90–1.21] | N/Ag | 0.59 | 0.45 |
CYP2D6*10 | (IM) | 0.02 | 1.07 [0.85–1.36] | N/Ag | 0.55 | 0.55 |
CYP2D6*UMb | (UM) | 0.08 | N/Ac | 1.25 [0.79–1.98] | 0.35 | 0.35 |
CYP2D6*PM Model 1d | (PM) | 1N/A | Trend test only:h | 0.99 [0.89–1.10] | N/A | 0.89 |
CYP2D6*PM Model 2e | (PM) | 1N/A | Trend test only:h | 0.98 [0.89–1.07] | N/A | 0.65 |
TagSNPs | ||||||
CYP2D6 01t | N/A | 0.32 | 0.92 [0.86–0.99] | 1.04 [0.93–1.18] | 0.03 | 0.53 |
CYP2D6 02t | N/A | 0.46 | 1.04 [0.96–1.13] | 1.06 [0.96–1.17] | 0.43 | 0.2 |
CYP2D6 03t | N/A | 0.24 | 0.95 [0.89–1.02] | 1.03 [0.88–1.21] | 0.31 | 0.48 |
CYP2D6_04t | N/A | 0.08 | 0.99 [0.90–1.09] | 0.98 [0.65–1.48] | 0.98 | 0.83 |
CYP2D6_05t | N/A | 0.22 | 0.98 [0.91–1.06] | 0.86 [0.73–1.01] | 0.19 | 0.15 |
SNPs . | Metabolizer status . | MAFf . | h OR (95% Cl; heterozygotes vs. common homozygotes) . | OR (95% Cl; rare homozygotes vs. common homozygotes) . | Heterogeneity P (2df) . | Trend test P (1df) . |
---|---|---|---|---|---|---|
Functional SNPs | ||||||
CYP2D6*41 | (IM) | 0.09 | 0.98 [0.89–1.07] | 0.95 [0.71–1.28] | 0.85 | 0.58 |
CYP2D6*04 | (PM) | 0.2 | 1.00 [0.92–1.07] | 0.86 [0.73–1.02] | 0.23 | 0.26 |
CYP2D6*05a | (PM) | 0.04 | N/Ac | 3.08 [0.83–11.37] | 0.07 | 0.09 |
CYP2D6*06 | (PM) | 0.01 | 1.33 [1.05–1.69] | N/Ag | 0.02 | 0.02 |
CYP2D6*09 | (IM) | 0.03 | 1.04 [0.90–1.21] | N/Ag | 0.59 | 0.45 |
CYP2D6*10 | (IM) | 0.02 | 1.07 [0.85–1.36] | N/Ag | 0.55 | 0.55 |
CYP2D6*UMb | (UM) | 0.08 | N/Ac | 1.25 [0.79–1.98] | 0.35 | 0.35 |
CYP2D6*PM Model 1d | (PM) | 1N/A | Trend test only:h | 0.99 [0.89–1.10] | N/A | 0.89 |
CYP2D6*PM Model 2e | (PM) | 1N/A | Trend test only:h | 0.98 [0.89–1.07] | N/A | 0.65 |
TagSNPs | ||||||
CYP2D6 01t | N/A | 0.32 | 0.92 [0.86–0.99] | 1.04 [0.93–1.18] | 0.03 | 0.53 |
CYP2D6 02t | N/A | 0.46 | 1.04 [0.96–1.13] | 1.06 [0.96–1.17] | 0.43 | 0.2 |
CYP2D6 03t | N/A | 0.24 | 0.95 [0.89–1.02] | 1.03 [0.88–1.21] | 0.31 | 0.48 |
CYP2D6_04t | N/A | 0.08 | 0.99 [0.90–1.09] | 0.98 [0.65–1.48] | 0.98 | 0.83 |
CYP2D6_05t | N/A | 0.22 | 0.98 [0.91–1.06] | 0.86 [0.73–1.01] | 0.19 | 0.15 |
NOTE: Abbreviation: N/A, not applicable.
Cases n = 6,640; Controls = 6,832.
aCYP2D6*05 - gene deletion
bCYP2D6*UM - gene duplication
cCYP2D6*05 - gene deletion and CYP2D6*UM expression levels assessed using real-time PCR, but expression levels could not adequately distinguish between heterozygotes and common homozygotes.
dThe PM group are classified as carriers of 2 variant alleles (rare homozygote alleles) for at least 1 of the functional SNPs associated with PM or IM status. Where IM = CYP2D6*41; CYP2D6*9; CYP2D6*10; and PM = CYP2D6*4; CYP2D6*5; CYP2D6*6. In CYP2D6* PM model 1: The PM group (classified as stated) is compared with individuals who carried 2 copies of the wild-type allele (EM) at all SNPs (i.e., common homozygotes) or individuals who carried a single variant allele at a single SNP (heterozygotes).
eThe PM group was classified as individuals carrying at least 1 variant allele at 1 or more of the functional SNPs (heterozygotes and rare homozygotes), associated with PM or IM status. Where IM = CYP2D6*41; CYP2D6*9; CYP2D6*10; and PM = CYP2D6*4; CYP2D6*5; CYP2D6*6. In CYP2D6* PM model 2: The PM group (classified as stated) is compared with individuals who carried 2 copies of the wild-type allele (EM) at all SNPs (common homozygotes).
fMAF shown for controls
gN/A - Either no or very few homozygotes for this variant within study population.
hTrend test only given as these are combined variants.
Discussion
There is little evidence to suggest that common or functional variants are associated with breast cancer risk. CYP2D6*6 was associated with a modestly increased breast cancer risk (P = 0.02). Given that there are many potential low to moderate penetrance breast cancer susceptibility genes and the prior probability of an underlying association for any particular gene is low, this finding is likely to represent a false positive. The SNPs in this study are in pairwise correlation (r2 ≥ 0.8) with 77% of all the common variants present in the resequenced CYP2D6 region available in the 1000 Genomes Project (8).
To assess if any, as yet undetermined, CYP2D6 variant could be associated with susceptibility, a global haplotype analysis was conducted (Table 2). No significant differences in haplotype frequencies were seen between cases and controls. The combined rare haplotype group was marginally statistically significant; however, this study does not have sufficient power to investigate rare variants.
Global haplotype analysis of CYP2D6 functional and tag SNPs
Genea . | Haplotypeb . | Frequency . | P . | dOR . | d95% Cls . | Global P . |
---|---|---|---|---|---|---|
CYP2D6 | h00100000 | 0.22 | 0.55 | 0.06 | ||
Block 1 | h00110000 | 0.09 | 0.47 | |||
5t, 2t, 1t, *41, *9, *6, *4, *10 | h01001000 | 0.03 | 0.62 | |||
h10000010 | 0.1 | 0.19 | ||||
Combined rarec | 0.04 | 0.01 | 1.18 | 1.04–1.35 | ||
CYP2D6+ | h01 | 0.08 | 0.73 | 0.74 | ||
Block 2 | h10 | 0.23 | 0.46 | |||
3t, 4t |
Genea . | Haplotypeb . | Frequency . | P . | dOR . | d95% Cls . | Global P . |
---|---|---|---|---|---|---|
CYP2D6 | h00100000 | 0.22 | 0.55 | 0.06 | ||
Block 1 | h00110000 | 0.09 | 0.47 | |||
5t, 2t, 1t, *41, *9, *6, *4, *10 | h01001000 | 0.03 | 0.62 | |||
h10000010 | 0.1 | 0.19 | ||||
Combined rarec | 0.04 | 0.01 | 1.18 | 1.04–1.35 | ||
CYP2D6+ | h01 | 0.08 | 0.73 | 0.74 | ||
Block 2 | h10 | 0.23 | 0.46 | |||
3t, 4t |
aSNPs used for haplotype analysis are arranged in chromosome order; 0 = major allele; 1 = minor allele.
bThe baseline for each test is the frequency of all the other haplotypes combined.
cRare haplotypes (≤ 5%) were pooled.
dOR and 95% CI only stated if haplotypes show a difference between cases and controls at a significance level P ≤ 0.05.
In conclusion, there is limited evidence to suggest that CYP2D6 plays a significant role in breast cancer susceptibility. However, the role of rarer variants, such as CYP2D6*6, requires further investigation.
Disclosure of Potential Conflicts of Interest
No potential conflicts of interest were disclosed.
Acknowledgments
The authors thank the women participating in this study; SEARCH study team; Eastern Cancer Registration and Information Centre; clinicians throughout East Anglia who recruited patients; Professor Ulrich M. Zanger and team (Dr. Margarete Fischer-Bosch Institute of Clinical Pharmacology, Germany) for providing DNA samples with known CYP2D6 genotypes.
Grant Support
This research was funded by Cancer Research UK (CRUK) grant refs: C10097/A7484 and C507/A6306 (J.E. Abraham during the period of this research); C490/A10119 (P.D.P. Pharoah and K.E. Driver); C8197/A10123 (A.M. Dunning and C. Luccarini); RHAG109 (M.J. Maranian); and C507/A7325 (C. Caldas).