Hepatic gene expression profiling using Genechips in zebrafish exposed to 17alpha-ethynylestradiol

A review of the influences of microplastics on toxicity and transgenerational effects of pharmaceutical and personal care products in aquatic environment

PPCPs (pharmaceutical and personal care products) and microplastics (MPs) are two types of emerging pollutants that are ubiquitous and widely concerned in the environment. Both of them can accumulate in fish or aquatic invertebrates and transfer to offspring, thereby producing toxic effects on both parents and offspring, in which the characteristics of MPs also enable them to adsorb PPCPs thus producing carrier effects. In this study, we have conducted a comprehensive review of MPs and PPCPs and found that MPs can act as a carrier of PPCPs to influence the bioaccumulation of PPCPs. MPs and PPCPs have toxicity and transgenerational effects on both fish and aquatic invertebrates in many aspects, and MPs can also affect the toxicity and transgenerational effects of PPCPs due to their carrier effects. This paper revealed that MPs may have an important impact on the bioavailability of PPCPs and the interaction between MPs and PPCPs is a hot topic in future research. This study also puts forward the shortcomings of the current research and related suggestions, and relevant research should be carried out as soon as possible to provide the basis for the prevention and treatment of fresh water.

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.

A Self-Organizing Map of the Fathead Minnow Liver Transcriptome to Identify Consistent Toxicogenomic Patterns across Chemical Fingerprints

Lack of consistent findings in different experimental settings remains a major challenge in toxicogenomics. The present study investigated whether consistency between findings of different microarray experiments can be improved when the analysis is based on a common reference frame ("toxicogenomic universe"), which can be generated using the machine learning algorithm of the self-organizing map (SOM). This algorithm arranges and clusters genes on a 2-dimensional grid according to their similarity in expression across all considered data. In the present study, 19 data sets, comprising of 54 different adult fathead minnow liver exposure experiments, were retrieved from Gene Expression Omnibus and used to train a SOM. The resulting toxicogenomic universe aggregates 58 872 probes to 2500 nodes and was used to project, visualize, and compare the fingerprints of these 54 different experiments. For example, we could identify a common pattern, with 14% of significantly regulated nodes in common, in the data sets of an interlaboratory study of ethinylestradiol exposures. Consistency could be improved compared with the 5% total overlap in regulated genes reported before. Furthermore, we could determine a specific and consistent estrogen-related pattern of differentially expressed nodes and clusters in the toxicogenomic universe by applying additional clustering steps and comparing all obtained fingerprints. Our study shows that the SOM-based approach is useful for generating comparable toxicogenomic fingerprints and improving consistency between results of different experiments. Environ Toxicol Chem 2020;39:526-537. © 2019 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals, Inc. on behalf of SETAC.

Experimental Approaches for Characterizing the Endocrine-Disrupting Effects of Environmental Chemicals in Fish

Increasing industrial and agricultural activities have led to a disturbing increase of pollutant discharges into the environment. Most of these pollutants can induce short-term, sustained or delayed impacts on developmental, physiological, and behavioral processes that are often regulated by the endocrine system in vertebrates, including fish, thus they are termed endocrine-disrupting chemicals (EDCs). Physiological impacts resulting from the exposure of these vertebrates to EDCs include abnormalities in growth and reproductive development, as many of the prevalent chemicals are capable of binding the receptors to sex steroid hormones. The approaches employed to investigate the action and impact of EDCs is largely dependent on the specific life history and habitat of each species, and the type of chemical that organisms are exposed to. Aquatic vertebrates, such as fish, are among the first organisms to be affected by waterborne EDCs, an attribute that has justified their wide-spread use as sentinel species. Many fish species are exposed to these chemicals in the wild, for either short or prolonged periods as larvae, adults, or both, thus, studies are typically designed to focus on either acute or chronic exposure at distinct developmental stages. The aim of this review is to provide an overview of the approaches and experimental methods commonly used to characterize the effects of some of the environmentally prevalent and emerging EDCs, including 17 α-ethinylestradiol, nonylphenol, BPA, phthalates, and arsenic; and the pervasive and potential carriers of EDCs, microplastics, on reproduction and growth. In vivo and in vitro studies are designed and employed to elucidate the direct effects of EDCs at the organismal and cellular levels, respectively. In silico approaches, on the other hand, comprise computational methods that have been more recently applied with the potential to replace extensive in vitro screening of EDCs. These approaches are discussed in light of model species, age and duration of EDC exposure.

In vivo assessment of gonad status, secondary sex characteristics and spawning in transparent Casper zebrafish

Important aspects of vertebrate reproduction, such as gametogenesis, involve changes in organs found deep internally and thus not easily studied directly in most living vertebrates due to obscuring pigment and overlying tissues. Transparent lines of zebrafish, especially the Casper double mutant, allow direct observation and analysis of reproductive events in the gonads in vivo. The natural production of fertilized eggs in zebrafish is a complex process involving oogenesis, spermatogenesis, mating behavior, endocrine and neurological processes with inputs from the environment including light, temperature and nutrition. While these factors play important roles, the hypothalamic-pituitary-gonadal axis (HPGA) is central in the regulation of embryo output. Endocrine disrupting compounds (EDCs) include a variety of pollutants often present in the environment. EDCs may have direct effects on the HPGA or indirect effects through toxic action on supporting organs such as the liver or kidney. Estrogenic compounds such as diethylstilbestrol (DES) have been reported to affect reproduction in a variety of species including man. In this study, the effects of DES on reproduction were determined in a novel way by using transparent Casper zebrafish that allow direct visualization of gonad status over time. Changes in gonad status with DES treatment were correlated with effects on secondary sex characteristics (i.e., genital vent size and breeding tubercles) spawning and embryo production. The results suggest that the Casper zebrafish is a useful model for studying dynamics of reproductive events in vertebrate gonads in vivo and for determining effects of EDCs on zebrafish reproduction.

Estrogen levels in surface sediments from a multi-impacted Brazilian estuarine system

Estrogen levels were assessed in surface sediments from one of the most industrialized and urbanized estuarine systems in Latin America (SSES, Santos and São Vicente estuarine system). Estriol (E3) presented quantifiable levels in all sampled sites, ranging from 20.9 ng g^-1 to 694.2 ng g^-1. 17β-estradiol (E2) and 17α-ethinylestradiol (EE2) were also detected in almost all sampled sites. The highest concentration of E2 was 23.9 ng g^-1, whereas high levels of EE2 86.3 ng g^-1. The occurrence of estrogens in SSES was diffuse and partially related to a domestic sewage outfall. Estrogens were also found in areas with substantial contribution of sanitary effluents from domiciles not covered by sanitation services. Our results reinforce that studies on environmental contamination by estrogens should not be spatially limited to the vicinities of point sources. These results contribute to raise awareness on the need of a formal approach to assess ecological risks of estrogens in the SSES.

Transcriptomic analysis of short-term 17α-ethynylestradiol exposure in two Californian sentinel fish species sardine (Sardinops sagax) and mackerel (Scomber japonicus)

Endocrine disrupting chemicals (EDCs) are substances which disrupt normal functioning of the endocrine system by interfering with hormone regulated physiological pathways. Aquatic environments provide the ultimate reservoir for many EDCs as they enter rivers and the ocean via effluent discharges and accumulate in sediments. One EDC widely dispersed in municipal wastewater effluent discharges is 17α-ethynylestradiol (EE2), which is one of the most widely prescribed medicines. EE2 is a bio-active estrogen employed in the majority of oral contraceptive pill formulations. As evidence of the health risks posed by EDCs mount, there is an urgent need to improve diagnostic tools for monitoring the effects of pollutants. As the cost of high throughput sequencing (HTS) diminishes, transcriptional profiling of an organism in response to EDC perturbation presents a cost-effective way of screening a wide range of endocrine responses. Coastal pelagic filter feeding fish species analyzed using HTS provide an excellent tool for EDC risk assessment in the marine environment. Unfortunately, there are limited genome sequence data and annotation for many of these species including Pacific sardine (Sardinops sagax) and chub mackerel (Scomber japonicus), which limits the utility of molecular tools such as HTS to interrogate the effects of endocrine disruption. In this study, we carried out RNA sequencing (RNAseq) of liver RNA harvested from wild sardine and mackerel exposed for 5 h under laboratory conditions to a concentration of 12.5 pM EE2 in the tank water. We developed an analytical framework for transcriptomic analyses of species with limited genomic information. EE2 exposure altered expression patterns of key genes involved in important metabolic and physiological processes. The systems approach presented here provides a powerful tool for obtaining a comprehensive picture of endocrine disruption in aquatic organisms.

Early-life exposure to the endocrine disruptor 17-α-ethinylestradiol induces delayed effects in adult brain, liver and ovotestis proteomes of a self-fertilizing fish

Early-life represents a critically sensitive window to endocrine disrupting chemicals, potentially leading to long-term repercussions on the phenotype later in life. The mechanisms underlying this phenomenon, referred to as the Developmental Origins of Health and Disease (DOHaD), are still poorly understood. To gain molecular understanding of these effects, we exposed mangrove rivulus (Kryptolebias marmoratus) for 28 days post hatching (dph) to 4 and 120 ng/L 17-α-ethinylestradiol, a model xenoestrogen. After 28 days, fish were raised for 140 days in clean water and we performed quantitative label-free proteomics on brain, liver and ovotestis of 168 dph adults. A total of 820, 888 and 420 proteins were robustly identified in the brain, liver and ovotestis, respectively. Effects of 17-α-ethinylestradiol were tissue and dose-dependent: a total of 31, 51 and 18 proteins were differentially abundant at 4 ng/L in the brain, liver and ovotestis, respectively, compared to 20, 25 and 39 proteins at 120 ng/L. Our results suggest that estrogen-responsive pathways, such as lipid metabolism, inflammation, and the innate immune system were affected months after the exposure. In addition, the potential perturbation of S-adenosylmethionine metabolism encourages future studies to investigate the role of DNA methylation in mediating the long-term effects of early-life exposures. SIGNIFICANCE: The Developmental Origins of Health and Disease (DOHaD) states that early life stages of humans and animals are sensitive to environmental stressors and can develop health issues later in life, even if the stress has ceased. Molecular mechanisms supporting DOHaD are still unclear. The mangrove rivulus is a new fish model species naturally reproducing by self-fertilization, making it possible to use isogenic lineages in which all individuals are highly homozygous. This species therefore permits to strongly reduce the confounding factor of genetic variability in order to investigate the effects of environmental stress on the phenotype. After characterizing the molecular phenotype of brain, liver and ovotestis, we obtained true proteomic reaction norms of these three organs in adults after early life stages have been exposed to the common endocrine disruptor 17-α-ethinylestradiol (EE2). Our study demonstrates long-term effects of early-life endocrine disruption at the proteomic level in diverse estrogen-responsive pathways 5 months after the exposure. The lowest tested and environmentally relevant concentration of 4 ng/L had the highest impact on the proteome in brain and liver, highlighting the potency of endocrine disruptors at low concentrations.

Estrogen-responsive gene networks in the teleost liver: What are the key molecular indicators?

An overarching goal of environmental genomics is to leverage sensitive suites of markers that are robust and reliable to assess biological responses in a range of species inhabiting variable environments. The objective of this study was to identify core groups of transcripts and molecular signaling pathways that respond to 17alpha-ethylinestadiol (EE2), a ubiquitous estrogenic contaminant, using transcriptome datasets generated from six independent laboratories. We sought to determine which biomarkers and gene networks were those most robust and reliably detected in multiple laboratories. Six laboratories conducted microarray analysis in pieces of the same liver from male fathead minnows exposed to ∼15ng/L EE2 for 96h. There were common transcriptional networks identified in every dataset. These included down-regulation of gene networks associated with blood clotting, complement activation, triglyceride storage, and xenobiotic metabolism. Noteworthy was that more than ∼85% of the gene networks were suppressed by EE2. Leveraging both these data and those mined from the Comparative Toxicogenomics Database (CTD), we narrowed in on an EE2-responsive transcriptional network. All transcripts in this network responded ∼±5-fold or more to EE2, increasing reliability of detection. This network included estrogen receptor alpha, transferrin, myeloid cell leukemia 1, insulin like growth factor 1, insulin like growth factor binding protein 2, and methionine adenosyltransferase 2A. This estrogen-responsive interactome has the advantage over single markers (e.g. vitellogenin) in that these entities are directly connected to each other based upon evidence of expression regulation and protein binding. Thus, it represents an interacting functional suite of estrogenic markers. Vitellogenin, the gold standard for estrogenic exposures, can show high individual variability in its response to estrogens, and the use of a multi-gene approach for estrogenic chemicals is expected to improve sensitivity. In our case, the coefficient of variation was significantly lowered by the gene network (∼67%) compared to Vtg alone, supporting the use of this transcriptional network as a sensitive alternative for detecting estrogenic effluents and chemicals. We propose that screening chemicals for estrogenicity using interacting genes within a defined expression network will improve sensitivity, accuracy, and reduce the number of animals required for endocrine disruption assessments.

Global transcriptomic analysis of zebrafish in response to embryonic exposure to three antidepressants, amitriptyline, fluoxetine and mianserin

Antidepressants are among the most commonly detected pharmaceuticals in aqueous systems, and, as emerging organic pollutants, may exert negative effects on non-target aquatic organisms. Previously, it has been revealed that antidepressant exposure significantly inhibits the growth and development of fish during their early developmental stages. Thus, in the present study, we aimed to identify and compare the underlying mechanisms of action of different antidepressants at the transcriptional level using zebrafish (Danio rerio) embryos. Through high-throughput RNA sequencing (RNA-Seq) data analysis, 32, 34, and 130 differentially expressed genes (DEGs) were obtained from zebrafish larvae after 120h of embryonic exposure to sublethal concentrations of amitriptyline, fluoxetine, and mianserin, respectively. The expression profiles of the identified DEGs showed similar trends in response to the three antidepressant treatments, suggesting consistent toxic effects of low concentrations of these three drugs on the regulation of gene expression in fish. Several metabolic and signaling pathways, including glycolysis/gluconeogenesis and the insulin pathway, were affected in the exposed fish larvae. The expression profiles of selected DEGs were then verified by the qRT-PCR method, which indicated significant positive correlations with the RNA-Seq results. Next, we determined the concentration-dependent expression patterns of 6 selected DEGs in fish larvae exposed to three antidepressants at a series of environmentally relevant concentrations. The results revealed a significant concentration-dependent reduction in the levels of dual-specificity phosphatase 5 (dusp5) mRNA, as well as a non-concentration-dependent gene expression inhibition of prostaglandin D2 synthase b (ptgdsb); the circadian rhythm-related genes, i.e. those encoding nuclear receptor subfamily 1, group D, member 1 (nr1d1) and period 2 (per2); and genes encoding early growth response factors (egr1 and egr4), in the antidepressant-treated fish larvae. In summary, to our knowledge, our findings demonstrate, for the first time, that the three different categories of antidepressants have common effects on the gene expression involved in multiple biological processes and signaling pathways during the early development of fish and thus provide information for characterizing the adverse outcome pathways and on the ecological risk assessment of these pharmaceutical pollutants in the aquatic environment.

How consistent are we? Interlaboratory comparison study in fathead minnows using the model estrogen 17α-ethinylestradiol to develop recommendations for environmental transcriptomics

Fundamental questions remain about the application of omics in environmental risk assessments, such as the consistency of data across laboratories. The objective of the present study was to determine the congruence of transcript data across 6 independent laboratories. Male fathead minnows were exposed to a measured concentration of 15.8 ng/L 17α-ethinylestradiol (EE2) for 96 h. Livers were divided equally and sent to the participating laboratories for transcriptomic analysis using the same fathead minnow microarray. Each laboratory was free to apply bioinformatics pipelines of its choice. There were 12 491 transcripts that were identified by one or more of the laboratories as responsive to EE2. Of these, 587 transcripts (4.7%) were detected by all laboratories. Mean overlap for differentially expressed genes among laboratories was approximately 50%, which improved to approximately 59.0% using a standardized analysis pipeline. The dynamic range of fold change estimates was variable between laboratories, but ranking transcripts by their relative fold difference resulted in a positive relationship for comparisons between any 2 laboratories (mean R²  > 0.9, p < 0.001). Ten estrogen-responsive genes encompassing a fold change range from dramatic (>20-fold; e.g., vitellogenin) to subtle (∼2-fold; i.e., block of proliferation 1) were identified as differentially expressed, suggesting that laboratories can consistently identify transcripts that are known a priori to be perturbed by a chemical stressor. Thus, attention should turn toward identifying core transcriptional networks using focused arrays for specific chemicals. In addition, agreed-on bioinformatics pipelines and the ranking of genes based on fold change (as opposed to p value) should be considered in environmental risk assessment. These recommendations are expected to improve comparisons across laboratories and advance the use of omics in regulations. Environ Toxicol Chem 2017;36:2593-2601. © 2017 SETAC.

Characterizing combined effects of antiestrogenic chemicals on vitellogenin production in rainbow trout (Oncorhynchus mykiss) hepatocytes

Fish are exposed to a complex mixture of endocrine disrupting compounds (EDC), some of which display antiestrogenic activity leading to suppression of estrogen receptor (ER)- mediated reproductive processes. Although the main mode of action (MoA) of these antiestrogens is to directly interfere with natural ligand binding of the ER, several other MoA have been proposed. The aim of the present study was to characterize single and combined antiestrogenic effects of the aryl hydrocarbon receptor (AhR)-agonist β-naphthoflavone (BNF) and ER-antagonist 4-hydroxytamoxifen (OHT) on vitellogenin (Vtg) protein using primary rainbow trout (Oncorhynchus mykiss) hepatocytes. Supporting transcriptional analysis of ER-responsive genes (estrogen receptor-α (er-α), vitellogenin-1 (vtg-1), eggshell zona radiata protein (zrp)) and AhR-mediated genes (aryl hydrocarbon receptor-2β, cytochrome p450-1a (cyp1a)) was performed by qPCR to characterize the antiestrogenic influence on ER- and AhR-mediated responses. Data demonstrated that both BNF and OHT significantly reduced 17β-estradiol (E2)-induced Vtg protein expression in a concentration responsive manner, whereas exposure to a mixture of these produced an additive antiestrogenic effect. The results observed at the protein level were further supported by transcriptional analysis of ER-responsive genes (er-α, vtg-1, zrp), where only E2-induced vtg-1 gene expression was significantly decreased by OHT and the mixture of OHT and BNF. E2-induced er-α and zrp gene expression was not markedly altered. The significant reduction of E2-induced vtg-1 gene expression by OHT suggested that the antiestrogenic effect of this compound may be associated with ER signaling pathway. Specific genes involved in putative AhR-ER cross-talk were also investigated, however none were directly associated with the compound anti-estrogenic MoA. Although the MoA of the single compounds and mixture were not completely characterized, the present study enhanced our knowledge of the combined toxicity mediated by antiestrogens acting through different MoA.

Evaluation of ecotoxicological effects of benzophenone UV filters: Luminescent bacteria toxicity, genotoxicity and hormonal activity

The widespread use of organic ultraviolet (UV) filters in personal care products raises concerns about their potentially hazardous effects on human and ecosystem health. In this study, the toxicities of four commonly used benzophenones (BPs) UV filters including benzophenone (BP), 2-Hydroxybenzophenone (2HB), 2-Hydroxy-4-methoxybenzophenone (BP3), and 2-Hydroxy-4-methoxybenzophenone-5-sulfonicacid (BP4) in water were assayed in vitro using Vibrio fischeri, SOS/umu assay, and yeast estrogen screen (YES) assay, as well as in vivo using zebrafish larvae. The results showed that the luminescent bacteria toxicity, expressed as logEC₅₀, increased with the lipophilicity (logKow) of BPs UV filters. Especially, since 2HB, BP3 and BP4 had different substituent groups, namely -OH, -OCH₃ and -SO₃H, respectively, these substituent functional groups had a major contribution to the lipophilicity and acute toxicity of these BPs. Similar tendency was observed for the genotoxicity, expressed as the value of induction ratio=1.5. Moreover, all the target BPs UV filters showed estrogenic activity, but no significant influences of lipophilicity on the estrogenicity were observed, with BP3 having the weakest estrogenic efficiency in vitro. Although BP3 displayed no noticeable adverse effects in any in vitro assays, multiple hormonal activities were observed in zebrafish larvae including estrogenicity, anti-estrogenicity and anti-androgenicity by regulating the expression of target genes. The results indicated potential hazardous effects of BPs UV filters and the importance of the combination of toxicological evaluation methods including in vitro and in vivo assays.

Microcystin-LR retards gonadal maturation through disrupting the growth hormone/insulin-like growth factors system in zebrafish

Recent studies have documented that microcystins (MCs) have potential toxic effects on growth and reproduction in fish. However, no systematic data exist on whether MCs cause gonadal development retardation through disrupting the growth hormone/insulin-like growth factors (GH/IGFs) system. To this end, zebrafish hatchlings (5 d post-fertilization) were exposed to 0, 0.3, 3 and 30µg/L microcystin-LR (MC-LR) for 90 d until they reached sexual maturity. Life-cycle exposure to MC-LR caused delayed ovarian maturation and sperm development along with ultrapathological lesions in the brain and liver. Moreover, the retarded gonadal development was accompanied by an inhibition of the GH/IGFs system, which was characterized by significant decreases in the transcriptional levels of brain gh (males only), hepatic igf2a and igf2b as well as gonadal igf1 (males only), igf3 and igf2r. These findings for the first time point to the influence of MC-LR on fish gonadal development via the GH/IGFs system. Also, sex-differential impairments suggested that gonadal development of males is more vulnerable than that of female to MC-LR. Our results provide evidence that MC-LR at environmentally relevant concentrations is able to induce impairments on fish gonadal development.

Characterisation of the transcriptome of male and female wild-type guppy brains with RNA-Seq and consequences of exposure to the pharmaceutical pollutant, 17α-ethinyl estradiol

Waterways are increasingly being contaminated by chemical compounds that can disrupt the endocrinology of organisms. One such compound is 17α-ethinyl estradiol (EE2), a synthetic estrogen used in the contraceptive pill. Despite considerable research interest in the effects of EE2 on reproduction and gene expression, surprisingly, only a few studies have capitalised on technologies, such as next-generation sequencing (NGS), to uncover the molecular pathways related to EE2 exposure. Accordingly, using high-throughput sequencing technologies, the aim of our study was to explore the effects of EE2 on brain transcriptome in wild-type male and female guppy (Poecilia reticulata). We conducted two sets of experiments, where fish were exposed to EE2 (measured concentrations: 8ng/L and 38ng/L) in a flow-through system for 21days. The effects on the brain transcriptome on both males and females were assessed using Illumina sequencing (MiSeq and HiSeq) platform followed by bioinformatics analysis (edgeR, DESeq2). Here, we report that exposure to EE2 caused both up- and downregulation of specific transcript abundances, and affected transcript abundance in a sex-specific manner. Specifically, we found 773 transcripts, of which 60 were male-specific, 61 female-specific and 285 treatment-specific. EE2 affected expression of 165 transcripts in males, with 88 downregulated and 77 upregulated, while in females, 120 transcripts were affected with 62 downregulated and 58 upregulated. Finally, RT-qPCR validation demonstrated that expression of transcripts related to transposable elements, neuroserpin and heat shock protein were significantly affected by EE2-exposure. Our study is the first to report brain transcriptome libraries for guppies exposed to EE2. Not only does our study provide a valuable resource, it offers insights into the mechanisms underlying the feminizing effects on the brains of organisms exposed to environmentally realistic concentrations of EE2.

Persistent Effects of Developmental Exposure to 17α-Ethinylestradiol on the Zebrafish (Danio rerio) Brain Transcriptome and Behavior

The synthetic estrogen 17α-ethinylestradiol (EE₂) is an endocrine disrupting compound of concern due to its persistence and widespread presence in the aquatic environment. Effects of developmental exposure to low concentrations of EE₂ in fish on reproduction and behavior not only persisted to adulthood, but have also been observed to be transmitted to several generations of unexposed progeny. To investigate the possible biological mechanisms of the persistent anxiogenic phenotype, we exposed zebrafish embryos for 80 days post fertilization to 0, 3, and 10 ng/L EE₂ (measured concentrations 2.14 and 7.34 ng/L). After discontinued exposure, the animals were allowed to recover for 120 days in clean water. Adult males and females were later tested for changes in stress response and shoal cohesion, and whole-brain gene expression was analyzed with RNA sequencing. The results show increased anxiety in the novel tank and scototaxis tests, and increased shoal cohesion in fish exposed during development to EE₂. RNA sequencing revealed 34 coding genes differentially expressed in male brains and 62 in female brains as a result of EE₂ exposure. Several differences were observed between males and females in differential gene expression, with only one gene, sv2b, coding for a synaptic vesicle protein, that was affected by EE₂ in both sexes. Functional analyses showed that in female brains, EE₂ had significant effects on pathways connected to the circadian rhythm, cytoskeleton and motor proteins and synaptic proteins. A large number of non-coding sequences including 19 novel miRNAs were also differentially expressed in the female brain. The largest treatment effect in male brains was observed in pathways related to cholesterol biosynthesis and synaptic proteins. Circadian rhythm and cholesterol biosynthesis, previously implicated in anxiety behavior, might represent possible candidate pathways connecting the transcriptome changes to the alterations to behavior. Further the observed alteration in expression of genes involved in synaptogenesis and synaptic function may be important for the developmental modulations resulting in an anxiety phenotype. This study represents an initial survey of the fish brain transcriptome by RNA sequencing after long-term recovery from developmental exposure to an estrogenic compound.

Reproductive and transcriptional effects of the antiandrogenic progestin chlormadinone acetate in zebrafish (Danio rerio)

Chlormadinone acetate (CMA) is a frequently used progestin with antiandrogenic activity in humans. Residues may enter the aquatic environment but potential adverse effects in fish are unknown. While our previous work focused on effects of CMA in vitro and in zebrafish eleuthero-embryos, the present study reports on reproductive and transcriptional effects in adult female and male zebrafish (Danio rerio). We performed a reproductive study using breeding groups of zebrafish. After 15 days of pre-exposure, we exposed zebrafish to different measured concentrations between 6.4 and 53,745 ng/L CMA for 21 days and counted produced eggs daily to determine fecundity. Additionally, transcriptional effects of CMA in brains, livers, and gonads were analyzed. CMA induced a slight but statistically significant reduction in fecundity at 65 ng/L and 53,745 ng/L compared to pre-exposure. Furthermore, we observed differential expression for gene transcripts of steroid hormone receptors, genes related to the hypothalamic-pituitary-gonadal axis, and steroidogenesis. In particular, we found a significant decrease of transcript levels of vitellogenin (vtg1) in ovaries and liver, and of cyp2k7 in the liver of males, as well as a significant increase of transcripts of the progesterone receptor (pgr) in testes, and cyp2k1 in the liver of females. The observed effects were weaker than those of other very potent progestins, which is probably related to the lack of interaction of CMA with the zebrafish progesterone receptor.

Molecular initiating events of the intersex phenotype: Low-dose exposure to 17α-ethinylestradiol rapidly regulates molecular networks associated with gonad differentiation in the adult fathead minnow testis

Intersex, or the presence of oocytes in the testes, has been documented in fish following exposure to wastewater effluent and estrogenic compounds. However, the molecular networks underlying the intersex condition are not completely known. To address this, we exposed male fathead minnows to a low, environmentally-relevant concentration of 17alpha-ethinylestradiol (EE2) (15ng/L) and measured the transcriptome response in the testis after 96h to identify early molecular initiating events that may proceed the intersex condition. The short-term exposure to EE2 did not affect gonadosomatic index and proportion of gametes within the testes. However, the production of 11-ketotestosterone and testosterone from the testis in vitro was decreased relative to controls. Expression profiling using a 8×60K fathead minnow microarray identified 10 transcripts that were differentially expressed in the testes, the most dramatic change being that of coagulation factor XIII A chain (20-fold increase). Transcripts that included guanine nucleotide binding protein (Beta Polypeptide 2), peroxisome proliferator-activated receptor delta, and WNK lysine deficient protein kinase 1a, were down-regulated by EE2. Subnetwork enrichment analysis revealed that EE2 suppressed transcriptional networks associated with steroid metabolism, hormone biosynthesis, and sperm mobility. Most interesting was that gene networks associated with doublesex and mab-3 related transcription factor 1 (dmrt1) were suppressed in the adult testis, despite the fact that dmrt1 itself was not different in expression from control males. Transcriptional networks involving forkhead box L2 (foxl2) (transcript involved in ovarian follicle development) were increased in expression in the testis. Noteworthy was that a gene network associated to granulosa cell development was increased over 100%, suggesting that this transcriptome network may be important for monitoring estrogenic exposures. Other cell processes rapidly downregulated by EE2 at the transcript level included glucose homeostasis, response to heavy metal, amino acid catabolism, and the cyclooxygenase pathway. Conversely, lymphocyte chemotaxis, intermediate filament polymerization, glucocorticoid metabolism, carbohydrate utilization, and anterior/posterior axis specification were increased. These data provide new insight into the transcriptional responses that are perturbed prior to gonadal remodeling and intersex following exposure to estrogens. These data demonstrate that low concentrations of EE2 (1) rapidly suppresses male hormone production, (2) down-regulate molecular networks related to male sex differentiation, and (3) induce transcriptional networks related to granulosa cell development in the adult testis. These responses are hypothesized to be key molecular initiating events that occur prior to the development of the intersex phenotype following estrogenic exposures.

Assessment of multiple hormone activities of a UV-filter (octocrylene) in zebrafish (Danio rerio)

In this study, zebrafish (Danio rerio) were exposed to a UV-filter-octocrylene (OCT) with elevated concentrations for 28 d. The total body accumulation of OCT in zebrafish was found to reach 2321.01 ("L" level), 31,234.80 ("M" level), and 70,593.38 ng g(-1) ("H" level) when the average OCT exposure concentration was controlled at 28.61, 505.62, and 1248.70 μg L(-1), respectively. Gross and histological observations as well as RT-qPCR analysis were conducted to determine the effects of OCT accumulation on zebrafish. After exposure, the gonad-somatic index and percentage of vitellogenic oocytes were found to increase significantly in the ovaries of female zebrafish at the H accumulation level. Significant up-regulation of esr1 and cyp19b were observed in the gonads, as well as vtg1 in the livers for both female and male zebrafish. At M and H accumulation levels, apparent down-regulation of ar was observed in the ovaries and testis of the female and male zebrafish, respectively. Although the extent of the effects on zebrafish differed at different accumulation levels, the induction of vtg1 and histological changes in the ovaries are indications of estrogenic activity and the inhibition of esr1 and ar showed antiestrogenic and antiandrogenic activity, respectively. Thus, as OCT could easily accumulate in aquatic life such as zebrafish, one of its most of concern hazards would be the disturbance of the histological development and its multiple hormonal activities.

Activity of binary mixtures of drospirenone with progesterone and 17α-ethinylestradiol in vitro and in vivo

Despite potential exposure of aquatic organisms to mixtures of steroid hormones, very little is known on their joint activity in fish. Drospirenone (DRS) is a new synthetic progestin used in contraceptive pills in combination with 17α-ethinylestradiol (EE2). Here we systematically analyzed effects of DRS in binary mixtures with progesterone (P4) and EE2. First, we determined the in vitro activity of single compounds in recombinant yeast assays that express the human progesterone, androgen, or estrogen receptor, followed by determination of mixture activities of DRS and P4, DRS and EE2, as well as medroxyprogesterone acetate (MPA) and dydrogesterone (DDG). Mixtures of DRS and P4, as well as of DRS and EE2 showed additive progestogenic and androgenic activities. However, DDG and MPA showed non-additive progestogenic and androgenic activities. We then analyzed the in vivo activity of single compounds and mixtures of DRS and P4, as well as DRS and EE2, by assessing transcriptional changes of up to 14 selected target genes in zebrafish embryos at 48h post fertilization (hpf), and in eleuthero-embryos at 96hpf and 144hpf. DRS, P4, and EE2 led to significant transcriptional alteration of genes, including those encoding hormone receptors (pgr, esr1), a steroidogenic enzyme (hsd17b3), and estrogenic markers (vtg1, cyp19b), in particular at 144 hpf. In general, DRS showed stronger transcriptional changes than P4. In mixtures of DRS and P4, they were mainly non-additive (antagonistic interaction). In mixtures of DRS and EE2, transcriptional responses of esr1, vtg1 and cyp19b were dominated by EE2, suggesting an antagonistic interaction or independent action. Equi-effective mixtures of DRS and EE2, based on progesterone receptor transcripts, showed antagonistic interactions. Our data suggest that interactions in mixtures assessed in vitro in recombinant yeast cannot be translated to the in vivo situation. The receptor-based responses did not correspond well to the transcriptional responses in embryos which are much more complex due to the interplay between hormonal pathways, receptor crosstalk, and hormonal feedback loops.

17α-Ethinylestradiol (EE2) effect on global gene expression in primary rainbow trout (Oncorhynchus mykiss) hepatocytes

The potential impact of endocrine disrupting chemicals (EDCs) in the aquatic environment has driven the development of screening assays to evaluate the estrogenic properties of chemicals and their effects on aquatic organisms such as fish. However, obtaining full concentration-response relationships in animal (in vivo) exposure studies are laborious, costly and unethical, hence a need for developing feasible alternative (non-animal) methods. Use of in vitro bioassays such as primary fish hepatocytes, which retain many of the native properties of the liver, has been proposed for in vitro screening of estrogen receptor (ER) agonists and antagonists. The aim of present study was to characterize the molecular mode of action (MoA) of the ER agonist 17α-ethinylestradiol (EE2) in primary rainbow trout (Oncorhynchus mykiss) hepatocytes. A custom designed salmonid 60,000-feature (60k) oligonucleotide microarray was used to characterize the potential MoAs after 48h exposure to EE2. The microarray analysis revealed several concentration-dependent gene expression alterations including classical estrogen sensitive biomarker gene expression (e.g. estrogen receptor α, vitellogenin, zona radiata). Gene Ontology (GO) analysis displayed transcriptional changes suggesting interference of cellular growth, fatty acid and lipid metabolism potentially mediated through the estrogen receptor (ER), which were proposed to be associated with modulation of genes involved in endocrine function and reproduction. Pathway analysis supported the identified GOs and revealed modulation of additional genes associated with apoptosis and cholesterol biosynthesis. Differentially expressed genes (DEGs) related to impaired lipid metabolism (e.g. peroxisome proliferator-activated receptor α and γ), growth (e.g. insulin growth factor protein 1), phase I and II biotransformation (e.g. cytochrome P450 1A, sulfotransferase, UDP-glucuronosyltransferase and glutathione S-transferase) provided additional insight into the MoA of EE2 in primary fish hepatocytes. Results from the present study suggest that biotransformation, estrogen receptor-mediated responses, lipid homeostasis, growth and cancer/apoptosis in primary fish hepatocytes may be altered after short-term exposure to ER-agonists such as EE2. In many cases the observed changes were similar to those reported for estrogen-exposed fish in vivo. In conclusion, global transcriptional analysis demonstrated that EE2 affected a number of toxicologically relevant pathways associated with an estrogenic MoA in the rainbow trout hepatocytes.

The Effects of Synthetic Estrogen Exposure on the Sexually Dimorphic Liver Transcriptome of the Sex-Role-Reversed Gulf Pipefish

Species exhibiting sex-role reversal provide an unusual perspective on the evolution of sex roles and sex differences. However, the proximate effects of sex-role reversal are largely unknown. Endocrine disruptors provide an experimental mechanism to address hormonal regulation of sexually dimorphic gene expression in sex-role-reversed taxa. Here, we investigate gene expression patterns in the liver of the sex-role-reversed Gulf pipefish, because the liver is known to be sexually dimorphic and estrogen-regulated in species with conventional sex roles. Using next-generation RNA-sequencing technology (RNA-seq), we detected sexually dimorphic hepatic gene expression patterns, with a total of 482 differentially expressed genes between the sexes in Gulf pipefish. Two-thirds of these genes were over-expressed in females, and the sex-specific transcriptomes of this sex-role-reversed pipefish’s liver were superficially similar to those of fishes with conventional sex-roles. We exposed females, pregnant males, and non-pregnant males to 17α-ethinylestradiol (EE2) at ecologically relevant concentrations of 5ng/L and compared gene expression patterns in the livers of exposed fish to control fish. Several genes that were up-regulated in EE2-exposed males relative to control males were also found to be female-biased in control animals. These genes included several of the classic estrogen biomarkers, such as vitellogenin, choriogenin, and zona pellucida. Thus, estrogen exposure induced feminization of the male liver transcriptome in a sex-role-reversed pipefish. These results suggest that the ancestral state of estrogen-regulated female reproductive physiology has been retained in all sex-role-reversed vertebrates thus far studied, despite substantial evolution of the hormonal regulation of ornamentation and mating behavior in these interesting taxa.

Adaptation of a Bioinformatics Microarray Analysis Workflow for a Toxicogenomic Study in Rainbow Trout

Sex steroids play a key role in triggering sex differentiation in fish, the use of exogenous hormone treatment leading to partial or complete sex reversal. This phenomenon has attracted attention since the discovery that even low environmental doses of exogenous steroids can adversely affect gonad morphology (ovotestis development) and induce reproductive failure. Modern genomic-based technologies have enhanced opportunities to find out mechanisms of actions (MOA) and identify biomarkers related to the toxic action of a compound. However, high throughput data interpretation relies on statistical analysis, species genomic resources, and bioinformatics tools. The goals of this study are to improve the knowledge of feminisation in fish, by the analysis of molecular responses in the gonads of rainbow trout fry after chronic exposure to several doses (0.01, 0.1, 1 and 10 μg/L) of ethynylestradiol (EE2) and to offer target genes as potential biomarkers of ovotestis development. We successfully adapted a bioinformatics microarray analysis workflow elaborated on human data to a toxicogenomic study using rainbow trout, a fish species lacking accurate functional annotation and genomic resources. The workflow allowed to obtain lists of genes supposed to be enriched in true positive differentially expressed genes (DEGs), which were subjected to over-representation analysis methods (ORA). Several pathways and ontologies, mostly related to cell division and metabolism, sexual reproduction and steroid production, were found significantly enriched in our analyses. Moreover, two sets of potential ovotestis biomarkers were selected using several criteria. The first group displayed specific potential biomarkers belonging to pathways/ontologies highlighted in the experiment. Among them, the early ovarian differentiation gene foxl2a was overexpressed. The second group, which was highly sensitive but not specific, included the DEGs presenting the highest fold change and lowest p-value of the statistical workflow output. The methodology can be generalized to other (non-model) species and various types of microarray platforms.

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

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

Influence of short-term exposure to low levels of 17α-ethynylestradiol on expression of genes involved in immunity and on immune parameters in rainbow trout, Oncorhynchus mykiss

Fish are exposed to endocrine-disrupting chemicals (EDC), which are well known to disturb not only the reproductive system but also the immune system in vertebrates. However, the mechanisms by which these compounds are able to modify fish immunity are not well understood. In order to test the EE2 effects on immunity in selected organs, we exposed rainbow trout male juveniles for 3 weeks to EE2 concentrations ranging from 0.01 to 1 μg/L. The results of this study suggest that EE2 affects the immunity of rainbow trout in a tissue dependent manner. This molecule affects both cellular and humoral immune systems. Indeed, blood leukocyte populations, as well as hepatic and plasma lysozyme, plasma MPO and renal complement activities, are modulated by EE2. Moreover, EE2 alters the gene expression of some mucus compounds, hepatic expression of complement sub-unit and lysozyme, or genes involved in the hepatic phagocytosis and transport of immunoglobulin across the liver.

Dietary exposure of 17-alpha ethinylestradiol modulates physiological endpoints and gene signaling pathways in female largemouth bass (Micropterus salmoides)

17Alpha-ethinylestradiol (EE2), used for birth control in humans, is a potent estrogen that is found in wastewater at low concentrations (ng/l). EE2 has the ability to interfere with the endocrine system of fish, affecting reproduction which can result in population level effects. The objective of this study was to determine if dietary exposure to EE2 would alter gene expression patterns and key pathways in the liver and ovary and whether these could be associated with reproductive endpoints in female largemouth bass during egg development. Female LMB received 70ng EE2/g feed (administered at 1% of body weight) for 60 days. EE2 dietary exposure significantly reduced plasma vitellogenin concentrations by 70%. Hepatosomatic and gonadosomatic indices were also decreased with EE2 feeding by 38.5% and 40%, respectively. Transcriptomic profiling revealed that there were more changes in steady state mRNA levels in the liver compared to the ovary. Genes associated with reproduction were differentially expressed, such as vitellogenin in the liver and aromatase in the gonad. In addition, a set of genes related with oxidative stress (e.g. glutathione reductase and glutathione peroxidase) were identified as altered in the liver and genes associated with the immune system (e.g. complement component 1, and macrophage-inducible C-type lectin) were altered in the gonad. In a follow-up study with 0.2ng EE2/g feed for 60 days, similar phenotypic and gene expression changes were observed that support these findings with the higher concentrations. This study provides new insights into how dietary exposure to EE2 interferes with endocrine signaling pathways in female LMB during a critical period of reproductive oogenesis.

Gonadal development and transcript profiling of steroidogenic enzymes in response to 17α-methyltestosterone in the rare minnow Gobiocypris rarus

It is well known that natural and anthropogenic chemicals interfere with the hormonal system of vertebrate and invertebrate organisms. How these chemicals regulate gonadal steroidogenesis remains to be determined. The main objective of this study was to evaluate the effects of 17α-methyltestosterone (MT), a synthetic model androgen, on gene expression profiles of six key steroidogenic genes in adult rare minnow. The full-length cDNA encoding 11β-hydroxysteroid dehydrogenase-2 (11β-HSD2) was firstly isolated and characterized by RT-PCR and RACE methods. The gonadal transcript changes of StAR, cyp11a1, 3β-HSD, cyp17a1, 11β-HSD2 and cyp19a1a in 6-month adult Gobiocypris rarus exposed to MT and 17α-ethinylestradiol (EE2) for 7, 14 and 21 days were detected by qRT-PCR. To make an effort to connect the transcriptional changes of steroidogenic enzymes with effects on higher levels of biological organization and on VTG, one remarkable sensitive target of steroids, body and gonad weights, histology of gonads, and hepatic vtg mRNA level were measured. MT caused varying degree of abnormalities in ovaries and testes. The hepatic vtg mRNA level was highly inhibited in females and slightly altered in males by MT. Transcripts of several steroidogenic genes including StAR, cyp17a1, and cyp11a1 showed high responsiveness to MT exposure in G. rarus. The gene expression profiles of these steroidogenic genes in MT-treated groups were much distinct with the EE2-treated group.

Occurrence of 17α-ethynylestradiol (EE2) in the environment and effect on exposed biota: a review

17α-ethynylestradiol (EE2) is a synthetic hormone, which is a derivative of the natural hormone, estradiol (E2). EE2 is an orally bio-active estrogen, and is one of the most commonly used medications for humans as well as livestock and aquaculture activity. EE2 has become a widespread problem in the environment due to its high resistance to the process of degradation and its tendency to (i) absorb organic matter, (ii) accumulate in sediment and (iii) concentrate in biota. Numerous studies have reported the ability of EE2 to alter sex determination, delay sexual maturity, and decrease the secondary sexual characteristics of exposed organisms even at a low concentration (ng/L) by mimicking its natural analogue, 17β-estradiol (E2). Thus, the aim of this review is to provide an overview of the science regarding EE2, the concentration levels in the environment (water, sediment and biota) and summarize the effects of this compound on exposed biota at various concentrations, stage life, sex, and species. The challenges in respect of EE2 include the extension of the limited database on the EE2 pollution profile in the environment, its fate and transport mechanism, as well as the exposure level of EE2 for better prediction and definition revision of EE2 toxicity end points, notably for the purpose of environmental risk assessment.

Benzotriazole UV-stabilizers and benzotriazole: Antiandrogenic activity in vitro and activation of aryl hydrocarbon receptor pathway in zebrafish eleuthero-embryos

Benzotriazole UV-stabilizers (BUVs) are applied in materials for protection against UV-irradiation. They are widely used, bioaccumulate and share structural similarities to benzotriazole. Benzotriazole (1HBT) finds application as corrosion inhibitor in dishwashing detergents, antifreeze (vehicles) and aircraft de-icing agent. BUVs and 1HBT are persistent and ubiquitous in the aquatic environment, but there is little understanding of the ecotoxicological implications. Here, we comparatively analyze the hormonal activity in vitro and effects in zebrafish eleuthero-embryos in vivo. 2-(2-Hydroxy-5-methylphenyl)benzotriazole (UV-P), 2-(3-t-butyl-2-hydroxy-5-methylphenyl)-5-chlorobenzotriazole (UV-326), UV-327, UV-328, UV-329 and UV-320 showed no estrogenicity (YES assay) and androgenicity (YAS assay). However, UV-P and 1HBT showed significant antiandrogenic activity. We assessed the transcription profiles of up to 26 genes associated with different toxicological pathways in zebrafish eleuthero-embryos to elucidate potential modes of action of UV-P, UV-326 and 1HBT. Embryos were experimentally exposed for 144hpf to three measured concentrations of 15.8, 70.8, and 690μg/L UV-P, 7.5, 31.7, and 84.3μg/L UV-326 and 7.9, 97.3 and 1197.3μg/L 1HBT. Among the 26 transcripts, the induction of the aryl hydrocarbon receptor (AHR) pathway by UV-P and UV-326 was the most significant finding. UV-P led to dose-related induction of AHR1, ARNT2 and cyp1a1, as well as of phase II enzymes glutathione-S-transferase (gstp1) and ugt1a. UV-326 led to a significant induction of cyp1a1 and AHR2, but down-regulation of gstp1 at 84μg/L. Only little transcriptional alterations occurred in genes related to apoptosis, oxidative stress, hormone receptors, and steroidogenesis including aromatase. 1HBT led to only a few expressional changes at 1197μg/L. Our data lead to the conclusion that UV-P and UV-326 activate the AHR-pathway, whereas 1HBT shows only little transcriptional alterations. It should be noted, however, that effects have been observed at concentration much higher than those occurring in the environment. Forthcoming studies should show whether the observed antiandrogenic activities and transcriptional changes translate into physiological effects .

Combined effects of pharmaceuticals, personal care products, biocides and organic contaminants on the growth of Skeletonema pseudocostatum

Organisms in the environment are exposed to a number of pollutants from different compound groups. In addition to the classic pollutants like the polychlorinated biphenyls, polyaromatic hydrocarbons (PAHs), alkylphenols, biocides, etc. other compound groups of concern are constantly emerging. Pharmaceuticals and personal care products (PPCPs) can be expected to co-occur with other organic contaminants like biocides, PAHs and alkylphenols in areas affected by wastewater, industrial effluents and intensive recreational activity. In this study, representatives from these four different compound groups were tested individually and in mixtures in a growth inhibition assay with the marine algae Skeletonema pseudocostatum (formerly Skeletonema costatum) to determine whether the combined effects could be predicted by models for additive effects; the concentration addition (CA) and independent action (IA) prediction model. The eleven tested compounds reduced the growth of S. pseudocostatum in the microplate test in a concentration-dependent manner. The order of toxicity of these chemicals were irgarol>fluoxetine>diuron>benzo(a)pyrene>thioguanine>triclosan>propranolol>benzophenone 3>cetrimonium bromide>4-tert-octylphenol>endosulfan. Several binary mixtures and a mixture of eight compounds from the four different compound groups were tested. All tested mixtures were additive as model deviation ratios, the deviation between experimental and predicted effect concentrations, were within a factor of 2 from one or both prediction models (e.g. CA and IA). Interestingly, a concentration dependent shift from IA to CA, potentially due to activation of similar toxicity pathways at higher concentrations, was observed for the mixture of eight compounds. The combined effects of the multi-compound mixture were clearly additive and it should therefore be expected that PPCPs, biocides, PAHs and alkylphenols will collectively contribute to the risk in areas contaminated by such complex mixtures.

Accumulation and effects of the UV-filter octocrylene in adult and embryonic zebrafish (Danio rerio)

Wide application of the UV-filter octocrylene (OC) in cosmetics leads to contamination of the aquatic environment, but effects of OC remain unclear. Here we determine bioaccumulation and molecular effects of OC. Adult male zebrafish were exposed to 22, 209 and 383 μg/L and embryos to 69, 293 and 925 μg/L OC. OC accumulated in fish up to 17 μg/g. Calculated BCF varied between 41 and 136. Microarray analysis in brain and liver following exposure to 383 μg/L OC revealed alteration of 628 and 136 transcripts, respectively. Most prominent GO processes included developmental processes, organ development, hematopoiesis, formation of blood vessels, blood circulation, fat cell differentiation and metabolism. Validation by RT-qPCR in brain and liver of adult fish and embryos included a series of genes. Blood levels of 11-ketotestosterone were not altered. The transcriptomics data suggest that OC mainly affects transcription of genes related to developmental processes in the brain and liver as well as metabolic processes in the liver.

Exposure to 17α-ethinylestradiol decreases motility and ATP in sperm of male fighting fish Betta splendens

The synthetic estrogen 17α-ethinylestradiol (EE2) is an endocrine-disrupting chemical released into aquatic environments from sewage treatment facilities. We tested the effects of two environmentally relevant concentrations of waterborne EE2, 10 and 100 ng L(-1) , on reproductive endpoints in the teleost fish Betta splendens. In the first experiment, testes were removed from males and sperm were exposed to EE2 directly through the activation water. Direct exposure to EE2 had no effect on any measure of sperm swimming performance. In the second experiment, we exposed sexually mature male B. splendens to EE2 using a semi-static exposure protocol for 4 weeks. There were no significant treatment effects in the 10 ng L(-1) treatment group, but at the 100 ng L(-1) dose we found that fish had smaller gonads and reduced sperm swimming velocity. When allowed to interact freely with female conspecifics, males exposed to 100 ng L(-1) EE2 built smaller nests and showed a nonsignificant decrease in fertilization success. To investigate further the potential mechanism underlying the decrease in sperm quality, we repeated the chronic exposure experiment and analyzed the ATP content of sperm from fish in each treatment group. We found that males exposed to 100 ng L(-1) of EE2 had fewer moles of ATP per sperm than did fish in the other two treatment groups, suggesting that a decrease in intracellular ATP caused a reduction in sperm swimming velocity. The current study adds to the growing body of literature that indicates the risks to aquatic organisms of exposure to environmentally relevant concentrations of EE2.

Gene expression changes in aging zebrafish (Danio rerio) brains are sexually dimorphic

Background

Brain aging is a multi-factorial process due to both genetic and environmental factors. The zebrafish has recently become a popular model organism for examining aging and age-related diseases because as in humans they age gradually and exhibit cognitive decline. Few studies have examined the biological changes in the aging brain that may contribute to these declines and none have examined them within individuals with respect to gender. Our aim was to identify the main genetic pathways associated with zebrafish brain aging across gender. We chose males and females from specific age groups (young, 7.5-8.5 months and old, 31-36 months) based on the progression of cognitive decline in zebrafish. RNA was isolated from individual brains and subjected to microarray and qPCR analysis. Statistical analyses were performed using a two-way ANOVA and the relevant post-hoc tests.

Results

Our results demonstrated that in the brains of young and old male and female zebrafish there were over 500 differentially expressed genes associated with multiple pathways but most notably were those related to neurogenesis and cell differentiation, as well as brain and nervous system development.

Conclusions

The gene expression of multiple pathways is altered with age and differentially expressed in males and females. Future studies will be aimed at determining the causal relationships of age-related changes in gene expression in individual male and female brains, as well as possible interventions that counteract these alterations.

Effects of low concentrations of the antiprogestin mifepristone (RU486) in adults and embryos of zebrafish (Danio rerio): 2. Gene expression analysis and in vitro activity

Here, we analyzed the transcriptional effects of the antiprogestin mifepristone (MIF, RU486) and progesterone (P4) in zebrafish as well as their in vitro activities in yeast-based reporter gene assays. This study is associated with the reproduction study in adult zebrafish and embryos exposed for 21 days to 5, 39, 77 ng/L MIF, and 25 ng/L P4 (Blüthgen et al., 2013a). The in vitro activities of MIF and P4 were investigated using a series of recombinant yeast-based assays (YES, YAS, YPS) and compared to transcriptional alterations obtained in fish tissues and embryos from the exposure study. MIF elicited antiestrogenic, androgenic and progestogenic activities in recombinant yeast, similar to P4, and no antiprogestogenic activity in vitro. The transcriptional alterations of steroid hormone receptors were similar in adult males and females, and more pronounced in embryos. MIF tended to transcriptionally down-regulate the androgen (ar), progesterone (pgr) and glucocorticoid (gr) receptors in adult fish and embryos. Transcripts of the estrogen receptor (esr1) and vitellogenin (vtg1) were not significantly altered. A trend for down-regulation was observed for transcripts of genes belonging to steroidogenic enzymes including 17β-hydroxysteroid dehydrogenase type 3 (hsd17b3), 3 β-hydroxysteroid dehydrogenase (hsd3b), P450 aromatase A (cyp19a) and 11β-hydroxylase (cyp11b). P4 resulted in similar transcriptional alterations as MIF. The data indicate that gene expression changes (here and later gene expression is taken as synonym to gene transcription) and in vitro activities match only in part including the lack of antiprogestogenic activity of MIF. Additionally, effects on reproduction and gonad histology described in the associated report (Blüthgen et al., 2013a) can only partly be explained by gene expression data presented here.

High-throughput sequencing and pathway analysis reveal alteration of the pituitary transcriptome by 17α-ethynylestradiol (EE2) in female coho salmon, Oncorhynchus kisutch

Considerable research has been done on the effects of endocrine disrupting chemicals (EDCs) on reproduction and gene expression in the brain, liver and gonads of teleost fish, but information on impacts to the pituitary gland are still limited despite its central role in regulating reproduction. The aim of this study was to further our understanding of the potential effects of natural and synthetic estrogens on the brain-pituitary-gonad axis in fish by determining the effects of 17α-ethynylestradiol (EE2) on the pituitary transcriptome. We exposed sub-adult coho salmon (Oncorhynchus kisutch) to 0 or 12 ng EE2/L for up to 6 weeks and effects on the pituitary transcriptome of females were assessed using high-throughput Illumina(®) sequencing, RNA-Seq and pathway analysis. After 1 or 6 weeks, 218 and 670 contiguous sequences (contigs) respectively, were differentially expressed in pituitaries of EE2-exposed fish relative to control. Two of the most highly up- and down-regulated contigs were luteinizing hormone β subunit (241-fold and 395-fold at 1 and 6 weeks, respectively) and follicle-stimulating hormone β subunit (-3.4-fold at 6 weeks). Additional contigs related to gonadotropin synthesis and release were differentially expressed in EE2-exposed fish relative to controls. These included contigs involved in gonadotropin releasing hormone (GNRH) and transforming growth factor-β signaling. There was an over-representation of significantly affected contigs in 33 and 18 canonical pathways at 1 and 6 weeks, respectively, including circadian rhythm signaling, calcium signaling, peroxisome proliferator-activated receptor (PPAR) signaling, PPARα/retinoid x receptor α activation, and netrin signaling. Network analysis identified potential interactions between genes involved in circadian rhythm and GNRH signaling, suggesting possible effects of EE2 on timing of reproductive events.

Transcriptome alterations in zebrafish embryos after exposure to environmental estrogens and anti-androgens can reveal endocrine disruption

Exposure to environmental chemicals known as endocrine disruptors (EDs) is in many cases associated with an unpredictable hazard for wildlife and human health. The identification of endocrine disruptive properties of chemicals certain to enter the aquatic environment relies on toxicity tests with fish, assessing adverse effects on reproduction and sexual development. The demand for quick, reliable ED assays favored the use of fish embryos as alternative test organisms. We investigated the application of a transcriptomics-based assay for estrogenic and anti-androgenic chemicals with zebrafish embryos. Two reference compounds, 17α-ethinylestradiol and flutamide, were tested to evaluate the effects on development and the transcriptome after 48h-exposures. Comparison of the transcriptome response with other estrogenic and anti-androgenic compounds (genistein, bisphenol A, methylparaben, linuron, prochloraz, propanil) showed commonalities and differences in regulated pathways, enabling us to classify the estrogenic and anti-androgenic potencies. This demonstrates that different mechanism of ED can be assessed already in fish embryos.

Transgenic (cyp19a1b-GFP) zebrafish embryos as a tool for assessing combined effects of oestrogenic chemicals

Endocrine disrupting chemicals and especially oestrogen receptor (ER) agonists have been extensively studied over the years due to their potential effects on sexual development and reproduction in vertebrates, notably fish. As ER agonists can exist as complex mixtures in the aquatic environment, evaluating the impact of combined exposure on oestrogenic effects has become increasingly important. Use of predictive models such as concentration addition (CA) and independent action (IA) has allowed assessment of combined estrogenic effects of complex multi-compound mixtures of ER agonists in various fish in vitro and in vivo experimental models. The present work makes use of a transgenic zebrafish strain, tg(cyp19a1b-GFP), which expresses the green fluorescent protein (GFP) under the control of the cyp19a1b (brain aromatase or aromatase B) gene to determine the oestrogenic potency of ER agonists alone or in mixtures. In these studies, tg(cyp19a1b-GFP) zebrafish embryos were exposed for four days (from one to five days post fertilization) to five different oestrogenic chemicals; 17α-ethinylestradiol (EE2), 17β-estradiol (E2), estrone (E1), bisphenol A (BPA) and 4-tert-octylphenol (OP), and three mixtures of up to four of these compounds. The mixture of BPA, OP and E2 was also tested with primary cultures of rainbow trout hepatocytes by analysing the ER-mediated induction of the oestrogenic biomarker vitellogenin in order to compare the performance of the two methods for assessing oestrogenic effects of complex mixtures. The three tested mixtures were predominantly acting in an additive manner on the expression of GFP. Additivity was indicated by the overlap of the 95% confidence interval of the concentration response curves for the observed data with the CA and IA prediction models, and model deviation ratios within a factor of two for a majority of the mixture concentrations. However, minor deviations determined as more than additive effects for the mixture of EE2, E1 and E2 and less than additive effects for the mixture of BPA, OP, EE2 and E1 were observed at the higher mixture concentrations tested. The successful prediction of additivity by CA and IA in tg(cyp19a1b-GFP) zebrafish embryos and deviations at high mixture concentrations seemed to correspond well to results obtained in the rainbow trout hepatocyte assay. The present results clearly show the usefulness of combining predictive modelling and use of in vitro bioassays for rapid screening of oestrogenic effects of complex mixtures and environmental samples.

The effects of 17-α-ethinylestradiol (EE2) on molecular signaling cascades in mummichog (Fundulus heteroclitus)

Exposures to ≤10 ng/L of 17-α-ethinylestradiol (EE2) will reduce or shut down egg production in freshwater fish models, while mummichog (Fundulus heteroclitus), an estuarine species, are able to produce eggs at EE2 concentrations >3000 ng/L. The objective of this study was to gain mechanistic insight into how mummichog are able to produce eggs during exposures to high EE2. Mummichog were exposed to 0, 50 or 250 ng/L of EE2 for 14 d. There were no changes in gonadosomatic index, liversomatic index, gonad development, or plasma estradiol levels after exposure to EE2. However, testosterone significantly decreased with EE2 exposures (50, 250 ng/L). Microarray analysis in the liver revealed that cell processes associated with lipids were affected by EE2 at the transcriptome level. Based on the transcriptomics data, we hypothesize that mummichog are able to maintain lipid transport and uptake into the ovary and this may be associated with apolipoproteins, facilitating normal oocyte development. Novel gene regulatory networks for protein modification targets were also constructed to learn more about the potential roles of estrogens in the teleost liver. Although post-translational modifications (PTMs) are important regulatory mechanisms, the roles of PTMs in protein regulation in fish and the susceptibility of PTMs to aquatic pollutants are largely unexplored and may offer novel insight into mechanisms of endocrine disruption.

Molecular crosstalk between a chemical and a biological stressor and consequences on disease manifestation in rainbow trout

The aim of the present study was to examine the molecular and organism reaction of rainbow trout, Oncorhynchus mykiss, to the combined impact of two environmental stressors. The two stressors were the myxozoan parasite, Tetracapsuloides bryosalmonae, which is the etiological agent of proliferative kidney disease (PKD) and a natural stressor to salmonid populations, and 17β-estradiol (E2) as prototype of estrogen-active chemical stressors in the aquatic environment. Both stressors, the parasite and estrogenic contaminants, co-exist in Swiss rivers and are discussed as factors contributing to the decline of Swiss brown trout populations over the last decades. Using a microarray approach contrasting parasite-infected and non-infected rainbow trout at low or high estrogen levels, it was observed that molecular response patterns under joint exposure differed from those to the single stressors. More specifically, three major response patterns were present: (i) expression responses of gene transcripts to one stressor are weakened by the presence of the second stressor; (ii) expression responses of gene transcripts to one stressor are enhanced by the presence of the second stressor; (iii) expression responses of gene transcripts at joint treatment are dominated by one of the two stressors. Organism-level responses to concurrent E2 and parasite treatment - assessed through measuring parasite loads in the fish host and cumulative mortalities of trout - were dominated by the pathogen, with no modulating influence of E2. The findings reveal function- and level-specific responses of rainbow trout to stressor combinations, which are only partly predictable from the response to the single stressors.

Current metabolomics: practical applications

The field of metabolomics continues to grow rapidly over the last decade and has been proven to be a powerful technology in predicting and explaining complex phenotypes in diverse biological systems. Metabolomics complements other omics, such as transcriptomics and proteomics and since it is a 'downstream' result of gene expression, changes in the metabolome is considered to best reflect the activities of the cell at a functional level. Thus far, metabolomics might be the sole technology capable of detecting complex, biologically essential changes. As one of the omics technology, metabolomics has exciting applications in varied fields, including medical science, synthetic biology, medicine, and predictive modeling of plant, animal and microbial systems. In addition, integrated applications with genomics, transcriptomics, and proteomics provide greater understanding of global system biology. In this review, we discuss recent applications of metabolomics in microbiology, plant, animal, food, and medical science.

Gene expression profiles in the testis associated with testis-ova in adult Japanese medaka (Oryziaslatipes) exposed to 17α-ethinylestradiol

The occurrence of oocytes in the testis (testis-ova) of several fish species is often associated with exposure of estrogenic chemicals. However, induction mechanisms of the testis-ova remain to be elucidated. To develop marker genes for detecting testis-ova in the testis, adult male medaka were exposed to nominal concentration of 100 ng L(-1) of 17α-ethinylestradiol (EE2) for 3-5 weeks, and 800 ng estradiol benzoate (EB) for 3 weeks (experiment I), and a measured concentration of 20 ng L(-1) EE2 for 1-6 weeks (experiment II). Histological analysis was performed for the testis, and microarray analyses were performed for the testis, liver and brain. Microarray analysis in the estrogen-exposed medaka liver showed vitellogenin and choriogenin as estrogen responsive genes. Testis-ova were induced in the testis after 4 weeks of exposure to 100 ng L(-1) EE2, 3 weeks of exposure to 800 ng EB, and 6 weeks of exposure to 20 ng L(-1) EE2. Microarray analysis of estrogen-exposed testes revealed up-regulation of genes related to zona pellucida (ZP) and the oocytes marker gene, 42Sp50. Using quantitative RT-PCR we confirmed that Zpc5 gene can be used as a marker for the detection of testis-ova in male medaka.

Progestins and antiprogestins affect gene expression in early development in zebrafish (Danio rerio) at environmental concentrations

Progesterone (P4) and synthetic progestins (gestagens) from contraceptives and hormone therapy occur in treated wastewater and surface water, and they may have adverse effects on aquatic organisms. Little is known about the molecular and reproductive effects of P4 and synthetic progestins in fish, and effects of the antiprogestin mifepristone (RU486, an abortive) are unknown. We aimed at elucidating effects on the hormone system by quantitative determination of transcriptional changes of target genes induced by 2, 20, and 200 ng/L P4, RU486, norethindrone (NET), and levonorgestrel (LNG). We exposed zebrafish embryos for 144 h post fertilization (hpf) to these compounds and analyzed expressional changes of ar, esr1, vtg1, hsd17ß3, and progesterone (pgr), mineralo- (mr), and glucocorticoid (gr) receptors, each at 48, 96, and 144 hpf. Concentrations of NET and LNG were constant during exposure, while P4 and RU486 decreased. P4 and RU486 were the most potent steroids. Significant up to 4-fold induction of pgr, ar, mr, and hsd17b3 occurred at 2 ng/L P4 and higher, while RU484 inhibited pgr expression. NET and LNG modulated some transcripts mainly above 2 ng/L. The expressional chances occurring at environmental levels may translate to negative interference with differentiation of brain and gonads, and consequently reproduction.

Insight into molecular pathways of retinal metabolism, associated with vitellogenesis in zebrafish

Retinal is the main retinoid stored in oviparous eggs of fish, amphibians, and reptiles, reaching the oocytes in association with vitellogenins, the yolk precursor proteins. During early presegmentation stages of zebrafish embryos, retinal is metabolized to retinoic acid (RA), which regulates genes involved in cell proliferation, differentiation, and tissue function and is therefore essential for normal embryonic development. While synthesis of vitellogenin and its regulation by 17β-estradiol (E(2)) were extensively investigated, pathways for retinal synthesis remain obscure. We determined the expression pattern of 46 candidate genes, aiming at identifying enzymes associated with retinal synthesis, ascertaining whether they were regulated by E(2), and finding pathways that could fulfill the demand for retinoids during vitellogenesis. Genes associated with retinal synthesis were upregulated in liver (rdh10, rdh13, sdr) and surprisingly also in intestine (rdh13) and ovary (rdh1, sdr), concomitantly with higher gene expression and synthesis of vitellogenins in liver but also in extrahepatic tissues, shown here for the first time. Vitellogenin synthesis in the ovary was regulated by E(2). Gene expression studies suggest that elevated retinal synthesis in liver, intestine, and ovary also depends on cleavage of carotenoids (by Bcdo2 or Bmco1), but in the ovary it may also be contingent on higher uptake of retinol from the circulatory system (via Stra6) and retinol synthesis from retinyl esters (by Lpl). Decrease in oxidation (by Raldh2 or Raldh3) of retinal to RA and/or degradation of RA (by Cyp26a1) may also facilitate higher hepatic retinal levels. Together, these processes enable meeting the putative demands of retinal for binding to vitellogenins. Bioinformatic tools reveal multiple hormone response elements in the studied genes, suggesting complex and intricate regulation of these processes.