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

Multiomics analysis reveal the impact of 17α-Ethinylestradiol on mortality in juvenile zebrafish

17α-Ethinylestradiol (EE2) is known for its endocrine-disrupting effects on embryonic and adult fish. However, its impact on juvenile zebrafish has not been well established. In this study, juvenile zebrafish were exposed to EE2 at concentrations of 5 ng/L (low dose, L), 10 ng/L (medium dose, M), and 50 ng/L (high dose, H) from 21 days post-fertilization (dpf) to 49 dpf. We assessed their growth, development, behavior, transcriptome, and metabolome. The findings showed that the survival rate in the EE2-H group was 66.8 %, with all surviving fish displaying stunted growth and swollen, transparent abdomens by 49 dpf. Moreover, severe organ deformities were observed in the gills, kidneys, intestines, and heart of fish in both the EE2-H and EE2-M groups. Co-expression analysis of mRNA and lncRNA revealed that EE2 downregulated the transcription of key genes involved in the cell cycle, DNA replication, and Fanconi anemia signaling pathways. Additionally, metabolomic analysis indicated that EE2 influenced metabolism and development-related signaling pathways. These pathways were also significantly identified based on the genes regulated by lncRNA. Consequently, EE2 induced organ deformities and mortality in juvenile zebrafish by disrupting signaling pathways associated with development and metabolism. The results of this study offer new mechanistic insights into the adverse effects of EE2 on juvenile zebrafish based on multiomics analysis. The juvenile zebrafish are highly sensitive to EE2 exposure, which is not limited to adult and embryonic stages. It is a potential model for studying developmental toxicity.

Influence of EE2 exposure, age and sex on telomere length in European long-snouted seahorse (Hippocampus guttulatus)

After a Telomere Lengthening in juvenile stage, a progressive telomere shortening occurs with age despite higher telomerase level. Telomere Length (TL) may also reflect past physiological state such as a chronic chemical stress. Several studies have revealed a correlation between TL, ageing and/or sex in vertebrates, including teleosts; however, the patterns of telomere dynamics with telomerase mRNA expression, sex, lifespan or chemical stress in teleosts are unclear. The first aim of this study is to verify if telomere length is age and sex-dependent. The second aim is to consider if TL is a useful indicator of stress response in European long-snouted seahorse, Hippocampus guttulatus, an ectothermic and non-model system. We showed that after telomere lengthening during the juvenile stage, a telomeric attrition became significant in sexually mature individuals (p = 0.042). TL decreased in older seahorses despite the presence of somatic telomerase mRNA expression at all life stages studied. There was no difference in TL between males and females, but telomerase mRNA expression was consistently higher in females than males. Exposure to EE2 had no effect on TL in young seahorses, but was correlated with a significant increase in telomerase mRNA expression and various physiological disruptions. Here, a growth retardation of -10 % for body length (p = 0.016) and approximately -45 % for mass (p = 0.001) compared to healthy juvenile seahorses was observed. Our data suggest that telomere dynamics alone should not be used as a marker of EE2 exposure in juvenile seahorses.

Multiple biomarker responses in female Clarias gariepinus exposed to acetaminophen

Several authors have documented the presences of acetaminophen (APAP) in both surface and groundwater and have received attention from government agencies and basic authorities across the globe. The impacts of such pharmaceutical products on non-target organism like fish are underestimated as a result of selected investigation using few biomarkers. We evaluated the sub-chronic impacts of APAP in female catfish (Clarias gariepinus) using multiple biomarkers. The exposure of female catfish to APAP induced oxidative stress. Markers such as superoxide dismutase (SOD), glutathione peroxidase (GPx), and total antioxidant capacity (TAC) were significantly higher in all exposed groups. Exposure of Clarias gariepinus to APAPA caused histological alterations in the gills (fusion and shortening of some filaments, hyperplasia of the epithelial gill cells, aneurism, congestion, and epithelial rupture of the gills), liver (apoptotic hyperplasia, sinusoidal congestion, and necrosis of the hepatocytes), and gonad (degenerated follicles and ovarian apoptosis). Furthermore, multivariate results indicated that there was a distinct response from the acetaminophen-exposed female catfish, with over 95% of the biomarkers significantly contributing to the discrimination between the acetaminophen-exposed female catfish and the control groups. Our research provides evidence supporting the use of a multiple biomarker approach to evaluate the impacts of drugs on the health status of exposed fish.

Bioinformatics approaches and big data analytics opportunities in improving fisheries and aquaculture

Aquaculture has witnessed an excellent growth rate during the last two decades and offers huge potential to provide nutritional as well as livelihood security. Genomic research has contributed significantly toward the development of beneficial technologies for aquaculture. The existing high throughput technologies like next-generation technologies generate oceanic data which requires extensive analysis using appropriate tools. Bioinformatics is a rapidly evolving science that involves integrating gene based information and computational technology to produce new knowledge for the benefit of aquaculture. Bioinformatics provides new opportunities as well as challenges for information and data processing in new generation aquaculture. Rapid technical advancements have opened up a world of possibilities for using current genomics to improve aquaculture performance. Understanding the genes that govern economically relevant characteristics, necessitates a significant amount of additional research. The various dimensions of data sources includes next-generation DNA sequencing, protein sequencing, RNA sequencing gene expression profiles, metabolic pathways, molecular markers, and so on. Appropriate bioinformatics tools are developed to mine the biologically relevant and commercially useful results. The purpose of this scoping review is to present various arms of diverse bioinformatics tools with special emphasis on practical translation to the aquaculture industry.

Synthetic estrogen bioaccumulates and changes the behavior and biochemical biomarkers in adult zebrafish

Estrogen is considered to be an endocrine disrupter and is becoming increasingly more prevalent in the daily life of humans. In some cases, estrogen is not fully metabolized by organisms and may be excreted in either its original form or in organic complex forms, into water residue systems reaching concentrations of 0.05 ng.L^-1 to 75 ng.L^-1. However, estrogen 17α-ethinylestradiol (EE2), which is used in oral contraceptives, is very difficult to remove from water. Here, we evaluated whether the synthetic hormone, EE2, affects the nervous system and the behavior of adult zebrafish. We established a range of concentrations (0.05, 0.5, 5, 50, and 75 ng.L^-1), in addition to the control, to evaluate the effect of this compound and its bioaccumulation in zebrafish tissues. Here we show that EE2 bioaccumulates in fish and can change its behavior with an increased time in the upper zone (novel tank test) and far from the shoal segment (social preference test), demonstrating a clear anxiolytic pattern. The anxiolytic effect of EE2 can be harmful as it can affect the stress response of the species. The results presented herein reinforce the idea that the presence of EE2 in environmental water can be dangerous for non-target animals.

Global expression response of genes in sex-undifferentiated Nile tilapia gonads after exposure to trace letrozole

The aromatase inhibitor letrozole can be found in rivers, effluents, and even drinking water. Studies have demonstrated that letrozole affects various metabolic pathways and may cause reproductive toxicity, especially in fish exposed during development. However, studies on the effect of a low concentration of letrozole at the whole-gonad transcriptomic level in the early stage of fish sexual development have not been investigated. The aim of our study was to explore the potential effects of a low concentration of letrozole on the gonad transcriptome of Nile tilapia at an early stage of sexual development. In this study, 9 dpf (days postfertilization) Nile tilapia were exposed to trace letrozole for 12 days. Letrozole exposure from 9 dpf to 21 dpf persistently altered phenotypic sex development and induced the male-biased sex ratio. The transcriptome results showed that 1173 differentially expressed genes (DEGs) were present in the female control vs 1.5 μg/L letrozole-treated female comparison group and that 1576 DEGs were present in the 1.5 μg/L letrozole-treated female vs male control comparison group. Differentially expressed gene enrichment analysis revealed several crucial pathways, including the drug metabolism-cytochrome P450 pathway, the ErbB-PI3K/Akt/mTOR pathway, and the calcium signalling pathway. Further analysis of these identified DEGs indicated that some key genes correlated with metabolism and epigenetic regulation were significantly affected by letrozole, such as UDP-glucuronosyltransferase (Ugt), glutathione S-transferase omega-1 (Gsto1), lysine-specific demethylase 6bb (Kdm6bb, original name is Kdm6a), jumonji and AT-rich interaction domain containing 2 (Jarid2b, original name is Jarid2), growth arrest and DNA damage inducible gamma (Gadd45g), and chromobox protein 7 (Cbx7). The qRT-PCR validation results for twelve DEGs showed that the Pearson's correlation of the log₁₀^{fold change} values between the qPCR and RNA-Seq results was 0.90, indicating the accuracy and reliability of the RNA-Seq results. Our study is the first to report the effect of letrozole on the transcriptome of gonads from fish during early-stage sexual development. These findings will be useful for understanding the toxic effects and molecular mechanisms of letrozole exposure at the early stage of gonad development on the sexual development of aquatic organisms.

Cross-species physiological interactions of endocrine disrupting chemicals with the circadian clock

Endocrine disrupting chemicals (EDCs) are endocrine-active chemical pollutants that disrupt reproductive, neuroendocrine, cardiovascular and metabolic health across species. The circadian clock is a transcriptional oscillator responsible for entraining 24-hour rhythms of physiology, behavior and metabolism. Extensive bidirectional cross talk exists between circadian and endocrine systems and circadian rhythmicity is present at all levels of endocrine control, from synthesis and release of hormones, to sensitivity of target tissues to hormone action. In mammals, a range of hormones directly alter clock gene expression and circadian physiology via nuclear receptor (NR) binding and subsequent genomic action, modulating physiological processes such as nutrient and energy metabolism, stress response, reproductive physiology and circadian behavioral rhythms. The potential for EDCs to perturb circadian clocks or circadian-driven physiology is not well characterized. For this reason, we explore evidence for parallel endocrine and circadian disruption following EDC exposure across species. In the reviewed studies, EDCs dysregulated core clock and circadian rhythm network gene expression in brain and peripheral organs, and altered circadian reproductive, behavioral and metabolic rhythms. Circadian impacts occurred in parallel to endocrine and metabolic alterations such as impaired fertility and dysregulated metabolic and energetic homeostasis. Further research is warranted to understand the nature of interaction between circadian and endocrine systems in mediating physiological effects of EDC exposure at environmental levels.

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

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

Multiple Xenosteroid Pollutants Biomarker Changes in Xultured Nile Tilapia Using Wastewater Effluents as Their Primary Water Source

Simple Summary

Estrogenic endocrine disruptive chemicals (E-EDCs) are important types of pollutants in fish farms worldwide and a globally concerned problem. In this study, Nile tilapia fish farms receiving wastewater effluents in Egypt were selected as highly, moderately polluted fish farms; besides, a putative control site was deemed low in contamination. Levels of E-EDCs (natural and synthetic steroids, and industrial phenolic compound/bisphenol A (BPA)) was recorded in farm water, and fish tissues at all sites under consideration, mainly, lower levels of testosterone, progesterone, zeranol, and 17β-estradiol were detected compared to the higher level of BPA. Moreover, the effects of these pollutants on fish biometric, reproductive genes, and hormonal biomarkers was evaluated along with the observed associated histopathological alterations. Our findings revealed the detection of some steroidal compounds with a higher level of the BPA. All analyzed biomarkers were reduced to a great extent in the highly polluted sites compared to others, and the histopathological alterations observed were supportive of other measurements. These observations warrant strict monitoring of aquatic pollution sources and the development of strategic plans to control aquaculture pollution.

Abstract

This study was undertaken to screen levels of xenosteroids (estrogenic endocrine disrupting chemicals/E-EDCs) in Nile tilapia (Oreochromis niloticus) fish farms subjected to water fill from the drain at three sites S1 (highly polluted), S2 (moderately polluted), and a putative reference site (RS). Biometric, hormonal, gene expression, and histopathological analysis were investigated. Testosterone, progesterone, and zeranol residues were detected at (0.12–3.44 µg/L) in water samples of different sites. Bisphenol-A (BPA) exhibited a very high concentration (6.5 µg/mL) in water samples from S1. Testosterone, 17β-estradiol residues were detected in fish tissues from all sites at (0.16–3.8 µg/Kg) and (1.05–5.01 µg/Kg), respectively. BPA residues were detected at a very high concentration in the liver and muscle of fish collected from S1 at higher levels of 25.9 and 48.07 µg/Kg, respectively. The detected E-EDCs, at different sites, particularly BPA, reduced the somatic and testicular growth among sites and oversampling time points. Meanwhile, hepatosomatic index (HSI) was significantly increased in S1 compared to S2. All analyzed genes estrogen receptor-type I (er-I, er-ɑ) and II (er-II, er-ß1), polypeptide 1a (cyp19a1), SRY-box containing gene 9 (sox9), and vitellogenin (vtg) and gonadotropin hormones (luteinizing hormone (LH), follicle-stimulating hormone (FSH)), testosterone, 17β-estradiol, and anti-Mullerian hormone (AMH) were significantly expressed at S1 compared to other sites. Histopathology was more evident in S1 than other sites. These findings warrant immediate strategies development to control aquatic pollution and maintain fish welfare and aquaculture sustainability.

Effects of new generation progestins, including as mixtures and in combination with other classes of steroid hormones, on zebrafish early life stages

Fish are exposed to progestins and steroid mixtures in contaminated waters but the ecotoxicological implications are not sufficiently known. Here we analyze effects of the new generation progestin dienogest (DNG) followed by investigating effects of mixtures of new generation progestins containing DNG, cyproterone acetate and drospirenone and the hormone progesterone. Furthermore, effects of this mixture were studied after adding 17β-estradiol (E2) and clobetasol propionate (CLO) in zebrafish embryos and larvae at concentrations between 0.01 and 10 μg/L. DNG showed only very minor transcriptional alterations among the 24 assessed genes with downregulation of the fshb transcript only. The progestin mixture caused weak induction of the lhb, cyp2k22 and sult2st3 transcripts. Addition of E2 to the mixture caused strong induction vtg1, cyp19b, esr1 and lhb, as well as downregulation of fshb from 0.01 μg/L onwards. Besides altering the same transcripts, addition of CLO altered glucocorticoid regulated genes mmp-9, mmp-13, g6pca, fkbp5 and irg1l. While each steroid class exhibited its specific activity independently in the mixture, sult2st3 and cyp2k22 were regulated by both E2 and CLO. Furthermore, CLO alone and in mixtures decreased spontaneous muscle contractions, increased heartrate and induced edema. Our study highlights the prominent effects of E2 and CLO in environmental steroid mixtures, while new generation progestins show relatively low activity.

Twenty years of transcriptomics, 17alpha-ethinylestradiol, and fish

In aquatic toxicology, perhaps no pharmaceutical has been investigated more intensely than 17alpha-ethinylestradiol (EE2), the active ingredient of the birth control pill. At the turn of the century, the fields of comparative endocrinology and endocrine disruption research witnessed the emergence of omics technologies, which were rapidly adapted to characterize potential hazards associated with exposures to environmental estrogens, such as EE2. Since then, significant advances have been made by the scientific community, and as a result, much has been learned about estrogen receptor signaling in fish from environmental xenoestrogens. Vitellogenin, the egg yolk precursor protein, was identified as a major estrogen-responsive gene, establishing itself as the premier biomarker for estrogenic exposures. Omics studies have identified a plethora of estrogen responsive genes, contributing to a wealth of knowledge on estrogen-mediated regulatory networks in teleosts. There have been ~40 studies that report on transcriptome responses to EE2 in a variety of fish species (e.g., zebrafish, fathead minnows, rainbow trout, pipefish, mummichog, stickleback, cod, and others). Data on the liver and testis transcriptomes dominate in the literature and have been the subject of many EE2 studies, yet there remain knowledge gaps for other tissues, such as the spleen, kidney, and pituitary. Inter-laboratory genomics studies have revealed transcriptional networks altered by EE2 treatment in the liver; networks related to amino acid activation and protein folding are increased by EE2 while those related to xenobiotic metabolism, immune system, circulation, and triglyceride storage are suppressed. EE2-responsive networks in other tissues are not as comprehensively defined which is a knowledge gap as regulated networks are expected to be tissue-specific. On the horizon, omics studies for estrogen-mediated effects in fish include: (1) Establishing conceptual frameworks for incorporating estrogen-responsive networks into environmental monitoring programs; (2) Leveraging in vitro and computational toxicology approaches to identify chemicals associated with estrogen receptor-mediated effects in fish (e.g., male vitellogenin production); (3) Discovering new tissue-specific estrogen receptor signaling pathways in fish; and (4) Developing quantitative adverse outcome pathway predictive models for estrogen signaling. As we look ahead, research into EE2 over the past several decades can serve as a template for the array of hormones and endocrine active substances yet to be fully characterized or discovered.

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 vitro exposure of vitellogenic rainbow trout ovarian follicles to endocrine disrupting chemicals can alter basal estradiol-17β production and responsiveness to a gonadotropin challenge

Endogenous estrogens play major roles in many aspects of female reproductive development in fish. In order to develop a relatively high-throughput assay to determine the potential impact on reproductive development, vitellogenic rainbow trout ovarian follicles were exposed to a suite of contaminants in vitro and then assessed for the ability to produce estradiol-17β (E2) after a 500 ng/ml salmon gonadotropin (sGTH) challenge. There was a positive correlation between ovarian follicle size and E2 production, but an inverse correlation between size and responsiveness to sGTH. Significant impacts on E2 levels were observed following treatment with different endocrine disrupting chemicals, such as 17α-ethinylestradiol (EE2), prochloraz, or trenbolone. EE2 was remarkably potent and significantly reduced ovarian follicle responsiveness to sGTH at concentrations as low as 0.1 nM. Of the other contaminants tested, only tamoxifen impacted E2 levels, and only at concentrations near the limits of solubility. Flutamide, fluoxetine, 4-hydroxy tamoxifen, hydroxyflutamide, and norfluoxetine had little or no impact. Quantitative PCR analyses of steroidogenesis-related genes were carried out on EE2 treated ovarian follicles, but significant transcriptional responses to EE2 were not observed. Overall, this study suggests that xenoestrogens and anti-estrogens are more likely to interfere with ovarian E2 synthesis than other classes of EDCs. This also provides a template for further testing of the effects of EDCs on ovarian function.

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.

Initiation of sex change and gonadal gene expression in black sea bass (Centropristis striata) exposed to exemestane, an aromatase inhibitor

Many teleost fishes exhibit sequential hermaphroditism, where male or female gonads develop first and later undergo sex change. Model sex change species are characterized by social hierarchies and coloration changes, which enable experimental manipulations to better understand these processes. However, other species such as the protogynous black sea bass (Centropristis striata) do not exhibit these characteristics and instead receive research attention due to their importance in fisheries or aquaculture. Black sea bass social structure is unknown, which makes sex change sampling difficult, and few molecular resources are available. The purpose of the present study was to induce sex change using exemestane, an aromatase inhibitor, and assess gonadal gene expression using sex markers (amh, zpc2) and genes involved in steroidogenesis (cyp19a1a, cyp11b), estrogen signaling (esr1, esr2b), and apoptosis or atresia (aen, casp9, fabp11, parg, pdcd4, rif1). Overall, dietary exemestane treatment was effective, and most exposed females exhibited early histological signs of sex change and significantly higher rates of ovarian atresia relative to control females. Genes associated with atresia did not reflect this, however, as expression patterns in sex changing gonads were overall similar to those of ovaries, likely due to a whole ovary dilution effect of the RNA. Still, small but insignificant expression decreases during early sex change were detected for ovary-related genes (aen, casp9, fabp11, zpc2) and anti-apoptotic factors (parg, rif1). Exemestane treatment did not impact spermatogenesis or testicular gene expression, but testes were generally characterized by elevated steroidogenic enzyme and estrogen receptor mRNAs. Further research will be needed to understand these processes in black sea bass, using isolated ovarian follicles and multiple stages of sex change.

In vivo effects of 17α-ethinylestradiol, 17β-estradiol and 4-nonylphenol on insulin-like growth-factor binding proteins (igfbps) in Atlantic salmon

Feminizing endocrine disrupting compounds (EDCs) affect the growth and development of teleost fishes. The major regulator of growth performance, the growth hormone (Gh)/insulin-like growth-factor (Igf) system, is sensitive to estrogenic compounds and mediates certain physiological and potentially behavioral consequences of EDC exposure. Igf binding proteins (Igfbps) are key modulators of Igf activity, but their alteration by EDCs has not been examined. We investigated two life-stages (fry and smolts) of Atlantic salmon (Salmo salar), and characterized how the Gh/Igf/Igfbp system responded to waterborne 17α-ethinylestradiol (EE₂), 17β-estradiol (E₂) and 4-nonylphenol (NP). Fry exposed to EE₂ and NP for 21 days had increased hepatic vitellogenin (vtg) mRNA levels while hepatic estrogen receptor α (erα), gh receptor (ghr), igf1 and igf2 mRNA levels were decreased. NP-exposed fry had reduced body mass and total length compared to controls. EE₂ and NP reduced hepatic igfbp1b1, -2a, -2b1, -4, -5b2 and -6b1, and stimulated igfbp5a. In smolts, hepatic vtg mRNA levels were induced following 4-day exposures to all three EDCs, while erα only responded to EE₂ and E₂. EDC exposures did not affect body mass or fork length; however, EE₂ diminished plasma Gh and Igf1 levels in parallel with reductions in hepatic ghr and igf1. In smolts, EE₂ and E₂ diminished hepatic igfbp1b1, -4 and -6b1, and stimulated igfbp5a. There were no signs of compromised ionoregulation in smolts, as indicated by unchanged branchial ion pump/transporter mRNA levels. We conclude that hepatic igfbps respond (directly and/or indirectly) to environmental estrogens during two key life-stages of Atlantic salmon, and thus may modulate the growth and development of exposed individuals.

Early phosphoproteomic changes for adverse outcome pathway development in the fathead minnow (Pimephales promelas) brain

Adverse outcome pathways (AOPs) are conceptual frameworks that organize and link contaminant-induced mechanistic molecular changes to adverse biological responses at the individual and population level. AOPs leverage molecular and high content mechanistic information for regulatory decision-making, but most current AOPs for hormonally active agents (HAAs) focus on nuclear receptor-mediated effects only despite the overwhelming evidence that HAAs also activate membrane receptors. Activation of membrane receptors triggers non-genomic signaling cascades often transduced by protein phosphorylation leading to phenotypic changes. We utilized label-free LC-MS/MS to identify proteins differentially phosphorylated in the brain of fathead minnows (Pimephales promelas) aqueously exposed for 30 minutes to two HAAs, 17α-ethinylestradiol (EE2), a strong estrogenic substance, and levonorgestrel (LNG), a progestin, both components of the birth control pill. EE2 promoted differential phosphorylation of proteins involved in neuronal processes such as nervous system development, synaptic transmission, and neuroprotection, while LNG induced differential phosphorylation of proteins involved in axon cargo transport and calcium ion homeostasis. EE2 and LNG caused similar enrichment of synaptic plasticity and neurogenesis. This study is the first to identify molecular changes in vivo in fish after short-term exposure and highlights transduction of rapid signaling mechanisms as targets of HAAs, in addition to nuclear receptor-mediated pathways.

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.

Transcriptional alterations induced by binary mixtures of ethinylestradiol and norgestrel during the early development of zebrafish (Danio rerio)

Synthetic estrogens and progestins are commonly used in human and veterinary medicine. After use, they reach aquatic environments via discharge of wastewaters from human and animals, thus posing potential risks to organisms. So far, very little is known about their combined effects in aquatic organisms. The aim of this study was to investigate the effects of binary mixtures of ethinylestradiol (EE2) and norgestrel (NGT) on embryonic zebrafish (Danio rerio) by measuring transcriptional alterations. Zebrafish embryos were exposed to EE2 and NGT alone or in combination at concentrations between 36 and 5513ngL^-1 for 96h post-fertilization (hpf). The results showed that most of gene transcriptions of hypothalamic-pituitary-gonadal axis (e.g., Pgr, Mprα, Esr1, Esr2a, Vtg1, Ar, Cyp11b, Star, Gnrh3 and Fshb) and circadian rhythm signaling (e.g., Cry1a, Cry2a, Cry2b, Per3, Arntl1b, Arntl2, Clock1a, Cry3 and Cry4) displayed most pronounced alterations in the mixtures as compared to single EE2 and NGT exposures. This finding suggests exposure to the binary mixtures of EE2 and NGT produced significantly enhanced effects in fish as compared to single chemical exposures, and their coexistence could have significant environmental implications.

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.

Transcriptomic alterations in the brain of painted turtles (Chrysemys picta) developmentally exposed to bisphenol A or ethinyl estradiol

Developmental exposure of turtles and other reptiles to endocrine-disrupting chemicals (EDCs), including bisphenol A (BPA) and ethinyl estradiol (EE), can stimulate partial to full gonadal sex-reversal in males. We have also recently shown that in ovo exposure to either EDC can induce similar sex-dependent behavioral changes typified by improved spatial learning and memory or possibly feminized brain responses. Observed behavioral changes are presumed to be due to BPA- and EE-induced brain transcriptomic alterations during development. To test this hypothesis, we treated painted turtles (Chrysemys picta) at developmental stage 17, incubated at 26°C (male-inducing temperature), with 1) BPA (1 ng/µl), 2) EE (4 ng/µl), or 3) vehicle ethanol (control group). Ten months after hatching and completion of the behavioral tests, juvenile turtles were euthanized, brains were collected and frozen in liquid nitrogen, and RNA was isolated for RNA-Seq analysis. Turtles exposed to BPA clustered separately from EE-exposed and control individuals. More transcripts and gene pathways were altered in BPA vs. EE individuals. The one transcript upregulated in both BPA- and EE-exposed individuals was the mitochondrial-associated gene, ND5, which is involved in oxidative phosphorylation. Early exposure of turtles to BPA increases transcripts linked with ribosomal and mitochondrial functions, especially bioenergetics, which has been previously linked with improved cognitive performance. In summary, even though both BPA and EE resulted in similar behavioral alterations, they diverge in the pattern of neural transcript alterations with early BPA significantly upregulating several genes involved in oxidative phosphorylation, mitochondrial activity, and ribosomal function, which could enhance cognitive performance.

Effects of antiandrogenic progestins, chlormadinone and cyproterone acetate, and the estrogen 17α-ethinylestradiol (EE2), and their mixtures: Transactivation with human and rainbowfish hormone receptors and transcriptional effects in zebrafish (Danio rerio) eleuthero-embryos

Synthetic progestins act as endocrine disrupters in fish but their risk to the environment is not sufficiently known. Here, we focused on an unexplored antiandrogenic progestin, chlormadinone acetate (CMA), and the antiandrogenic progestin cyproterone acetate (CPA). The aim was to evaluate whether their in vitro interaction with human and rainbowfish (Melanotaenia fluviatilis) sex hormone receptors is similar. Furthermore, we investigated their activity in zebrafish (Danio rerio) eleuthero-embryos. First, we studied agonistic and antagonistic activities of CMA, CPA, and 17α-ethinylestradiol (EE2), in recombinant yeast expressing either the human progesterone (PGR), androgen (AR), or estrogen receptor. The same compounds were also investigated in vitro in a stable transfection cell system expressing rainbowfish nuclear steroid receptors. For human receptors, both progestins exhibited progestogenic, androgenic and antiestrogenic activity with no antiandrogenic or estrogenic activity. In contrast, interactions with rainbowfish receptors showed no progestogenic, but antiandrogenic, antiglucocorticoid, and some antiestrogenic activity. Thus, interaction with and transactivation of human and rainbowfish PGR and AR were distinctly different. Second, we analyzed transcriptional alterations in zebrafish eleuthero-embryos at 96 and 144h post fertilization after exposure to CPA, CMA, EE2, and binary mixtures of CMA and CPA with EE2, mimicking the use in oral contraceptives. CMA led to slight down-regulation of the ar transcript, while CPA down-regulated ar and pgr transcripts. EE2 exposure resulted in significant transcriptional alterations of several genes, including esr1, pgr, vtg1, cyp19b, and gonadotropins (fshb, lhb). The mixture activity of CMA and EE2 followed the independent action model, while CPA and EE2 mixtures showed additive action in transcriptional alterations. Third, we analyzed the interactions of binary mixtures of CMA and CPA, and of CMA and EE2 for their joint activity in vitro and in eleuthero-embryos. Both mixtures behaved according to the concentration addition model in their in vitro interaction with human and rainbowfish receptors, often showing antagonism. In zebrafish eleuthero-embryos, binary mixtures of CMA and EE2 showed the same expression patterns as EE2 alone, indicating an independent action in vivo. Our study demonstrates that CMA and CPA interact distinctly with human and rainbowfish receptors, suggesting that activities of these and possibly additional environmental steroids determined with yeast expressing human receptors cannot simply be translated to fish. The lack of agonistic activities of both progestins to rainbowfish PGR and AR is the probable reason for the low activity found in zebrafish eleuthero-embryos.

Wastewater treatment plant effluent alters pituitary gland gonadotropin mRNA levels in juvenile coho salmon (Oncorhynchus kisutch)

It is well known that endocrine disrupting compounds (EDCs) present in wastewater treatment plant (WWTP) effluents interfere with reproduction in fish, including altered gonad development and induction of vitellogenin (Vtg), a female-specific egg yolk protein precursor produced in the liver. As a result, studies have focused on the effects of EDC exposure on the gonad and liver. However, impacts of environmental EDC exposure at higher levels of the hypothalamic-pituitary-gonad axis are less well understood. The pituitary gonadotropins, follicle-stimulating hormone (Fsh) and luteinizing hormone (Lh) are involved in all aspects of gonad development and are subject to feedback from gonadal steroids making them a likely target of endocrine disruption. In this study, the effects of WWTP effluent exposure on pituitary gonadotropin mRNA expression were investigated to assess the utility of Lh beta-subunit (lhb) as a biomarker of estrogen exposure in juvenile coho salmon (Oncorhynchus kisutch). First, a controlled 72-h exposure to 17α-ethynylestradiol (EE2) and 17β-trenbolone (TREN) was performed to evaluate the response of juvenile coho salmon to EDC exposure. Second, juvenile coho salmon were exposed to 0, 20 or 100% effluent from eight WWTPs from the Puget Sound, WA region for 72h. Juvenile coho salmon exposed to 2 and 10ng EE2L(-1) had 17-fold and 215-fold higher lhb mRNA levels relative to control fish. Hepatic vtg mRNA levels were dramatically increased 6670-fold, but only in response to 10ng EE2L(-1) and Fsh beta-subunit (fshb) mRNA levels were not altered by any of the treatments. In the WWTP effluent exposures, lhb mRNA levels were significantly elevated in fish exposed to five of the WWTP effluents. In contrast, transcript levels of vtg were not affected by any of the WWTP effluent exposures. Mean levels of natural and synthetic estrogens in fish bile were consistent with pituitary lhb expression, suggesting that the observed lhb induction may be due to estrogenic activity of the WWTP effluents. These results suggest that lhb gene expression may be a sensitive index of acute exposure to estrogenic chemicals in juvenile coho salmon. Further work is needed to determine the kinetics and specificity of lhb induction to evaluate its utility as a potential indicator of estrogen exposure in immature fish.

High-throughput transcriptome sequencing reveals the combined effects of key e-waste contaminants, decabromodiphenyl ether (BDE-209) and lead, in zebrafish larvae

PBDEs and heavy metals are two major contaminants at e-waste disposal sites, but their combined effects remain largely unexplored. In the present study, the transcriptomic profiles of zebrafish larvae were examined after acute exposure of embryos to 200 μg/L BDE-209, 20 μg/L lead (Pb) or their mixture (Mix). Stimulation of steroidogenic pathway and vitellogenesis in the BDE-209 and Mix treatments indicated the estrogenic activities of BDE-209, while Pb antagonized those estrogenic effects in the Mix treatment. Increased heart rates were observed in zebrafish exposed to the Pb and Mix treatments. The cardiac dysfunction probably resulted from the promotion of angiogenesis, increased adrenergic drive and induction of the formation of blood clot. Furthermore, the Pb and Mix treatments activated neuroendocrine regulation of the pituitary in a positive feedback loop, via the thyrotropin-releasing hormone receptor, thus increasing thyroid hormone production self-adaptively. Overall, the interaction between BDE-209 and Pb led to synergistic and antagonistic effects on gene transcriptions, with concerted contribution from their individual toxicological properties.

De novo characterization of the liver transcriptome of javelin goby Synechogobius hasta and analysis of its transcriptomic profile following waterborne copper exposure

Previous studies have investigated the physiological responses to chronic copper (Cu) exposure in the liver of Synechogobius hasta; however, little information is available on the underlying molecular mechanisms. In an effort to better understand the mechanisms of Cu toxicity and to illuminate global gene expression patterns modulated by Cu exposure, we obtained the liver transcriptome information of S. hasta by RNA sequencing (RNA-seq) technology and also investigated the differential expression of genes following waterborne Cu exposure. Using the Illumina sequencing platform, as many as 60,217 unigenes were generated, with 815 bp of average length and 1298 bp of unigene N50 after filtering and assembly. For functional annotation analysis, 34,860, 31,526, 31,576, 25,808, 11,542, and 21,721 unigenes were annotated to the NR, NT, Swiss-Prot, KEGG, COG, and GO databases, respectively, and total annotation unigenes were 37,764. After 30 days of exposure to 55 μg Cu/l, a total of 292 and 1076 genes were significantly up- and down-regulated, respectively. By KEGG analysis, 660 had a specific pathway annotation. Subsequent bioinformatics analysis revealed that the differentially expressed genes were mainly related to lipid metabolism, immune system, apoptosis, and signal transduction, suggesting that these signaling pathways may be regulated by Cu exposure. The present study provides comprehensive sequence information for subsequent gene expression studies regarding S. hasta, and the transcriptome profiling after Cu exposure is also expected to improve our understanding of the molecular toxicology of Cu.

Transcriptomic response to water accommodated fraction of crude oil exposure in the gill of Japanese flounder, Paralichthys olivaceus

Illumina-based RNA-seq was used to determine the short-term transcriptomic responses of Paralichthys olivaceus gill to an environmentally relevant level of water accommodated fraction (WAF) of crude oil. 213,979 transcripts and 128,482 unigenes were obtained. Differential expression analysis revealed that 1641 and 2142 genes were significantly up- and down-regulated. Enrichment analysis identified a set of GO terms and putative pathways involved in the response of P. olivaceus to WAF exposure. Analysis of the transcripts revealed the effective protective mechanisms of P. olivaceus to reduce the toxic effects of WAF. Moreover, WAF exposure induced the metabolism of energy substrates, and downstream pathway genes were modified to provide protection against toxic damage. Transcripts analysis demonstrated that the genes involved in circadian rhythm signaling were regulated in gills of P. olivaceus exposed to WAF. These results provide insights into the mechanisms of WAF-induced toxicity in fishes and into the WAF-sensitive biomarkers in P. olivaceus.

Actions of Bisphenol A and Bisphenol S on the Reproductive Neuroendocrine System During Early Development in Zebrafish

Bisphenol A (BPA) is a well-known environmental, endocrine-disrupting chemical, and bisphenol S (BPS) has been considered a safer alternative for BPA-free products. The present study aims to evaluate the impact of BPA and BPS on the reproductive neuroendocrine system during zebrafish embryonic and larval development and to explore potential mechanisms of action associated with estrogen receptor (ER), thyroid hormone receptor (THR), and enzyme aromatase (AROM) pathways. Environmentally relevant, low levels of BPA exposure during development led to advanced hatching time, increased numbers of GnRH3 neurons in both terminal nerve and hypothalamus, increased expression of reproduction-related genes (kiss1, kiss1r, gnrh3, lhβ, fshβ, and erα), and a marker for synaptic transmission (sv2). Low levels of BPS exposure led to similar effects: increased numbers of hypothalamic GnRH3 neurons and increased expression of kiss1, gnrh3, and erα. Antagonists of ER, THRs, and AROM blocked many of the effects of BPA and BPS on reproduction-related gene expression, providing evidence that those three pathways mediate the actions of BPA and BPS on the reproductive neuroendocrine system. This study demonstrates that alternatives to BPA used in the manufacture of BPA-free products are not necessarily safer. Furthermore, this is the first study to describe the impact of low-level BPA and BPS exposure on the Kiss/Kiss receptor system during development. It is also the first report of multiple cellular pathways (ERα, THRs, and AROM) mediating the effects of BPA and BPS during embryonic development in any species.

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.

Development of approaches to induce puberty in cultured female sablefish (Anoplopoma fimbria)

Efforts to establish sustainable and efficient aquaculture production of sablefish (Anoplopoma fimbria) have been constrained by delayed puberty in cultured females. This study integrates a series of experiments aimed at gaining an understanding of the reproductive physiology of puberty in female sablefish. We detected transcripts for the dopamine D2 receptor (drd2) in brain, pituitary and ovary of sablefish, and prepubertal females exhibited significantly elevated brain and pituitary drd2 expression relative to wild maturing females. Treatments with sustained-release cholesterol pellets containing testosterone (T) and the dopamine D2 receptor antagonist, metoclopramide (Met), stimulated expression of pituitary luteinizing hormone beta subunit (lhb) and follicle-stimulating hormone beta subunit (fshb), respectively, in prepubertal females, whereas a combination of T and gonadotropin-releasing hormone agonist (GnRHa) had a strong synergistic effect on lhb expression (2000-fold higher than control). Although T induced a significant increase in the maximum ovarian follicle volume, none of the treatments tested stimulated onset of vitellogenesis. Using liquid chromatography/tandem mass spectrometry, we demonstrated that Met stimulated production of T by previtellogenic ovarian follicles in vitro, whereas gonadotropin preparations enhanced 17α-hydroxyprogesterone, androstenedione (A4), T and 17β-estradiol (E2) production. Treatment with T increased production of A4, 11β-hydroxyandrostenedione, 11β-hydroxytestosterone, E2, 11-ketotestosterone, and 5α-dihydrotestosterone (DHT). Interestingly, in the presence of high doses of T the previtellogenic ovary preferentially produced A4 and DHT over any other metabolite. Our data suggest the existence of dopamine inhibition of the reproductive axis in female sablefish. Treatments with Met and T elevated gonadotropin mRNAs in prepubertal females but failed to stimulate the transition into vitellogenic growth, suggesting a possible failure in pituitary gonadotropin protein synthesis/release. Previtellogenic ovarian follicles of sablefish are equipped to synthesize steroids, including those required for vitellogenic growth, and DHT, a steroid hormone whose role in reproduction of fishes remains unknown.

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.

Transcriptome analysis of the Capra hircus ovary

BACKGROUND

Capra hircus is an important economic livestock animal, and therefore, it is necessary to discover transcriptome information about their reproductive performance. In this study, we performed de novo transcriptome sequencing to produce the first transcriptome dataset for the goat ovary using high-throughput sequencing technologies. The result will contribute to research on goat reproductive performance.

METHOD AND RESULTS

RNA-seq analysis generated more than 38.8 million clean paired end (PE) reads, which were assembled into 80,069 unigenes (mean size = 619 bp). Based on sequence similarity searches, 64,824 (60.6%) genes were identified, among which 29,444 and 11,271 unigenes were assigned to Gene Ontology (GO) categories and Clusters of Orthologous Groups (COG), respectively. Searches in the Kyoto Encyclopedia of Genes and Genomes pathway database (KEGG) showed that 27,766 (63.4%) unigenes were mapped to 258 KEGG pathways. Furthermore, we investigated the transcriptome differences of goat ovaries at two different ages using a tag-based digital gene expression system. We obtained a sequencing depth of over 5.6 million and 5.8 million tags for the two ages and identified a large number of genes associated with reproductive hormones, ovulatory cycle and follicle. Moreover, many antisense transcripts and novel transcripts were found; clusters with similar differential expression patterns, enriched GO terms and metabolic pathways were revealed for the first time with regard to the differentially expressed genes.

CONCLUSIONS

The transcriptome provides invaluable new data for a functional genomic resource and future biological research in Capra hircus, and it is essential for the in-depth study of candidate genes in breeding programs.

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.

A novel model for development, organization, and function of gonadotropes in fish pituitary

The gonadotropins follicle-stimulating hormone (FSH) and luteinizing hormone (LH) are key regulators of the reproductive axis in vertebrates. Despite the high popularity of zebrafish as a model organism for studying reproductive functions, to date no transgenic zebrafish with labeled gonadotropes have been introduced. Using gonadotropin regulatory elements from tilapia, we generated two transgenic zebrafish lines with labeled gonadotropes. The tilapia and zebrafish regulatory sequences were highly divergent but several conserved elements allowed the tilapia promoters to correctly drive the transgenes in zebrafish pituitaries. FSH cells reacted to stimulation with gonadotropin releasing hormone by proliferating and showing increased transgene fluorescence, whereas estrogen exposure caused a decrease in cell number and transgene fluorescence. Transgene fluorescence reflected the expression pattern of the endogenous fshb gene. Ontogenetic expression of the transgenes followed typical patterns, with FSH cells appearing early in development, and LH cells appearing later and increasing dramatically in number with the onset of puberty. Our transgenic lines provide a powerful tool for investigating the development, anatomy, and function of the reproductive axis in lower vertebrates.