Carla Miranda1
Macarena Galleguillos1Roberto Torres2Karla Tardón2
Dante D. Cáceres3
Kuen Lee1,4
María A. Redal5,6
Nelson M. Varela1,6*
Luis A. Quiñones1,6*- 1Laboratory of Chemical Carcinogenesis and Pharmacogenetics, Department of Basic-Clinical Oncology (DOBC), Faculty of Medicine, University of Chile, Santiago, Chile
- 2National Cancer Institute, Santiago, Chile
- 3Institute of Population Health, Faculty of Medicine, University of Chile, Santiago, Chile
- 4Faculty of Medicine, University of Chile, Santiago, Chile
- 5Genetic Division, Department of Medicine, Hospital de Clínicas José de San Martín, Buenos Aires, Argentina
- 6Latin American Network for Implementation and Validation of Clinical Pharmacogenomics Guidelines (RELIVAF-CYTED), Madrid, Spain
Tamoxifen (TAM), a selective oestrogen receptor modulator, is one of the most used treatments in oestrogen receptor-positive (ER+) early and metastatic breast cancer (BC) patients. The response to TAM has a high degree of inter-individual variability. This is mainly due to genetic variants in CYP2D6 gene, as well as other genes encoding proteins involved in the TAM pharmacokinetic and/or pharmacodynamic. Therefore, prediction of the TAM response using these genetic factors together with other non-genetic variables may be relevant to improve breast cancer treatment. Thus, in this work, we used genetic polymorphisms and clinical variables for TAM response modelling. One hundred sixty-two ER + BC patients with 2 years of TAM treatment were retrospectively recruited, and the genetic polymorphisms CYP2D6*4, CYP3A4*1B (CYP3A4*1.001), CYP3A5*3, UGT2B7*2, UGT2B15*2, SULT1A1*2, and ESRA V364E were analyzed by PCR-RFLP. Concomitantly, the therapeutic response was obtained from clinical records for association with genotypes using univariate and multivariate biostatistical models. Our results show that UGT2B15*1/*2 genotype protects against relapse (OR = 0.09; p = 0.02), CYP3A5*3/*3 genotype avoids endometrial hyperplasia (OR = 0.07; p = 0.01), SULT1A1*1/*2 genotype avoids vaginal bleeding (OR = 0.09; p = 0.03) and ESRA 364E/364E genotype increases the probability of vaginal bleeding (OR = 5.68; p = 0.02). Logistic regression models, including genomic and non-genomic variables, allowed us to obtain preliminary predictive models to explain relapse (p = 0.010), endometrial hyperplasia (p = 0.002) and vaginal bleeding (p = 0.014). Our results suggest that the response to TAM treatment in ER + BC patients might be associated with the presence of the studied genetic variants in UGT2B15, CYP3A5, SULT1A1 and ESRA genes. After clinical validation protocols, these models might be used to help to predict a percentage of BC relapse and adverse reactions, improving the individual response to TAM-based treatment.
1 Introduction
Breast cancer (BC) is the second leading cause of cancer death in women (American Cancer Society, 2017) and is caused by many factors, such as age, breast density, heredity, and exposure to oestrogens (Kumar et al., 2010). Breast cancer (BC) patients with oestrogen receptor positivity (ER+) are treated in addition to surgery, chemotherapy, and/or radiation therapy with hormone therapy (Maughan et al., 2010; (McDonnell and Wardell 2010; NCI 2012). Tamoxifen (TAM), a type of hormonal therapy known as a selective oestrogen receptor modulator, behaves as an oestrogen antagonist in the breast and as an agonist in the endometrium (Lewis and Jordan 2005; Terán and Teppa 2005; Hertz et al., 2009; NCI 2012), and despite the new hormonal agents, TAM remains the drug of choice for any BC stage (Davies et al., 2011) since it increases the disease-free time in pre- or postmenopausal women, reducing the annual death rate by 34% (MINSAL, 2005; Untch and Thomssen, 2010). However, the response among patients is variable in terms of efficacy and side effects.
TAM undergoes extensive biotransformation in the liver by phase I (CYP) and phase II enzymes. Phase I generates, among other metabolites, N-desmethyltamoxifen (N-desmethylTAM), 4-hydroxy tamoxifen (4-hydroxyTAM), and 4-hydroxy-N-desmethyltamoxifen (endoxifen) (Brauch et al., 2009). Endoxifen and 4-hydroxyTAM are the main and unique active metabolites of TAM due to their high affinity to ERs and their activity being 30- to 100-fold greater than TAM itself (Wakeling and Slater 1980; Johnson et al., 2004; Fernández-Santander et al., 2007; Brauch et al., 2009). In vitro studies have determined that these metabolites are more effective in reducing cell proliferation (Jordan et al., 1977; Clarke et al., 2003; Borges et al., 2006). The two most abundant metabolites in plasma are N-desmethylTAM and endoxifen (Clarke et al., 2003). Subsequently, these metabolites undergo biotransformation by phase II enzymes, the main enzymes involved being SULT1A1, UGT2B7, and UGT2B15 (Nishiyama et al., 2002). TAM has a long half-life, and steady-state concentrations are obtained at 4 weeks for TAM and at 8 weeks for N-desmethylTAM (Buckley and Goa, 1989; Lien et al., 1992; Kisanga et al., 2004).
Differences have been found between the response to TAM and the histological type, degree of differentiation of the breast tumor cells, age, and state of menopause in the patient. TAM has an overall response rate of almost 75% in patients with ER+ and PR + BC. In the case of adjuvant treatment, TAM reduces the risk of relapse by 25% and mortality by 17% (Fisher et al., 1998; Davies et al., 2011). On the other hand, the long-term safety of TAM is well elucidated, and possible harmful effects on bone metabolism, endometrial cancer, thromboembolic diseases, and cognitive disorders due to chronic oestrogen deprivation in brain tissue have been described (Fisher et al., 1998; Cuzick et al., 2002; Cuzick, 2005). A 1% increase in the incidence of endometrial cancer has been observed due to its oestrogen-inducing effect in the uterus and thromboembolic events in 1–2% of cases.
Despite several studies having been conducted, after 3 decades, there are still differences in the TAM response that have not been explained. It is known that these variations in the response to drugs are multifactorial, the result of the interaction between multiple and complex genetic, physiological and environmental factors (Ma and Lu 2011). A potential explanation is the presence of genetic variants in the genes encoding biotransformation enzymes affecting their efficacy and safety (Zhou et al., 2008). In this respect, undoubtedly genetic variation in the CYP2D6 driving to ultrarapid (XN), intermediate metabolizer or poor metabolizer phenotypes (e.g. *3, *4, *5, *9, *10, *17 haplotypes) may affects plasma concentrations of TAM and its metabolites, as CPIC has stablished in the updated guidelines (Goetz et al, 2018). CYP2D6*4 (rs3892097) is associated with decreased side effects, such as hot flashes, disease-free time, and survival rate (Goetz et al., 2005; Jin et al., 2005; Schroth et al., 2007; Bijl et al., 2009; Ramón y Cajal et al., 2010; Rae et al., 2012), and decreased serum levels of TAM metabolites (Gjerde et al., 2007; Kiyotani et al., 2010; Lim et al., 2011; Madlensky et al., 2011). In addition, CYP3A4*1B (currently named CYP3A4*1.001, rs2740574) is associated with a 3-fold higher risk of endometrial cancer in TAM-treated BC patients (Chu et al., 2007). CYP3A5*3 (rs776746) is associated with characteristics of the tumor in BC postmenopausal women treated with TAM (Tucker et al., 2005). On the other hand, CYP2C9 *2 and *3 variants reduce the concentration of TAM metabolites, affecting the therapeutic response (Mürdter et al., 2011).
Phase II enzymes also affect TAM response. SULT1A1*2 (rs9282861) increases the risk of relapse in TAM-treated BC patients (Wegman et al., 2005). UGT2B15*2 (rs1902023) decreases the risk of relapse of BC (Nowell et al., 2005). Furthermore, the presence of UGT2B15*2 and SULT1A1*2 is associated with a lower risk of relapse and a significant reduction in survival time in patients with BC treated with TAM (Nowell et al., 2005).
Finally, several mutations have been identified in the oestrogen receptor α gene (ESRA) in BC patients (Herynk et al., 2004); however, their effect on the efficacy and safety of TAM treatment has not been elucidated. In this respect, using SIFT and PolyPhen bioinformatic tools, it was established that the single nucleotide polymorphism (SNP) ESRA (OMIM* 133430) V364E (rs121913044; 1453T>A; Val364Glu) generates a deleterious change with possible damage to the N-terminus of the hormone-binding domain in the ER, giving rise to a 40 times lower affinity for oestrogen (Flicek et al., 2011). This SNP also exhibits a dependence on oestrogen for binding to an ERE, although it maintains its negative dominant activity entirely. Therefore, ESRA V364E is highly active and capable of repressing ER-mediated transcription, both when ESRA V364E and normal ESRA proteins are present together in cells and even without DNA binding (Wrenn and Katzenellenbogen 1993; Ince et al., 1995; McInerney et al., 1996).
Although in the last decade several studies have investigated genetic variants in TAM-metabolizing enzymes that might determine the differences in the response to treatment in patients with BC, there are still controversies about their relationship with the response (Brauch and Schwab, 2014; Brewer et al., 2014; Province et al., 2014; Binkhorst et al., 2015). To address this question, we aimed to study the association among response (relapse and ADR) and 7 genetic variants in genes encoding proteins involved in the pharmacokinetics and pharmacodynamics of TAM (CYP2D6∗4, CYP3A4∗1B (CYP3A4∗1.001), CYP3A5∗3, SULT1A1∗2, UGT2B7∗2, UGT2B15∗2 and ESRA V364E) in women with hormone-dependent BC and under adjuvant treatment with TAM. The main goal is to generate predictive models to approximate the response of patients according to their genetic-metabolic characteristics.
2 Patients and Methods
2.1 Study Design
A retrospective case-control study was carried out from January 2014 to January 2015 at the Polyclinic of Oncology of the National Institute of Cancer (INC). The sample size was determined according to the frequency of carriers with the variant allele carriers in the population under study, using PS Power and Sample Size Calculations Version 3.0, January 2009, considering an 80% power; α = 0.05; OR = 2.0; and the less frequent CYP3A4∗1B (CYP3A4∗1.001, rs2740574), according to the literature (Fuentes and Silveyra, 2019).
2.2 Patients
One hundred sixty-two (162) patients were enrolled. The inclusion criteria were as follows: 1) patients with histologically confirmed BC, 2) in adjuvant treatment with 20 mg of tamoxifen (Novadex®) daily for at least 24 months; 2) > 18 years-old, 3) positive ER and negative HER-2 status, 4) cell differentiation degree I to III. The exclusion criteria were as follows: 1) patients diagnosed with in situ BC; 2) negative ER, negative PR, triple-negative or Her-2 positive patients; 3) concomitant treatments with vitamin K antagonists (warfarin, acenocoumarol or dicoumarol), serotonin reuptake antagonist antidepressants (fluoxetine or paroxetine), mitomycin, nitonavir, primidone, or chemotherapy (fluorouracil, methotrexate or cyclophosphamide), 4) chronic unbalanced systemic pathology or other active cancers and, 5) aromatase inhibitor treatment or LHRH agonist. The events (relapses and ADR) were evaluated after 6 months of TAM treatment. All the patients signed a written consent form and agreed to be included in this study.
Toxicity was assessed based on the Common Terminology Criteria for Adverse Events, v4.03. Vaginal bleeding was defined as any vaginal bleeding in postmenopausal patients (more than 1 year without menstruation or oophorectomized) or any abnormal bleeding, such as unexpected bleeding or in quantity or habitual bleeding in premenopausal patients. Endometrial hyperplasia was defined as an endometrial thickness greater than 5 mm detected by transvaginal echotomography in postmenopausal patients or an abnormal thickness of the endometrium according to the phase of the menstrual cycle in premenopausal patients.
The ethnicity of the study group was approached by the percentage of Amerindian-Caucasian admixture (%Ma-c) based on the ABO system (Valenzuela et al., 1987; Acuña et al., 2000).
The patients were gathered into two groups: the control group include patients who do not present the analyzed ADRs (relapse, hot flashes, etc.), and the cases group include patients who have the ADRs.
2.3 Genotyping Analyses
Potentially functional SNPs encoding the proteins related to the TAM response were obtained from the PharmGKB database (Whirl-Carrillo et al., 2012), the NCBI dbSNP database (https://www.ncbi.nlm.nih.gov/snp/), the SNPinfo Web Server (https://snpinfo.niehs.nih.gov), the Ensembl® genome database project (https://www.ensembl.org/index.html) and the level of evidence for each SNP (Supplementary Table S1). Genomic DNA was isolated from the peripheral blood samples of the subjects using a Genomic DNA Extraction Blood DNA Kit (FAVORGEN®, BIOTECH CORP, Headquarters, Taiwan, China) and from buccal mucosa cells using a MasterAmp™ Buccal Swab Kit (Epicentre®, an Illumina company, Madison, USA). DNA samples were quantified at 260/280 nm using a Nanodrop spectrophotometer (model DS-11, FX Series Spectrophotometer, USA). CYP2D6∗4 (rs3892097), CYP3A4∗1B (CYP3A4∗1.001, rs2740574), CYP3A5∗3 (rs776746), SULT1A1∗2 (rs9282861), UGT2B7∗2 (rs7439366), UGT2B15∗2 (rs1902023), and ESRA V364E (rs121913044) were analyzed using polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) (Cavalli et al., 2001; Schur et al., 2001; Hajdinjak and Zagradišnik 2004; Lee et al., 2005; Kagaya et al., 2007; Arslan et al., 2011). The primers and restriction enzymes used are presented in Supplementary Table S2. Each assay contained four controls: one sample for each genotype as a positive control and one negative sample with nuclease-free pure water to volume. For quality assurance purposes, we randomly chose 20% of the samples for repetition of the analysis. When analyses were not coincident, we excluded the samples.
2.4 Statistical Analyses
Statistical analyses were performed using GraphPad Prism 4.0 and STATA 11.1, with p < 0.05 considered statistically significant. The results are expressed as the mean ± standard deviation (SD), number, percentage, or frequency where appropriate. The Shapiro-Wilk test was used to determine quantitative variable distributions. Mean values between the two groups were compared using the unpaired t-test for variables with a normal distribution, while the Mann-Whitney U test was used for variables with non-normal distributions. Frequencies of qualitative variables were also calculated. The chi-square test with Fischer’s exact test was used to investigate differences in genotypic and allelic frequencies between the two groups.
We checked for Hardy-Weinberg equilibrium (HWE) in our sample even though the conditions for HWE are not applicable because it is not a random sampling in a random-mating population, a control population or the general population (Namipashaki et al., 2015) and is a group with a selection bias based on the disease (i.e., SNPs can also be related to cancer). The polymorphic variants CYPD6∗4, CYP3A4∗1B (CYP3A4∗1.001), CYP3A5∗3, SULT1A1∗2, and accomplish HWE, while UGT2B15∗2 and ESRA V364E do not (cut-off chi2: 3.84).
The polymorphisms were evaluated using co-dominant (wild type vs heterozygote vs variant), dominant (wild type vs. heterozygote/variant), and recessive (wild type/heterozygote vs. variant) inheritance models. Bi-variable analyses were performed to determine the association between the events (relapses and ADR), considered the dependent variables, and the characteristics of the patients, considered the independent variables. Thus, the independent variables considered were BMI, blood group, smoking habit, socioeconomic status, treatment with oral contraceptives, hormone replacement therapies, age at menarche, menopausal status, family history of cancer, age at diagnosis, cancer stage, and presence of genetic variants (exposure variable).
Multivariate logistic regression models and multivariate linear regression analyses were employed to investigate associations between genetic or non-genetic characteristics and relapse or ADR. The logistic multivariate models were adjusted stepwise using a forward and backward procedure to include potentially relevant variables to derive statistical association models characterized by pseudo R2. All association studies were assayed by choosing parameters with a better statistical association for each analysis. The odds ratio (OR), or coefficients (Coef.), and 95% confidence intervals (CIs) are reported in the multivariate logistic regression models.
2.5 Ethics
The study was carried out under strict ethical procedures recommended by the Ethics Committee of the Faculty of Medicine of the University of Chile (July 24, 2013) and the Research Ethics Committee of the Northern Metropolitan Health Service (June 3, 2013), in accordance with the procedures suggested in the Declaration of Helsinki (WMA, 2013) and according to Chilean Laws 20.120, 20.584, and 19.628 and the guidelines of Good Clinical Practices. The treatment schedule for patients was according to the cancer stage (I-III), histology, cell differentiation degree, and presence and relative abundance of 3 selected differentially abundant proteins (ER+, PR+, Her2-), which involved surgery followed by radiotherapy and/or chemotherapy according to the Breast Cancer Clinical Guide, 2nd Ed. Santiago, Chile.
3 Results
The baseline characteristics of the patients (58 ± 13 years; BMI: 29 ± 6) are shown in Table 1, and the characteristics of the patients according to genotype are shown in Supplementary Table S3. Sociogenetic gradients, risk factors, gynecological characteristics, and pathological features for descriptive analyses can also be observed. Of the total patients, 68 were premenopausal, and 94 were postmenopausal. In relation to the cancer stage at diagnosis, the majority of patients were in stage II (51.23%), the predominant histology type was invasive ductal carcinoma (86.88%), and the cell differentiation degree was predominantly G2 (54.17%). A 49.3% Ma-c was found.
TABLE 1. Characteristics of patients (n = 162).
Table 2 shows the therapeutic response characteristics of the patients. Relapse was found in 9 patients (5.73%). The most severely observed ADRs among patients were endometrial hyperplasia (9.20%) and vaginal bleeding (4.91%), and the most frequent were hot flashes (61.96%), bone pain (13.50%), and cramps (12.27%).
TABLE 2. Relapse and Adverse drug reactions (ADRs) in patients.
The genotypic and allelic frequencies for the analyzed polymorphisms are shown in Table 3.
TABLE 3. Genotype and allele frequencies of the studied polymorphisms in patients.
We performed univariable logistic regression of the risk of relapse (Supplementary Tables S4, S5), endometrial hyperplasia (Supplementary Tables S6, S7), and vaginal bleeding (Supplementary Tables S8, S9) in association with genetic and non-genetic variables. In Table 4, only statistically significant results for the univariable logistic regression analysis (risk of relapse, endometrial hyperplasia, and vaginal bleeding) are shown. The results show that UGT2B15∗2 A/C (rs1902023) in a dominant model of inheritance was associated with relapse; CYP3A5∗3 A/G (rs776746) in a recessive model of inheritance was associated with endometrial hyperplasia; and SULT1A1∗2 G/A (rs9282861) in a dominant model of inheritance and ESRA V364E (rs121913044) in a recessive model of inheritance were associated with vaginal bleeding.
TABLE 4. Univariable logistic regression analysis of risk of severe ADRs according to genotypes.
After using stepwise forward and backward procedures, multivariate logistic regression analyses for the risk of relapse and ADRs, including genetic and non-genetic variables, were performed. We obtained significant models for endometrial hyperplasia, vaginal bleeding, phlebitis, headache, nausea, hot flash, cramps, bone pain, and urticaria (Figure 1). Supplementary Tables S10–S12 show significant multivariate logistic regression analyses or logit models for relapse (Supplementary Table S10), endometrial hyperplasia (Supplementary Table S11) and vaginal bleeding (Supplementary Table S12) with Pseudo R2 = 0.495 (p = 0.01), Pseudo R2 = 0.42 (p = 0.002) and Pseudo R2 = 0.34 (p = 0.014), respectively.
FIGURE 1. Research Scheme and obtained multivariate predictive models. ADRs, adverse drug reactions; BC: Breast Cancer; SULT: Sulfotransferase; UGT: UDP-glucuronosyl transferase; CYP: Cytochrome P-450; ESRA: Estrogen Receptor 1; SE: Socioeconomic status (per capita income), I: <$135.000 CLP (U$ 200); II: $135.001-$500.000 CLP (U$ 200–750); III: $500.001-$1.000.000 CLP (U$ >750–1,450); IV: >2.000.000 CLP >$CLP 2,000,000 (U$ >2,900); HBC: Family History of Breast Cancer; HBOC: Family History of breast or ovary cancer; M: Menarche age; G: Number of Gestations; A: Number of Abortions; P: Number of deliveries; L: Breastfeeding time; Age diag: Age of diagnosis, years; OCT: Oral Contraceptive Treatment; HRT: Treatment with HRT for menopause; Menop: postmenopause; S: Cancer stage at diagnosis (Stage I-IV); DCis: Ductal carcinoma in situ; IDC: Invasive Ductal Carcinoma; ILC: Invasive Lobular Carcinoma; OH: Others Histology (OH).
Associations among genetic variants with any adverse effect versus no adverse effect or with some grades of adverse effects were also assessed, but no statistical significance was observed (Supplementary Tables S13, S14).
4 Discussion
Patient response to TAM has been investigated for a long time. In this respect, differences of approximately 25–50% in patients are observed (Fisher et al., 1998; Davies et al., 2011). This is an important issue due to the harmful effects on bone metabolism, cancer of the endometrium, thromboembolic diseases, and cognitive disorders caused by chronic oestrogen deprivation in brain tissue (Fisher et al., 1998; Cuzick et al., 2002; Cuzick, 2005). Moreover, various side effects similar to those seen in menopause may occur, such as hot flashes, weight gain, vaginal dryness, nausea, and elevated liver transaminases (Cuzick, 2005; Carpenter and Miller 2005).
Although it is well known that genetic variants in CYP2D6 could explain between 10–20% of these cases and its use as a biomarker for TAM is recommended by the FDA, there is still no agreement on its clinical utility for predicting results in BC, the studies reveal contradictory results and the role of other pharmacokinetic and pharmacodynamic proteins have been poorly studied (Brauch and Schwab, 2014; Brewer et al., 2014; Province et al., 2014; Binkhorst et al., 2015; Goetz et al., 2018). Therefore, the present research focused on studying the association between 7 genetic variants in several genes encoding proteins involved in the pharmacokinetics and pharmacodynamics of TAM and response (relapse) and ADRs in TAM-treated hormone-dependent BC women and contributing to the field, generating predictive models that approximate the response of patients according to their genetic-metabolic characteristics. CYP2D6∗3 and ∗17 alleles were previously found to be absent in different sub-groups of the Chilean population (Roco et al., 2012; Varela et al., 2015) therefore we exclude them from this study and other potentially relevant variants could not be studied, unfortunately.
In this study, we recorded ADRs after 6 months of TAM treatment. This is because several other factors influence TAM ADRs, especially before TAM active metabolites reach steady state, which is after 2 months of treatment (Buckley and Goa, 1989; Lien et al., 1992; Kisanga et al., 2004). Therefore, we considered that after 6 months of treatment, the ADRs might be associated mainly with TAM. Similarly, after 2 years of treatment, relapses can be properly evaluated (Cheng et al., 2012).
The study patients have a 49.3% Ma-c, which agrees with that described by Valenzuela et al. in 1987, in relation to admixture and socioeconomic stratum, because main of the enrolled population of patients (94%) belong to a low socioeconomic stratum with higher Mapuche ancestry; this was an expected result according to the place of recruitment of the patients.
The study of the association between genetic variants and clinical parameters showed that the presence of UGT2B15∗1/∗2 protects against relapse (OR = 0.09; p = 0.008), CYP3A5∗3/∗3 protects against presenting endometrial hyperplasia (OR = 0.07, p = 0.01), and SULT1A1∗1/∗2 protects against vaginal bleeding (OR = 0.09; p = 0.03) and that ESRA 364E/364E is a risk factor for vaginal bleeding (OR = 5.68; p = 0.03). In this respect, it seems to be relevant to perform more studies regarding the ESRA V364E variant, a very poorly studied variant that can be determined in patients with ER + BC to prevent vaginal bleeding caused by TAM therapy, which is also a suspicious sign of endometrial cancer.
In our study, some analyses showed associations without statistical significance, which may be due to the low number of occurrences of the event in each subgroup, lack of clinical data in any of the contrasting groups, or because the variant explains a smaller part of the response, which could be clarified by increasing the sample size in future studies. Thus, in small groups, Fisher’s exact test was used to obtain proper results.
The biostatistical tools used in this field have established that non-significant trends, but with p < 0.2 and physiological significance, can be included in predictive models and in multi-variable logistic regression analysis to help predict the response (relapses) and adverse reactions to treatment (Núñez et al., 2011). Therefore, preliminary predictive models to predict the response (efficacy and safety) in TAM-treated BC patients were generated. The efficacy result, analyzed as relapse, led to a model that explains 49.5% of the possibility of presenting a relapse (p = 0.01) by including SULT1A1∗2 and UGT2B15∗2 and a series of relevant non-genomic variables. These results were expected since both the SULT1A1 and UGT2B15 enzymes have specificity for 4-hydroxyTAM, even though the genetic variant UGT2B15∗2 presents a nucleotide change located in a substrate binding site that is associated with a decrease in catalytic activity and alterations in the Km and Vmax values (Lévesque et al., 1997). Similarly, SULT1A1∗2 has a lower catalytic activity than the wild-type allele (Arslan et al., 2011). Thus, these variants would produce a decrease in the elimination of 4-hydroxyTAM.
Moreover, our results are similar to those published by Nowell and his colleagues (Nowell et al., 2002 and Nowell et al., 2005), who conducted a retrospective study in 162 TAM-treated BC patients from Arkansas who received TAM to determine whether genetic variability (CYP2D6∗4, SULT1A1∗2, and UGT2B15∗2) in the TAM metabolic pathway was associated with recurrence. They found that high activity genotypes UGT2B15 ∗1/∗1 and combined genotypes UGT2B15∗1/∗1 and SULT1A1∗2/∗2 lead to an increased risk of disease recurrence and that SULT1A1∗2/∗2 causes an approximately 3 times higher risk of death than SULT1A1∗1/∗1 and ∗1/∗2. Conversely, Wegman et al. (2005) found that the genetic variants CYP2D6∗4/∗4, and SULT1A1∗1/∗1 cause an increase in the disease-free relapse time in patients with premenopausal and postmenopausal BC treated with 40 mg TAM using proportional hazards Cox regression adjusted for age, tumor size, and lymph nodes. However, Ahern et al. (2011) and Lash et al. (2011) in homologous models found that the genetic variants CYP2D6∗4, UGT1A8∗3, UGT2B7∗2, and UGT2B15∗2 did not modify the relapse rate of the disease in patients with ER+/TAM + BC. On the other hand, Sensorn et al., 2013 investigated genetic variant CYP3A4∗1B (−392A>G, actually CYP3A4∗1.001 G/A)/∗18 (878T>C), CYP3A5∗3 (6986 G>A), ABCB1 3435 C>T, ABCC2∗1C (−24C > T) and ABCC2 68231 A>G to evaluate the risk of recurrence within 3 years among thirty Thai women after receiving tamoxifen adjuvant therapy but they could not find CYP3A4∗1B/∗18, and did not observe a statistical association between CYP3A5∗3, ABCC2∗1C, and ABCC2 68231 A>G with clinical outcome.
Additionally, preliminary predictive models were generated to predict ADRs. The endometrial hyperplasia model explained 41.6% of the possibility of presenting the event (p = 0.002). The model includes the CYP2D6∗4, CYP3A4∗1B (CYP3A4∗1.001), CYP3A5∗3, SULT1A1∗2, UGT2B15∗2 and ESRA V364E genotypes plus several relevant non-genomic factors. Likewise, a preliminary predictive model was generated that explains 33.7% of the possibility of presenting vaginal bleeding (p = 0.014) by including the CYP2D6∗4, CYP3A4∗1B (CYP3A4∗1.001), SULT1A1∗2, and ESRA V364E genotypes and relevant non-genomic factors. Currently, there is no evidence regarding the association of genetic variants with TAM-induced vaginal bleeding or endometrial thickening. Chu et al. (2007) found that women with BC who carried CYP3A4∗1B (CYP3A4∗1.001), had a 3-fold increased risk of developing endometrial cancer after treatment with TAM. The phlebitis model explained 66% of the possibility of presenting this ADR (p = 0.036), a significant finding that has not been previously reported. In the same way, preliminary predictive models were generated to predict the appearance of headache (42.2%; p = 0.028), nausea (36.3%; p = 0.013), and hot flashes (12.1%; p = 0.023) in patients. Our results differ from those reported by Goetz et al. (2005), where CYP2D6∗4/∗4 tends to contribute to a higher risk of disease relapse and a lower incidence of hot flashes in a similar study. No other studies were found in the literature describing associations with headache and nausea. Finally, predictive models to predict the appearance of cramps (24.6%; p = 0.053), bone pain (14.1%; p = 0.017) and urticaria (34%; p = 0.039) were obtained. We were not able to find any similar study in the scientific literature.
Altogether, our results show that some specific genetic variants in genes encoding proteins involved in the TAM pharmacokinetic and/or pharmacodynamic can influence the efficacy of TAM therapy in BC.
Since there are still discrepancies between the findings, which may be due to the genetic differences present in the subjects analyzed, it is crucial to continue with similar and larger studies that confirm the association of CYP2D6∗4, CYP2C9∗2 and ∗3, CYP3A4∗1B (CYP3A4∗1.001), CYP3A5∗3, SULT1A1∗2, UGT2B15∗2 and ESR1 V364E and response to treatment with TAM. Differences in allele frequencies of variants influencing drug’s pharmacokinetics across populations are considered to result in the interethnic differences in pharmacokinetics observed. For example, the lower mean body weights of populations of East Asian in relation to European populations are well established to contribute to differences in the clearance or volume of distribution of some drugs (Lu et al., 2018). On the other hand, in relation to that scarce information is available in Latin American populations, but allele frequencies heterogeneously differ from Asian, Caucasian and African populations (Mendoza et al., 2001; López et al., 2005; Roco et al., 2012; Friedrich et al., 2014), giving rise an idea of the phenotypic response and patient outcome.
Despite our analysis, the present study has some limitations. Even though we had a good sample size for combinatorial analyses, a small number of patients examined in each sub-group comparison could mask potential associations, especially for low-frequency polymorphisms, particularly in the multivariate analyses. The small number of variants analyzed could restrict the statistical power and predictive potential of models, considering that some other candidate genes/polymorphisms were not evaluated in this study (e.g. CYP2D6 ∗2, ∗5, ∗6, ∗10 and ∗41, CYP2C9∗2, CYP2C19 ∗2 and ∗17, UGT1A10 139Lys, UGT1A8 173Gly/277Cys and ABCC2), which could still be relevant. In some cases, missing clinical data was relevant, giving rise to a possible differential misclassification bias affecting the estimated associations between potentially relevant combinations of risk factors and adverse reactions. Other environmental factors, such as contaminants and stress, or social determinants were not explored in this study, which could potentially affect the associations.
However, our findings and preliminary predictive models proposed can be used to prevent, in a good percentage of patients, the negative effects of TAM therapy. Additionally, they could help in the decision-making of a change in therapy to aromatase inhibitors (AIs) (Bernhard et al., 2015).
Certainly, for routine clinical application, the models must be validated in terms of their specificity, sensitivity, and predictive value, so it is necessary to carry out a second stage of research that involves treatment intervention in the patient according to their genotypic profile and relevant non-genomic variables derived from this study. It should be considered that the predictive models generated in this study after being validated could be used in the Chilean population. They could be used in another population with the same percentage of Caucasian aboriginal mix.
5 Conclusion
In our study, UGT2B15*2 A/C is likely to have a clinically significant impact on recurrence risk, CYP3A5∗3 A/G on endometrial hyperplasia, and SULT1A1∗2 G/A, and ESRA V364E on vaginal bleeding in Chilean patients with breast cancer who receive tamoxifen adjuvant therapy. The obtained predictive models showed that the response to TAM treatment in BC patients is partly associated with some of the genetic variants studied. This, after a clinical validation protocol, might help to predict a percentage of disease relapse and negative effects of TAM therapy, improving the efficacy and safety of TAM pharmacotherapeutic treatment.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation, to any qualified researcher.
Ethics Statement
The studies involving human participants were reviewed and approved by the Ethics Committee of the University of Chile (July 24, 2013), the Northern Metropolitan Health Service, National Cancer Institute (May 8, 2013), and Research Ethics Committee of the North Metropolitan Health Service (June 3, 2012. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
CM: The conception of the research, enrolment of patients, experimental analyses, analysis of data and, writing the manuscript. MG: Experimental analyses and analyses of data. RT: Enrolment of patients. KT: Enrolment of patients. DC: Analysis of data. KL: Analysis of data. MR: Analysis of data and writing the manuscript. NV: Analysis of data and writing the manuscript. LQ: The conception of the research, analysis of data and, writing the manuscript.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
The authors wish to thank the breast cancer patients from the “Instituto Nacional del Cáncer” for their altruistic collaboration in pursuit of the common welfare. We also thank the Latin American Society of Pharmacogenomics and Personalized Medicine (SOLFAGEM) by sponsoring this article.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2021.661443/full#supplementary-material
References
Acuña, M., Llop, E., and Rothhammer, F. (2000). Composición genética de la población chilena: las comunidades urales de los valles de Elqui, Limarí y Choapa. Rev. Méd. Chile 128 (6), 593–600. doi:10.4067/S0034-98872000000600004
Ahern, T. P., Christensen, M., Cronin-Fenton, D. P., Lunetta, K. L., Søiland, H., Gjerde, J., et al. (2011). Functional Polymorphisms in UDP-Glucuronosyl Transferases and Recurrence in Tamoxifen-Treated Breast Cancer Survivors. Cancer Epidemiol. Biomarkers Prev. 20 (9), 1937–1943. doi:10.1158/1055-9965.EPI-11-0419
American Cancer Society (2017). “What Is Breast Cancer?” American Cancer Society. Cancer Facts and Figures. Atlanta, GA: American Cancer Society. Available at: http://www.cancer.org/cancer/breast-cancer/about/what-is-breast-cancer.html (Accessed June 20, 2021).
Arslan, S., Silig, Y., and Pinarbasi, H. (2011). Sulfotransferase 1A1 Arg(213)His Polymorphism and Prostate Cancer Risk. Exp. Ther. Med. 2 (6), 1159–1162. doi:10.3892/etm.2011.334
Bernhard, J., Luo, W., Ribi, K., Colleoni, M., Burstein, H. J., Tondini, C., et al. (2015). Patient-reported Outcomes with Adjuvant Exemestane versus Tamoxifen in Premenopausal Women with Early Breast Cancer Undergoing Ovarian Suppression (TEXT and SOFT): a Combined Analysis of Two Phase 3 Randomised Trials. Lancet Oncol. 16, 848–858. doi:10.1016/S1470-2045(15)00049-2
Bijl, M. J., van Schaik, R. H., Lammers, L. A., Hofman, A., Vulto, A. G., van Gelder, T., et al. (2009). The CYP2D6*4 Polymorphism Affects Breast Cancer Survival in Tamoxifen Users. Breast Cancer Res. Treat. 118 (1), 125–130. doi:10.1007/s10549-008-0272-2
Binkhorst, L., Mathijssen, R. H., Jager, A., and van Gelder, T. (2015). Individualization of Tamoxifen Therapy: Much More Than Just CYP2D6 Genotyping. Cancer Treat. Rev. 41 (3), 289–299. doi:10.1016/j.ctrv.2015.01.002
Borges, S., Desta, Z., Li, L., Skaar, T. C., Ward, B. A., Nguyen, A., et al. (2006). Quantitative Effect of CYP2D6 Genotype and Inhibitors on Tamoxifen Metabolism: Implication for Optimization of Breast Cancer Treatment. Clin. Pharmacol. Ther. 80 (1), 61–74. doi:10.1016/j.clpt.2006.03.013
Brauch, H., Mürdter, T. E., Eichelbaum, M., and Schwab, M. (2009). Pharmacogenomics of Tamoxifen Therapy. Clin. Chem. 55 (10), 1770–1782. doi:10.1373/clinchem.2008.121756
Brauch, H., and Schwab, M. (2014). Prediction of Tamoxifen Outcome by Genetic Variation of CYP2D6 in Post-Menopausal Women with Early Breast Cancer. Br. J. Clin. Pharmacol. 77 (4), 695–703. doi:10.1111/bcp.12229
Brewer, N. T., Defrank, J. T., Chiu, W. K., Ibrahim, J. G., Walko, C. M., Rubin, P., et al. (2014). Patients' Understanding of How Genotype Variation Affects Benefits of Tamoxifen Therapy for Breast Cancer. Public Health Genomics 17 (1), 43–47. doi:10.1159/000356565
Buckley, M. M., and Goa, K. L. (1989). Tamoxifen. A Reappraisal of its Pharmacodynamic and Pharmacokinetic Properties, and Therapeutic Use. Drugs 37 (4), 451–490. doi:10.2165/00003495-198937040-00004
Carpenter, R., and Miller, W. R. (2005). Role of Aromatase Inhibitors in Breast Cancer. Br. J. Cancer 93 (Suppl. 1), S1–S5. doi:10.1038/sj.bjc.6602688
Cavalli, S. A., Hirata, M. H., and Hirata, R. D. (2001). Detection of MboII Polymorphism at the 5' Promoter Region of CYP3A4. Clin. Chem. 47 (2), 348–351. doi:10.1093/clinchem/47.2.348
Cheng, L., Swartz, M. D., Zhao, H., Kapadia, A. S., Lai, D., Rowan, P. J., et al. (2012). Hazard of Recurrence Among Women After Primary Breast cancer Treatment--a 10-Year Follow-Up Using Data from SEER-Medicare. Cancer Epidemiol Biomarkers Prev. 21 (5), 800–809. doi:10.1158/1055-9965.EPI-11-1089
Chu, W., Fyles, A., Sellers, E. M., McCready, D. R., Murphy, J., Pal, T., et al. (2007). Association between CYP3A4 Genotype and Risk of Endometrial Cancer Following Tamoxifen Use. Carcinogenesis 28 (10), 2139–2142. doi:10.1093/carcin/bgm087
Clarke, R., Liu, M. C., Bouker, K. B., Gu, Z., Lee, R. Y., Zhu, Y., et al. (2003). Antiestrogen Resistance in Breast Cancer and the Role of Estrogen Receptor Signaling. Oncogene 22 (47), 7316–7339. doi:10.1038/sj.onc.1206937
Cuzick, J. (2005). Aromatase Inhibitors for Breast Cancer Prevention. J. Clin. Oncol. 23 (8), 1636–1643. doi:10.1200/JCO.2005.11.027
Cuzick, J., ForbesEdwards, J. R., Edwards, R., Baum, M., Cawthorn, S., Coates, A., et al. (2002). First Results from the International Breast Cancer Intervention Study (IBIS-I): A Randomised Prevention Trial. Lancet 360 (9336), 817–824. doi:10.1016/s0140-6736(02)09962-2
Davies, C., Davies, C., Godwin, J., Gray, R., Clarke, M., Cutter, D., et al. (2011). Relevance of Breast Cancer Hormone Receptors and Other Factors to the Efficacy of Adjuvant Tamoxifen: Patient-Level Meta-Analysis of Randomised Trials. Lancet 378 (9793), 771–784. doi:10.1016/S0140-6736(11)60993-8
Fernández-Santander, A., Gómez, A. T., Bernal, A. T., Gómez-Gallego, F., Santiago Dorrego, C., and Moya, F. B. (2007). Aspectos Farmacogenéticos En El Tratamiento Del Cáncer de Mama Con Tamoxifeno. Biological Sciences. Madrid, Spain: ADEMAS Com, FundaciÓn Tejerina. Available at: http://books.google.com/books?hl=en& lr=& id=GanArxPeUeAC& oi=fnd& pg=PA5& dq=Aspectos+farmacogenéticos+en+el+tratamiento+del+cáncer+de+mama+con+Tamoxifeno& ots=P2jmC77GWg& sig=GWxp_taJbt5AJiCdE3MHLyfUZUg.
Fisher, B., Costantino, J. P., Wickerham, D. L., Redmond, C. K., Kavanah, M., Cronin, W. M., et al. (1998). Tamoxifen for Prevention of Breast Cancer: Report of the National Surgical Adjuvant Breast and Bowel Project P-1 Study. J. Natl. Cancer Inst. 90 (18), 1371–1388. doi:10.1093/jnci/90.18.1371
Flicek, P., Amode, M. R., Barrell, D., Beal, K., Brent, S., Chen, Y., et al. (2011). Ensembl 2011. Nucleic Acids Res. 39 (Database issue), D800–D806. doi:10.1093/nar/gkq1064
Friedrich, D. C., Genro, J. P., Sortica, V. A., Suarez-Kurtz, G., de Moraes, M. E., Pena, S. D., et al. (2014). Distribution of CYP2D6 Alleles and Phenotypes in the Brazilian Population. PLoS One 9 (10), e110691. doi:10.1371/journal.pone.0110691
Fuentes, N., and Silveyra, P. (2019). Estrogen Receptor Signaling Mechanisms. Adv. Protein Chem. Struct. Biol. 116, 135–170. doi:10.1016/bs.apcsb.2019.01.001
Gjerde, J., Hauglid, M., Breilid, H., Lundgren, S., Varhaug, J. E., Kisanga, E. R., et al. (2007). Effects of CYP2D6 and SULT1A1 Genotypes Including SULT1A1 Gene Copy Number on Tamoxifen Metabolism. Ann. Oncol. 19 (1), 56–61. doi:10.1093/annonc/mdm434
Goetz, M. P., Rae, J. M., Suman, V. J., Safgren, S. L., Ames, M. M., Visscher, D. W., et al. (2005). Pharmacogenetics of Tamoxifen Biotransformation Is Associated with Clinical Outcomes of Efficacy and Hot Flashes. J. Clin. Oncol. 23 (36), 9312–9318. doi:10.1200/JCO.2005.03.3266
Goetz, M. P., Sangkuhl, K., Guchelaar, H. J., Schwab, M., Province, M., Whirl-Carrillo, M., et al. (2018). Clinical Pharmacogenetics Implementation Consortium (CPIC) Guideline for CYP2D6 and Tamoxifen Therapy. Clin. Pharmacol. Ther. 103 (5), 770–777. doi:10.1002/cpt.1007
Hajdinjak, T., and Zagradišnik, B. (2004). Prostate Cancer and Polymorphism D85Y in Gene for Dihydrotestosterone Degrading Enzyme UGT2B15: Frequency of DD Homozygotes Increases with Gleason Score. Prostate 59 (4), 436–439. doi:10.1002/pros.20024
Hertz, D. L., McLeod, H. L., and Hoskins, J. M. (2009). Pharmacogenetics of Breast Cancer Therapies. Breast 18 (Suppl. 3), S59–S63. doi:10.1016/S0960-9776(09)70275-9
Herynk, M. H., and Fuqua, S. A. (2004). Estrogen Receptor Mutations in Human Disease. Endocr. Rev. 25 (6), 869–898. doi:10.1210/er.2003-0010
Ince, B. A., Schodin, D. J., Shapiro, D. J., and Katzenellenbogen, B. S. (1995). Repression of Endogenous Estrogen Receptor Activity in MCF-7 Human Breast Cancer Cells by Dominant Negative Estrogen Receptors. Endocrinology 136 (8), 3194–3199. doi:10.1210/endo.136.8.7628351
Jin, Y., Desta, Z., Stearns, V., Ward, B., Ho, H., Lee, K. H., et al. (2005). CYP2D6 Genotype, Antidepressant Use, and Tamoxifen Metabolism during Adjuvant Breast Cancer Treatment. J. Natl. Cancer Inst. 97 (January), 30–39. doi:10.1093/jnci/dji005
Johnson, M. D., Zuo, H., Lee, K. H., Trebley, J. P., Rae, J. M., Weatherman, R. V., et al. (2004). Pharmacological Characterization of 4-Hydroxy-N-Desmethyl Tamoxifen, a Novel Active Metabolite of Tamoxifen. Breast Cancer Res. Treat. 85 (2), 151–159. doi:10.1023/B:BREA.0000025406.31193.e8
Jordan, V. C., Collins, M. M., Rowsby, L., and Prestwich, G. (1977). A Monohydroxylated Metabolite of Tamoxifen with Potent Antioestrogenic Activity. J. Endocrinol. 75 (2), 305–316. doi:10.1677/joe.0.0750305
Kagaya, H., Inoue, K., Miura, M., Satoh, S., Saito, M., Tada, H., et al. (2007). Influence of UGT1A8 and UGT2B7 Genetic Polymorphisms on Mycophenolic Acid Pharmacokinetics in Japanese Renal Transplant Recipients. Eur. J. Clin. Pharmacol. 63 (3), 279–288. doi:10.1007/s00228-006-0248-2
Kisanga, E. R., Gjerde, J., Guerrieri-Gonzaga, A., Pigatto, F., Pesci-Feltri, A., Robertson, C., et al. (2004). Tamoxifen and Metabolite Concentrations in Serum and Breast Cancer Tissue during Three Dose Regimens in a Randomized Preoperative Trial. Clin. Cancer Res. 10 (7), 2336–2343. doi:10.1158/1078-0432.ccr-03-0538
Kiyotani, K., Mushiroda, T., Imamura, C. K., Hosono, N., Tsunoda, T., Kubo, M., et al. (2010). Significant Effect of Polymorphisms in CYP2D6 and ABCC2 on Clinical Outcomes of Adjuvant Tamoxifen Therapy for Breast Cancer Patients. J. Clin. Oncol. 28 (8), 1287–1293. doi:10.1200/JCO.2009.25.7246
Lash, T. L., Cronin-Fenton, D., Ahern, T. P., Rosenberg, C. L., Lunetta, K. L., Silliman, R. A., et al. (2011). CYP2D6 Inhibition and Breast Cancer Recurrence in a Population-Based Study in Denmark. J. Natl. Cancer Inst. 103 (6), 489–500. doi:10.1093/jnci/djr010
Lee, S. J., Bell, D. A., Coulter, S. J., Ghanayem, B., and Goldstein, J. A. (2005). Recombinant CYP3A4*17 Is Defective in Metabolizing the Hypertensive Drug Nifedipine, and the CYP3A4*17 Allele May Occur on the Same Chromosome as CYP3A5*3, Representing a New Putative Defective CYP3A Haplotype. J. Pharmacol. Exp. Ther. 313 (1), 302–309. doi:10.1124/jpet.104.078758
Lévesque, E., Beaulieu, M., Green, M. D., Tephly, T. R., Bélanger, A., and Hum, D. W. (1997). Isolation and Characterization of UGT2B15(Y85): A UDP-Glucuronosyltransferase Encoded by a Polymorphic Gene. Pharmacogenetics 7 (4), 317–325. doi:10.1097/00008571-199708000-00007
Lewis, J. S., and Jordan, V. C. (2005). Selective Estrogen Receptor Modulators (SERMs): Mechanisms of Anticarcinogenesis and Drug Resistance. Mutat. Res. 591 (1–2), 247–263. doi:10.1016/j.mrfmmm.2005.02.028
Lien, E. A., Johannessen, D. C., Aakvaag, A., and Lønning, P. E. (1992). Influence of Tamoxifen, Aminoglutethimide and Goserelin on Human Plasma IGF-I Levels in Breast Cancer Patients. J. Steroid Biochem. Mol. Biol. 41 (3–8), 541–543. doi:10.1016/0960-0760(92)90380-2
Lim, J. S., Chen, X. A., Singh, O., Yap, Y. S., Ng, R. C., Wong, N. S., et al. (2011). Impact of CYP2D6, CYP3A5, CYP2C9 and CYP2C19 Polymorphisms on Tamoxifen Pharmacokinetics in Asian Breast Cancer Patients. Br. J. Clin. Pharmacol. 71 (5), 737–750. doi:10.1111/j.1365-2125.2011.03905.x
López, M., Guerrero, J., Jung-Cook, H., and Alonso, M. E. (2005). CYP2D6 Genotype and Phenotype Determination in a Mexican Mestizo Population. Eur. J. Clin. Pharmacol. 61 (10), 749–754. doi:10.1007/s00228-005-0038-2
Lu, S., Nand, R. A., Yang, J. S., Chen, G., and Gross, A. S. (2018). Pharmacokinetics of CYP2C9, CYP2C19, and CYP2D6 Substrates in Healthy Chinese and European Subjects. Eur. J. Clin. Pharmacol. 74, 285–296. doi:10.1007/s00228-017-2375-3
Ma, Q., and Lu, A. Y. (2011). Pharmacogenetics, Pharmacogenomics, and Individualized Medicine. Pharmacol. Rev. 63 (2), 437–459. doi:10.1124/pr.110.003533
Madlensky, L., Natarajan, L. S. T., Tchu, S., Pu, M. S. F., Mortimer, J., Flatt, S. W., et al. (2011). Tamoxifen Metabolite Concentrations, CYP2D6 Genotype, and Breast Cancer Outcomes. Clin. Pharmacol. Ther. 89 (5), 718–725. doi:10.1038/clpt.2011.32
Maughan, K. L., Lutterbie, M. A., and Ham, P. S. (2010). Treatment of Breast Cancer. Am. Fam. Physician 81 (June), 1339–1346.
McDonnell, D. P., and Wardell, S. E. (2010). The Molecular Mechanisms Underlying the Pharmacological Actions of ER Modulators: Implications for New Drug Discovery in Breast Cancer. Curr. Opin. Pharmacol. 10 (6), 620–628. doi:10.1016/j.coph.2010.09.007
McInerney, E. M., Ince, B. A., Shapiro, D. J., and Katzenellenbogen, B. S. (1996). A Transcriptionally Active Estrogen Receptor Mutant Is a Novel Type of Dominant Negative Inhibitor of Estrogen Action. Mol. Endocrinol. 10 (12), 1519–1526. doi:10.1210/mend.10.12.8961262
Mendoza, R., Wan, Y. J., Poland, R. E., Smith, M., Zheng, Y., Berman, N., et al. (2001). CYP2D6 Polymorphism in a Mexican American Population. Clin. Pharmacol. Ther. 70 (6), 552–560. doi:10.1067/mcp.2001.120675
MINSAL (2005). Guía Clínica Cáncer de Mama En Personas de 15 Años y Más. 1st Ed. Santiago: Ministerio de Salud de Chile.
Mürdter, T. E., Schroth, W., Bacchus-Gerybadze, L., Winter, S., Heinkele, G., Simon, W., et al. (2011). Activity Levels of Tamoxifen Metabolites at the Estrogen Receptor and the Impact of Genetic Polymorphisms of Phase I and II Enzymes on Their Concentration Levels in Plasma. Clin. Pharmacol. Ther. 89 (5), 708–717. doi:10.1038/clpt.2011.27
Namipashaki, A., Razaghi-Moghadam, Z., and Ansari-Pour, N. (2015). The Essentiality of Reporting Hardy-Weinberg Equilibrium Calculations in Population-Based Genetic Association Studies. Cel J. 17, 187. doi:10.22074/cellj.2016.3711
NCI (2012). Terapia Hormonal Para El Cáncer de Seno. Available at: http://www.cancer.gov/espanol/recursos/hojas-informativas/tratamiento/terapia-hormonal-seno (Accessed June 20, 2021).
Nishiyama, T., Ogura, K., Nakano, H., Ohnuma, T., Kaku, T., Hiratsuka, A., et al. (2002). Reverse Geometrical Selectivity in Glucuronidation and Sulfation of Cis- and Trans-4-hydroxytamoxifens by Human Liver UDP-Glucuronosyltransferases and Sulfotransferases. Biochem. Pharmacol. 63 (10), 1817–1830. doi:10.1016/s0006-2952(02)00994-2
Nowell, S., Sweeney, C., Winters, M., Stone, A., Lang, N. P., Hutchins, L. F., et al. (2002). Association between Sulfotransferase 1A1 Genotype and Survival of Breast Cancer Patients Receiving Tamoxifen Therapy. J. Natl. Cancer Inst. 94 (21), 1635–1640. doi:10.1093/jnci/94.21.1635
Nowell, S. A., Ahn, J., Rae, J. M., Scheys, J. O., Trovato, A., Sweeney, C., et al. (2005). Association of Genetic Variation in Tamoxifen-Metabolizing Enzymes with Overall Survival and Recurrence of Disease in Breast Cancer Patients. Breast Cancer Res. Treat. 91 (3), 249–258. doi:10.1007/s10549-004-7751-x
Núñez, E., Steyerberg, E. W., and Núñez, J. (2011). Regression Modeling Strategies. Rev. Esp Cardiol. 64 (6), 501–507. doi:10.1016/j.recesp.2011.01.019
Province, M. A., Goetz, M. P., Brauch, H., Flockhart, D. A., Hebert, J. M., Whaley, R., et al. (2014). CYP2D6 Genotype and Adjuvant Tamoxifen: Meta-Analysis of Heterogeneous Study Populations. Clin. Pharmacol. Ther. 95 (2), 216–227. doi:10.1038/clpt.2013.186
Rae, J. M., Drury, S., Hayes, D. F., Stearns, V., Thibert, J. N., Haynes, B. P., et al. (2012). CYP2D6 and UGT2B7 Genotype and Risk of Recurrence in Tamoxifen-Treated Breast Cancer Patients. J. Natl. Cancer Inst. 104 (6), 452–460. doi:10.1093/jnci/djs126
Ramón y Cajal, T., Altés, A., Paré, L., del Rio, E., Alonso, C., Barnadas, A., et al. (2010). Impact of CYP2D6 Polymorphisms in Tamoxifen Adjuvant Breast Cancer Treatment. Breast Cancer Res. Treat. 119 (1), 33–38. doi:10.1007/s10549-009-0328-y
Kumar, V., Abbas, A. K., and Aster, J. C. (2010). Robbins y Cotran. Patología estructural y funcional. 9th Ed. Barcelona: Elsevier España SA.
Roco, A., Quiñones, L., Agúndez, J. A., García-Martín, E., Squicciarini, V., Miranda, C., et al. (2012). Frequencies of 23 Functionally Significant Variant Alleles Related with Metabolism of Antineoplastic Drugs in the Chilean Population: Comparison with Caucasian and Asian Populations. Front. Genet. 3, 229. doi:10.3389/fgene.2012.00229
Schroth, W., Antoniadou, L., Fritz, P., Schwab, M., Muerdter, T., Zanger, U. M., et al. (2007). Breast Cancer Treatment Outcome with Adjuvant Tamoxifen Relative to Patient CYP2D6 and CYP2C19 Genotypes. J. Clin. Oncol. 25 (33), 5187–5193. doi:10.1200/JCO.2007.12.2705
Schur, B. C., Bjerke, J., Nuwayhid, N., and Wong, S. H. (2001). Genotyping of Cytochrome P450 2D6*3 and *4 Mutations Using Conventional PCR. Clin. Chim. Acta 308 (1–2), 25–31. doi:10.1016/s0009-8981(01)00422-3
Sensorn, I., Sirachainan, E., Chamnanphon, M., Pasomsub, E., Trachu, N., Supavilai, P., et al. (2013). Association of CYP3A4/5, ABCB1 and ABCC2 Polymorphisms and Clinical Outcomes of Thai Breast Cancer Patients Treated with Tamoxifen. Pharmgenomics Pers Med. 6 (1), 93–98. doi:10.2147/PGPM.S44006
Terán, J., and Teppa, A. (2005). Moduladores Selectivos de Los Receptores Estrogénicos (SERMs): Bioquímica, Farmacología y Aplicación En Guinecología. Ginecol. Obstet. Mex 73, 424–435.
Tucker, A. N., Tkaczuk, K. A., Lewis, L. M., Tomic, D., Lim, C. K., and Flaws, J. A. (2005). Polymorphisms in Cytochrome P4503A5 (CYP3A5) May Be Associated with Race and Tumor Characteristics, but Not Metabolism and Side Effects of Tamoxifen in Breast Cancer Patients. Cancer Lett. 217 (1), 61–72. doi:10.1016/j.canlet.2004.08.027
Untch, M., and Thomssen, C. (2010). Clinical Practice Decisions in Endocrine Therapy. Cancer Invest. 28 (January), 4–13. doi:10.3109/07357907.2010.501637
Valenzuela, C., Acuña, M., and Harb, Z. (1987). Gradiente Sociogenetico En La Poplación Chilena. Rev. Med. Chil 115, 295–299.
Varela, N., Quiñones, L. A., Stojanova, J., Garay, J., Cáceres, D., Cespedes, S., et al. (2015). Characterization of the CYP2D6 Drug Metabolizing Phenotypes of the Chilean Mestizo Population through Polymorphism Analyses. Pharmacol. Res. 101, 124–129. doi:10.1016/j.phrs.2015.07.020
Wakeling, A. E., and Slater, S. R. (1980). Estrogen-Receptor Binding and Biologic Activity of Tamoxifen and its Metabolites. Cancer Treat. Rep. 64 (6–7), 741–744.
Wegman, P., Vainikka, L., Stål, O., Nordenskjöld, B., Skoog, L., Rutqvist, L. E., et al. (2005). Genotype of Metabolic Enzymes and the Benefit of Tamoxifen in Postmenopausal Breast Cancer Patients. Breast Cancer Res. 7 (3), R284–R290. doi:10.1186/bcr993
Whirl-Carrillo, M., McDonagh, E. M., Hebert, J. M., Gong, L., Sangkuhl, K., Thorn, C. F., et al. (2012). Pharmacogenomics Knowledge for Personalized Medicine. Clin. Pharmacol. Ther. 92 (4), 414–417. doi:10.1038/clpt.2012.96
WMA (2013). World Medical Association Declaration of Helsinki: Ethical Principles for Medical Research Involving Human Subjects. JAMA 310 (20), 2191–2194. doi:10.1001/jama.2013.281053
Wrenn, C. K., and Katzenellenbogen, B. S. (1993). Structure-Function Analysis of the Hormone Binding Domain of the Human Estrogen Receptor by Region-specific Mutagenesis and Phenotypic Screening in Yeast. J. Biol. Chem. 268 (32), 24089–24098. doi:10.1016/s0021-9258(20)80497-9
Keywords: breast cancer, pharmacokinetics, pharmacodynamics, polymorphism, ADR, Relapse
Citation: Miranda C, Galleguillos M, Torres R, Tardón K, Cáceres DD, Lee K, Redal MA, Varela NM and Quiñones LA (2021) Preliminary Pharmacogenomic-Based Predictive Models of Tamoxifen Response in Hormone-dependent Chilean Breast Cancer Patients. Front. Pharmacol. 12:661443. doi: 10.3389/fphar.2021.661443
Received: 30 January 2021; Accepted: 25 October 2021;
Published: 25 November 2021.
Edited by:
Chonlaphat Sukasem, Mahidol University, ThailandReviewed by:
Umamaheswaran Gurusamy, University of California, San Francisco, United StatesChalirmporn Atasilp, Thammasat University, Thailand
Monpat Chamnanphon, Srinakharinwirot University, Thailand
Copyright © 2021 Miranda, Galleguillos, Torres, Tardón, Cáceres, Lee, Redal, Varela and Quiñones. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Luis A. Quiñones, lquinone@uchile.cl; Nelson M. Varela, nvarela@med.uchile.cl