<figure><img src="image_1.png" /></figure>

HAL
open science

# Glyceollins trigger anti-proliferative effects through estradiol-dependent and independent pathways in breast cancer cells

Sylvain Lecomte, Frédéric Chalmel, François Ferrière, Frédéric Percevault,
Nicolas Plu, Christian Saligaut, Claire Surel, Marie Lelong, Théo Efstathiou,
Farzad Pakdel

<figure><img src="image_2.png" /></figure>

## To cite this version:

Sylvain Lecomte, Frédéric Chalmel, François Ferrière, Frédéric Percevault, Nicolas Plu, et al.. Glyceollins trigger anti-proliferative effects through estradiol-dependent and independent pathways in breast cancer cells. Cell Communication and Signaling, 2017, 15 (7), pp.26. 10.1186/s12964-017-0182-1. hal-01560288

HAL Id: hal-01560288

https://univ-rennes.hal.science/hal-01560288

Submitted on 21 Jun 2018

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_3.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH

Open Access

<figure><img src="image_4.png" /></figure>

CrossMark

# Glyceollins trigger anti-proliferative effects through estradiol-dependent and independent pathways in breast cancer cells

Sylvain Lecomte<fnref n="1" />,<fnref n="2" />, Frederic Chalmel<fnref n="1" />,³, François Ferriere<fnref n="1" />,<fnref n="2" />, Frederic Percevault<fnref n="1" />,<fnref n="2" />, Nicolas Plu⁴, Christian Saligaut<fnref n="1" />,<fnref n="2" />, Claire Surel⁴, Marie Lelong<fnref n="1" />,<fnref n="2" />, Theo Efstathiou⁴ and Farzad Pakdel<fnref n="1" />,<fnref n="2" />*

**Abstract**

**Background:** Estrogen receptors (ER) α and β are found in both women and men in many tissues, where they have different functions, including having roles in cell proliferation and differentiation of the reproductive tract. In addition to estradiol (E2), a natural hormone, numerous compounds are able to bind ERs and modulate their activities. Among these compounds, phytoestrogens such as isoflavones, which are found in plants, are promising therapeutics for several pathologies. Glyceollins are second metabolites of isoflavones that are mainly produced in soybean in response to an elicitor. They have potentially therapeutic actions in breast cancer by reducing the proliferation of cancer cells. However, the molecular mechanisms driving these effects remain elusive.

**Methods:** First, to determine the proliferative or anti-proliferative effects of glyceollins, in vivo and in vitro approaches were used. The length of epithelial duct in mammary gland as well as uterotrophy after treatment by E2 and glyceollins and their effect on proliferation of different breast cell line were assessed. Secondly, the ability of glyceollin to activate ER was assessed by luciferase assay. Finally, to unravel molecular mechanisms involved by glyceollins, transcriptomic analysis was performed on MCF-7 breast cancer cells.

**Results:** In this study, we show that synthetic versions of glyceollin I and II exert anti-proliferative effects in vivo in mouse mammary glands and in vitro in different ER-positive and ER-negative breast cell lines. Using transcriptomic analysis, we produce for the first time an integrated view of gene regulation in response to glyceollins and reveal that these phytochemicals act through at least two major pathways. One pathway involving FOXM1 and ERα is directly linked to proliferation. The other involves the HIF family and reveals that stress is a potential factor in the anti-proliferative effects of glyceollins due to its role in increasing the expression of REDD1, an mTORC1 inhibitor.

**Conclusion:** Overall, our study clearly shows that glyceollins exert anti-proliferative effects by reducing the expression of genes encoding cell cycle and mitosis-associated factors and biomarkers overexpressed in cancers and by increasing the expression of growth arrest-related genes. These results reinforce the therapeutic potential of glyceollins for breast cancer.

**Keywords:** Glyceollins, Estrogen receptors, Cell proliferation, Transcriptomic, Gene expression, Breast cancer

<footnote>* Correspondence: farzad.pakdel@univ-rennes1.fr</footnote>

<footnote n="1">¹Institut de Recherche en Santé-Environnement-Travail (IRSET), University of Rennes 1, 9 Avenue du Pr Léon Bernard, 35000 Rennes, France</footnote>

<footnote n="2">²Inserm U1085, Team Transcription, Environment and Cancer, 9 Avenue du Pr Léon Bernard, 35000 Rennes, France</footnote>

<footnote>Full list of author information is available at the end of the article</footnote>
## Background

Breast cancer is the most prevalent cancer in women worldwide and has an incidence of 89.7 per 100,000 women in Western Europe. The WHO estimates that breast cancer was responsible for over 500,000 deaths in 2011 [1]. Among the different types of breast cancer, the most common is estrogen receptor (ER)-positive cancer, which represents approximately 80% of diagnosed cases of cancer. The ER belongs to the nuclear receptor super-family and is divided into two subtypes, ERα and ERβ. The ER acts in cells by directly binding to DNA on a responsive element called the estrogen-responsive element (ERE) or by interacting with other transcription factors, such as stimulating protein 1 (Sp1) or activator protein 1 (AP1), which are already bound to responsive elements in promoter regions. The ER also modulates signaling pathways such as the MAPK and PI3K/AKT pathways [2]. Thus, ERα, the major isoform in breast tissue, plays an essential role in normal mammary gland development and function as well as in breast cancer initiation and growth. ERs are bound by the natural hormone estradiol (E2), which has a pleiotropic effect and is responsible for the proliferation and survival of breast epithelial cells. Therefore, E2 plays an important role in breast cancer growth. However, E2 is also essential for maintaining cell differentiation, which consequently limits metastatic potential. Hence, ERα is a good prognostic marker and a prime target for therapy. Endocrine therapies such as tamoxifen or fulvestrant are effective for ER-positive cancer, but frequent relapses are observed [3]. Currently, the focus is on the discovery of new compounds with selective ER-modulator (SERM) activities. Phytoestrogens appear to be promising candidates and have been well studied [4]. Phytoestrogens could be ERα agonist or antagonist depending of cellular type and phenotype studied such as cell proliferation and differentiation [5]. Among this family, pterocarpans, which are second metabolites of isoflavones and the best-known members are glyceollins, have been studied since the 2000s.

Glyceollins are phytoalexins produced mainly by soybeans after elicitation by different types of stressors, such as UV, low temperatures or microorganisms. The glyceollin family includes three compounds: glyceollin I (GI), glyceollin II (GII) and glyceollin III (GIII). In plants, glyceollins are involved in host defense against pathogens such as fungi [6, 7] or nematodes [8]. Glyceollins are promising therapeutic compounds for numerous human pathologies, including breast cancer [9]. Interactions between glyceollins and ERα or ERβ were described for the first time in 2000 [10]. Glyceollins act as antiestrogenic compounds that directly interact with both ER isoforms [11, 12], and they have the capacity to suppress tumorigenesis of breast and ovarian cancer [13]. Glyceollin I is the most potent antiestrogenic molecule, and docking

experiments have shown that glyceollin I can interact with the ER in a similar manner to tamoxifen to exert antagonist activity [12]. Thus, chemically synthesized glyceollin I was generated and assessed. It was shown that the natural enantiomer exerted anti-proliferative activities against numerous cell lines, including ER-positive breast cancer cells [14], and that the compound is a potent inhibitor of ER activation [15].

However, the precise antiestrogenic mechanisms associated with glyceollins in ER-positive breast cancer remain elusive. In this work, we synthesized natural enantiomers of glyceollin I and II to determine their impact in vivo on the growth of galactophore ducts in mouse mammary glands as well as in vitro in different breast cell lines. The ability of glyceollins to bind and activate ERs as well as their effects on the expression of endogenous E2-dependent genes were assessed. Glyceollins showed surprising effects on gene expression, which led us to perform transcriptomic analysis of the ER-positive breast cell line MCF-7 to better elucidate the mechanisms underlying the actions of these compounds. We found that glyceollins exert their effects through both ER-dependent and ER-independent pathways involving different transcription factors.

## Methods

### Animals

Ovariectomized (ovx) (at 4.5 weeks) and intact non-ovariectomized (non-ovx) RjOrl SWISS female mice were purchased from Janvier Labs. The animals were housed at 20 ± 2 °C under 60 ± 10% humidity with a 12 h light/dark cycle. After 3 weeks of acclimation, 9-week-old mice were weighed and randomly assigned to one of 11 groups (Table 1) of at least six mice each with similar average body weights. Treatment (Table 1) was delivered by subcutaneous injection (sc) at 24 h intervals for 3 consecutive days with sesame oil as a vehicle. On

<table><caption>**Table 1** Description of randomized mouse groups and their treatments</caption><thead><tr><td>Mice groups</td><td>Treatment</td></tr></thead><tbody><tr><td>Non-Ovariectomized (non-Ovx)</td><td>Vehicle</td></tr><tr><td>Ovariectomized (Ovx)</td><td>Vehicle</td></tr><tr><td>Estradiol (E2)</td><td>E2 10 µg/kg</td></tr><tr><td>Glyceollin I Low (GI L)</td><td>GI 50 mg/kg</td></tr><tr><td>Glyceollin II Low (GII L)</td><td>GII 50 mg/kg</td></tr><tr><td>Estradiol + Glyceollin I Low (E2 + GI L)</td><td>E2 10 µg/kg + GI 50 mg/kg</td></tr><tr><td>Estradiol + Glyceollin II Low (E2 + GII L)</td><td>E2 10 µg/kg + GII 50 mg/kg</td></tr><tr><td>Glyceollin I High (GI H)</td><td>GI 100 mg/kg</td></tr><tr><td>Glyceollin II High (GI H)</td><td>GII 100 mg/kg</td></tr><tr><td>Estradiol + Glyceollin I High (E2 + GI H)</td><td>E2 10 µg/kg + GI 100 mg/kg</td></tr><tr><td>Estradiol + Glyceollin II High (E2 + GII H)</td><td>E2 10 µg/kg + GII 100 mg/kg</td></tr></tbody></table>
day 3, the animals were weighed and then euthanized by cervical dislocation. Uteri and mammary glands were removed.

## Uterotrophic measurements

The uterus was carefully dissected at the level of the vaginal fornix, trimmed of fascia and fat, gently blotted on moistened filter paper and weighed.

## Mammary whole-mount preparations and immunostaining

The first inguinal mammary fat pads were removed and stained as described by Tian et al. [16]. Briefly, mammary fat pads were spread as flat as possible on a glass surface and fixed with 4% paraformaldehyde. For assessment of epithelial duct length, the mammary glands were stained overnight in carmine alum (0.2% carmine, 0.5% aluminum potassium sulfate), dehydrated in an ethanol gradient and clarified overnight in xylene. Tissues were then photographed under a SteREO Discovery V8 microscope (Zeiss, original magnification, ×1). For Ki-67 and Epcam immunostaining, mammary fat pads were embedded in Tissue Tek mounting medium (Sakura) and sliced with a cryostat. The slices were then incubated at room temperature with a rabbit anti-Ki-67 antibody (Abcam) and a rat anti-Epcam (sc53532, Santa Cruz) for 1 h in PBS supplemented with 0.3% Triton ×100 and 0.5% milk. A dye-conjugated secondary antibody was then incubated with the sections at room temperature for 1 h in PBS supplemented with 0.3% Triton ×100 and 0.5% milk. Images were obtained with an Imager.Z1 ApoTome AxioCam (Zeiss) microscope and processed with Axio Vision Software. The percentage of Ki-67-positive cells was determined by counting the total numbers of ductal epithelial cells and Ki-67-positive cells using ImageJ software.

## Cell culture and reagents

MCF-7 cells were maintained in DMEM, 4.5 g/L glucose supplemented with non-essential amino acids (NEAA) (Invitrogen) and 10% fetal bovine serum (FBS) (Biowest). T47D cells were maintained in RPMI 1640 supplemented with NEAA, sodium pyruvate (Invitrogen) and 10% FBS (Biowest). HC-11 cells were maintained in RPMI 1640 supplemented with 2 mM L-glutamine, 5 µg/mL insulin (Invitrogen), 0.01 µg/mL epidermal growth factor (EGF) (Abcys) and 10% FBS (Biowest). MCF10-A cells were maintained in DMEM/F12 supplemented with 0.5 µg/mL hydrocortisone, 10 µg/mL insulin (Invitrogen), 20 ng/mL EGF (Abcys), 100 ng/mL cholera toxin (Sigma) and 5% horse serum (Invitrogen). All cell lines were cultured with penicillin/streptomycin (Invitrogen) at 37 °C under 5% CO₂. For steroid treatments, cells were cultured for at least 24 h in steroids

and serum-free DMEM without phenol red and with 2.5% or 5% charcoal/dextran-stripping FBS (Biowest) for MCF-7 and T47D cells, respectively. E2 was purchased from Sigma. Natural enantiomers (6aS and 11aS asymmetric carbon configurations) of glyceollin I and glyceollin II were chemically synthesized by HPC Pharma adapting the synthesis method described by Khupse et al. [17] and Luniwal et al. [18]. The purity was determined at 98% and 99% for glyceollin I and glyceollin II, respectively.

## Proliferation assay

Cells (7500 cells/well for HC-11, 20,000 cells/well for MCF-7 and MCF10-A, and 40,000 cells/well for T47D) were plated in 24-well plates and then deprived of steroids and serum for 72 h. The cells were treated with different doses of glyceollin I or II with or without 10⁻⁹ M E2 for 6 days with renewal of the treatment mixture on day 3. After treatment, the cells were trypsinized, and the cell number was determined using a TC10 Automated Cell Counter (Bio-Rad).

## Luciferase assay

MCF-7 cells (30,000 cells/well) were plated in 24-well plates. After serum and steroid deprivation, the cells were transfected overnight with 100 ng of an ERE-TK-luciferase vector, which encodes luciferase under the control of one ERE, and with 20 ng of a CMV-β galactosidase vector, which served as a control of transfection efficiency control. JetPEI was used as a transfection reagent (Polyplus transfection). Next, the cells were treated with 10⁻⁹ M E2 and/or with different doses of glyceollin I or II. ICI₁₈₂.₇₈₀ (Tocris) was used as ER-inhibitor. The cells were lysed in Passive Lysis Buffer (Promega), and luciferase activity was determined using a commercial luciferase assay system (Promega).

## RNA extraction and real-time PCR

MCF-7 cells (250,000 cells/well) were plated in 6-well plates. After 30 h of serum and steroid deprivation, the cells were treated with solvent as a control, with 10⁻⁹ M E2 or with different concentrations of glyceollin I or II. Total RNA was extracted using a RNeasy mini kit (Qiagen) according to the manufacturer's instructions. Then, the RNA was reverse-transcribed using an M-MLV RT kit (Invitrogen) according to the manufacturer's instructions. For real-time PCR, 5 ng of cDNA was used with 150 nM primers (Table 2) and iTaq Universal SYBR Green Supermix (BioRad). Real-time PCR was performed on a CFX 384 apparatus, and the results were analyzed with CFX Manager software (BioRad).
<table><caption>**Table 2** Gene names and primer sequences used in real-time PCR experiments</caption><thead><tr><th>Gene name and symbol</th><th>Forward primer</th><th>Reverse primer</th></tr></thead><tbody><tr><td>Progesterone receptor (PgR)</td><td>CCCGCCGTCGTAACTTTGG</td><td>GTGCCTATCCTGCCTCTCAATC</td></tr><tr><td>Growth regulation in breast cancer 1 (GREB1)</td><td>GAGGATGTGGAGTGGAGACC</td><td>CAGTACCTCAAAGACCTCGGC</td></tr><tr><td>Trefoil Factor 1(TFFI/pS2)</td><td>ACCATGGAGAACAAGGTGA</td><td>CCGAGCTCTGGGACTAATCA</td></tr><tr><td>Amphiregulin (AREG)</td><td>GTATTTCACTTTCCGTCTTGTTTTG</td><td>CCTGGCTATATTGTCGATTCA</td></tr><tr><td>Forkhead box M1 (FOXM1)</td><td>AGCGAGACCCATCAAAGTGG</td><td>GGTCTTGGGGTGGGAGATTG</td></tr><tr><td>Estrogen receptor 1 (ESR1/ERα)</td><td>TTTATGGGAAAAGGCTCAAA</td><td>GACAAAACCGAGTCACATCA</td></tr><tr><td>Estrogen receptor 2 (ESR2/ERβ)</td><td>AGAGTCCCTGGTGTGAAGCAAG</td><td>GACAGCGCAGAAGTGAGCATC</td></tr><tr><td>FBJ murine osteosarcoma viral oncogene homolog (FOS)</td><td>GAATTAACCTGGTGCTGGAT</td><td>GAACACTATTGCCAGGAA</td></tr><tr><td>Peroxisome proliferator-activated receptor gamma (PPARG)</td><td>GCAATCAAAGTGGAGCCTGC</td><td>CCCTTGCATCCTTCACAAGC</td></tr><tr><td>Hypoxia inducible factor 1, alpha subunit (HIF1A)</td><td>CTGCCACCACTGATGAATTA</td><td>GTATGTGGGTAGGAGATGGA</td></tr><tr><td>Endothelial PAS domain protein 1 (EPAS1/HIF2α)</td><td>GCGCTAGACTCCGAGAACAT</td><td>TGGCCACTTACTACCTGACCCTT</td></tr><tr><td>Vascular endothelial growth factor A (VEGFA)</td><td>AGGAGGAGGGCAGAATCATCA</td><td>CTCGATTGGATGGCAGTAGCT</td></tr><tr><td>DNA-damage-inducible transcript 4 (DDIT4/REDD1)</td><td>AGGAAGCTCATTGAGTTGTG</td><td>GGTACATGCTACACACACAT</td></tr><tr><td>Nuclear Receptor Subfamily 2 Group F member 1 (NR2F1/COUP-TFI)</td><td>TACGTGAGGAGCCAGTACCC</td><td>CGATGGGGGTTTTACCTACC</td></tr><tr><td>Chemokine (C-X-C motif) receptor 4 (CXCR4)</td><td>GCCTTATCCTGCCTGGTATTGTC</td><td>GCGAAGAAAGCCAGGATGAGGA</td></tr><tr><td>Atypical chemokine receptor 3 (ACKR3/CXCR7)</td><td>ACAGGCTATGACACGCACTG</td><td>ACGAGACTGACCACCCAGAC</td></tr><tr><td>Chemokine (C-X-C motif) ligand 12 (CXCL12)</td><td>CTCCTGGGGATGTGTAATGG</td><td>GCCTCCATGGCATACATAGG</td></tr><tr><td>Glyceraldehyde-3-phosphate dehydrogenase (GAPDH)</td><td>GGGCATCCTGGGCTACACTG</td><td>GGGCATCCTGGGCTACACTG</td></tr><tr><td>TATA box binding protein (TBP)</td><td>TGCACAGGAGCCAAGAGTGAA</td><td>CACATCACAGCTCCCCACCA</td></tr></tbody></table>

### Transcriptomic analysis

MCF-7 cells (250,000 cells/well) were plated in 6-well plates. After 30 h of serum and steroid deprivation, the cells were treated for 24 h with solvent as a control, with $10^{-9}$ M E2, with $10^{-5}$ M glyceollin I or II, or with an E2 and glyceollin I or II co-treatment. Total RNA was prepared as described above. RNA quantity and purity were determined using a Nanodrop (Thermo Fisher). Only RNA samples with 260/280 and 260/230 ratios >1.8 were selected. RNA quality was analyzed with a bioanalyzer (Agilent), and RNA samples with a RIN > 8.5 and an 18S/28S ratio > 1.7 were selected for spotting on a SurePrint G3 Human Gene Expression v2 8x60K Microarray (Agilent Technologies). Total RNA was reverse-transcribed and labeled according to the manufacturer's instructions. All samples were prepared and spotted in quadruplicate. Sample hybridization, microarray scanning and results extraction were performed by the GeT-Biopuces Platform in Toulouse, France.

### Microarray data analysis and gene filtration

Data analysis was performed using the AMEN suite of tools [19]. Briefly, probes showing a signal higher than a given background cutoff (median of the normalized dataset, cutoff 5.48) and at least a 2-fold change in at least one pairwise comparison were selected. To define a set of 1852 transcripts displaying significant statistical changes across comparisons, the LIMMA (linear models

for microarray data) package was used (F-value adjusted with the false discovery rate method, $p \le 0.05$) [20]. The resulting probes were then partitioned into eight expression patterns (termed P1-P8) using the k-means algorithm.

### Functional data mining

The enrichment analysis module implemented in AMEN [19] was employed to identify human diseases, biological processes, molecular pathways and subcellular components significantly over-represented in each expression pattern by calculating Fisher's exact probability using the Gaussian hypergeometric function (FDR-adjusted $p$-value $\le 0.01$, number of probes in a given group associated with a given annotation term $\ge 5$).

### Regulatory network analysis

Protein-gene regulation data were downloaded from the Transcription Factor Encyclopedia database [21]. A network representation showing all known protein-gene interactions between transcripts differentially expressed in the current project was drawn using AMEN software.

### Statistical analysis

Mann-Whitney tests were performed using the Bios-taTGV website (http://marne.u707.jussieu.fr/biostatgv/), and significant $p$-values were adjusted with Bonferroni correction.
# Results

## Glyceollins have anti-proliferative properties

To test the in vivo effects of synthetic glyceollins I and II on mammary gland growth and uterotrophy, we first compared the epithelial duct length (Fig. 1a) and the uterine weight (Additional file 1: Figure S1a) from controls, E2-exposed and glyceollin-exposed mice. We used ovx animals at 4.5 weeks in age, which allowed us to work with animals possessing estrogen-sensitive tissues but having never been under estrogenic influence since birth. This ensured reduction of the endogenous hormonal background, low baseline uterine weights and a maximum range of response to administered estrogens. This also allowed efficient monitoring of antiestrogenic effects after exposure to estrogens [22]. As controls, ovx mice were injected daily for three days with vehicle or 10 µg/kg E2. For the treatment groups, ovx mice were injected with 50 mg/kg glyceollin I (GI L) or II (GII L) or with 100 mg/kg GI H or GII H either alone or in combination with E2. Epithelial duct length was measured around the lymphatic ganglion on a surface representing nearly 25% of the total surface of the mammary gland (Fig. 1a cartridge). In the ovx mice, the length of the epithelial duct was markedly reduced (median, 1.17 µm/µm²) compared to non-ovx mice (median, 2.65 µm/µm²). E2 treatment partially restored the length and the ramification of the epithelial duct compared to the non-ovx animals and significantly increased them compared to the ovx mice (median, 1.75 µm/µm², *p* < 0.001). Although treatment with glyceollin I or II alone seemed to enhance epithelial duct proliferation at a low dose (GI L and GII L) (median, 1.35 µm/µm² and 1.52 µm/µm², respectively) no significant differences were found between the ovx vehicle- and E2-treated groups. In fact, epithelial duct length in the glyceollin-treated mice varied widely, ranging from 0.805 µm/µm² to 2.07 µm/µm². At a high dose (GI H), glyceollin I seemed to enhance the growth of the epithelial duct (median, 1.52 µm/µm²). Unlike glyceollin I, high-dose glyceollin II (GII H) did not exert any estrogenic activity on the growth of galactophore ducts (median, 1.04 µm/µm²), which was significantly less than that in the E2-treated mice (*p* < 0.01). In contrast, glyceollins exhibited antiestrogenic activity when combined with E2. As shown in Fig. 1a, co-treatment with E2 and glyceollin I or II significantly reduced epithelial duct length (median, 1.29 µm/µm² and 1.28 µm/µm², respectively, *p* < 0.05 at the low dose and median, 1.18 µm/µm² and 1.29 µm/µm², respectively, *p* < 0.001, at the high dose) compared to the E2-treated mice. In contrast, neither glyceollin showed significant estrogenic or antiestrogenic effects on uterine weight (Additional file 1: Figure S1a). To verify that the differences observed in epithelial duct length between the E2 treatment and the combination treatment of E2 with the two

doses of glyceollin I and II were linked to proliferation, Ki-67 immunostaining was performed (Fig. 1b). Three animals were analyzed from the ovx group and for each treatment group. The immunostaining was performed on slices taken from around the lymphatic ganglion, and the total number of epithelial duct cells was counted. Ki-67 immunostaining colocalized with Epcam immunostaining suggesting that Ki-67 expression is restrained to mammary epithelial cells (Fig. 1b cartridge). In ovx mice, the percentage of Ki-67-positive cells was very low at 4%. As described above, E2 treatment, either with or without concomitant glyceollin treatment, significantly increased (*p* < 0.001) the percentage of Ki-67-positive cells compared to the ovx group. However, at the two doses tested, the glyceollins significantly reduced the number of Ki-67-positive cells compared to the E2 treatment alone.

We next tested the in vitro effects of glyceollin I and II on the growth of ER-positive breast cell lines (MCF-7, T47D and HC-11) and an ER-negative nonmalignant breast cell line (MCF10-A). The cells were treated with vehicle and different doses of glyceollin I and II either alone or in combination with E2 for six days. As the proliferation of ER-positive cells is controlled by estrogens, the MCF-7 and T47D cell lines were also treated with 10⁻⁹ M E2, which respectively led to 5- and 3.5-fold increases in cell number (Fig. 1c and d). Interestingly, when cells were treated with E2 in combination with glyceollin I or II, a significant anti-proliferative effect was observed for a doses of 10⁻⁶ M glyceollin I and 5 × 10⁻⁶ M glyceollin II compared to the E2-treated cells. Thus, glyceollin I showed stronger anti-proliferative effects than glyceollin II. However, globally, glyceollin treatment alone did not augment cell number, except for glyceollin I at a dose of 10⁻⁶ M, which showed a very low proliferative effect in MCF-7 cells. It should also be noted that both glyceollins showed a weak but significant anti-proliferative effect at a dose of 10⁻⁵ M compared to the vehicle-treated control cells. E2 did not have any effect on the proliferation of the mouse mammary gland cell line HC-11 (data not shown). However, the glyceollins showed a strong anti-proliferative effect that was statistically significant starting from 5 × 10⁻⁶ M (Fig. 1e). Unlike the other ER-positive cell lines used above, glyceollin I and II showed the same dose effect. Because HC-11 cells require EGF and insulin for growth, it is possible that the glyceollins affected the signaling pathways involved by these two growth factors. To assess the precise role of the ER in this phenotype, ER-negative MCF10-A breast cancer cells were treated with different doses of glyceollin I and II for six days (Fig. 1f). A significant decrease in proliferation was observed with both glyceollins at 10⁻⁵ M compared to vehicle-treated control cells. However, the decrease was less evident in
<figure><img src="image_5.png" /><figcaption>b</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>c</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>d</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>e</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>f</figcaption></figure>

<figure><img src="image_10.png" /></figure>

**Fig. 1** Effects of glyceollin I and II on epithelial duct growth in mouse mammary glands and on ER-positive and ER-negative breast cell lines. Mammary glands were obtained from non-ovariectomized (non-ovx) or ovariectomized (ovx) mice treated with vehicle, 10 µg/kg E2, 50 mg/kg glyceollin I or glyceollin II (GI L or GII L, respectively), or 100 mg/kg glyceollin I or glyceollin II (GI H or GII H, respectively) either alone or in combination with E2 for 72 h (a). Epithelial duct length was assessed by delimiting an area around the lymphatic ganglion and measuring the length of this area using FUJI software (Cartridge). The results are represented in box-and-whisker plots, where the top and the bottom of the box correspond to the 75th and the 25th percentile, respectively. The horizontal bar in the box is the median, and the points outside the box correspond to extreme values. The results are expressed in µm/µm² +/- SEM and were obtained from four independent experiments with at least 5 mice per group. \*p-value <0.05 and \*\*p-value <0.001 with a Mann-Whitney test followed by Bonferroni correction for comparisons of ovx vs treated mice. \#p-value <0.05, \#\#p-value <0.01 and \#\#\#p-value <0.001 with a Mann-Whitney test followed by Bonferroni correction for comparisons of E2 vs the other treatments. The anti-proliferative effect of glyceollin was assessed using Ki-67 and Epcam immunostaining of mammary gland tissue (b and cartridge). The immunostaining was performed on frozen slices of mammary glands around the lymphatic ganglion for ovx and E2-treated mice in combination with the two glyceollin concentrations. Then, the number of Ki-67-positive epithelial cells was determined with FUJI software. The results are expressed as the percentage of Ki-67-positive epithelial cells and represent the mean +/- SEM of three independent mice per treatment. \*\*\#p-value <0.001 with a Mann-Whitney test followed by Bonferroni correction for comparisons of ovx vs treated mice. \#\#\#p-value <0.001 with a Mann-Whitney test followed by Bonferroni correction for comparisons of E2 vs the other treatments. The effects of glyceollin on cell proliferation were also assessed in the ER-positive breast cell lines MCF-7 (c), T47D (d), and HC-11 (e) and the ER-negative breast cell line MCF10A (f). Cells were grown for 6 days with or without 10⁻⁹ M estradiol (E2) in the presence or absence of different concentrations of glyceollins. The numbers of cells measured in the presence of glyceollin I (black hatched line) or II (grey hatched line), either alone or in combination with E2 (full line), are expressed as percentages of E2-treated cells (c and d) or untreated control cells (e and f). The results are presented as the mean +/- SEM of four independent experiments
these cells compared to the MCF-7 or T47D cell lines. Altogether, these data show that the anti-proliferative effects of glyceollins are primarily produced by the ER but also occur through other ER-independent pathways. To establish whether the glyceollins are cytostatic or cytotoxic, cell cycle and apoptosis were assessed in MCF-7 cells after treatment with 10⁻⁵ M glyceollins I or II either alone or in combination with E2 (Additional files 1: Figure S1b and c, and 2 respectively). Cell cycle was analyzed by flow cytometry and showed that E2 induced MCF-7 cells to enter the cell cycle by significantly increasing the percentage of cells in S and G2/M phases compared to control cells. In contrast, glyceollins I and II only weakly induced cells to enter S phase and blocked their passage to G2/M phase. In combination with E2, glyceollins I and II reduced cell entry into S phase compared to treatment with E2 alone, but the combination clearly blocked cells from entering G2/M phase. Apoptosis was analyzed by TUNEL assay. As described in numerous previous studies, E2 reduced the percentage of apoptotic cells. In contrast, 10⁻⁵ M glyceollin, whether alone or in combination with E2, did not significantly increase the percentage of apoptotic cells. Thus, glyceollins appear to be more cytostatic than cytotoxic.

## Glyceollins interact with the ER

To verify ER activation by glyceollins, MCF-7 cells were transfected with an ERE-TK-luciferase reporter gene containing an ERE that has been classically used to assay the estrogenic potencies of xenoestrogens (Fig. 2). The cells were treated with 10⁻⁹ M E2, which served as a positive control, and different concentrations of glyceollin I or II with or without 10⁻⁶ M ICI₁₈₂,₇₈₀ (Fig. 2a and b). As shown in Fig. 2a, both glyceollins weakly activated the ER, as a significant increase in luciferase activity was observed in cells treated with 10⁻⁶ M and 10⁻⁵ M glyceollin I or II compared to cells treated with solvent (*p* < 0.05). Nevertheless, even if these glyceollin concentrations induced expression of the luciferase reporter gene, they did not reach the maximal transactivation efficiency observed with 10⁻⁹ M E2. Indeed, glyceollin I and II showed approximately 7.5- and 4-fold less induction, respectively, than that obtained with E2. This activation was clearly linked to ER as shown in Fig. 2b where a co-treatment with ER-inhibitor, ICI₁₈₂,₇₈₀, abolished the increase of luciferase activity. Surprisingly, when MCF-7 cells were treated with 10⁻⁹ M E2 and different doses of glyceollin I or II (Fig. 2c), neither glyceollin I nor glyceollin II had an antiestrogenic effect on E2

<figure><img src="image_11.png" /></figure>

<figure><img src="image_12.png" /></figure>

<figure><img src="image_13.png" /></figure>

<figure><img src="image_14.png" /><figcaption>**Fig. 2** Effect of glyceollins I and II on ER activation by luciferase assay. MCF-7 cells were transfected with an ERE-TK-luciferase reporter plasmid and a CMV-β-galactosidase plasmid as a control of transfection efficiency. The cells were then treated with 10⁻⁹ M E2 (white) and different doses of glyceollin I (light grey) or II (hard grey) either alone (a), in combination with ICI₁₈₂,₇₈₀ (b) or in combination with E2 (c). Moreover, cells were pretreated with 10⁻⁵ M G1 or GII at different time indicated (3 h, 1 h and 0 h) before 10⁻⁹ M E2 adding (d). The results are expressed as the percentage of luciferase activity obtained with E2 treatment and are represented as the mean ± SEM of four independent experiments. \* *p*-value < 0.05, \*\* *p*-value < 0.01 and \*\*\* *p*-value < 0.001 with a Mann-Whitney test followed by Bonferroni correction for comparisons of the control vs the treatments (a) or for treatment with E2 vs the other compounds (c and d)</figcaption></figure>
induction. In contrast, it appears that the glyceollins acted additively with E2 to significantly increase luciferase activity; these increases were observed with 10⁻⁶ M and 5 × 10⁻⁶ M glyceollin I (*p* < 0.05) and with 5 × 10⁻⁵ M and 10⁻⁵ M glyceollin II (*p* < 0.01 and *p* < 0.001, respectively). In a recent study [23], it was shown an antagonist effect of glyceollins on ER activation in presence of E2 when cells were pretreated 1 h by glyceollins before E2 adding. Thus, we have pretreated cells with 10⁻⁵ M glyceollin I or II for 1 h or 3 h before adding 10⁻⁹ M E2 (Fig. 2d). Neither 3 h nor 1 h of pretreatment with GI affected E2-induction of luciferase activity. Like for GI, GII pretreatment did not antagonize E2 effect but the additive effect observed in absence of pretreatment was abolished.

## Glyceollins modulate E2-related gene expression

We also performed real-time PCR to evaluate the effects of glyceollins on different endogenous E2-dependent genes in MCF-7 cells. As shown in Fig. 3, expression levels of the progesterone receptor (PgR) (Fig. 3a), growth regulation by estrogen in breast cancer 1 (GREB1) (Fig. 3b), trefoil factor 1 (TFF1, also known as pS2) (Fig. 3c) and amphiregulin (AREG) (Fig. 3d) were measured after 24 h of treatment with 10⁻⁹ M E2 and different doses of glyceollin I or II either alone or in combination with E2. Different profiles were observed

depending on the considered gene. Concerning PgR expression, the glyceollins did not show any effects when used alone. When cells were co-treated with E2 and glyceollins, a 50% decrease in PgR expression was observed with 10⁻⁵ M glyceollins compared to that obtained in cells treated with E2 alone. Concerning GREB1 expression, we observed a 2-fold increase in expression with glyceollin I treatment compared to the control cells. Conversely, glyceollin II did not show any effects. Co-treatment of cells with E2 and glyceollins did not considerably affect E2-mediated induction of GREB1 expression, except for a slight repression with 10⁻⁵ M glyceollin II. For TFF1, we observed a similar profile to that obtained with GREB1: neither glyceollin I nor glyceollin II significantly affected TFF1 gene expression regardless of the presence of E2. A third type of expression profile was obtained with the AREG gene, which is expressed at very low basal levels in MCF-7 cells. E2 treatment strongly induced AREG gene expression, increasing it by approximately 40-fold. Surprisingly, 10⁻⁵ M glyceollin I or II alone also induced AREG gene expression. In addition, at 10⁻⁵ M, both glyceollins showed an additive effect with E2. To assess if the difference observed in the transcription of these genes could be linked to the recruitment of ERα on ERE, chromatin immunoprecipitation (ChIP) experiments were performed on GREB1 and PgR gene promoters. Additional file 2. As shown in the Additional file 3: Figure S2, E2-treatment increased ERα

<figure><img src="image_15.png" /></figure>

<figure><img src="image_16.png" /></figure>

<figure><img src="image_17.png" /></figure>

<figure><img src="image_18.png" /></figure>

**Fig. 3** Effects of glyceollins I and II on the expression of endogenous E2 target genes. MCF-7 cells were treated with different doses of glyceollin I (white) or II (black) with or without 10⁻⁹ M E2. The relative expression of PgR (a), GREB1 (b), TFF1 (c) and AREG (d) was assessed by real-time PCR and normalized to the expression of the housekeeping genes GAPDH and TBP. The results are expressed as the percentage of relative expression of each transcript obtained in the E2-treated cells and are represented as the mean ± SEM of four independent experiments. \*p-value < 0.05, \*\*p-value < 0.01 and \*\*\*p-value < 0.001 with a Mann-Whitney test followed by Bonferroni correction for comparisons of the control vs the treatments. \#p-value < 0.05 with a Mann-Whitney test followed by Bonferroni correction for comparisons of the E2 treatment vs the other treatments
recruitment on GREB1 promoter and PgR enhancer about 3 and 10 times, respectively. MCF-7 cell treated with GI or GII alone showed a slight or no effect on ERα recruitment on GREB1 and PgR promoter, respectively. However, when cells were co-treated with E2 and glyceollins, ERα recruitment on GREB1 promoter was not modified whereas it was significantly reduced on PgR enhancer. Together, these results are in good agreement with the expression profiles of these genes observed from the Q-PCR data (Fig. 3a and b).

## Genome-wide analysis of glyceollin's effects

To more precisely understand the actions of glyceollins and to particularly assess how glyceollins affect E2 gene regulation, we performed transcriptomic analysis of MCF-7 cells treated with glyceollins either alone or in combination with E2. The method used to select differentially expressed genes is described in Fig. 4a and in materials and methods. After this initial sorting, all the selected transcripts were combined, for a total of 1965 transcripts. Then, LIMMA tests were applied with a *p*-value ≤0.05, which led to the selection of 1852 transcripts representing 1550 genes. These differentially expressed genes were clustered, and 8 expression patterns emerged (P1-P8) (Fig. 4b). All the genes from the different patterns are referenced in Additional file 4: Table S1. P1, P2, P3 and P4 corresponded to estrogen-regulated genes that were upregulated (P1, P2 and P4) or downregulated (P3). In these groups, the glyceollins exerted different effects. Glyceollin I and glyceollin II both exerted antiestrogenic activity in P1 and to a lesser extent in P2 when cells were co-treated with E2. Moreover, in P1, glyceollin reduced gene expression independently of E2 co-treatment. P4 contained genes that were induced by E2 and also by glyceollin I or II. Interestingly, when cells were co-treated with E2 and the glyceollins, an additive effect was observed between E2 and glyceollin on the induction of P4 genes. P3 corresponded to genes repressed by E2 and to a lesser extent by glyceollin I or II. In co-treated cells, the glyceollins partially restored the expression of these genes compared to E2 treatment alone. P5, P6, P7 and P8 corresponded to genes unaffected or slightly affected by E2 but regulated by the glyceollins. P5 and P6 included genes upregulated by glyceollin I or II. The difference between these patterns resided in the effect of glyceollin II on gene induction, which was globally less important in P6. P7 and P8 corresponded to genes downregulated by both glyceollins. The difference between these patterns was the effect of glyceollin II, which was a less powerful inhibitor in P7. Then, differentially expressed genes from each pattern were subjected to GO and pathway analysis (Fig. 4b and Additional file 5: Figure S3). Fig. 4b shows the most relevant significant biological processes, cellular components and pathways for patterns P1, P3, P5 and P6. For

each term, the number of genes associated with the term is indicated and compared to the number of genes expected by chance. Interestingly, P1 essentially included genes linked to the cell cycle (68 genes were associated, 7 genes were expected by chance; 68/7) and cell proliferation (32/11). The effects of the glyceollins on the genes from this group could partially explain the anti-proliferative action of these compounds. In P3, GO terms related to cell communication (64/42), cell surface (16/6) and secretion pathways prevailed. In P5, as described previously [24], genes linked to monocarboxylic acid metabolism (23/7) and long fatty-acyl-CoA biosynthesis (5/0) were found. Finally, genes associated with cholesterol (10/0) and steroid biosynthesis (5/0) were found in P6. A Venn diagram is shown in Additional file 6: Figure S4.

Regulatory network analysis combining the differentially expressed gene data as well as the protein-gene interaction data showed a small network with two principal communities represented by three TFs, FOXM1, HIF1α and EPAS1/HIF2α (Fig. 5). In the FOXM1 group, genes from P1 were overrepresented, as shown by light red circles, whereas in the community with HIFα and EPAS1/HIF2α, E2-independent genes from P5, P6 and P7 were more important. However, glyceollins might regulate cell functions through two pathways, one involving ER and FOXM1 and the other involving the HIF pathway.

## Genes linked to proliferation and growth arrest were differentially expressed

As described above, glyceollins I and II affected FOXM1 expression and, consequently, downstream genes involved in G2/M transition and mitosis, such as CENPA, CENPF, AURKB and PLK1 [25, 26] (Fig. 6 and Additional file 4: Table S1). However, to validate the results from our transcriptomic analysis, genes not involved in G2/M transition but linked to FOXM1 were chosen (Fig. 6). Expression levels of FOXM1, ERα, ERβ, FOS and PPARγ were assessed by real-time PCR. Cells were treated under the same conditions as those used for the transcriptomic analysis, and the results are expressed as the percentage of normalized expression following E2 treatment. FOXM1 (Fig. 6a) was not induced by E2 at 10⁻⁹ M but was significantly reduced by glyceollin I and II both in the presence (*p* < 0.01 and *p* < 0.001, respectively) and absence of E2 (*p* < 0.001). ERα was affected and linked to FOXM1 in our network analysis (Fig. 5). To further understand the antiestrogenic properties of the glyceollins, ERα expression was also measured. ERα levels were significantly reduced (*p* < 0.01) in E2-treated cells compared to control cells. Glyceollin I and II were both strong inhibitors of ERα expression. As shown in Fig. 6b, ERα expression was significantly reduced in glyceollin I- and glyceollin II-treated cells compared to E2-treated cells
<figure><img src="image_19.png" /></figure>

<figure><img src="image_20.png" /><figcaption>**Fig. 4** Transcriptomic analysis, selection and clustering of differentially expressed genes. MCF-7 cells were treated with vehicle, 10⁻⁹ M E2, 10⁻⁵ M glyceollin I and II, or a combination of E2 and each of the glyceollins. Total RNA was extracted, reverse-transcribed, labeled and spotted onto a DNA chip. To select differentially expressed genes, all pairwise conditions were performed. All intensity signals above the overall median and with a fold change ≥2 were selected from each comparison. Then, the selected transcripts were joined and submitted to LIMMA tests. Only transcripts with a *p*-value ≤0.05 were selected, for a total of 1550 genes (**a**). These genes were clustered into eight different expression patterns (P1-P8) (**b**). Each pattern shows the number of genes in parentheses. The conditions are classified according to control, E2, glyceollin I, glyceollin II, E2 and glyceollin I, and E2 and glyceollin II. The most relevant GO terms for each pattern are indicated at the right side of the panel. For example, we found a significant enrichment of genes associated with cell division in P1 (68 genes, 7 expected by chance), as indicated (68/7) in the figure</figcaption></figure>

regardless of whether they were treated with (*p* < 0.01) or
without E2 (*p* < 0.01 and *p* < 0.001, respectively). In con-
trast, ERβ expression was not affected by either E2 or gly-
ceollin I or II (Fig. 6c). FOS is involved in tumorigenesis
and ER signaling. FOS expression was significantly induced
(*p* < 0.01) by E2 treatment but not by treatments with

glyceollin I or II. Moreover, glyceollin II significantly
inhibited the E2-mediated induction of FOS expres-
sion (*p* < 0.001) (Fig. 6d). Glyceollin I also tended to
inhibit this E2-mediated induction, but the effect was
not statistically significant. Finally, as glyceollins have
been shown to regulate lipid metabolism [24], we also
<figure><img src="image_21.png" /></figure>

measured the expression of PPARγ, a member of the FOXM1 community. As shown in Fig. 6e, PPARγ expression was significantly downregulated in E2-treated cells compared to controls (*p* < 0.001), whereas it was significantly upregulated when cells were treated with glyceollin I or II either in the presence or absence of E2 (*p* < 0.001) compared to E2 alone.

Next, the expression of genes linked to the hub cre-
ated by HIF1α and EPAS1/HIF2α was also validated by
real-time PCR (Fig. 7). HIF1α was slightly but signifi-
cantly induced following E2 treatment compared to the
control (*p* < 0.01). Unlike glyceollin II, glyceollin I was
1.5-fold more potent in inducing HIF1α expression than
E2 (*p* < 0.001). In addition, when cells were co-treated
with E2 and glyceollin I, we observed a synergistic effect
resulting in 2-fold higher E2-mediated induction of
HIF1α expression (*p* < 0.001). Cells co-treated with E2
and glyceollin II showed significantly induced HIF1α ex-
pression (*p* < 0.01) compared to those treated with E2
alone (Fig. 7a). EGLN3, part of the HIF1α community,
encodes a protein involved in the degradation of HIF1α
under normoxic conditions and was downregulated by
glyceollin I and II (Additional file 4: Table S1). Further-
more, the expression of EPAS1/HIF2α was strongly re-
pressed by E2 by 5-fold (*p* < 0.001) compared to the
control. Glyceollin I and II did not affect EPAS1/HIF2α

expression compared to the control. In co-treated cells,
glyceollin I and II partially restored the expression of
this gene (Fig. 7b). A classical target of HIF1α is VEGFA.
Therefore, it was not surprising that this gene emerged
in our analysis. Compared to control cells, VEGFA ex-
pression was not induced by E2 treatment, but it was in-
duced by treatments with glyceollin I and II (3.5-fold
and 2.5-fold, respectively; *p* < 0.001) independently of E2
co-treatment (Fig. 7c). Finally, it is interesting to note
that our network analysis identified a few genes associ-
ated with stress response and growth arrest, such as
GADD45B or DDIT4 L (REDD2) and DDIT4 (REDD1).
This last gene was chosen to validate the transcriptomic
analysis results. DDIT4 expression was significantly re-
duced (*p* < 0.001) following E2 treatment compared to
the control and was significantly induced by glyceollin I
and II (3.5-fold and 3-fold, respectively; *p* < 0.001). In
co-treated cells, E2 had a slight impact on glyceollin I
and II-mediated induction of DDIT4 (6-fold and 5.5-
fold, respectively, compared to E2 alone; *p* < 0.001).

The final pathway from the transcriptomic data is the CXCR4/CXCR7/CXCL12 axis, which is associated with the orphan receptor COUP-TFI (Fig. 8). In previous studies [27, 28], we have reported on the importance of the role of COUP-TFI in this signaling axis, which modulates the proliferation and migration of breast cancer cells. Indeed, breast tumor aggressiveness has been associated with increased COUP-TFI expression, decreased CXCL12 and CXCR7 expression, and increased CXCR4 expression. However, it is interesting to note that treatment with glyceollin I or II significantly reduced the expression of COUP-TFI and CXCR4 by 2-fold (*p* < 0.001 and *p* < 0.01, respectively, compared to E2) with or without E2 co-treatment (Fig. 8a and b). As previously described, E2 significantly decreased CXCR7 expression compared to the control (*p* < 0.001) (Fig. 8c). Glyceollin I and II also repressed CXCR7 expression; however, glyceollin I was significantly less potent than glyceollin II. Furthermore, neither of the glyceollins had an effect on E2-mediated repression of CXCR7 expression when used in co-treatments. E2 increased the expression of CXCL12 by approximately 5-fold, whereas glyceollin I and II did not modify CXCL12 expression compared to the control. In contrast, both glyceollins markedly reduced the E2-mediated induction of CXCL12 expression (Fig. 8d).

Discussion

Estrogens are involved in multiple physiological pro-
cesses and act on various tissues. In particular, they par-
ticipate in the development and maintain the function of
reproductive organs such as the gonads or the mammary
gland through their binding to ERα and ERβ. ERα is the
major isoform in the mammary gland and has a role in
<figure><img src="image_22.png" /></figure>

<figure><img src="image_23.png" /></figure>

<figure><img src="image_24.png" /></figure>

<figure><img src="image_25.png" /></figure>

<figure><img src="image_26.png" /><figcaption>**Fig. 6** Validation of differentially expressed genes corresponding to the FOXM1 community. MCF-7 cells were treated with vehicle (black), 10⁻⁹ M E2 (white), 10⁻⁵ M glyceollin I (light grey) and II (hard grey), or a combination of E2 and each of the glyceollins (hatched squares). The relative expression of FOXM1 (a), ERα (b), ERβ (c), FOS (d) and PPARG (e) was assessed by real-time PCR and normalized to the expression of the housekeeping genes GAPDH and TBP. The results are expressed as the percentage of relative expression of each transcript obtained in E2-treated cells and are represented as the mean +/- SEM of 10 independent experiments. \*p-value <0.05, \*\*p-value <0.01 and \*\*\*p-value <0.001 with a Mann-Whitney test followed by Bonferroni correction for comparisons of the control vs the treatments. \#p-value <0.01 and \##p-value <0.001 with a Mann-Whitney test followed by Bonferroni correction for comparisons of E2 vs the other treatments</figcaption></figure>

epithelial duct proliferation and differentiation [29]. E2 also promotes cell survival, and due to its proliferative effect, this hormone has been linked to breast cancer [3]. Many natural and synthetic chemicals in the environment have been reported to exhibit hormonal activity, particularly estrogenic potency [30]. This is the case for the well-known compound bisphenol A, which has well-documented effects on breast cancer [31]. Phytoestrogens, which are estrogenic compounds from plants, are also found in food, particularly in soy, and have been reported to decrease the risk of breast cancer at high doses

[32, 33]. Among the phytoestrogens, glyceollins emerged as promising compounds in the 2000s.

A previous study from Burow et al. [11] reported that glyceollins bind the ER and act as antiestrogenic compounds by inhibiting cell proliferation. However, the detailed molecular mechanisms driving the anti-proliferative actions of these phytochemicals remain elusive and appear to be more complex than those based only on ER interactions [14]. In the present work, we utilized synthetic glyceollins I and II to better delineate the modes of action exhibited by these compounds. We showed that glyceollin
<figure><img src="image_27.png" /></figure>

<figure><img src="image_28.png" /></figure>

<figure><img src="image_29.png" /></figure>

<figure><img src="image_30.png" /></figure>

**Fig. 7** Validation of differentially expressed genes corresponding to the HIF1α/EPAS1 community. MCF-7 cells were treated with vehicle (black), 10⁻⁹ M E2 (white), 10⁻⁵ M glyceollin I (light grey) and II (hard grey), or a combination of E2 and each of the glyceollins (hatched squares). The relative expression of HIF1α (a), EPAS1 (b), VEGFA (c) and DDIT4 (d) was assessed by real-time PCR and normalized to the expression of the housekeeping genes GAPDH and TBP. The results are expressed as the percentage of relative expression of each transcript obtained in E2-treated cells and are represented as the mean ± SEM of 10 independent experiments. **p**-value <0.01 and ***p**-value <0.001 with a Mann-Whitney test followed by Bonferroni correction for comparisons of the control vs the treatments. ##p**-value <0.01 and ###p**-value <0.001 with a Mann-Whitney test followed by Bonferroni correction for comparisons of E2 vs the other treatments

I and glyceollin II inhibit the trophic action of E2 in vivo in mouse mammary glands, but not in uteri. Our data differ from those obtained by Salvo et al., who reported that natural glyceollins antagonize the trophic effects of E2 in uteri [13]. One explanation for these different observations regarding uterotrophy is that the referenced study used nude mice that were treated daily with the natural glyceollins for twenty days, whereas our current study treated mice for only three days. Nevertheless, our data clearly show that mammary epithelium growth is not influenced by glyceollins when they are administered alone; however, when administered together with E2, they are capable of inhibiting the stimulatory effect of E2 on ductal epithelium growth. It should be noted that since proliferation is partially blocked by glyceollins during ductal elongation, further in vivo experiments, such as TUNEL assays, would be required to determine other mechanisms of action by which glyceollins may block E2-mediated ductal elongation. However, to the best of our knowledge, this is the first in vivo study showing the anti-proliferative effect of glyceollins on epithelial ductal extension. In this way, glyceollins appear to be selective estrogen receptor modulators (SERM). In accordance with our results, a previous study investigated the effects of various SERM such as raloxifen, on mammary gland development [34]. The authors showed that raloxifen alone induced a slight ductal

tree invasion in the fat pad, whereas, in combination with conjugated estrogens, raloxifen had a clear antagonistic effect [34]. Next, we determined the in vitro effects of glyceollins on cell proliferation. Interestingly, we found that glyceollins I and II exert anti-proliferative effects in both ER-positive and ER-negative breast cancer cells in accordance with the study of Rhodes et al. [35]. This suggests that glyceollins do not act as conventional antiestrogens, such as tamoxifen, but rather act through both ER-dependent and ER-independent pathways, although the ER-dependent pathway seems to be predominant. Burow and collaborators reported that glyceollins could antagonize the E2-mediated stimulation of an ERE-luciferase reporter plasmid [11, 12, 15]. Although this was not observed in our study, the use of different ERE sequences and cell lines and differences in the duration of the treatment could account for this discrepancy. For example, in our experiments, an ERE-luciferase reporter plasmid was used that contains only a single ERE sequence upstream of the luciferase gene, whereas Burow and collaborators used a luciferase reporter with two ERE motifs in addition to pre-treating cells with glyceollins before E2 was added. Based on these observations, we suggest that glyceollins may prevent cooperative effects between ER dimers on ERE sequences, which could explain the decreased luciferase activity. Nevertheless,
<figure><img src="image_31.png" /></figure>

<figure><img src="image_32.png" /></figure>

<figure><img src="image_33.png" /></figure>

<figure><img src="image_34.png" /></figure>

**Fig. 8** Validation of differentially expressed genes corresponding to the CXCL12/CXCR4/CXCR7 axis. MCF-7 cells were treated with vehicle (black), 10⁻⁹ M E2 (white), 10⁻⁵ M glyceollin I (light grey) and II (hard grey), or a combination of E2 and each of the glyceollins (hatched squares). The relative expression of COUP-TFI (a), CXCR4 (b), CXCR7 (c) and CXCL12 (d) was assessed by real-time PCR and normalized to the expression of the housekeeping genes GAPDH and TBP. The results are expressed as the percentage of relative expression of each transcript obtained in E2-treated cells and are represented as the mean +/- SEM of 10 independent experiments. \*p-value <0.05, \*\*p-value <0.01 and \*\*\*p-value <0.001 with a Mann-Whitney test followed by Bonferroni correction for comparisons of the control vs the treatments. \#p-value <0.01 and \####p-value <0.001 with a Mann-Whitney test followed by Bonferroni correction for comparisons of E2 vs the other treatments

glyceollin I and II inhibited the expression of the endogenous PgR gene induced by E2 by over 50% in MCF-7 cells.

To further explore the molecular mechanisms underlying how glyceollins exert their anti-proliferative effects, we performed transcriptomic analysis of MCF-7 breast cancer cells exposed to glyceollins and created a gene regulatory network of differentially expressed genes. This integrative genomic approach was followed by the quantification of several key genes, which allowed us to identify, for the first time, two major pathways involving the ER and FOXM1 factors and the other including the hypoxia inducible factor (HIF) family (HIF1α and EPAS1/ HIF2α). The first hub highlighted in our gene regulatory network is represented by the forkhead transcription factor FOXM1 and the ER. FOXM1 is a well-known key regulator of the cell cycle and is involved in G1/S and G2/M transition [36]. Thus, the downregulation of this gene could explain the effects of glyceollins on cellular proliferation. FOXM1 gene expression involves ERα [37], and in return, expression of ERα involves FOXM1 [38]. By targeting ERα, glyceollins could affect this auto-regulatory loop. The mechanisms by which glyceollins act through ERβ are not fully defined. Competition

binding assays showed that glyceollins are able to bind both ERα and ERβ with, however, a greater sensitivity of glyceollins for ERα vs. ERβ [11]. In addition, since MCF-7 cells express mainly ERα (the ratio ERα/ERβ is 8/1) [5], the effects of glyceollins are likely mediated by ERα signaling. Nevertheless, glyceollins may affect ERα/FOXM1 regulatory loop by another pathway that may potentially involve ERβ. Indeed, ERβ represses FOXM1 expression by displacing ERα from the FOXM1 promoter [39]. Our data showed that glyceollins markedly decreased ERα expression but do not affect ERβ expression in MCF-7 cells. A change in the equilibrium of the ERα / ERβ ratio could then contribute to the antiestrogenic activity exerted by the glyceollins and could reinforce the possibility of an involvement of the ERβ. In addition to the auto-regulatory loop that exists between ERα and FOXM1, decreased expression and activity of these two factors could explain the downregulation of GREB1, at least with glyceollin II. Indeed, FOXM1 and ERα co-bind DNA in breast cancer cells and modulate the expression of specific genes [40]. In the referenced work, the authors showed that FOXM1 knockout affected GREB1 expression. Therefore, one could easily
hypothesize that glyceollins inhibit E2-related gene expression via this pathway. Moreover, overexpression of FOXM1 is a hallmark of many cancers and a sign of poor prognosis. In ER-positive breast cancer, overexpression of FOXM1 is associated with endocrine resistance and invasiveness because it favors the expansion of stem-like cancer cells [26]. A recent study showed that the FOXM1 cistrome is a powerful index to predict breast cancer outcomes [41]. Thus, it will be very interesting to test the plasticity of the binding interaction that exists between FOXM1 and ERα in response to glyceollin treatment.

The second hub highlighted in our gene regulatory network is centered on the HIF family. The HIF family is composed of three O₂-regulated members (HIF1α, EPAS1/HIF2α and HIF3α) that become stabilized under hypoxic conditions. To accomplish this, they heterodimerize with the constitutively expressed HIF1β (also known as ARNT) to regulate genes necessary for adaptation to low-oxygen conditions [42]. It was surprising that the glyceollins in this study induced HIF1α expression under normoxic conditions because they have been previously described as inhibitors of this factor at both the synthesis and stability levels under hypoxic conditions [43]. Under normoxic conditions, HIF1α is controlled by numerous stimuli, including reactive oxygen species (ROS) [44]. Recently, it was shown that glyceollin at a concentration of approximately 18 µM induced ROS generation in a hepatic cell line [45]. In our experiments utilizing 10 µM glyceollin, it is possible that moderate ROS production was induced that consequently induced HIF family activity. We identified DDIT4 (also known as REDD1) and DDIT4L (also known as REDD2), both inhibitors of mTORC1 [46], in the HIF family community. REDD1 and REDD2 are stress-responsive genes induced by different stimuli, such as DNA damage or hypoxia. mTORC1 is a member of the PI3K/AKT signaling pathway and acts downstream of AKT. It participates in protein synthesis by promoting the phosphorylation of p70S6K. Therefore, glyceollins alter the phosphorylation of p70S6K in ER-positive breast cancer [23]. Thus, inhibition of mTORC1 could be a factor involved in the anti-proliferative effects produced by glyceollins due to perturbations in the PI3K/AKT/mTOR pathway. PI3K mutations are frequently observed in ER-positive breast cancer. Many inhibitors of PI3K pathway are under clinical trials or approved as therapeutics such as everolimus which is a mTORC1 inhibitor [47]. Thus, this observation reinforces the therapeutic potential of glyceollin in ER-positive breast cancer. Moreover, another study reported that loss of the REDD1 gene leads to an increase in HIF1 level and consequently an increase in tumorigenicity. The authors also showed that REDD1 localizes to mitochondria to regulate ROS production [48]. Overall, REDD1 appears to act as a tumor suppressor that

works through different levels, reinforcing the therapeutic potential of glyceollins.

Glyceollins are studied in part for their ability to inhibit E2-related gene expression [12]. Thus, we were surprised to note that, unlike the other genes that were assessed, AREG expression was induced by 10⁻⁵ M glyceollins, and this effect was increased by E2 co-treatment. AREG is regulated by numerous transcription factors, including the ER [49]. Recently, a role for EPAS1/HIF2α in the induction of AREG expression was described in MCF-7 cells [50]. Thus, glyceollins might induce AREG expression through ERs and EPAS1/HIF2α, which would explain the synergistic effect. Moreover, high expression of EPAS1/HIF2α, AREG and WISP2 is linked to improved survival in breast cancer [50]; glyceollin treatment does not affect EPAS1/HIF2α expression in the absence of E2 treatment and even partially restores expression in E2-treated cells. AREG and WISP2 (Fig. 3 and Additional file 4: Table S1) are overexpressed in glyceollin-treated cells.

Finally, our transcriptomic analysis identified differential effects on the expression of the orphan receptor COUP-TFI and the chemokine CXCL12, as well as its receptors CXCR4 and CXCR7, following glyceollin treatment. Considering the important role of COUP-TFI in CXCL12 expression and the importance of the CXCL12 signaling axis in tumor growth and metastasis, the effects produced by glyceollins seem very important. Indeed, the chemokine CXCL12 plays critical roles in cell migration, angiogenesis, proliferation, and survival in many types of cancer, including breast cancer, by interacting with the transmembrane receptors CXCR4 and CXCR7 [51, 52]. CXCR4 is often overexpressed in metastatic tumors and promotes the migration of invasive cells to tissues where local CXCL12 secretion is increased, such as bone, liver, brain and lung [53, 54]. We recently reported that E2 controls the activity of the CXCL12/CXCR4/CXCR7 signaling axis in breast tumor cells and influences the proliferation and migration of breast cancer cells [27, 28]. Furthermore, COUP-TFI and the CXCL12 signaling axis are dysregulated in breast tumor biopsies compared to normal epithelium. Indeed, primarily in ER-positive invasive ductal cancer, we observed significant upregulation of COUP-TFI and CXCR4 and downregulation of CXCR7 and CXCL12, and the levels of these changes showed correlations with tumor grade [28]. Downregulation of CXCL12 in cancer cells is frequently associated with promoter methylation, which encourages cells to migrate toward a CXCL12 gradient and establish metastases [55]. Interestingly, glyceollins repress the expression of CXCR4 and do not affect CXCR7. However, they exert antiestrogenic activity in E2-mediated induction of CXCL12 and they maintain the expression of this gene; thus, they might limit
the metastatic potential of tumor cells. In accordance with our observation, a previous study showed that glyceollins could reverse the epithelial to mesenchymal transition of letrozole resistant cells and thus decreased their invasion and migration [56]. It would therefore be interesting to test the ability of glyceollins to limit the loss of expression of key genes in cancer and, in particular, the activation of enzymes involved in epigenetic modifications, as was previously described for genistein [57].

## Conclusion

In conclusion, glyceollins I and II did not show any effect on mouse uterotrophy, whereas they did exert antiproliferative effects on mammary gland epithelial duct growth. This antagonistic activity was confirmed in different ER-positive and ER-negative breast cell lines. Our mechanistic studies revealed that glyceollins are more cytostatic than cytotoxic. Moreover, they have some similarity to SERMs which have partial agonist and antagonist properties depending on E2-target genes. For the first time, a genome-wide microarray was performed on an ER-positive breast cell line to identify pathways involved in the anti-proliferative effects of glyceollins. We identified two major pathways, centered on FOXM1/ERα and HIF1α/HIF2α, which could explain the activity of glyceollins on ER-positive and ER-negative cell lines. These results confirm and reinforce the therapeutic potential of glyceollins for managing breast cancer.

## Additional files

**Additional file 1: Figure S1.** Effect of glyceollin I and II on ovariectomized mouse uterotrophy and on cell cycle and apoptosis in MCF-7 cells. Uteri were obtained from ovariectomized (ovx) or intact (non-ovx) mice treated with vehicle, 10 µg/kg E2, 50 mg/kg glyceollin I or glyceollin II (GII) (GI L or GII L), or 100 mg/kg glyceollin I or glyceollin II (GI H or GII H) either alone or in combination with E2 for 72 h (Additional file 2: Figure S1a). The mice were then sacrificed, and their uteri were removed and weighted. The results are represented in box-and-whisker plots, where the top and the bottom of the box correspond to the 75th and the 25th percentile, respectively. The horizontal bar in the box is the median, and the points outside the box correspond to extreme values. The results are expressed as relative uteri weight (g per g of body weight) and were taken from 4 independent experiments with at least 5 mice per group. \*\*\*p-value <0.001 with a Mann-Whitney test followed by Bonferroni correction for comparisons of the control vs the treatments. \#\#p-value <0.01 with a Mann-Whitney test followed by Bonferroni correction for comparisons of E2 vs the other treatments. For analyses of cell cycle (Additional file 2: Figure S1b) and apoptosis (Additional file 2: Figure S1c), cells were treated for 3 days with 10⁻⁹ E2 with or without 10⁻⁵ M glyceollin I or II. To analyze cell cycle, cells were stained with propidium iodide and subjected to flow cytometry analysis. The results are expressed as the percentages of cells in each cell cycle phase and are represented as the mean of 4 independent experiments +/- SEM. \*p-value <0.05, \*\*p-value <0.01 and \*\*\*p-value <0.001 with a Mann-Whitney test followed by Bonferroni correction for comparisons of the control vs the treatments. For analysis of apoptosis, cells were stained using a TUNEL assay, and the percentage of apoptotic cells was assessed with an Array Scan VTI. The results are expressed as the percentage of TUNEL-positive cells compared to total cells. (TIFF 2098 kb)

**Additional file 2: Supplemental material. (PDF 157 kb)**

**Additional file 3: Figure S2.** Effect of glyceollins on ERα recruitment on GREB1 promoter and PgR enhancer. MCF-7 cells were treated with vehicle

(black), 10⁻⁹ M E2 (white), 10⁻⁵ M glyceollin I (light grey) and II (hard grey), or a combination of E2 and each of the glyceollins (hatched squares). The recruitment of ERα on GREB1 promoter (a) and PgR enhancer (b) was assessed by chromatin immunoprecipitation followed by real time PCR. Results are expressed in fold recruitment compared to control and are the mean of two independent experiments. (TIFF 940 kb)

**Additional file 4: Table S1.** All differentially expressed genes. The ID of all genes are indicated with their functions and expression levels. *Cell cycle analysis.* MCF-7 cells (100,000 cells/well) were plated in 6-well plates. After 72 h of serum and steroid deprivation, the cells were treated for 72 h with solvent as control, 10⁻⁹ M E2, 10⁻⁵ M glyceollin I or II, or a combination of E2 and glyceollin I or II. Then, the cells were trypsinized and fixed in 70% ethanol before staining with propidium iodide. The percentage of cells in each cell cycle phase was assessed by flow cytometry with a FACS Calibur (BD Biosciences). *Measurement of apoptosis* MCF-7 cells (4000 cells/well) were plated in 96-well plates. After 72 h of serum and steroid deprivation, the cells were treated for 72 h with solvent as control, 10⁻⁹ M E2, 10⁻⁵ M glyceollin I or II, or a combination of E2 and glyceollin I or II. TUNEL staining was assessed with an In Situ Cell Death Detection Kit, Fluorescein (Roche) according to the manufacturer's instructions. The fluorescence and percentage of TUNEL-positive cells were determined with an Array Scan VTI (Thermo Fisher Scientific) on the ImPACcell platform (Rennes, France). *Chromatin Immunoprecipitation (ChIP)* MCF-7 cells (2,000,000 cells/dishes) were plated in 10 cm dishes and then deprived of steroids and serum for 72 h. The cells were treated for 1 h with 10⁻⁹ M E2, with 10⁻⁵ M Gl or GII with or without 10⁻⁹ M E2. Then, cells were cross-linked for 10 min with 1.5% of formaldehyde (Sigma). Cells were lysed in lysis buffer (50 mM Tris-HCl, pH 8.1, 10 mM EDTA, 0.5% Empigen BB and 1% SDS). Chromatin was sonicated 10 min (15 s on/off cycles) on Bioruptor (Diagenode) at highest intensity. Soluble chromatin was diluted in IP buffer (20 mM Tris-HCl, pH 8.1, 2 mM EDTA, 0.1% Triton X-100) with 2 µg of ERα antibody (E115, Abcam) and yeast RNA as non-specific competitor and incubated overnight at 4 °C on rocking platform. Then, protein G coupled sepharose beads were added to the samples and were incubated 4 h à 4 °C. Immune complexes were washed one time in washing buffer 1 (20 mM Tris-HCl, pH 8.1, 2 mM EDTA, 150 mM NaCl, 1% Triton X-100 and 0.1% SDS), one time in washing buffer 2 (20 mM Tris-HCl, pH 8.1, 2 mM EDTA, 500 mM NaCl, 1% Triton X-100 and 0.1% SDS), one time in washing buffer 3 (10 mM Tris-HCl, pH 8.1, 1 mM EDTA, 250 mM LiCl, 1% Deoxycholate and 1% NP-40) and finally two times in washing buffer 4 (10 mM Tris-HCl, pH 8.1, 1 mM EDTA). After washing, immune complexes were extracted with 100 µl of extraction buffer (0.1 M NaHCO₃ and 1% SDS). Cross-linking was reverse by incubation of samples overnight at 65 °C and DNA was purified using the Nucleospin Gel and PCR cleanup kit (Macherey Nagel). Enrichment analysis on the ERE proximal of GREB1 (Fwd: CACTTTGAGCAAAAGCCACA and Rev.: GACCCAGTTGCCACACTTTT) and on an enhancer 1 of PgR described in [58] was normalized using an irrelevant region on the chromosome 10 (Fwd: AGGTGACAAGCCAAGTGTCC and Rev.: GCCTGGTGGCATACTAAAGG). Analysis was performed by real time PCR on a CFX 384 apparatus (BioRad) on 2 µL of immunoprecipitation or 0.2 µL of input with 500 nM of primers and iTaq Universal SYBR Green Supermix (BioRad). (XLSX 590 kb)

**Additional file 5: Figure S3.** GO enrichment analysis of different treatment-related expression patterns. Eight expression patterns are matched with a selection of GO terms from the ontology "phenotypes," "biological process," "cellular component" and "pathways." The numbers of genes associated with each GO term are indicated in the first column. Enrichment is indicated by bolded rectangles, where the first number indicates the number of genes found in our analysis and the second the number expected with a random list of genes. Overrepresented genes in a specific GO term are shown in red, and underrepresented genes are shown in blue. (TIFF 2724 kb)

**Additional file 6: Figure S4.** Venn diagram. A Venn diagram was created from the list of differentially expressed genes obtained from comparisons of the control and E2 (red), Gl (yellow), GII (green), E2 + Gl (blue) and E2 + GII (purple) treatments. (TIFF 3761 kb)

## Abbreviations

E2: Estradiol; ER: Estrogen receptor; ERE: Estrogen responsive element; LIMMA: Linear model for microarray data; OVX: Ovariectomized; SERM: Selective estrogen receptor modulator
## Acknowledgments

We kindly thank Catherine Martin, Adrien Alusse and Charly Jehano for their technical assistance with mouse experiments and Rémy Le Guevel for his assistance with TUNEL analysis. We also thank the University of Rennes 1, Inserm and CNRS for supporting our research programs.

## Availability and data materials

All data generated or analyzed during this study are included in this published article [and its additional information files].

## Funding

This work was funded by FUI mVolio, Région Bretagne and Rennes Métropole.

## Authors' contributions

SL conducted, analyzed and interpreted all cellular experiments and was a major contributor in writing the manuscript. FC analyzed the transcriptomic data. FF and FPE conducted the animal experiments. FF, CSa and ML participated to the in vivo data analyses. NP, CSu and TE produced and characterized the glyceollins. FPA supervised and coordinated the study, participated to the conception and design of data and was a major contributor in writing the manuscript. All authors read and approved the final manuscript.

## Ethics approval and consent to participate

All animal procedures were performed according to the guidelines for animal models in research defined by the Ethics Committee and approved by the Ministry of France (reference project number, 2015061812074056_V2). All experiments were conducted by FF and FPE, who are qualified in laboratory animal care and use procedures.

## Consent for publication

Not applicable

## Competing interests

The authors declare that they have no conflicts of interest.

## Author details

¹Institut de Recherche en Santé-Environnement-Travail (IRSET), University of Rennes 1, 9 Avenue du Pr Léon Bernard, 35000 Rennes, France. ²Inserm U1085, Team Transcription, Environment and Cancer, 9 Avenue du Pr Léon Bernard, 35000 Rennes, France. ³Inserm U1085, Team Viral and Chemical Environment & Reproduction, 9 Avenue du Pr Léon Bernard, 35000 Rennes, France. ⁴Laboratoire Nutrinov, Technopole Atalante Champeaux, 8 rue Jules Maillard de la Gournerie, 35012 Rennes Cedex, France.

Received: 10 March 2017 Accepted: 20 June 2017

Published online: 30 June 2017

## References

1. WHO, Globocan 2012. Fact Sheets by Cancer. Available from: http://globocan.iarc.fr/Pages/fact_sheets_cancer.aspx

2. Levin ER, Pietras RJ. Estrogen receptors outside the nucleus in breast cancer. Breast Cancer Res Treat. 2008;108:351-61.

3. Kerdivel G, Flouriot G, Pakdel F. Modulation of estrogen receptor alpha activity and expression during breast cancer progression. Vitam Horm. 2013; 93:135-60.

4. Sirokin AV, Harrath AH. Phytoestrogens and their effects. Eur J Pharmacol. 2014;741:230-6.

5. Lecomte S, Lelong M, Bourgine G, Efstathiou T, Saligaut C, Pakdel F. Assessment of the potential activity of major dietary compounds as selective estrogen receptor modulators in two distinct cell models for proliferation and differentiation. Toxicol Appl Pharmacol. 2017;325:61-70.

6. Lozovaya VV, Lygin AV, Zernova OV, Li S, Hartman GL, Widholm JM. Isoflavonoid accumulation in soybean hairy roots upon treatment with Fusarium solani. Plant Physiol Biochem. 2004;42:671-9.

7. Lygin AV, Zernova OV, Hill CB, Kholina NA, Widholm JM, Hartman GL, et al. Glyceollin is an important component of soybean plant defense against Phytophthora sojae and Macrophomina phaseolina. Phytopathology. 2013; 103:984-94.

8. Edens RM, Anand SC, Bolla RI. Enzymes of the Phenylpropanoid pathway in soybean infected with Meloidogyne Incognita or *Heterodera glycines*. J Nematol. 1995;27:292-303.

9. Ng TB, Ye XJ, Wong JH, Fang EF, Chan YS, Pan W, et al. Glyceollin, a soybean phytoalexin with medicinal properties. Appl Microbiol Biotechnol. 2011;90:59-68.

10. Nikov GN, Hopkins NE, Boue S, Alworth WL. Interactions of dietary estrogens with human estrogen receptors and the effect on estrogen receptor-estrogen response element complex formation. Environ Health Perspect. 2000;108:867-72.

11. Burow ME, Boue SM, Collins-Burow BM, Melnik LI, Duong BN, Carter-Wientjes CH, et al. Phytochemical glyceollins, isolated from soy, mediate antihormonal effects through estrogen receptor alpha and beta. J Clin Endocrinol Metab. 2001;86:1750-8.

12. Zimmermann MC, Tilghman SL, Boué SM, Salvo VA, Elliott S, Williams KY, et al. Glyceollin I, a novel antiestrogenic phytoalexin isolated from activated soy. J Pharmacol Exp Ther. 2010;332:35-45.

13. Salvo VA, Boué SM, Fonseca JP, Elliott S, Corbitt C, Collins-Burow BM, et al. Antiestrogenic glyceollins suppress human breast and ovarian carcinoma tumorigenesis. Clin Cancer Res. 2006;12:7159-64.

14. Khupse RS, Sarver JG, Trendel JA, Bearss NR, Reese MD, Wiese TE, et al. Biomimetic syntheses and antiproliferative activities of racemic, natural (-), and unnnatural (+) glyceollin I. J Med Chem. 2011;54:3506-23.

15. Payton-Stewart F, Khupse RS, Boué SM, Elliott S, Zimmermann MC, Skripnikova EV, et al. Glyceollin I enantiomers distinctly regulate ER-mediated gene expression. Steroids. 2010;75:870-8.

16. Tian J, Berton TR, Shirley SH, Lambertz I, Gimenez-Conti IB, DiGiovanni J, et al. Developmental stage determines estrogen receptor alpha expression and non-genomic mechanisms that control IGF-1 signaling and mammary proliferation in mice. J Clin Invest. 2012;122:192-204.

17. Khupse RS, Erhardt PW. Total syntheses of Racemic, natural (-) and unnnatural (+) Glyceollin I. Org Lett. 2008;10:5007-10.

18. Luniwal A, Khupse R, Reese M, Liu J, El-Dakdouki M, Malik N, et al. Multigram synthesis of Glyceollin I. Org Process Res Dev. 2011;15:1149-62.

19. Chalmel F, Primig M. The annotation, mapping, expression and network (AMEN) suite of tools for molecular systems biology. BMC Bioinformatics. 2008;9:86.

20. Smyth GK. Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat. Appl. Genet. Mol. Biol. 2004;3:1-25.

21. Yusuf D, Butland SL, Swanson MI, Bolotin E, Ticoll A, Cheung WA, et al. The transcription factor encyclopedia. Genome Biol. 2012;13:R24.

22. Takahashi O, Oishi S. Male reproductive toxicity of four bisphenol antioxidants in mice and rats and their estrogenic effect. Arch Toxicol. 2006; 80:225-41.

23. Bratton MR, Martin EC, Elliott S, Rhodes LV, Collins-Burow BM, McLachlan JA, et al. Glyceollin, a novel regulator of mTOR/p70S6 in estrogen receptor positive breast cancer. J Steroid Biochem Mol Biol. 2015;150:17-23.

24. Wood CE, Boue SM, Collins-Burow BM, Rhodes LV, Register TC, Cline JM, et al. Glyceollin-elicited soy protein consumption induces distinct transcriptional effects as compared to standard soy protein. J Agric Food Chem. 2012;60:81-6.

25. Stoll SW, Stuart PE, Swindell WR, Tsoi LC, Li B, Gandarillas A, et al. The EGF receptor ligand amphiregulin controls cell division via FoxM1. Oncogene. 2015;35:2075-86.

26. Bergamaschi A, Madak-Erdogan Z, Kim YJ, Choi Y-L, Lu H, Katzenellenbogen BS. The forkhead transcription factor FOXM1 promotes endocrine resistance and invasiveness in estrogen receptor-positive breast cancer by expansion of stem-like cancer cells. Breast Cancer Res. 2014;16:436.

27. Boudot A, Kerdivel G, Habauzit D, Eeckhoute J, Le Dily F, Flouriot G, et al. Differential estrogen-regulation of CXCL12 chemokine receptors, CXCR4 and CXCR7, contributes to the growth effect of estrogens in breast cancer cells. PLoS One. 2011;6:e20898.

28. Boudot A, Kerdivel G, Lecomte S, Flouriot G, Desille M, Godey F, et al. COUP-TFI modifies CXCL12 and CXCR4 expression by activating EGF signaling and stimulates breast cancer cell migration. BMC Cancer. 2014;14:407.

29. Couse JF, Korach KS. Reproductive phenotypes in the estrogen receptor-alpha knockout mouse. Ann Endocrinol. 1999;60:143-8.

30. Kerdivel G, Habauzit D, Pakdel F. Assessment and molecular actions of endocrine-disrupting chemicals that interfere with estrogen receptor pathways. Int J Endocrinol. 2013;2013:501851.
31. Seachrist DD, Bonk KW, Ho S-M, Prins GS, Soto AM, Keri RA. A review of the carcinogenic potential of bisphenol a. Reprod Toxicol. 2015;59:167-82.

32. Wu AH, Yu MC, Tseng C-C, Pike MC. Epidemiology of soy exposures and breast cancer risk. Br J Cancer. 2008;98:9-14.

33. Caan BJ, Natarajan L, Parker B, Gold EB, Thomson C, Newman V, et al. Soy food consumption and breast cancer prognosis. Cancer Epidemiol Biomark Prev. 2011;20:854-8.

34. Peano BJ, Crabtree JS, Komm BS, Winneker RC, Harris HA. Effects of various selective Estrogen receptor modulators with or without conjugated estrogens on mouse mammary gland. Endocrinology. 2009;150:1897-903.

35. Rhodes LV, Tilghman SL, Boue SM, Wang S, Khalili H, Muir SE, et al. Glyceollins as novel targeted therapeutic for the treatment of triple-negative breast cancer. Oncol Lett. 2012;3:163-71.

36. Wang I-C, Chen Y-J, Hughes D, Petrovic V, Major ML, Park HJ, et al. Forkhead box M1 regulates the transcriptional network of genes essential for mitotic progression and genes encoding the SCF (Skp2-Cks1) ubiquitin ligase. Mol Cell Biol. 2005;25:10875-94.

37. Millour J, Constantinidou D, Stavropoulou AV, Wilson MSC, Myatt SS, Kwok JM-M, et al. FOXM1 is a transcriptional target of ERalpha and has a critical role in breast cancer endocrine sensitivity and resistance. Oncogene. 2010; 29:2983-95.

38. Madureira PA, Varshochi R, Constantinidou D, Francis RE, Coombes RC, Yao K-M, et al. The Forkhead box M1 protein regulates the transcription of the estrogen receptor alpha in breast cancer cells. J Biol Chem. 2006;281:25167-76.

39. Horimoto Y, Hartman J, Millour J, Pollock S, Olmos Y, Ho K-K, et al. ERβ1 represses FOXM1 expression through targeting ERα to control cell proliferation in breast cancer. Am J Pathol. 2011;179:1148-56.

40. Sanders DA, Ross-Innes CS, Beraldi D, Carroll JS, Balasubramanian S. Genome-wide mapping of FOXM1 binding reveals co-binding with estrogen receptor alpha in breast cancer cells. Genome Biol. 2013;14:R6.

41. Yau C, Meyer L, Benz S, Vaske C, Scott G, Egan B, et al. FOXM1 cistrome predicts breast cancer metastatic outcome better than FOXM1 expression levels or tumor proliferation index. Breast Cancer Res Treat. 2015;154:23-32.

42. Gilkes DM, Semenza GL. Role of hypoxia-inducible factors in breast cancer metastasis. Future Oncol. 2013;9:1623-36.

43. Lee S-H, Jee J-G, Bae J-S, Liu K-H, Lee YM. A group of novel HIF-1α inhibitors, glyceollins, blocks HIF-1α synthesis and decreases its stability via inhibition of the PI3K/AKT/mTOR pathway and Hsp90 binding. J Cell Physiol. 2015;230:853-62.

44. Ke Q, Costa M. Hypoxia-inducible factor-1 (HIF-1). Mol Pharmacol. 2006; 70:1469-80.

45. Kim HJ, Jung CL, Jeong YS, Kim J-S. Soybean-derived glyceollins induce apoptosis through ROS generation. Food Funct. 2014;5:688-95.

46. DeYoung MP, Horak P, Sofer A, Sgroi D, Ellisen LW. Hypoxia regulates TSC1/2 mTOR signaling and tumor suppression through REDD1-mediated 14 3 3 shuttling. Genes Dev. 2008;22:239-51.

47. Dey N, De P, Leyland-Jones B. PI3K-AKT-mTOR inhibitors in breast cancers: From tumor cell signaling to clinical trials. Pharmacol. Ther. 2017;175:91-106. [Epub ahead of print].

48. Horak P, Crawford AR, Vadysirisack DD, Nash ZM, DeYoung MP, Sgroi D, et al. Negative feedback control of HIF-1 through REDD1-regulated ROS suppresses tumorigenesis. Proc Natl Acad Sci U S A. 2010;107:4675-80.

49. McBryan J, Howlin J, Napoletano S, Martin F. Amphiregulin: role in mammary gland development and breast cancer. J Mammary Gland Biol Neoplasia. 2008;13:159-69.

50. Stiehl DP, Bordoli MR, Abreu-Rodríguez I, Wollenick K, Schraml P, Gradin K, et al. Non-canonical HIF-2α function drives autonomous breast cancer cell growth via an AREG-EGFR/ErbB4 autocrine loop. Oncogene. 2012; 31:2283-97.

51. Müller A, Homey B, Soto H, Ge N, Catron D, Buchanan ME, et al. Involvement of chemokine receptors in breast cancer metastasis. Nature. 2001;410:50-6.

52. Sun X, Cheng G, Hao M, Zheng J, Zhou X, Zhang J, et al. CXCL12 / CXCR4 / CXCR7 chemokine axis and cancer progression. Cancer Metastasis Rev. 2010; 29:709-22.

53. Kang H, Watkins G, Parr C, Douglas-Jones A, Mansel RE, Jiang WG. Stromal cell derived factor-1: its influence on invasiveness and migration of breast cancer cells in vitro, and its association with prognosis and survival in human breast cancer. Breast Cancer Res. 2005;7:R402-10.

54. Burger JA. CXCR4: a key receptor in the crosstalk between tumor cells and their microenvironment. Blood. 2006;107:1761-7.

55. Wendt MK, Cooper AN, Dwinell MB. Epigenetic silencing of CXCL12 increases the metastatic potential of mammary carcinoma cells. Oncogene. 2008;27:1461-71.

56. Carriere P, Llopis S, Naiki A, Nguyen G, Phan T, Nguyen M, et al. Glyceollin I reverses epithelial to Mesenchymal transition in Letrozole resistant breast cancer through ZEB1. Int J Environ Res Public Health. 2015;13:10.

57. Xie Q, Bai Q, Zou L-Y, Zhang Q-Y, Zhou Y, Chang H, et al. Genistein inhibits DNA methylation and increases expression of tumor suppressor genes in human breast cancer cells. Genes Chromosomes Cancer. 2014;53:422-31.

58. Carroll JS, Meyer CA, Song J, Li W, Geistlinger TR, Eeckhoute J, et al. Genome-wide analysis of estrogen receptor binding sites. Nat Genet. 2006; 38:1289-97.

Submit your next manuscript to BioMed Central and we will help you at every step:

* We accept pre-submission inquiries

* Our selector tool helps you to find the most relevant journal

* We provide round the clock customer support

* Convenient online submission

* Thorough peer review

* Inclusion in PubMed and all major indexing services

* Maximum visibility for your research

Submit your manuscript at
www.biomedcentral.com/submit

BioMed Central