- Research article
- Open Access
Polymorphisms in the estrogen receptor alpha gene (ESR1), daily cycling estrogen and mammographic density phenotypes.
BMC Cancer volume 16, Article number: 776 (2016)
Single nucleotide polymorphisms (SNPs) involved in the estrogen pathway and SNPs in the estrogen receptor alpha gene (ESR1 6q25) have been linked to breast cancer development, and mammographic density is an established breast cancer risk factor. Whether there is an association between daily estradiol levels, SNPs in ESR1 and premenopausal mammographic density phenotypes is unknown.
We assessed estradiol in daily saliva samples throughout an entire menstrual cycle in 202 healthy premenopausal women in the Norwegian Energy Balance and Breast Cancer Aspects I study. DNA was genotyped using the Illumina Golden Gate platform. Mammograms were taken between days 7 and 12 of the menstrual cycle, and digitized mammographic density was assessed using a computer-assisted method (Madena). Multivariable regression models were used to study the association between SNPs in ESR1, premenopausal mammographic density phenotypes and daily cycling estradiol.
We observed inverse linear associations between the minor alleles of eight measured SNPs (rs3020364, rs2474148, rs12154178, rs2347867, rs6927072, rs2982712, rs3020407, rs9322335) and percent mammographic density (p-values: 0.002–0.026), these associations were strongest in lean women (BMI, ≤23.6 kg/m2.). The odds of above-median percent mammographic density (>28.5 %) among women with major homozygous genotypes were 3–6 times higher than those of women with minor homozygous genotypes in seven SNPs. Women with rs3020364 major homozygous genotype had an OR of 6.46 for above-median percent mammographic density (OR: 6.46; 95 % Confidence Interval 1.61, 25.94) when compared to women with the minor homozygous genotype. These associations were not observed in relation to absolute mammographic density. No associations between SNPs and daily cycling estradiol were observed. However, we suggest, based on results of borderline significance (p values: 0.025–0.079) that the level of 17β-estradiol for women with the minor genotype for rs3020364, rs24744148 and rs2982712 were lower throughout the cycle in women with low (<28.5 %) percent mammographic density and higher in women with high (>28.5 %) percent mammographic density, when compared to women with the major genotype.
Our results support an association between eight selected SNPs in the ESR1 gene and percent mammographic density. The results need to be confirmed in larger studies.
Genetic factors are believed to account for 30–60 % of the variance in mammographic density [1, 2] while established breast cancer risk factors such as age, body mass index (BMI), parity, age at first birth, the use of hormone therapy and physical activity account for the remainder of its variability . Substantial evidence supports the role of ovarian steroid hormones in breast cancer , and estrogens are associated with breast cancer development in both pre- and postmenopausal women [5–7]. Estrogens increase cellular proliferation in breast tissue, which may increase mammographic density . Mammographic density reflects the proportion of fibroglandular cells in the breast tissue; higher density indicates increased potential for proliferative activity  and is an established risk factor for breast cancer [10, 11]. Recently, studies of associations between breast cancer risk factors and mammographic density has focused not only on percent mammographic density [12, 13], but on various mammographic density phenotypes [14–17] as absolute mammographic density is believed to represent the actual target tissue for tumor development [18–21].
The genetic determinants of mammographic density have yet to be identified. A recent combined meta analysis of data from five genome wide studies (GWAS) among women of European decent suggests that multiple loci might be involved [22, 23]. Given the associations between estrogens and breast cancer, it is plausible that genetic variation in estrogen receptors may be important.
Both estrogen receptors, alpha (ERα) and beta (ERβ), are members of the nuclear receptor family of ligand-inducible transcription factors, but they are believed to have different transcriptional activation properties . It has been assumed for a long time that the action of estrogens in carcinogenesis is via the ERα signaling and the proliferation induced by the interaction between estrogens and estrogen receptors. [4, 24] Thus, the estrogen receptor alpha gene (ESR1), which encodes ERα, is of special interest in relation to breast cancer initiation, development, and therapeutics. Specific polymorphisms (SNPs) in ESR1 may directly or indirectly lead to variations in its activity, and may have an effect on breast cancer risk. Although a number of studies have examined SNPs located in estrogen related genes and possible associations with mammographic density [25–36], few have focused on premenopausal women [25, 27, 28, 33, 36], and only a handful have considered SNPs in the ESR1 (rs2234693 and rs9340799) [25, 26, 29–31]. Further, none have chacterized participants by menstrual cycle hormonal milieu.
Previously, in the Norwegian Energy Balance and Breast Cancer Aspects (EBBA)-I study, we observed a positive association between daily circulating ovarian sex hormones and both absolute and percent mammographic density [15, 37]. In addition, daily 17β-estradiol profiles were positively associated with traditional breast cancer risk factors, such as early age at menarche , short time since last birth , and unfavourable metabolic profile [38, 40, 41]. We also observed associations between two SNPs in the estrogen pathway (rs7172156 and rs749292 in CYP19A1), daily cycling 17β-estradiol, and absolute and percent mammographic density  as well as between CYP17 (rs2486758) and metabolic risk factors  These associations also point to the need for further studies of estrogen, mammographic density phenotypes and susceptibility genes in combination. Here, we expand upon those results to consider the ESR1 region, examining whether SNPs are associated with daily cycling estradiol and premenopausal mammographic density phenotypes.
Participants and study design
The Norwegian Energy Balance and Breast cancer Aspects I study (EBBA-I) included a total of 204 healthy premenopausal women, aged 25–35, recruited from the general population by announcements in local newspapers, and public meeting places. The study was conducted at the Department of Clinical Research at the University Hospital of Northern Norway (UNN), Tromsø, between 2000 and 2002 and has been described in detail elsewhere [14, 17, 38], we briefly summarize the methods. The participating women had to meet the following eligibility criteria: self-reported regular menstruation (cycle length: 22–38 days within the previous 3 months), no use of steroid contraceptives, pregnancy or lactation in the previous 6 months, no infertility, no history of gynecological disorders, and no chronic disorders (e.g. diabetes, hypo-/hyperthyroidism). At recruitment, subjects completed questionnaires and were interviewed by a trained nurse. Recall and memory-probing aids, including a lifetime calendar, were used to date specific life events [14, 17, 38]. These interviews including items on demographics, reproductive history, and lifestyle factors including: age at menarche, marital status, education, ethnicity, parity, physical activity, previous use of hormonal contraceptives, family history of cancer, smoking, and alcohol use [17, 38]. Birth weight obtained from the questionnaire and interview was also obtained by a linkage to the national Birth Registry. Two women were excluded from the current analyses due to missing mammographic data, resulting in 202 participants .
All participants underwent clinical examinations within three specified intervals during a single menstrual cycle: (1) between days 1–5 after onset of bleeding; (2) days 7–12; and (3) days 21–25. Height was measured to the nearest 0.5 cm and weight to the nearest 0.1 kg on a regularly calibrated electronic scale. Body mass index (BMI) for our analysis was calculated as weight in kilograms per height in square meter (kg/m2)  using data from the first visit. Waist circumference (WC) was measured to the nearest 0.5 cm, 2.5 cm above the umbilicus. During the second visit a whole body scan was obtained for the estimation of the total percentage of fat tissue, using dual-energy X-ray absorptiometry (DEXA; DPLX-L 2288, Lunar Radiation Corporation, Madison, WI, USA). The percentage of fat tissue was estimated using Lunar software [14, 17, 37].
Assessment of 17β-estradiol in serum and in saliva
At all three scheduled visits overnight fasting serum concentrations of 17β-estradiol were measured in fresh sera using a direct immunometric assay (Immuno-1; Bayer Diagnostics, Norway) at the Department of Clinical Chemistry, UNN . The sensitivity for estradiol was 0.01 nmol/L, and the coefficient of variation (CV) was 3.9 %.
The participants self-collected daily morning saliva samples during the course of a whole menstrual cycle in order to assess the daily bioavailable fraction of 17β-estradiol. Sampling started on the first morning of menstrual bleeding and was conducted according to previously established and validated collection protocols [38, 43]. Samples were sent to the Reproductive Ecology Laboratory at Harvard University where they were stored at −70 °C until analysis. 17β-estradiol concentrations were measured in each saliva sample using 125I- labeled RIA kits (#39100, Diagnostic Systems Laboratories, Webster, TX, USA). All samples were run in duplicate, and samples from a single participant were run within the same assay batch. CVs were calculated based on the high and low value pools included in each assay [14, 17, 37].
Following estradiol assay, all cycles were aligned based on the identification of the mid-cycle drop in salivary 17 β-estradiol concentration (hereafter designated cycle day 0), which provides a reasonable estimate of the day of ovulation [17, 43, 44]. Two measures of 17β-estradiol were calculated for all participants: mean cycle-long concentration and mean mid-cycle (day −7 to +6) concentrations, as well as using daily levels of salivary 17β-estradiol. The mid-cycle 17β-estradiol drop could not be identified for 14 women, hence their cycles could not be aligned and they were omitted from the statistical analysis .
Single-nucleotide polymorphism selection and genotyping
We analysed genetic polymorphisms in ESR1, 6q25, that encodes the estrogen receptor α. DNA was extracted from frozen whole blood using a MagAttract DNA Blood Mini M48 kit (QIAGEN, Valencia, CA, USA) by the Department of Medical Genetics, UNN. DNA was genotyped on the Golden Gate Platform (Illumina, San Diego, CA, USA) at the Fred Hutchinson Cancer Research Centre (Makar Lab), using the manufacturer’s protocol . These methods have previously been described in detail . In brief, 250 ng of genomic DNA was divided into aliquots in 96-well plates, processed accordingly and scanned on the Illumina iScan reader using BeadStudio software [17, 42].
We conducted a series of quality control procedures . SNP call rates exceeded 99 % for this study, with 100 % concordance of blinded duplicates. The linkage disequilibrium select algorithm was employed to choose the tag SNPs via the Genome Variation Server [46, 47]. The SNPs were selected using an r2 threshold of 0.8 and a minor allele frequency >5 %, representing variability in the white European population. Tag SNP coverage extended 2 kilobases (kb) upstream and 1 kb downstream of the gene, and 76 SNPs were covered. We further reduced the number of SNPs using power calculations and ended up with a final selection of 34 common SNPs with minor allele frequency >0.2. (Additional file 1: Table S1) None of the selected SNPs were monomorphic or significantly out of Hardy–Weinberg equilibrium .
Mammograms and mammographic density
We obtained bilateral two-view mammograms from our subjects during the second scheduled visit (between cycle days 7 and 12) at the Centre for Breast Imaging, UNN, using a standard protocol [17, 38, 48]. The left craniocaudal mammograms were digitised and imported into a computerised mammographic density assessment programme (Madena) developed at the University of Southern California School of Medicine (Los Angeles, CA, USA) [38, 49, 50]. A single trained reader (G. Ursin) conducted all density measurements as follows: A region of interest (ROI) that included the entire breast was identified, (excluding the pectoralis muscle, prominent veins and fibrous strands). A tinting tool was used to highlight pixels representing dense areas of the mammograms within the ROI. The size of these dense areas (in square centimetres) was automatically calculated by the Madena software, giving a measure of absolute mammographic breast density [15, 38]. We then calculated percentage mammographic density as the ratio of absolute mammographic breast density to total breast area multiplied by 100 . The mammograms were read in four batches, with an equal number of mammograms included in each batch. A duplicate reading of 26 randomly selected mammograms from two of the batches showed a Pearson’s correlation coefficient of 0.97. The reader was blinded to any characteristics of the study population [14, 17, 37, 38].
Descriptive characteristics were calculated by means (standard deviation) for continuous variables and percent for binary data. On the basis of the plausible biological mechanisms related to the estrogen metabolic pathway, we selected 34 SNPs in the ESR1, 6q25 gene for further analysis. These SNPs were coded as AA = 0 (major homozygous), Aa = 1 (heterozygous) and aa = 2 (minor homozygous) .
Mammographic density (absolute and percent) was considered continuously in our first set of models. We used multivariable linear regression models to assess the association between mammographic density phenotypes (absolute and percent mammographic density) as dependent variables and ESR1 SNPs as ordinal independent variables.
Percent mammographic density and absolute mammographic density were used as both continuous and dichotomized variables, representing lower and higher density, using median values as cut-off points; Percent mammographic density (28.5 %), and absolute mammographic density (32.4 cm2) . Previous studies of pre- and postmenopausal women have observed a 2–3 fold increase in breast cancer risk in women with percent mammographic density > 25 % [8, 51] and absolute mammographic density > 32 cm2. These observations along with our own observations related to estradiol levels [15, 37] support the comparison of women with above versus below median percent and absolute mammographic density . Thus, mammographic density outcome variables were also used as dichotomized variables in logistic regression models where indicator variables of each SNP was included using aa as the reference level.
We considered several potential confounding factors known to be associated with mammographic density phenotypes, estrogen concentrations, and/or ESR1 variant. These included age (continuous), BMI (continuous), birth weight (continuous), age at menarche (continuous), parity (categorical), previous oral contraceptive use (categorical) and current smoking habits (categorical) . Age, BMI, parity, current smoking habits, and previous oral contraceptive use were included as covariates in the final models. Salivary 17β-estradiol (continuous) and birth weight [13, 40], did not influence our estimates, and these were left out of our final model.
Based on plausible biological mechanisms and results from the first set of models, we selected eight SNPs (rs3020364, rs2474148, rs12154178, rs2347867, rs6927072, rs2982712, rs3020407, rs9322335) for further analyses in which we stratified the subjects based on median BMI (23.6 kg/m2). We again fitted multivariable linear regression models using the same set of covariates as previously specified, with the exception of BMI, to show how the relationship between percent mammographic density and the SNP genotypes might vary among different strata of women. More detailed stratification (i.e. to tertiles of body mass index) gave no additional information. Thus, BMI was used both as a continuous variable as well as a dichotomized one when we performed stratified analysis (comparing low BMI vs high BMI).
We used linear mixed models for repeated measures to study variations of daily salivary 17β-estradiol across the menstrual cycle, for subgroups of women with either major, minor homozygous or heterozygous genotypes for all eight SNPs, and adjusted for current smoking habits (yes/no), previous use of oral contraceptives (yes/no) as well as the same confounders as in the initial linear regression models, and stratified our data by median percent mammographic density (28.5 %). This revealed a suggested pattern for three of our eight SNPs (rs3020364, rs2474148, rs2982712).
All P-values were two-tailed and considered significant when the value was <0.05. The analyses were conducted with SPSS version 22.0 software (IBM, Armonk, NY, USA).
Table 1 provides selected general characteristics of the study participants. The participating women were on average 30.7 years old and had a mean BMI of 24.4 kg/m2. We found a mean salivary 17β-estradiol concentration of 17.9 pmol/l, mean percent mammographic density of 29.8 % (median 28.5 %), and mean absolute mammographic density of 34.7 cm2 (median 32.4 cm2). Common risk factors such as age, parity and body composition (BMI, WC, and total tissue fat) were inversely associated with mammographic density phenotypes. (Results not shown).
In our multiple linear regression models we found that having the minor homozygous allele was associated with decreased percent mammographic density for eight of our 34 SNPs: rs3020364 (p-value 0.01), rs2474148 (p-value 0.023), rs12154178 (p-value 0.026), rs2347867 (p-value 0.002), rs6927072 (p-value 0.007), rs2982712 p-value 0.006), rs3020407 (p-value 0.020), rs9322335 (p-value 0.024) (Table 2). For all but two (rs2347867, rs9322335) of these eight SNPs this association was seen in lean, but not in heavier women when we dichotomized our data by median split of BMI (Lean: ≤23.6, Heavy > 23.6). With the exception of rs9322335 (p-value 0.031), we did not observe similar associations in relation to absolute mammographic density (Additional file 2: Table S2).
The frequencies of genotypes of our selected SNPs in the study population were similar to those recorded in HapMap (Table 3). For seven out of eight SNPs, the odds of above-median percent mammographic density (>28.5 %) were 3–6 times higher among women with major homozygous genotypes compared to women with minor homozygous genotypes: rs3020364: OR 6.46 (p-value 0.009), rs2474148: OR 4.23 (p-value 0.028), rs12154178: OR 5.44 (p-value 0.014), rs2347867: OR 3.38 (p-value 0.030), rs6927072: OR 3.44 (p-value 0.028), rs2982712: OR 3.97 (p-value 0.013), rs3020407: OR 3.97 (p-value 0.024). When comparing the minor homozygous genotypes to the heterozygous, the odds of above median mammographic were higher for the heterozygous, but only significant for rs2982712 (OR: 3.48, p-value 0.014) (Table 3).
When examining the 17β-estradiol concentrations throughout the menstrual cycle for women with high and low percent mammographic density separately (median split of percent mammographic density) using linear mixed models we found a striking, albeit only borderline significant pattern for 3 of our 8 SNPs; the level of 17β-estradiol for women with the minor aa genotype for rs3020364, rs24744148 and rs2982712 were lower throughout the cycle for women with low (<28.5 %) percent mammographic density and higher throughout the cycle for women with high (>28.5 %) percent mammographic density, when compared to women with the major AA genotype (Fig. 1) We observed a 24.3 % lower level of mean 17β- estradiol throughout a menstrual cycle in women with low mammographic density (<28.5 %) with minor genotype aa of rs3020364 compared to those with major genotype AA (Fig. 1a, p-value 0.076) and in women with high mammographic density and minor genotype aa of the same SNP we observed a 58 % higher level of mean 17β-estradiol throughout the cycle compared with women with major genotype AA. (Fig. 1b: p-value 0.079). Fig. 1c-f show a similar pattern for rs2474148 and rs2982714, results in Fig. 1f are significant (p-value 0.025).
We have assessed the association between 34 polymorphisms in the ESR1 gene with mammographic density phenotypes and daily cycling estradiol in healthy, premenopausal, Norwegian women. Our results support an association between 8 (rs3020364, rs2474148, rs12154178, rs2347867, rs6927072, rs2982712, rs3020407, rs9322335) studied SNPs and percent mammographic density, and one SNP (rs9322335) also associated with absolute mammographic density. For all but two of these SNPs (rs2347867, rs9322335) this inverse linear association was maintained in lean (BMI ≤23.6 kg/m2), but not in heavier women (BMI >23.6 kg/m2). We also suggest, based on results of borderline significance that the level of 17β-estradiol for women with the minor aa genotype for rs3020364, rs24744148 and rs2982712 were lower throughout the cycle for women with low (<28.5 %) percent mammographic density and higher throughout the cycle for women with high (>28.5 %) percent mammographic density, when compared to women with the major AA genotype.
All 8 SNPs associated with mammographic density in the present study are intronic. It is not yet clear why and how noncoding SNPs influence the gene activity , but one suggested mechanism is that intronic SNPs may regulate gene expression through endogenous transacting factors, epigenetics and chromosome conformation . Previous GWAS studies have shown intronic SNPs to be important breast cancer risk loci . This does not imply causality, but may aid in the identification of novel susceptibility loci . The ERα-receptor mediates estrogen action by regulating gene transcription, triggering the mRNA-synthesis and consequently the production of proteins that produce the physiological hormone effects . The expression of ERα seems to be linked to breast cell proliferation as well as to be increased during tumorigenesis .
Since one of the mechanisms by which estrogen promotes the proliferation of both normal and neoplastic breast epithelial cells is through its binding to ERα, and that the content of ERα in breast tissue is associated with increased breast cell proliferation [56, 57], a modification of the expression of this gene could influence mammographic density or breast cancer risk. Indeed, a study among postmenopausal women found a statistically significant increased protein expression of ERα in normal epithelium of women with high mammographic density compared to those with low density . To our knowledge, there are no studies related to ESR1 SNPs, daily levels of estrogens throughout an entire menstrual cycle and mammographic density phenotypes in premenopausal women, and published results on the specific SNPs analysed in our current study do not exist.
The possible associations between SNPs in the ER genes and mammographic density are still under debate [25, 26, 29–31]. Some in vitro studies have shown that the expression of ERα in breast cancer cells may be down regulated by estradiol . Recently published results from a study of 3,872 common genetic variants across three ESR1 loci with subjects from three international consortia found evidence for at least five independent causal variants associated with phenotype sets, mammographic density being one. We found no evidence of correspondence between our SNPs and the causal variants identified, but this may be explained by the fact that due to our small sample size we identified 34 SNPs with a sufficient MAF >0.2, while the newly published article could include the genotypes of MAF > 0.02 .
Thus far, research has mostly focused on SNPs rs2234693 (PvuII) or rs9340799 (XbaI). The impact of PvuII and XbaI variant alleles on ESR1 transcription is not well established. A recent meta-analysis including 11 studies identified reduced breast cancer risk for those carrying homozygous the minor allele of SNP rs2234693, while no association was observed for SNP rs9340799 . Others have suggested that these variants are associated with reduced breast cancer risk among women with lower levels of estrogen [62, 63], and these findings complement those emerging from some of the studies on mammographic density: Among a population of pre- and postmenopausal women, the minor allele of the rs9340799 SNP was associated with decreased mammographic density, while no association was found for ERα rs2234693 . By contrast, others have found that premenopausal women with the homozygous minor allele of SNP rs2234693 (ERα) have higher mammographic density. Similar associations have been suggested for rs9340799 . It is, therefore, plausible that estrogen concentrations may influence the relationship between susceptibility genes and mammographic density. This suggests a biological rationale for our findings of different estrogen levels throughout the menstrual cycle of genotype aa compared to AA when our population was split by percent mammographic density (Fig. 1). Our results suggest that having the minor homozygous allele is associated with decreased percent mammographic density, and the associations are stronger in lean women. However, no association has been observed within strata of estrogen-related factors (parity, hormonal derivatives used, age at menarche, BMI) among premenopausal women , nor in populations of premenopausal and postmenopausal women ([26, 30, 31]. Other investigated SNPs (rs2228480, rs728524, rs3798577 and rs2077647) have shown no association with mammographic density in premenopausal and/or postmenopausal populations [25, 26].
Several of our findings are novel. To our knowledge, no other studies have focused on the associations between mammographic density phenotypes and these 8 SNPs in the ESR1 gene. In a large case–control study, 41 independent breast cancer susceptibility variants were discovered using a custom genotyping array designed, in part, by the Breast Cancer Association Consortium . Of these variants, a recent report found novel associations between breast cancer SNPs in 6q25 including rs9383938 in the ESR1 region with a volumetric measure of mammographic density in 5000 Swedish women . The marker of density consortium (MODE) recently identified nine loci associated with an area-based mammographic density one of which is ESR1 [66, 67].
Specific SNPs located in genes of this pathway could directly or indirectly lead to variations that may have effects on breast cancer risk. Also, estrogens are known to regulate the activity of several enzymes and receptors in the estrogen pathway [59, 68–70] Therefore, the relation between SNPs of such genes and breast cancer risk could vary among strata of women based on their estradiol levels . However, no overall linear association was found between our SNPs and estradiol.
The strengths of our study include daily measurements of salivary concentrations of unbound, bioactive 17β-estradiol collected across an entire menstrual cycle , as well as following strict procedures , and validated methods . Among ovulating premenopausal women, the estrogen levels vary considerably throughout the menstrual cycle . Using daily salivary samples, we were able to measure the free biologically active form of estrogen, which is considered to be the ideal among ovulating premenopausal women [14, 44, 73]. Previous research has indicated that multiple measurements of unbound bioavailable levels probably give us a picture of the real and cumulative estrogen exposure over time [7, 74]. Thus, we were able to capture the continuous estrogen exposure of the included women. We also included standardized repeated serum hormone levels. We have previously observed, based on the same study population, that daily levels of salivary estradiol may be associated with known breast cancer risk factors including mammographic density phenotypes [15, 37]. Moreover, mammographic density measures were obtained within a narrow time frame during the late follicular phase (days 7–12), thereby avoiding the bias of variation in mammographic density throughout the menstrual cycle . The validated computer-assisted method used to quantify mammographic density has been shown to give a superior prediction of breast cancer risk compared to qualitative methods . A single experienced blinded reader read all mammograms, and the assessed mammographic density was negatively associated with age, BMI, and parity .
The relatively small sample size of the current study must be considered as one of the limitations warranting attention, and underlines the need for further and larger studies. The study design was cross-sectional, and population stratification can be of concern in this type of study , but our study population consisted of Norwegian/Caucasian women with only a ten-year age span.
Based on the biological hypothesis that polymorphisms in the ESR1 gene may influence 17β-estradiol levels and mammographic phenotypes, we examined a limited number of SNPs . Even though there is a risk of false-positive results/type I errors in multiple testing, our analysis are based on carefully considered plausible biological mechanisms, and an a priori hypothesis of difference rather than a universal null-hypothesis. Hence, we considered the increased likelihood of type II errors and the risk of deeming truly interesting difference non-significant when introducing Bonferroni-adjustments just as important.
In recent studies, the immunoassay methods used in the present study is most often replaced by Liquid Chromatography (LC) - Mass Spectrometry (MS/MS), which is considered a more efficient method of analysing salivary hormones with higher specificity and sensitivity. However, studies on estradiol measurements have demonstrated a high correlation between LC/MS-MS and immunoassays of 0.969 [15, 17, 78].
The 6q25 region is important in the etiology of breast cancer among women with dense breasts , but its putative functions are still undefined. Identifying genetic variants that are associated with both breast cancer risk and mammographic density measures that predict breast cancer has the potential to reveal underlying biological pathways that explain the associations between those mammographic measures and cancer. Our data provide additional evidence for the relationship between genetic variants, mammographic density phenotypes and estradiol. We hypothesize that ESR1 may play a role in the relationship between dense breast tissue and estradiol levels, and therefore, the risk of breast cancer.
Body mass index
Coefficient of variance
Norwegian Energy Balance and Breast cancer Aspects
Estrogen receptor alpha gene
Genome wide association study
Minor allele frequency
Region of interest
Single nucleotide polymorphism
University hospital of North Norway
Ursin G, Lillie EO, Lee E, Cockburn M, Schork NJ, Cozen W, Parisky YR, Hamilton AS, Astrahan MA, Mack T. The relative importance of genetics and environment on mammographic density. Cancer Epidemiol Biomarkers Prev. 2009;18(1):102–12.
Boyd NF, Guo H, Martin LJ, Sun L, Stone J, Fishell E, Jong RA, Hislop G, Chiarelli A, Minkin S, et al. Mammographic density and the risk and detection of breast cancer. N Engl J Med. 2007;356(3):227–36.
Vachon CM, Kuni CC, Anderson K, Anderson VE, Sellers TA. Association of mammographically defined percent breast density with epidemiologic risk factors for breast cancer (United States). Cancer Causes Control. 2000;11(7):653–62.
Russo J, Russo IH. Genotoxicity of steroidal estrogens. Trends Endocrinol Metab. 2004;15(5):211–4.
Folkerd E, Dowsett M. Sex hormones and breast cancer risk and prognosis. Breast. 2013;22 Suppl 2:S38–43.
Key T, Appleby P, Barnes I, Reeves G, Endogenous H, Breast Cancer Collaborative G. Endogenous sex hormones and breast cancer in postmenopausal women: reanalysis of nine prospective studies. J Natl Cancer Inst. 2002;94(8):606–16.
Endogenous H, Breast Cancer Collaborative G, Key TJ, Appleby PN, Reeves GK, Travis RC, Alberg AJ, Barricarte A, Berrino F, Krogh V, et al. Sex hormones and risk of breast cancer in premenopausal women: a collaborative reanalysis of individual participant data from seven prospective studies. Lancet Oncol. 2013;14(10):1009–19.
Yaghjyan L, Colditz GA, Rosner B, Tamimi RM. Mammographic breast density and subsequent risk of breast cancer in postmenopausal women according to the time since the mammogram. Cancer Epidemiol Biomarkers Prev. 2013;22(6):1110–7.
Boyd NF, Jensen HM, Cooke G, Han HL, Lockwood GA, Miller AB. Mammographic densities and the prevalence and incidence of histological types of benign breast disease. Reference Pathologists of the Canadian National Breast Screening Study. Eur J Cancer Prev. 2000;9(1):15–24.
McCormack VA, dos Santos SI. Breast density and parenchymal patterns as markers of breast cancer risk: a meta-analysis. Cancer Epidemiol Biomarkers Prev. 2006;15(6):1159–69.
Pettersson A, Graff RE, Ursin G, Santos Silva ID, McCormack V, Baglietto L, Vachon C, Bakker MF, Giles GG, Chia KS et al. Mammographic density phenotypes and risk of breast cancer: a meta-analysis. J Natl Cancer Inst. 2014;106(5). doi:10.93/jnci/dju078.Review.
Greendale GA, Reboussin BA, Slone S, Wasilauskas C, Pike MC, Ursin G. Postmenopausal hormone therapy and change in mammographic density. J Natl Cancer Inst. 2003;95(1):30–7.
McTiernan A, Martin CF, Peck JD, Aragaki AK, Chlebowski RT, Pisano ED, Wang CY, Brunner RL, Johnson KC, Manson JE, et al. Estrogen-plus-progestin use and mammographic density in postmenopausal women: Women’s Health Initiative randomized trial. J Natl Cancer Inst. 2005;97(18):1366–76.
Frydenberg H, Flote VG, Iversen A, Finstad SE, Furberg AS, Torjesen PA, Wilsgaard T, Schlichting E, Ellison PT, Ursin G, et al. Insulin-like growth factor-1, growth hormone, and daily cycling estrogen are associated with mammographic density in premenopausal women. Cancer Causes Control. 2014;25(7):891–903.
Frydenberg H, Flote VG, Larsson IM, Barrett ES, Furberg AS, Ursin G, Wilsgaard T, Ellison PT, McTiernan A, Hjartaker A, et al. Alcohol consumption, endogenous estrogen and mammographic density among premenopausal women. Breast Cancer Res. 2015;17:103.
Flote VG, Frydenberg H, Ursin G, Iversen A, Fagerland MW, Ellison PT, Wist EA, Egeland T, Wilsgaard T, McTiernan A, et al. High-density lipoprotein-cholesterol, daily estradiol and progesterone, and mammographic density phenotypes in premenopausal women. Cancer Prev Res (Phila). 2015;8(6):535–44.
Flote VG, Furberg A, McTiernan A, Frydenberg H, Ursin G, Iversen A, Lofteroed T, Ellison PT, Wist EA, Egeland T, et al. Gene variations in oestrogen pathways, CYP19A1, daily 17ss-estradiol and mammographic density phenotypes in premenopausal women. Breast Cancer Res. 2014;16(6):499.
Ursin G, Hovanessian-Larsen L, Parisky YR, Pike MC, Wu AH. Greatly increased occurrence of breast cancers in areas of mammographically dense tissue. Breast Cancer Res. 2005;7(5):R605–8.
Gill JK, Maskarinec G, Pagano I, Kolonel LN. The association of mammographic density with ductal carcinoma in situ of the breast: the Multiethnic Cohort. Breast Cancer Res. 2006;8(3):R30.
Rossouw JE, Anderson GL, Prentice RL, LaCroix AZ, Kooperberg C, Stefanick ML, Jackson RD, Beresford SA, Howard BV, Johnson KC, et al. Risks and benefits of estrogen plus progestin in healthy postmenopausal women: principal results From the Women’s Health Initiative randomized controlled trial. JAMA. 2002;288(3):321–33.
Boyd NF, Melnichouk O, Martin LJ, Hislop G, Chiarelli AM, Yaffe MJ, Minkin S. Mammographic density, response to hormones, and breast cancer risk. J Clin Oncol. 2011;29(22):2985–92.
Ahmed S, Thomas G, Ghoussaini M, Healey CS, Humphreys MK, Platte R, Morrison J, Maranian M, Pooley KA, Luben R, et al. Newly discovered breast cancer susceptibility loci on 3p24 and 17q23.2. Nat Genet. 2009;41(5):585–90.
Fletcher O, Johnson N, Orr N, Hosking FJ, Gibson LJ, Walker K, Zelenika D, Gut I, Heath S, Palles C, et al. Novel breast cancer susceptibility locus at 9q31.2: results of a genome-wide association study. J Natl Cancer Inst. 2011;103(5):425–35.
Yamamoto-Ibusuki M, Yamamoto Y, Fujiwara S, Sueta A, Yamamoto S, Hayashi M, Tomiguchi M, Takeshita T, Iwase H. C6ORF97-ESR1 breast cancer susceptibility locus: influence on progression and survival in breast cancer patients. Eur J Hum Genet. 2014;23(7):949–56.
Crandall CJ, Sehl ME, Crawford SL, Gold EB, Habel LA, Butler LM, Sowers MR, Greendale GA, Sinsheimer JS. Sex steroid metabolism polymorphisms and mammographic density in pre- and early perimenopausal women. Breast Cancer Res. 2009;11(4):R51.
Dumas I, Diorio C. Polymorphisms in genes involved in the estrogen pathway and mammographic density. BMC Cancer. 2010;10:636.
Haiman CA, Hankinson SE, De Vivo I, Guillemette C, Ishibe N, Hunter DJ, Byrne C. Polymorphisms in steroid hormone pathway genes and mammographic density. Breast Cancer Res Treat. 2003;77(1):27–36.
Maskarinec G, Lurie G, Williams AE, Le Marchand L. An investigation of mammographic density and gene variants in healthy women. Int J Cancer. 2004;112(4):683–8.
van Duijnhoven FJ, Bezemer ID, Peeters PH, Roest M, Uitterlinden AG, Grobbee DE, van Gils CH. Polymorphisms in the estrogen receptor alpha gene and mammographic density. Cancer Epidemiol Biomarkers Prev. 2005;14(11 Pt 1):2655–60.
van Duijnhoven FJ, Peeters PH, Warren RM, Bingham SA, Uitterlinden AG, van Noord PA, Monninkhof EM, Grobbee DE, van Gils CH. Influence of estrogen receptor alpha and progesterone receptor polymorphisms on the effects of hormone therapy on mammographic density. Cancer Epidemiol Biomarkers Prev. 2006;15(3):462–7.
Warren R, Skinner J, Sala E, Denton E, Dowsett M, Folkerd E, Healey CS, Dunning A, Doody D, Ponder B, et al. Associations among mammographic density, circulating sex hormones, and polymorphisms in sex hormone metabolism genes in postmenopausal women. Cancer Epidemiol Biomarkers Prev. 2006;15(8):1502–8.
Barzan D, Veldwijk MR, Herskind C, Li Y, Zhang B, Sperk E, Du WD, Zhang XJ, Wenz F. Comparison of genetic variation of breast cancer susceptibility genes in Chinese and German populations. Eur J Hum Genet. 2013;21(11):1286–92.
Haiman CA, Bernstein L, Berg D, Ingles SA, Salane M, Ursin G. Genetic determinants of mammographic density. Breast Cancer Res. 2002;4(3):R5.
Li J, Eriksson L, Humphreys K, Czene K, Liu J, Tamimi RM, Lindstrom S, Hunter DJ, Vachon CM, Couch FJ, et al. Genetic variation in the estrogen metabolic pathway and mammographic density as an intermediate phenotype of breast cancer. Breast Cancer Res. 2010;12(2):R19.
Stone J, Gurrin LC, Byrnes GB, Schroen CJ, Treloar SA, Padilla EJ, Dite GS, Southey MC, Hayes VM, Hopper JL. Mammographic density and candidate gene variants: a twins and sisters study. Cancer Epidemiol Biomarkers Prev. 2007;16(7):1479–84.
Takata Y, Maskarinec G, Le Marchand L. Breast density and polymorphisms in genes coding for CYP1A2 and COMT: the Multiethnic Cohort. BMC Cancer. 2007;7:30.
Iversen A, Frydenberg H, Furberg AS, Flote VG, Finstad SE, McTiernan A, Ursin G, Wilsgaard T, Ellison PT, Jasienska G, et al. Cyclic endogenous estrogen and progesterone vary by mammographic density phenotypes in premenopausal women. Eur J Cancer Prev. 2015;25(1):9–18.
Furberg AS, Jasienska G, Bjurstam N, Torjesen PA, Emaus A, Lipson SF, Ellison PT, Thune I. Metabolic and hormonal profiles: HDL cholesterol as a plausible biomarker of breast cancer risk. The Norwegian EBBA Study. Cancer Epidemiol Biomarkers Prev. 2005;14(1):33–40.
Iversen A, Thune I, McTiernan A, Emaus A, Finstad SE, Flote V, Wilsgaard T, Lipson SF, Ellison PT, Jasienska G, et al. Ovarian hormones and reproductive risk factors for breast cancer in premenopausal women: the Norwegian EBBA-I study. Hum Reprod. 2011;26(6):1519–29.
Emaus A, Veierod MB, Furberg AS, Espetvedt S, Friedenreich C, Ellison PT, Jasienska G, Andersen LB, Thune I. Physical activity, heart rate, metabolic profile, and estradiol in premenopausal women. Med Sci Sports Exerc. 2008;40(6):1022–30.
Finstad SE, Emaus A, Tretli S, Jasienska G, Ellison PT, Furberg AS, Wist EA, Thune I. Adult height, insulin, and 17beta-estradiol in young women. Cancer Epidemiol Biomarkers Prev. 2009;18(5):1477–83.
Iversen A, Thune I, McTiernan A, Makar KW, Wilsgaard T, Ellison PT, Jasienska G, Flote V, Poole EM, Furberg AS. Genetic polymorphism CYP17 rs2486758 and metabolic risk factors predict daily salivary 17beta-estradiol concentration in healthy premenopausal Norwegian women. The EBBA-I study. J Clin Endocrinol Metab. 2012;97(5):E852–7.
Lipson SF, Ellison PT. Comparison of salivary steroid profiles in naturally occurring conception and non-conception cycles. Hum Reprod. 1996;11(10):2090–6.
Ellison PT, Lipson SF. Salivary estradiol--a viable alternative? Fertil Steril. 1999;72(5):951–2.
Passarelli MN, Phipps AI, Potter JD, Makar KW, Coghill AE, Wernli KJ, White E, Chan AT, Hutter CM, Peters U, et al. Common single-nucleotide polymorphisms in the estrogen receptor beta promoter are associated with colorectal cancer survival in postmenopausal women. Cancer Res. 2013;73(2):767–75.
Carlson CS, Eberle MA, Rieder MJ, Yi Q, Kruglyak L, Nickerson DA. Selecting a maximally informative set of single-nucleotide polymorphisms for association analyses using linkage disequilibrium. Am J Hum Genet. 2004;74(1):106–20.
Thorisson GA, Smith AV, Krishnan L, Stein LD. The International HapMap Project Web site. Genome Res. 2005;15(11):1592–3.
Bjurstam N, Bjorneld L, Warwick J, Sala E, Duffy SW, Nystrom L, Walker N, Cahlin E, Eriksson O, Hafstrom LO, et al. The Gothenburg Breast Screening Trial. Cancer. 2003;97(10):2387–96.
Ursin G, Astrahan MA, Salane M, Parisky YR, Pearce JG, Daniels JR, Pike MC, Spicer DV. The detection of changes in mammographic densities. Cancer Epidemiol Biomarkers Prev. 1998;7(1):43–7.
Ursin G, Ma H, Wu AH, Bernstein L, Salane M, Parisky YR, Astrahan M, Siozon CC, Pike MC. Mammographic density and breast cancer in three ethnic groups. Cancer Epidemiol Biomarkers Prev. 2003;12(4):332–8.
van Gils CH, Hendriks JH, Otten JD, Holland R, Verbeek AL. Parity and mammographic breast density in relation to breast cancer risk: indication of interaction. Eur J Cancer Prev. 2000;9(2):105–11.
Robbez-Masson LJ, Bodor C, Jones JL, Hurst HC, Fitzgibbon J, Hart IR, Grose RP. Functional analysis of a breast cancer-associated FGFR2 single nucleotide polymorphism using zinc finger mediated genome editing. PLoS One. 2013;8(11):e78839.
Milne RL, Burwinkel B, Michailidou K, Arias-Perez JI, Zamora MP, Menendez-Rodriguez P, Hardisson D, Mendiola M, Gonzalez-Neira A, Pita G, et al. Common non-synonymous SNPs associated with breast cancer susceptibility: findings from the Breast Cancer Association Consortium. Hum Mol Genet. 2014;23(22):6096–111.
Bai Z, Gust R. Breast cancer, estrogen receptor and ligands. Arch Pharm (Weinheim). 2009;342(3):133–49.
Brandenberger AW, Tee MK, Jaffe RB. Estrogen receptor alpha (ER-alpha) and beta (ER-beta) mRNAs in normal ovary, ovarian serous cystadenocarcinoma and ovarian cancer cell lines: down-regulation of ER-beta in neoplastic tissues. J Clin Endocrinol Metab. 1998;83(3):1025–8.
Russo J, Ao X, Grill C, Russo IH. Pattern of distribution of cells positive for estrogen receptor alpha and progesterone receptor in relation to proliferating cells in the mammary gland. Breast Cancer Res Treat. 1999;53(3):217–27.
Russo J, Russo IH. The role of estrogen in the initiation of breast cancer. J Steroid Biochem Mol Biol. 2006;102(1–5):89–96.
Lundstrom E, Sahlin L, Skoog L, Hagerstrom T, Svane G, Azavedo E, Sandelin K, von Schoultz B. Expression of Syndecan-1 in histologically normal breast tissue from postmenopausal women with breast cancer according to mammographic density. Climacteric. 2006;9(4):277–82.
Borras M, Hardy L, Lempereur F, el Khissiin AH, Legros N, Gol-Winkler R, Leclercq G. Estradiol-induced down-regulation of estrogen receptor. Effect of various modulators of protein synthesis and expression. J Steroid Biochem Mol Biol. 1994;48(4):325–36.
Dunning AM, Michailidou K, Kuchenbaecker KB, Thompson D, French JD, Beesley J, Healey CS, Kar S, Pooley KA, Lopez-Knowles E, et al. Breast cancer risk variants at 6q25 display different phenotype associations and regulate ESR1, RMND1 and CCDC170. Nat Genet. 2016;48(4):374–86.
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.
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.
Onland-Moret NC, van Gils CH, Roest M, Grobbee DE, Peeters PH. The estrogen receptor alpha gene and breast cancer risk (The Netherlands). Cancer Causes Control. 2005;16(10):1195–202.
Michailidou K, Hall P, Gonzalez-Neira A, Ghoussaini M, Dennis J, Milne RL, Schmidt MK, Chang-Claude J, Bojesen SE, Bolla MK, et al. Large-scale genotyping identifies 41 new loci associated with breast cancer risk. Nat Genet. 2013;45(4):353–61. 361e351-352.
Brand JS, Humphreys K, Thompson DJ, Li J, Eriksson M, Hall P, Czene K. Volumetric mammographic density: heritability and association with breast cancer susceptibility loci. J Natl Cancer Inst. 2014;106(12). doi:10.093/jnci/dju334.Print2014dec.
Lindstrom S, Thompson DJ, Paterson AD, Li J, Gierach GL, Scott C, Stone J, Douglas JA, Dos-Santos-Silva I, Fernandez-Navarro P, et al. Genome-wide association study identifies multiple loci associated with both mammographic density and breast cancer risk. Nat Commun. 2014;5:5303.
Lindstrom S, Vachon CM, Li J, Varghese J, Thompson D, Warren R, Brown J, Leyland J, Audley T, Wareham NJ, et al. Common variants in ZNF365 are associated with both mammographic density and breast cancer risk. Nat Genet. 2011;43(3):185–7.
Tsuchiya Y, Nakajima M, Kyo S, Kanaya T, Inoue M, Yokoi T. Human CYP1B1 is regulated by estradiol via estrogen receptor. Cancer Res. 2004;64(9):3119–25.
Tsuchiya Y, Nakajima M, Yokoi T. Cytochrome P450-mediated metabolism of estrogens and its regulation in human. Cancer Lett. 2005;227(2):115–24.
Jiang H, Xie T, Ramsden DB, Ho SL. Human catechol-O-methyltransferase down-regulation by estradiol. Neuropharmacology. 2003;45(7):1011–8.
Dumas I, Diorio C. Estrogen pathway polymorphisms and mammographic density. Anticancer Res. 2011;31(12):4369–86.
Gann PH, Giovanazzi S, Van Horn L, Branning A, Chatterton Jr RT. Saliva as a medium for investigating intra- and interindividual differences in sex hormone levels in premenopausal women. Cancer Epidemiol Biomarkers Prev. 2001;10(1):59–64.
Bellem A, Meiyappan S, Romans S, Einstein G. Measuring estrogens and progestagens in humans: an overview of methods. Gend Med. 2011;8(5):283–99.
Schoemaker MJ, Folkerd EJ, Jones ME, Rae M, Allen S, Ashworth A, Dowsett M, Swerdlow AJ. Combined effects of endogenous sex hormone levels and mammographic density on postmenopausal breast cancer risk: results from the Breakthrough Generations Study. Br J Cancer. 2014;110(7):1898–907.
Morrow M, Chatterton Jr RT, Rademaker AW, Hou N, Jordan VC, Hendrick RE, Khan SA. A prospective study of variability in mammographic density during the menstrual cycle. Breast Cancer Res Treat. 2010;121(3):565–74.
Byrne C, Schairer C, Wolfe J, Parekh N, Salane M, Brinton LA, Hoover R, Haile R. Mammographic features and breast cancer risk: effects with time, age, and menopause status. J Natl Cancer Inst. 1995;87(21):1622–9.
Rothman KJ, Greenland S. Modern epidemiology. 2nd ed. Philadelphia, PA: Lippincott-Raven; 1998.
Holst JP, Soldin OP, Guo T, Soldin SJ. Steroid hormones: relevance and measurement in the clinical laboratory. Clin Lab Med. 2004;24(1):105–18.
Holst F, Stahl PR, Ruiz C, Hellwinkel O, Jehan Z, Wendland M, Lebeau A, Terracciano L, Al-Kuraya K, Janicke F, et al. Estrogen receptor alpha (ESR1) gene amplification is frequent in breast cancer. Nat Genet. 2007;39(5):655–60.
We give special thanks to each woman who participated in the EBBA-I study and acknowledge the efforts and professional assistance of Gunn Kristin Knudsen, Heidi Jakobsen, Anna-Kirsti Kvitnes, and Sissel Andersen from the Clinical Research Department at the University Hospital of North Norway, and Karen Makar for her invaluable contribution in the analysis of the single-nucleotide polymorphism analysis.
This work was supported by the University Hospital of Oslo; Foundation for the Norwegian Health and Rehabilitation Organization (grant number 59010–2000, 59010–2001, 59010–2002); and Norwegian Cancer Society (grant number 05087, TP 49 258). The funders did not influence the conduct of the study, the preparation of this paper, or the decision to publish.
Availability of data and materials
Because of the small sample size and the risk of patient set identification the data set supporting the conclusion of this article will not be presented within the manuscript. However, information can be made available at request from the corresponding author.
IT, GJ and PTE conceived of and designed the study and collected the clinical data. PTE carried out the radioimmunoassay and measured salivary hormone levels. GU digitized and measured the mammographic density data. AM performed the single-nucleotide polymorphism analysis. FNF, VGF, HF, TW, IT and AM performed statistical analyses. FNF, VGF, HF, ASF, AM, GU, PTE, EAW, EBS, TW, GJ, and IT interpreted the results. FNF drafted the manuscript in cooperation with IT. All authors contributed to critical revision of the manuscript, editing of the final version of the manuscript and agree to be accountable for the accuracy and integrity of the work. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethical approval and consent to participate
All participants underwent informed consent procedures and signed a consent form. The Norwegian Data Inspectorate and the Regional Committee for Medical Research Ethics approved the study.
About this article
Cite this article
Fjeldheim, F.N., Frydenberg, H., Flote, V.G. et al. Polymorphisms in the estrogen receptor alpha gene (ESR1), daily cycling estrogen and mammographic density phenotypes.. BMC Cancer 16, 776 (2016). https://doi.org/10.1186/s12885-016-2804-1
- Mammographic density