Endocrine-disrupting chemical concentrations in follicular fluid and follicular reproductive hormone levels.

PURPOSE: To determine correlations between chemicals in follicular fluid (FF) and follicular reproductive hormone levels. METHODS: The analysis was part of a larger cohort study to determine associations between exposure to EDCs and in vitro fertilization (IVF) outcomes. FF was aspirated from a single leading follicle per participant. Demographics and data on exposure to EDCs were self-reported by the participants using a questionnaire. The concentrations of estradiol (E2), progesterone (PG), anti-Mullerian hormone (AMH), and inhibin B, as well as that of 12 phthalate metabolites and 12 phenolic chemicals were measured in each FF sample. Multivariate linear regression model was used to identify the drivers of hormone levels based on participant's age, BMI, smoking status, and chemical exposure for the monitored chemicals detected in more than 50% of the samples. Benjamini-Hochberg false discovery rate (FDR) correction was applied on the resulting p values (q value). RESULTS: FF samples were obtained from 72 women (mean age 30.9 years). Most of the phthalates and phenolic substances monitored (21/24, 88%) were identified in FF. Ten compounds (7 phthalate metabolites, 3 phenols) were found in more than 50% of samples. In addition, there were positive associations between E2 levels and mono-n-butyl phthalate (MnBP) (beta = 0.01) and mono-isobutyl phthalate (MiBP) (beta = 0.03) levels (q value < 0.05). CONCLUSION: Higher concentrations of several phthalate metabolites, present among others in personal care products, were associated with increased E2 levels in FF. The results emphasize the need to further investigate the mechanisms of action of such EDCs on hormonal cyclicity and fertility in women.

Corrigendum: Development of new approach methods for the identification and characterization of endocrine metabolic disruptors-a PARC project.

[This corrects the article DOI: 10.3389/ftox.2023.1212509.].

The Conflict between Regulatory Agencies over the 20,000-Fold Lowering of the Tolerable Daily Intake (TDI) for Bisphenol A (BPA) by the European Food Safety Authority (EFSA).

BACKGROUND: The European Food Safety Authority (EFSA) recommended lowering their estimated tolerable daily intake (TDI) for bisphenol A (BPA) 20,000-fold to 0.2 ng/kg body weight (BW)/day. BPA is an extensively studied high production volume endocrine disrupting chemical (EDC) associated with a vast array of diseases. Prior risk assessments of BPA by EFSA as well as the US Food and Drug Administration (FDA) have relied on industry-funded studies conducted under good laboratory practice protocols (GLP) requiring guideline end points and detailed record keeping, while also claiming to examine (but rejecting) thousands of published findings by academic scientists. Guideline protocols initially formalized in the mid-twentieth century are still used by many regulatory agencies. EFSA used a 21st century approach in its reassessment of BPA and conducted a transparent, but time-limited, systematic review that included both guideline and academic research. The German Federal Institute for Risk Assessment (BfR) opposed EFSA's revision of the TDI for BPA. OBJECTIVES: We identify the flaws in the assumptions that the German BfR, as well as the FDA, have used to justify maintaining the TDI for BPA at levels above what a vast amount of academic research shows to cause harm. We argue that regulatory agencies need to incorporate 21st century science into chemical hazard identifications using the CLARITY-BPA (Consortium Linking Academic and Regulatory Insights on BPA Toxicity) nonguideline academic studies in a collaborative government-academic program model. DISCUSSION: We strongly endorse EFSA's revised TDI for BPA and support the European Commission's (EC) apparent acceptance of this updated BPA risk assessment. We discuss challenges to current chemical risk assessment assumptions about EDCs that need to be addressed by regulatory agencies to, in our opinion, become truly protective of public health. Addressing these challenges will hopefully result in BPA, and eventually other structurally similar bisphenols (called regrettable substitutions) for which there are known adverse effects, being eliminated from all food-related and many other uses in the EU and elsewhere. https://doi.org/10.1289/EHP13812.

Innovative tools and methods for toxicity testing within PARC work package 5 on hazard assessment.

New approach methodologies (NAMs) have the potential to become a major component of regulatory risk assessment, however, their actual implementation is challenging. The European Partnership for the Assessment of Risks from Chemicals (PARC) was designed to address many of the challenges that exist for the development and implementation of NAMs in modern chemical risk assessment. PARC's proximity to national and European regulatory agencies is envisioned to ensure that all the research and innovation projects that are initiated within PARC agree with actual regulatory needs. One of the main aims of PARC is to develop innovative methodologies that will directly aid chemical hazard identification, risk assessment, and regulation/policy. This will facilitate the development of NAMs for use in risk assessment, as well as the transition from an endpoint-based animal testing strategy to a more mechanistic-based NAMs testing strategy, as foreseen by the Tox21 and the EU Chemical's Strategy for Sustainability. This work falls under work package 5 (WP5) of the PARC initiative. There are three different tasks within WP5, and this paper is a general overview of the five main projects in the Task 5.2 'Innovative Tools and methods for Toxicity Testing,' with a focus on Human Health. This task will bridge essential regulatory data gaps pertaining to the assessment of toxicological prioritized endpoints such as non-genotoxic carcinogenicity, immunotoxicity, endocrine disruption (mainly thyroid), metabolic disruption, and (developmental and adult) neurotoxicity, thereby leveraging OECD's and PARC's AOP frameworks. This is intended to provide regulatory risk assessors and industry stakeholders with relevant, affordable and reliable assessment tools that will ultimately contribute to the application of next-generation risk assessment (NGRA) in Europe and worldwide.

Development of new approach methods for the identification and characterization of endocrine metabolic disruptors-a PARC project.

In past times, the analysis of endocrine disrupting properties of chemicals has mainly been focused on (anti-)estrogenic or (anti-)androgenic properties, as well as on aspects of steroidogenesis and the modulation of thyroid signaling. More recently, disruption of energy metabolism and related signaling pathways by exogenous substances, so-called metabolism-disrupting chemicals (MDCs) have come into focus. While general effects such as body and organ weight changes are routinely monitored in animal studies, there is a clear lack of mechanistic test systems to determine and characterize the metabolism-disrupting potential of chemicals. In order to contribute to filling this gap, one of the project within EU-funded Partnership for the Assessment of Risks of Chemicals (PARC) aims at developing novel in vitro methods for the detection of endocrine metabolic disruptors. Efforts will comprise projects related to specific signaling pathways, for example, involving mTOR or xenobiotic-sensing nuclear receptors, studies on hepatocytes, adipocytes and pancreatic beta cells covering metabolic and morphological endpoints, as well as metabolism-related zebrafish-based tests as an alternative to classic rodent bioassays. This paper provides an overview of the approaches and methods of these PARC projects and how this will contribute to the improvement of the toxicological toolbox to identify substances with endocrine disrupting properties and to decipher their mechanisms of action.

Ranking food products based on estimating and combining their microbiological, chemical and nutritional risks: Method and application to Ready-To-Eat dishes sold in France.

Microbiological (M), chemical (C), and nutritional (N) risks associated with food products are usually assessed and managed independently by experts in public services or food companies. This can render difficult the comparison of food products in term of overall risk for the consumer. The objective of this study was to suggest a relatively simple method to (i) classify food products based on their M, C and N risks, and (ii) aggregate these risks and rank the food products accordingly. The method was developed and applied to 17 ready-to-eat (RTE) dishes available on the French market. With regard to food safety, the individual M and C risks were characterized considering likelihood and severity as recommended by the Codex Alimentarius. With regard to nutrition/health, the N risk was estimated based on the tendency of the dish to contribute to nutrient adequacy and to a healthy eating pattern. Finally, the outranking method PROMETHEE was applied to aggregate the three M, C, N risks and rank the food dishes. Food products were ranked relatively to each other, not in absolute terms. When we attributed the same weight to M, C and N risks, the RTE dish "Duck Parmentier" had the highest risk score while "Papillote of chicken, potatoes and small vegetables" and "Vegetarian plate vegetables and quinoa" had the lowest. However, this overall ranking changed according to the weight assigned to individual M, C and N risks, at least for food products whose scores varied according to risk types, such as"sushi discovery" (high M and C risks, low N risk). Since the risk ranking method developed here was built with assumptions and hypotheses related to the specific case study, more applications are needed to assess whether it can be generic. Nevertheless, this method is well grounded, objective, transparent, relatively fast and easy to set up. It might lead to further development of decision tools, particularly for consumers. This study paves the way towards food product multi-risk ranking.

The fate of bisphenol A, bisphenol S, and their respective glucuronide metabolites in ovarian cells.

Ovarian cells are critical for reproduction and steroidogenesis, which are functions that can be impacted by exposure to xenobiotics. As in other extra-hepatic tissues, biotransformation events may occur at the ovarian level. Such metabolic events deserve interest, notably as they may modulate the overall exposure and toxicity of xenobiotics. In this study, the comparative metabolic fate of two bisphenols was investigated in ovarian cells. Bisphenol A (BPA), a model endocrine disruptor, and its major substitute bisphenol S (BPS) were selected. Bovine granulosa cells (primary cultures) and theca explants (ex vivo tissue) were exposed for 24 hr to tritium-labeled BPA, BPS and their respective glucuronides (i.e. their major circulating forms), at concentrations consistent with low-dose exposure scenarios. Mass balance studies were performed, followed by radio-HPLC profiling. The capability of both cell compartments to biotransform BPA and BPS into their respective sulfo-conjugates was demonstrated, with sulfation being the predominant metabolic route. In theca, there was a significantly higher persistence of BPA (compared to BPS) residues over 24 hr. Moreover, only theca explants were able to deconjugate inactive BPA-glucuronide and BPS-glucuronide back into their biologically active aglycone forms. Deconjugation rates were demonstrated to be higher for BPS-G than for BPA-G. These findings raise concerns about the in situ direct release of bisphenols at the level of the ovary and demonstrate the relevance of exploring the biotransformation of bisphenols and their circulating metabolites in different ovarian cells with specific metabolic capabilities. This work also provides essential knowledge for the improved risk assessment of bisphenols.

Importance of relative binding of bisphenol A and bisphenol S to plasma proteins for predicting their in vivo potencies.

Many studies suggest that the potential impact of bisphenol S (BPS) as an endocrine disruptor is comparable to that of bisphenol A (BPA). However, in vitro-to-in vivo and from animal to human extrapolations require knowledge of the plasma free fraction of the active endocrine compounds. The present study aimed to characterise BPA and BPS binding to plasma proteins both in humans and different animal species. The plasma protein binding of BPA and BPS was assessed by equilibrium dialysis in plasma from adult female mice, rats, monkeys, early and late pregnant women as well as paired cord blood, early and late pregnant sheep and foetal sheep. The fraction of free BPA was independent of plasma concentrations and ranged between 4% and 7% in adults. This fraction was 2 to 3.5 times lower than that of BPS in all species except sheep, ranging from 3% to 20%. Plasma binding of BPA and BPS was not affected by the stage of pregnancy, BPA and BPS free fractions representing about 4% and 9% during early and late human pregnancy, respectively. These fractions were lower than the free fractions of BPA (7%) and BPS (12%) in cord blood. Our results suggest that similarly to BPA, BPS is extensively bound to proteins, mainly albumin. The higher fraction of free BPS compared to BPA may have implications for human exposure assessment since BPS free plasma concentrations are expected to be 2 to 3.5 times higher than that of BPA for similar plasma concentration.

Gestational exposure to bisphenol A induces region-specific changes in brain metabolomic fingerprints in sheep.

Fetal brain development depends on maternofetal thyroid function. In rodents and sheep, perinatal BPA exposure is associated with maternal and/or fetal thyroid disruption and alterations in central nervous system development as demonstrated by metabolic modulations in the encephala of mice. We hypothesized that a gestational exposure to a low dose of BPA affects maternofetal thyroid function and fetal brain development in a region-specific manner. Pregnant ewes, a relevant model for human thyroid and brain development, were exposed to BPA (5 µg/kg bw/d, sc). The thyroid status of ewes during gestation and term fetuses at delivery was monitored. Fetal brain development was assessed by metabolic fingerprints at birth in 10 areas followed by metabolic network-based analysis. BPA treatment was associated with a significant time-dependent decrease in maternal TT4 serum concentrations. For 8 fetal brain regions, statistical models allowed discriminating BPA-treated from control lambs. Metabolic network computational analysis revealed that prenatal exposure to BPA modulated several metabolic pathways, in particular excitatory and inhibitory amino-acid, cholinergic, energy and lipid homeostasis pathways. These pathways might contribute to BPA-related neurobehavioral and cognitive disorders. Discrimination was particularly clear for the dorsal hippocampus, the cerebellar vermis, the dorsal hypothalamus, the caudate nucleus and the lateral part of the frontal cortex. Compared with previous results in rodents, the use of a larger animal model allowed to examine specific brain areas, and generate evidence of the distinct region-specific effects of fetal BPA exposure on the brain metabolome. These modifications occur concomitantly to subtle maternal thyroid function alteration. The functional link between such moderate thyroid changes and fetal brain metabolomic fingerprints remains to be determined as well as the potential implication of other modes of action triggered by BPA such as estrogenic ones. Our results pave the ways for new scientific strategies aiming at linking environmental endocrine disruption and altered neurodevelopment.

Transgenerational metabolomic fingerprints in mice ancestrally exposed to the obesogen TBT.

BACKGROUND: Endocrine disrupting chemicals (EDCs) contribute to the etiology of metabolic disorders such as obesity, insulin resistance and hepatic dysfunction. Concern is growing about the consequences of perinatal EDC exposure on disease predisposition later in life. Metabolomics are promising approaches for studying long-term consequences of early life EDC exposure. These approaches allow for the identification and characterization of biomarkers of direct or ancestral exposures that could be diagnostic for individual susceptibility to disease and help to understand mechanisms through which EDCs act. OBJECTIVES: We sought to identify metabolomic fingerprints in mice ancestrally exposed to the model obesogen tributyltin (TBT), to assess whether metabolomics could discriminate potential trans-generational susceptibility to obesity and recognize metabolic pathways modulated by ancestral TBT exposure. METHODS: We used non-targeted 1H NMR metabolomic analyses of plasma and liver samples collected from male and female mice ancestrally exposed to TBT in two independent transgenerational experiments in which F3 and F4 males became obese when challenged with increased dietary fat. RESULTS: Metabolomics confirmed transgenerational obesogenic effects of environmentally relevant doses of TBT in F3 and F4 males, in two independent studies. Although females never became obese, their specific metabolomic fingerprint evidenced distinct transgenerational effects of TBT in female mice consistent with impaired capacity for liver biotransformation. DISCUSSION: This study is the first application of metabolomics to unveil the transgenerational effects of EDC exposure. Very early, significant changes in the plasma metabolome were observed in animals ancestrally exposed to TBT. These changes preceded the onset of obesogenic effects elicited by increased dietary fat in the TBT groups, and which ultimately resulted in significant changes in the liver metabolome. Development of metabolomic fingerprints could facilitate the identification of individuals carrying the signature of ancestral obesogen exposure that might increase their susceptibility to other risk factor such as increased dietary fat.

Safety assessment of cosmetics by read across applied to metabolomics data of in vitro skin and liver models.

As a result of the cosmetics testing ban, safety evaluations of cosmetics ingredients must now be conducted using animal-free methods. A common approach is read across, which is mainly based on structural similarities but can also be conducted using biological endpoints. Here, metabolomics was used to assess biological effects to enable a read across between a candidate cosmetic ingredient, DIV665, only studied using in vitro assays, and a structurally similar reference compound, PA102, previously investigated using traditional in vivo toxicity methods. The (1) cutaneous distribution after topical application, (2) skin metabolism, (3) liver metabolism and (4) effect on the intracellular metabolomic profiles of in vitro skin and hepatic models, SkinEthic®RHE model and HepaRG® cells were investigated. The compounds exhibited similar skin penetration and skin and liver metabolism, with small differences attributed to their physicochemical properties. The effects of both compounds on the metabolome of RHE and HepaRG® cells were similarly small, both in terms of the metabolites modulated and the magnitude of changes. The patterns of metabolome changes did not fit with any known signature relating to a mode of action known to be linked to liver toxicity e.g. modification of the Krebs cycle, urea synthesis and lipid metabolism, were more reflective of transient adaptive responses. Overall, these studies indicate that PA102 is biologically similar to DIV665, allowing read across of safety endpoints, such as in vivo sub-chronic (but not reproduction toxicity) studies, for the former to be applied to DIV665. Based on this, in the absence of animal data (which is prohibited for new chemicals), it could be concluded that DIV665 applied according to the consumer topical use scenario, is similar to PA102, and is predicted to exhibit low local skin and systemic toxicity.

The pregnane X receptor drives sexually dimorphic hepatic changes in lipid and xenobiotic metabolism in response to gut microbiota in mice.

BACKGROUND: The gut microbiota-intestine-liver relationship is emerging as an important factor in multiple hepatic pathologies, but the hepatic sensors and effectors of microbial signals are not well defined. RESULTS: By comparing publicly available liver transcriptomics data from conventional vs. germ-free mice, we identified pregnane X receptor (PXR, NR1I2) transcriptional activity as strongly affected by the absence of gut microbes. Microbiota depletion using antibiotics in Pxr+/+ vs Pxr-/- C57BL/6J littermate mice followed by hepatic transcriptomics revealed that most microbiota-sensitive genes were PXR-dependent in the liver in males, but not in females. Pathway enrichment analysis suggested that microbiota-PXR interaction controlled fatty acid and xenobiotic metabolism. We confirmed that antibiotic treatment reduced liver triglyceride content and hampered xenobiotic metabolism in the liver from Pxr+/+ but not Pxr-/- male mice. CONCLUSIONS: These findings identify PXR as a hepatic effector of microbiota-derived signals that regulate the host's sexually dimorphic lipid and xenobiotic metabolisms in the liver. Thus, our results reveal a potential new mechanism for unexpected drug-drug or food-drug interactions. Video abstract.

An Optimized Dual Extraction Method for the Simultaneous and Accurate Analysis of Polar Metabolites and Lipids Carried out on Single Biological Samples.

The functional understanding of metabolic changes requires both a significant investigation into metabolic pathways, as enabled by global metabolomics and lipidomics approaches, and the comprehensive and accurate exploration of specific key pathways. To answer this pivotal challenge, we propose an optimized approach, which combines an efficient sample preparation, aiming to reduce the variability, with a biphasic extraction method, where both the aqueous and organic phases of the same sample are used for mass spectrometry analyses. We demonstrated that this double extraction protocol allows working with one single sample without decreasing the metabolome and lipidome coverage. It enables the targeted analysis of 40 polar metabolites and 82 lipids, together with the absolute quantification of 32 polar metabolites, providing comprehensive coverage and quantitative measurement of the metabolites involved in central carbon energy pathways. With this method, we evidenced modulations of several lipids, amino acids, and energy metabolites in HepaRG cells exposed to fenofibrate, a model hepatic toxicant, and metabolic modulator. This new protocol is particularly relevant for experiments involving limited amounts of biological material and for functional metabolic explorations and is thus of particular interest for studies aiming to decipher the effects and modes of action of metabolic disrupting compounds.

Synergic toxic effects of food contaminant mixtures in human cells.

Humans are exposed to multiple exogenous substances, notably through food consumption. Many of these compounds are suspected to impact human health, and their combination could exacerbate their harmful effects. We previously observed in human cells that, among the six most prevalent food contaminant complex mixtures identified in the French diet, synergistic interactions between component appeared in two mixtures compared with the response with the chemicals alone. In the present study, we demonstrated in human cells that these properties are driven only by two heavy metals in each mixture: tellurium (Te) with cadmium (Cd) and Cd with inorganic arsenic (As), respectively. It appeared that the predicted effects for these binary mixtures using the mathematical model of Chou and Talalay confirmed synergism between these heavy metals. Based on different cell biology experiments (cytotoxicity, genotoxicity, mutagenesis and DNA repair inhibition experiments), a detailed mechanistic analysis of these two mixtures suggests that concomitant induction of oxidative DNA damage and decrease of their repair capacity contribute to the synergistic toxic effect of these chemical mixtures. Overall, these results may have broad implications for the fields of environmental toxicology and chemical mixture risk assessment.

The GOLIATH Project: Towards an Internationally Harmonised Approach for Testing Metabolism Disrupting Compounds.

The purpose of this project report is to introduce the European "GOLIATH" project, a new research project which addresses one of the most urgent regulatory needs in the testing of endocrine-disrupting chemicals (EDCs), namely the lack of methods for testing EDCs that disrupt metabolism and metabolic functions. These chemicals collectively referred to as "metabolism disrupting compounds" (MDCs) are natural and anthropogenic chemicals that can promote metabolic changes that can ultimately result in obesity, diabetes, and/or fatty liver in humans. This project report introduces the main approaches of the project and provides a focused review of the evidence of metabolic disruption for selected EDCs. GOLIATH will generate the world's first integrated approach to testing and assessment (IATA) specifically tailored to MDCs. GOLIATH will focus on the main cellular targets of metabolic disruption-hepatocytes, pancreatic endocrine cells, myocytes and adipocytes-and using an adverse outcome pathway (AOP) framework will provide key information on MDC-related mode of action by incorporating multi-omic analyses and translating results from in silico, in vitro, and in vivo models and assays to adverse metabolic health outcomes in humans at real-life exposures. Given the importance of international acceptance of the developed test methods for regulatory use, GOLIATH will link with ongoing initiatives of the Organisation for Economic Development (OECD) for test method (pre-)validation, IATA, and AOP development.

Evaluation of the genotoxic potential of apoptosis inducers with the γH2AX assay in human cells.

Human risk assessment of genotoxic chemicals is an important area of research. However, the specificity of in vitro mammalian genotoxicity assays is sometime low, as they yield to misleading positive results that are not observe in in vivo studies. Apoptosis can be a confounding factor in the interpretation of the results. Recently, a new strategy for genotoxicity screening, based on the combined analysis of phosphorylated histones H2AX (γH2AX) and H3 (pH3), was proposed to discriminate efficiently aneugenic from clastogenic compounds. However, γH2AX biomarker could also be induce by apoptosis. The aim of the present study was to investigate the specificity of this genotoxic biomarker. For this purpose, we analyzed 26 compounds inducing apoptosis by different mechanism of action, with the γH2AX assay in three human cell lines after 24 h treatment. Most of the tested chemicals were negative in the assay, whatever the cell line tested. The few compounds that generated positive data have also been report positive in other genotoxicity assays. The data presented here demonstrate that the γH2AX assay is not vulnerable to the generation of misleading positive results by apoptosis inducers. Currently, no formal guidelines have been approve for the γH2AX assay for regular genotoxicity studies, but we suggest that this biomarker could be used as a new standard genotoxicity assay.

Large-Scale Modeling Approach Reveals Functional Metabolic Shifts during Hepatic Differentiation.

Being able to explore the metabolism of broad metabolizing cells is of critical importance in many research fields. This article presents an original modeling solution combining metabolic network and omics data to identify modulated metabolic pathways and changes in metabolic functions occurring during differentiation of a human hepatic cell line (HepaRG). Our results confirm the activation of hepato-specific functionalities and newly evidence modulation of other metabolic pathways, which could not be evidenced from transcriptomic data alone. Our method takes advantage of the network structure to detect changes in metabolic pathways that do not have gene annotations and exploits flux analyses techniques to identify activated metabolic functions. Compared to the usual cell-specific metabolic network reconstruction approaches, it limits false predictions by considering several possible network configurations to represent one phenotype rather than one arbitrarily selected network. Our approach significantly enhances the comprehensive and functional assessment of cell metabolism, opening further perspectives to investigate metabolic shifts occurring within various biological contexts.

An Untargeted Metabolomics Approach to Investigate the Metabolic Modulations of HepG2 Cells Exposed to Low Doses of Bisphenol A and 17ß-Estradiol.

The model xeno-estrogen bisphenol A (BPA) has been extensively studied over the past two decades, contributing to major advances in the field of endocrine disrupting chemicals research. Besides its well documented adverse effects on reproduction and development observed in rodents, latest studies strongly suggest that BPA disrupts several endogenous metabolic pathways, with suspected steatogenic and obesogenic effects. BPA's adverse effects on reproduction are attributed to its ability to activate estrogen receptors (ERs), but its effects on metabolism and its mechanism(s) of action at low doses are so far only marginally understood. Metabolomics based approaches are increasingly used in toxicology to investigate the biological changes induced by model toxicants and chemical mixtures, to identify markers of toxicity and biological effects. In this study, we used proton nuclear magnetic resonance (1H-NMR) based untargeted metabolite profiling, followed by multivariate statistics and computational analysis of metabolic networks to examine the metabolic modulation induced in human hepatic cells (HepG2) by an exposure to low and very low doses of BPA (10-6M, 10-9M, and 10-12M), vs. the female reference hormone 17ß-estradiol (E2, 10-9M, 10-12M, and 10-15M). Metabolomic analysis combined to metabolic network reconstruction highlighted different mechanisms at lower doses of exposure. At the highest dose, our results evidence that BPA shares with E2 the capability to modulate several major metabolic routes that ensure cellular functions and detoxification processes, although the effects of the model xeno-estrogen and of the natural hormone can still be distinguished.