Expressomal approach for comprehensive analysis and visualization of ligand sensitivities of xenoestrogen responsive genes

Novel Approach to Screen Endocrine-Disrupting Chemicals via Endocrine-Enhanced Reduced Human Transcriptome

Endocrine-disrupting chemicals (EDCs) can interfere with multiple pathways and trigger different modes of action. Thus, the traditional EDC in vitro screening processes often require a battery of bioassays to cover multiple target pathways. Here we developed an endocrine-enhanced reduced human transcriptome (ERHT) focused on hormone receptor signaling induced by the EDCs regulating specific genes. ERHT was developed based on 1200 prioritized genes covering 110 endocrine-related biological pathways across eight potential adverse outcomes. The ability of this approach to identify EDCs was derived from machine learning of 1068 dose-dependent transcriptome profiles and enhanced by quantifying chemical-induced critical pathway responses, and thus, it demonstrated excellent classification performance (AUC = 0.84 ± 0.03) in internal cross-validation. We ultimately applied this approach to known EDCs and inactive substances to validate the reliability of this approach. Through external validation on 210 chemicals, the extrapolation accuracy exceeded 80%, demonstrating the outstanding practical performance of this approach. Meanwhile, the pathway responses induced by the same chemical were consistent with the experimental results reported by multiple sequencing platforms, highlighting the robustness of this approach. The above results demonstrate that this approach can provide novel insights for EDCs' high-throughput screening and comprehensive toxic mechanisms through biological pathways.

A transcriptomic biomarker predictive of cell proliferation for use in adverse outcome pathway-informed testing and assessment

High-throughput transcriptomics (HTTr) is increasingly being used to identify molecular targets of chemicals that can be linked to adverse outcomes. Cell proliferation (CP) is an important key event in chemical carcinogenesis. Here, we describe the construction and characterization of a gene expression biomarker that is predictive of the CP status in human and rodent tissues. The biomarker was constructed from 30 genes known to be increased in expression in prostate cancers relative to surrounding tissues and in cycling human MCF-7 cells after estrogen receptor (ER) agonist exposure. Using a large compendium of gene expression profiles to test utility, the biomarker could identify increases in CP in (i) 308 out of 367 tumor vs. normal surrounding tissue comparisons from 6 human organs, (ii) MCF-7 cells after activation of ER, (iii) after partial hepatectomy in mice and rats, and (iv) the livers of mice and rats after exposure to nongenotoxic hepatocarcinogens. The biomarker identified suppression of CP (i) under conditions of p53 activation by DNA damaging agents in human cells, (ii) in human A549 lung cells exposed to therapeutic anticancer kinase inhibitors (dasatinib, nilotnib), and (iii) in the mouse liver when comparing high levels of CP at birth to the low background levels in the adult. The responses using the biomarker were similar to those observed using conventional markers of CP including PCNA, Ki67, and BrdU labeling. The CP biomarker will be a useful tool for interpretation of HTTr data streams to identify CP status after exposure to chemicals in human cells or in rodent tissues.

A 50-gene biomarker identifies estrogen receptor-modulating chemicals in a microarray compendium

High throughput transcriptomics (HTTr) profiling has the potential to rapidly and comprehensively identify molecular targets of environmental chemicals that can be linked to adverse outcomes. We describe here the construction and characterization of a 50-gene expression biomarker designed to identify estrogen receptor (ER) active chemicals in HTTr datasets. Using microarray comparisons, the genes in the biomarker were identified as those that exhibited consistent directional changes when ER was activated (4 ER agonists; 4 ESR1 gene constitutively active mutants) and opposite directional changes when ER was suppressed (4 antagonist treatments; 4 ESR1 knockdown experiments). The biomarker was evaluated as a predictive tool using the Running Fisher algorithm by comparison to annotated gene expression microarray datasets including those evaluating the transcriptional effects of hormones and chemicals in MCF-7 cells. Depending on the reference dataset used, the biomarker had a predictive accuracy for activation of up to 96%. To demonstrate applicability for HTTr data analysis, the biomarker was used to identify ER activators in a set of 15 chemicals that are considered potential bisphenol A (BPA) alternatives examined at up to 10 concentrations in MCF-7 cells and analyzed by full-genome TempO-Seq. Using benchmark dose (BMD) modeling, the biomarker genes stratified the ER potency of BPA alternatives consistent with previous studies. These results demonstrate that the ER biomarker can be used to accurately identify ER activators in transcript profile data derived from MCF-7 cells.

Angiogenic Activity and Mechanism for Bisphenols on Endothelial Cell and Mouse: Evidence of a Structural-Selective Effect

Increased epidemiological evidence indicates the association of bisphenol exposure with human vascular disorders, while the underlying mechanism has not been clarified. Here, we sought to unveil the potential angiogenic effect and the underlying mechanism of bisphenols with different structural features using endothelial cells treated with an environmentally relevant concentration of bisphenols (range: 1 nM to 10 μM) and a C57BL/6 mouse model fed with doses of 0.002, 0.02, 2, and 20 mg/kg BW/day for 5 weeks. Bisphenol A (BPA) and bisphenol S (BPS) at a 1 nM level significantly increased tube formation by 45.1 and 30.2% and induced the microvessel sprouting, while tube length and microvessel sprouting were significantly inhibited by 37.2 and 55.7% after exposure to tetrabromobisphenol S (TBBPS) at 1 μM, respectively. Mechanistically, TBBPA and TBBPS significantly inhibited the interaction between phosphatidylinositol 3-kinase (PI3K) and thyroid receptor (TR), while BPA and BPS favored the interaction between PI3K and estrogen receptor (ER), resulting in abnormal PI3K signaling with consequent distinct angiogenic activity. BPA- and BPS-induced pro-angiogenic effects and TBBPS showed anti-angiogenic effects due to their distinct disruption on the TR/ER-PI3K pathway. Our work provided new evidence and mechanistic insight on the angiogenic activity of bisphenols and expanded the scope of endocrine disruptors with interference in vascular homeostasis.

Global gene expression analysis reveals novel transcription factors associated with long-term low-level exposure of EA.hy926 human endothelial cells to bisphenol A

Bisphenol A (BPA) is an endocrine disruptor that binds to estrogen receptors (ER); however, studies have shown that the ER pathway was not always the primary molecular mechanism of BPA's action in cells and that gene transcription could be altered by different exposure times and doses. Here, we sought to understand the correlation between the BPA-responsive genes that have associated biological functions and the transcription factors (TFs) involved in their regulation by repeatedly exposing human endothelial cells EA.hy926 to three nanomolar concentrations of BPA (10^-9 M, 10^-8 M, and 10^-7 M) for 14 weeks, after which changes in global gene expression were determined by RNA sequencing. Cytoscape plug-in iRegulon was used to infer TFs involved in the control of BPA-deregulated genes. The results show a minimal overlap in deregulated genes between three concentrations of BPA, with 10^-9 M BPA having the highest number of deregulated genes. TF analysis suggests that all three concentrations of BPA were active in the absence of an ER-mediated pathway. A unique set of TFs (NES≥4) has been identified for each BPA concentration, including the NFκB family and CEBPB for 10^-9 M BPA, MEF family, AHR/ARNT, and ZBTB33 for 10^-8 M BPA, and IRF1-7 and OVOL1/OVOL2 for 10^-7 M BPA, whereas STAT1/STAT2 were common TFs for 10^-9 M and 10^-7 M BPA. Overall, our data suggest that long-term low-level exposure of EA.hy926 cells to BPA leads to concentration-specific changes in gene expression that are not controlled by the ER-mediated signaling but rather by other mechanisms.

Toxicity testing and endocrine disrupting chemicals

Regulatory agencies around the world depend on standardized testing approaches to evaluate environmental chemicals for endocrine disrupting properties. The US Environmental Protection Agency (EPA) has developed a two-tiered testing approach within its Endocrine Disruptor Screening Program (EDSP). The eleven Tier 1 and three Tier 2 EDSP assays can be used to identify chemicals that act as agonists or antagonists of estrogen receptor, androgen receptor, or thyroid hormone receptor, or chemicals that interfere with steroidogenesis. Additional assays have been developed in the context of Tox21, and others have been validated by the OECD. In spite of the availability of validated toxicity tests, problems have been identified with the approaches and methods used to identify endocrine disrupting chemicals (EDCs). This chapter will provide an overview of several of these issues including: (1) The way an EDC is defined by an agency impacts whether a specific test can be used to determine if a chemical is an EDC. This is especially important when considering which assays examine outcomes that are considered "adverse effects." (2) Some assumptions about the validated studies used to identify EDCs may not be true (e.g., their reproducibility has been questioned). (3) Many of the validated assays are less sensitive than other methods that have not yet been validated. Ultimately, these and other problems contribute to the current landscape, where testing approaches have failed to protect the public from known EDCs. The chapter concludes with a review of approaches that have been taken to improve current guideline studies.

Similarities and Differences in Gene Expression Networks Between the Breast Cancer Cell Line Michigan Cancer Foundation-7 and Invasive Human Breast Cancer Tissues

Failure to adequately characterize cell lines, and understand the differences between in vitro and in vivo biology, can have serious consequences on the translatability of in vitro scientific studies to human clinical trials. This project focuses on the Michigan Cancer Foundation-7 (MCF-7) cells, a human breast adenocarcinoma cell line that is commonly used for in vitro cancer research, with over 42,000 publications in PubMed. In this study, we explore the key similarities and differences in gene expression networks of MCF-7 cell lines compared to human breast cancer tissues. We used two MCF-7 data sets, one data set collected by ARCHS4 including 1032 samples and one data set from Gene Expression Omnibus GSE50705 with 88 estradiol-treated MCF-7 samples. The human breast invasive ductal carcinoma (BRCA) data set came from The Cancer Genome Atlas, including 1212 breast tissue samples. Weighted Gene Correlation Network Analysis (WGCNA) and functional annotations of the data showed that MCF-7 cells and human breast tissues have only minimal similarity in biological processes, although some fundamental functions, such as cell cycle, are conserved. Scaled connectivity—a network topology metric—also showed drastic differences in the behavior of genes between MCF-7 and BRCA data sets. Finally, we used canSAR to compute ligand-based druggability scores of genes in the data sets, and our results suggested that using MCF-7 to study breast cancer may lead to missing important gene targets. Our comparison of the networks of MCF-7 and human breast cancer highlights the nuances of using MCF-7 to study human breast cancer and can contribute to better experimental design and result interpretation of study involving this cell line.

Bisphenol A promotes breast cancer cell proliferation by driving miR-381-3p-PTTG1-dependent cell cycle progression

Bisphenol A (BPA) is a high-production-volume industrial chemical that facilitates the development of breast cancer. However, the molecular mechanism associated with BPA-induced breast cancer cell proliferation and migration remains elusive. In our study, we exposed MCF-7 cells to different concentrations of BPA (0.1, 1 and 10 μM) for 24, 48, or 72 h. We found that BPA exposure significantly promoted MCF-7 cell proliferation and migration but not invasion. To elucidate the mechanisms, the differentially expressed genes between the BPA and control groups were investigated with the Gene Expression Omnibus (GEO) database through GEO2R. Kyoto Encyclopedia of Genes and Genomes (KEGG) and pathway action network analyses demonstrated the important role of the cell cycle pathway in the effects of BPA exposure on MCF-7 cells. Importantly, analysis with the cytoHubba plugin of Cytoscape software coupled with analysis of enriched genes in the cell cycle pathway identified PTTG1 and CDC20 (two hub genes) as key targets associated with BPA-induced MCF-7 cell proliferation and migration. Interestingly, BPA significantly increased the protein expression levels of PTTG1 but not CDC20. Knockdown of PTTG1 inhibited the BPA-induced increase in proliferation and maintained cell cycle progression. In addition, we confirmed that the increased expression of PTTG1 upon BPA exposure was caused by miR-381-3p inhibition. Moreover, we verified that miR-381-3p expression was low and inversely correlated with PTTG1 expression in breast cancer tissues. Together, these findings demonstrate that BPA promotes high PTTG1 expression and alters the cell cycle to enhance MCF-7 cell proliferation by inhibiting miR-381-3p expression.

Visceral obesity and cell cycle pathways serve as links in the association between bisphenol A exposure and breast cancer

It has been identified that bisphenol A (BPA) exposure causes developmental toxicity in breast cells. However, the exact molecular mechanisms underlying the association between exposure to BPA and breast cancer remain unclear. The aim of the present study was to investigate the BPA-regulated signaling pathways associated with the aggressiveness and the development of breast cancer. Microarray technology and functional gene set analyses were used to evaluate BPA and breast cancer-associated biomarkers and pathways in a discovery-driven manner. Using individual dataset analyses, it was indicated that two BPA-associated gene sets, the visceral obesity pathway, involved in visceral fat deposits and the metabolic syndrome, and the cell cycle pathway, involved in cyclins and cell cycle regulation, were significantly associated with a high grade of aggressiveness and the development of estrogen receptor (ER)-positive breast cancer (between P<0.05 and 0.0001). The pooled analysis indicated that the most significant pathway was G₁/S checkpoint regulation, and the cyclin and cell cycle regulation pathway for BPA-associated ER-positive cancer. Cancer cell signaling pathways were associated with healthy breast cells developing into breast cancer. The visceral obesity and the cell cycle pathways were indicated to link BPA exposure to breast cancer. The results of the present study demonstrate a significant association between breast cancer and BPA-regulated gene pathways.

PD-1 Imposes Qualitative Control of Cellular Transcriptomes in Response to T Cell Activation

Targeted blockade of programmed cell death 1 (PD-1), an immune-checkpoint receptor that inhibits T cell activation, provides clinical benefits in various cancers. However, how PD-1 modulates gene expression in T cells remains enigmatic. Here we investigated how PD-1 affects transcriptome changes induced by T cell receptor (TCR) activation. Intriguingly, we identified a huge variance in PD-1 sensitivity among TCR-inducible genes. When we quantified the half maximal effective concentration (EC₅₀) as the relationship between change in gene expression and TCR signal strength, we found that genes associated with survival and proliferation were efficiently expressed upon TCR activation and resistant to PD-1-mediated inhibition. Conversely, genes encoding cytokines and effector molecules were expressed less efficiently and sensitive to PD-1-mediated inhibition. We further demonstrated that transcription factor binding motifs and CpG frequency in the promoter region affect EC₅₀ and thus the PD-1 sensitivity of genes. Our findings explain how PD-1, dependent on the TCR signal strength, calibrates cellular transcriptomes to shape functional properties of T cell populations.

Dichotomy of bisphenol A-induced expression of peroxisome proliferator-activated receptors in hepatic and testicular tissues in mice

Environmental and dietary exposure to bisphenol A (BPA) and its toxicological consequences are extensively reported. BPA has multiple cellular targets. One of the mechanisms of action of BPA involves interaction with and activation of nuclear receptors (NRs) including peroxisome proliferator activated-receptors (PPARs). PPARs regulate genes involved in adipogenesis, and metabolism of glucose, lipid and cholesterol. Study of tissue and dose specific PPAR expression may decipher the toxicity outcome of BPA exposure. We studied expression of three forms of PPARs in mouse liver and testes exposed to BPA for 14 days. mRNA and protein expression of all forms of PPAR increased linearly (monotonic) with the dose in the liver while non-monotonic but dose specific effects were observed in the testes showing a differential pattern of expression. However, histopathological study showed a dose-dependent pattern of changes in liver as well as testes demonstrating a monotonic effect. These findings imply that other PPAR-independent mechanisms may play a role in BPA-induced pathological changes. The present study warrants exploration of the role of PPARs in BPA-induced effects on male reproductive functions and offers an insight into the peculiar response of BPA at low subchronic levels which may be helpful in designing appropriate risk assessment framework.

Weighted Gene Correlation Network Analysis (WGCNA) Reveals Novel Transcription Factors Associated With Bisphenol A Dose-Response

Despite Bisphenol-A (BPA) being subject to extensive study, a thorough understanding of molecular mechanism remains elusive. Here we show that using weighted gene correlation network analysis (WGCNA), which takes advantage of a graph theoretical approach to understanding correlations amongst genes and grouping genes into modules that typically have co-ordinated biological functions and regulatory mechanisms, that despite some commonality in altered genes, there is minimal overlap between BPA and estrogen in terms of network topology. We confirmed previous findings that ZNF217 and TFAP2C are involved in the estrogen pathway, and are implicated in BPA as well, although for BPA they appear to be active in the absence of canonical estrogen-receptor driven gene expression. Furthermore, our study suggested that PADI4 and RACK7/ZMYNDB8 may be involved in the overlap in gene expression between estradiol and BPA. Lastly, we demonstrated that even at low doses there are unique transcription factors that appear to be driving the biology of BPA, such as SREBF1. Overall, our data is consistent with other reports that BPA leads to subtle gene changes rather than profound aberrations of a conserved estrogen signaling (or other) pathways.

Alteration of mammary gland development by bisphenol a and evidence of a mode of action mediated through endocrine disruption

The development and function of the mammary gland are endocrine-dependent processes, depending on the stage of development. Foetal and/or postnatal exposure to low doses of BPA alters tissue organisation through epithelial proliferation and stroma-epithelial interactions. BPA also alters the expression of E2-dependent epithelial and stroma transcriptomes. Several signalling pathways are consistent with the observed phenotype: proliferation and apoptosis, a focal adhesion pathway indicating changes in biomechanical properties of the extracellular matrix, and immune function. Some of BPA's effects are reversed by oestrogen and/or GPER inhibitors. BPA also alters the expression of epigenetic marks (EZH2, HOTAIR), which would explain the delayed effect of foetal BPA exposure. In conclusion, experimental evidence shows that pre- or postnatal BPA exposure consistently causes endocrine modifications in the mammary tissue of different animal species, disrupting stromal-epithelial interactions and ultimately increasing its susceptibility to carcinogens. An interspecies comparison highlights why and how these effects apply to humans.

Endocrine Disruption by Mixtures in Topical Consumer Products

Endocrine disruption has been gathering increasing attention in the past 25 years as a possible new threat for health and safety. Exposure to endocrine disruptor has been progressively linked with a growing number of increasing disease in the human population. The mechanics through which endocrine disruptors act are not yet completely clear, however a number of pathways have been identified. A key concern is the cumulative and synergic effects that endocrine disruptors could have when mixed in consumer products. We reviewed the available literature to identify known or potential endocrine disruptors, as well as endocrine active substances that could contribute to cumulative effects, in topical consumer products. The number of endocrine actives used daily in consumer products is staggering and even though most if not all are used in concentrations that are considered to be safe, we believe that the possibility of combined effects in mixtures and non-monotonic dose/response is enough to require further precautions. A combined in vitro approach based on existing, validated OECD test methods is suggested to screen consumer products and mixtures for potential interaction with estrogen and androgen hormone receptors, in order to identify products that could have cumulative effects or support their safety concerning direct endocrine disruption capabilities.

Novel application of normalized pointwise mutual information (NPMI) to mine biomedical literature for gene sets associated with disease: use case in breast carcinogenesis

Advances in technology within biomedical sciences have led to an inundation of data across many fields, raising new challenges in how best to integrate and analyze these resources. For example, rapid chemical screening programs like the US Environmental Protection Agency's ToxCast and the collaborative effort, Tox21, have produced massive amounts of information on putative chemical mechanisms where assay targets are identified as genes; however, systematically linking these hypothesized mechanisms with in vivo toxicity endpoints like disease outcomes remains problematic. Herein we present a novel use of normalized pointwise mutual information (NPMI) to mine biomedical literature for gene associations with biological concepts as represented by Medical Subject Headings (MeSH terms) in PubMed. Resources that tag genes to articles were integrated, then cross-species orthologs were identified using UniRef50 clusters. MeSH term frequency was normalized to reflect the MeSH tree structure, and then the resulting GeneID-MeSH associations were ranked using NPMI. The resulting network, called Entity MeSH Co-occurrence Network (EMCON), is a scalable resource for the identification and ranking of genes for a given topic of interest. The utility of EMCON was evaluated with the use case of breast carcinogenesis. Topics relevant to breast carcinogenesis were used to query EMCON and retrieve genes important to each topic. A breast cancer gene set was compiled through expert literature review (ELR) to assess performance of the search results. We found that the results from EMCON ranked the breast cancer genes from ELR higher than randomly selected genes with a recall of 0.98. Precision of the top five genes for selected topics was calculated as 0.87. This work demonstrates that EMCON can be used to link in vitro results to possible biological outcomes, thus aiding in generation of testable hypotheses for furthering understanding of biological function and the contribution of chemical exposures to disease.

Editor's Highlight: Transcriptome Profiling Reveals Bisphenol A Alternatives Activate Estrogen Receptor Alpha in Human Breast Cancer Cells

Plasticizers with estrogenic activity, such as bisphenol A (BPA), have potential adverse health effects in humans. Due to mounting evidence of these health effects, BPA is being phased out and replaced by other bisphenol variants in “BPA-free” products. We have compared estrogenic activity of BPA with 6 bisphenol analogues [bisphenol S (BPS); bisphenol F (BPF); bisphenol AP (BPAP); bisphenol AF (BPAF); bisphenol Z (BPZ); bisphenol B (BPB)] in 3 human breast cancer cell lines. Estrogenicity was assessed (10⁻¹¹–10⁻⁴ M) by cell growth in an estrogen receptor (ER)-mediated cell proliferation assay, and by the induction of estrogen response element-mediated transcription in a luciferase assay. BPAF was the most potent bisphenol, followed by BPB > BPZ ∼ BPA > BPF ∼ BPAP > BPS. The addition of ICI 182,780 antagonized the activation of ERs. Data mining of ToxCast high-throughput screening assays confirm our results but also show divergence in the sensitivities of the assays. Gene expression profiles were determined in MCF-7 cells by microarray analysis. The comparison of transcriptome profile alterations resulting from BPA alternatives with an ERα gene expression biomarker further indicates that all BPA alternatives act as ERα agonists in MCF-7 cells. These results were confirmed by Illumina-based RNA sequencing. In conclusion, BPA alternatives are not necessarily less estrogenic than BPA in human breast cancer cells. BPAF, BPB, and BPZ were more estrogenic than BPA. These findings point to the importance of better understanding the risk of adverse effects from exposure to BPA alternatives, including hormone-dependent breast cancer.

BCScreen: A gene panel to test for breast carcinogenesis in chemical safety screening

Targeted gene lists have been used in clinical settings to specify breast tumor type, and to predict breast cancer prognosis and response to treatment. Separately, panels have been curated to predict systemic toxicity and xenoestrogen activity as a part of chemical screening strategies. However, currently available panels do not specifically target biological processes relevant to breast development and carcinogenesis. We have developed a gene panel called the Breast Carcinogen Screen (BCScreen) as a tool to identify potential breast carcinogens and characterize mechanisms of toxicity. First, we used four seminal reviews to identify 14 key characteristics of breast carcinogenesis, such as apoptosis, immunomodulation, and genotoxicity. Then, using a hybrid data and knowledge-driven framework, we systematically combined information from whole transcriptome data from genomic databases, biomedical literature, the CTD chemical-gene interaction database, and primary literature review to generate a panel of 500 genes relevant to breast carcinogenesis. We used normalized pointwise mutual information (NPMI) to rank genes that frequently co-occurred with key characteristics in biomedical literature. We found that many genes identified for BCScreen were not included in prognostic breast cancer or systemic toxicity panels. For example, more than half of BCScreen genes were not included in the Tox21 S1500+ general toxicity gene list. Of the 230 that did overlap between the two panels, representation varied across characteristics of carcinogenesis ranging from 21% for genes associated with epigenetics to 82% for genes associated with xenobiotic metabolism. Enrichment analysis of BCScreen identified pathways and processes including response to steroid hormones, cancer, cell cycle, apoptosis, DNA damage and breast cancer. The biologically-based systematic approach to gene prioritization demonstrated here provides a flexible framework for creating disease-focused gene panels to support discovery related to etiology. With validation, BCScreen may also be useful for toxicological screening relevant to breast carcinogenesis.

Evidence of Absence: Estrogenicity Assessment of a New Food-Contact Coating and the Bisphenol Used in Its Synthesis

Consumer concerns about exposure to substances found in food contact materials with estrogenic activity (EA) have created substantial demand for alternatives. We assessed the potential EA of both a new bisphenol monomer used to synthesize polymeric coatings for metal food-contact applications and the nonintentionally added substances (NIAS) that may migrate into food. We evaluated tetramethyl bisphenol F (TMBPF) using in vitro and in vivo assays. We extracted the polymeric coating using food simulants ethanol (50% v/v) and acetic acid (3% w/v) and measured migration using tandem liquid chromatography (LC)/mass spectrometry (MS) and LC time-of-flight MS for TMBPF and NIAS, respectively. We also tested migrants for EA using the E-SCREEN assay. TMBPF did not show estrogenic activity in the uterotrophic assay and did not alter puberty in male and female rats or mammary gland development in female rats. Neither TMBPF nor the migrants from the final polymeric coating increased proliferation of estrogen-sensitive MCF7 cells. TMBPF did not show estrogen-agonist or antagonist activity in the estrogen receptor-transactivation assay. TMBPF migration was below the 0.2 parts per billion detection limit. Our findings provide compelling evidence for the absence of EA by TMBPF and the polymeric coating derived from it and that human exposure to TMBPF would be negligible.

Combination of Metabolomics with Cellular Assays Reveals New Biomarkers and Mechanistic Insights on Xenoestrogenic Exposures in MCF-7 Cells

The disruptive potential of xenoestrogens like bisphenol A (BPA) lies in their 17β-estradiol (E2)-like binding to estrogen receptors (ERs) followed by concomitant modulation of ER target gene expression. Unsurprisingly, most endocrine testing systems focus on the quantification of canonical transcripts or ER-sensitive reporters. However, only little information is available about the corresponding metabolomic changes in vitro. This knowledge gap becomes particularly relevant in the context of potential mixture effects, for example, as a consequence of coexposure to potentially estrogenically active pollutants (e.g., Cd²⁺). Such effects are often difficult to dissect with molecular tools, especially with regard to potential physiological relevance. Metabolomic biomarkers are well-suited to address this latter aspect as they provide a comprehensive readout of whole-cell physiology. Applying a targeted metabolomics approach (FIA-MS/MS), this study looked for biomarkers indicative of xenoestrogenic exposure in MCF-7 cells. Cells were treated with E2 and BPA in the presence or absence of Cd²⁺. Statistical analysis revealed a total of 11 amino acids and phospholipids to be related to the compound's estrogenic potency. Co-exposure to Cd²⁺ modulated the estrogenic profile. However, the corresponding changes were found to be moderate with cellular assays such as the E-screen failing to record any Cd²⁺-specific estrogenic effects. Overall, metabolomics analysis identified proline as the most prominent estrogenic biomarker. Its increase could clearly be related to estrogenic exposure and concomitant ERα-mediated induction of proliferation. Involvement of the latter was confirmed by siRNA-mediated knockdown studies as well as by receptor inhibition. Further, the underlying signaling was also found to involve the oncoprotein MYC. Taken together, this study provides insights into the underlying mechanisms of xenoestrogenic effects and exemplify the strength of the complementary use of metabolomics and cellular and molecular assays.

Dose- and Time-Dependent Transcriptional Response of Ishikawa Cells Exposed to Genistein

To further define the utility of the Ishikawa cells as a reliable in vitro model to determine the potential estrogenic activity of chemicals of interest, transcriptional changes induced by genistein (GES) in Ishikawa cells at various doses (10 pM, 1 nM, 100 nM, and 10 μM) and time points (8, 24, and 48 h) were identified using a comprehensive microarray approach. Trend analysis indicated that the expression of 5342 unique genes was modified by GES in a dose- and time-dependent manner (P ≤ 0.0001). However, the majority of gene expression changes induced in Ishikawa cells were elicited by the highest dose of GES evaluated (10 μM). The GES' estrogenic activity was identified by comparing the Ishikawa cells' response to GES versus 17 α-ethynyl estradiol (EE, at equipotent doses, ie, 10 μM vs 1 μM, respectively) and was defined by changes in the expression of 284 unique genes elicited by GES and EE in the same direction, although the magnitude of the change for some genes was different. Further, comparing the response of the Ishikawa cells exposed to high doses of GES and EE versus the response of the juvenile rat uterus exposed to EE, we identified 66 unique genes which were up- or down regulated in a similar manner in vivo as well as in vitro Genistein elicits changes in multiple molecular pathways affecting various biological processes particularly associated with cell organization and biogenesis, regulation of translation, cell proliferation, and intracellular transport; processes also affected by estrogen exposure in the uterus of the rat. These results indicate that Ishikawa cells are capable of generating a biologically relevant estrogenic response and offer an in vitro model to assess this mode of action.

Fish connectivity mapping: linking chemical stressors by their mechanisms of action-driven transcriptomic profiles

Background

A very large and rapidly growing collection of transcriptomic profiles in public repositories is potentially of great value to developing data-driven bioinformatics applications for toxicology/ecotoxicology. Modeled on human connectivity mapping (Cmap) in biomedical research, this study was undertaken to investigate the utility of an analogous Cmap approach in ecotoxicology. Over 3500 zebrafish (Danio rerio) and fathead minnow (Pimephales promelas) transcriptomic profiles, each associated with one of several dozen chemical treatment conditions, were compiled into three distinct collections of rank-ordered gene lists (ROGLs) by species and microarray platforms. Individual query signatures, each consisting of multiple gene probes differentially expressed in a chemical condition, were used to interrogate the reference ROGLs.

Results

Informative connections were established at high success rates within species when, as defined by their mechanisms of action (MOAs), both query signatures and ROGLs were associated with the same or similar chemicals. Thus, a simple query signature functioned effectively as an exposure biomarker without need for a time-consuming process of development and validation. More importantly, a large reference database of ROGLs also enabled a query signature to cross-interrogate other chemical conditions with overlapping MOAs, leading to novel groupings and subgroupings of seemingly unrelated chemicals at a finer resolution. This approach confirmed the identities of several estrogenic chemicals, as well as a polycyclic aromatic hydrocarbon and a neuro-toxin, in the largely uncharacterized water samples near several waste water treatment plants, and thus demonstrates its future potential utility in real world applications.

Conclusions

The power of Cmap should grow as chemical coverages of ROGLs increase, making it a framework easily scalable in the future. The feasibility of toxicity extrapolation across fish species using Cmap needs more study, however, as more gene expression profiles linked to chemical conditions common to multiple fish species are needed.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2406-y) contains supplementary material, which is available to authorized users.

Nodes-and-connections RNAi knockdown screening: identification of a signaling molecule network involved in fulvestrant action and breast cancer prognosis

Although RNA interference (RNAi) knockdown screening of cancer cell cultures is an effective approach to predict drug targets or therapeutic/prognostic biomarkers, interactions among identified targets often remain obscure. Here, we introduce the nodes-and-connections RNAi knockdown screening that generates a map of target interactions through systematic iterations of in silico prediction of targets and their experimental validation. An initial RNAi knockdown screening of MCF-7 human breast cancer cells targeting 6560 proteins identified four signaling molecules required for their fulvestrant-induced apoptosis. Signaling molecules physically or functionally interacting with these four primary node targets were computationally predicted and experimentally validated, resulting in identification of four second-generation nodes. Three rounds of further iterations of the prediction–validation cycle generated third, fourth and fifth generation of nodes, completing a 19-node interaction map that contained three predicted nodes but without experimental validation because of technical limitations. The interaction map involved all three members of the death-associated protein kinases (DAPKs) as well as their upstream and downstream signaling molecules (calmodulins and myosin light chain kinases), suggesting that DAPKs play critical roles in the cytocidal action of fulvestrant. The in silico Kaplan–Meier analysis of previously reported human breast cancer cohorts demonstrated significant prognostic predictive power for five of the experimentally validated nodes and for three of the prediction-only nodes. Immunohistochemical studies on the expression of 10 nodal proteins in human breast cancer tissues not only supported their prognostic prediction power but also provided statistically significant evidence of their synchronized expression, implying functional interactions among these nodal proteins. Thus, the Nodes-and-Connections approach to RNAi knockdown screening yields biologically meaningful outcomes by taking advantage of the existing knowledge of the physical and functional interactions between the predicted target genes. The resulting interaction maps provide useful information on signaling pathways cooperatively involved in clinically important features of the malignant cells, such as drug resistance.

Regulation of epithelial-mesenchymal transition in breast cancer cells by cell contact and adhesion

Epithelial-mesenchymal transition (EMT) is a physiological program that is activated during cancer cell invasion and metastasis. We show here that EMT-related processes are linked to a broad and conserved program of transcriptional alterations that are influenced by cell contact and adhesion. Using cultured human breast cancer and mouse mammary epithelial cells, we find that reduced cell density, conditions under which cell contact is reduced, leads to reduced expression of genes associated with mammary epithelial cell differentiation and increased expression of genes associated with breast cancer. We further find that treatment of cells with matrix metalloproteinase-3 (MMP-3), an inducer of EMT, interrupts a defined subset of cell contact-regulated genes, including genes encoding a variety of RNA splicing proteins known to regulate the expression of Rac1b, an activated splice isoform of Rac1 known to be a key mediator of MMP-3-induced EMT in breast, lung, and pancreas. These results provide new insights into how MMPs act in cancer progression and how loss of cell-cell interactions is a key step in the earliest stages of cancer development.

Reproductive performance: at the cross-road of genetics, technologies and environment

Sexual reproduction depends on a negotiation between the sexes at the level of the cells (gametes), tissue (trophectoderm of the blastocyst and endometrium in the uterus) and organisms (to allow sexual intercourse). This review evaluates new questions linked to sexual reproduction in the biosphere in the context of the 21st century, in light of current knowledge in genetics and epigenetics. It presents the challenge of ‘forcing reproductive efficiency’ using ineffective gametes, or despite other fertility problems, through medically assisted reproduction and presents the reproductive challenge of high production farm animals, which are in a situation of chronically negative energy balance. It also analyses the situation created by the release of endocrine disruptors into the environment and discusses the possible transgenerational consequences of environmental modifications linked to these compounds.

A path forward in the debate over health impacts of endocrine disrupting chemicals

Several recent publications reflect debate on the issue of "endocrine disrupting chemicals" (EDCs), indicating that two seemingly mutually exclusive perspectives are being articulated separately and independently. Considering this, a group of scientists with expertise in basic science, medicine and risk assessment reviewed the various aspects of the debate to identify the most significant areas of dispute and to propose a path forward. We identified four areas of debate. The first is about the definitions for terms such as "endocrine disrupting chemical", "adverse effects", and "endocrine system". The second is focused on elements of hormone action including "potency", "endpoints", "timing", "dose" and "thresholds". The third addresses the information needed to establish sufficient evidence of harm. Finally, the fourth focuses on the need to develop and the characteristics of transparent, systematic methods to review the EDC literature. Herein we identify areas of general consensus and propose resolutions for these four areas that would allow the field to move beyond the current and, in our opinion, ineffective debate.

The effects of soy supplementation on gene expression in breast cancer: a randomized placebo-controlled study

BACKGROUND

There are conflicting reports on the impact of soy on breast carcinogenesis. This study examines the effects of soy supplementation on breast cancer-related genes and pathways.

METHODS

Women (n = 140) with early-stage breast cancer were randomly assigned to soy protein supplementation (n = 70) or placebo (n = 70) for 7 to 30 days, from diagnosis until surgery. Adherence was determined by plasma isoflavones: genistein and daidzein. Gene expression changes were evaluated by NanoString in pre- and posttreatment tumor tissue. Genome-wide expression analysis was performed on posttreatment tissue. Proliferation (Ki67) and apoptosis (Cas3) were assessed by immunohistochemistry.

RESULTS

Plasma isoflavones rose in the soy group (two-sided Wilcoxon rank-sum test, P < .001) and did not change in the placebo group. In paired analysis of pre- and posttreatment samples, 21 genes (out of 202) showed altered expression (two-sided Student's t-test, P < .05). Several genes including FANCC and UGT2A1 revealed different magnitude and direction of expression changes between the two groups (two-sided Student's t-test, P < .05). A high-genistein signature consisting of 126 differentially expressed genes was identified from microarray analysis of tumors. This signature was characterized by overexpression (>2-fold) of cell cycle transcripts, including those that promote cell proliferation, such as FGFR2, E2F5, BUB1, CCNB2, MYBL2, CDK1, and CDC20 (P < .01). Soy intake did not result in statistically significant changes in Ki67 or Cas3.

CONCLUSIONS

Gene expression associated with soy intake and high plasma genistein defines a signature characterized by overexpression of FGFR2 and genes that drive cell cycle and proliferation pathways. These findings raise the concerns that in a subset of women soy could adversely affect gene expression in breast cancer.