|Home | About | Journals | Submit | Contact Us | Français|
Cytochrome P450 2D6 (CYP2D6) inhibition reduces the concentration of 4-hydroxylated tamoxifen metabolites, but the clinical relevance remains uncertain.
We conducted a large case–control study nested in the population of 11251 women aged 35–69 years at diagnosis of stage I–III breast cancer between 1985 and 2001 on Denmark’s Jutland Peninsula and registered with the Danish Breast Cancer Cooperative Group. We identified 541 recurrent or contralateral breast cancers among women with estrogen receptor–positive (ER+) disease treated with tamoxifen for at least 1 year and 300 cancers in women with ER-negative (ER−) disease never treated with tamoxifen. We matched one control subject per case patient on ER status, menopausal status, stage, calendar time, and county, genotyped the CYP2D6*4 allele to assess genetic inhibition, and ascertained prescription history to assess drug–drug inhibition. We estimated the odds ratio (OR), associating CYP2D6 inhibition with breast cancer recurrence and adjusted for potential confounding with logistic regression. To address bias from incomplete information on CYP2D6 function, we used Monte Carlo simulation to complete a record-level probabilistic bias analysis. All statistical tests were two-sided.
The frequency of the CYP2D6*4 minor allele was 24% in case patients with ER+ tumors, 23% in case patients with ER− tumors, and 22% each in control subjects with ER+ and ER− tumors. In women with ER+ tumors, the associations of one functional allele with recurrence (OR = 0.99; 95% confidence interval = 0.76 to 1.3) and no functional allele with recurrence (OR = 1.4; 95% confidence interval = 0.84 to 2.3) were near null, as were those for women with ER− tumors. The near-null associations persisted when evaluated by intake of medications, by combining genotype with medication history, in the probabilistic bias analysis, or by restricting the analysis to women with ER expression confirmed by re-assay.
The association between CYP2D6 inhibition and recurrence in tamoxifen-treated patients is likely null or small.
The cytochrome P450 2D6 (CYP2D6) enzyme, which metabolizes tamoxifen, is inhibited by the selective serotonin reuptake inhibitor paroxetine, but it is not known whether women with fewer than two functional CYP2D6*4 alleles or those who take selective serotonin reuptake inhibitors are poor candidates for tamoxifen therapy.
Five hundred and forty -one women in the Danish Breast Cancer Cooperative Group registry with recurrent or contralateral estrogen receptor–positive breast cancer who were treated with tamoxifen and 300 women with estrogen receptor–negative breast cancer who were never treated with tamoxifen were matched on clinical and tumor characteristics, CYP2D6 genotype, and selective serotonin reuptake inhibitor prescription history with control subjects from the same registry who had no recurrent or contralateral breast cancer.
There was no statistically significant association between CYP2D6 inhibition and breast cancer recurrence in tamoxifen-treated women. The near-null association persisted regardless of whether CYP2D6 inhibition was assessed by genotype, by intake of medications that inhibit CYP2D6 function, or by a combination of genotype and medication history.
Tamoxifen treatment can be effective in women with estrogen receptor –positive breast cancer who have fewer than two functional CYP2D6 alleles or who take medications, such as selective serotonin reuptake inhibitors, that inhibit the CYP2D6 enzyme.
Genotyping data for only one CYP2D6 allele were available, so the association between other CYP2D6 alleles and breast cancer recurrence was not ascertained. There was no information on tamoxifen adherence by case patients and control subjects, so the full extent of tamoxifen treatment was unknown.
From the Editors
Tamoxifen (TAM), a selective estrogen receptor (ER) modulator, halves the risk of breast cancer recurrence in patients with nonmetastatic ER-positive (ER+) disease and is also a potent therapy in women with metastatic ER+ disease (1). The effectiveness of tamoxifen therapy is, however, incomplete. Some women relapse and others do not respond at all. Mechanisms of resistance to tamoxifen therapy and predictive markers of susceptibility to resistance other than lack of ER expression have been widely researched (2–4). Accurate markers are clinically important, allowing prediction of tamoxifen response, adverse effects, and personalization of combined therapies (5,6). Tamoxifen and its primary metabolite (N-desmethyl tamoxifen) are metabolized mostly by the gene product of cytochrome P450 2D6 (CYP2D6) (7,8) to 4-hydroxytamoxifen (9) and 4-hydroxy-N-desmethyltamoxifen (10,11) (now often called endoxifen). These secondary metabolites bind the receptor about 100-fold more readily than tamoxifen and are thus the most important modulators of the ER in the tamoxifen pathway (12). Women who inherit two functional CYP2D6 alleles have higher steady-state concentrations of 4-hydroxytamoxifen (8,13) and 4-hydroxy-N-desmethyltamoxifen (7,13,14) than women who inherit no functional alleles when treated with tamoxifen. Women who inherit one functional allele have intermediate concentrations (13). Similarly, women who inherit two functional alleles and take the potent CYP2D6 inhibitor paroxetine (a selective serotonin reuptake inhibitor [SSRI]) have lower concentrations of 4-hydroxy-N-desmethyltamoxifen (7,14) when treated with tamoxifen. It has been suggested that ER+ breast cancer patients with nonfunctional alleles of CYP2D6 or those who take other medications that inhibit CYP2D6 function may be poor candidates for adjuvant tamoxifen therapy (15–17) because their lower concentrations of the potent metabolites may place them at higher risk for relapse or failure to respond.
Although the molecular and pharmacological bases for this hypothesis are compelling, earlier clinical epidemiology studies focusing on associations between CYP2D6 inhibition and breast cancer outcomes have had widely heterogeneous results (18). Thirteen published studies have examined the association between inheriting a variant CYP2D6 allele and risk of breast cancer recurrence or mortality (19–31), and several have been recently updated (32,33). Relative risks reported in these studies range from 0.52 to 6.7, with six reporting associations below the null and seven reporting associations above the null (Phomogeneity < .001). A similarly heterogeneous pattern has been found in studies of the association between drug-induced CYP2D6 inhibition and breast cancer recurrence (34–41). Some of these studies have been criticized on the grounds of small sample size, survivor and other selection biases, potential for uncontrolled confounding by prognostic markers, and information bias arising from retrospective or absent information on use of CYP2D6 inhibiting medications or from noncentralized testing of ER expression [see (18) for a review]. Genotype frequencies were not in Hardy–Weinberg equilibrium in many of the earlier studies (20,22,26,29,31,32). The inconsistent pattern of associations from earlier studies, combined with limitations offering only a partial explanation for the heterogeneity, cautions against any strong inference based on results available to date.
To address limitations of earlier research and to provide a precise estimate obtained from a large well-identified study population, we conducted a study of the association between genetic and pharmacological markers of CYP2D6 inhibition and breast cancer recurrence nested in a population-based clinical registry of breast cancer patients in Denmark (42).
The source population consisted of 11251 female residents of the Jutland Peninsula in Denmark aged 35–69 years between 1985 and 2001, who were diagnosed with stage I, II, or III breast cancer as defined by the Union for International Cancer Control (UICC) (43), and who were registered with the Danish Breast Cancer Cooperative Group (DBCG) (42). Since 1977, the DBCG has enrolled into its clinical database nearly all Danish breast cancer patients younger than 70 years at diagnosis. The same 10-year follow-up protocol is used for all patients registered with the DBCG (44), regardless of their participation in clinical trials. Approximately half of the DBCG patients are enrolled in clinical trials (44). Thus, studies nested in the DBCG registry combine the data quality advantages of a clinical trial with the generalizability benefits of a representative population. The study was approved by the Regional Committee on Biomedical Research Ethics of Aarhus County in Denmark and by the Boston University Medical Campus Institutional Review Board.
We divided the source population into three groups: women whose tumor expressed the ER and who were treated with tamoxifen for at least 1 year (ER+/TAM+, n = 1826); women whose tumor did not express the ER, who were not treated with tamoxifen, and who survived for at least 1 year (ER−/TAM−, n = 1808); and all other breast cancer patients (n = 7617, such as ER+ patients not treated with tamoxifen and ER− patients who were treated with tamoxifen) who were excluded (Figure 1). ER+/TAM+ women were assigned to tamoxifen- therapy protocols lasting 1, 2, or 5 years, depending on the guideline prevailing in Denmark at the time of diagnosis (44). As discussed below, most women assigned to 1- or 2-year protocols actually were likely to receive tamoxifen therapy for much longer. We included ER−/TAM− women to estimate the direct association between CYP2D6 inhibition and the breast cancer recurrence rate. Follow-up time began 1 year after the date of breast cancer diagnosis and continued until the date of the first breast cancer recurrence, death from any cause, loss to follow-up (eg, emigration), 10 years of follow-up, or until September 1, 2006.
Case patients were women with local or distant breast cancer recurrence or contralateral breast cancer occurrence during their follow-up time. We designed the ER+/TAM+ sample size to achieve statistical power of 90% to detect an odds ratio (OR) of 1.5 associating reduced CYP2D6 function with recurrence risk. We therefore included all 541 ER+/TAM+ case patients. The number of ER−/TAM− case patients was greater than needed to achieve adequate statistical power. We therefore selected ER−/TAM− case patients at random (n = 300), after frequency matching as close as possible to the distribution of stage and calendar period of diagnosis among ER+/TAM+ case patients (Figure 1).
For each case patient, we selected without replacement one control subject from members of the source population who were alive and had no recurrence or contralateral breast cancer after the same amount of follow-up time. We matched control subjects to case patients according to group membership (ER+/TAM+ or ER−/TAM−), menopausal status at diagnosis (premenopausal or postmenopausal), date of breast cancer surgery (caliper matched ± 12 months), county of residence at time of diagnosis, and cancer stage at diagnosis (UICC stage I, II, or III).
We used the Danish Civil Personal Registration number to link datasets. The number is a unique identifier assigned to all Danish residents who were alive on April 1, 1968, born thereafter, or upon immigration.
We collected data from the DBCG registry on demographic (age, menopausal status, and hospital of diagnosis), tumor (UICC stage, histological grade, and ER expression), and therapy characteristics (primary surgical tumor management, receipt of radiation therapy, receipt of chemotherapy, and receipt of tamoxifen therapy). We collected information from the Danish National Registry of Patients on the conditions included in the Charlson comorbidity index (46) that were present at the time of breast cancer diagnosis.
We obtained data on receipt of prescriptions for SSRIs and other potential CYP2D6 inhibitors by linking to the Registry of Medicinal Products, which is maintained by Statistics Denmark as part of the Danish national health-care system. Because this registry’s records are only complete from 1995 forward , analyses incorporating prescription information are limited to case patients and control subjects diagnosed with breast cancer that year or later (Figure 1).
Laboratory personnel were blinded to all clinical information, including case or control status, ER status, and receipt or nonreceipt of tamoxifen therapy.
We collected formalin-fixed paraffin-embedded tissue blocks from the pathology department archives of treating hospitals. We reviewed hematoxylin- and eosin-stained glass slides and the original pathology reports to identify those blocks to be processed. All tissue blocks were processed using the laboratory’s established sterile protocols designed to avoid case contamination. For DNA extractions, three to six 10-μm paraffin sections were cut from each sample and placed in a 1.5-mL microtube. Parallel whole sections were cut from tumor blocks for confirmatory ER expression assays.
Before DNA extraction, tissue samples were deparaffinized by repeated treatment with xylene. After deparaffinization, DNA was extracted by repeated ethanol treatment, proteinase K digestion, and the QIAmp DNA FFPE Tissue kit (Qiagen AB, Dusseldorf, Germany) according to the manufacturer’s protocol.
From each tissue sample, 50 ng of extracted DNA were amplified with 25 μL polymerase chain reaction with 50 denaturation cycles at 92°C for 15 seconds, followed by annealing and extension at 60°C for 90 seconds, using primers and reagents supplied with TaqMan genotyping kits (ABI kit: C-27102431-D0, Applied BioSystems [ABI], Foster City, CA; https://products.appliedbiosystems.com/ab/en/US/adirect/ab?cmd=ABAssayDetailDisplay&assayID=C__27102431_D0&Fs=y) to assay CYP2D6*4 (rs3892097, the minor allele with no function). All samples were assayed in duplicate using the MX3000P real-time polymerase chain reaction system (Stratagene, Cedar Creek, TX). Positive and negative control assays for the variant were identified by sequencing peripheral blood DNA from 30 healthy individuals and were included with each set of assays.
For 106 participants, we paired lymph node tissue blocks with tumor tissue blocks and processed them according to the protocols for sterile technique, DNA extraction, and genotyping described above.
ER expression at diagnosis was the basis for inclusion in the study and segregation into ER+/TAM+ or ER−/TAM− groups. During the study period, the DBCG recommended an evolving series of assays for ER expression for all Danish breast cancer patients (47). Because assay methods have improved over time, and to reduce the potential effect of variability across diagnosing hospitals, we re-assayed ER expression using whole sections from the original diagnostic paraffin-embedded tissues and primary antibody against ER alpha (clone 6F11; Novocastra, Newcastle Upon Tyne, UK). Heat-induced epitope retrieval for ER was achieved by incubation in a Tris–EDTA buffer, pH 9 (Target Retrieval Solution, pH 9; Dako, Glostrup, Denmark) using a microwave oven. Sections were stained on a Lab Vision Autostainer (Thermo Fisher Scientific, Fremont, CA) using the EnVision+ detection system (Dako) and visualized with horseradish peroxidase and diaminobenzadine. Slides were scored positive for ER when there was distinct nuclear staining of neoplastic cells. The percentage of positive cells was recorded and for the purposes of this study, a cutoff point of at least 10% positive tumor was chosen for ER- positivity in accordance with previous DBCG recommendations (47).
We selected 20 ER+/TAM+ and 10 ER−/TAM− patients from one participating hospital. An investigator (J.P.G.) or a surgical colleague under his supervision reviewed their medical records blinded to the DBCG registry data. Half of the women in each group were diagnosed during the period 1985–1993 and half were diagnosed during the period 1994–2001. To guide the review of medical records, we adapted a standardized medical abstraction form and accompanying codebook from similar research tools used in earlier validation (48,49) and data collection (50) studies of breast cancer patients. The abstraction ascertained demographic information, tumor characteristics, therapy characteristics, recurrence status, and occurrence of a second primary breast cancer.
We used the DBCG definition of breast cancer recurrence, that is, any type of breast cancer or distant metastases diagnosed subsequent to the initial course of therapy (51). Given the follow-up time in the source population, all recurrences occurred between 1 and 10 years after the primary breast cancer diagnosis.
We classified case patients and control subjects as having two functional CYP2D6*4 alleles, one functional CYP2D6*4 allele, or no functional CYP2D6*4 allele.
Prescriptions in the Registry of Medicinal Products are coded using the Anatomical Therapeutic Chemical (ATC) classification system of the World Health Organization. We defined SSRI antidepressants as all medications included in ATC group N06AB. This group comprises the following drugs: zimeldine, fluoxetine, citalopram, paroxetine, sertraline, alaproclate, fluvoxamine, etoperidone, and escitalopram.
We classified case patients and control subjects as those with no record of a SSRI prescription during their follow-up time (never SSRI) and those with any record of a prescription for a SSRI during their follow-up time (ever SSRI). We used a similar procedure to classify case patients and control subjects as ever or never users of another prescription medication that is a CYP2D6 inhibitor or substrate, as defined previously (37). In earlier investigations using this study population, we observed a near-null association between use of SSRIs or other CYP2D6 inhibitors and recurrence risk (35–37).
To combine genetic and pharmacological information on CYP2D6 inhibition, we used the following index: 1) high inhibition: no functional allele or paroxetine or fluoxetine use during 30% or more of time on tamoxifen, 2) high intermediate inhibition: one functional allele and use of any CYP2D6 inhibitor (SSRI or other) during 30% or more of time on tamoxifen, 3) low intermediate inhibition: one functional allele and use of any CYP2D6 inhibitor during less than 30% of time on tamoxifen or no history of CYP2D6 inhibitor use, 4) low inhibition: two functional alleles and some history of CYP2D6 inhibitor use not previously classified, and 5) no inhibition: two functional alleles and no history of CYP2D6 inhibitor use. Women with missing genotype or medication history were excluded from this analysis.
We defined the following set of covariates: time of breast cancer diagnosis, age at diagnosis, Charlson comorbidity index score at diagnosis, menopausal status at diagnosis, county of residence at diagnosis, UICC stage at diagnosis, histological grade, surgery type, and receipt of systemic adjuvant chemotherapy.
Analyses were conducted within strata of the two groups (ER+/TAM+ or ER−/TAM−). We computed the frequency and proportion of case patients and control subjects within categories of assigned protocol of tamoxifen duration, genotype, SSRI use, use of other CYP2D6 inhibitors or substrates, the index of CYP2D6 inhibition, and the covariates. We tested whether the genotypes at the CYP2D6*4 locus observed among control subjects were in Hardy–Weinberg equilibrium by calculating the χ2 test statistic, with expected genotype frequencies based on the observed prevalence of major and minor alleles (52). We performed these calculations in all control subjects combined and within group strata.
We estimated the rate ratio associating CYP2D6 inhibition with breast cancer recurrence as the odds ratio in a conditional logistic regression, including genetic information, use of CYP2D6-inhibiting medications, or the index of CYP2D6 inhibition as the prognostic variable, conditioned on the matched factors. We adjusted for potential confounding in a logistic regression model, which included the marker of CYP2D6 inhibition, time to recurrence or control subject selection, menopausal status, stage, receipt of chemotherapy, receipt of radiation therapy, and type of surgery. We repeated this analysis restricting it to the ER+/TAM+ women who received no chemotherapy. All estimates of association are accompanied by a 95% confidence interval (95% CI) calculated by the profile likelihood method. We computed the two-sided P value for a test of homogeneity of the adjusted odds ratios in the ER+/TAM+ women vs the corresponding adjusted odds ratios in the ER−/TAM− women.
As noted above, prescription data were missing for women diagnosed before 1995. In addition, genotyping results and ER re-assay results were missing for a small proportion of subjects (Figure 1), due to unavailable tumor blocks or indeterminate assay results. Case patients and control subjects with missing data were excluded from analyses that required the variable with a missing result (45).
We examined the concordance of CYP2D6 genotypes in the 106 DNA samples extracted from normal tissue with the genotypes in the DNA samples extracted from the paired tumor tissues. We calculated the crude odds ratio associating failure to confirm the original ER status with recurrence risk. In a subanalysis, we calculated genotype associations with recurrence rate after excluding participants for whom re-assay of the tumor showed ER expression discordant with the assay result at diagnosis. We evaluated the concordance of the information from the DBCG registry with the information from the gold standard medical records. All analyses were performed using SAS version 9 (SAS Institute Inc., Cary, NC).
To address bias from incomplete information on CYP2D6 function, we conducted a quantitative bias analysis (53) informed by classifications based on comprehensive genotyping of the CYP2D6 gene in another study population of German breast cancer patients (54). We used this study to estimate the sensitivity (s = 187/303) and specificity (t = 186/189) of CYP2D6 functional classification based on genotyping only the CYP2D6*4 allele, as well as the positive predictive value (PPVco = 187/190) and negative predictive value (NPVco = 186/302) of reduced CYP2D6 function in control subjects based on genotyping only the CYP2D6*4 allele. We set the prevalence of reduced CYP2D6 function in control subjects (pco = 303/492) equal to the prevalence observed in the German study population (54) based on the comprehensive genotype. For the probabilistic bias analysis, we assigned beta distributions to all of these classification parameters using standard methods (α = numerator + 1, β = total − numerator + 1) (53). The prevalence of reduced CYP2D6 function in case patients (pca) is a function of the prevalence in control subjects and the association between reduced function and breast cancer recurrence (OR):
In the ER−/TAM− strata, we set OR = 1. In the ER+/TAM+ group, OR was the parameter we wished to estimate, so we substituted in its place a beta distribution with α = 1.83 and β = 4.54 scaled to the interval ln|1| to ln|2.77|. This distribution has minimum OR = 1 (null association), maximum OR = 2.77 [the strongest association reported in the study of German breast cancer patients (54)], and mode of OR = 1.22 [the result of a recent meta-analysis (55) of this topic among Caucasians (OR = 1.22, 95% CI = 0.88 to 1.68)]. Finally, we calculated the positive- predictive value in case patients (PPVca) and negative- predictive value in case patients (NPVca) as
where s is sensitivity and t is specificity.
so, when OR is greater than 1.0, NPVca is less than NPVco, which means that a case patient classified as having normal function based only on having no *4 allele is more likely to be reclassified as truly having reduced function than an analogous control subject. With this model and assigned probability distributions, we used Monte Carlo simulation to complete a record-level probabilistic bias analysis following established methods (53). The bias analysis was performed using SAS version 9.
More than 90% of the women in the ER+/TAM+ group were postmenopausal and had stage II or III disease (Table 1), a consequence of the DBCG criteria for assignment to standard tamoxifen protocols during the era when the study population was diagnosed with breast cancer (44). More than 80% of the breast cancer patients underwent mastectomy, which is consistent with treatment patterns in Denmark during this era (56), and a substantial proportion (34%–47%) received radiation therapy (Table 1). As expected, the prevalence of chemotherapy was much higher in the ER−/TAM− group (63% in control subjects and 83% in case patients) than in the ER+/TAM+ group (12% in control subjects and 13% in case patients) (Table 1).
The frequency of the CYP2D6*4 minor allele was 24% in ER+ case patients, 23% in ER− case patients, and 22% among both the ER+/TAM+ control subjects and the ER−/TAM− control subjects. The minor allele frequency among control subjects compares well with the minor allele frequency reported in reference populations (57,58). The CYP2D6*4 allele was in Hardy–Weinberg equilibrium among all control subjects (P = .21 overall, P = .13 in ER+/TAM+ control subjects, and P = .97 in ER−/TAM− control subjects).
In the ER+/TAM+ group, the associations of one functional allele with breast cancer recurrence (adjusted OR = 0.99, 95% CI = 0.76 to 1.3) and no functional allele with breast cancer recurrence (adjusted OR = 1.4, 95% CI = 0.84 to 2.3) were near null (Table 2), as was the association of genetically reduced function (zero or one functional CYP2D6*4 allele) with breast cancer recurrence (adjusted OR = 1.1, 95% CI = 0.81 to 1.4). Similar near-null associations were found for the ER−/TAM− group (adjusted OR = 0.91, 95% CI = 0.62 to 1.3; adjusted OR = 1.3, 95% CI = 0.60 to 2.9; and adjusted OR = 0.95, 95% CI = 0.66 to 1.4, respectively; Table 2). When we restricted the ER+/TAM+ group to women who received no chemotherapy, the associations of one functional allele with breast cancer recurrence (adjusted OR = 1.0, 95% CI = 0.76 to 1.4) and no functional allele with breast cancer recurrence (OR = 1.5, 95% CI = 0.88 to 2.6) were similarly small or null. No estimate of association in the ER+/TAM+ group was statistically– significantly different from the corresponding estimate of association in the ER−/TAM− group (Table 2).
As we reported previously (37), neither SSRI medications nor other CYP2D6 inhibitors increased the risk of breast cancer recurrence (Table 2). When we included information on medications that inhibit CYP2D6 function with the genetic information, the near-null results persisted in all categories of the index of CYP2D6 inhibition and in both groups (Table 2). These associations were less precisely measured because prescription information was not available throughout the study period.
The probabilistic bias analysis revealed that the results were unlikely to have been substantially biased by the absence of genetic information at alleles other than CYP2D6*4, assuming a valid bias model. The estimate of the association between reduced function and breast cancer recurrence in the ER+/TAM+ group changed from 1.1 (95% CI = 0.81 to 1.4) in the conventional analysis reported above to 1.3 (95% simulation interval 0.87 to 1.9) in the bias analysis. In the ER−/TAM− group, the estimate of association changed from 0.95 (95% CI = 0.66 to 1.4) in the conventional analysis reported above to 0.97 (95% CI = 0.71 to 1.3) in the bias analysis. Substituting genetic information from a second cohort of persons without breast cancer (57) yielded comparable results.
DNA extracted from normal tissue vs DNA extracted from paired tumor samples showed perfect concordance between CYP2D6 genotypes (Figure 2, A), so all DNA extractions used for genotyping in this study were from tumor samples. Concordance between positive ER expression at diagnosis and the centralized immunohistochemistry re-assay results (94%, Figure 2, B) was better than the concordance between negative ER expression at diagnosis and the re-assay results (74%, Figure 2, B). Further stratification of these results by case and control status was also informative. In the ER+/TAM+ group, 96% of tumors from control participants expressed the ER when re-assayed, whereas only 92% of tumors from case patients with breast cancer recurrence expressed the ER when re-assayed (Table 1). The crude odds ratio associating failure to express the ER in the re-assay with recurrence equaled 2.0 (95% CI = 1.2 to 3.6). In the ER−/TAM− group, the crude odds ratio associating failure to confirm the absence of ER expression with recurrence was null (OR = 1.0, 95% CI = 0.71 to 1.5). Limiting the genotype analyses to women whose ER expression at re-assay was concordant with ER status at diagnosis did not appreciably change the adjusted estimates of association in the ER+/TAM+ group (one functional allele OR = 1.0, 95% CI = 0.78 to 1.4 and no functional allele OR = 1.5, 95% CI = 0.87 to 2.4) or in the ER−/TAM− group (one functional allele OR = 0.86, 95% CI = 0.55 to 1.4 and no functional allele OR = 1.1, 95% CI = 0.40 to 3.0).
The 30 patients selected for the registry validation substudy were typical of the parent study population, except that they were all treated at Aalborg Hospital. Of the 22 patients categorized as having a recurrence in the DBCG registry, 19 had a recurrence and three had a second primary breast cancer according to the medical record review. One patient had a recurrence that was not registered in the DBCG (Figure 2, C).
All patients who were categorized as receiving tamoxifen by the DBCG had at least 1 year of treatment, and no patient categorized as not receiving tamoxifen by the DBCG had evidence in the medical record of receiving tamoxifen. Most importantly, medical record review showed that most patients assigned to 1- or 2-year treatment protocols at diagnosis actually received tamoxifen treatment for a longer duration (Figure 2, D).
In this population-based study, we observed no substantial association between CYP2D6 inhibition and breast cancer recurrence. The adjusted odds ratio associating no functional CYP2D6*4 allele with breast cancer recurrence in the ER+/TAM+ group (OR = 1.4, 95% CI = 0.84 to 2.3) approximately equaled the corresponding adjusted odds ratio in the ER−/TAM− group (OR = 1.3, 95% CI = 0.60 to 2.9). Taken together, these results should caution against overinterpretation of the former association because no association via the tamoxifen pathway is plausible to explain the latter association. The near-null association persisted regardless of whether we assessed CYP2D6 inhibition only by *4 genotype, by intake of medications that inhibit CYP2D6 function, by combination of *4 genotype with medication history, by limiting the ER+/TAM+ group to women who received no chemotherapy, by limiting the dataset to women with confirmed ER expression, or by analysis of the potential bias arising from misclassification of CYP2D6 function because we genotyped only the *4 allele.
The study design took advantage of the long-standing high-quality clinical database maintained by the DBCG (42) and the tissue archives maintained by the pathology departments of Danish hospitals (59). Linking these resources yielded a study effectively immune to selection bias and, to our knowledge, the largest study to date of the association between CYP2D6 inhibition and breast cancer recurrence. The 184 recurrent cases with at least one *4 allele in the ER+/TAM+ group represent approximately one-third of the accumulated cases with genetically identified reduced function reported in studies thus far.
Data quality was high, with perfect concordance between genotypes in DNA extracted from tumor tissue and lymph node tissue in the validation subset and with alleles at the CYP2D6*4 locus in Hardy–Weinberg equilibrium among all control subjects and consistent with the expected allele frequencies. Concordance of re-assayed ER expression with ER status at diagnosis was consistent with earlier reports (60,61) and showed the expected associations with recurrence risk. The association between failure to express the ER at re-assay and recurrence risk in the ER+/TAM+ group provides the study design with face validity and demonstrates the study’s ability to detect associations between rare events and recurrence risk.
All recurrences in the validation subset were confirmed by medical record review, which is consistent with an earlier and much larger validation study that reported a positive predictive value of 99.4% for breast cancer recurrence recorded by the DBCG (62). We observed perfect agreement for all of the matched factors except for one patient’s menopausal status. Only duration of tamoxifen therapy differed frequently between the registry data and medical records, with many patients who were originally assigned to tamoxifen therapy of short duration receiving tamoxifen for longer durations. This discrepancy likely resulted from modifications to the tamoxifen protocol occurring after the patient’s diagnosis and initial assignment to tamoxifen durations of 1 or 2 years. The discrepancy actually strengthens our study’s results because it suggests that the null association did not arise from some patients receiving short-duration tamoxifen therapy.
The major limitation of this study was absence of genotyping data for CYP2D6 alleles other than *4. Complete genotyping of the DNA samples extracted from paraffin-embedded tissues required more resources than were available to us. To address this limitation, we implemented a quantitative bias analysis, which assumed that case patients with recurrence were more likely to carry alleles with reduced function than were control subjects. All of the parameters of the bias model were informed by published external data sources. Assuming an accurate bias model, this bias analysis showed that the near-null results would have changed very little with complete genotyping data. A second limitation was the absence of information on tamoxifen adherence by case patients and control subjects. About half of tamoxifen-treated patients fail to complete the intended duration of their tamoxifen therapy (63). Among the 20 ER+/TAM+ patients included in our medical record review, six did not complete their intended duration of tamoxifen therapy, two because of recurrent breast cancer. Tamoxifen adherence is related to recurrence risk (64) and may be predicted by CYP2D6 genotype (65), in which case, adherence to the intended duration would be a causal intermediate between CYP2D6 genotype and recurrence. Results adjusted for adherence would be more biased than without adjustment, usually toward the null (66).
Earlier studies of the association between CYP2D6 inhibition and risk of breast cancer recurrence or mortality in tamoxifen-treated women have reported widely heterogeneous results (55), and there is no adequate explanation for this heterogeneity (18). Although the heterogeneity presents an important barrier in interpreting a quantitative synthesis (18,67), a recent meta-analysis concluded that the effect of CYP2D6 inhibition on recurrence risk may be relatively small (55). This conclusion is consistent with our results, with reasonable bounds that can be placed on the expected association (18), and with the well-understood pharmacology of tamoxifen therapy (18,68), which requires that tamoxifen and its metabolites overwhelm estrogen in competition for binding to the ER. It is common for genetic research to initially show strong associations, sometimes in more than one study (69), only to find that the strength of association diminishes or disappears as evidence accumulates (70–72). Despite early enthusiasm for the potential to identify poor candidates for tamoxifen therapy by CYP2D6 genotyping (15), the true association between CYP2D6 inhibition and breast cancer recurrence risk in tamoxifen-treated patients is likely to be null or small.
US National Cancer Institute at the National Institutes of Health (R01 CA118708 to T.L.L.); Danish Cancer Society (DP06117 to S.H.-D.); Karen Elise Jensen Foundation (H.T.S.).
The authors thank the Danish Breast Cancer Cooperative Group for access to its registry data and for preparing the initial dataset.
The granting agencies had no role in the design of the study; the collection, analysis, and interpretation of the data; the writing of the article; or the decision to submit the article for publication.