Cross-Ancestry Genome-Wide Association Study Defines the Extended CYP2D6 Locus as the Principal Genetic Determinant of Endoxifen Plasma Concentrations.

The therapeutic efficacy of tamoxifen is predominantly mediated by its active metabolites 4-hydroxy-tamoxifen and endoxifen, whose formation is catalyzed by the polymorphic cytochrome P450 2D6 (CYP2D6). Yet, known CYP2D6 polymorphisms only partially determine metabolite concentrations in vivo. We performed the first cross-ancestry genome-wide association study with well-characterized patients of European, Middle-Eastern, and Asian descent (n = 497) to identify genetic factors impacting active and parent metabolite formation. Genome-wide significant variants were functionally evaluated in an independent liver cohort (n = 149) and in silico. Metabolite prediction models were validated in two independent European breast cancer cohorts (n = 287, n = 189). Within a single 1-megabase (Mb) region of chromosome 22q13 encompassing the CYP2D6 gene, 589 variants were significantly associated with tamoxifen metabolite concentrations, particularly endoxifen and metabolic ratio (MR) endoxifen/N-desmethyltamoxifen (minimal P = 5.4E-35 and 2.5E-65, respectively). Previously suggested other loci were not confirmed. Functional analyses revealed 66% of associated, mostly intergenic variants to be significantly correlated with hepatic CYP2D6 activity or expression (ρ = 0.35 to -0.52), and six hotspot regions in the extended 22q13 locus impacting gene regulatory function. Machine learning models based on hotspot variants (n = 12) plus CYP2D6 activity score (AS) increased the explained variability (~ 9%) compared with AS alone, explaining up to 49% (median R2 ) and 72% of the variability in endoxifen and MR endoxifen/N-desmethyltamoxifen, respectively. Our findings suggest that the extended CYP2D6 locus at 22q13 is the principal genetic determinant of endoxifen plasma concentration. Long-distance haplotypes connecting CYP2D6 with adjacent regulatory sites and nongenetic factors may account for the unexplained portion of variability.

Loss of heterozygosity at the CYP2D6 locus in breast cancer: implications for germline pharmacogenetic studies.

BACKGROUND: Controversy exists regarding the impact of CYP2D6 genotype on tamoxifen responsiveness. We examined loss of heterozygosity (LOH) at the CYP2D6 locus and determined its impact on genotyping error when tumor tissue is used as a DNA source. METHODS: Genomic tumor data from the adjuvant and metastatic settings (The Cancer Genome Atlas [TCGA] and Foundation Medicine [FM]) were analyzed to characterize the impact of CYP2D6 copy number alterations (CNAs) and LOH on Hardy Weinberg equilibrium (HWE). Additionally, we analyzed CYP2D6 *4 genotype from formalin-fixed paraffin-embedded (FFPE) tumor blocks containing nonmalignant tissue and buccal (germline) samples from patients on the North Central Cancer Treatment Group (NCCTG) 89-30-52 tamoxifen trial. All statistical tests were two-sided. RESULTS: In TCGA samples (n =627), the CYP2D6 LOH rate was similar in estrogen receptor (ER)-positive (41.2%) and ER-negative (35.2%) but lower in HER2-positive tumors (15.1%) (P < .001). In FM ER+ samples (n = 290), similar LOH rates were observed (40.8%). In 190 NCCTG samples, the agreement between CYP2D6 genotypes derived from FFPE tumors and FFPE tumors containing nonmalignant tissue was moderate (weighted Kappa = 0.74; 95% CI = 0.63 to 0.84). Comparing CYP2D6 genotypes derived from buccal cells to FFPE tumor DNA, CYP2D6*4 genotype was discordant in six of 31(19.4%). In contrast, there was no disagreement between CYP2D6 genotypes derived from buccal cells with FFPE tumors containing nonmalignant tissue. CONCLUSIONS: LOH at the CYP2D6 locus is common in breast cancer, resulting in potential misclassification of germline CYP2D6 genotypes. Tumor DNA should not be used to determine germline CYP2D6 genotype without sensitive techniques to detect low frequency alleles and quality control procedures appropriate for somatic DNA.