The effect of 17α-ethynylestradiol on steroidogenesis and gonadal cytokine gene expression is related to the reproductive stage in marine hermaphrodite fish

Embryo developmental toxicity in marine medaka (Oryzias melastigma) due to parental and embryonic 17α-ethinylestradiol exposure

The synthetic estrogen 17α-ethinylestradiol (EE2) is a common component of hormone therapy and oral contraceptives and has been widely used for nearly 60 years. Numerous studies have shown that exposure to EE2 can affect embryonic development in a number of fish species. The effects of parental and embryonic EE2 exposure on embryo developmental toxicity and the underlying molecular mechanisms, however, have rarely been examined. In this study, embryos collected from parental EE2-exposed adult fish were examined to assess EE2-induecd toxicity during embryo development. The rate of embryo development including heart rate, hatching rate, and larval locomotion were measured to assess embryo developmental toxicity. The embryonic transcriptome was used to delineate the related developmental toxicity pathways. Our results suggest that parental and embryonic EE2 exposure resulted in growth retardation including a reduction in embryo heart rate, a delay in the appearance eye pigmentation, decreased hatching rate and impaired larval locomotion. In addition, gene ontology (GO) enrichment analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, and Ingenuity Pathway Analysis (IPA) of transcriptome revealed that these impairments are controlled by estrogen receptor and related to eye structure, neuronal and synaptic structure, and behaviour. The key factors identified, including PRKAA2, APOB, EPHB2, OXTR, NR2E3, and POU4F2, could serve as biomarkers for assessing EE2-induced embryo developmental toxicity. For the first time, our results show that eye pigmentation is a potentially sensitive marker of EE2-induced embryo developmental toxicity.

The Effect of 17α-Ethynilestradiol and GPER1 Activation on Body and Muscle Growth, Muscle Composition and Growth-Related Gene Expression of Gilthead Seabream, Sparus aurata L

Endocrine-disrupting chemicals include natural and synthetic estrogens, such as 17α-ethynilestradiol (EE2), which can affect reproduction, growth and immunity. Estrogen signalling is mediated by nuclear or membrane estrogen receptors, such as the new G-protein-coupled estrogen receptor 1 (GPER1). The present work studies the effect of EE2 and G1 (an agonist of GPER1) on body and muscle parameters and growth-related genes of 54 two-year-old seabreams. The fish were fed a diet containing EE2 (EE2 group) and G1 (G1 group) for 45 days and then a diet without EE2 or G1 for 122 days. An untreated control group was also studied. At 45 days, the shortest body length was observed in the G1 group, while 79 and 122 days after the cessation of treatments, the shortest body growth was observed in the EE2 group. Hypertrophy of white fibers was higher in the EE2 and G1 groups than it was in the control group, whereas the opposite was the case with respect to hyperplasia. Textural hardness showed a negative correlation with the size of white fibers. At the end of the experiment, all fish analyzed in the EE2 group showed a predominance of the gonadal ovarian area. In addition, the highest expression of the mafbx gene (upregulated in catabolic signals) and mstn2 (myogenesis negative regulator) was found in EE2-exposed fish.

17α-ethynylestradiol prevents the natural male-to-female sex change in gilthead seabream (Sparus aurata L.)

Exposure to 17α-ethynylestradiol (EE2, 5 μg/g food) impairs some reproductive events in the protandrous gilthead seabream and a short recovery period does not allow full recovery. In this study, spermiating seabream males in the second reproductive cycle (RC) were fed a diet containing 5 or 2.5 μg EE2/g food for 28 days and then a commercial diet without EE2 for the remaining RC. Individuals were sampled at the end of the EE2 treatment and then at the end of the RC and at the beginning of the third RC, 146 and 333 days after the cessation of treatment, respectively. Increased hepatic transcript levels of the gene coding for vitellogenin (vtg) and plasma levels of Vtg indicated both concentrations of EE2 caused endocrine disruption. Modifications in the histological organization of the testis, germ cell proliferation, plasma levels of the sex steroids and pituitary expression levels of the genes coding for the gonadotropin β-subunits, fshβ and lhβ were detected. The plasma levels of Vtg and most of the reproductive parameters were restored 146 days after treatments. However, although 50% of the control fish underwent sex reversal as expected at the third RC, male-to female sex change was prevented by both EE2 concentrations.

Role of sox30 in regulating testicular steroidogenesis of common carp

Expression of transcription factors is crucial for the regulation of steroidogenesis and gonadal development in fish. SRY-related box (SOX) proteins regulate gene expression of various events related to vertebrate reproduction. This study reports the role of sox30 and its influence on sox9a/b in regulating testicular steroidogenesis of the common carp, Cyprinus carpio. Tissue distribution showed predominant expression of sox30 in gonads, while gonadal ontogeny indicated significant dimorphic expression of sox30 from 120 days post hatch. Higher sox30 transcripts during the spawning season, an elevation of sox30 after human chorionic gonadotropin induction, and 11-ketotestosterone (11-KT) treatment authenticate gonadotropin dependency. Treatment of 17α-methyl-di-hydroxy-testosterone to juvenile common carp for mono-sex induction, vis-à-vis elevated sox30 expression. Sox30 protein was detected abundantly in spermatocytes and spermatid/sperm of carp testis. Transient silencing of sox30 using small interfering RNAs decreased sox9a/b expression, lead to downregulation of certain molecule/factor, transcription factor, germ/stem cell marker, and steroidogenesis-related enzyme genes. Serum testosterone and 11-KT decreased significantly upon transient silencing of sox30, in vivo. Concomitantly, a reduction in testicular microsomal 11-β hydroxysteroid dehydrogenase activity was observed. These results demonstrate the influence of sox30 as well as sox9a/b in the regulation of testicular steroidogenesis in common carp.

Quantitative transcriptomics, and lipidomics in evaluating ovarian developmental effects in Atlantic cod (Gadus morhua) caged at a capped marine waste disposal site

In the present study, a previously capped waste disposal site at Kollevåg (Norway) was selected to study the effects of contaminant leakage on biomarkers associated with Atlantic cod (Gadus morhua) reproductive endocrinology and development. Immature cod were caged for 6 weeks at 3 locations, selected to achieve a spatial gradient of contamination, and compared to a reference station. Quantitative transcriptomic, and lipidomic analysis was used to evaluate the effects of the potential complex contaminant mixture on ovarian developmental and endocrine physiology. The number of expressed transcripts, with 0.75 log2-fold differential expression or more, varied among stations and paralleled the severity of contamination. Particularly, significant bioaccumulation of ∑PCB-7, ∑DDTs and ∑PBDEs were observed at station 1, compared to the other station, including the reference station. Respectively 1416, 698 and 719 differentially expressed genes (DEGs), were observed at stations 1, 2 and 3, compared to the reference station, with transcripts belonging to steroid hormone synthesis pathway being significantly upregulation. Transcription factors such as esr2 and ahr2 were increased at all three stations, with highest fold-change at Station 1. MetaCore pathway maps identified affected pathways that are involved in ovarian physiology, where some unique pathways were significantly affected at each station. For the lipidomics, sphingolipid metabolism was particularly affected at station 1, and these effects paralleled the high contaminant burden at this station. Overall, our findings showed a novel and direct association between contaminant burden and ovarian toxicological and endocrine physiological responses in cod caged at the capped Kollevåg waste disposal site.

Gonad-Specific Transcriptomes Reveal Differential Expression of Gene and miRNA Between Male and Female of the Discus Fish (Symphysodon aequifasciatus)

The discus fish (Symphysodon aequifasciatus) is an ornamental fish that is well-known around the world. In artificial reproduction, they must be matched by one male and one female, whereas phenotype investigation indicated that there are no significant differences in appearance between males and females, which causes great difficulties in the mating during artificial reproduction. So, it is of great importance to establish artificial sex identification methods for the discus fish. The molecular mechanism of the sexual dimorphism of the discus fish was previously unknown. In this study, we constructed six cDNA libraries from three adult testes and three adult ovaries and performed RNA sequencing for identifying sex-biased candidate genes and microRNAs (miRNAs). A total of 50,082 non-redundant genes (unigenes) were identified, of which 18,570 unigenes were significantly overexpressed in testes, and 11,182 unigenes were significantly overexpressed in ovaries. A total of 551 miRNAs were identified, of which 47 miRNAs were differentially expressed between testes and ovaries. Eight differentially expressed unigenes, seven differentially expressed miRNAs and one non-differential miRNA were validated by quantitative real-time polymerase chain reaction. Twenty-four of these differentially expressed miRNAs and their 15 predicted target genes constituted 41 miRNA–mRNA interaction pairs, and some of vital sex-related metabolic pathways were also identified. These results revealed these differentially expressed genes and miRNAs between ovary and testis might be involved in regulating gonadal development, sex determination, gametogenesis, and physiological function maintenance, and there are complex regulatory networks between genes and miRNAs. It can help us understand the molecular mechanism of the sexual dimorphism and obtain a high-efficiency sex identification method in the artificial reproduction process of the discus fish.

Endocrine disrupter chemicals affect the humoral antimicrobial activities of gilthead seabream males even upon the cease of the exposure

17α-ethynilestradiol (EE₂) and tamoxifen (Tmx) are pollutants world-wide distributed in aquatic environments. Gilthead seabream, Sparus aurata L., is highlighted as a species model of intensively culture in anthropogenic disturbed environments. The effects of these pollutants on gilthead seabream reproduction and some immune responses have been described but, the humoral innate antimicrobial activities have never received attention. In this work we analysed the latest in the plasma of gilthead seabream males of different ages and reproductive stages treated with 0, 2.5, 5 or 50 μg EE₂ or 100 μg Tmx g^-1 food during different times of exposure and of reverting to commercial diet (recovery). The peroxidase and protease activities decreased as the spermatogenesis of the first reproductive cycle (RC) proceeded in control fish. However, only protease and antiprotease activities showed different level at different stages of the second RC in control fish, but showed scarce disruption in fish treated with EE₂ or Tmx. Peroxidase and bactericide activities are more sensitive to EE₂, than to Tmx. The effects induced by EE₂ varied depending on the activity analyzed, the dose and the time of exposure and the reproductive stage and the age of the specimens.

Effects of Dietary Bisphenol A on the Reproductive Function of Gilthead Sea Bream (Sparus aurata) Testes

Bisphenol A (BPA), a known endocrine disrupting chemical (EDC), was administered by diet to gilthead sea bream (Sparus aurata) in order to study its effects on the endocannabinoid system (ECS) and gonadal steroidogenesis. 2-year-old male gilthead sea bream were fed with two different concentrations of BPA (LOW at 4 and HIGH at 4000 µg/kg body weight for 21 days during the reproductive season. Exposure to 4000 µg BPA/kg bw/day (BPA HIGH) reduced sperm motility and altered the straight-line velocity (VSL) and linearity (LIN). Effects on steroidogenesis were evident, with testosterone (T) being up-regulated by both treatments and 11-ketotestosterone (11-KT) down-regulated by BPA HIGH. Plasma levels of 17β-estradiol (E₂) were not affected. The Gonadosomatic Index (GSI) increased in the BPA HIGH group. Interestingly, the levels of endocannabinoids and endocannabinoid-like compounds were significantly reduced after both treatments. Unpredictably, a few changes were noticed in the expression of genes coding for ECS enzymes, while the receptors were up-regulated depending on the BPA dose. Reproductive markers in testis (leptin receptor (lepr), estrogen receptors (era, erb), progesterone receptors (pr) and the gonadotropin releasing hormone receptor (gnrhr)) were up-regulated. BPA induced the up-regulation of the hepatic genes involved in oogenesis (vitellogenin (vtg) and zona pellucida 1 (zp1)).

GPER-independent inhibition of adrenocortical cancer growth by G-1 involves ROS/Egr-1/BAX pathway

We previously demonstrated that treatment of the H295R adrenocortical cancer cell line with the non-steroidal, high-affinity GPER (G protein-coupled estrogen receptor 1) agonist G-1 reduced tumor growth in vitro and in vivo through a GPER independent action. Moreover, we observed that G-1 treatment induces cell-cycle arrest and apoptosis following a sustained ERK1/2 activation. However, the precise mechanisms causing these effects were not clarified. Starting from our preliminary published results, we performed a microarray study that clearly evidenced a strong and significative up-regulation of EGR-1 gene in H295R cells treated for 24h with micromolar concentration of G-1. The microarray findings were confirmed by RT-PCR and Western-blot analysis as well as by immunofluorescence that revealed a strong nuclear staining for EGR-1 after G-1 treatment. EGR-1 is a point of convergence of many intracellular signaling cascades that control tumor cell growth and proliferation as well as others that relate to cell death machinery. Here we found that the increased Egr-1 expression was a consequence of G-1-mediated ROS-dependent ERK activation that were promptly reversed by the presence of the antioxidant n-acetyl-cysteine. Finally, we observed that silencing EGR-1 gene expression reversed the main effects induced by G-1 in ACC cells, including upregulation of the negative regulator of cell cycle, p21^Waf1/Cip1 and the positive regulator of mitochondrial apoptotic pathway, BAX, as well as the cell growth inhibition. The identified ROS/MAPK/Egr-1/BAX pathway as a potential off-target effect of the G-1 could be useful in implementing the pharmacological approach for ACC therapy.

17α-Ethynylestradiol alters the peritoneal immune response of gilthead seabream

17α-Ethynylestradiol (EE₂), a synthetic estrogen used in most oral contraceptives pills and hormone replacement therapies, is found in many water bodies, where it can modulate the fish immune response. EE₂ acts as an endocrine disruptor in gilthead seabream, Sparus aurata L., a marine teleost fish of great economic value in Mediterranean aquaculture, as it induces hepatic vitellogenin gene (vtg) expression. Moreover, EE₂ also alters the capacity of gilthead seabream to appropriately respond to infection although it does not behave as an immunosuppressor. Nevertheless, these previous studies have mainly focused on the head kidney leukocytes and no information exists on peritoneal leukocytes, including mast cells. In the present work, juvenile gilthead seabream fish were fed a pellet diet supplemented with EE₂ for 76 days and intraperitoneally injected with hemocyanin plus imject alum adjuvant at the end of EE₂ treatment and 92 days later, and the peritoneal immune response was analyzed. EE₂ supplementation induced vtg expression but returned to basal levels by 3 months post-treatment. Interestingly, gilthead seabream peritoneal leukocytes express the genes encoding for the nuclear estrogen receptor α and the G protein-coupled estrogen receptor 1 and the dietary intake of EE₂ induced these expression. Moreover, EE₂ induced an inflammatory response in the peritoneal cavity in unvaccinated fish, which was largely maintained for several months after the cessation of the treatment. However, the impact of EE₂ in vaccinated fish was rather minor and transient. Taken together, the study provides fresh information about endocrine immune disruption, focusing on peritoneal leukocytes.

Estrogens Promote the Production of Natural Neutralizing Antibodies in Fish through G Protein-Coupled Estrogen Receptor 1

Natural antibodies play crucial roles in pathogen elimination, B-cell survival and homeostasis, and inflammatory and autoimmune diseases. Although estrogens are able to regulate both innate and adaptive immune responses, their role in the production of natural antibodies is unknown. Here, we show that the dietary intake of the synthetic estradiol analog, 17α-ethinylestradiol (EE₂), one of the most potent pharmaceutical estrogens and intensively used in human therapeutics as a component of most oral contraceptives, regulates the abundance and proliferation of T and IgM⁺ B lymphocytes in the teleost fish gilthead seabream (Sparus aurata L.). Furthermore, for the first time in vertebrates, it is shown that estrogen signaling through G protein-coupled estrogen receptor 1 (GPER1) induces the production of polyreactive natural antibodies, which are able to crossreact with unrelated antigens and commensal and pathogenic bacteria. In addition, the serum from fish treated with EE₂ or the GPER1 agonist G1 shows higher complement-dependent bactericidal activity than that from non-treated specimens. These results demonstrate that estrogens and GPER1 are the key regulators of natural antibody production and pathogen clearance in fish, paving the way for future studies in other vertebrate classes.

Physiological fidelity or model parsimony? The relative performance of reverse-toxicokinetic modeling approaches

Background

Physiologically-based toxicokinetic (PBTK) models are often developed to facilitate in vitro to in vivo extrapolation (IVIVE) using a top-down, compartmental approach, favoring architectural simplicity over physiological fidelity despite the lack of general guidelines relating model design to dynamical predictions. Here we explore the impact of design choice (high vs. low fidelity) on chemical distribution throughout an animal’s organ system.

Results

We contrast transient dynamics and steady states of three previously proposed PBTK models of varying complexity in response to chemical exposure. The steady states for each model were determined analytically to predict exposure conditions from tissue measurements. Steady state whole-body concentrations differ between models, despite identical environmental conditions, which originates from varying levels of physiological fidelity captured by the models. These differences affect the relative predictive accuracy of the inverted models used in exposure reconstruction to link effects-based exposure data with whole-organism response thresholds obtained from in vitro assay measurements.

Conclusions

Our results demonstrate how disregarding physiological fideltiy in favor of simpler models affects the internal dynamics and steady state estimates for chemical accumulation within tissues, which, in turn, poses significant challenges for the exposure reconstruction efforts that underlie many IVIVE methods. Developing standardized systems-level models for ecological organisms would not only ensure predictive consistency among future modeling studies, but also ensure pragmatic extrapolation of in vivo effects from in vitro data or modeling exposure-response relationships.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-017-0407-3) contains supplementary material, which is available to authorized users.

Estrogen regulation of gene expression in the teleost fish immune system

Elucidating the mechanisms of estrogens-induced immunomodulation in teleost fish is of great importance due to the observed worldwide continuing decrease in pristine environments. However, little is know about the immunotoxicological consequences of exposure to these chemicals in fish, or of the mechanisms through which these effects are mediated. In this review, we summarize the results showing estrogens (natural or synthetic) acting through estrogen receptors and regulating specific target genes, also through microRNAs (miRNAs), leading to modulation of the immune functioning. The identification and characterization of miRNAs will provide new opportunities for functional genome research on teleost immune system and can also be useful when screening for novel molecule biomarkers for environmental pollution.

Molecular mechanism of endocrine system impairment by 17α-methyltestosterone in gynogenic Pengze crucian carp offspring

The effects of synthetic androgen 17α-methyltestosterone (MT) on endocrine impairment were examined in crucian carp. Immature 7-month old mono-female Pengze crucian carp (Pcc) F2 offspring were exposed to 50 and 100 μg/L of MT (week 2, 4, and 8). Gonadosomatic index, hepatosomatic index and intestine weight altered considerably and oocyte development was repressed. In the treatment groups, ovarian 11-ketotestosterone decreased, whereas 17β-estradiol and testosterone increased, and ovarian aromatase activities increased at week 4. However, in the brain tissue, those values significantly decreased. Quantitative RT-PCR analysis demonstrated changes in steroid receptor genes and upregulation of steroidogenic genes (Pcc-3bhsd, Pcc-11bhsd2 Pcc-cyp11a1), while the other three steroidogenic genes (Pcc-cyp17a1, Pcc-cyp19a1a and Pcc-star) decreased from week 4 to week 8. Ovarian, hepatic Pcc-vtg B and vitellogenin concentration increased in both 50 and 100 μg/L of MT exposure groups. This study adds further information regarding the effects of androgens on the development of previtellogenic oocytes, which suggests that MT could directly target estrogen signaling pathway, or indirectly affect steroidogenesis and vitellogenesis.

Selective estrogen receptor modulators differentially alter the immune response of gilthead seabream juveniles

17α-ethynylestradiol (EE2), a synthetic estrogen used in oral contraceptives and hormone replacement therapy, tamoxifen (Tmx), a selective estrogen-receptor modulator used in hormone replacement therapy, and G1, a G protein-coupled estrogen receptor (GPER) selective agonist, differentially increased the hepatic vitellogenin (vtg) gene expression and altered the immune response in adult gilthead seabream (Sparus aurata L.) males. However, no information exists on the effects of these compounds on the immune response of juveniles. This study aims, for the first time, to investigate the effects of the dietary intake of EE2, Tmx or G1 on the immune response of gilthead seabream juveniles and the capacity of the immune system of the specimens to recover its functionality after ceasing exposures (recovery period). The specimens were immunized with hemocyanin in the presence of aluminium adjuvant 1 (group A) or 120 (group B) days after the treatments ceased (dpt). The results indicate that EE2 and Tmx, but not G1, differentially promoted a transient alteration in hepatic vtg gene expression. Although all three compounds did not affect the production of reactive oxygen intermediates, they inhibited the induction of interleukin-1β (il1b) gene expression after priming. Interestingly, although Tmx increased the percentage of IgM-positive cells in both head kidney and spleen during the recovery period, the antibody response of vaccinated fish varied depending on the compound used and when the immunization was administered. Taken together, our results suggest that these compounds differentially alter the capacity of fish to respond to infection during ontogeny and, more interestingly, that the adaptive immune response remained altered to an extent that depends on the compound.

Tamoxifen disrupts the reproductive process in gilthead seabream males and modulates the effects promoted by 17α-ethynylestradiol

17α-Ethynylestradiol (EE2), which is used in oral contraceptives and hormone replacement therapy, is a well documented estrogenic endocrine disruptor and an aquatic contaminant. In the present study, adult male specimens of the marine hermaphrodite teleost gilthead (Sparus aurata L.) were fed a diet containing tamoxifen (Tmx), an estrogen receptor ligand used in cancer therapy, alone or combined with EE2, for 25 days and then fed a commercial diet for a further 25 days (recovery period). The effects of short (5days) and long (25 days) treatments on several reproductive and gonad immune parameters and the reversibility of the disruptive effects after the recovery period were examined. Our data showed that Tmx acted as an estrogenic endocrine disruptor as revealed by the increase in the hepatic transcription of the vitellogenin gene in males, the serum levels of 17β-estradiol and the gonad expression levels of the estrogen receptor α and G protein-coupled estrogen receptor genes, and the recruitment of leukocytes into the gonad, a well known estrogenic-dependent process in gilthead seabream males. On the other hand, Tmx also increased sperm concentration and motility as well as the serum levels of androgens and the expression levels of genes that codify for androgenic enzymes, while decreasing the expression levels of the gene that code for gonadal aromatase. When applied simultaneously, Tmx and EE2 could act in synergy or counteract, each other, depending on the parameter measured. The disruptive effect of EE2 and/or Tmx was not reversible after a 25 day recovery period.

Tamoxifen persistently disrupts the humoral adaptive immune response of gilthead seabream (Sparus aurata L.)

There is increasing concern about the possible effect of pharmaceutical compounds may have on the fish immune system. Bath exposition of 17α-ethynylestradiol (EE2), a synthetic estrogen used in oral contraceptives, altered the immune response of the gilthead seabream (Sparus aurata L.), a marine hermaphrodite teleost. Tamoxifen (Tmx) is a selective estrogen-receptor modulator used in hormone replacement therapy, the effects of which are unknown in fish immunity. This study aims to investigate the effects of dietary administration of EE2 (5 μg/g food) and Tmx (100 μg/g food) on the immune response of gilthead seabream, and the capacity of the immune system to recover its functionality after a recovery period. The results show for the first time the reversibility of the effect of EE2 and Tmx on the fish immune response. Tmx promoted a transient alteration in hepatic vitellogenin gene expression of a different magnitude to that produced by EE2. Both, EE2 and Tmx inhibited the induction of interleukin-1β gene expression while reversed the inhibition of ROI production in leukocytes following vaccination. However, none of these effects were observed after ceasing EE2 and Tmx exposure. EE2 and Tmx stimulated the antibody response of vaccinated fish although Tmx, but not EE2, altered the antibody response and modulated the percentage of IgM(+) B lymphocytes of vaccinated fish during the recovery phase. Taken together, our results suggest that EE2 and Tmx might alter the capacity of fish to appropriately respond to infection and show that Tmx has a long-lasting effect on humoral adaptive immunity.

In Silico analysis of perturbed steroidogenesis and gonad growth in fathead minnows (P. promelas) exposed to 17α-ethynylestradiol

The multi-factorial nature of adverse reproductive effects mediated by endocrine disrupting compounds (or EDCs) makes understanding the mechanistic basis of reproductive dysfunction a highly pertinent area of research. As a consequence, a main motivator for continued research is to integrate 'multi-leveled' complexity (i.e., from genes to phenotype) using mathematical methods capable of encapsulating properties of physiological relevance. In this study, an in silico stoichiometric model of piscine steroidogenesis was augmented with a 'biomass' reaction associating the underlying stoichiometry of steroidogenesis with a reaction representative of gonad growth. The ability of the in silico model to predict perturbed steroidogenesis and subsequent effects on gonad growth was tested by exposing reproductively active male and female fathead minnows (Pimephales promelas) to 88 ng/L of the synthetic estrogen, 17α-ethynylestradiol (EE2). The in silico model was parameterized (or constrained) with experimentally quantified concentrations of selected steroid hormones (using mass spectrometry) and fold changes in gene expression (using RT-qPCR) for selected steroidogenic enzyme genes, in gonads of male and female fish. Once constrained, the optimization framework of flux balance analysis (FBA) was used to calculate an optimal flux through the biomass reaction (analogous to gonad growth) and associated steroidogenic flux distributions required to generate biomass. FBA successfully predicted effects of EE2 exposure on fathead minnow gonad growth (%gonadosomatic index or %GSI) and perturbed production of steroid hormones. Specifically, FBA accurately predicted no effects of exposure on male %GSI and a significant reduction for female %GSI. Furthermore, in silico simulations accurately identified disrupted reaction fluxes catalyzing productions of androgens (in male fish) and progestogens (in female fish), an observation which agreed with in vivo experimentation. The analyses presented is the first-ever to successfully associate underlying flux properties of the steroidogenic network with gonad growth in fish, an approach which can incorporate in silico predictions with toxicological risk assessments.