Genome-Wide Analysis of Low Dose Bisphenol-A (BPA) Exposure in Human Prostate Cells.

Endocrine disrupting compounds (EDCs) have the potential to cause adverse effects on wild-life and human health. Two important EDCs are the synthetic estrogen 17α-ethynylestradiol (EE2) and bisphenol-A (BPA) both of which are xenoestrogens (XEs) as they bind the estrogen receptor and dis-rupt estrogen physiology in mammals and other vertebrates. In the recent years the influence of XEs on oncogenes, specifically in relation to breast and prostate cancer has been the subject of considerable study.
METHODOLOGY: In this study, healthy primary human prostate epithelial cells (PrECs) were exposed to environmentally relevant concentrations of BPA (5nM and 25nM BPA) and interrogated using a whole genome microarray.
RESULTS: Exposure to 5 and 25nM BPA resulted in 7,182 and 7,650 differentially expressed (DE) genes, respectively in treated PrECs. Exposure to EE2 had the greatest effect on the PrEC transcriptome (8,891 DE genes).
CONCLUSION: We dissected and investigated the nature of the non-estrogenic gene signature associated with BPA with a focus on transcripts relevant to epigenetic modifications. The expression of transcripts encoding nuclear hormone receptors as well as histone and DNA methylation, modifying enzymes were significantly perturbed by exposure to BPA.

INTRODUCTION

Endocrine disrupting compounds (EDCs) are chemicals found in many consumer products that mimic natural hormones [1]. EDCs have been shown to disrupt normal endocrine function by binding nuclear hormone receptors including estrogen receptors and activating estrogen related signaling pathways, thus they represent an emerging problem for ecosystem and human health. A specific subset of EDCs, the xenoestrogens (XEs), can mimic 17β-estradiol (E2), the female estrogen [2]. One such XE is bisphenol-A (BPA), a chemical building block used to strengthen polycarbonate plastic products [3]. It is a global contaminant found in the environment at concentrations up to 92 nM in surface waters [4]. The primary source of BPA exposure is dietary. Canned foods, beer and sodas are frequently coated with BPA-containing epoxy films and BPA monomers can leach into the contents of the can. BPA has also been found in fast-food hamburgers [5].

Even though the prostate gland is androgen-dependent and develops during puberty under androgen control, estrogens also play a part in normal and malignant prostate development [6]. During the 3rd trimester of pregnancy, maternal estrogens allow the growth of the prostate gland in men [7], and it has previously been shown that embryos exposed to the exogenous estrogen diethylstilbestrol (DES) in that critical window present prostate hyperplasia as newborns [8].

In humans, estrogen receptor alpha (ERα) is encoded by the gene ESR1 (EStrogen Receptor 1) and estrogen receptor beta (ERβ) by the gene ESR2 (EStrogen Receptor 2). Prostate cells express ERα and ERβ, which are localized characteristically in stroma and epithelium, respectively [9]. Pasquali et al. [10] have provided evidence for the co-expression of both ESR1 and ESR2 transcripts in normal epithelial cells in primary cultures. In both human and rodent prostate, ERβ is more highly expressed than ERα, while ERβ is found in both epithelial and (to a lesser extent) stromal cells, ERα is generally localized to the stroma [9].

Estrogen used to be the foundation of prostate cancer therapy prior to the understanding that it was also associated with a high risk of serious cardiovascular complications, which ultimately limited its clinical use [11]. The indirect therapeutic effect of estrogen was primarily mediated by suppression of the hypothalamic-pituitary-gonadal axis and by direct effects on the Leydig cells in testis, which together lead to decreased serum levels of testosterone and castration like effects in prostate [12]. These studies also reported that estrogen had direct “toxic” effects on prostate.

In rodents, perinatal or neonatal exposure to 17β-estradiol (estradiol) or other estrogens can lead to “imprinting” of prostate and disrupt normal morphogenesis and inhibit prostate growth via ERα, resulting in aberrant prostate function and morphology in the adult prostate, as well as increased proliferation, inflammation and dysplasia, and cause epithelial hyperplasia including drastic alterations in the expression pattern of steroid receptors [13-15]. Prolonged treatment of adult rodents with estrogens also leads to epithelial metaplasia, prostatic intraepithelial neoplasia (PIN)-like lesions and even adenocarcinoma of the prostate highlighting the role of estrogen in prostate cancer development [9]. Chronic exposure to chlordecone, an estrogenic organochloride pesticide commonly used in the farming of bananas, is also associated with increased risk of prostate cancer [16].

BPA binds and activates the estrogen receptors (ERα and ERβ) [17, 18]. Additionally, BPA can stimulate other nuclear hormone receptors including estrogen-related receptor gamma (ESRRG) [19], the androgen receptor (NR3C4) [20], thyroid hormone receptors (THRs) [21], G-protein-coupled receptors (GPER) [22], glucocorticoid and mineralocorticoid receptors (NR3C1 and NR3C2) [23, 24], and the pregnane X receptor (NR1/2) [25]. BPA has also been shown to interfere with endocannabinoid receptors CNR1 and CNR2 [26, 27]. This promiscuous affinity for multiple nuclear hormone and G-protein coupled receptors results in BPA impacting multiple physiological processes in both humans and wildlife.

In the United States, BPA has been detected in 93% of 2,517 urine samples collected from both adults and children, indicating pervasive exposure to this contaminant [3, 28, 29]. Data from multiple sources indicate that its ubiquitous presence and thereby continuous exposure in humans cause adverse health effects including premature puberty [30], reduced fertility [31-33], obesity [34-37], metabolic diseases [38] and cancer [39-42]. The initial study to assess BPA’s estrogenic effects demonstrated that it increased the growth rate of MCF-7 cells [43]. Low (nM) concentrations of BPA exhibit estrogenic activity [44, 45]. Furthermore, a BPA metabolite, 4-methyl-2,4-bis(4-hydroxyphenyl)pent-1-ene (MBP) exhibits transcriptional activity at low nM concentrations and suggests that metabolism of BPA to MBP may explain some of the adverse effects observed [46]. Furthermore, a comparison of 3D structural models of human ERα and ERβ with both MBP and BPA revealed that MBP, but not BPA, had key interactions with amino acid residues important in E2 binding in human ERα and ERβ [47].

Since the late 1980s, prostate cancer has been the second leading cause of death due to cancer in men [48]. BPA has since been linked to increased risks for prostate cancer (PCa) by affecting cellular proliferation and prostate cancer cell migration [49]. Additionally, BPA has been shown to promote PCa by disrupting centrosome amplification and microtubule dynamics that contribute to neoplastic transformation of prostate cells [50, 51]. Moreover, a correlation exists between men affected by PCa and higher levels of BPA in urine than those without the disease, suggesting that urinary BPA levels may have prognostic value for PCa [51].

Interestingly, early-life exposure to BPA during critical developmental windows has also been shown to cause epigenetic alterations in the prostate, and, in so doing, promote prostate disease in men as they get older [52, 53]. BPA alters the epigenome by interfering with the expression of histone modifying enzymes for example, the histone deacetylase SIRT1, and the histone methyltransferase SET8, which cause gene expression modifications in PCa cells [54]. Cheong et al. [55] suggested that DNA methylation changes may be an epigenetic signature associated with increased cancer susceptibility in the adult gland due to early-life exposure.

The objective of this study was to perform a systems level analysis of the effects of environmentally relevant doses of BPA on the transcriptome of healthy primary human prostate epithelial cells (PrECs). As a positive estrogenic control for these experiments, we selected 17α-ethinylestradiol (EE2), an artificial analogue of E2, commonly used in birth control pills and a well-characterized XE [56, 57]. To our knowledge, this is the first in vitro study that (1) examines the effects of low dose BPA exposure on the healthy human prostate cell transcriptome using a system level approach, and (2) characterizes non-estrogenic effects of BPA by distinguishing how BPA exposure differs from EE2, (3) defines the difference between a low (5 nM) and higher dose (25 nM) of BPA on primary prostate epithelial cells (PrECs), and (4) specifically examines the effects of BPA on the expression of genes encoding histone and DNA methylation modifying enzymes.

MATERIALS AND METHODS

Primary Human Cells

Primary prostate epithelial cells (PrECs) derived from a healthy 23-year old male with no history of prostate disease were obtained from CloneticsTM (San Diego, CA). The PrEC used in this study were cryopreserved in the second passage as proliferating cultures and were not an immortalized cell line. Telomerase reverse transcriptase (TERT) expression negative immortalizes various normal cells, including prostate epithelial and stromal-derived cells in culture. The cells used in this study were TERT expression negative. PrECs were cultured in Clonetics™ PrEGM™ Prostate Epithelial Cell Growth Medium, with a confluence of 70-90% in 0.2 mL/well 96 well plates- and approximately 0.1 x 106 cells for each well - and supplemented with 10% fetal bovine serum, 2mM glutamine, and 1% penicillin/streptomycin. Within the first passage, PrECs were exposed to either 5 nM BPA, 25 nM BPA, or 0.1 nM EE2 for 24 hours. A control group received no treatment.

RNA Extraction and Microarray Analysis

RNA was extracted using TRIzol reagent (Invitrogen) and the extracted RNA was further purified using the RNeasy Mini kit (Qiagen, Valencia, CA). RNA was assessed via absorbance readings (OD) at 260nm using an ND-1000 (Nanodrop, Wilmington, DE). RNA integrity was examined with the 6000 Nano LabChip assay from Agilent, (Santa Clara, CA). Only RNA samples with an RNA integrity number (RIN) score of >7.0 indicating intact RNA were used for microarray experiments (RIN is a measure of the degradation of RNA sample) [58, 59]. The samples were then sent to the UCSD BIOGEM core for Illumina Beadarray processing. HumanHT-12 v4 Expression BeadChip Kits were used for the array experiments (Illumina, San Diego). RNA labeling and microarray hybridization have been described in detail elsewhere [60]. Microarray experiments were performed with three biological replicates. The arrays were scanned on the Illumina BeadArray Reader, a confocal-type imaging system with 532 (cye3) nm laser illuminations. Image analysis and data extraction was carried out as in accordance with Illumina specifications. Preliminary data analysis and QC was carried out using GenomeStudio version 1.8.0 (Illumina Inc.). All array data have been deposited in the EBI ArrayExpress Database accession number E-MTAB-7576.

Gene and Systems Level Analysis

Microarray data were analyzed (Fig. ) using the Pipeline for Integrated Microarray Expression and Normalization Toolkit (PIMENTo) [61, 62]. The tool follows best practices for array data analyses and streamlines the processes necessary for gene expression analysis. It is built with the R programming language and leverages several open-source packages available through CRAN and Bioconductor, routinely used for microarray analyses. PIMENTo enables researchers to perform complex tasks with a minimal number of operations. It allows the user to review necessary data inputs, provides (1) data pre-processing and quality control, (2) normalization to remove the technical variability across arrays data [63, 64], (3) data visualization and (4) differential expression (DE) analysis. Three biological replicate cultures exposed to BPA or EE2 were used in the analyses. Expression level data from the Illumina Bead Studio software were normalized using a loess algorithm [42]. Probes, whose expression level that exceeded a threshold value in at least one sample, are called detected. The threshold value was found by inspection from the distribution plots of (log) expression levels. As one of the biological replicates was performed over 9 months after the other two we used the Bioconductor package 'limma' to remove technical batch effects. Limma allows significant differentially expressed genes to be determined based on the p-values, adjusted p-values or false discovery rate (FDR) estimates. The methodology follows the premise that transcripts are sorted according to their q-value, which is the smallest false discovery rate (FDR) at which the transcript is called significant. FDR is the expected fraction of false positive tests among significant tests. We set FDR at 0.1. We then search for groups of genes that collectively enrich GO terms or pathways. It should also be noted that the systems level analysis we performed (GO and Pathway) was itself subjected to FDR testing thereby adding rigor to the data analysis.

This was conducted using iPathwayGuide (Advaita Bioinformatics) [65]. Venn comparisons were performed using VENNY 2.1 online tool [66].

RESULTS

The Transcriptome of PrECs Reveals an Altered Signature in Response to BPA (Low and High) and EE2 Exposures

Using Venn diagrams, the list of significant differentially expressed (DE) microarray probes (q ≤ 0.1) after exposures to 5 and 25 nM BPA as well as 0.1 nM EE2 (positive estrogenic control) relative to control PrECs were compared (Fig. ). Exposure to EE2 had the strongest effect on the PrEC transcriptome and resulted in 11,241 DE probes corresponding to 8,891 transcripts. Exposure to 5 and 25 nM BPA perturbed the expression patterns of 8,876 and 9,525 probes corresponding to 7,182 and 7,650 unique transcripts respectively, indicating that BPA also impacts the PrECs transcriptome. The number of array probes with an absolute log2 fold change differences of 1 (linear fold-change of 2) in PrECs exposed to BPA and EE2 compared to control cells were 821, 755 and 715 for EE2, 5nM and 25 nM BPA respectively (Fig. S1). All three exposures have 7,011 DE probes in common (Fig. , intersection) while the 5 nM BPA exposure has 548 unique DE probes (orange), the 25 nM BPA has 593 unique DE probes (green) and EE has 2,389 unique DE probes (purple). Together, the low and high doses of BPA affected 1,839 DE probes not shared with the EE2 signature.

Meta-analysis (Advaita iPathwayGuide)

Meta-analysis carried out using iPathwayGuide using unique significant DE genes for all three exposures revealed that 24 pathways were commonly enriched (Fig. ). Amongst the 24 commonly enriched pathways all pathways were related to metabolism (Table , orange), DNA replication and cell cycle (light green), cancer (light yellow), inflammation and immune response (light blue), protein biogenesis (purple), biological adhesion (pink). The TGF-β signaling pathway (dark green), which is involved in many cellular processes including cell growth, cell differentiation, apoptosis, cellular homeostasis and other cellular functions, was also enriched in all three exposures.

Meta-analysis also revealed 276 Biological Process (BP) terms that were significantly enriched as part of the Gene Ontology (GO) analysis. Table displays the top 30 BP terms that relate to inflammation and immune response (light blue), metabolic pathways (orange), embryonic development (dark yellow), angiogenesis and blood circulation (turquoise), DNA replication and cell cycle (light green), biological adhesion (pink), cancer (light yellow), and protein biogenesis (purple). This suggests that BPA has a strong estrogenic effect similar to EE2 on several biological pathways and processes.

Additionally, the meta-analysis revealed that each of the three exposures affected PrECs differently with each yielding an exclusive gene expression signature. Exposure to 5 nM BPA perturbed 15 biological pathways (Fig. , orange), including apelin signaling (Table , dark green) a key regulator of angiogenesis/cell proliferation/energy metabolism regulation, cell metabolism (orange), fat digestion/absorption and lysosome (white), focal adhesion (pink), neurotransmission and taste transduction (grey) and inflammation/immune response (light blue). Exposure to 25 nM BPA enriched 5 unique pathways (Fig. , green) related to protein synthesis (Table , purple) and metabolism (orange) as well as antifolate resistance and the Fanconi anemia pathway (white). Exposure to EE2, the positive estrogenic control, affected 7 unique pathways (Fig. , purple), including circadian rhythm (Table 3C, white), metabolism (orange) and Wnt/HIF-1/PPAR signaling (dark green).

Together these results suggest that each exposure has a unique impact on specific biological pathways. For instance, 5 nM BPA has a strong impact on neurological pathways such as gap junction, glutamatergic synapse and serotonergic synapse, a signature that is not observed in 25 nM BPA and 0.1 nM EE2 exposures. These findings highlight the non-estrogenic signature of BPA exposure.

Expression Levels of BPA-activated Receptors in PrECs: Treated Versus Basal Expression

As BPA has been previously shown to bind and activate several receptors other than ERα and ERβ (listed in the introduction), we examined their expression after exposure to BPA and EE2 and compared all exposures to the control group in the meta-analysis. We determined that ESR1 was upregulated only by the 25nM BPA exposure (Fig. ), and ESR2 remained unaffected (not present in meta-analysis). Both ESRRA and THRA were upregulated after exposure to BPA and EE2, with 25nM BPA causing a robust increase in expression similar to the positive control EE2 (Fig. ). At the other end of the spectrum, BPA and EE2 exposures profoundly decreased GPER, NR3C1 and APLN expression levels and even the low-level BPA exposure caused similar decreases compared to EE2 (Fig. ).

Effects of BPA Exposure on Histone and DNA Methylation Modifying Enzymes

To assess the effects of BPA on epigenetic factors, we compiled a list of all histone and DNA methylation modifying enzymes present in our DE analysis and generated a bar plot (Fig. and Supplementary Table ). Two DNA methyltransferases (DNMT1 and DNMT3B), as well as 2 histone deacetylases (HDAC4 and SIRT5), 5 histone methyltransferases (EZH2, SETD2, SETD6, SETMAR, SUV39H1) and UHRF1 were present in our DE gene list and upregulated in all 3 exposures. In contrast, 2 histone deacetylases (HDAC9 and SIRT2) and 3 histone methyltransferases (MLL3, SETD5, SMYD3) were downregulated in all 3 exposures. Certain genes were only significantly upregulated (EHMT2, SETD8) and downregulated (HDAC3, HDAC5, HDAC8, SIRT1) by the positive estrogenic control EE2, while other genes were upregulated (SETD1A) or downregulated (MLL5) by both EE2 and 25nM BPA exposures. HDAC7A was only upregulated in 25nM BPA exposure.

DISCUSSION

BPA Induces Estrogen-like Response in PrECs

We demonstrated that the synthetic estrogen EE2, used here as estrogenic control, and BPA both have strong impacts on transcripts encoding proteins belonging to fundamental signaling pathways and biological processes in PrECs comparable to the effects of estrogen reported above, confirming that they can both be categorized as xenoestrogens. BPA (5 and 25 nM) and EE2 (0.1nM) shared approximately 54% of the DE genes and 24 biological pathways in our in vitro study. There is evidence from rodent models and human prostate cell lines that BPA can influence carcinogenesis by modulating prostate cancer cell proliferation, and for some tumors, stimulating progression [67, 68]. Our data showed that all three exposures enriched cell cycle, DNA replication and mismatch repair biological pathways, suggesting that BPA and EE2 induce at a molecular level the potential for cancer development and progression in prostate cells. This data is consistent with previous studies showing that BPA induces prostate cancer (PCa) cells migration via modulation of the ion channel protein expression involved in calcium entry and in cancer cell migration [49].

Additionally, all three exposures perturbed several biological pathways related to angiogenesis and blood circulation. In fact, new growth in the vascular network is important since the proliferation, as well as metastatic spread, of cancer cells depends on an adequate supply of oxygen and nutrients and the removal of waste products [69]. BPA has been shown previously to promote angiogenesis later in life following prenatal exposure [70], suggesting that early life exposure to BPA can imprint the epigenome and increase cancer risk later in life. Prins et al. [71] reported that perinatal exposure to BPA or to estradiol increases the susceptibility to develop prostate cancer either spontaneously or after a second estrogenic exposure during adulthood.

Our data showed a strong genomic signature of metabolic pathways with all exposures. Prenatal exposure to BPA exhibited important effects in adult males including food intake, an increase in body weight and liver weight, abdominal adipocyte mass, number and volume, and in serum leptin and insulin, and a decrease in serum adiponectin and in glucose tolerance [72], highlighting that BPA is both an endocrine and metabolic disruptor [7].

Inflammation and immune response pathways were also enriched by all three exposures. This is consistent with several studies demonstrating that pro- and anti-inflammatory response, as well as the production of natural killer cells, vary according to estrogen levels [73]. More specifically, our system level analysis revealed that all three exposures commonly affected cytokine-cytokine receptor interaction, leukocyte transendothelial migration and natural killer cell mediated cytotoxicity biological pathways as well as BP GO terms related to acute inflammatory response. BPA has been shown to induce an inflammation-like response in human adipocytes by inducing the release of inflammatory factors, such as IL-6 and IFN-γ, and activate pro-inflammatory pathways including JNK, JAK/STAT and NFκB pathways [74].

Non-estrogenic Effects of BPA in PrECs

The main goal of this study was to tease out the non-estrogenic effects of BPA in PrECs by examining the unique genomic signature BPA has in contrast to the positive estrogenic control EE2. Additionally, we also aimed to compare the 5 and 25 nM exposures. Our in vitro study revealed that BPA exposure (both at 5 and 25 nM) results in a BPA specific signature and that approximately 14% of the DE genes are strictly affected by BPA (and not by EE2), suggesting that BPA also induces a non-estrogenic response in PrECs. PrECs treated with 5 and 25 nM BPA had approximately the same number of DE genes (7,182 and 7,650 DE genes respectively) but 5 nM BPA perturbed about 3 times as many biological pathways and 4 times as many predicted miRNAs than 25 nM BPA did, suggesting that at low doses BPA has a greater biological impact on the PrEC transcriptome. The higher dose of BPA affected pathways related to protein biogenesis and metabolism-hormone biosynthesis whereas the lower dose had more far reaching effects including perturbation of the apelin signaling pathway, focal adhesion, glycolysis/gluconeogenesis biological pathways.

The apelin signaling pathway was the most significantly enriched pathway in the gene list that is unique to the 5 nM BPA exposure, suggesting that BPA interferes with apelin expression or binds to the G-protein-coupled apelin receptor (APJ). Apelin and APJ are implicated in different key physiological processes such as angiogenesis, cell proliferation and energy metabolism regulation [75, 76]. Wan et al. [77] have recently described a miR-224/apelin axis that regulates cell invasion and migration in prostate cancer. Additionally, a role in the regulation of blood pressure, cardiac function, and fluid homeostasis via the renin-angiotensin-aldosterone

system (RAS) has been suggested [78]. Given that a locally based RAS exists in the prostate [79], further examination of the effect of BPA on the apelin signaling pathway and how apelin interferes with the RAS in the prostate is required to gain a better understanding of the impact BPA can have on this organ. No direct evidence exists yet of direct interaction between BPA and APJ but a recent study demonstrated that 1 nM BPA exposure increased apelin expression and secretion in ovarian carcinoma cell lines via activation of the peroxisome proliferator-activated receptor γ (PPARγ) [80]. Of

note is the fact it has been shown previously that PPARγ is expressed on prostate-derived cells [81] and that BPA can bind to PPARγ [82].

Focal adhesion was another biological pathway significantly enriched in PrECs exposed to 5 nM BPA. This is consistent with another study in which BPA caused alterations in the focal adhesion pathway following fetal exposure in the peri-ductal stroma of mice [83]. This pathway is important for cell motility, cell proliferation, cell differentiation, regulation of gene expression and cell survival as these biological processes are regulated by cell-matrix adhesion [84, 85]. A recent study has shown that BPA promotes migration, invasion, and an increase in the number of focal contacts in breast cancer cells and induces activation of FAK, Src, and ERK2 migration-promoting kinases [86].

Exposure to 5 nM BPA also enriched the gap junction pathway. Gap junctions contain intercellular channels that allow direct communication between the cytosolic compartments of adjacent cells and the direct transfer of small molecules including ions, amino acids, nucleotides, second messengers and other metabolites essential for many physiological events, including embryonic development, electrical coupling, metabolic transport, apoptosis, and tissue homeostasis [87, 88]. It has previously been shown that compared to normal cells, gap‐junctional communication in malignant cells is either reduced (or undetectable) via a cAMP‐dependent signal transduction pathway [89]. Normal PrECs express mRNA transcripts for connexin-32 (GJB1), connexin-40 (GJA5) and connexin-43 (GJA1) whereas in malignant cells, the expression of GJA1 and GJB1 are decreased [90], signifying that as cancer cells become autonomous, a decrease in intercellular communicative activity occurs. In our datasets, expression of both GJB1 and GJA5 were not detected and GJA1 was significantly down-regulated in PrECs exposed to BPA and EE2 relative to control (Supplementary Fig ). This suggests that exposures to BPA and EE2 shifted the profile of normal PrECs towards a malignant signature associated with perturbation of gap junction communication.

To further examine this signature, we have to consider other BPA-activated receptors expressed in PrECs. BPA has been shown to bind to the androgen receptor (AR) present in prostate cells, causing an anti-androgenic effect [91-93]. Given that androgens and AR play a pivotal role in the expression of the male phenotype, this BPA-induced anti-androgen effect could have far-reaching impacts on male reproduction and physiology.

The nuclear estrogen-related receptor-gamma (ERRγ) receptor also has a good affinity for BPA and is expressed in PrECs [94-96]. ERRγ is an inducible transcription factor and

a downstream mediator of endocrine and metabolic signals such as glucagon, insulin, and endocannabinoids [97]. Nuclear receptors mediate gene transcription through interactions with coactivators and corepressors, a group of proteins that modulates chromatin structure and the recruitment of the basal transcription machinery [98, 99]. ERRγ is an in vivo receptor of BPA. In zebrafish during development ERRγ mediates BPA-induced malformations in otoliths, structures found in the saccule and utricle of the inner ear [19]. BPA strongly binds to ERRγ in a dose-dependent manner with an IC50 value of 13.1 nM [95]. We did not observe any perturbations in ERRγ expression levels upon EE2 or BPA treatments suggesting that BPA exposure does not influence the expression of ERRγ. Few studies have examined the function of ERR in the prostate but recent studies have shown that it is a crucial mediator of multiple endocrine and metabolic signals and it plays important roles in pathological conditions such as insulin resistance, alcoholic liver injury, and cardiac hypertrophy, and controls energy metabolism in the heart, skeletal muscle, and pancreatic β cells [97]. ERRγ has been shown to regulate mitochondrial-oxidative metabolism and promote angiogenesis, mitochondrial biogenesis, and oxidative remodeling [100, 101].

Effects of BPA on Histone and DNA Methylation Modifying Enzymes

The epigenetic effects of BPA were first described by Dolinoy et al. (2007) who showed that BPA exposure during pregnancy in Agouti viable (Avy) yellow mice causes hypomethylation of DNA on the epigenome of the offspring and a change in coat color [102]. Doherty et al. (2011) demonstrated that BPA promotes translation of the histone methyltransferase EZH2 (Enhancer of Zeste Homolog 2), which increases the histone modification H3K27me3, a mark of transcriptional repression [39, 103, 104]. Doshi et al. (2011) reported that neonatal BPA exposure in male rats led to hypermethylation of the promoter region of ESR1 and ESR2 and downregulation of both ESRs expression in adult testes, and increased DNMT3a and 3b expression by 2 fold, highlighting that BPA induces aberrant DNA methylation [105]. Santangeli et al. (2016) also described enrichment of H3K4me3 and H3K27me3 histone modifications in genes of interest and increased DNMT1 and DNMT3 expression levels in BPA exposed zebrafish [106]. Our data is consistent with these findings since we observed a significant increase in EZH2, DNMT1 and DNMT3B after exposure to BPA and EE2.

The gene encoding the transcription factor ubiquitin-like with PHD and ring finger domains 1 (UHRF1) controls the expression of topoisomerase IIα (TOP2A) and is involved in a wide range of physiological and pathological phenomena, including cancer development and metastasis, by controlling gene expression through regulating epigenetic mechanisms, including DNA methylation, histone deacetylation, histone methylation, and histone ubiquitination [107, 108]. In our datasets, UHRF1 is upregulated in all three exposures, a signature that is observed in many types of cancer, and contributes to cancer cell activation through hyper-methylation of tumor-suppressor genes such as BRCA1, CDKN2A, p73, and RASSF1, making UHRF1 a progression marker in cancer [109].

CONCLUSION

Low levels of BPA and EE2, comparable to those routinely detected in the environment, critically affect global transcriptomic responses in prostate epithelial cells. Comparing and contrasting BPA to EE2, we were able to tease out the estrogenic from non-estrogenic transcriptomic responses that BPA induces. Similar to EE2, (1) BPA impacts expression of genes belonging to fundamental pathways that regulates cell proliferation, cell cycle, DNA replication, mismatch repair and angiogenesis, all of which are relevant for cancer development and progression in prostate cells, (2) BPA impacts endocrine and metabolic transcripts in prostate cells, and (3) BPA perturbs expression of transcripts belonging to inflammation and immune response pathways.

Interestingly, BPA also has a non-estrogenic signature in PrECs and we showed in this study that BPA has profound effects on transcripts related to biological pathways involved in regulation of energy metabolism, protein biogenesis and metabolism-hormone biosynthesis, focal adhesion and gap junctions. Additionally, BPA perturbed the expression of the receptors AR and ERRγ, suggesting that its effects are driven via other receptors beside estrogen receptors.

We also show here that BPA can also deregulate expression of key epigenetic factors including EZH2, DNMT1, DNMT3B and UHRF1, highlighting that BPA-induced transcriptional perturbations have the potential to have an influence at the level of the epigenome.

In conclusion, we note that this study focused on an in vitro model using primary prostate cells and differences may exist with the in vivo BPA exposure response. The results presented here provide a compelling rationale for the need for future studies on BPA exposure in human tissues, using newer technologies such as RNAseq and mass spectrometry.

Table 1

Pathways that are commonly enriched in PrECs exposed to 5 nM BPA, 25 nM BPA and 0.1 nM EE2. Orange: pathways related to metabolism pathways, green: pathways related to DNA replication and cell cycle, yellow: pathways related to cancer, light blue: pathways related to inflammation and immune response, purple: pathways related to protein biogenesis, dark green: signaling pathways, pink: pathways related to biological adhesion. (

Pathways Common to All 3 Exposures	5 nM BPA (p-value)	25 nM BPA (p-value)	0.1 nM EE2 (p-value)
Alanine, aspartate and glutamate metabolism	5.11E-03	1.07E-02	9.81E-03
Biosynthesis of amino acids	1.36E-04	8.12E-04	1.70E-03
Carbon metabolism	2.98E-05	1.82E-04	1.04E-03
Cell cycle	9.63E-04	1.82E-03	1.37E-02
Central carbon metabolism in cancer	6.05E-04	1.01E-02	6.84E-03
Chemical carcinogenesis	2.68E-03	7.53E-03	3.01E-02
Cytokine-cytokine receptor interaction	5.00E-03	2.97E-03	1.22E-02
DNA replication	3.88E-06	1.56E-08	1.46E-07
ECM-receptor interaction	5.28E-06	6.99E-06	4.75E-04
HTLV-I infection	2.67E-03	1.69E-03	2.86E-03
Leukocyte transendothelial migration	4.34E-02	1.81E-02	3.71E-02
Metabolic pathways	5.32E-03	1.39E-03	4.28E-04
Mismatch repair	3.63E-02	4.50E-04	7.97E-05
Natural killer cell mediated cytotoxicity	1.60E-02	3.43E-02	4.87E-02
Protein digestion and absorption	4.46E-03	3.00E-03	3.44E-02
Proteoglycans in cancer	1.12E-02	1.06E-02	6.08E-03
Purine metabolism	1.20E-02	2.28E-04	6.75E-04
Pyrimidine metabolism	6.38E-03	1.06E-03	3.82E-03
Regulation of actin cytoskeleton	3.32E-02	2.24E-02	2.58E-02
Ribosome biogenesis in eukaryotes	6.62E-04	1.01E-03	1.28E-02
RNA transport	2.28E-02	3.28E-02	2.45E-02
Small cell lung cancer	4.79E-02	1.74E-02	1.70E-02
TGF-beta signaling pathway	4.30E-02	3.21E-02	4.19E-02

Table 2

Top 30 GO biological process terms that are commonly enriched in PrECs exposed to 5 nM BPA, 25 nM BPA and 0.1 nM EE2. Orange: BP terms related to metabolic pathways, green: BP terms related to DNA replication and cell cycle, light yellow: BP terms related to cancer, dark yellow: BP terms related to embryonic development, light blue: BP terms related to inflammation and immune response, purple: BP terms related to protein biogenesis, turquoise: BP terms related to angiogenesis and blood circulation, pink: pathways related to biological adhesion. (

Biological Process Terms Common to All 3 Exposures	5 nM BPA (p-value)	25 nM BPA (p-value)	0.1 nM EE2 (p-value)
Acute inflammatory response	3.81E-02	1.29E-02	2.62E-02
Acylglycerol metabolic process	3.20E-02	2.12E-02	7.64E-03
ADP metabolic process	5.60E-04	5.75E-03	3.30E-02
Alditol metabolic process	1.48E-03	1.67E-02	2.13E-03
Alpha-amino acid metabolic process	1.83E-02	6.40E-03	4.39E-02
Ameboidal-type cell migration	2.64E-03	1.49E-02	2.58E-02
Aminoglycan catabolic process	3.25E-03	3.59E-03	1.58E-02
Anatomical structure development	8.20E-07	1.90E-04	1.00E-04
Anatomical structure formation involved in morphogenesis	2.20E-04	2.69E-03	4.29E-03
Anatomical structure homeostasis	8.87E-03	2.34E-02	1.78E-02
Anatomical structure morphogenesis	1.40E-04	9.77E-03	1.01E-02
Angiogenesis	2.50E-05	9.50E-05	3.10E-05
Animal organ development	6.00E-07	3.90E-06	1.00E-04
Animal organ morphogenesis	5.76E-03	3.09E-03	4.65E-03
Animal organ regeneration	1.57E-03	9.07E-03	1.85E-02
ATP generation from ADP	3.70E-04	8.19E-03	3.75E-02
Attachment of spindle microtubules to kinetochore	2.08E-02	1.60E-03	2.64E-02
Biological adhesion	6.60E-04	2.60E-04	4.23E-02
Blood circulation	2.46E-03	2.50E-05	2.63E-03
Blood vessel development	1.40E-04	1.20E-04	1.10E-04
Blood vessel morphogenesis	9.30E-05	1.30E-04	5.00E-04
C-terminal protein amino acid modification	2.50E-03	4.89E-02	2.02E-02
C-terminal protein lipidation	1.34E-03	3.43E-02	2.22E-02
Calcium-dependent cell-cell adhesion via plasma membrane cell adhesion molecules	2.61E-02	3.80E-02	1.24E-02
Carbohydrate biosynthetic process	2.13E-02	4.60E-04	1.39E-03
Carbohydrate metabolic process	3.50E-05	2.40E-05	9.25E-03
Carboxylic acid biosynthetic process	9.28E-03	3.70E-05	9.00E-04
Carboxylic acid metabolic process	4.40E-05	1.60E-05	6.80E-04
Cardiovascular system development	3.10E-04	9.30E-05	1.60E-04
Cell adhesion	4.80E-04	2.80E-04	3.80E-02

Table 3

Pathways that are uniquely enriched in PrECs exposed to 5 nM BPA, 25 nM BPA and 0.1 nM EE2. Orange: pathways related to metabolic and hormone biosynthesis, purple: pathways related to protein biogenesis, pink: pathways related to biological adhesion, grey: pathways related to neurotransmission and taste transduction, light blue: pathways related to inflammation and immune response, dark green: cell signaling pathways. (

A) Pathways Unique to 5 nM BPA	p-value
Apelin signaling pathway	4.52E-02
Drug metabolism - cytochrome P450	3.71E-02
Fat digestion and absorption	2.34E-02
Focal adhesion	9.05E-03
Fructose and mannose metabolism	4.72E-03
Gap junction	3.78E-02
Glutamatergic synapse	4.73E-02
Glycolysis / Gluconeogenesis	1.10E-03
Lysosome	2.46E-02
Metabolism of xenobiotics by cytochrome P450	1.62E-02
Phenylalanine metabolism	4.25E-02
Pyruvate metabolism	8.55E-03
Rheumatoid arthritis	4.93E-02
Serotonergic synapse	2.73E-02
Taste transduction	2.69E-02
B) Pathways unique to 25 nM BPA	p-value
Antifolate resistance	1.90E-02
Fanconi anemia pathway	2.60E-02
Ribosome	9.87E-03
Thiamine metabolism	3.52E-02
Ubiquinone and other terpenoid-quinone biosynthesis	3.52E-02
C) Pathway name unique to 0.1 nM EE2	p-value
Circadian rhythm	4.26E-02
Glycosaminoglycan degradation	3.68E-02
Glyoxylate and dicarboxylate metabolism	9.81E-03
HIF-1 signaling pathway	3.11E-02
PPAR signaling pathway	4.89E-02
Steroid hormone biosynthesis	2.44E-02
Wnt signaling pathway	1.91E-02