Expression signatures for a model androgen and antiandrogen in the fathead minnow (Pimephales promelas) ovary

Identification of Hub Genes and Biomarkers between Hyperandrogen and Normoandrogen Polycystic Ovary Syndrome by Bioinformatics Analysis

BACKGROUND

The common and divergent genetic mechanisms of hyperandrogen (HA) and normoandrogen (NA) polycystic ovary syndrome (PCOS) are currently unknown.

OBJECTIVE

This study aimed to explore the hub genes and potential mechanisms of HA and NA PCOS through bioinformatics analysis.

METHODS

The GSE137684 dataset was downloaded from the Gene Expression Omnibus (GEO) database. The co-expressed genes and differentially expressed genes (DEGs) between HA and NA PCOS samples were functionally annotated by gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses. The protein-protein interaction (PPI) network of the DEGs was constructed and visualized using STRING and Cytoscape, respectively, and the hub genes were screened using the Cytohubba plug-in. The transcription factors (TFs) of these hub genes were identified with the JASPAR database, and the hub gene-TF regulatory network was constructed.

RESULTS

A total of 327 DEGs, including 191 upregulated and 136 downregulated genes, were identified in HA PCOS relative to NA PCOS. Ten hub genes were screened, of which MYC, CAV1, and HGF were mainly enriched in the Proteoglycans in the cancer pathway. In addition, 47 TFs were identified that were found to be involved in the regulation of hub genes.

CONCLUSION

MYC, CAV1, and HGF are potential diagnostic biomarkers and therapeutic targets for HA PCOS.

New Insights Into the Evolution of Corticotropin-Releasing Hormone Family With a Special Focus on Teleosts

Corticotropin-releasing hormone (CRH) was discovered for its role as a brain neurohormone controlling the corticotropic axis in vertebrates. An additional crh gene, crh2, paralog of crh (crh1), and likely resulting from the second round (2R) of vertebrate whole genome duplication (WGD), was identified in a holocephalan chondrichthyan, in basal mammals, various sauropsids and a non-teleost actinopterygian holostean. It was suggested that crh2 has been recurrently lost in some vertebrate groups including teleosts. We further investigated the fate of crh1 and crh2 in vertebrates with a special focus on teleosts. Phylogenetic and synteny analyses showed the presence of duplicated crh1 paralogs, crh1a and crh1b, in most teleosts, resulting from the teleost-specific WGD (3R). Crh1b is conserved in all teleosts studied, while crh1a has been lost independently in some species. Additional crh1 paralogs are present in carps and salmonids, resulting from specific WGD in these lineages. We identified crh2 gene in additional vertebrate groups such as chondrichthyan elasmobranchs, sarcopterygians including dipnoans and amphibians, and basal actinoperygians, Polypteridae and Chondrostei. We also revealed the presence of crh2 in teleosts, including elopomorphs, osteoglossomorphs, clupeiforms, and ostariophysians, while it would have been lost in Euteleostei along with some other groups. To get some insights on the functional evolution of the crh paralogs, we compared their primary and 3D structure, and by qPCR their tissue distribution, in two representative species, the European eel, which possesses three crh paralogs (crh1a, crh1b, crh2), and the Atlantic salmon, which possesses four crh paralogs of the crh1-type. All peptides conserved the structural characteristics of human CRH. Eel crh1b and both salmon crh1b genes were mainly expressed in the brain, supporting the major role of crh1b paralogs in controlling the corticotropic axis in teleosts. In contrast, crh1a paralogs were mainly expressed in peripheral tissues such as muscle and heart, in eel and salmon, reflecting a striking subfunctionalization between crh1a and b paralogs. Eel crh2 was weakly expressed in the brain and peripheral tissues. These results revisit the repertoire of crh in teleosts and highlight functional divergences that may have contributed to the differential conservation of various crh paralogs in teleosts.

Towards regulation of Endocrine Disrupting chemicals (EDCs) in water resources using bioassays - A guide to developing a testing strategy

Endocrine disrupting chemicals (EDCs) are found in every environmental medium and are chemically diverse. Their presence in water resources can negatively impact the health of both human and wildlife. Currently, there are no mandatory screening mandates or regulations for EDC levels in complex water samples globally. Bioassays, which allow quantifying in vivo or in vitro biological effects of chemicals are used commonly to assess acute toxicity in water. The existing OECD framework to identify single-compound EDCs offers a set of bioassays that are validated for the Estrogen-, Androgen-, and Thyroid hormones, and for Steroidogenesis pathways (EATS). In this review, we discussed bioassays that could be potentially used to screen EDCs in water resources, including in vivo and in vitro bioassays using invertebrates, fish, amphibians, and/or mammalians species. Strengths and weaknesses of samples preparation for complex water samples are discussed. We also review how to calculate the Effect-Based Trigger values, which could serve as thresholds to determine if a given water sample poses a risk based on existing quality standards. This work aims to assist governments and regulatory agencies in developing a testing strategy towards regulation of EDCs in water resources worldwide. The main recommendations include 1) opting for internationally validated cell reporter in vitro bioassays to reduce animal use & cost; 2) testing for cell viability (a critical parameter) when using in vitro bioassays; and 3) evaluating the recovery of the water sample preparation method selected. This review also highlights future research avenues for the EDC screening revolution (e.g., 3D tissue culture, transgenic animals, OMICs, and Adverse Outcome Pathways (AOPs)).

Corticotropin-releasing hormone reduces basal estradiol production in zebrafish follicular cells

Corticotropin-releasing hormone (CRH) plays a key regulatory role in coordinating the regulation of endocrine, autonomic nervous, immune, and reproductive systems. Two CRH (CRHα and CRHβ) and their receptors (CRHR1 and CRHR2) had been identified in zebrafish. However, their functions remained uncovered in the ovary of zebrafish. Therefore, this study aimed to determine whether CRH acts directly on the ovary to regulate steroidogenesis in cultured zebrafish follicular cells. Firstly, CRH and its receptors are expressed in the zebrafish ovary. The expression profile of CRHβ fluctuated during ovarian development in zebrafish, and the highest CRHα mRNA levels were observed at the mature follicle. The highest CRHR1 and CRHR2 mRNA levels existed in mid-vitellogenic (MV) and early vitellogenic (EV) stages, respectively. In primary cultured zebrafish follicular cells, both of the CRHα and CRHβ inhibited expression of hsd17b3 mRNA levels and decreased content of estradiol (E2) in the medium. Furthermore, CRH activated p38 MAPK and p38 MAPK inhibitor SB203580 attenuated the phosphorylation of p38 MAPK induced by CRHα. Simultaneously, SB203580 changed the effect of CRH on cyp19a1a expression but not hsd17b1 and hsd17b3. SB203580 alone or combined with CRH inhibited the E2 content. Finally, the CRHR inhibitor α-helical 9-41 also blocked the phosphorylation of p38 MAPK induced by CRHα but did not change the inhibitory effect of CRH on the mRNA expression of the steroidogenic gene and the content of E2 in the culture medium. Taken together, our findings suggest that the anti-steroidogenic effects of CRH may be mediated partly through activation of the p38 MAPK signaling pathway.

Effects of environmental steroid mixtures are regulated by individual steroid receptor signaling

Fish are exposed to steroids of different classes in contaminated waters, but their effects are not sufficiently understood. Here we employed an anti-sense technique using morpholino oligonucleotides to knockdown the glucocorticoid receptors (GRs, GRα and GRβ) and androgen receptor (AR) to investigate their role in physiological and transcriptional responses. To this end, zebrafish embryos were exposed to clobetasol propionate (CLO), androstenedione (A4) and mixtures containing different classes of steroids. CLO caused a decrease of spontaneous muscle contraction and increase of heart rate, as well as transcriptional induction of pepck1, fkbp5, sult2st3 and vitellogenin (vtg1) at 24 and/or 48 h post fertilization (hpf). Knockdown of GRs eliminated these effects, while knockdown of AR decreased the ar transcript but caused no expressional changes, except induction of sult2st3 after exposure to A4 at 24 hpf. Exposure to a mixture of 6 steroids comprising progesterone (P4) and three progestins, cyproterone acetate, dienogest, drospirenone, 17β-estradiol (E2) and CLO caused a significant induction of pepck1, sult2st3, vtg1 and per1a. Knockdown of GRs eliminated the physiological effects and the up-regulation of vtg1, sult2st3, pepck1, fkbp5 and per1a. Thus, as with CLO, responses in mixtures were regulated by GRs independently from the presence of other steroids. Exposure to a mixture comprising A4, CLO, E2 and P4 caused induction of vtg1, cyp19b, sult2st3 and fkbp5. Knockdown of AR had no effect, indicating that regulation of these genes occurred by the GRs and estrogen receptor (ER). Our findings show that in early embryos GRs cause vtg1 and sult2st3 induction in addition to known glucocorticoid target genes. Each steroid receptor regulated its own target genes in steroid mixtures independently from other steroids. However, enhanced expressional induction occurred for vtg1 and fkbp5 in steroid mixtures, indicating an interaction/cross-talk between GRs and ER. These findings have importance for the understanding of molecular effects of steroid mixtures.

Mummichog (Fundulus heteroclitus) are less sensitive to 17α-ethinylestradiol (EE2) than other common model teleosts: A comparative review of reproductive effects

The environmental estrogen 17α-ethinylestradiol (EE₂) will depress or completely inhibit egg production in many common model teleosts at low concentrations (≤0.5 ng/L; Runnalls et al., 2015). This inhibition is not seen in the estuarine killifish, or mummichog (Fundulus heteroclitus), even when exposed to 100 ng/L EE₂. This relative insensitivity to EE₂ exposure indicates species-specific mechanisms for compensating for exogenous estrogenic exposure. This review compares various reproductive responses elicited by EE₂ in mummichog to other common model teleosts, such as zebrafish (Danio rerio) and fathead minnow (Pimephales promelas), identifying key endpoints where mummichog differ from other studied fish. For example, EE₂ accumulates primarily in the liver/gall bladder of mummichog, which is different than zebrafish and fathead minnow in which accumulation is predominantly in the carcass. Despite causing species-specific differences in fecundity, EE₂ has been shown to consistently induce hepatic vitellogenin in males and cause feminization/sex reversal during gonadal differentiation in larval mummichog, similar to other species. In addition, while gonadal steroidogenesis and plasma steroid levels respond to exogenous EE_2, it is generally at higher concentrations than observed in other species. In mummichog, production of 17β-estradiol (E₂) by full grown ovarian follicles remains high; unlike other teleost models where E₂ synthesis decreases as 17α,20β-dihydroxy-4-prenen-3-on levels increase to induce oocyte maturation. New evidence in mummichog indicates some dissimilarity in gonadal steroidogenic gene expression responses compared to gene expression responses in zebrafish and fathead minnow exposed to EE₂. The role of ovarian physiology continues to warrant investigation regarding the tolerance of mummichog to exogenous EE₂ exposure. Here we present a comprehensive review, highlighting key biological differences in response to EE₂ exposure between mummichog and other commonly used model teleosts.

Effects of 17β-trenbolone exposure on sex hormone synthesis and social behaviours in adolescent mice

17β-Trenbolone (17β-TBOH) is an endocrine disruptor that has been widely reported in aquatic organisms. However, little is known about the effect of 17β-TBOH on mammals, particularly on the development of adolescents. Through a series of behavioural experiments, exposure to at 80 μg kg ^-1 d ^-1 and 800 μg kg ^-1 d ^-1 17β-TBOH during puberty (from PND 28 to 56, male mice) increased anxiety-like behaviours. Exposure to the low dose of 80 μg kg ^-1 d ^-1 resulted in a clear social avoidance behaviour in mice. The two doses affected testicular development and endogenous androgen synthesis in male mice. In addition, 17β-TBOH exposure altered the differentiation of oligodendrocytes and the formation of the myelin sheath in the medial prefrontal cortex (mPFC). These results reveal the effects of 17β-TBOH on the behaviours, gonadal and neurodevelopment of adolescent mammals. In addition, the inhibition of the secretion of endogenous hormones and decrease in the formation of the myelin sheath in mPFC may be associated with the 17β-TBOH-induced behavioural changes in mice.

Adverse Outcome Pathway Network-Based Assessment of the Interactive Effects of an Androgen Receptor Agonist and an Aromatase Inhibitor on Fish Endocrine Function

Predictive approaches to assessing the toxicity of contaminant mixtures have been largely limited to chemicals that exert effects through the same biological molecular initiating event. However, by understanding specific pathways through which chemicals exert effects, it may be possible to identify shared "downstream" nodes as the basis for forecasting interactive effects of chemicals with different molecular initiating events. Adverse outcome pathway (AOP) networks conceptually support this type of analysis. We assessed the utility of a simple AOP network for predicting the effects of mixtures of an aromatase inhibitor (fadrozole) and an androgen receptor agonist (17β-trenbolone) on aspects of reproductive endocrine function in female fathead minnows. The fish were exposed to multiple concentrations of fadrozole and 17β-trenbolone individually or in combination for 48 or 96 h. Effects on 2 shared nodes in the AOP network, plasma 17β-estradiol (E2) concentration and vitellogenin (VTG) production (measured as hepatic vtg transcripts) responded as anticipated to fadrozole alone but were minimally impacted by 17β-trenbolone alone. Overall, there were indications that 17β-trenbolone enhanced decreases in E2 and vtg in fadrozole-exposed fish, as anticipated, but the results often were not statistically significant. Failure to consistently observe hypothesized interactions between fadrozole and 17β-trenbolone could be due to several factors, including lack of impact of 17β-trenbolone, inherent biological variability in the endpoints assessed, and/or an incomplete understanding of interactions (including feedback) between different pathways within the hypothalamic-pituitary-gonadal axis. Environ Toxicol Chem 2020;39:913-922. © 2020 SETAC.

Physiological and Transcriptional Effects of Mixtures of Environmental Estrogens, Androgens, Progestins, and Glucocorticoids in Zebrafish

Fishes are exposed to mixtures of different classes of steroids, but ecotoxicological implications are not sufficiently known. Here, we systematically analyze effects of different combinations of steroid mixtures in zebrafish embryos to assess their joint activities on physiology and transcriptional alterations of steroid-specific target genes at 96 and 120 h post fertilization. In binary mixtures of clobetasol propionate (CLO) with estradiol (E2) or androstenedione (A4), each steroid exhibited its own expression profile. This was also the case in mixtures of 5-, 8-, and 13-different classes of steroids in exposure concentrations of 10-10,000 ng/L. The transcriptional expression of most genes in different mixtures was steroid-specific except for genes encoding aromatase (cyp19b), sulfotransferase (sult2st3), and cyp2k22 that were induced by androgens, progestins, and glucocorticoids. Marked alterations occurred for sult2st3 in binary mixtures of CLO + E2 and CLO + A4. Glucocorticoids increased the heart rate and muscle contractions. In mixtures containing estrogens, induction of the cyp19b transcript occurred at 10 ng/L and protc from the anticoagulation system at 100 ng/L. Our study demonstrates that steroids can act independently in mixtures; the sum of individual steroid profiles is expressed. However, some genes, including cyp19b, sult2st3, and cyp2k22, are regulated by several steroids. This joint effect on different pathways may be of concern for fish development.

Specificity of Biogenic Selenium Nanoparticles for Prostate Cancer Therapy With Reduced Risk of Toxicity: An in vitro and in vivo Study

Selenium deficiency is associated with many physiological disorders including the high risk of cancer. The rehabilitation of selenium with different selenium supplements, however, fails due to their low therapeutic index. Therefore, it is advantageous to have a less toxic form of selenium for supplementation with potentially high anticancer activity. Here we show Bacillus licheniformis derived biogenic selenium nanoparticles at a minimal concentration of 2 μg Se/ml induce necroptosis in LNCaP-FGC cells, without affecting the RBC integrity. Real-time gene expression analysis indicated the overexpression of tumor necrotic factor (TNF) and interferon regulatory factor (IRF1) and decreased expression of androgen receptor (AR) and prostate-specific antigen (PSA). Furthermore, histopathological analysis showed the subsequent oral administrations of 10 times higher concentration of these endotoxin free selenium nanoparticles in C3H/HeJ mice (50 mg Se/kg of body weight), induce significantly lower toxicity compared to the L-selenomethionine (5 mg Se/kg). Our study suggested that the biogenic SeNP could emerge as the safest form of selenium supplementation with potent anticancer activity.

A critical review of the environmental occurrence and potential effects in aquatic vertebrates of the potent androgen receptor agonist 17β-trenbolone

Trenbolone acetate is widely used in some parts of the world for its desirable anabolic effects on livestock. Several metabolites of the acetate, including 17β-trenbolone, have been detected at low nanograms per liter concentrations in surface waters associated with animal feedlots. The 17β-trenbolone isomer can affect androgen receptor signaling pathways in various vertebrate species at comparatively low concentrations/doses. The present article provides a comprehensive review and synthesis of the existing literature concerning exposure to and biological effects of 17β-trenbolone, with an emphasis on potential risks to aquatic animals. In vitro studies indicate that, although 17β-trenbolone can activate several nuclear hormone receptors, its highest affinity is for the androgen receptor in all vertebrate taxa examined, including fish. Exposure of fish to nanograms per liter water concentrations of 17β-trenbolone can cause changes in endocrine function in the short term, and adverse apical effects in longer exposures during development and reproduction. Impacts on endocrine function typically are indicative of inappropriate androgen receptor signaling, such as changes in sex steroid metabolism, impacts on gonadal stage, and masculinization of females. Exposure of fish to 17β-trenbolone during sexual differentiation in early development can greatly skew sex ratios, whereas adult exposures can adversely impact fertility and fecundity. To fully assess ecosystem-level risks, additional research is warranted to address uncertainties as to the degree/breadth of environmental exposures and potential population-level effects of 17β-trenbolone in sensitive species. Environ Toxicol Chem 2018;37:2064-2078. Published 2018 Wiley Periodicals Inc. on behalf of SETAC. This article is a US government work and, as such, is in the public domain in the United States of America.

Testosterone-induced modulation of peroxisomal morphology and peroxisome-related gene expression in brown trout (Salmo trutta f. fario) primary hepatocytes

Disruption of androgenic signaling has been linked to possible cross-modulation with other hormone-mediated pathways. Therefore, our objective was to explore effects caused by testosterone - T (1, 10 and 50μM) in peroxisomal signaling of brown trout hepatocytes. To study the underlying paths involved, several co-exposure conditions were tested, with flutamide - F (anti-androgen) and ICI 182,780 - ICI (anti-estrogen). Molecular and morphological approaches were both evaluated. Peroxisome proliferator-activated receptor alpha (PPARα), catalase and urate oxidase were the selected targets for gene expression analysis. The vitellogenin A gene was also included as a biomarker of estrogenicity. Peroxisome relative volumes were estimated by immunofluorescence, and transmission electron microscopy was used for qualitative morphological control. The single exposures of T caused a significant down-regulation of urate oxidase (10 and 50μM) and a general up-regulation of vitellogenin. A significant reduction of peroxisome relative volumes and smaller peroxisome profiles were observed at 50μM. Co-administration of T and ICI reversed the morphological modifications and vitellogenin levels. The simultaneous exposure of T and F caused a significant and concentration-dependent diminishing in vitellogenin expression. Together, the findings suggest that in the tested model, T acted via both androgen and estrogen receptors to shape the peroxisomal related targets.

Estrogen-dependent seasonal adaptations in the immune response of fish

Clinical and experimental evidence shows that estrogens affect immunity in mammals. Less is known about this interaction in the evolutionary older, non-mammalian, vertebrates. Fish form an excellent model to identify evolutionary conserved neuroendocrine-immune interactions: i) they are the earliest vertebrates with fully developed innate and adaptive immunity, ii) immune and endocrine parameters vary with season, and iii) physiology is constantly disrupted by increasing contamination of the aquatic environment. Neuro-immuno-endocrine interactions enable adaption to changing internal and external environment and are based on shared signaling molecules and receptors. The presence of specific estrogen receptors on/in fish leukocytes, implies direct estrogen-mediated immunoregulation. Fish leukocytes most probably are also capable to produce estrogens as they express the cyp19a and cyp19b - genes, encoding aromatase cytochrome P450, the enzyme critical for conversion of C19 steroids to estrogens. Immunoregulatory actions of estrogens, vary among animal species, and also with dose, target cell type, or physiological condition (e.g., infected/non-infected, reproductive status). They moreover are multifaceted. Interestingly, season-dependent changes in immune status correlate with changes in the levels of circulating sex hormones. Whereas E2 circulating in the bloodstream is perhaps the most likely candidate to be the physiological mediator of systemic immune-reproductive trade-offs, leukocyte-derived hormones are hypothesized to be mainly involved in local tuning of the immune response. Contamination of the aquatic environment with estrogenic EDCs may violate the delicate and precise allostatic interactions between the endogenous estrogen system and the immune system. This has negative effects on fish health, but will also affect the physiology of its consumers.

Transcriptional responses in male Japanese medaka exposed to antiandrogens and antiandrogen/androgen mixtures

The occurrence of androgenic endocrine disrupting chemicals (EDCs) in water is thought to be linked to deviation from normal male developmental and reproductive functions in exposed aquatic organisms. Because aquatic environments represent a chemically complex medium, the combined effects of androgenic EDCs require urgent attention. In the present study, the effects of two model androgen receptor (AR) antagonists, flutamide (FLU), and vinclozolin (VIN), were first determined individually in male Japanese medaka using the transcriptional response for genes associated with the hypothalamic-pituitary-gonadal axis. The fish were further exposed to binary mixtures of VIN and 17β-trenbolone (TRE, AR agonist) to confirm the theoretical opposing effects of the AR antagonist and agonist. The results showed that exposure to FLU or VIN alone induced very similar transcriptional responses, demonstrating that gene transcription analysis could be successfully employed in identifying the action of single chemicals. For example, both exposures increased the transcription of cyp17b but decreased that of cyp19b in the gonad, demonstrating the compensatory response for AR blockage. However, in the case of exposure to mixtures, although the joint antagonistic action of TRE and VIN affected the most genes, the transcription profiles after exposure to mixtures were not consistent with expectations based on the results for individual chemicals, such as hepatic vtg, and star or cyp19a in gonads. Therefore, the limitation of gene transcription analyses in exposures to mixtures, as well as the potential for the extrapolation of single chemicals, should be considered in future studies. © 2015 Wiley Periodicals, Inc. Environ Toxicol 31: 1591-1599, 2016.

A study of temporal effects of the model anti-androgen flutamide on components of the hypothalamic-pituitary-gonadal axis in adult fathead minnows

The aim of this study was to investigate temporal changes in the hypothalamic-pituitary-gonadal (HPG) axis of fathead minnows (Pimephales promelas) treated with the model androgen receptor (AR) antagonist flutamide. Reproductively-mature fish were exposed in a flow-through test to analytically-confirmed concentrations of either 50 or 500μg flutamide/L for 8 d, followed by an 8-d recovery period in clean water. Fish were sampled at 1, 2, 4 and 8days during each phase of the experiment. Flutamide (500μg/L) caused significant reductions in relative gonad size of the females on day 8 of the exposure and day 1 of the recovery, and reduced expression of secondary sex characteristics in males during the exposure phase of the experiment. Ex vivo gonadal synthesis of testosterone in both sexes (and 17β-estradiol in females) was reduced in the 500μg/L treatment within 2 d of exposure; however, steroid synthesis returned to levels comparable to controls by the end of the exposure portion of the test. Ex vivo testosterone synthesis in males exposed to 50μg flutamide/L was greater than in controls on days 4 and 8 of the exposure. Both the enhanced steroid production in the low treatment males, and return to control levels in the high treatment males and females during chemical exposure are indicative of a compensatory HPG response. One contributor to this response could be increased expression of genes responsible for enzymes involved in steroid synthesis; for example, transcripts for both cytochrome P450 side- chain cleavage and 11β-hydroxysteroid dehydrogenase were significantly elevated in flutamide-exposed males. Overall, responses of the HPG axis in adult male and female fathead minnows exposed to flutamide were both dynamic and comparatively rapid during exposure and recovery. These observations have ramifications both for the development of short-term fish assays to detect endocrine-active chemicals, and the derivation of robust adverse outcome pathways for AR antagonists in fish.

Fish Reproduction Is Disrupted upon Lifelong Exposure to Environmental PAHs Fractions Revealing Different Modes of Action

Polycyclic aromatic hydrocarbons (PAHs) constitute a large family of organic pollutants emitted in the environment as complex mixtures, the compositions of which depend on origin. Among a wide range of physiological defects, PAHs are suspected to be involved in disruption of reproduction. In an aquatic environment, the trophic route is an important source of chronic exposure to PAHs. Here, we performed trophic exposure of zebrafish to three fractions of different origin, one pyrolytic and two petrogenic. Produced diets contained PAHs at environmental concentrations. Reproductive traits were analyzed at individual, tissue and molecular levels. Reproductive success and cumulative eggs number were disrupted after exposure to all three fractions, albeit to various extents depending on the fraction and concentrations. Histological analyses revealed ovary maturation defects after exposure to all three fractions as well as degeneration after exposure to a pyrolytic fraction. In testis, hypoplasia was observed after exposure to petrogenic fractions. Genes expression analysis in gonads has allowed us to establish common pathways such as endocrine disruption or differentiation/maturation defects. Taken altogether, these results indicate that PAHs can indeed disrupt fish reproduction and that different fractions trigger different pathways resulting in different effects.

Transcriptional networks associated with the immune system are disrupted by organochlorine pesticides in largemouth bass (Micropterus salmoides) ovary

Largemouth bass (Micropterus salmoides) inhabiting Lake Apopka, Florida are exposed to high levels of persistent organochlorine pesticides (OCPs) and dietary uptake is a significant route of exposure for these apex predators. The objectives of this study were to determine the dietary effects of two organochlorine pesticides (p, p'-dichlorodiphenyldichloroethylene; p, p' DDE and methoxychlor; MXC) on the reproductive axis of largemouth bass. Reproductive bass (late vitellogenesis) were fed one of the following diets: control pellets, 125ppm p, p'-DDE, or 10ppm MXC (mg/kg) for 84days. Due to the fact that both p,p' DDE and MXC have anti-androgenic properties, the anti-androgenic pharmaceutical flutamide was fed to a fourth group of largemouth bass (750ppm). Following a 3 month exposure, fish incorporated p,p' DDE and MXC into both muscle and ovary tissue, with the ovary incorporating 3 times more organochlorine pesticides compared to muscle. Endpoints assessed were those related to reproduction due to previous studies demonstrating that these pesticides impact the reproductive axis and we hypothesized that a dietary exposure would result in impaired reproduction. However, oocyte distribution, gonadosomatic index, plasma vitellogenin, and plasma sex steroids (17β-estradiol, E2 and testosterone, T) were not different between control animals and contaminant-fed largemouth bass. Moreover, neither p, p' DDE nor MXC affected E2 or T production in ex vivo oocyte cultures from chemical-fed largemouth bass. However, both pesticides did interfere with the normal upregulation of androgen receptor that is observed in response to human chorionic gonadotropin in ex vivo cultures, an observation that may be related to their anti-androgenic properties. Transcriptomics profiling in the ovary revealed that gene networks related to cell processes such as leukocyte cell adhesion, ossification, platelet function and inhibition, xenobiotic metabolism, fibrinolysis, and thermoregulation were altered by p, p' DDE, MXC, and flutamide. Interestingly, immune-related gene networks were suppressed by all three chemicals. The data suggest that p, p' DDE and flutamide affected more genes in common with each other than either chemical with MXC, consistent with studies suggesting that p, p' DDE is a more potent anti-androgen than MXC. These data demonstrate that reproductive health was not affected by these specific dietary treatments, but rather the immune system, which may be a significant target of organochlorine pesticides. The interaction between the reproductive and immune systems should be considered in future studies on these legacy and persistent pesticides.

Transcriptional networks associated with 5-alpha-dihydrotestosterone in the fathead minnow (Pimephales promelas) ovary

Androgens play a significant role in regulating oogenesis in teleost fishes. The androgen dihydrotestosterone (DHT) is a potent non-aromatizable androgen involved in sexual differentiation in mammals; however, its actions are not well understood in teleost fish. To better characterize the physiological role of DHT in the fathead minnow (FHM) ovary on a temporal scale, in vitro assays for 17β-estradiol (E2) production were conducted in parallel with microarray analysis. Ovarian explants were incubated at different concentrations of DHT (10(-6), 10(-7), and 10(-8)M DHT) in three separate experiments conducted at 6, 9, and 12h. DHT treatment resulted in a rapid and consistent increase in E2 production from the ovary at all three time points. Therefore, DHT may act to shift the balance of metabolites in the steroidogenic pathway within the ovary. Major biological themes affected by DHT in the ovary in one or more of the time points included those related to blood (e.g. vasodilation, blood vessel contraction, clotting), lipids (e.g. lipid storage, cholesterol metabolism, lipid degradation) and reproduction (e.g. hormone and steroid metabolism). Gene networks related to immune responses and calcium signaling were also affected by DHT, suggesting that this androgen may play a role in regulating these processes in the ovary. This study detected no change in mRNA levels of steroidogenic enzymes (cyp19a1, star, 11βhsd, 17βhsd, srd5a isoforms), suggesting that the observed increase in E2 production is likely more dependent on the pre-existing gene or protein complement in the ovary rather than the de novo expression of transcripts. This study increases knowledge regarding the roles of DHT and androgens in general in the teleost ovary and identifies molecular signaling pathways that may be associated with increased E2 production.

Effects of 17 α-methyltestosterone on transcriptome, gonadal histology and sex steroid hormones in rare minnow Gobiocypris rarus

The 17α-methyltestosterone (MT), a synthetic androgen, is known for its interference effects on the endocrine system. Aiming to investigate the transcriptome profiling of gonads induced by MT and to understand the molecular mechanism by which MT causes adverse effects in fish, transcriptome profiling of gonads, gonadal histology and the sex steroid hormones in response to MT were analyzed in Gobiocypris rarus. Eight libraries, 4 from the ovary and 4 from the testis, were constructed and sequenced and then a total number of clean reads per sample ranging from 7.03 to 9.99 million were obtained. In females, a total of 191 transcripts were differentially regulated by MT, consisting of 102 up-regulated transcripts and 89 down-regulated transcripts. In males, 268 differentially expressed genes with 108 up-regulated and 160 down-regulated were detected upon MT exposure. Testosterone serves as the major sex steroid hormone content in G. rarus of both sexes. The concentrations of 17β-estradiol, testosterone and 11-ketotestosterone were significantly increased in females and decreased in males after MT exposure. Interestingly, MT caused a decreased number of vitellogenic oocytes in the ovary and spermatozoa in the testis. After MT exposure, four differentially expressed genes (ndufa4, slc1a3a, caskin-2 and rpt3) were found in G. rarus of both sexes. Overall, we suggest that MT seemed to affect genes involved in pathways related to physiological processes in the gonads of G. rarus. These processes include the electron transfer of Complex IV, endothelial cell activation, axon growth and guidance, and proteasome assembly and glutamate transport metabolic.

Endocrinology: advances through omics and related technologies

The rapid development of new omics technologies to measure changes at genetic, transcriptomic, proteomic, and metabolomics levels together with the evolution of methods to analyze and integrate the data at a systems level are revolutionizing the study of biological processes. Here we discuss how new approaches using omics technologies have expanded our knowledge especially in nontraditional models. Our increasing knowledge of these interactions and evolutionary pathway conservation facilitates the use of nontraditional species, both invertebrate and vertebrate, as new model species for biological and endocrinology research. The increasing availability of technology to create organisms overexpressing key genes in endocrine function allows manipulation of complex regulatory networks such as growth hormone (GH) in transgenic fish where disregulation of GH production to produce larger fish has also permitted exploration of the role that GH plays in testis development, suggesting that it does so through interactions with insulin-like growth factors. The availability of omics tools to monitor changes at nearly any level in any organism, manipulate gene expression and behavior, and integrate data across biological levels, provides novel opportunities to explore endocrine function across many species and understand the complex roles that key genes play in different aspects of the endocrine function.

Transcriptional and physiological responses induced by binary mixtures of drospirenone and progesterone in zebrafish (Danio rerio)

Drospirenone (DRS) is a synthetic progestin increasingly used in oral contraceptives with similar effects to progesterone (P4). Wild fish are exposed to DRS and P4 through wastewater. However, the effects of DRS on fish, both as an individual compound and in mixtures, have not been extensively studied. Therefore, in this study, global gene expression profiles of ovary and brain of female zebrafish (Danio rerio) were characterized after exposure to 55, 553, and 5442 ng/L DRS for 14 days. The effects were then compared to the observed responses after exposure to mixtures of DRS and P4 (DRS+P4: 27 + 0.8, 277 + 8 and 3118 + 123 ng/L). Transcriptomics findings were related to the changes in vitellogenin protein concentrations in the blood, morphology, and histology of gonads. Multivariate analysis indicated tissue-, dose-, and treatment-dependent expression profiles. Genes involved in steroid hormone receptor activity and circadian rhythm were enriched in DRS and mixture groups, among other pathways. In mixtures, the magnitude of response was dose- and transcript-dependent, both at the molecular and physiological levels. Effects of DRS and P4 were additive for most of the investigated parameters and occurred at environmentally relevant concentrations. They may translate to adverse reproductive effects in fish.

Gene expression of fathead minnows (Pimephales promelas) exposed to two types of treated municipal wastewater effluents

Contaminants of emerging concern (CECs) in treated municipal effluents have the potential to adversely impact exposed organisms prompting elevated public concern. Using transcriptomic tools, we investigated changes in gene expression and cellular pathways in the liver of male fathead minnows (Pimephales promelas) exposed to 5% concentrations of full secondary-treated (HTP) or advanced primary-treated (PL) municipal wastewater effluents containing CECs. Gene expression changes were associated with apical end points (plasma vitellogenin and changes in secondary sexual characteristics). Of 32 effluent CECs analyzed, 28 were detected including pharmaceuticals, personal care products, hormones, and industrial compounds. Exposure to both effluents produced significantly higher levels of plasma VTG and changes in secondary sexual characteristics (e.g., ovipositor development). Transcript patterns differed between effluents, with <10% agreement in the detected response (e.g., altered production of transcripts involved in xenobiotic detoxification, oxidative stress, and apoptosis were observed following exposure to both effluents). Exposure to PL effluent caused changes in transcription of genes involved in metabolic pathways (e.g., lipid transport and steroid metabolism). Exposure to HTP effluent affected transcripts involved in signaling pathways (e.g., focal adhesion assembly and extracellular matrix). The results suggest a potential association between some transcriptomic changes and physiological responses following effluent exposure. This study identified responses in pathways not previously implicated in exposure to complex chemical mixtures containing CECs, which are consistent with effluent exposure (e.g., oxidative stress) in addition to other pathway responses specific to the effluent type.

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.

Anti-androgens act jointly in suppressing spiggin concentrations in androgen-primed female three-spined sticklebacks - prediction of combined effects by concentration addition

Increasing attention is being directed at the role played by anti-androgenic chemicals in endocrine disruption of wildlife within the aquatic environment. The co-occurrence of multiple contaminants with anti-androgenic activity highlights a need for the predictive assessment of combined effects, but information about anti-androgen mixture effects on wildlife is lacking. This study evaluated the suitability of the androgenised female stickleback screen (AFSS), in which inhibition of androgen-induced spiggin production provides a quantitative assessment of anti-androgenic activity, for predicting the effect of a four component mixture of anti-androgens. The anti-androgenic activity of four known anti-androgens (vinclozolin, fenitrothion, flutamide, linuron) was evaluated from individual concentration-response data and used to design a mixture containing each chemical at equipotent concentrations. Across a 100-fold concentration range, a concentration addition approach was used to predict the response of fish to the mixture. Two studies were conducted independently at each of two laboratories. By using a novel method to adjust for differences between nominal and measured concentrations, good agreement was obtained between the actual outcome of the mixture exposure and the predicted outcome. This demonstrated for the first time that androgen receptor antagonists act in concert in an additive fashion in fish and that existing mixture methodology is effective in predicting the outcome, based on concentration-response data for individual chemicals. The sensitivity range of the AFSS assay lies within the range of anti-androgenicity reported in rivers across many locations internationally. The approach taken in our study lays the foundations for understanding how androgen receptor antagonists work together in fish and is essential in informing risk assessment methods for complex anti-androgenic mixtures in the aquatic environment.

Hepatic protein expression networks associated with masculinization in the female fathead minnow (Pimephales promelas)

Endocrine disruptors that act via the androgen receptor (AR) are less well studied than environmental estrogens, and there is evidence that treatment with AR agonists can result in masculinization of female fish. In this study, female fathead minnows (FHM) were exposed to the model nonaromatizable androgen 5-alpha dihydrotestosterone (DHT) (100 μg/L), the ureic-based herbicide linuron (LIN) (100 μg/L), and a mixture of DHT and LIN (100 μg/L each) to better characterize androgen action in females. LIN was used because of reports that this chemical has an antiandrogenic mode of action in fish. After 21d, DHT and LIN treatments resulted in a significant depression of plasma vitellogenin (Vtg) and DHT and DHT+LIN increased the prevalence of nuptial tubercles in female FHMs indicating masculinization. Using iTRAQ and an LTQ Orbitrap Velos, ∼2000 proteins were identified in the FHM liver and the number of proteins quantified after exposures was >1200. Proteins that significantly and consistently changed in abundance across biological replicates included prostaglandin E synthase 3, programmed cell death 4a, glutathione S transferases, canopy, selenoprotein U, and ribosomal proteins. Subnetwork enrichment analysis identified that interferon and epidermal growth factor signaling were regulated by DHT and LIN, suggesting that these signaling pathways are correlated to depressed plasma vitellogenin. These data provide novel insight into hepatic protein networks that are associated with the process of masculinization in teleosts.

Exploring androgen-regulated pathways in teleost fish using transcriptomics and proteomics

In the environment, there are aquatic pollutants that disrupt androgen signaling in fish. Laboratory and field-based experiments have utilized omics technologies to characterize the molecular mechanisms underlying androgen-receptor agonism/antagonism. Transcriptomics and proteomics studies with 17β-trenbolone, a growth-promoting pharmaceutical found in water systems surrounding cattle feed lots, and androgens such as 17α-methyltestosterone and 17α-methyldihydrotestosterone, have been conducted in ovary and liver of fish that include the fathead minnow (FHM) (Pimephales promelas), common carp (Cyprinus carpio), Qurt medaka (Oryzias latipes), and zebrafish (Danio rerio). In this mini-review, we survey recent omics studies in fish and reveal that, despite the diversity of species and tissues examined, there are common cellular responses that are observed with waterborne androgenic treatments. Recurring themes in gene ontology include apoptosis, transport and oxidation of lipids, synthesis and transport of hormones, immune response, protein metabolism, and cell proliferation. However, we also discuss other mechanisms other than androgen receptor (AR) activation, such as responses to toxicant stress, estrogen receptor agonism, aromatization of androgens into estrogens, and inhibitory feedback mechanisms by high levels of androgens that may also explain molecular responses in fish. To further explore androgen-responsive protein networks, a sub-network enrichment analysis was performed on protein data collected from the livers of female FHMs exposed to 17β-trenbolone. We construct a putative AR-regulated protein/cell process network in the liver that includes B-lymphocyte differentiation, xenobiotic clearance, low-density lipoprotein oxidation, proliferation of smooth muscle cells, and permeability of blood vessels. We demonstrate that construction of protein networks can offer insight into cell processes that are potentially regulated by androgens.

Taking microarrays to the field: differential hepatic gene expression of caged fathead minnows from Nebraska watersheds

This study aimed to evaluate the utility of microarrays as a biomonitoring tool in field studies. A 15,000-oligonucleotide microarray was used to measure the hepatic gene expression of fathead minnows (Pimephales promelas) caged in four Nebraska, USA watersheds - the Niobrara and Dismal Rivers (low-impact agricultural sites) and the Platte and Elkhorn Rivers (high-impact agricultural sites). Gene expression profiles were site specific and fish from the low- and high-impact sites aggregated into distinct groups. Over 1500 genes were differentially regulated between fish from the low- and high-impact sites. Many gene expression differences (1218) were also noted when the Platte and Elkhorn minnows were compared to one another and Platte fish experienced a higher degree of transcript alterations than Elkhorn fish. These findings indicate that there are differences between the low-impact and high-impact sites, as well as between the two high-impact sites. Historical water quality data support these results as only trace levels of agrichemicals have been detected at the low-impact sites, while substantial levels of agrichemicals have been reported at the high-impact sites with agrichemical loads at the Platte generally exceeding those at the Elkhorn. Overall, this study demonstrates that microarrays can be utilized to discriminate sites with different contaminant loads from one another.

DIGE and iTRAQ as biomarker discovery tools in aquatic toxicology

Molecular approaches in ecotoxicology have greatly enhanced mechanistic understanding of the impact of aquatic pollutants in organisms. These methods have included high throughput Omics technologies, including quantitative proteomics methods such as 2D differential in-gel electrophoresis (DIGE) and isobaric tagging for relative and absolute quantitation (iTRAQ). These methods are becoming more widely used in ecotoxicology studies to identify and characterize protein bioindicators of adverse effect. In teleost fish, iTRAQ has been used successfully in different fish species (e.g. fathead minnow, goldfish, largemouth bass) and tissues (e.g. hypothalamus and liver) to quantify relative protein abundance. Of interest for ecotoxicology is that many proteins commonly utilized as bioindicators of toxicity or stress are quantifiable using iTRAQ on a larger scale, providing a global baseline of biological effect from which to assess changes in the proteome. This review highlights the successes to date for high throughput quantitative proteomics using DIGE and iTRAQ in aquatic toxicology. Current challenges for the iTRAQ method for biomarker discovery in fish are the high cost and the lack of complete annotated genomes for teleosts. However, the use of protein homology from teleost fishes in protein databases and the introduction of hybrid LTQ-FT (Linear ion trap-Fourier transform) mass spectrometers with high resolution, increased sensitivity, and high mass accuracy are able to improve significantly the protein identification rates. Despite these challenges, initial studies utilizing iTRAQ for ecotoxicoproteomics have exceeded expectations and it is anticipated that the use of non-gel based quantitative proteomics will increase for protein biomarker discovery and for characterization of chemical mode of action.

Progress and promises in toxicogenomics in aquatic toxicology: is technical innovation driving scientific innovation?

In the last decade, new technologies have been invented to analyze large amounts of information such as gene transcripts (transcriptomics), proteins (proteomics) and small cellular molecules (metabolomics). Many studies have been performed in the last few years applying these technologies to aquatic toxicology, mainly in fish. In this article, we summarize the current state of knowledge and question whether the application of modern technology for descriptive purposes truly represents scientific advancement in aquatic toxicology. We critically discuss the advantages and disadvantages of these technologies and emphasize the importance of these critical aspects. To date, these techniques have been used mainly as a proof of principle, demonstrating effects of model compounds. The potential to use these techniques to better analyze the mode-of-action of a toxicant or the effects of a compound within organisms has rarely been met. This is partly due to a lack of baseline data and the fact that the expression of mRNA and protein profiles is rarely linked to physiology or toxicologically meaningful outcomes. It seems premature to analyze mixtures or environmental samples until more is known about the expression profiles of individual toxicants. Gene transcription, protein, or metabolic data give only a partial view of these effects. Thus, we emphasize that data obtained by these technologies must be linked to physiological changes to fully understand their significance. The use of these techniques in aquatic toxicology is still in its infancy, data cannot yet be applied to environmental risk assessment or regulation until more emphasis is placed on interpreting the data within their physiological and toxicological contexts.

New insights into placozoan sexual reproduction and development

Unraveling animal life cycles and embryonic development is basic to understanding animal biology and often sheds light on phylogenetic relationships. A key group for understanding the evolution of the Metazoa is the early branching phylum Placozoa, which has attracted rapidly increasing attention. Despite over a hundred years of placozoan research the life cycle of this enigmatic phylum remains unknown. Placozoa are a unique model system for which the nuclear genome was published before the basic biology (i.e. life cycle and development) has been unraveled. Four organismal studies have reported the development of oocytes and one genetic study has nourished the hypothesis of sexual reproduction in natural populations at least in the past. Here we report new observations on sexual reproduction and embryonic development in the Placozoa and support the hypothesis of current sexual reproduction. The regular observation of oocytes and expressed sperm markers provide support that placozoans reproduce sexually in the field. Using whole genome and EST sequences and additional cDNA cloning we identified five conserved sperm markers, characteristic for different stages in spermatogenesis. We also report details on the embryonic development up to a 128-cell stage and new ultrastructural features occurring during early development. These results suggest that sperm and oocyte generation and maturation occur in different placozoans and that clonal lineages reproduce bisexually in addition to the standard mode of vegetative reproduction. The sum of observations is best congruent with the hypothesis of a simple life cycle with an alternation of reproductive modes between bisexual and vegetative reproduction.

A computational model of the hypothalamic: pituitary: gonadal axis in female fathead minnows (Pimephales promelas) exposed to 17α-ethynylestradiol and 17β-trenbolone

Background

Endocrine disrupting chemicals (e.g., estrogens, androgens and their mimics) are known to affect reproduction in fish. 17α-ethynylestradiol is a synthetic estrogen used in birth control pills. 17β-trenbolone is a relatively stable metabolite of trenbolone acetate, a synthetic androgen used as a growth promoter in livestock. Both 17α-ethynylestradiol and 17β-trenbolone have been found in the aquatic environment and affect fish reproduction. In this study, we developed a physiologically-based computational model for female fathead minnows (FHM, Pimephales promelas), a small fish species used in ecotoxicology, to simulate how estrogens (i.e., 17α-ethynylestradiol) or androgens (i.e., 17β-trenbolone) affect reproductive endpoints such as plasma concentrations of steroid hormones (e.g., 17β-estradiol and testosterone) and vitellogenin (a precursor to egg yolk proteins).

Results

Using Markov Chain Monte Carlo simulations, the model was calibrated with data from unexposed, 17α-ethynylestradiol-exposed, and 17β-trenbolone-exposed FHMs. Four Markov chains were simulated, and the chains for each calibrated model parameter (26 in total) converged within 20,000 iterations. With the converged parameter values, we evaluated the model's predictive ability by simulating a variety of independent experimental data. The model predictions agreed with the experimental data well.

Conclusions

The physiologically-based computational model represents the hypothalamic-pituitary-gonadal axis in adult female FHM robustly. The model is useful to estimate how estrogens (e.g., 17α-ethynylestradiol) or androgens (e.g., 17β-trenbolone) affect plasma concentrations of 17β-estradiol, testosterone and vitellogenin, which are important determinants of fecundity in fish.

The effects of estrogenic and androgenic endocrine disruptors on the immune system of fish: a review

During the last decade, a number of studies have shown that, in addition to their classically described reproductive function, estrogens and androgens also regulate the immune system in teleosts. Today, several molecules are known to interfere with the sex-steroid signaling. These chemicals are often referred to as endocrine disrupting contaminants (EDCs). We review the growing evidence that these compounds interfere with the fish immune system. These studies encompass a broad range of approaches from field studies to those at the molecular level. This integrative overview improves our understanding of the various endocrine-disrupting processes triggered by these chemicals. Furthermore, the research also explains why fish that have been exposed to EDCs are more sensitive to pathogens during gametogenesis. In this review, we first discuss the primary actions of sex-steroid-like endocrine disruptors in fish and the specificity of the fish immune system in comparison to mammals. Then, we review the known interactions between the immune system and EDCs and interpret the primary effects of sex steroids (estrogens and androgens) and their related endocrine disruptors on immune modulation. The recent literature suggests that immune parameters may be used as biomarkers of contamination by EDCs. However, caution should be used in the assessment of such immunotoxicity. In particular, more attention should be paid to the specificity of these biomarkers, the external/internal factors influencing the response, and the transduction pathways induced by these molecules in fish. The use of the well-known mammalian models provides a useful guide for future research in fish.

Systems biology: leading the revolution in ecotoxicology

The rapid development of new technologies such as transcriptomics, proteomics, and metabolomics (Omics) are changing the way ecotoxicology is practiced. The data deluge has begun with genomes of over 65 different aquatic species that are currently being sequenced, and many times that number with at least some level of transcriptome sequencing. Integrating these top-down methodologies is an essential task in the field of systems biology. Systems biology is a biology-based interdisciplinary field that focuses on complex interactions in biological systems, with the intent to model and discover emergent properties of the system. Recent studies demonstrate that Omics technologies provide valuable insight into ecotoxicity, both in laboratory exposures with model organisms and with animals exposed in the field. However, these approaches require a context of the whole animal and population to be relevant. Powerful approaches using reverse engineering to determine interacting networks of genes, proteins, or biochemical reactions are uncovering unique responses to toxicants. Modeling efforts in aquatic animals are evolving to interrelate the interacting networks of a system and the flow of information linking these elements. Just as is happening in medicine, systems biology approaches that allow the integration of many different scales of interaction and information are already driving a revolution in understanding the impacts of pollutants on aquatic systems.

Relationships between the transcriptome and physiological indicators of reproduction in female rainbow trout over an annual cycle

Normal transcriptomic patterns along the brain-pituitary-gonad-liver (BPGL) axis should be better characterized if endocrine-disrupting compound-induced changes in gene expression are to be understood. Female rainbow trout were studied over a complete year-long reproductive cycle. Tissue samples from pituitary, ovary, and liver were collected for microarray analysis using the 16K Genomic Research on Atlantic Salmon Project (GRASP) microarray and for quantitative polymerase chain reaction measures of estrogen receptor (ER) isoform messenger RNA (mRNA) levels. Plasma was collected to determine levels of circulating estradiol-17β (E2), follicle-stimulating hormone (FSH), and luteinizing hormone (LH). As an a priori hypothesis, changes in gene expression were correlated to either circulating levels of E2, FSH, and LH, or ER mRNAs quantified by quantitative polymerase chain reaction. In the liver, most transcriptomic patterns correlated to levels of either E2, LH, or ERs. Fewer ovarian transcripts could be correlated to levels of E2, ERα, or FSH. No significant associations were obvious in the pituitary. As a post hoc hypothesis, changes in transcript abundance were compared with microarray features with known roles in gonadal maturation. Many altered transcripts in the ovary correlated to transcript levels of estradiol 17-beta-dehydrogenase 8 or 17 B HSD12, or to glycoprotein alpha chain 1 or 2. In the pituitary, genes involved with the growth axis (e.g., growth hormone, insulin-related growth factor binding protein) correlated with the most transcripts. These results suggest that transcriptional networks along the BPGL axis may be regulated by factors other than circulating steroid hormones.

Investigations of transcript expression in fathead minnow (Pimephales promelas) brain tissue reveal toxicological impacts of RDX exposure

Production, usage and disposal of the munitions constituent (MC) cyclotrimethylenetrinitramine (RDX) has led to environmental releases on military facilities. The chemical attributes of RDX are conducive for leaching to surface water which may put aquatic organisms at risk of exposure. Because RDX has been observed to cause aberrant neuromuscular effects across a wide range of animal phyla, we assessed the effects of RDX on central nervous system (CNS) functions in the representative aquatic ecotoxicological model species, fathead minnow (Pimephales promelas). We developed a fathead minnow brain-tissue cDNA library enriched for transcripts differentially expressed in response to RDX and trinitrotoluene (TNT) exposure. All 4,128 cDNAs were sequenced, quality filtered and assembled yielding 2230 unique sequences and 945 significant blastx matches (E ≤10(-5)). The cDNA library was leveraged to create custom-spotted microarrays for use in transcript expression assays. The impact of RDX on transcript expression in brain tissue was examined in fathead minnows exposed to RDX at 0.625, 2.5, 5, 10mg/L or an acetone-spike control for 10 days. Overt toxicity of RDX in fathead minnow occurred only at the highest exposure concentration resulting in 50% mortality and weight loss. Conversely, Bayesian analysis of microarray data indicated significant changes in transcript expression at concentrations as low as 0.625 mg/L. In total, 154 cDNAs representing 44 unique transcripts were differentially expressed in RDX exposures, the majority of which were validated by reverse transcriptase-quantitative PCR (RT-qPCR). Investigation of molecular pathways, gene ontology (GO) and individual gene functions affected by RDX exposures indicated changes in metabolic processes involved in: oxygen transport, neurological function, calcium binding/signaling, energy metabolism, cell growth/division, oxidative stress and ubiquitination. In total, our study indicated that RDX exposure affected molecular processes critical to CNS function in fathead minnow.

Behavioral and genomic impacts of a wastewater effluent on the fathead minnow

Rivers containing effluents from water treatment plants are complex soups of compounds, ranging from pharmaceuticals to natural hormones. Male fathead minnows (Pimephales promelas) were exposed for 3 weeks to effluent waters from the Metropolitan Wastewater Treatment Plant in St. Paul, MN. Fish were tested for their competitive nest holding behavior. Changes in vitellogenin were measured and these were correlated to changes in gene expression using a 22,000 gene microarray developed specifically for fathead minnows. Significant changes in gene expression were observed in both liver and testis, which correlate to phenotypic changes of vitellogenin induction and reduced competitive behavior. We also compared by real-time PCR the expression changes in key genes related to steroid biosynthesis and metabolism in fish exposed to the effluent as well as in fish exposed to a model estrogen and a model androgen. While the gene expression signature from effluent-exposed fish shared some elements with estrogen and androgen signatures, overall it was different, underscoring the complexity of compounds present in sewage and their different modes of action.

Reverse engineering adverse outcome pathways

The toxicological effects of many stressors are mediated through unknown, or incompletely characterized, mechanisms of action. The application of reverse engineering complex interaction networks from high dimensional omics data (gene, protein, metabolic, signaling) can be used to overcome these limitations. This approach was used to characterize adverse outcome pathways (AOPs) for chemicals that disrupt the hypothalamus-pituitary-gonadal endocrine axis in fathead minnows (FHM, Pimephales promelas). Gene expression changes in FHM ovaries in response to seven different chemicals, over different times, doses, and in vivo versus in vitro conditions, were captured in a large data set of 868 arrays. Potential AOPs of the antiandrogen flutamide were examined using two mutual information-based methods to infer gene regulatory networks and potential AOPs. Representative networks from these studies were used to predict network paths from stressor to adverse outcome as candidate AOPs. The relationship of individual chemicals to an adverse outcome can be determined by following perturbations through the network in response to chemical treatment, thus leading to the nodes associated with the adverse outcome. Identification of candidate pathways allows for formation of testable hypotheses about key biological processes, biomarkers, or alternative endpoints that can be used to monitor an AOP. Finally, the unique challenges facing the application of this approach in ecotoxicology were identified and a road map for the utilization of these tools presented.

Neuropeptides in the gonads: from evolution to pharmacology

Vertebrate gonads are the sites of synthesis and binding of many peptides that were initially classified as neuropeptides. These gonadal neuropeptide systems are neither well understood in isolation, nor in their interactions with other neuropeptide systems. Further, our knowledge of the control of these gonadal neuropeptides by peripheral hormones that bind to the gonads, and which themselves are under regulation by true neuropeptide systems from the hypothalamus, is relatively meager. This review discusses the existence of a variety of neuropeptides and their receptors which have been discovered in vertebrate gonads, and the possible way in which such systems could have evolved. We then focus on two key neuropeptides for regulation of the hypothalamo-pituitary-gonadal axis: gonadotropin-releasing hormone (GnRH) and gonadotropin-inhibitory hormone (GnIH). Comparative studies have provided us with a degree of understanding as to how a gonadal GnRH system might have evolved, and they have been responsible for the discovery of GnIH and its gonadal counterpart. We attempt to highlight what is known about these two key gonadal neuropeptides, how their actions differ from their hypothalamic counterparts, and how we might learn from comparative studies of them and other gonadal neuropeptides in terms of pharmacology, reproductive physiology and evolutionary biology.

In vivo assessment and potential diagnosis of xenobiotics that perturb the thyroid pathway: Proteomic analysis of Xenopus laevis brain tissue following exposure to model T4 inhibitors

As part of a multi-endpoint systems approach to develop comprehensive methods for assessing endocrine stressors in vertebrates, differential protein profiling was used to investigate expression patterns in the brain of the amphibian model (Xenopus laevis) following in vivo exposure to a suite of T4 synthesis inhibitors. We specifically address the application of Two Dimensional Polyacrylamide Gel Electrophoresis (2D PAGE), Isobaric Tags for Relative and Absolute Quantitation (iTRAQ) and LC-MS/MS to assess changes in relative protein expression levels. 2D PAGE and iTRAQ proved to be effective complementary techniques for distinguishing protein changes in the developing amphibian brain in response to T4 synthesis inhibition. This information served to evaluate the use of distinctive protein profiles as a potential mechanism to screen chemicals for endocrine activity in anurans. Regulatory pathways associated with proteins expressed as a result of chemical effect are reported. To our knowledge, this is also the first account of the anuran larvae brain proteome characterization using proteomic technologies. Correlation of protein changes to other cellular and organism-level responses will aid in the development of a more rapid and cost-effective, non-mammalian screening assay for thyroid axis-disrupting chemicals.

Site-specific impacts on gene expression and behavior in fathead minnows (Pimephales promelas) exposed in situ to streams adjacent to sewage treatment plants

BACKGROUND

Environmental monitoring for pharmaceuticals and endocrine disruptors in the aquatic environment traditionally employs a variety of methods including analytical chemistry, as well as a variety of histological and biochemical endpoints that correlate with the fish fitness. It is now clear that analytical chemistry alone is insufficient to identify aquatic environments that are compromised because these measurements do not identify the biologically available dose. The biological endpoints that are measured are important because they relate to known impairments; however, they are not specific to the contaminants and often focus on only a few known endpoints. These studies can be enhanced by looking more broadly at changes in gene expression, especially if the analysis focuses on biochemical pathways. The present study was designed to obtain additional information for well-characterized sites adjacent to sewage treatment plants in MN that are thought to be impacted by endocrine disruptors.

RESULTS

Here we examine five sites that have been previously characterized and examine changes in gene expression in fathead minnows (Pimephales promelas) that have been caged for 48 h in each of the aquatic environments. We find that the gene expression changes are characteristic and unique at each of the five sites. Also, fish exposed to two of the sites, 7 and 12, present a more aggressive behavior compared to control fish.

CONCLUSION

Our results show that a short-term exposure to sewage treatment plant effluents was able to induce a site-specific gene expression pattern in the fathead minnow gonad and liver. The short-term exposure was also enough to affect fish sexual behavior. Our results also show that microarray analysis can be very useful at determining potential exposure to chemicals, and could be used routinely as a tool for environmental monitoring.