- Primary research
- Open Access
Genetic polymorphisms of estrogen receptor genes are associated with breast cancer susceptibility in Chinese women
Cancer Cell Internationalvolume 19, Article number: 11 (2019)
Estrogen exposure is a widely known risk factor for BC. And the interaction of estrogen with estrogen receptor (ER) plays an important role in breast cancer development. This case–control study aims to assess the association of genetic polymorphisms in the estrogen receptor genes with breast cancer (BC) susceptibility in Chinese Han women.
Four polymorphisms (rs2881766, rs9383951, rs9340799 in ESR1 and rs3020449 in ESR2) were genotyped in 459 patients and 549 healthy controls using the Sequenom MassARRAY method. Odds ratio (OR) and 95% confidence intervals (95% CI) were calculated to evaluate the associations. False-positive report probability (FPRP) was utilized to examine the noteworthiness of significant findings.
We observed that rs2881766 was associated with a decreased BC risk (GG vs. TT: OR = 0.63, 95% CI = 0.44–0.91; GG vs. TT/GT: OR = 0.68, 95% CI = 0.49–0.95), while rs3020449 was associated with an increased risk of BC (CT vs. TT: OR = 1.58, 95% CI = 1.21–2.06; CT/CC vs. TT: OR = 1.54, 95% CI = 1.20–1.98; TT/CC vs. CT: OR = 1.48, 95% CI = 1.15–1.90). The other two polymorphisms have no relation with BC susceptibility. In addition, rs2881766 was correlated with lymph node metastasis and ER expression, and rs3020449 was related to tumor size, histological grade and ER expression. The values of false-positive report probability indicated that the significant associations of BC risk with both rs2881766 and rs3020449 were noteworthy.
Our study suggests that polymorphisms rs2881766 and rs3020449 in estrogen receptor genes were associated with BC susceptibility as well as clinical features in Chinese women. These findings need further validation in a large population.
Breast cancer (BC), which had approximately 1.7 million new cases and 0.5 million deaths in 2012, has become the most common cancer type and the leading cause of global cancer death among females . In Chinese women, BC also holds the highest incidence, accounting for 15% of all new cancers and 7% malignancy deaths . BC is a complex heterogeneous disease. Both genetic and environmental factors are involved in the occurrence of BC .
Estrogen exposure is a widely known risk factor for BC. And the interaction of estrogen with estrogen receptor (ER), which can alter the expression of downstream genes, plays an important role in breast cancer development . ER alpha (ERα) and beta (ERβ) are the main forms of ER. They were encoded by the gene estrogen receptor 1 (ESR1) and estrogen receptor 2 (ESR2) respectively. And evidence showed that genetic variants in these two genes were associated with BC susceptibility [5,6,7,8]. Polymorphisms in ESR1 have been investigated in a number of studies. The pooled results demonstrated that the AA genotype of rs2228480 was correlated with a lower risk of BC in Caucasians. The C allele of the rs3798577 decreased BC risk in Asians while increased BC risk in Caucasians . The TT/TC genotype of rs2234693 was related to high risk of BC . In addition, meta-analysis of the relationship between ESR2 polymorphisms and BC susceptibility showed that the polymorphism rs4986938 was associated with reduced BC risk in overall population under both dominant and heterozygous models .
In the present study, four single nucleotide polymorphisms (SNPs) in the two estrogen receptor genes were selected to be studied. The effects of these SNPs (rs2881766, rs9383951, rs9340799 in ESR1 and rs3020449 in ESR2) in BC risk in Chinese population were either seldom being explored or the conclusions were still controversial. Therefore, we conducted a hospital-based case–control study to explore the association of the four polymorphisms with BC susceptibility in Chinese Han women.
BC patients treated at the Department of Oncology, the Second Affiliated Hospital of Xi’an Jiaotong University (Xi’an, China) were enrolled as cases . Women who came to the hospital for a checkup during the same period of time were recruited to form the control group. The cases were newly diagnosed and confirmed by histology or pathology. None of them had received chemotherapy, endocrinotherapy or radiotherapy before surgery. Patients with other types of cancer were excluded. The controls were healthy individuals without any history of a tumor or chronic diseases. All of the participants were Han population from Northwest China and have no relation to each other. Clinical information was collected from the medical records of the study subjects.
Sample collection and genotyping
Peripheral blood samples were collected in tubes coated with EDTA and were stored at − 80 °C. Genomic DNA was extracted from whole blood using the Universal Genomic DNA Extraction Kit (TaKaRa Bio Inc., Japan). DNA concentration was measured with the UV/VIS spectrophotometer (DU530, Beckman Instruments, USA). Four tag-SNPs (rs2881766, rs9383951, rs9340799 and rs3020449) were selected from The Single Nucleotide Polymorphism database (dbSNP). A multiplexed SNP MassEXTEND assay was designed by Sequenom MassARRAY Assay Design V3.0 Software (Agena Bioscience Inc., USA). And SNP genotypes was detected using Sequenom MassARRAY RS1000 [10, 11]. The primers of the selected SNPs were listed in Additional file 1: Table S1. Data was analyzed with Sequenom Typer Software (version 4.0, USA).
The Student t test and χ2-test were adopted to compare the differences in basic characteristics between cases and controls. Hardy–Weinberg equilibrium (HWE) for selected SNPs in controls and the differences in genotypes distribution between cases and controls were examined by χ2-test. We evaluated the associations of the four SNPs with BC risk in codominant, dominant, recessive, and overdominant models. Odds ratios (ORs) with 95% confidence intervals (CIs) were calculated by multivariate logistic regression, with adjustment for age and body mass index (BMI). All the tests were two-tailed and P-value < 0.05 was considered statistically significant. Akaike’s information criterion (AIC) and Bayesian information criterion (BIC) values were used to determine the best-fit model for each SNP. We then conducted False-positive report probability (FPRP) analysis to examine whether the significant results were noteworthy. The prior probability of 0.1 was set to detect the noteworthiness for OR of 1.50/0.67 and 0.25 was determined as a FPRP cut-off value as described in previous studies [12, 13]. All data were analyzed using SPSS software (version 22.0, IBM Corporation, USA).
Basic information of the study population and SNPs
A total of 1008 subjects including 459 BC patients and 549 healthy individuals were enrolled in this study. As shown in Table 1, no significant difference was observed between the case and control groups in the distribution of age, menopausal status, procreative times and BMI. It suggested that the cases and controls in this study were matched adequately on basic characteristics. The basic information of the four SNPs was presented in Table 2. The genotypic frequencies of all the selected SNPs in controls were in accordance with HWE. The frequency distribution of clinicopathological features in BC patients was shown in Additional file 1: Table S2.
Association between ESR1 and ESR2 SNPs with BC risk
Among the three ESR1 polymorphisms, rs2881766 was found to reduce BC risk under homozygous model (GG vs. TT: OR = 0.63, 95% CI = 0.44–0.91) and recessive model (GG vs. TT/GT: OR = 0.68, 95% CI = 0.49–0.95). The model with the lowest AIC and BIC values was considered as the best-fit model. For rs2881766, the recessive model might serve as the best-fit model. The other two SNPs (rs9383951 and rs9340799) were not related to BC susceptibility. As for the ESR2 polymorphism rs3020449, an increased risk of BC was found in heterozygous model (CT vs. TT: OR = 1.58, 95% CI = 1.21–2.06), dominant model (CT/CC vs. TT: OR = 1.54, 95% CI = 1.20–1.98), and overdominant model (TT/CC vs. CT: OR = 1.48, 95% CI = 1.15–1.90). And the AIC and BIC values suggested the dominant model may be the best-fit model for rs3020449 (Table 3).
FPRP analysis results
FPRP was adopted to assess the noteworthiness of the significant associations between the selected SNPs and BC risk. At the prior probability of 0.1 and FPRP cut-off value of 0.25, the associations between ESR2 rs3020449 and BC risk remained noteworthy in the three models (FPRP = 0.018, 0.016, and 0.034 respectively). In addition, the significant decrease of BC risk in carrier of ESR1 rs2881766 GG genotype was noteworthy (FPRP = 0.250 under GG vs.TT) (Table 4).
Relationship of rs2881766 and rs3020449 with clinicopathological features of BC
We further explored the relationship of ESR1 rs2881766 and ESR2 rs3020449 with BC clinicopathological features. The results showed that the GG genotype of rs2881766 was negatively correlated with lymph node metastasis (GG vs.TT: OR = 0.47, 95% CI = 0.27–0.83; GG vs.TT/GT: OR = 0.54, 95% CI = 0.32–0.91) and ER+ status (GG vs.TT: OR = 0.52, 95% CI = 0.29–0.91; GG vs.TT/GT: OR = 0.53, 95% CI = 0.32–0.90) (Table 5).
The rs3020449 was related to a larger tumor size under heterozygous and dominant model (CT vs. TT: OR = 1.60, 95% CI = 1.06–2.43; CT/CC vs. TT: OR = 1.55, 95% CI = 1.05–2.31). In addition, the CC genotype of rs3020449 was associated with higher histological grade (CC vs. TT: OR = 1.90, 95% CI = 1.02–3.54) as well as ER positive status (CC vs. TT: OR = 2.16, 95% CI = 1.12–4.17; CC vs. TT/CT: OR = 1.96, 95% CI = 1.05–3.63) (Table 6).
Genetic variation is an important risk factor for BC . Single nucleotide polymorphism is the most frequent variation in the genome. Genome-wide association study (GWAS) has identified numerous BC susceptibility loci in tumor-related genes such as FGFR2, TOX3, TP53, PTEN, MAP3K1, c-MYC, LSP1, and CASP8 [14, 15]. Estrogen plays a critical role in the development of breast cancer . Its effect on the breast epithelium is mediated by ER. The two major forms of ER (ERα and ERβ) are encoded by two separate genes, ESR1 located on Chr6 and ESR2 located on Chr14. And genetic polymorphisms in these two genes were reported to associate with BC susceptibility, including the extensively studied SNPs such as rs2228480, rs3798577, rs2077647, rs2234693 in ESR1 and rs4986938, rs1256049 in ESR2 as well as less common SNPs such as rs3020314, rs1514348, rs1062577, rs1271572, rs1256054 and rs2987983 [5,6,7,8, 16,17,18,19].
In this study, we investigated the associations of BC risk with four SNPs in the estrogen receptor genes ESR1 and ESR2. SNP rs2881766 is located in the promoter region of ESR1. This polymorphism was reported to increase BC risk in Korean women in a previously . However, another study conducted in Chinese population showed that the GG genotype of this SNP decreases BC risk in menarche > 13-year-old while increases BC risk in menarche ≤ 13-year-old subgroup . Our results indicated that the carriers of rs2881766 GG genotype had a lower risk of BC compared with the TT and GT genotype carriers. The discrepancy between our results and previous studies may caused by different sample sizes and the effects of confounding factors such as age, ethnicity and environmental effects. Moreover, the GG genotype may be a protective factor of lymph node metastasis for BC patients. And patients with rs2881766 GG genotype are more likely to have less expression of ER. SNP rs9340799 and rs9383951 are located in the first and fifth intron of ESR1 respectively [8, 22]. The association between rs9340799 (also known as Xbal) and BC risk has been evaluated in several studies before, but the conclusions were debatable. A few studies suggested this polymorphism can affect BC risk [23,24,25], whereas others did not find any relationship between this polymorphism and BC susceptibility [8, 18, 26, 27]. There is only one previous study of rs9383951and BC risk, but no significant association was found. Our study indicated that neither rs9340799 nor rs9383951 were related to BC risk, which is in line with most of previous studies. SNP rs3020449 is located in the promoter region of the ESR2 gene . Two previous studies investigated the effect of this SNP on BC risk, but neither one observed significant association [19, 28]. In contrast, our results suggested that the CT genotype of rs3020449 increased the risk of BC. In addition, BC patients with CT genotype tend to have a larger tumor size compared with TT genotype carriers. Furthermore, patients with CC genotype may have a greater tumor grade and higher expression of ER.
FPRP analysis is an effective approach to verify the noteworthiness of significant findings. In this study, we adopted a relatively stringent cut-off value for FPRP. The FPRP value of the significant association between ESR2 rs3020449 and BC risk was much lower than the threshold, suggesting that our findings of this SNP were noteworthy and authentic. The significant association of ESR1 rs2881766 with BC risk was noteworthy only in homozygous model. But actually, FPRP value of significant association in another model was still quite small (< 0.3). Hence, we believe our findings for this polymorphism were credible to some extent.
In conclusion, the current study suggests that the ESR1 rs2881766 decreases BC risk while ESR2 rs3020449 increases BC risk in Chinese women. Future large studies are required to validate these findings. The possible mechanisms underlying the associations also need to be explored.
single nucleotide polymorphism
Akaike’s information criterion
Bayesian information criterion
false-positive report probability
Torre LA, Bray F, Siegel RL, Ferlay J, Lortet-Tieulent J, Jemal A. Global cancer statistics, 2012. Cancer J Clin. 2015;65(2):87–108.
Chen W, Zheng R, Baade PD, Zhang S, Zeng H, Bray F, Jemal A, Yu XQ, He J. Cancer statistics in China. Cancer J Clin. 2016;66(2):115–32.
Fletcher O, Dudbridge F. Candidate gene-environment interactions in breast cancer. BMC Med. 2014;12:195.
Yager JD, Davidson NE. Estrogen carcinogenesis in breast cancer. N Engl J Med. 2006;354(3):270–82.
Li T, Zhao J, Yang J, Ma X, Dai Q, Huang H, Wang L, Liu P. A meta-analysis of the association between ESR1 genetic variants and the risk of breast cancer. PLoS ONE. 2016;11(4):e0153314.
Yu KD, Rao NY, Chen AX, Fan L, Yang C, Shao ZM. A systematic review of the relationship between polymorphic sites in the estrogen receptor-beta (ESR2) gene and breast cancer risk. Breast Cancer Res Treat. 2011;126(1):37–45.
Jahandoost S, Farhanghian P, Abbasi S. The effects of sex protein receptors and sex steroid hormone gene polymorphisms on breast cancer Risk. J Natl Med Assoc. 2017;109(2):126–38.
Zhang Y, Zhang M, Yuan X, Zhang Z, Zhang P, Chao H, Jiang L, Jiang J. Association between ESR1 PvuII, XbaI, and P325P polymorphisms and breast cancer susceptibility: a meta-analysis. Med Sci Monit. 2015;21:2986–96.
Liu D, Wang M, Tian T, Wang XJ, Kang HF, Jin TB, Zhang SQ, Guan HT, Yang PT, Liu K, et al. Genetic polymorphisms (rs10636 and rs28366003) in metallothionein 2A increase breast cancer risk in Chinese Han population. Aging. 2017;9(2):547–55.
Gabriel S, Ziaugra L, Tabbaa D: SNP genotyping using the Sequenom MassARRAY iPLEX platform. Curr Prot Hum Genet 2009, Chapter 2:Unit 2.12.
Wang M, Tian T, Ma X, Zhu W, Guo Y, Duan Z, Fan J, Lin S, Liu K, Zheng Y, et al. Genetic polymorphisms in caveolin-1 associate with breast cancer risk in Chinese Han population. Oncotarget. 2017;8(53):91654–61.
Wacholder S, Chanock S, Garcia-Closas M, El Ghormli L, Rothman N. Assessing the probability that a positive report is false: an approach for molecular epidemiology studies. J Natl Cancer Inst. 2004;96(6):434–42.
Lai MC, Yang Z, Zhou L, Zhu QQ, Xie HY, Zhang F, Wu LM, Chen LM, Zheng SS. Long non-coding RNA MALAT-1 overexpression predicts tumor recurrence of hepatocellular carcinoma after liver transplantation. Med Oncol (Northwood, London, England). 2012;29(3):1810–6.
Mavaddat N, Antoniou AC, Easton DF, Garcia-Closas M. Genetic susceptibility to breast cancer. Mol Oncol. 2010;4(3):174–91.
Mahdi KM, Nassiri MR, Nasiri K. Hereditary genes and SNPs associated with breast cancer. Asian Pac J Cancer Prev. 2013;14(6):3403–9.
Lipphardt MF, Deryal M, Ong MF, Schmidt W, Mahlknecht U. ESR1 single nucleotide polymorphisms predict breast cancer susceptibility in the central European Caucasian population. Int J Clin Exp Med. 2013;6(4):282–8.
Dehghan Z, Sadeghi S, Tabatabaeian H, Ghaedi K, Azadeh M, Fazilati M, Bagheri F. ESR1 single nucleotide polymorphism rs1062577 (c.*3804T > A) alters the susceptibility of breast cancer risk in Iranian population. Gene. 2017;611:9–14.
Ghali R, Al-Mutawa MA, Al-Ansari AK, Zaied S, Bhiri H, Mahjoub T, Almawi WY. Differential association of ESR1 and ESR2 gene variants with the risk of breast cancer and associated features: a case–control study. Gene. 2018;651:194–9.
Chen L, Liang Y, Qiu J, Zhang L, Chen X, Luo X, Jiang J. Significance of rs1271572 in the estrogen receptor beta gene promoter and its correlation with breast cancer in a southwestern Chinese population. J Biomed Sci. 2013;20:32.
Son BH, Kim MK, Yun YM, Kim HJ, Yu JH, Ko BS, Kim H, Ahn SH. Genetic polymorphism of ESR1 rs2881766 increases breast cancer risk in Korean women. J Cancer Res Clin Oncol. 2015;141(4):633–45.
Chen L, Kang H, Jin GJ, Chen X, Zhang QY, Lao WT, Li R. The association between a novel polymorphism (rs1062577) in ESR1 and breast cancer susceptibility in the Han Chinese women. Gynecol Endocrinol. 2016;32(7):553–6.
Long J, Cai Q, Sung H, Shi J, Zhang B, Choi JY, Wen W, Delahanty RJ, Lu W, Gao YT, et al. Genome-wide association study in east Asians identifies novel susceptibility loci for breast cancer. PLoS Genet. 2012;8(2):e1002532.
Shin A, Kang D, Nishio H, Lee MJ, Park SK, Kim SU, Noh DY, Choe KJ, Ahn SH, Hirvonen A, et al. Estrogen receptor alpha gene polymorphisms and breast cancer risk. Breast Cancer Res Treat. 2003;80(1):127–31.
Wang J, Higuchi R, Modugno F, Li J, Umblas N, Lee J, Lui LY, Ziv E, Tice JA, Cummings SR, et al. Estrogen receptor alpha haplotypes and breast cancer risk in older Caucasian women. Breast Cancer Res Treat. 2007;106(2):273–80.
Madeira KP, Daltoe RD, Sirtoli GM, Carvalho AA, Rangel LB, Silva IV. Estrogen receptor alpha (ERS1) SNPs c454-397T > C (PvuII) and c454-351A > G (XbaI) are risk biomarkers for breast cancer development. Mol Biol Rep. 2014;41(8):5459–66.
Sun MY, Du HY, Zhu AN, Liang HY, de Garibay GR, Li FX, Li M, Yang XX. Genetic polymorphisms in estrogen-related genes and the risk of breast cancer among Han Chinese women. Int J Mol Sci. 2015;16(2):4121–35.
Li N, Dong J, Hu Z, Shen H, Dai M. Potentially functional polymorphisms in ESR1 and breast cancer risk: a meta-analysis. Breast Cancer Res Treat. 2010;121(1):177–84.
Treeck O, Elemenler E, Kriener C, Horn F, Springwald A, Hartmann A, Ortmann O. Polymorphisms in the promoter region of ESR2 gene and breast cancer susceptibility. J Steroid Biochem Mol Biol. 2009;114(3–5):207–11.
ZJD conceived and designed the study. MW, HTL, SL, and QH collected samples. TLY provided technical support. MW, PX, YJD, and YD collected and processed data. TT, ZJD and YD analyzed data. LHZ and NL prepared tables. TT drafted the manuscript. ZJD and YD revised the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Availability of data and materials
All data generated and analyzed during this study are included in this published article and its additional file.
Consent for publication
Ethics approval and consent to participate
The current study was approved by the Ethics Committee of Xi’an Jiaotong University. A written consent from each participant was obtained after they were informed the purpose of this study. The research protocol was in accordance with the approved guidelines.
This study was supported by National Natural Science Foundation, People’s Republic of China (No. 81471670), the Key research and development plan, Shaanxi Province, People’s Republic of China (2017ZDXM-SF-066), Science and technology branch project of Xinjiang Uygur Autonomous Region, People’s Republic of China (2017E0262).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.