Genetic polymorphisms of estrogen receptor genes are associated with breast cancer susceptibility in Chinese women

Background 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. Methods 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. Results 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. Conclusions 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. Electronic supplementary material The online version of this article (10.1186/s12935-019-0727-z) contains supplementary material, which is available to authorized users.


Background
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 [1]. In Chinese women, BC also holds the highest incidence, accounting for 15% of all new cancers and 7% malignancy deaths [2].
BC is a complex heterogeneous disease. Both genetic and environmental factors are involved in the occurrence of BC [3].
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 [4]. 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 [5]. The TT/TC genotype of rs2234693 was related to high risk of BC [8]. 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 [6].
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.

Study population
BC patients treated at the Department of Oncology, the Second Affiliated Hospital of Xi'an Jiaotong University (Xi'an, China) were enrolled as cases [9]. 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).

Statistical analysis
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 Falsepositive 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  (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).

Discussion
Genetic variation is an important risk factor for BC [14]. 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 [4]. 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 [20]. 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 [21]. 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 [7]. 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.

Conclusions
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

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
Not applicable.

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