Neonatal Genistein Exposure and Glucocorticoid Signaling in the Adult Mouse Uterus

Background: Female reproductive tract development is sensitive to the endocrine-disrupting potential of environmental estrogens. Early-life exposure to the dietary phytoestrogen genistein impairs fertility and persistently alters the transcriptome in the oviduct and uterus of rodents. Glucocorticoid signaling, which has recently been shown to be essential for normal fertility in the female mouse uterus, is antagonized by genistein. Objective: Our goal was to determine whether early-life exposure to genistein disrupts glucocorticoid signaling in the mouse uterus, which may contribute to infertility. Methods: Female C57Bl/6 mice were exposed to either 50mg/kg per day genistein, 10μg/kg per day estradiol, or vehicle (corn oil) on postnatal days 1–5 (PND1–5), and then treated with the synthetic glucocorticoid dexamethasone (Dex: 1mg/kg) or vehicle (saline) on PND5, at weaning on PND21, or as adults on PND56 following adrenalectomy and ovariectomy to evaluate glucocorticoid responsiveness. Uteri were isolated following treatment for gene expression or chromatin immunoprecipitation. Results: Neonatal exposure to genistein altered the uterine transcriptome of adult mice and caused substantial changes to the transcriptional response to glucocorticoids. Although expression of the glucocorticoid receptor was not affected, genistein exposure disrupted glucocorticoid receptor recruitment to specific regulatory sites in target genes. Many genes involved in chromatin remodeling were dysregulated in genistein-exposed mice, suggesting that epigenetic reprograming may contribute to the altered glucocorticoid response of the uterus following early-life exposure to genistein. These changes affected the biological activity of glucocorticoids within the uterus, as glucocorticoids antagonized the proliferative effects of estradiol in the uterus of control mice but not genistein-exposed mice. Conclusions: Our findings suggest that disruption of glucocorticoid signaling due to early-life exposure to environmental estrogens may in part render the uterus unable to support implantation. https://doi.org/10.1289/EHP1575


Introduction
Environmentally derived compounds with estrogenic structures are recognized endocrine disruptors. The female reproductive tract is particularly sensitive to the effects of such compounds, so much so that toxicologists use the uterotrophic assay to screen for health risk (O'Connor et al. 1996). Exposure to these environmental estrogens, which are present in household and cosmetic products, pesticides and herbicides, food additives, groundwater, plastics, and plants, can impair reproductive function in a number of species. For this reason, the effects of soy consumption on human health have increasingly been the subject of much debate. Soy contains high levels of isoflavones, a class of phytoestrogens that can mimic endogenous estradiol (E 2 ) activity by binding to estrogen receptors (ERs) (Choi et al. 2008). The endocrine-disrupting properties of these compounds present a potential threat to fertility and reproduction in mammals (Caserta et al. 2008).
Although isoflavones have well-described health benefits in cancer (Mohamed et al. 2017;Spagnuolo et al. 2015), the mechanisms underlying these benefits also result in adverse effects on the proliferative nature of the estrogen-sensitive endometrium (Plaza-Parrochia et al. 2017). Genistein is the most abundant of the soybean isoflavones, accounting for approximately 50% of the total soybean isoflavone content (Murphy et al. 2002). Reproductive disturbances have been reported in a number of species fed soy as a significant portion of their diet, including rats, mice, rabbits, sheep, cattle, and cheetahs (Bennetts et al. 1946;Carter et al. 1955;Kendall et al. 1950;Setchell et al. 1987;Thain 1966). A randomized study described an increased incidence of endometrial hyperplasia in women receiving soy supplements long term (Unfer et al. 2004). Serum genistein levels in women consuming a nonvegetarian diet fall within a range of 2:6-22:6 nM, whereas levels are reported to be between 148 and 360 nM in vegetarians and likely higher in those consuming soy supplements (Elorinne et al. 2016;Peeters et al. 2007).
Early-life exposures to exogenous compounds that mimic the activity of endogenous hormones have the potential to permanently alter developing organs and tissues. Therefore, developmental exposure to genistein is of particular concern given that about 12% of formula-fed infants in the United States are fed soy-based formula during their first year of life (Rossen et al. 2016). Serum genistein levels in these infants occur in the range of 1-5 lM, which is several-fold higher than serum levels experienced in adults (vegetarian or nonvegetarian diet) and the dose reported to compete with E 2 for estrogen receptor binding (Cao et al. 2009;Rossen et al. 2016;Wang et al. 1996). The reported serum concentrations in infants fed soy formula also overlap with the concentration range shown in rodents to produce persistent adverse reproductive effects (approximately 3-7 lM serum genistein) (Doerge et al. 2002). In rodents, neonatal genistein exposure results in significant disruptions to the structure and function of the female reproductive tract that manifest in adults Newbold et al. 2001). Adult female rodents exposed to genistein as neonates exhibit sub-to complete infertility, resulting from altered estrous cyclicity, disrupted development of the oviduct, and an insufficient uterine environment (Awoniyi et al. 1998;Carter et al. 1955;Jefferson et al. 2009Jefferson et al. , 2012Nagao et al. 2001). Global gene analysis of the adult female oviduct following neonatal genistein exposure revealed substantial changes to basal gene expression, as well as the transcriptional response to pregnancy (Jefferson et al. , 2012. Interestingly, marked changes in immune response genes were reported following neonatal genistein exposure, some specific to the onset of pregnancy. Neonatal exposure to diethylstilbestrol (DES), a potent synthetic estrogen, also causes changes to female reproductive tract gene expression (Newbold et al. 2007). Neonatal DES exposure temporarily altered the expression of many chromatin-modifying proteins and persistently altered epigenetic marks in specific genomic regulatory regions in the adult mouse uterus (Jefferson et al. 2013). Soy-formula-fed girls demonstrate site-specific differences in DNA methylation in vaginal epithelial cells compared with girls fed cow formula (Harlid et al. 2017). These studies suggest that the mechanism of infertility following neonatal estrogen exposure during the developmental period is altered transcriptional activity in the adult female reproductive tract.
Classically, glucocorticoids are considered the primary mediators of pro-and anti-inflammatory actions within the immune response (Busillo and Cidlowski 2013). Glucocorticoids exert their effects through the glucocorticoid receptor (GR), a member of the nuclear receptor superfamily of transcription factors (Baxter and Tomkins 1970). Elevated levels of glucocorticoids due to stress or exogenous administration have deleterious effects on the reproductive tract, impacting the hypothalamus, pituitary, and gonads (Whirledge and Cidlowski 2017). Glucocorticoids regulate uterine biology through antagonizing the biological effects of E 2 , as well as independently coordinating uterine gene expression (Rhen et al. 2003). During pregnancy, glucocorticoids are thought to regulate many key processes, including the maternal immune response to pregnancy (Whirledge and Cidlowski 2017). In vitro studies have shown that glucocorticoids regulate thousands of genes within human uterine cell types, including genes involved in the inflammatory response that are also commonly regulated by E 2 (Whirledge et al. 2012(Whirledge et al. , 2013. In a mouse model, disruption of uterine GR signaling results in a profound subfertile phenotype, characterized by reduced blastocyst implantation, defects in stromal cell decidualization, dysregulation of immune response genes, and altered immune cell recruitment during early pregnancy (SD . When coadministered with glucocorticoids, genistein disrupts GR-dependent transcription of target genes in human uterine endometrial cells in vitro (S . These data, in addition to the overlap in reproductive and immune-modulatory functions of genistein and glucocorticoids, suggest that genisteininduced infertility may in part be driven by antagonism of glucocorticoid signaling. To better understand the contribution of glucocorticoid signaling to genistein-induced infertility, we investigated the long-term effects of early-life genistein exposure on glucocorticoid action in the uterus by using a mouse model of subcutaneous neonatal genistein exposure. Glucocorticoid-induced gene expression in the uterus was evaluated in the early postnatal period, prior to puberty, and in adult female mice that were adrenalectomized and ovariectomized postpuberty. Dexamethasone (Dex;1,17,dione; ≥98%) by thin layer chromatography (TLC) and E 2 (17-b-estradiol; ≥98% TLC) were purchased from Steraloids. Genistein (4 0 ,5,7-trihydroxyisoflavone; ≥98% HPLC) was purchased from Sigma Aldrich.

Animal Experiment Study Design
Timed pregnant C57Bl/6 mice were obtained from Charles River at 14.5 d following the presence of a vaginal plug, where the detection of the vaginal plug was considered day 0.5 of pregnancy. Pregnant mice were individually housed, maintained on a controlled 12:12 h light-dark cycle (lights on 0700-1900 hours), and provided ad libitum access to water and food. All mice were provided irradiated Teklad global soy protein-free extruded rodent diet (2920X). Isoflavone levels were measured in this diet to be 12:6 ppm, compared with commercial diets that contain approximately 830 ppm (Naaz et al. 2003). Litters were randomly assigned to experimental groups, and pups were injected subcutaneously on PND1-5 with either vehicle (corn oil) or genistein (50 mg=kg per day dissolved in corn oil). The dose of genistein chosen for animal injections produces circulating serum concentrations in mice that correspond to the serum genistein concentration of infants fed soy-based infant formulas (1-5 lM) (Cao et al. 2009;Doerge et al. 2002). For female pups, 50 mg=kg injected subcutaneously results in a maximum serum genistein concentration of 6:8 ± 1:4 lM with a halflife of 19 h (Doerge et al. 2002). The mouse and human estrogen receptor display relatively similar binding affinities for genistein, 8:1 ± 0:4 × 10 −7 M and 6:3 ± 0:7 × 10 −7 M respectively, and genistein has a reported half maximal inhibitory concentration (IC 50 ) of 5 × 10 −7 M for the estrogen receptor (Matthews et al. 2000;Wang et al. 1996). On the last day of control or genistein exposure (PND5), a subset of female mice pooled from several litters was randomly assigned to receive either vehicle (saline) or Dex (1 mg=kg dissolved in saline) by subcutaneous injection. The dose of Dex chosen for experiments has been validated to produce a robust transcriptional response in the mouse uterus and falls within the range of doses employed clinically in obstetrics (Boomsma et al. 2012;SD Whirledge et al. 2015). Tissues were harvested after 4 h, immediately snap-frozen with liquid nitrogen, and stored at −80 C until RNA extraction. The remaining pups were weaned on PND21 and litters from the same treatment groups were collectively housed at five females per cage. At weaning, a subset of female mice from pooled litters was randomly assigned to receive vehicle (saline) or Dex (1 mg=kg dissolved in saline) by intraperitoneal injection, and tissues were harvested after 4 h to isolate mRNA. Tissue collected for RNA extraction was snap-frozen with liquid nitrogen and stored at −80 C until RNA extraction was performed.
In previous studies employing the postnatal genistein exposure regimen, adult female mice were found to lack a regular estrous cycle . Moreover, variations in endogenous glucocorticoids among mice contribute to heterogeneity in basal gene expression (see Figure S1). In order to evaluate the uterine response to glucocorticoids in adult females without confounding effects of hormonal fluctuations from irregular cycles and endogenous glucocorticoids, the remaining mice underwent bilateral adrenalectomy (ADX) and ovariectomy (OVX) on PND42 to remove endogenous hormones. ADX/OVX animals were maintained on the same diet and provided 0.9% saline water during the 14-d recovery period.
To determine whether neonatal exposure to other estrogens could impact the uterine response to glucocorticoids, a subset of pups were exposed to 10 lg=kg E 2 dissolved in corn oil by subcutaneous injection on PND1-5. The dose of E 2 was selected to produce in vivo concentrations that model serum concentrations of E 2 during the estrous cycle and a robust uterotrophic response (Modder et al. 2004;Walmer et al. 1992) (see Figure S2). Adult mice underwent OVX/ADX on PND42 as previously described and were maintained on the same diet and provided 0.9% saline water during a 14-d recovery period.
On PND56, female mice from all treatment groups received either vehicle (saline), Dex (1 mg=kg dissolved in saline), E 2 (10 lg=kg dissolved in ethanol and diluted in saline), or Dex + E 2 by intraperitoneal injection. Uteri were collected 1.5 h later for chromatin immunoprecipitation and immunofluorescence, 4 h later for gene expression analysis, or 24 h later for cell proliferation. Tissue collected for RNA extraction was snap-frozen with liquid nitrogen and stored at −80 C until RNA extraction was performed. All animal procedures complied with the National Institutes of Health directives for the care and use of laboratory animals. All animals were treated humanely with regard to the alleviation of suffering. These studies were approved by the Institutional Animal Care and Use committees at the National Institute of Environmental Health Sciences.

Quantitative Real-Time PCR
Total RNA was extracted from the whole mouse uterus using the Qiagen RNeasy mini kit (Qiagen, Valencia, CA) according to the manufacturer's protocol with the deoxyribonuclease (DNase) treatment performed on the column. RNA purity and yield was assessed be evaluating the A260/A280 and concentration using the NanoDrop One Spectrophotometer (ThermoFisher Scientific, Waltham, MA). One hundred nanograms total RNA was used to synthesize cDNA using the One-Step RT-PCR Universal Master Mix reagent (ThermoFisher Scientific). Quantitative real-time polymerase chain reaction (qRT-PCR) was performed with the CFX Connect or CFX384 thermocycler (BioRad) using predesigned primer-probe sets (see Table S1) (ThermoFisher Scientific) in a 10-lL reaction volume. The thermocycling parameters for each reaction were 48 C for 30 min, 95 C for 10 min, followed by 40 cycles of 95 C for 15 s and 60 C for 60 s. Each gene primer-probe set was evaluated in technical duplicates from a standard curve and normalized to that of the reference gene peptidylprolyl isomerase B (PPIB), which was unaffected by treatment. Data are reported as fold change, and analyzed by ANOVA or Student's t-test using three to five biological replicates per treatment group.

Microarray Analysis
Gene expression analysis was conducted with three independent biological replicates using Agilent Whole Mouse Genome 4x44 multiplex format oligo arrays (product number G4122F; design ID 014868; Agilent Technologies) following the Agilent one-color microarray-based gene expression analysis protocol. The array is formatted to contain 43,379 gene probes. Starting with 400 ng of total RNA, cyanine 3 (Cy3)-labeled cRNA was produced according to the manufacturer's protocol. For each sample, 1:65 lg of Cy3-labeled cRNAs were fragmented and hybridized for 17 h in a rotating hybridization oven. Slides were washed and then scanned with an Agilent Scanner. Data were obtained using the Agilent Feature Extraction software (version 9.5), using the onecolor defaults for all parameters. The Agilent Feature Extraction Software performed error modeling, adjusting for additive and multiplicative noise. Preliminary analyses were performed with OmicSoft Array Studio software (version 7.0). In accordance with Minimum Information About a Microarray Experiment (MIAME) guidelines, the raw microarray data will be available in the Gene Expression Omnibus repository at the National Center for Biotechnology Information (https://www.ncbi.nlm. nih.gov/geo/query/acc.cgi?acc=GSE107788).
Differential gene expression was analyzed using the Partek Genomics Suite (version 6.6; (Partek Inc.). Significant changes in gene expression were defined based on an analysis of variance (ANOVA) p-value of <0:01 and a false discovery rate of p < 0:05. The lists of significantly regulated probes were visually sorted by Venn diagram (http://bioinformatics.lu/venn.php). Heat maps were generated in Partek Genomics Suite and provided hierarchical clustering dendrograms. The statistically significant probes were also analyzed with Ingenuity Pathway Analysis (IPA; Ingenuity Systems) to evaluate their functional annotation. Gene set enrichment was determined by IPA using Fisher's exact test (p < 0:05). Significantly regulated genes were also analyzed by gene ontology using Gene Annotation Tool to Help Explain Relationships (GATHER) (Chang and Nevins 2006). Gene ontology GATHER analysis used a Bayes factor cutoff of 5.0.

Corticosterone Assay
Blood samples from experimental animals (11-12 animals per experimental group) were obtained by submandibular vein bleed between 0930 and 1100 hours. Submandibular vein bleed was performed in a procedure room, separate from where the experimental mice were housed. Care was taken to minimize the amount of time each mouse was handled, and mice were placed in a new cage separate from cage-mates immediately following the procedure. At the completion of the procedure, mice were returned to their original cage. Corticosterone was measured by the commercially available DetectX Corticosterone Enzyme Immunoassay Kit (Arbor Assays). All samples were run in duplicate. The limit of detection for the assay was determined to be 16:9 pg=mL. The sensitivity was determined to be 18:6 pg=mL. The intraassay precision was determined to be 6.5-8.1%. The interassay precision was determined to be 9.9-16.4%.

Immunofluorescence
Cryopreserved frozen uteri from three mice per treatment group were sectioned to 5-lm serial sections and stored at −20 C. Prior to immunostaining, slides were thawed at room temperature. Three sections per uterus were stained for histological analysis with antibodies to GR (1:500; Cell Signaling Technologies) or MKi67 (1:400; Cell Signaling Technologies) overnight at 4°C (see Table  S2). Sections were then washed with phosphate buffered saline (PBS) and incubated with secondary antibodies (see Table S2). Hoescht 33342 (1:5,000; Molecular Probes) stain was applied prior to mounting to visualize the nucleus. Images were obtained on a Zeiss LSM780 confocal microscope equipped with a 63 × (oil) objective and processed using the Zen 2012 software. Epithelial cell proliferation was scored for each image taken by counting the number of Ki67-positive epithelial cells and dividing by the total number of epithelial cells in the image. For each treatment group, three to four images were counted.

Western Blotting
Tissues from control (n = 5) and genistein-treated (n = 4) mice were lysed in Tris glycine SDS sample buffer supplemented with 2-mercaptoethanol (BME). Equal amounts of protein from tissue extracts were separated via SDS-PAGE and transferred to a nitrocellulose membrane. Membranes were blocked with 7.5% skim milk in Tris-buffered saline (TBS) and incubated overnight with mouse monoclonal anti-b-actin antibodies (1:10,000; Millipore) and rabbit monoclonal anti-GR antibodies (1:1,000; Cell Signaling Technologies) in 5% milk in TBS-Tween (0.1%) (see Table S2). Membranes were washed and incubated with goat anti-rabbit IRDye 680-conjugated secondary antibody and goat antimouse IRDye 800conjugated secondary antibody (LI-COR Biosciences) for 1 h at room temperature. The Odyssey LI-COR imaging system (LI-COR Biosciences) was used to visualize protein expression. GR protein levels were normalized to b-actin and expressed relative to control.

Chromatin Immunoprecipitation
Uteri were pooled together (five per sample) and pulverized by mortar and pestle in liquid nitrogen. Proteins were fixed to chromatin with 1% paraformaldehyde for 10 min at room temperature.

047002-3
Cells were then washed in ice-cold PBS and cross-linking stopped with 125 nM glycine. Cells were resuspended in cell lysis buffer containing 50 mM Hepes-KOH pH 8.0, 1 mM EDTA, 140 mM NaCl, 10% glycerol, 0.5% NP-40, 0.25% Triton X-100, and protease inhibitors and dounce homogenized. Samples were centrifuged for 10 min at 5,000 rpm at 4°C and resuspended in shearing buffer containing 10 mM Tris-HCl pH 8.0, 1 mM EDTA, 140 mM NaCl, 1.0% SDS, 0.1% Na deoxycholate, and 1% Triton X-100 with protease inhibitors. Samples were then sonicated using the Fisher Scientific Model 120 Sonic Dismembrator (ThermoFisher Scientific) at 35% for 6 min to obtain 200-to 500-bp fragments as confirmed by agarose gel electrophoresis. Sheared chromatin was precleared with protein A agarose/salmon sperm DNA (Millipore) and immunoprecipitated overnight with monoclonal GR antibody (1:175; Cell Signaling Technologies) or normal rabbit IgG (1:700; Millipore) (see Table S2). The protein-DNA complexes were then precipitated with protein A magnetic beads (BioRad). After washing and elution, cross-links were reversed and DNA was purified using the QiaQuick PCR Purification kit (Qiagen). Immunoprecipitated DNA was quantified from four to five independent pooled samples per treatment group using qRT-PCR with custom-designed primer-probe sets ordered from Integrated DNA Technologies (see Table S3). Relative expression values for each primer were calculated using the DDCt method and set relative to control vehicle IgG values. Previously reported glucocorticoid response elements (GREs) in glucocorticoidinduced leucine zipper (Gilz) and Krüppel-like factor 13 (Klf13) were used to design primers (Cruz-Topete et al. 2016). Regions of GR binding in the mouse serum/glucocorticoid regulated kinase 1 (Sgk1) and period circadian clock 1 (Per1) genes have been previously described (Yu et al. 2010

Statistical Analysis
The data are presented as mean ± SEM of a minimum of three biological replicates. Statistical significance was determined by ANOVA with Tukey's post hoc analysis (control vs. treatment group exposed to vehicle or Dex) or Student's t-test (control vs. genistein) using StatTools software. Statistical significance is reported as p < 0:05 (*) or p < 0:01 (**).

Effect of Neonatal Genistein Exposure on Glucocorticoid-Regulated Gene Expression in the Uterus
To assess the impact of neonatal genistein exposure on the uterine response to glucocorticoids, female mice were injected subcutaneously on PND1-5 with 50 mg=kg genistein or corn oil control ( Figure 1A). At the completion of the treatment regimen on PND5, a subset of the control and genistein-exposed female mice were treated with the synthetic glucocorticoid Dex (1 mg=kg) or vehicle (saline) by subcutaneous injection and uteri were collected 4 h post injection. Body and uterine weights were measured and mRNA was isolated from the uterus. Developmental exposure to genistein did not alter body weight, but the estrogenic activity of genistein transiently increased the uterine/body weight ratio at PND5 ( Figure 1A). Expression of Gilz, FK506 binding protein 5 (Fkbp5), Sgk1, Klf13, and Per1, and interleukin 13 receptor subunit alpha 2 (Il13ra2) and 11b-hydroxysteroid dehydrogenase II (11bhsdII), classic target genes of GR signaling, was quantified by qRT-PCR.
On PND5, genistein exposure did not alter basal mRNA expression of Gilz, Fkbp5, Klf13, or Per1, 11bhsdII, or the reference gene peptidylprolyl isomerase B (Ppib) ( Figure 1B), although basal expression of Sgk1 was significantly higher and basal expression of Il13ra2 was significantly lower in the genistein-exposed mice. Dex treatment induced all of these genes equivalently in both treatment groups, suggesting that genistein-exposed mice had normal glucocorticoid responsiveness at that time ( Figure 1C). Basal expression and glucocorticoid-mediated induction of Ppib, Gilz, Fkbp5, Sgk1, Klf13, Per1, Il13ra2, and 11bhsdII were also evaluated at weaning (PND21) in a subset of the control and genisteinexposed female mice treated for 4 h with Dex (1 mg=kg) or vehicle (saline) (Figure 2A). On PND21, prior genistein exposure did not affect basal expression of Ppib, Gilz, Fkbp5, Sgk1, Klf13, Per1, Il13ra2, or 11bhsdII mRNA ( Figure 2B). Furthermore, the induction of Gilz, Fkbp5, Sgk1, Klf13, Per1, and 11bhsdII by Dex treatment was indistinguishable in the genistein-exposed and control uteri, suggesting that the uterine response to glucocorticoids was not altered prior to puberty ( Figure 2C). Long-term effects of the neonatal genistein exposure were evaluated in adult female mice that underwent the postnatal genistein dosing regimen ( Figure 3A). Adult control and genistein-exposed female mice were adrenalectomized/ovariectomized (ADX/OVX) on PND42 to remove endogenous hormones. Following a 2-wk recovery period, ADX/OVX mice were treated with 1 mg=kg Dex or vehicle (saline) by intraperitoneal injection and uteri were harvested after 4 h. At this time point, body weight and uterine/body weight ratio were equivalent between treatment groups ( Figure  3A). Basal mRNA expression of Ppib, Gilz, Fkbp5, Sgk1, Klf13, Per1, IL13ra2, and 11bhsdII was not different in adult mice exposed to genistein on PND1-5 compared with control mice ( Figure 3B). However, induction of the glucocorticoid-responsive genes Gilz, Fkbp5, Sgk1, and Klf13 was significantly blunted in adult genistein-exposed mice compared with control mice ( Figure  3C). Glucocorticoid-mediated induction of Per1, 11bhsdII, and Il13ra2 was not affected by early-life genistein exposure ( Figure  3D). These data indicate that glucocorticoid signaling is selectively altered in the mature adult uterus but not the neonatal uterus by early-life exposure to genistein (see Table S4).
Neonatal genistein exposure caused lasting effects to basal gene expression in the oviduct . In order to determine the long-term effects on basal gene expression in the uterus, whole-genome microarray analysis was performed on uteri from three biological replicates of adult ADX/OVX mice previously exposed to corn-oil control or genistein as neonates (Figure 4). A comparison of the significantly regulated genes identified 943 genes that demonstrated different basal expression levels in the uterus of adult mice exposed to genistein as neonates compared with controls ( Figure 4A). Of the significantly regulated genes, 350 genes were induced, whereas 593 genes were repressed. Shown in Table S5 are the top 10 induced and repressed genes determined by fold change. Interestingly, IPA performed on the genes with significantly different basal expression identified many molecular and biological pathways related to regulation of cell cycle and cancer ( Figure 4B). Among IPA diseases and biological functions, cancer was the top ranked and cell cycle the fifth ranked annotation most significantly associated with the dysregulated genes ( Figure 4C).

Genome-Wide Gene Expression Analysis
Attenuated induction of Gilz, Fkbp5, Sgk1, and Klf13 by glucocorticoids in the uterus of adult mice exposed to genistein during development suggested that the transcriptional response to glucocorticoids was altered, thus whole-genome microarray analysis was performed on uteri from adult control and genistein-exposed mice treated with vehicle or 1 mg=kg Dex for 4 h. Three biological replicates were evaluated for each treatment group. A heat map of sample replicates and hierarchical clustering analysis of glucocorticoid-responsive genes indicated samples separated into distinct groups based on exposure and treatment ( Figure 5A). Glucocorticoids significantly altered the expression of 3,832 genes (1,903 genes induced and 1,929 genes repressed) in the uteri of control mice compared with 2,353 genes (1,335 genes induced and 1,018 genes repressed) in the uteri of genistein-exposed mice ( Figure 5B). A comparison of the genes significantly regulated by glucocorticoids in both the control and genistein-exposed groups revealed only 1,033 genes in common, including 77 genes that were differentially regulated by glucocorticoids (see Table S6). Neonatal genistein exposure resulted in the loss of glucocorticoid regulation of 2,799 genes and the unique regulation of 1,320 genes by glucocorticoids. To identify the molecular and biological pathways associated with the glucocorticoid-regulated uterine transcriptome in control and genistein-exposed mice, the annotation of significantly regulated genes was evaluated in IPA software. IPA identified GR signaling as the top canonical pathway regulated by glucocorticoid treatment in both control and genistein-exposed mice ( Figure 5C). Evaluation of this pathway revealed significantly fewer genes regulated in genistein-exposed mice (see Figure S3). All of the top canonical pathways regulated in genistein-exposed mice had a smaller percent of significantly regulated genes than the top five canonical pathways regulated in control mice, reflecting the overall suppression of glucocorticoid signaling following neonatal genistein exposure. For example, the top canonical pathway regulated by glucocorticoids in control mice was GR signaling, where 85 of 287 genes were regulated (29.6%). In genistein-exposed mice, only 46 of 287 genes were regulated by glucocorticoids, representing 16% of the total genes in the pathway. The top predicted diseases and biological functions by IPA analysis and a separate gene ontology analysis revealed similar annotated functions in the control-and genistein-treated mice, though the number of significantly regulated genes in each function was diminished in the genistein-treated mice (see Tables S7-S9).
The genes commonly regulated by glucocorticoids in the control and genistein-exposed mice were compared to determine if early-life exposure to genistein altered the extent or directionality of the uterine transcriptional response to glucocorticoids. Of the 1,033 genes, 13 genes were regulated in the opposite direction (anticorrelated) (see Table S10), and a change in the magnitude of induction or repression greater than 2-fold was found for 64 genes. A comparison analysis of the genes common to the control and genistein-exposed mice indicated the predicted activation of the top pathways and functions was analogous between treatment groups, indicating that the differential regulation of common genes did not affect activation of predicted functions (see Figure S4). Interestingly, a comparison analysis of the top canonical pathways and diseases and biological functions in the uniquely regulated genes (2,799 unique genes in control mice; 1,320 unique genes in genistein-exposed mice) showed numerous differences in the predicted activation, including several pathways/functions that were predicted to have opposing activation scores ( Figure 5D). These data indicate that disruption of glucocorticoid signaling by neonatal exposure to genistein leads to the dysregulation of many genes and pathways with potentially important roles in uterine biology.
In order to assess whether the altered transcriptional response to glucocorticoids in adult mice following neonatal genistein exposure was specific to genistein or broadly represents early-life exposure to estrogens, subcutaneous injection of 10 lg=kg per day E 2 dissolved in corn oil was also included as a treatment arm in the PND1-5 injection regimen. Adult ADX/OVX female mice exposed to control, genistein, or E 2 treatment on PND1-5 were treated with 1 mg=kg Dex or vehicle (saline) for 4 h on PND56. The differential expression of select genes identified by the whole-genome microarray analysis was validated in an independent cohort of genistein-exposed mice and compared with mice exposed to E 2 ( Figure 6). Genes were chosen for analysis based on their differential patterns of glucocorticoidresponsiveness determined by whole-genome microarray analysis. Specifically, genes with changes in the magnitude of induction or Whole genome microarray analysis was performed on total RNA isolated from the uterus of adult ADX/OVX PND56 mice previously exposed to control or genistein. Three biological replicates were evaluated for each treatment group. The number of probes that were significantly different (p < 0:01) between the control and genistein-exposed group were graphed as induced (red) or repressed (green). (B) Significantly different genes were evaluated in IPA software and the top 10 regulated canonical pathways are reported. The total number of genes in the IPA-annotated pathway is listed on the right side of the graph. The total percentage of significantly regulated genes is depicted as the sum of induced and repressed genes. (C) The top 10 diseases and biological functions are listed. The p-value range and number of genes within the listed functions is provided. ADX/OVX, adrenalectomized/ovariectomized; IPA, Ingenuity Pathway Analysis; PND, postnatal day. exhibiting a novel regulation pattern following neonatal genistein exposure were evaluated. Expression of a disintegrin-like and metallopeptidase with thrombospondin type 1 motif (ADAMTS)-like 2 (Adamtsl2) mRNA was significantly increased in response to glucocorticoid treatment in control mice. The magnitude of induction by glucocorticoids was blunted in genistein-exposed mice and further decreased in E 2 -exposed mice. Induction of Indian hedgehog (Ihh) mRNA by glucocorticoids was validated in the uterus of control and genistein-exposed mice. However, upregulation of Ihh was absent following Dex treatment in E 2 -exposed mice. Microarray results indicated that early-life exposure to genistein impeded glucocorticoid-mediated induction Figure 5. Effect of genistein exposure on genome-wide regulation by glucocorticoids. (A) Adult ADX/OVX (PND56) mice, exposed as neonates to control or genistein, were treated with vehicle (saline) or 1 mg=kg Dex for 4 h. Whole genome microarray analysis was performed on total RNA isolated from the uterus, with three biological replicates evaluated for each treatment group. Heat map and hierarchical clustering of significantly regulated genes (p < 0:01), where red represents induced genes and green represents repressed genes. (B) The number of significantly regulated genes by treatment (p < 0:01) were compared between the control and genistein groups by Venn diagram. (C) Significantly regulated genes were annotated with IPA software. The top five regulated canonical pathways are reported for the control and genistein groups. The total number of genes in the IPA-annotated pathway is listed on the right side of the graph. The total percentage of significantly regulated genes is depicted as the sum of induced and repressed genes. (D) Genes uniquely regulated by Dex in the control (2,799) and genistein (1,320) groups. The top five canonical pathways and diseases and biological functions were ranked and assigned an activation z-score (orange indicates induced; blue indicates repressed). ADX/OVX, adrenalectomized/ovariectomized; Con, control; Dex, dexamethasone; IPA, Ingenuity Pathway Analysis; Neo, neonatal; Veh, vehicle.
of Resistin (Retn), which was validated in independent samples. Interestingly, early-life exposure to E 2 did not prevent the induction of Retn by glucocorticoids. Oviductal glycoprotein 1 (Ovgp1) was uniquely regulated by glucocorticoids following neonatal genistein exposure. E 2 exposure produced a comparable effect. These results suggest that the ability of estrogens to specifically alter glucocorticoid signaling in the uterus may be applicable to many other estrogenic compounds.

Effect of Neonatal Genistein Exposure on GR Recruitment to Target Genes
To define the mechanism by which neonatal genistein exposure alters uterine glucocorticoid signaling, we evaluated circulating levels of glucocorticoids, expression of GR, and nuclear translocation of GR in response to ligand. There were no differences in serum levels of corticosteroids in adult control and genisteinexposed mice ( Figure 7A). The expression of uterine GR mRNA (Nr3c1) and protein was similar in control and genistein-treated mice, indicating that differences in glucocorticoid responsiveness were not due to insufficient ligand or receptor expression ( Figure  7B). The ability of GR to undergo translocation to the nucleus following ligand binding was assessed by immunostaining and confocal microscopy. ADX/OVX adult control and genisteintreated mice were treated with vehicle (saline) or 1 mg=kg Dex for 1.5 h. Glucocorticoid treatment led to nuclear translocation of GR, which was otherwise localized to the cytoplasm, in luminal epithelial and stromal cells of both control and genistein-treated mice, suggesting that altered gene expression in the genisteinexposed mice was not a consequence of deficient nuclear translocation ( Figure 7C). These data demonstrate the presence of an intact glucocorticoid signaling system in the uterus of genisteinexposed mice.
The genomic actions of GR are mediated through a physical interaction with DNA response elements or through associations with chromatin-bound proteins. Chromatin immunoprecipitation (ChIP) assays were performed on the uteri of control and genisteintreated ADX/OVX adults 1.5 h following vehicle or Dex treatment to determine whether altered glucocorticoid signaling following neonatal genistein exposure resulted from dysregulated GR recruitment to Figure 6. Effect of neonatal genistein exposure on glucocorticoid-responsive genes differentially regulated in genistein-exposed mice. Adult ADX/OVX (PND56) mice, exposed as neonates to control, 50 mg=kg genistein or 10 lg=kg E 2 on PND1-5, were treated with vehicle (saline) or 1 mg=kg Dex for 4 h. Relative mRNA expression in the uterus of Adamtsl2, Ihh, Retn, and Ovgp1 were measured by qRT-PCR and values were normalized to the reference gene Ppib. Dex treatment for each neonatal exposure group (vehicle, genistein, or E 2 ) was set relative to vehicle treatment for that group (i.e., E 2 -exposed treated with Dex set relative to E 2 -exposed vehicle). The reported fold change from microarray analysis is listed below the qRT-PCR results. The results represent the mean of 4-5 animals ± SEM. ** p < 0:01 according to ANOVA with Turkey's post hoc analysis. ADX/OVX, adrenalectomized/ovariectomized; Dex, dexamethasone; E 2 , estradiol; Gen, genistein; N/D, not determined to be significantly regulated by microarray analysis; PND, postnatal day; qRT-PCR, quantitative real-time polymerase chain reaction; Veh, vehicle. regulatory regions of target genes. For this experiment, GR recruitment was evaluated for Gilz, Sgk1, and Klf13, which demonstrated blunted induction of mRNA expression by glucocorticoid treatment following neonatal genistein exposure, and Per1, which did not demonstrate differences in glucocorticoid-induced mRNA expression. A functional GRE has been reported for Gilz and Klf13 (Cruz-Topete Figure 7. Effect of neonatal genistein exposure on ligand and GR expression. (A) Serum levels of corticosterone were measured in adult mice with intact ovaries and adrenal glands from serum collected between 0930 and 1130 hours in the morning. Data represent the mean of 11-12 animals ± SEM. (B) Relative mRNA expression of GR (Nr3c1) was measured by qRT-PCR in adult ADX/OVX control and genistein-exposed mice. Expression was normalized to the reference gene Ppib and set relative to expression in control mice. GR protein levels, quantified by western blot analysis, were normalized to levels of the reference protein b-actin, which was not altered by treatment group (see Figure S5). GR protein levels were set relative to control mice. The results and image represent the mean of 4-5 animals ± SEM. (C) Representative images of GR expression and localization (red) in luminal epithelial and endometrial stromal cells. Hoescht 33342 was used to visualize nuclei (blue). Images were taken at 630 × . ADX/OVX, adrenalectomized/ovariectomized; Dex, dexamethasone; GR, glucocorticoid receptor; qRT-PCR, quantitative real-time polymerase chain reaction; Veh, vehicle. et al. 2016), and in silico analysis of the Sgk1 and Per1 genes identified GREs in the introns of these genes with complete homology to the critical nucleotides of the consensus sequence. Similar to our findings of reduced induction of Gilz mRNA in response to glucocorticoids in mice exposed to genistein, recruitment of GR to a previously identified GRE in Gilz was abrogated in the uterus of genisteinexposed mice ( Figure 8A). Glucocorticoid receptors were strongly recruited to the Sgk1 intron-3 GRE in uteri from control mice exposed to glucocorticoids ( Figure 8A). However, Dex-induced recruitment of GR to the GRE of Sgk1 was weak in genistein-treated mice and did not obtain significance. Activated GRs were recruited to the Klf13 intron-1 GRE in both control and genistein-treated mice. However, recruitment of GR was significantly less in the genistein-treated mice compared with control mice. In contrast to our findings for Gilz, Sgk1, and Klf13, glucocorticoid-induced GR recruitment to the Per1 gene was not affected by neonatal genistein exposure ( Figure 8B). These results suggest that neonatal genistein exposure alters GR recruitment to regulatory regions in a gene-specific manner.
The recruitment of GR to accessible regions of chromatin is partially dependent on the maintenance of an open state by chromatinand DNA-modifying enzymes (Mundade et al. 2014). mRNA expression of genes involved in DNA and histone modifications was evaluated to explore the underlying mechanism by which neonatal genistein exposure results in changes to the profile of GR recruitment. Compared with control mice, basal mRNA levels of DNA (cytosine-5)-methyltransferases 1 (Dnmt1), 3A (Dnmt3a), and 3B (Dnmt3b) and histone deacetylases 1 (Hdac1) and 2 (Hdac2) were lower in adult ADX/OVX mice exposed to genistein as neonates ( Figure 9A). Relative mRNA expression of Methyl-CpGbinding domain protein 2 (Mbd2) and Methyl-CpG-binding protein 2 (Mecp2), transcriptional repressors that bind to methylated DNA, did not significantly change ( Figure 9B). These data suggest that neonatal genistein exposure causes long-term effects on the expression of specific genes involved in DNA and histone modifications.

Effect of Neonatal Genistein Exposure on the Functions of Glucocorticoids in the Adult Mouse Uterus
The ability of exogenous glucocorticoids to antagonize the uterotrophic effects of E 2 has been well described (Howe et al. 1990;Rabin et al. 1990;Rhen et al. 2003). To evaluate whether changes in glucocorticoid signaling resulting from neonatal genistein exposure disrupted the biological functions of glucocorticoids in the uterus, morphology was evaluated in control and genistein-exposed mice 24 h following vehicle (saline), 10 lg=kg E 2 , or 1 mg=kg Dex and 10 lg=kg E 2 injection in 3 mice per treatment group ( Figure 10A). In control and genistein-exposed mice, E 2 treatment induced the classic uterotrophic response, which is characterized by increased uterine weight through water imbibition and cell proliferation (see Figure S2). Glucocorticoid administration with E 2 blocked estrogeninduced uterine growth in control mice but was only able to partially restrict estrogen-induced uterine growth in genistein-exposed mice. Cell proliferation is partially responsible for increased uterine weight following E 2 exposure (Quarmby and Korach 1984). Uterine sections from control and genistein-exposed mice treated with vehicle, E 2 , or Dex + E 2 were stained for the proliferation marker Ki-67 to assess whether glucocorticoid antagonism of cell proliferation was altered ( Figure 10B). E 2 significantly increased uterine epithelial cell proliferation in control and genistein-exposed mice compared with vehicle-treated mice. In control mice, glucocorticoids administered with E 2 blocked the increase in proliferation. However, the number of Ki67 positive cells was not reduced in the genistein-exposed mice treated with Dex + E 2 . Because we observed significant differences Figure 8. Effect of neonatal genistein exposure on GR recruitment to specific GREs. Adult ADX/OVX PND56 mice, exposed neonatally to control or genistein, were treated with vehicle (saline) or 1 mg=kg Dex for 1.5 h. ChIP assays were performed with rabbit IgG or rabbit anti-GR antibody. Coimmunoprecipitated DNA was analyzed by qRT-PCR using primers to GREs and plotted relative to input DNA. (A) Relative fold enrichment of GR recruitment to Gilz, Sgk1, and Klf13 as determined by ChIP assays. (B) Relative fold enrichment of GR recruitment to Per1. The results represent the mean of 4-5 animals ± SEM. * p < 0:05, ** p < 0:01 as determined by Student's t-test. ADX/OVX, adrenalectomized/ovariectomized; ChIP, chromatin immunoprecipitation; Dex, dexamethasone; GR, glucocorticoid receptor; GRE, glucocorticoid response element; IgG, immunoglobulin G; PND, postnatal day; Veh, vehicle. in the regulation of cell proliferation, we evaluated genes significantly regulated by glucocorticoids in the cellular growth and proliferation pathway by IPA ( Figure 10C). In control mice, glucocorticoids regulated 175 genes related to cellular growth and proliferation. In genistein-exposed mice, glucocorticoid regulation of this pathway was limited, where only 64 of 175 genes were significantly regulated. Therefore, the failure of glucocorticoids to antagonize the uterotrophic effects of E 2 in genistein-exposed mice are likely the result of altered glucocorticoid signaling.

Discussion
Disruptions to the endocrine system by environmental compounds have varying consequences depending on the timing of exposure. Early-life exposures can result in persistent adverse effects due to the sensitive nature of the developing organs. For example, the long-term consequences of in utero exposure to the xenoestrogen DES are well-documented in humans (Reed and Fenton 2013). Human epidemiological data linking developmental exposure to phytoestrogens and adult pathologies remain limited, but animal models have clearly described lasting effects, some of which parallel the effects resulting from neonatal DES exposure (Adgent et al. 2012;D'Aloisio et al. 2012;Harlid et al. 2017;Nagao et al. 2001;Newbold et al. 2001;Padilla-Banks et al. 2006;Upson et al. 2015). In rodents, disruptions to the reproductive tract following neonatal genistein exposure correspond to alterations in the expression of genes key to the morphology and function of the reproductive tract . Global gene expression has been evaluated in the rodent oviduct following neonatal genistein exposure, but the assessment of the long-term consequences of neonatal genistein exposure in the rodent uterus has been limited to candidate genes (Begum et al. 2006;Jefferson et al. 2012;Tang et al. 2008). The results presented here demonstrate that early-life exposure to genistein at an environmentally relevant dose, which is comparable to exposure in infants consuming soy-based formula, causes substantial genome-wide changes in uterine gene expression and glucocorticoid responsiveness in adult mice (Jefferson and Williams 2011). , and Hdac2 in adult ADX/OVX mice exposed previously to control or genistein. (B) Relative mRNA expression of Mbd2 and Mecp2 in adult ADX/OVX mice exposed previously to control or genistein. Values were normalized to the reference gene Ppib and set relative to expression in control mice. **. p < 0:01 as determined by Student's t-test. ADX/OVX, adrenalectomized/ovariectomized. Neonatal genistein exposure, at concentrations equivalent in estrogenic activity to DES, resulted in uterine adenocarcinoma in adult mice (Newbold et al. 2001). Our microarray analysis indicated that neonatal genistein exposure altered the transcriptional profile of the adult uterus at baseline, and analysis of these genes suggested a potential mechanism for the observed uterine adenocarcinoma. IPA software ranked canonical pathways related to cell cycle regulation and cancer as the top pathways for the set of genes significantly altered by neonatal genistein exposure. In agreement with the identified canonical pathways, the top regulated disease was cancer. These data support genome-wide transcriptional reprograming as a possible mechanism leading to uterine adenocarcinoma.
Neonatal genistein exposure also renders the uterus unable to support implantation and pregnancy. Uterine glucocorticoid signaling plays an essential role in establishing endometrial receptivity and implantation (SD . Moreover, aberrant glucocorticoid signaling driven by stress-induced or exogenous elevations in glucocorticoid levels impairs reproductive success (Whirledge and Cidlowski 2017). In order to understand whether long-term effects on glucocorticoid signaling contribute to the deficient uterine environment, the glucocorticoid-regulated transcriptome was evaluated in the uterus of adult mice previously exposed to vehicle or genistein as neonates. The transcriptional response to glucocorticoids was persistently altered by neonatal exposure to genistein, and this differential pattern of regulation was validated in independent samples for the majority of tested genes. Gene ontology analysis indicated that the uniquely altered genes reflect preferentially activated and repressed canonical pathways and biological functions in the uterus. The pathways and functions activated by those genes uniquely regulated by glucocorticoids in the control group compared with those in the genistein-exposed group suggest long-lasting changes to uterine function may in part be mediated by selectively altered glucocorticoid signaling. For example, the cAMP signal transduction pathway mediates hormonal signals in the endometrium and plays a critical role in endometrial stromal cell decidualization (Brar et al. 1997;Telgmann et al. 1997). The activation status of the cAMP-mediated signaling pathway in response to glucocorticoids was predicted to be reversed in the uterus of adult mice exposed to genistein early in life. Furthermore, the glucocorticoid-mediated activation status of four of the top five diseases and biological functions were also predicted to be inversely correlated in the genisteinexposed mice compared with controls. These differentially responsive functions include cell viability, cell survival, organismal death, and differentiation of cells, which are essential functions critical to uterine remodeling in support of pregnancy. Similar functions were also found to be dysregulated in the subfertile uterine-GR deficient mouse in response to hormone priming for implantation, suggesting that disruption to glucocorticoid regulation of these functions could impair uterine competence (SD . The altered regulation of glucocorticoid target genes could be accounted for by differential recruitment of GR to regulatory sequences upon activation. Indeed, enrichment of GR at GREs within classic glucocorticoid target genes Gilz, Sgk1, and Klf13 was reduced in the uterus of adult female mice following neonatal genistein exposure, which corresponded to blunted mRNA induction of these genes in response to glucocorticoids. Deficient GR transactivation in this limited set of candidate genes may have implications for uterine biology and fertility. Sgk1 null mice have significantly fewer offspring related to spontaneous fetal loss, and lower endometrial SGK1 expression was reported in women with recurrent pregnancy loss (Salker et al. 2011). Although the actions of Gilz and Klf13 in the uterus are still under investigation, early studies suggest the expression of these genes may be important for endometrial function. Targeted reduction of Klf13 mRNA in human endometrial stromal cells decreased the expression of Bmp2, a morphogenic factor essential for stromal cell differentiation in early human pregnancy (Heard et al. 2012). Gilz is expressed in the mouse endometrium during pregnancy and levels distinctly change near parturition (Zhao et al. 2006). The biological significance of this is not clear, but GILZ is an important mediator of inflammation, which is a hallmark of parturition.
Modulation of transcription factor recruitment by genistein has been demonstrated in vitro (Chang et al. 2008;Chen et al. 2013;Li et al. 2008), although these findings are the first to show the persistent effects of genistein on recruitment in vivo. The mechanism leading to altered GR recruitment in a gene-specific manner subsequent to neonatal genistein exposure is likely reprograming of the epigenetic landscape during development, which has been documented for several xenoestrogens (Bromer et al. 2010;Jefferson et al. 2013;Strakovsky et al. 2014). Neonatal exposure to genistein in rodents resulted in site-specific hypo-and hypermethylation and associated changes in gene expression, which was also reported in vaginal epithelial cells from soy-formula-fed infants (Harlid et al. 2017;Strakovsky et al. 2014;Tang et al. 2008). Such site-specific changes to the epigenetic landscape may have contributed to the latent differences in glucocorticoid responsiveness of Gilz, Sgk1, and Klf13 compared with Per1, 11b-hsdII, and Il13ra2 in the uterus of adult mice following neonatal genistein exposure. Epigenetic changes have been associated with changes in the expression of chromatin-modifying enzymes (Greathouse et al. 2012). We found that mRNA expression of certain DNA methyltransferases and histone deacetylases was repressed in the uterus of adult genisteinexposed mice, possibly promoting global changes to basal gene expression and facilitating altered glucocorticoid responses via changes in the expression of specific co-regulators and/or enzymes that influence GR binding to target genes. In addition to changes to the chromatin environment, other mechanisms of transcriptional regulation, including posttranslational modifications to GR or cofactor availability, may also contribute to altered glucocorticoid signaling following early-life exposure to genistein.
It is notable that genistein can act as both an estrogen receptor ligand and inhibitor of receptor tyrosine kinases (Akiyama et al. 1987). However, neonatal genistein exposure with the estrogen receptor antagonist ICI 182780 blocked the transcriptional and functional consequences in the oviduct, suggesting that the effects of genistein on the female reproductive tract are mediated by its estrogenic activities (Jefferson et al. , 2012. In agreement with this, we found that early-life exposure to E 2 also altered glucocorticoidmediated gene expression in the uterus of adult mice. The expression of Ovgp1 was uniquely induced by Dex following postnatal treatment with genistein or E 2 . Early-life exposure to genistein or E 2 also diminished the magnitude of induction of Adamtsl2 by Dex. For some genes evaluated, the response to prenatal E 2 was unique from that of genistein, suggesting these molecules can have gene-specific effects. Prenatal exposure to the xenoestrogens bisphenol A and DES induce opposite effects on the methylation status of the promoter and intron 1 of homeobox A10, which suggests that xenoestrogens can exert distinct developmental effects potentially related to their estrogenic activity (Bromer et al. 2009(Bromer et al. , 2010. Moreover, these findings indicate that disrupted uterine glucocorticoid signaling may occur following exposure to a broad range of estrogenic chemicals and expand the implications for these findings.

Conclusion
The data presented herein demonstrate that early-life exposure to genistein disrupts both basal gene expression in the uterus and the transcriptional response to glucocorticoids in adult mice. We found that these differences occurred in a gene-specific manner. Although expression of GR was not altered, GR recruitment to transcriptional response elements was diminished in genes where neonatal genistein exposure blunted the glucocorticoid response. Furthermore, the ability of glucocorticoids to function as biological antagonists to the uterotrophic actions of E 2 was blocked by early-life exposure to genistein. Selectively altered glucocorticoid signaling may in part contribute to the subfertile/infertile phenotype following environmental estrogen exposure, though further study is warranted to determine the mechanisms by which disrupted glucocorticoid signaling impairs uterine receptivity.