rccqf Revista Colombiana de Ciencias Químico - Farmacéuticas Rev. colomb. cienc. quim. farm. 0034-7418 Departamento de Farmácia, Facultad de Ciencias, Universidade Nacional da Colombia 10.15446/rcciquifa.v47n2.73971 Artículos de Investigación Científica Frequency of 14 genetic variants associated with breast cancer risk and treatment in a Colombian population Frecuencia de 14 variantes genéticas asociadas con el riesgo y el tratamiento del cáncer de mama en una población colombiana Ariza Márquez Yeimy Viviana 1 * Aristizábal Gutiérrez Fabio Ancízar 2 Gómez-Gutiérrez Alberto 3 Gómez Piedad Elena 4 Casas Gómez María Consuelo 4 Balcázar Ignacio Briceño 5 Programa de doctorado en Biotecnología, Universidad Nacional de Colombia, sede Bogotá, Colombia, Universidad Nacional de Colombia Programa de doctorado en Biotecnología Universidad Nacional de Colombia Bogotá Colombia Departamento de Farmacia e Instituto de Biotecnología, Universidad Nacional de Colombia, sede Bogotá, Colombia. Universidad Nacional de Colombia artamento de Farmacia e Instituto de Biotecnología Universidad Nacional de Colombia Bogotá Colombia Instituto de Genética Humana, Facultad de Medicina, Pontificia Universidad Javeriana, Bogotá, Colombia. Pontificia Universidad Javeriana Instituto de Genética Humana Facultad de Medicina Pontificia Universidad Javeriana Bogotá Colombia Instituto de Referencia Andino, Bogotá, Colombia. Instituto de Referencia Andino Bogotá Colombia Facultad de Medicina, Universidad de la Sabana, Chía, Colombia. Universidad de la Sabana Facultad de Medicina Universidad de la Sabana Chía Colombia E-mail: yvarizam@unal.edu.co. May-Aug 2018 47 2 277 288 21 04 2017 24 04 2018 This is an open-access article distributed under the terms of the Creative Commons Attribution License SUMMARY <italic>Introduction:</italic>

Genetic variations have been related to risk and treatment efficacy. Many polymorphisms in breast cancer are known to influence susceptibility, breast cancer risk and treatment outcome. Polymorphisms vary among populations; therefore, local studies are necessary.

 <italic>Objective:</italic>

To establish the frequency of polymorphisms associated to breast cancer risk and treatment pharmacogenomics in a group of Colombian individuals.

 <italic>Methods:</italic>

Data from microarray profiles including gene polymorphisms associated with breast cancer treatment were retrospectively collected (Pathway Genomics®). The frequency of marker CYP2D6 rs3892097 and a breast cancer panel (CAS8 rs1045485, CHEK21100delC, ESR1 rs2046210, FGFR2 rs1219648, intergenic_2q35rs13387042, intergenic_8q24 rs13281615, MSRP30 rs10941679, TNRC9 rs3803662, AKAP9 rs6964587, LSP1 rs3817198, MAP3K1rs889312, PALBS1592 delT, ESR1 rs3020314) were studied.

 <italic>Results:</italic>

Microarray data from 68 men and 92 women were analyzed. All polymorphisms were in Hardy-Weinberg equilibrium. Genotypic frequencies of CYP2D6 rs3892097 C/T, CAS8 rs1045485 G/C, and those of genes included in a breast cancer panel (CAS8 rs1045485, CHEK21100delC, FGFR2rs1219648, intergenic_2q35rs13387042, intergenic_8q24 rs13281615, MSRP30 rs10941679, TNRC9 rs3803662, LSP1 rs3817198, MAP3K1rs889312, PALBS1592 del T, ESR1rs3020314) did not significantly differ from previously published data. ESR1 rs2046210, with allele frequencies of C=0.04 and T=0.02, and AKAP9 rs6964587, with a frequency of A=0.005, were determined as rare.

 <italic>Conclusions:</italic>

The population studied was not significantly different in allele distribution from previously reported data at HapMap. Genotypes in Colombian population are similar to other previously studied groups of healthy subjects. Extended use of genotyping pharmacogenetic polymorphisms will prevent toxicity and adverse effects in tamoxifen treatment (for example in CYP2D6 rs3892097). Therefore, therapeutic alternatives should be evaluated based on individual pharmacogenetic studies.

 RESUMEN <italic>Introducción:</italic>

las variaciones genéticas se han relacionado con el riesgo y la eicacia del tratamiento. Es sabido que muchos polimorfismos en cáncer de mama influyen en la susceptibilidad, el riesgo de cáncer y el resultado del tratamiento. Los polimorfismos varían entre las poblaciones, y por tanto, es necesario realizar estudios locales.

 <italic>Objetivo:</italic>

establecer la frecuencia de polimorismos asociados al riesgo de cáncer de mama y la farmacogenómica del tratamiento en un grupo de individuos colombianos.

 <italic>Métodos:</italic>

los datos de los perfiles de microarreglos, incluidos los polimorismos genéticos asociados con el tratamiento del cáncer de mama, se obtuvieron de forma retrospectiva (Pathway Genomics®). Se estudiaron la frecuencia del marcador CYP2D6 rs3892097 y un panel de cáncer de mama (CAS8 rs1045485, CHEK21100delC, ESR1 rs2046210,FGFR2 rs1219648, intergenic_2q35rs13387042, intergenic_8q24 rs13281615, MSRP30 rs10941679, TNRC9 rs3803662, AKAP9 rs6964587, LSP1 rs3817198, MAP3K1rs889312, PALBS1592 delT, ESR1 rs3020314).

 <italic>Resultados:</italic>

se analizaron los datos de microarreglos de 68 hombres y 92 mujeres. Todos los polimorfismos siguieron el equilibrio Hardy-Weinberg. Las frecuencias fenotípicas de CYP2D6 rs3892097 C/T, CAS8 rs1045485 G/C, y aquellas de los genes incluidos en un panel de cáncer de mama (CAS8 rs1045485, CHEK21100delC, FGFR2rs1219648, intergenic_2q35rs13387042, intergenic_8q24 rs13281615, MSRP30 rs10941679, TNRC9 rs3803662, LSP1 rs3817198, MAP3K1rs889312, PALBS1592 del T, ESR1rs3020314) no difirieron significativamente de los datos publicados previamente. ESR1 rs2046210, con frecuencias alélicas de C = 0,04 y T = 0,02, y AKAP9 rs6964587, con una frecuencia de A = 0,005, se determinaron como raras.

 <italic>Key words:</italic> Colombia pharmacogenomics breast Cancer <italic>Palabras clave:</italic> Colombia farmacogenómica cáncer de mama INTRODUCTION

Advances in technology have made multigene testing, or "panel testing", a practical option when looking for genetic variants that may be associated with a risk of breast cancer and treatment outcome. Use of pharmacogenomics will optimize treatment and lower health costs. In this context, it is necessary to study disease associations and pharmacogenetic polymorphisms in populations to optimize disease treatment. In the case of breast cancer, tamoxifen metabolism related polymorphisms are associated to treatment, but their frequency is poorly known in Colombian population.

Using different platforms to process and analyze samples can simultaneously identify genetic variants associated with susceptibility, clinical outcomes and drug response. For example, the spectrum of the frequency of the MTHFR gene polymorphism C677T in a group of healthy Colombian individuals was previously determined to evaluate the genetic background in relation to disease susceptibility and pharmacogenetic applications; the results for the C/C and C/T genotypes in a Colombian population were similar to other previously studied groups of healthy subjects. Subjects from Colombian populations might be at risk of developing diseases associated with MTHFR polymorphisms and might present toxicity and adverse effects if treated with MTX, which suggests the need to evaluate therapeutic alternatives based on individual pharmacogenetic studies [1].

In the case of breast cancer, the response to several hormone therapeutics, such as tamoxifen are mediated by the action of CYP2D6, which is highly polymorphic, and genotypes has been identified allowing classifying patients as ultra-rapid, normal and poor metabolizers.

Genetic variants in certain genes have been associated with early cancer risk [2], risk of developing cancer [3] and susceptibility to cancer [4]. In this way, CYP2D6*4 variant (rs3892097) is associated with response to tamoxifen, since patients carrying a T allele will present slow metabolism of this drug and would therefore need a dose adjustment or an alternative therapy [5]. In the case of CAS8 rs1045485, its C allele is associated with a lower risk of breast cancer [6,7]. On the other hand, ESR1 rs2046210 has been associated with risk of developing breast and endometrial cancer in Chinese women [8-10], and variant C has been associated with ER-negative but not ER-positive breast cancer (P > 0.05). These findings provide further evidence for distinct etiological pathways associated with invasive ER-positive and ER-negative breast cancer in younger women and women of African ancestry [11]. FGFR2rs1219648 showed a similar association with the risk of developing breast cancer in Chinese [12] and Indian populations [13]. Polymorphism rs13387042-2q35(G/A), has been associated with breast cancer risk only among postmenopausal women who never used hormone therapy [14]. Furthermore, the GG genotype of SNP rs13281615 plays a role in breast cancer through PVT1 expression during oncogenesis, and "protective" mutations could occur [15]. MSRP30 rs10941679 was significantly associated with BCIS (Breast cancer in situ) risk [16]. As for TNRC9 rs3803662, allele GG increased the risk of breast cancer in codominant inheritance (OR=2.19, 95% CI:1.19-4.02), and recessive genetic models (OR=2.06, 95% CI: 1.15-3.70) [17]. Moreover, the variant rs6964587 was associated with increased breast edema due to radiotherapy after the following 5 years (Beta, 0.22; 95% confidence interval, 0.09-0.34; P = 7 x 10-4) [18] and increased to risk of breast cancer [19]. LSP1 rs3817198 has been associated with mammographic density in genome-wide studies [20]. Likewise, MAP3K1rs889312, has been associated with breast cancer risk in women of European ancestry [21] and ESR1rs3020314 showed stronger associations in ER-positive [22]. Finally, in AKAP9, the minor T allele of 7q21-rs6964587 was associated with breast cancer risk under a recessive model (OR 1.07, 95% CI 1.00 to 1.13, p = 0.04) [23].

Knowledge of Colombian disease associated polymorphism frequencies will allow to establish National Health policies to improve health and treatment outcome.

 METHODS Population

A retrospective cross-sectional study in 160 healthy Colombian individuals over 18 years was conducted using molecular genetic profiling, age ranges between 18 and 74 years with a mean of 40.07 years. DNA samples were collected from saliva samples in Instituto de Referencia Andino, as well as all corresponding applications from medical orders and/or individuals received between the years 2012 and 2014. All individuals signed an informed consent.

The Healthy Woman DNA Insight™ panel from Pathway Genomics® evaluates panel of breast cancer (CAS8 rs1045485, CHEK21100delC, ESR1 rs2046210, FGFR2rs1219648, intergenic_2q35rs13387042, intergenic_8q24 rs13281615, MSRP30 rs10941679, TNRC9 rs3803662, AKAP9 rs6964587, LSP1 rs3817198, MAP3K1rs889312, PALBS1592 delT, ESR1rs3020314). The pharmacogenetics (CYP2D6*4), in addition to Abacavir® hypersensitivity, aminoglycoside-induced ototoxicity, response to beta blockers, carbamazepine hypersensitivity, clopidogrel metabolism, estrogen supplementation, alpha interferon, metoprolol metabolism, metabolism and hypersensitivity to phenytoin, proton-pump inhibitor, simvastatin induced myopathy, metabolism of voriconazole and warfarin.

DNA microarrays were used to measure expression levels of probes, corresponding to alleles in CAS8 rs1045485, CHEK21100delC, ESR1 rs2046210,FGFR2rs1219648, intergenic_2q35rs13387042, intergenic_8q24 rs13281615, MSRP30 rs10941679, TNRC9 rs3803662,AKAP9 rs6964587, LSP1 rs3817198, MAP3K1rs889312, PALBS1592 delT, ESR1rs3020314 genes. Each microarray position contains 10-12 picomoles of a specific DNA sequence under high stringency. Probe-target hybridization was detected and quantified by fluorophores attached to the probe.

The frequency of the polymorphism substitution in the genes markers was analyzed. Each participant was identified by a code used in a database, where the descriptive analysis for continuous variables such as age was performed with measures of central tendency and Shapiro-Wilk normality test. The description of categorical variables (origin, sex, genotype and allele) was performed by frequency analysis. Given the variable nature of the study association levels, they were evaluated by a Chi square test (X2) with a confidence level of 95%. In addition, Hardy-Weinberg equilibrium was tested. Data base designed in Access includes the following variables that were Collected: Code anonymity of each individual, age, gender, date of birth, maternal ancestry, paternal ancestry, nationality, CYP2D6 rs3892097, panel for breast cancer (CAS8 rs1045485, CHEK21100delC, ESR1 rs2046210, FGFR2rs1219648, intergenic_2q35rs13387042, intergenic_8q24 rs13281615, MSRP30 rs10941679, TNRC9 rs3803662,AKAP9 rs6964587, LSP1 rs3817198, MAP3K1rs889312, PALBS1592 delT, ESR1rs3020314).

 Statistical Methods

Available results were tabulated to perform the respective calculation of frequencies of each polymorphism and Hardy-Weinberg equilibrium was calculated using the online Curt Lab - HW calculator, and then a Genpop software was used in order to define genetic frequencies. Later on, contingency tables were established with a double entry and Chi squared distribution calculated. Results were transferred to the SPSS statistical program for Windows V20. They were not stratified by ethnic group or socioeconomic levels. Available frequencies were collected in HapMap Genome Browser release #28 (Phases 1, 2 and 3 - merged genotypes and frequencies) for the general population and also for the Colombian population for each of the variants tested in this review.

 Ethical considerations

This study was approved by the Ethics Committee of Universidad Nacional de Colombia.

 RESULTS

No deviation from Hardy-Weinberg equilibrium is present; therefore, the population included in this study showed no significant migration, subgroup stratification, inbreeding or selection processes. Additionally, by contrasting data reported in Hap-Map for general population, Colombians of the present study have no significant differences with the sample of Colombian population reported in HapMap.

After evaluatingeach particular polymorphism, CYP2D6*4 rs3892097, CAS8rs1045485, MSRP30 rs10941679, LSP1 rs3817198 and MAP3K1 rs889312, data found in the present study are similar to data report for general and Colombian population previously evaluated in HapMap. Allele Frequency Comparison between Present Study and Hap-Map is shown in Table 1.

Allele frequency comparison between a Colombian population in the present study and HapMap database.

Bar chart for gene frequencies.

Polymorphisms FGFR2 rs1219648, intergenic 2q35-rs13387042, intergenic 8q24-rs13281615, TNRC9 rs3803662 and ESR1 rs3020314, showed higher frequencies than those reported previously for the Colombian population in HapMap, and has greater similarity to those reported for the global population.

In addition, we identified in ESR1 rs2046210 variants C and G have frequencies of 0.04 and 0.02, respectively; similarly, AKAP9 rs6964587 allele A was found in a frequency of0.005. The presence of these variants had not been previously reported in the Colombian population or the HapMap database.

 DISCUSSION

As we move from gene disease discovery into the application of genetic knowledge, allele frequency studies in populations become more important for the design, implementation of studies and for determining treatment according to personal polymorphisms. Studies in Colombian population are scarce, and the present study determined the allele frequencies of 14 polymorphisms, previously associated with breast cancer status. Allele frequencies were compared with dbSNP and other data sources in HapMap.

Allele frequencies for minor alleles relevant for Colombian population are here reported, although the number of individuals within the categories is small. Local differences in allele frequencies should be considered when designing association studies in the future. One weakness of our analysis is the lack of treatment follow-up due to the diversity of places where patients are treated.

The number of useful biomarkers for predicting drug response will rapidly increase for improving treatment and settle the basis for personalized medicine. Pharmacogenetic profiles will be useful for prevention programs, and personalized therapies will reduce health costs and ensure a more efficient health care.

 Conclusion

The population studied was not significantly different in allele distribution with previously reported data from HapMap. Genotypes in Colombian population are similar to other previously studied groups of healthy subjects. The use of genotyping pharmacogenetic polymorphisms will prevent toxicity and adverse effects in tamoxifen treatment (for example in CYP2D6 rs3892097), and therapeutic alternatives should be evaluated based on individual pharmacogenetic studies.

 REFERENCES 1. C. Romero-Sánchez, A. Gómez-Gutiérrez, P.E. Gómez, M.C. Casas-Gómez, I. Briceño, C677T (RS1801133) MTHFR gene polymorphism frequency in a colombian population, Colombia Médica, 46(2), 75-79 (2015). Romero-Sánchez C. Gómez-Gutiérrez A. Gómez P.E. Casas-Gómez M.C. Briceño I. C677T (RS1801133) MTHFR gene polymorphism frequency in a colombian population Colombia Médica 46 2 75 79 2015 2. O. Fletcher, F. Dudbridge, Candidate gene-environment interactions in breast cancer, BMC Medicine, 12, 195 (2014). Fletcher O. Dudbridge F. Candidate gene-environment interactions in breast cancer BMC Medicine 12 195 195 2014 3. Y. Sun, C. Ye, X. Guo et al ., Evaluation of potential regulatory function of breast cancer risk locus at 6q25.1, Carcinogenesis, 37(2), 163-168 (2016). Sun Y. Ye C. Guo X. Evaluation of potential regulatory function of breast cancer risk locus at 6q25.1 Carcinogenesis 37 2 163 168 2016 4. F. Cui, D. Wu, W. Wang, X. He, M. Wang, Variants of FGFR2 and their associations with breast cancer risk: a HUGE systematic review and meta-analysis. Breast Cancer Res. Treat ., 155(2), 313-335 (2016). Cui F. Wu D. Wang W. X. He Wang M. Variants of FGFR2 and their associations with breast cancer risk: a HUGE systematic review and meta-analysis Breast Cancer Res. Treat 155 2 313 335 2016 5. E. Sirachainan, S. Jaruhathai, N. Trachu et al ., CYP2D6 polymorphisms influence the efficacy of adjuvant tamoxifen in Thai breast cancer patients, Pharmgenomics Pers. Med ., 5, 149-153 (2012). Sirachainan E. Jaruhathai S. Trachu N. CYP2D6 polymorphisms influence the efficacy of adjuvant tamoxifen in Thai breast cancer patients Pharmgenomics Pers. Med 5 149 153 2012 6. A. Cox, A.M. Dunning, M. García-Closas et al ., A common coding variant in CASP8 is associated with breast cancer risk, Nat. Genet ., 39(3), 352-358 (2007). Cox A. Dunning A.M. García-Closas M. A common coding variant in CASP8 is associated with breast cancer risk Nat. Genet 39 3 352 358 2007 7. T.N. Sergentanis, K.P. Economopoulos, Association of two CASP8 polymorphisms with breast cancer risk: A meta-analysis, Breast Cancer Res. Treat ., 120(1), 229-234 (2010). Sergentanis T.N. Economopoulos K.P. Association of two CASP8 polymorphisms with breast cancer risk: A meta-analysis Breast Cancer Res. Treat 120 1 229 234 2010 8. W. Zheng, J. Long, Y.T. Gao et al ., Genome-wide association study identifies a new breast cancer susceptibility locus at 6q25.1, Nat. Genet ., 41(3), 324-328 (2009). Zheng W. Long J. Gao Y.T. Genome-wide association study identifies a new breast cancer susceptibility locus at 6q25.1 Nat. Genet 41 3 324 328 2009 9. Q. Cai, W. Wen, S. Qu et al ., Replication and functional genomic analyses of the breast cancer susceptibility locus at 6q25.1 generalize its importance in women of chinese, Japanese, and European ancestry, Cancer Res ., 71(4), 1344-1355 (2011). Cai Q. Wen W. Qu S. Replication and functional genomic analyses of the breast cancer susceptibility locus at 6q25.1 generalize its importance in women of chinese, Japanese, and European ancestry Cancer Res 71 4 1344 1355 2011 10. G. Li, Y.B. Xiang, R. Courtney et al ., Association of a single nucleotide polymorphism at 6q25.1,rs2046210, with endometrial cancer risk among Chinese women. Chin. J. Cancer, 30(2), 138-143 (2011). Li G. Xiang Y.B. Courtney R. Association of a single nucleotide polymorphism at 6q25.1,rs2046210, with endometrial cancer risk among Chinese women Chin. J. Cancer 30 2 138 143 2011 11. M. García-Closas , F.J. Couch, S. Lindstrom et al ., Genome-wide association studies identify four ER negative-specific breast cancer risk loci, Nat. Genet ., 45(4), 392-398 (2013). García-Closas M. Couch F.J. Lindstrom S. Genome-wide association studies identify four ER negative-specific breast cancer risk loci Nat. Genet 45 4 392 398 2013 12. H. Liang, X. Yang, L. Chen et al ., Heterogeneity of breast cancer associations with common genetic variants in FGFR2 according to the intrinsic subtypes in Southern Han Chinese women, BiomedRes. Int ., 2015, 626948 (2015). Liang H. Yang X. Chen L. Heterogeneity of breast cancer associations with common genetic variants in FGFR2 according to the intrinsic subtypes in Southern Han Chinese women BiomedRes. Int 2015 626948 626948 2015 13. S. Siddiqui, S. Chattopadhyay, M.S. Akhtar et al ., A study on genetic variants of Fibroblast growth factor receptor 2 (FGFR2) and the risk of breast cancer from North India, PLoS One, 9(10), e110426 (2014). Siddiqui S. Chattopadhyay S. Akhtar M.S. A study on genetic variants of Fibroblast growth factor receptor 2 (FGFR2) and the risk of breast cancer from North India PLoS One 9 10 e110426 2014 14. L. Fejerman, M.C. Stern, E.M. John et al ., Interaction between common breast cancer susceptibility variants, genetic ancestry, and nongenetic risk factors in Hispanic women, Cancer Epidemiol. Biomarkers Prev ., 24(11), 1731-1738 (2015). Fejerman L. Stern M.C. John E.M. Interaction between common breast cancer susceptibility variants, genetic ancestry, and nongenetic risk factors in Hispanic women Cancer Epidemiol. Biomarkers Prev 24 11 1731 1738 2015 15. Z. Zhang, A. Zhu, B. Zhang et al ., Frequent mutation of rs13281615 and its association with PVT1 expression and cell proliferation in breast cancer, J. Genet. Genomics, 41(4), 187-195 (2014). Zhang Z. Zhu A. Zhang B. Frequent mutation of rs13281615 and its association with PVT1 expression and cell proliferation in breast cancer J. Genet. Genomics 41 4 187 195 2014 16. D. Campa, M. Barrdahl, M.M. Gaudet et al ., Genetic risk variants associated with in situ breast cancer, Breast Cancer Res ., 17, 82 (2015). Campa D. Barrdahl M. Gaudet M.M. Genetic risk variants associated with in situ breast cancer Breast Cancer Res 17 82 82 2015 17. Q. Chen, Y. He, C. Liu et al ., [Associations between TOX3 rs3803662 polymorphisms and immunological markers of breast cancer], Zhonghua YiXueZaZhi, 95(34), 2783-2786 (2015). Chen Q. He Y. Liu C. Associations between TOX3 rs3803662 polymorphisms and immunological markers of breast cancer Zhonghua YiXueZaZhi 95 34 2783 2786 2015 18. L. Dorling, G.C. Barnett, K. Michailidou et al ., Patients with a high polygenic risk of breast cancer do not have an increased risk of radiotherapy toxicity, Clin. Cancer Res ., 22(6), 1413-1420 (2016). Dorling L. Barnett G.C. Michailidou K. Patients with a high polygenic risk of breast cancer do not have an increased risk of radiotherapy toxicity Clin. Cancer Res 22 6 1413 1420 2016 19. R.L. Milne, B. Burwinkel, K. Michailidou et al ., Common non-synonymous SNPs associated with breast cancer susceptibility: Findings from the Breast Cancer Association Consortium, Hum. Mol. Genet ., 23(22), 6096-6111 (2014). Milne R.L. Burwinkel B. Michailidou K. Common non-synonymous SNPs associated with breast cancer susceptibility: Findings from the Breast Cancer Association Consortium Hum. Mol. Genet 23 22 6096 6111 2014 20. A. Rudolph, P.A. Fasching, S. Behrens et al ., A comprehensive evaluation of interaction between genetic variants and use of menopausal hormone therapy on mammographic density, Breast Cancer Res ., 17, 110 (2015). Rudolph A. Fasching P.A. Behrens S. A comprehensive evaluation of interaction between genetic variants and use of menopausal hormone therapy on mammographic density Breast Cancer Res 17 110 110 2015 21. D.M. Glubb, M.J. Maranian, K. Michailidou et al ., Fine-scale mapping of the 5q11.2 breast cancer locus reveals at least three independent risk variants regulating MAP3K1, Am. J. Hum. Genet ., 96(1), 5-20 (2015). Glubb D.M. Maranian M.J. Michailidou K. Fine-scale mapping of the 5q11.2 breast cancer locus reveals at least three independent risk variants regulating MAP3K1 Am. J. Hum. Genet 96 1 5 20 2015 22. L. Quan, C.C. Hong, G. Zirpoli et al ., Variants of estrogen-related genes and breast cancer risk in European and African American women, Endocr. Relat. Cancer, 21(6), 853-864 (2014). Quan L. Hong C.C. Zirpoli G. Variants of estrogen-related genes and breast cancer risk in European and African American women Endocr. Relat. Cancer 21 6 853 864 2014 23. R.L. Milne , J. Lorenzo-Bermejo, B. Burwinkel et al ., 7q21-rs6964587 and breast cancer risk: an extended case-control study by the Breast Cancer Association Consortium, J. Med. Genet ., 48(10), 698-702 (2011). Milne R.L. Lorenzo-Bermejo J. Burwinkel B. 7q21-rs6964587 and breast cancer risk: an extended case-control study by the Breast Cancer Association Consortium J. Med. Genet 48 10 698 702 2011

No potential conflict of interest was reported by the authors.

Y.V. Ariza-Márquez, F.A. Aristizábal-Gutiérrez, A. Gómez-Gutiérrez, P.E. Gómez, M.C. Casas-Gómez, I. Briceño-Balcázar, Frequency of 14 genetic variants associated with breast cancer risk and treatment in a Colombian population, Rev. Colomb. Cienc. Quím. Farm., 47(2), 277-288 (2018).