Non-CYP2D6 Variants Selected by a GWAS Improve the Prediction of Impaired Tamoxifen Metabolism in Patients with Breast Cancer

A Genome-Wide Association Study of Endoxifen Serum Concentrations and Adjuvant Tamoxifen Efficacy in Early-Stage Breast Cancer Patients

Tamoxifen is part of the standard of care of endocrine therapy for adjuvant treatment of breast cancer. However, survival outcomes with tamoxifen are highly variable. The concentration of endoxifen, the 30-100 times more potent metabolite of tamoxifen and bioactivated by the CYP2D6 enzyme, has been described as the most relevant metabolite of tamoxifen metabolism. A genome-wide association study (GWAS) was performed with the objective to identify genetic polymorphisms associated with endoxifen serum concentration levels and clinical outcome in early-stage breast cancer patients receiving tamoxifen. A GWAS was conducted in 608 women of the CYPTAM study (NTR1509/PMID: 30120701). Germline DNA and clinical and survival characteristics were readily available. Genotyping was performed on Infinium Global Screening Array (686,082 markers) and single nucleotide polymorphism (SNP) imputation by using 1000 Genomes. Relapse-free survival during tamoxifen (RFSt) was defined the primary clinical outcome. Endoxifen serum concentration was analyzed as a continuous variable. Several genetic variants reached genome-wide significance (P value: ≤5 × 10-8). Endoxifen concentrations analysis identified 430 variants, located in TCF20 and WBP2NL genes (chromosome 22), which are in strong linkage disequilibrium with CYP2D6 variants. In the RFSt analysis, several SNP were identified (LPP gene: rs77693286, HR 18.3, 95% CI: 15.2-21.1; rs6790761, OR 18.2, 95% CI: 15.5-21.1). Endoxifen concentrations have a strong association with the chromosome 22, which contains the CYP2D6 gene.

Novel Independent Trans- and Cis-Genetic Variants Associated with CYP2D6 Expression and Activity in Human Livers

Cytochrome P450 2D6 (CYP2D6) is a critical hepatic drug-metabolizing enzyme in humans, responsible for metabolizing approximately 20%-25% of commonly used medications such as codeine, desipramine, fluvoxamine, paroxetine, and tamoxifen. The CYP2D6 gene is highly polymorphic, resulting in substantial interindividual variability in its catalytic function and the pharmacokinetics and therapeutic outcomes of its substrate drugs. Although many functional CYP2D6 variants have been discovered and validated, a significant portion of the variability in the expression and activity of CYP2D6 remains unexplained. In this study, we performed a genome-wide association study (GWAS) to identify novel variants associated with CYP2D6 protein expression in individual human livers, followed by a conditional analysis to control for the effect of functional CYP2D6 star alleles. We also examined their impact on hepatic CYP2D6 activity. Genotyping on a genome-wide scale was achieved using the Illumina Multi-Ethnic Genotyping Array (MEGA). A data-independent acquisition (DIA)-based proteomics method was used to quantify CYP2D6 protein concentrations. CYP2D6 activity was determined by measuring the dextromethorphan O-demethylation in individual human liver s9 fractions. The GWAS identified 44 single nuclear polymorphisms (SNPs) that are significantly associated with CYP2D6 protein expressions with a P value threshold of 5.0 × 10-7 After the conditional analysis, five SNPs, including the cis-variants rs1807493 and rs1062753 and the trans-variants rs4073010, rs729559, and rs80274432, emerged as independent variants significantly correlated with hepatic CYP2D6 protein expressions. Notably, four of these SNPs, except for rs80274432, also exhibited a significant association with CYP2D6 activities in human livers, suggesting their potential as novel and independent cis- and trans-variants regulating CYP2D6. SIGNIFICANT STATEMENT: Using individual human livers, we identified four novel cis- and trans-pQTLs/aQTLs (protein quantitative trait loci/activity quantitative trait loci) of Cytochrome P450 2D6 (CYP2D6) that are independent from known functional CYP2D6 star alleles. This study connects the CYP2D6 gene expression and activity, enhancing our understanding of the genetic variants associated with CYP2D6 protein expression and activity, potentially advancing our insight into the interindividual variability in CYP2D6 substrate medication response.

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 R² ) 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.

GWAS Links New Variant in Long Non-Coding RNA LINC02006 with Colorectal Cancer Susceptibility

Simple Summary

Identifying risk factors for cancer development can allow for appropriate stratification and surveillance of individuals at risk, increasing their chances of benefiting from early disease detection; however, most of the genetic factors contributing to the risk of colorectal cancer (CRC) remain undetermined. Here, we adopted a new approach for selecting index polymorphism for further validation in combination with a genome-wide association study of pooled DNA samples for CRC susceptibility variants in the Polish population. This study, including 2013 patients and controls, uncovered five susceptibility loci not previously reported for CRC. Four of identified variants were located within genes likely involved in tumor invasiveness and metastasis, suggesting that they could be markers of poor prognosis in CRC patients. Our results provide evidence that conducting association studies on small but homogenous populations can help us discover new common risk variants specific to the studied population.

Abstract

Despite great efforts, most of the genetic factors contributing to the risk of colorectal cancer (CRC) remain undetermined. Including small but homogenous populations in genome-wide association studies (GWAS) can help us discover new common risk variants specific to the studied population. In this study, including 465 CRC patients and 1548 controls, a pooled DNA samples-based GWAS was conducted in search of genetic variants associated with CRC in a Polish population. Combined with a new method of selecting single-nucleotide polymorphisms (SNPs) for verification in individual DNA samples, this approach allowed the detection of five new susceptibility loci not previously reported for CRC. The discovered loci were found to explain 10% of the overall risk of developing CRC. The strongest association was observed for rs10935945 in long non-coding RNA LINC02006 (3q25.2). Three other SNPs were also located within genes (rs17575184 in NEGR1, rs11060839 in PIWIL1, rs12935896 in BCAS3), while one was intergenic (rs9927668 at 16p13.2). An expression quantitative trait locus (eQTL) bioinformatic analysis suggested that these polymorphisms may affect transcription factor binding sites. In conclusion, four of the identified variants were located within genes likely involved in tumor invasiveness and metastasis. Therefore, they could possibly be markers of poor prognosis in CRC patients.

Recent advances in the personalized treatment of estrogen receptor-positive breast cancer with tamoxifen: a focus on pharmacogenomics

Introduction: Tamoxifen is still an important drug in hormone-dependent breast cancer therapy. Personalization of its clinical use beyond hormone receptor positivity could improve the substantial variability of the treatment response.Areas covered: The overview of the current evidence for the treatment personalization using therapeutic drug monitoring, or using genetic biomarkers including CYP2D6 is provided. Although many studies focused on the PK aspects or the impact of CYP2D6 variability the translation into clinical routine is not clearly defined due to the inconsistent clinical outcome data.Expert opinion: We believe that at least the main candidate factors, i.e. CYP2D6 polymorphism, CYP2D6 inhibition, endoxifen serum levels may become important predictors of clinical relevance for tamoxifen treatment personalization in the future. To achieve this aim, however, further research should take into consideration more precise characterization of the disease, epigenetic factors and also utilize an appropriately powered multifactorial approach instead of a single gene evaluating studies.

Pharmacogenomics of breast cancer: highlighting CYP2D6 and tamoxifen

PURPOSE

To review recent pharmacogenomics studies on breast cancer patients undergoing tamoxifen therapy, highlighting how our knowledge on cytochrome P450 2D6 (CYP2D6) can help to guide the development of adjuvant therapies for these patients.

METHODS

A comprehensive literature search was conducted. Articles reporting findings pertaining to the effect of CYP2D6 on the therapeutic efficacy of tamoxifen, those reporting how targeting CYP2D6 could inform tamoxifen-based therapy development, and those on the tamoxifen effects on cell lines and animal models were included in the review.

RESULTS

With CYP2D6 being the primary enzyme for tamoxifen metabolism, single-nucleotide polymorphisms (SNPs) in this gene were one of the determinants in the rate of tamoxifen metabolism, thereby potentially having an effect on the efficacy of tamoxifen-based therapies. Our review indicates the potential effectiveness of targeting these SNPs, including those for the CYP2D6*10 allele (c. 100C > T), in modifying the level of tamoxifen metabolism. These findings suggest the importance of pharmacogenomics research in our understanding of the efficacy of adjuvant therapies. However, the involvement of multiple enzymes in tamoxifen metabolism, dietary factors, ethnic differences in gene frequencies, and patients' compliance to tamoxifen therapies in studies do present challenges in pharmacogenomics research.

CONCLUSIONS

Pharmacogenomics could play important roles in mediating the advancement of tamoxifen-based adjuvant therapies. Research efforts should be directed towards the exploration of further SNPs of CYP2D6 that affect tamoxifen metabolism, as well as epigenetic changes in CYP2D6, enabling the design of precision medicine and confirming clinical validity in the use of pharmacogenomics for tamoxifen.