Influence of ovarian stage on transcript profiles in fathead minnow (Pimephales promelas) ovary tissue

Reaction thermodynamics as a constraint on piscine steroidogenesis flux distributions

While a considerable amount is known of the dynamics of piscine steroidogenesis during reproduction, the influence of thermodynamics constraints on its control has not been studied. In this manuscript, Gibbs free energy change of reactions was calculated for piscine steroidogenesis using the in silico eQuilibrator thermodynamics calculator. The analysis identified cytochrome P450 (cyp450) oxidoreductase reactions to have more negative Gibbs free energy changes relative to hydroxysteroid (HSD) and transferase reactions. In addition, a more favorable Gibbs free energy change was predicted for the Δ5 (cyp450 catalyzed) vs. Δ4 (HSD catalyzed) steroidogenesis branch-point, which converts pregnenolone to 17α-hydroxypregnenolone or progesterone respectively. Comparison of in silico predictions with in vivo experimentally measured flux across the Δ5 vs. Δ4 branch-point showed higher flux through the thermodynamically more favorable Δ5 pathway in reproducing or spawning vs. non-spawning fathead minnows (Pimephales promelas). However, the exposure of fish to endocrine stressors such as hypoxia or the synthetic estrogen 17α-ethinylestradiol (EE2), resulted in increased flux through both Δ5 and Δ4 pathways, indicating an adaptive response to increase steroidogenic redundancy. The correspondence of elevated flux through the Δ5 branch-point in spawning fish indicated the use of a thermodynamically favorable pathway to optimize steroid hormone productions during reproduction. We hypothesize that such selective use of a thermodynamically favorable steroidogenesis pathway may conserve reduced equivalents or transcriptional costs for investment to other biosynthetic or catabolic reactions to support reproduction. If generalizable, such an approach can provide novel insights into the structural principles and regulation of steroidogenesis or other metabolic pathways.

Transcriptome analysis of flathead grey mullet (Mugil cephalus) ovarian development induced by recombinant gonadotropin hormones

Background: Treatment with recombinant gonadotropin hormones (rGths), follicle-stimulating hormone (rFsh) and luteinizing hormone (rLh), was shown to induce and complete vitellogenesis to finally obtain viable eggs and larvae in the flathead grey mullet (Mugil cephalus), a teleost arrested at early stages of gametogenesis in intensive captivity conditions. This study aimed to investigate the transcriptomic changes that occur in the ovary of females during the rGths-induced vitellogenesis. Methods: Ovarian samples were collected through biopsies from the same five females at four stages of ovarian development. RNASeq libraries were constructed for all stages studied, sequenced on an Illumina HiSeq4000, and a de novo transcriptome was constructed. Differentially expressed genes (DEGs) were identified between stages and the functional properties of DEGs were characterized by comparison with the gene ontology and Kyoto Encyclopedia. An enrichment analysis of molecular pathways was performed. Results: The de novo transcriptome comprised 287,089 transcripts after filtering. As vitellogenesis progressed, more genes were significantly upregulated than downregulated. The rFsh application induced ovarian development from previtellogenesis to early-to-mid-vitellogenesis with associated pathways enriched from upregulated DEGs related to ovarian steroidogenesis and reproductive development, cholesterol metabolism, ovarian growth and differentiation, lipid accumulation, and cell-to-cell adhesion pathways. The application of rFsh and rLh at early-to-mid-vitellogenesis induced the growth of oocytes to late-vitellogenesis and, with it, the enrichment of pathways from upregulated DEGs related to the production of energy, such as the lysosomes activity. The application of rLh at late-vitellogenesis induced the completion of vitellogenesis with the enrichment of pathways linked with the switch from vitellogenesis to oocyte maturation. Conclusion: The DEGs and enriched molecular pathways described during the induced vitellogenesis of flathead grey mullet with rGths were typical of natural oogenesis reported for other fish species. Present results add new knowledge to the rGths action to further raise the possibility of using rGths in species that present similar reproductive disorders in aquaculture, the aquarium industry as well as the conservation of endangered species.

miR-133b targets tagln2 and functions in tilapia oogenesis

microRNAs (miRNAs) are important components of non-coding RNAs that participate in diverse life activities by regulating gene expression at the post transcriptional level through base complementary pairing with 3'UTRs of target mRNAs. miR-133b is a member of the miR-133 family, which play important roles in muscle differentiation and tumorigenesis. Recently, miR-133b was reported to affect estrogen synthesis by targeting foxl2 in mouse, while its role in fish reproduction remains to be elucidated. In the present study, we isolated the complete sequence of miR-133b, which was highly expressed in tilapia ovary at 30 and 90 dah (days after hatching) and subsequently decreased at 120 to 150 dah by qPCR. Interestingly, only a few oogonia were remained in the antagomir-133b treated tilapia ovary, while phase I and II oocytes were observed in the ovaries of the control group. Unexpectedly, the expression of foxl2 and cyp19a1a, as well as estradiol levels in serum were increased in the treated group. Furthermore, tagln2, an important factor for oogenesis, was predicted as the target gene of miR-133b, which was confirmed by dual luciferase reporter vector experiments. miR-133b and tagln2 were co-expressed in tilapia ovaries. Taken together, miR-133b may be involved in the early oogenesis of tilapia by regulating tagln2 expression. This study enriches the understanding of miR-133b function during oogenesis and lays a foundation for further study of the regulatory network during oogenesis.

Comparative Toxicological Evaluation of UV-Degraded versus Parent-Insensitive Munition Compound 1-Methyl-3-Nitroguanidine in Fathead Minnow

The US Army is replacing traditional munitions with insensitive munitions resistant to accidental detonation. Although the parent insensitive munition compound nitroguanidine (NQ) is generally not acutely toxic at concentrations >1000 mg/L in aquatic exposures, products formed by intensive ultraviolet (UV) degradation resulted in multiple-order of magnitude increases in toxicity. A methylated congener of NQ, 1-methyl-3-nitroguanidine (MeNQ), is also being assessed for potential use in insensitive munition explosive formulations; therefore, the present study investigated the hazard of parent versus UV-degraded MeNQ using fathead minnows (Pimephales promelas). Although up to 716 mg/L parent MeNQ caused no significant mortality or effects on growth in larval P. promelas fish in 7-d exposures, a similar concentration of MeNQ subjected to UV treatment resulted in 85% mortality. The UV treatment degraded only 3.3% of the MeNQ (5800 mg/L stock, UV-treated for 6 h), indicating that MeNQ degradation products have potentially high toxicity. The parent MeNQ exposure caused significantly decreased transcriptional expression of genes within the significantly enriched insulin metabolic pathway, suggesting antagonism of bioenergetics pathways, which complements observed, although nonsignificant, decreases in body weight. Significant differential transcriptional expression in the UV-degraded MeNQ treatments resulted in significant enrichment of pathways and functions related to the cell cycle, as well as erythrocyte function involved in O₂ /CO₂ exchange. These functions represent potential mechanistic sources of increased toxicity observed in the UV-degraded MeNQ exposures, which are distinct from previously observed mechanisms underlying increased toxicity of UV-degraded NQ in fish. Environ Toxicol Chem 2020;39:612-622. © 2019 SETAC.

Molecular Evaluation of Impacted Reproductive Physiology in Fathead Minnow Testes Provides Mechanistic Insights into Insensitive Munitions Toxicology

Previous toxicological investigations of the insensitive munition (IM), 3-nitro-1,2,4-triazol-5-one (NTO), demonstrated histopathological and physiological impacts in mammalian testes. The implications of these findings for fish was unknown, therefore we investigated the effects of chronic (21 day) exposures to NTO and an NTO-containing IM formulation called IMX-101 (composed of 2,4-dinitroanisole (DNAN), nitroguanidine (NQ), and NTO) in adult male fathead minnows to assess if impacts on testes were conserved. The NTO exposure caused no significant mortality through the maximum exposure concentration (720 mg/L, measured), however NTO elicited testicular impacts causing significant asynchrony in spermatogenesis and necrosis in secondary spermatocytes at the two highest exposure concentrations (383 mg/L and 720 mg/L) and testicular degeneration at the highest exposure. Microarray-based transcriptomics analysis identified significant enrichment of steroid metabolism pathways and mTORC-signal control of spermatogonia differentiation in NTO exposures each having logical connections to observed asynchronous spermatogenesis. Additionally, NTO impaired transcriptional expression for genes supporting sperm structural and flagellar development including sperm-associated antigen 6 (Spag6). These functional transcriptomic responses are hypothesized contributors to impacted reproductive physiology in NTO exposures that ultimately lead to reductions in spermatozoa. In contrast to NTO, the IMX-101 formulation elicited significant mortality at the two highest exposure concentrations of 25.2 and 50.9 mg/L (DNAN nominal + NTO measured + NQ measured). Unlike NTO and NQ, the DNAN component of the IMX-101 formulation underwent significant transformation in the 21d exposure. From previous investigations, neither NTO nor NQ caused mortality in fish at >1000 mg/L suggesting that mortality in the present study arose from DNAN / DNAN-attributable transformation products. The 12.6 mg/L IMX-101 exposure caused significant sublethal impacts on testes including sperm necrosis, interstitial fibrosis, and Sertoli-like cell hyperplasia. Transcriptional profiles for IMX-101 indicated significant enrichment on multiple signaling pathways supporting spermatogenesis, mitosis / meiosis, and flagellar structure, all logically connected to observed sperm necrosis. Additionally, pronounced transcriptional increases within the PPARα-RXRα pathway, a known DNAN target, has been hypothesized to correspond to Sertoli cell hyperplasia, presumably as a compensatory response to fulfill the nurse-function of Sertoli cells during spermatogenesis. Overall, the transcriptional results indicated unique molecular responses for NTO and IMX-101. Regarding chemical hazard, NTO impacted testes and impaired spermatogenesis, but at high exposure concentrations (≥ 192 mg/L), whereas the IMX-101 formulation, elicited mortality and impacts on reproductive physiology likely caused by DNAN and its transformation products present at concentrations well below the NTO-component concentration within the IMX-101 mixture formulation.

Alterations to the vision-associated transcriptome of zebrafish (Danio rerio) following developmental norethindrone exposure

Synthetic sex steroids, like the synthetic progestin norethindrone (NET), can affect a wide variety of biological processes via highly conserved mechanisms. NET is prevalent in surface waters, yet the sub-lethal effects of NET exposure are not are net yet well characterized in aquatic biota. A few targeted gene expression and behavioral studies have concluded that NET affects the vision of adult fish; however, early life stage (ELS) fish are often more sensitive to contaminants. Furthermore, many species of fish rely heavily on visual perception for survival during development. The goal of the present study was to characterize the effects of developmental exposure to environmentally relevant concentrations of NET on the visual system of ELS zebrafish, using transcriptomics and histological methods. Results indicate that exposure to relatively low levels of NET in aquatic systems may be sufficient to affect the visual function of developing fish.

Fer1l6 is essential for the development of vertebrate muscle tissue in zebrafish

The precise spatial and temporal expression of genes is essential for proper organismal development. Despite their importance, however, many developmental genes have yet to be identified. We have determined that Fer1l6, a member of the ferlin family of genes, is a novel factor in zebrafish development. We find that Fer1l6 is expressed broadly in the trunk and head of zebrafish larvae and is more restricted to gills and female gonads in adult zebrafish. Using both genetic mutant and morpholino knockdown models, we found that loss of Fer1l6 led to deformation of striated muscle tissues, delayed development of the heart, and high morbidity. Further, expression of genes associated with muscle cell proliferation and differentiation were affected. Fer1l6 was also detected in the C2C12 cell line, and unlike other ferlin homologues, we found Fer1l6 expression was independent of the myoblast-to-myotube transition. Finally, analysis of cell and recombinant protein-based assays indicate that Fer1l6 colocalizes with syntaxin 4 and vinculin, and that the putative C2 domains interact with lipid membranes. We conclude that Fer1l6 has diverged from other vertebrate ferlins to play an essential role in zebrafish skeletal and cardiac muscle development.

In silico predicted transcriptional regulatory control of steroidogenesis in spawning female fathead minnows (Pimephales promelas)

Oocyte development and maturation (or oogenesis) in spawning female fish is mediated by interrelated transcriptional regulatory and steroidogenesis networks. This study integrates a transcriptional regulatory network (TRN) model of steroidogenic enzyme gene expressions with a flux balance analysis (FBA) model of steroidogenesis. The two models were functionally related. Output from the TRN model (as magnitude gene expression simulated using extreme pathway (ExPa) analysis) was used to re-constrain linear inequality bounds for reactions in the FBA model. This allowed TRN model predictions to impact the steroidogenesis FBA model. These two interrelated models were tested as follows: First, in silico targeted steroidogenic enzyme gene activations in the TRN model showed high co-regulation (67-83%) for genes involved with oocyte growth and development (cyp11a1, cyp17-17,20-lyase, 3β-HSD and cyp19a1a). Whereas, no or low co-regulation corresponded with genes concertedly involved with oocyte final maturation prior to spawning (cyp17-17α-hydroxylase (0%) and 20β-HSD (33%)). Analysis (using FBA) of accompanying steroidogenesis fluxes showed high overlap for enzymes involved with oocyte growth and development versus those involved with final maturation and spawning. Second, the TRN model was parameterized with in vivo changes in the presence/absence of transcription factors (TFs) during oogenesis in female fathead minnows (Pimephales promelas). Oogenesis stages studied included: PreVitellogenic-Vitellogenic, Vitellogenic-Mature, Mature-Ovulated and Ovulated-Atretic stages. Predictions of TRN genes active during oogenesis showed overall elevated expressions for most genes during early oocyte development (PreVitellogenic-Vitellogenic, Vitellogenic-Mature) and post-ovulation (Ovulated-Atretic). Whereas ovulation (Mature-Ovulated) showed highest expression for cyp17-17α-hydroxylase only. FBA showed steroid hormone productions to also follow trends concomitant with steroidogenic enzyme gene expressions. General trends predicted by in silico modeling were similar to those observed in vivo. The integrated computational framework presented was capable of mechanistically representing aspects of reproductive function in fish. This approach can be extended to study reproductive effects under exposure to adverse environmental or anthropogenic stressors.

Hypoxia is involved in the reduction of HtrA3 in patients with endometrial hyperplasia and cancer

Endometrial cancer (EC) has recently become a major gynecological cancer and endometrial hyperplasia increases the risk for developing EC. Previous studies have reported that human high temperature requirement factor A3 (HtrA3), a member of ATP independent serine proteases family, is involved in endometrial carcinogenesis. However, the underlying mechanism of HtrA3 function is unclear in endometrial hyperplasia and cancer. In this study, we investigated that HtrA3 expression was reduced in endometrial hyperplasia as well as EC. The circulating levels of HtrA3 were also significantly reduced in both atypical hyperplasia and EC. Whether hypoxia is involved in the reduction of HtrA3 in EC was further investigated. Immunohistochemistry (IHC) scores of Glut1 and HtrA3 in type 1 and type 2 EC tissues showed the inverse correlation. And hypoxic condition reduced the expression of HtrA3. Furthermore, silencing HtrA3 promoted EC cell migration. Our study demonstrated the reduced levels of HtrA3 in endometrial hyperplasia including atypical hyperplasia which is a premalignant condition; and as the degree of hypoxia increases in EC, HtrA3 eventually loses its expression. Hypoxia is responsible for the reduction of HtrA3 which in turn promotes EC progression. These findings suggested that HtrA3 is an important adaptor in hypoxic regions that drives endometrial cancer development.

Transcriptomics provides mechanistic indicators of mixture toxicology for IMX-101 and IMX-104 formulations in fathead minnows (Pimephales promelas)

Within the US military, new insensitive munitions (IMs) are rapidly replacing conventional munitions improving safety from unintended detonation. Toxicity data for IM chemicals are expanding rapidly, however IM constituents are typically deployed in mixture formulations, and very little is known about their mixture toxicology. In the present study we sought to characterize the mixture effects and toxicology of the two predominant IM formulations IMX-101 and IMX-104 in acute (48 h) larval fathead minnow (Pimephales promelas) exposures. IMX-101 consists of a mixture of 2,4-dinitroanisole (DNAN), 3-nitro-1,2,4-triazol-5-one (NTO), and nitroguanidine (NQ) while IMX-104 is composed of DNAN, NTO, and hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX). DNAN was the most potent constituent in IMX-101 eliciting an LC50 of 36.1 mg/L, whereas NTO and NQ did not elicit significant mortality in exposures up to 1040 and 2640 mg/L, respectively. Toxic unit calculations indicated that IMX-101 elicited toxicity representative of the component concentration of DNAN within the mixture. Toxicogenomic responses for the individual constituents of IMX-101 indicated unique transcriptional expression and functional responses characteristic of: oxidative stress, impaired energy metabolism, tissue damage and inflammatory responses in DNAN exposures; impaired steroid biosynthesis and developmental cell-signaling in NQ exposures; and altered mitogen-activated protein kinase signaling in NTO exposures. Transcriptional responses to the IMX-101 mixture were driven by the effects of DNAN where expression and functional responses were nearly identical comparing DNAN alone versus the fractional equivalent of DNAN within IMX-101. Given that each individual constituent of the IMX-101 mixture elicited unique functional responses, and NTO and NQ did not interact with DNAN within the IMX-101 mixture exposure, the overall toxicity and toxicogenomic responses within acute exposures to the IMX-101 formulation are indicative of "independent" mixture toxicology. Alternatively, in the IMX-104 exposure both DNAN and RDX were each present at concentrations sufficient to elicit lethality (RDX LC50 = 28.9 mg/L). Toxic-unit calculations for IMX-104 mixture formulation exposures indicated slight synergistic toxicity (ΣTU LC50 = 0.82, 95% confidence interval = 0.73-0.90). Unique functional responses relative to DNAN were observed in the IMX-104 exposure including responses characteristic of RDX exposure. Based on previous transcriptomics responses to acute RDX exposures in fathead minnow larvae, we hypothesize that the potentially synergistic responses within the IMX-104 mixture are related to interactive effects of each DNAN and RDX on oxidative stress mitigation pathways.

Transcriptome analysis of three critical periods of ovarian development in Yellow River carp (Cyprinus carpio)

Ovary development is a complex process involving numerous genes; the molecular mechanism underlying the ovary development of carp is still unknown. Here we used Illumina HiSeq™ 2500 to explore the transcriptome of undifferentiated gland (PG), juvenile ovary (OJ) and adult ovary (OA) of Yellow River carp (Cyprinus carpio). A total of 58,749 unigenes were obtained, comprising 45,707 known genes and 13,042 new genes. We identified differentially-expressed genes (DEGs) during development and characterized the functional properties of DEGs by comparison with the gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes databases. qRT-PCR was used to analyze the expression of 22 DEGs and the results corresponded with those of RNA-Seq. Among DEGs between PG and OJ, some upstream regulators of gonad development were upregulated in PG, such as cyp19a and sox9, while some oocyte-specific genes were upregulated in OJ, such as nobox, bmp15 and zp2. Among DEGs between OJ and OA, many oocyte physiological function-related genes were upregulated in OA, such as fem-1 and foxl2. GO analysis showed a higher number of DEGs from PG-OJ analysis were assigned to reproduction terms. Furthermore, our investigation has also revealed DEGs identified from PG-OJ analysis were enriched in several important functional pathways, such as Fanconi anemia and the notch signal pathway. These data suggested a dynamic shift in gene expression during ovary development, and DEGs between PG and OJ provided crucial candidate gene data for the study of ovarian differentiation. Additionally, a total of 1,776,769 single nucleotide polymorphisms and 157,279 INDEs were revealed from transcriptome data. This result will contribute to knowledge of ovary differentiation of Yellow River carp.

Transgelin-2: A potential oncogenic factor

Actin-binding proteins are proteins that could bind to actin or actin fibers. As a member of actin-binding proteins, Transgelin-2 is expressed in smooth muscle cells and non-smooth muscle cells, and its gene, TAGLN2, is differently expressed in all cells and tissues. The deregulation of Transgelin-2 is considered to be correlated with progression of many kinds of diseases, especially the development of malignant tumors, such as invasion, metastasis, and resistance, yet the function and mechanism of action of Transgelin-2 remain elusive. Therefore, we reviewed the basic characteristics and function of Transgelin-2 and its biological role in various types of diseases in order to provide the theoretical basis for further research and new perspectives on cancer development.

Fine selection of up-regulated genes during ovulation by in vivo induction of oocyte maturation and ovulation in zebrafish

BACKGROUND

Two essential processes, oocyte maturation and ovulation, are independently induced, but proceed cooperatively as the final step in oogenesis before oocytes become fertilizable. Although these two processes are induced by the same maturation-inducing steroid, 17α, 20β-dihydroxy-4-pregnen-3-one (17, 20β-DHP), in zebrafish, it has been suggested that the receptor, and thus the signal transduction pathway is different for each process. Although much progress has been made in understanding the molecular mechanisms underlying the induction of oocyte maturation, the mechanisms for inducing ovulation remain under investigation. In the present study, in vivo induction techniques that permit the induction of oocyte maturation and ovulation in living zebrafish (in vivo assays) were used to select highly up-regulated genes (genes associated with ovulation). Using an in vivo assay, ovarian tissues that induced only oocyte maturation could be obtained. This made it possible for the first time to distinguish maturation-inducing genes from ovulation-inducing genes. Using a genome-wide microarray of zebrafish sequences, the gene expression levels were compared among an ethanol (EtOH)-treated group (non-activated group), a diethylstilbestrol (DES)- or testosterone (Tes)-treated group (maturation-induced group), and a 17, 20β-DHP-treated group (maturation- and ovulation-induced group). Ovulation-specific up-regulated genes were selected. The mRNA expression levels of the selected genes were measured by quantitative polymerase chain reaction (qPCR).

RESULTS

Among 34 genes identified, three that showed ovulation-specific increases were selected as candidates potentially associated with ovulation. The ovulation-specific up-regulation of three candidates, slc37a4a, zgc:65811 and zgc:92184 was confirmed by qPCR.

CONCLUSION

Our in vivo assay provides a new approach to precisely select genes associated with ovulation.

Reproduction impairment and endocrine disruption in adult zebrafish (Danio rerio) after waterborne exposure to TBOEP

Tris (2-butoxyethyl) phosphate (TBOEP) is widely used as a substitute of polybrominated diphenyl ethers (PBDEs). It has been frequently measured at concentrations of micrograms per liter (μg/L) in surface waters and waste water. However, limited information is available about the reproduction toxicology of TBOEP. In this study, adult zebrafish pairs were exposed to TBOEP at concentrations of 0, 5, 50, and 500μg/L for 21days. The effects on reproduction, hormone concentration, transcription of genes along the hypothalamic-pituitary-gonadal (HPG) axis, and gonadal development were investigated. After exposure to TBOEP, plasma concentrations of 17β-estradiol were significantly increased in both sexes of fish, while increase of testosterone was observed only in male fish. Transcription of genes along the HPG axis was significantly influenced by exposure to TBOEP in both male and female fish. Moreover, TBOEP decreases the average number of eggs production, as well as hatching success and survival rates in offspring. Histological examination shows inhibition of oocyte maturation in females and retardation spermiation in males, respectively. The results demonstrate that TBOEP could disturb the sex hormone balance by altering regulatory circuits of the HPG axis, affect gonadal development, eventually leading to disruption of reproductive performance and the development of progeny.

Characterizing Transcriptional Networks in Male Rainbow Darter (Etheostoma caeruleum) that Regulate Testis Development over a Complete Reproductive Cycle

Intersex is a condition that has been associated with exposure to sewage effluents in male rainbow darter (Etheostoma caeruleum). To better understand changes in the transcriptome that are associated with intersex, we characterized annual changes in the testis transcriptome in wild, unexposed fish. Rainbow darter males were collected from the Grand River (Ontario, Canada) in May (spawning), August (post-spawning), October (recrudescence), January (developing) and March (pre-spawning). Histology was used to determine the proportion of spermatogenic cell types that were present during each period of testicular maturation. Regression analysis determined that the proportion of spermatozoa versus spermatocytes in all stages of development (R² ≥ 0.58) were inversely related; however this was not the case when males were in the post-spawning period. Gene networks that were specific to the transition from developing to pre-spawning stages included nitric oxide biosynthesis, response to wounding, sperm cell function, and stem cell maintenance. The pre-spawning to spawning transition included gene networks related to amino acid import, glycogenesis, Sertoli cell proliferation, sperm capacitation, and sperm motility. The spawning to post-spawning transition included unique gene networks associated with chromosome condensation, ribosome biogenesis and assembly, and mitotic spindle assembly. Lastly, the transition from post-spawning to recrudescence included gene networks associated with egg activation, epithelial to mesenchymal transition, membrane fluidity, and sperm cell adhesion. Noteworthy was that there were a significant number of gene networks related to immune system function that were differentially expressed throughout reproduction, suggesting that immune network signalling has a prominent role in the male testis. Transcripts in the testis of post-spawning individuals showed patterns of expression that were most different for the majority of transcripts investigated when compared to the other stages. Interestingly, many transcripts associated with female sex differentiation (i.e. esr1, sox9, cdca8 and survivin) were significantly higher in the testis during the post-spawning season compared to other testis stages. At post-spawning, there were higher levels of estrogen and androgen receptors (esr1, esr2, ar) in the testis, while there was a decrease in the levels of sperm associated antigen 1 (spag1) and spermatogenesis associated 4 (spata4) mRNA. Cyp17a was more abundant in the testis of fish in the pre-spawning, spawning, and post-spawning seasons compared to those individuals that were recrudescent while aromatase (cyp19a) did not vary in expression over the year. This study identifies cell process related to testis development in a seasonally spawning species and improves our understanding regarding the molecular signaling events that underlie testicular growth. This is significant because, while there are a number of studies characterizing molecular pathways in the ovary, there are comparatively less describing transcriptomic patterns in the testis in wild fish.

Development and recovery of histopathological alterations in the gonads of zebrafish (Danio rerio) after single and combined exposure to endocrine disruptors (17α-ethinylestradiol and fadrozole)

Exposure of wildlife to endocrine disrupting chemicals (EDCs) is not necessarily continuous. Due to seasonal changes and variable industrial and agricultural activities it often occurs intermittently. Thus, it is possible that aquatic organisms may be more affected by periodic peak exposure than by chronic exposure. Therefore, an experimental scenario including an exposure from 2h to 90 days post-fertilization (dpf) and a subsequent recovery period until 150 dpf was chosen to assess the potential reversibility of the effects of sex steroids on sexual and gonad development of zebrafish (Danio rerio). The aim of this study was to investigate the persistence of the endocrine effects of an estrogen (EE2-17α-ethinylestradiol, 4ng/L), an inhibitor of estrogen synthesis (Fad-fadrozole, 50μg/L) or their binary mixture (Mix-EE2+ Fad, 4ng/L+50μg/L). Afterwards, a semi-quantitative histological assessment was used to investigate histopathological changes on gonad differentiation and development. The data showed that fadrozole, alone or in combination with EE2, permanently disrupts the sexual development, inducing masculinization and causing severe pathological alterations in testis, such as intersex associated to the enlargement of sperm ducts, interstitial changes, asynchronous development and detachment of basal membrane. After exposures to both EDCs and their mixture, the gonad histopathology revealed interstitial proteinaceous fluid deposits and, in ovaries, there were atretic oocytes, and presumably degenerative mineralization. On the other hand, the gonadal changes induced by EE2 alone seem to be partially reversible when the exposure regime changed to a recovery period. In addition, EE2 enhanced zebrafish growth in both genders, with male fish presenting signs of early obesity such as the presence of adipocytes in testis. Moreover, sex ratio was slightly skewed toward females, at 90 and 105 dpf, in zebrafish exposed to EE2. The data further indicate that long-term studies on impacts of single EDCs and their mixtures with recovery periods are crucial to reveal the possibility of sex reversal and pathological changes of gonads that can adversely affect breeding.

Defining the role of omics in assessing ecosystem health: Perspectives from the Canadian environmental monitoring program

Scientific reviews and studies continue to describe omics technologies as the next generation of tools for environmental monitoring, while cautioning that there are limitations and obstacles to overcome. However, omics has not yet transitioned into national environmental monitoring programs designed to assess ecosystem health. Using the example of the Canadian Environmental Effects Monitoring (EEM) program, the authors describe the steps that would be required for omics technologies to be included in such an established program. These steps include baseline collection of omics endpoints across different species and sites to generate a range of what is biologically normal within a particular ecosystem. Natural individual variability in the omes is not adequately characterized and is often not measured in the field, but is a key component to an environmental monitoring program, to determine the critical effect size or action threshold for management. Omics endpoints must develop a level of standardization, consistency, and rigor that will allow interpretation of the relevance of changes across broader scales. To date, population-level consequences of routinely measured endpoints such as reduced gonad size or intersex in fish is not entirely clear, and the significance of genome-wide molecular, proteome, or metabolic changes on organism or population health is further removed from the levels of ecological change traditionally managed. The present review is not intended to dismiss the idea that omics will play a future role in large-scale environmental monitoring studies, but rather outlines the necessary actions for its inclusion in regulatory monitoring programs focused on assessing ecosystem health.

Steroidogenesis and phase II conjugation during the gametogenesis of thicklip grey mullet (Chelon labrosus) from a population showing intersex condition

Steroidogenesis, the process by which steroid hormones are synthesized, involves a vast number of enzymes and biochemical pathways that are susceptible to chemical modulation. Endocrine disrupting compounds (EDCs) are of special concern since they can alter hormone homeostasis by interfering with synthesis, transport and elimination of hormones. It is important to understand gender differences and the natural variation in steroid balance through gamete development in fish exposed to EDCs. The aim of this study was to determine mRNA levels of genes encoding for Steroidogenic Acute Regulatory (star) protein; the steroidogenic enzymes P450 11β hydroxylase (cyp11b1) and P450 aromatase (cyp19a1a); as well as the phase II conjugation enzymes sulfotransferase (sult) and UDP-glucuronosyltransferase (ugt), together with the activity of P450 aromatase and plasma levels of 17β-estradiol (E2) and 11-ketotestosterone (11-KT), at different gametogenic stages and in intersex individuals of the thicklip grey mullet Chelon labrosus. Results demonstrated that the transcription levels of star, sult and ugt and levels of E2 and 11-KT in plasma significantly changed with the interaction between gender and reproductive stage. Cyp11b1 and cyp19a1a transcription levels were significantly different between genders while the activity of P450 aromatase varied significantly between genders and reproductive stages. Results from a multivariate assessment demonstrated that measured endpoints distinguished male, female and intersex mullets at immature gametogenic stage. Intersex distinction was based on sult, ugt and cyp19a1a transcript levels and P450 aromatase activity. The present work provides data to be used in future experimental designs with C. labrosus species, and gives new clues about the molecular events that lead to intersex occurrence in mullets.

The role of eicosanoids in 17α, 20β-dihydroxy-4-pregnen-3-one-induced ovulation and spawning in Danio rerio

This study employed a hormone bioassay to characterize the eicosanoids involved in zebrafish ovulation and spawning, in particular the prostaglandin (PG) products of cyclooxygenase (COX) metabolism and the leukotriene (LT) products of lipoxygenase (LOX) metabolism. Exposure to the teleost progestogen 17α, 20β-dihydroxy-4-pregnen-3-one (17,20βP) induced ovulation, but not spawning, in solitary females and both ovulation and spawning in male-female pairs. Transcription of the eicosanoid-synthesizing enzymes cytosolic phospholipase A2 (cPLA(2)) and COX-2 increased and LTC(4) synthase decreased in peri-ovulatory ovaries of 17,20βP-exposed fish. Ovarian PGF(2α) levels increased post-spawning in 17,20βP-exposed fish, but there was no difference in LTB(4) or LTC(4). Pre-exposure to cPLA(2) or LOX inhibitors reduced 17,20βP-induced ovulation rates, while a COX inhibitor had no effect on ovulation or spawning. Collectively, these findings suggest that eicosanoids, in particular LOX metabolites, mediate 17,20βP-induced ovulation in zebrafish. COX metabolites also appear to be involved in ovulation and spawning but their role remains undefined.

Effects of progesterone and norethindrone on female fathead minnow (Pimephales promelas) steroidogenesis

As knowledge of contaminants capable of adversely modulating endocrine functions increases, attention is focused on the effects of synthetic progestins as environmental endocrine disrupters. In the present study, effects of exposure to a synthetic progestin (norethindrone, 168 ± 7.5 ng/L) and endogenous progestogen (progesterone, 34 ± 4.1 ng/L) on steroidogenesis in adult female fathead minnows were examined. In vivo exposure to either compound lowered expression (nonsignificant) of luteinizing hormone (LHβ) levels in the brain along with significantly down-regulating the beta isoform of membrane progesterone receptor (mPRβ) in ovary tissue. The correspondence between lowered LHβ levels in the brain and mPRβ in the ovary is suggestive of a possible functional association as positive correlations between LHβ and mPR levels have been demonstrated in other fish species. In vitro exposure of ovary tissue to progesterone resulted in significantly elevated progestogen (pregnenolone, 17α-hydroxyprogesterone, and 17α,20β-dihydroxypregnenone) and androgen (testosterone) production. Whereas in vitro exposure to norethindrone did not significantly impact steroid hormone production but showed decreased testosterone production relative to solvent control (however this was not significant). Overall, this study showed that exposure to a natural progestogen (progesterone) and synthetic progestin (norethindrone), was capable of modulating LHβ (in brain) and mPRβ expression (in ovary).

Natural Variation in Fish Transcriptomes: Comparative Analysis of the Fathead Minnow (Pimephales promelas) and Zebrafish (Danio rerio)

Fathead minnow and zebrafish are among the most intensively studied fish species in environmental toxicogenomics. To aid the assessment and interpretation of subtle transcriptomic effects from treatment conditions of interest, better characterization and understanding are needed for natural variation in gene expression among fish individuals from lab cultures. Leveraging the transcriptomics data from a number of our toxicogenomics studies conducted over the years, we conducted a meta-analysis of nearly 600 microarrays generated from the ovary tissue of untreated, reproductively mature fathead minnow and zebrafish samples. As expected, there was considerable batch-to-batch transcriptomic variation; this “batch-effect” appeared to differentially impact subsets of fish transcriptomes in a nonsystematic way. Temporally more closely spaced batches tended to share a greater transcriptomic similarity among one another. The overall level of within-batch variation was quite low in fish ovary tissue, making it a suitable system for studying chemical stressors with subtle biological effects. The observed differences in the within-batch variability of gene expression, at the levels of both individual genes and pathways, were probably both technical and biological. This suggests that biological interpretation and prioritization of genes and pathways targeted by experimental conditions should take into account both their intrinsic variability and the size of induced transcriptional changes. There was significant conservation of both the genomes and transcriptomes between fathead minnow and zebrafish. The high degree of conservation offers promising opportunities in not only studying fish molecular responses to environmental stressors by a comparative biology approach, but also effective sharing of a large amount of existing public transcriptomics data for developing toxicogenomics applications.

Integrated approach to explore the mechanisms of aromatase inhibition and recovery in fathead minnows (Pimephales promelas)

Aromatase, a member of the cytochrome P450 superfamily, is a key enzyme in estradiol synthesis that catalyzes the aromatization of androgens into estrogens in ovaries. Here, we used an integrated approach to assess the mechanistic basis of the direct effects of aromatase inhibition, as well as adaptation and recovery processes in fish. We exposed female fathead minnows (Pimephales promelas) via the water to 30 μg/L of a model aromatase inhibitor, fadrozole, during 8 days (exposure phase). Fish were then held in clean water for 8 more days (recovery phase). Samples were collected at 1, 2, 4, and 8 days of both the exposure and the recovery phases. Transcriptomics, metabolomics, and network inference were used to understand changes and infer connections at the transcript and metabolite level in the ovary. Apical endpoints directly indicative of endocrine function, such as plasma estradiol, testosterone, and vitellogenin levels were also measured. An integrated analysis of the data revealed changes in gene expression consistent with increased testosterone in fadrozole-exposed ovaries. Metabolites such as glycogen and taurine were strongly correlated with increased testosterone levels. Comparison of in vivo and ex vivo steroidogenesis data suggested the accumulation of steroidogenic enzymes, including aromatase, as a mechanism to compensate for aromatase inhibition.

Endocrine disruption and reproduction impairment in zebrafish after long-term exposure to DE-71

The objective of the present study was to investigate the impact of polybrominated diphenyl ethers (PBDEs) on fish reproduction over 2 generations. Zebrafish (Danio rerio) embryos (F0) were exposed to low concentrations (3 µg/L, 10 µg/L, and 30 µg/L) of the PBDE mixture DE-71 until they were sexually mature, and steroid hormone production, expression of genes involved in steroidogenesis, gonadal development, and gamete characteristics were examined. Exposure of female zebrafish to DE-71 resulted in lower estradiol production and downregulation of cytochrome P450 aromatase mRNA. In males, exposure to DE-71 resulted in greater testosterone production and greater cytochrome P450 c17 α-hydroxylase,17,20-lase mRNA expression. Moreover, hepatic vitellogenin mRNA and estrogenic receptor β gene transcription were downregulated in females and males. Expression of the follicle-stimulating hormone β gene in the pituitary was upregulated, and the expression of luteinizing hormone β was downregulated in both sexes. Histological examination showed inhibition of oocyte maturation in females and retarded spermiation in males. The average number of eggs (F1) produced was also reduced. Additionally, exposure of F0 embryos to DE-71 did not result in developmental toxicity, whereas delayed hatching, reduced survival, and decreased growth were observed in the F1 embryos derived from parent fish exposed to DE-71. Therefore, long-term exposure to low concentrations of PBDEs in zebrafish could cause reproductive impairment, suggesting that PBDEs might have significant adverse effects on fish population in the highly PBDEs-contaminated aquatic environment.

Transcriptional signature of progesterone in the fathead minnow ovary (Pimephales promelas)

A growing number of studies have examined transcriptional responses to sex steroids along the hypothalamic-pituitary-gonadal axis in teleost fishes. However, data are lacking on the molecular cascades that underlie progesterone signaling. The objective of this study was to characterize the transcriptional response in the ovary of fathead minnows (Pimephales promelas) in response to progesterone (P4). Fathead minnow ovaries were exposed in vitro to 500 ng P4/L. Germinal vesicle migration and breakdown (GVBD) was observed and microarrays were used to identify gene cascades affected by P4. Microarray analysis identified 1702 differentially expressed transcripts after P4 treatment. Functional enrichment analysis revealed that transcripts involved in the molecular functions of protein serine/threonine kinase activity, ATP binding, and activity of calcium channels were increased after P4 treatment. There was an overwhelming decrease in levels of transcripts of genes that are structural constituents of ribosomes with P4 treatment. There was also evidence for gene expression changes in steroid and maturation-related transcripts. Pathway analyses identified cell cycle regulation, insulin action, hedgehog, and B cell activation as pathways containing an over-representation of highly regulated transcripts. Significant regulatory sub-networks of P4-mediated transcripts included genes regulated by tumor protein p53 and E2F transcription factor 1. These data provide novel insight into the molecular signaling cascades that underlie P4-signaling in the ovary and identify genes and processes that may indicate premature GVBD due to environmental pollutants that mimic progestins.

Assessing gene network stability and individual variability in the fathead minnow (Pimephales promelas) transcriptome

Transcriptomics is increasingly used to assess biological responses to environmental stimuli and stressors such as aquatic pollutants. However, fundamental studies characterizing individual variability in mRNA levels are lacking, which currently limits the use of transcriptomics in environmental monitoring assessments. To address individual variability in transcript abundance, we performed a meta-analysis on 231 microarrays that were conducted in the fathead minnow (FHM), a widely used toxicological model. The mean variability for gene probes was ranked from most to least variable based upon the coefficient of variation. Transcripts that were the most variable in individual tissues included NADH dehydrogenase flavoprotein 1, GTPase IMAP family member 7-like and v-set domain-containing T-cell activation inhibitor 1-like while genes encoding ribosomal proteins (rpl24 and rpl36), basic transcription factor 3, and nascent polypeptide-associated complex alpha subunit were the least variable in individuals across a range of microarray experiments. Gene networks that showed high variability (based upon the variation in expression of individual members within the network) included cell proliferation, metabolism (steroid, lipids, and glucose), cell adhesion, vascularization, and regeneration while those that showed low variability (more stability) included mRNA and rRNA processing, regulation of translational fidelity, RNA splicing, and ribosome biogenesis. Real-time PCR was conducted on a subset of genes for comparison of variability collected from the microarrays. There was a significant positive relationship between the two methods when measuring individual variability, suggesting that variability detected in microarray data can be used to guide decisions on sample sizes for measuring transcripts in real-time PCR experiments. A power analysis revealed that measuring estrogen receptor ba (esrba) requires fewer biological replicates than that of estrogen receptor bb (esrbb) in the gonad and samples sizes required to detect a 50% change for reproductive-related transcripts is between 12 and 20. Characterizing individual variability at the molecular level will prove necessary as efforts are made toward integrating molecular tools into environmental risk assessments.

Gene expression networks underlying ovarian development in wild largemouth bass (Micropterus salmoides)

Background

Oocyte maturation in fish involves numerous cell signaling cascades that are activated or inhibited during specific stages of oocyte development. The objectives of this study were to characterize molecular pathways and temporal gene expression patterns throughout a complete breeding cycle in wild female largemouth bass to improve understanding of the molecular sequence of events underlying oocyte maturation.

Methods

Transcriptomic analysis was performed on eight morphologically diverse stages of the ovary, including primary and secondary stages of oocyte growth, ovulation, and atresia. Ovary histology, plasma vitellogenin, 17β-estradiol, and testosterone were also measured to correlate with gene networks.

Results

Global expression patterns revealed dramatic differences across ovarian development, with 552 and 2070 genes being differentially expressed during both ovulation and atresia respectively. Gene set enrichment analysis (GSEA) revealed that early primary stages of oocyte growth involved increases in expression of genes involved in pathways of B-cell and T-cell receptor-mediated signaling cascades and fibronectin regulation. These pathways as well as pathways that included adrenergic receptor signaling, sphingolipid metabolism and natural killer cell activation were down-regulated at ovulation. At atresia, down-regulated pathways included gap junction and actin cytoskeleton regulation, gonadotrope and mast cell activation, and vasopressin receptor signaling and up-regulated pathways included oxidative phosphorylation and reactive oxygen species metabolism. Expression targets for luteinizing hormone signaling were low during vitellogenesis but increased 150% at ovulation. Other networks found to play a significant role in oocyte maturation included those with genes regulated by members of the TGF-beta superfamily (activins, inhibins, bone morphogenic protein 7 and growth differentiation factor 9), neuregulin 1, retinoid X receptor, and nerve growth factor family.

Conclusions

This study offers novel insight into the gene networks underlying vitellogenesis, ovulation and atresia and generates new hypotheses about the cellular pathways regulating oocyte maturation.

Identification of ovarian gene expression patterns during vitellogenesis in Atlantic cod (Gadus morhua)

Follicular maturational competence and ovulatory competence in teleost fish refer to the ability of the ovarian follicle to undergo final oocyte maturation and ovulation, respectively, in response to gonadotropin stimulation and other external cues. Some gene products related to competence acquisition are likely synthesized during vitellogenic growth, as these follicles gain in vivo responsiveness to exogenous gonadotropin stimulation and can be induced to undergo maturation and ovulation. In Atlantic cod (Gadus morhua), gonadotropin responsiveness has been shown to be oocyte size-dependent, and only ovaries containing late-stage vitellogenic follicles can be induced to ovulate. The purpose of the present study was to compare gene expression patterns between mid (unresponsive) and late (responsive) vitellogenic ovaries to identify genes involved in gonadotropin responsiveness and the acquisition of maturational and ovulatory competencies. Representational difference analysis was conducted in two reciprocal comparisons using intact ovarian fragments and follicle wall-enriched tissues, and genes of interest were used in real time quantitative PCR to confirm differential expression. Few differences were detected in intact ovarian fragments, but type IV ice-structuring protein and gephyrin were upregulated later in development and may be involved in lipid and sulfur metabolism, respectively. Candidate gene assays for luteinizing hormone receptor and aromatase also exhibited significant upregulation during vitellogenesis. Many genes were differentially expressed in follicle wall-enriched tissues, including endocrine maturational regulators and smooth muscle genes. Overall, maturational and ovulatory competencies during vitellogenesis in Atlantic cod are associated with up- and downregulation of many genes involved in lipid metabolism, endocrine regulation, and ovulatory preparation.

Quantitative proteomics in teleost fish: insights and challenges for neuroendocrine and neurotoxicology research

Neuroendocrine systems integrate both extrinsic and intrinsic signals to regulate virtually all aspects of an animal's physiology. In aquatic toxicology, studies have shown that pollutants are capable of disrupting the neuroendocrine system of teleost fish, and many chemicals found in the environment can also have a neurotoxic mode of action. Omics approaches are now used to better understand cell signaling cascades underlying fish neurophysiology and the control of pituitary hormone release, in addition to identifying adverse effects of pollutants in the teleostean central nervous system. For example, both high throughput genomics and proteomic investigations of molecular signaling cascades for both neurotransmitter and nuclear receptor agonists/antagonists have been reported. This review highlights recent studies that have utilized quantitative proteomics methods such as 2D differential in-gel electrophoresis (DIGE) and isobaric tagging for relative and absolute quantitation (iTRAQ) in neuroendocrine regions and uses these examples to demonstrate the challenges of using proteomics in neuroendocrinology and neurotoxicology research. To begin to characterize the teleost neuroproteome, we functionally annotated 623 unique proteins found in the fish hypothalamus and telencephalon. These proteins have roles in biological processes that include synaptic transmission, ATP production, receptor activity, cell structure and integrity, and stress responses. The biological processes most represented by proteins detected in the teleost neuroendocrine brain included transport (8.4%), metabolic process (5.5%), and glycolysis (4.8%). We provide an example of using sub-network enrichment analysis (SNEA) to identify protein networks in the fish hypothalamus in response to dopamine receptor signaling. Dopamine signaling altered the abundance of proteins that are binding partners of microfilaments, integrins, and intermediate filaments, consistent with data suggesting dopaminergic regulation of neuronal stability and structure. Lastly, for fish neuroendocrine studies using both high-throughput genomics and proteomics, we compare gene and protein relationships in the hypothalamus and demonstrate that correlation is often poor for single time point experiments. These studies highlight the need for additional time course analyses to better understand gene-protein relationships and adverse outcome pathways. This is important if both transcriptomics and proteomics are to be used together to investigate neuroendocrine signaling pathways or as bio-monitoring tools in ecotoxicology.

Applications for next-generation sequencing in fish ecotoxicogenomics

The new technologies for next-generation sequencing (NGS) and global gene expression analyses that are widely used in molecular medicine are increasingly applied to the field of fish biology. This has facilitated new directions to address research areas that could not be previously considered due to the lack of molecular information for ecologically relevant species. Over the past decade, the cost of NGS has decreased significantly, making it possible to use non-model fish species to investigate emerging environmental issues. NGS technologies have permitted researchers to obtain large amounts of raw data in short periods of time. There have also been significant improvements in bioinformatics to assemble the sequences and annotate the genes, thus facilitating the management of these large datasets.The combination of DNA sequencing and bioinformatics has improved our abilities to design custom microarrays and study the genome and transcriptome of a wide variety of organisms. Despite the promising results obtained using these techniques in fish studies, NGS technologies are currently underused in ecotoxicogenomics and few studies have employed these methods. These issues should be addressed in order to exploit the full potential of NGS in ecotoxicological studies and expand our understanding of the biology of non-model organisms.

Evaluation of water quality threats to the endangered Okaloosa darter (Etheostoma okaloosae) in East Turkey Creek on Eglin Air Force Base

The threatened Okaloosa darter (Etheostoma okaloosae) is found almost exclusively on the Eglin Air Force Base in the Choctawhatchee Bay watershed of Florida. Portions of this limited habitat are threatened with soil erosion, altered hydrology, and impaired water quality. In the present study, general water quality parameters (i.e., dissolved oxygen, specific conductance, pH, temperature, relative turbidity, and primary productivity) were characterized in East Turkey Creek, which is a body of water potentially impacted by treated wastewater sprayfields, and Long Creek, an adjacent reference stream that does not border the sprayfields. Water quality was assessed during a 30-day exposure using passive samplers for both non-polar and polar effluent parameters. Because the Okaloosa darter was listed as endangered at the time of sampling we chose a closely related species from the same creeks, the sailfin shiner (Pteronotropis hypseleotris) in which to measure metal body burdens. Additionally, fathead minnows (Pimephales promelas) were used for microarray analysis on gonad and liver tissues after 48 h exposures to water collected from the two creeks and brought into the laboratory. Waters from all sites, including reference sites, affected the expression of genes related to various biological processes including transcription and translation, cell cycle control, metabolism, and signaling pathways, suggesting that the sum of anthropogenic compounds in the site waters may cause a generalized stress response in both liver and testis, an effect that could be related to the generally low populations of the Okaloosa darter. Furthermore, effects of site waters on fish gene expression may be related to the impact of human activities other than the wastewater sprayfields, as nearby areas are closed to the public for military testing, training, and administrative activities and due to ordnance contamination.

A graphical systems model and tissue-specific functional gene sets to aid transcriptomic analysis of chemical impacts on the female teleost reproductive axis

Oligonucleotide microarrays and other 'omics' approaches are powerful tools for unsupervised analysis of chemical impacts on biological systems. However, the lack of well annotated biological pathways for many aquatic organisms, including fish, and the limited power of microarray-based analyses to detect low level differential expression of individual genes can hinder the ability to infer and understand chemical effects based on transcriptomic data. Here we report on the supervised assembly of a series of tissue-specific functional gene sets intended to aid transcriptomic analysis of chemical impacts on the female teleost reproductive axis. Gene sets were defined based on an updated graphical systems model of the teleost brain-pituitary-gonadal-hepatic axis. Features depicted in the model were organized into gene sets and mapped to specific probes on three zebrafish (Danio rerio) and two fathead minnow (Pimephales promelas) microarray platforms. Coverage of target genes on the microarrays ranged from 48% for the fathead minnow arrays to 88% for the most current zebrafish platform. Additionally, extended fathead minnow gene sets, incorporating first degree neighbors identified from a Spearman correlation network derived from a large compendium of fathead minnow microarray data, were constructed. Overall, only 14% of the 78 genes queried were connected in the network. Among those, over half had less than five neighbors, while two genes, cyclin b1 and zona pellucida glycoprotein 3, had over 100 first degree neighbors, and were neighbors to one another. Gene set enrichment analyses were conducted using microarray data from a zebrafish hypoxia experiment and fathead minnow time-course experiments conducted with three different endocrine-active chemicals. Results of these analyses demonstrate the utility of the approach for supporting biological inference from ecotoxicogenomic data and comparisons across multiple toxicogenomic experiments. The graphical model, gene mapping, and gene sets described are now available to the scientific community as tools to support ecotoxicogenomic research.