Microarray analysis in the zebrafish (Danio rerio) liver and telencephalon after exposure to low concentration of 17alpha-ethinylestradiol

The uses of zebrafish (Danio rerio) as an in vivo model for toxicological studies: A review based on bibliometrics

An in vivo model is necessary for toxicology. This review analyzed the uses of zebrafish (Danio rerio) in toxicology based on bibliometrics. Totally 56,816 publications about zebrafish from 2002 to 2023 were found in Web of Science Core Collection, with Toxicology as the top 6 among all disciplines. Accordingly, the bibliometric map reveals that "toxicity" has become a hot keyword. It further reveals that the most common exposure types include acute, chronic, and combined exposure. The toxicological effects include behavioral, intestinal, cardiovascular, hepatic, endocrine toxicity, neurotoxicity, immunotoxicity, genotoxicity, and reproductive and transgenerational toxicity. The mechanisms include oxidative stress, inflammation, autophagy, and dysbiosis of gut microbiota. The toxicants commonly evaluated by using zebrafish model include nanomaterials, arsenic, metals, bisphenol, and dioxin. Overall, zebrafish provide a unique and well-accepted model to investigate the toxicological effects and mechanisms. We also discussed the possible ways to address some of the limitations of zebrafish model, such as the combination of human organoids to avoid species differences.

A complex mixture of polycyclic aromatic compounds causes embryotoxic, behavioral, and molecular effects in zebrafish larvae (Danio rerio), and in vitro bioassays

Polycyclic aromatic compounds (PACs) are prevalent in the environment, typically found in complex mixtures and high concentrations. Our understanding of the effects of PACs, excluding the 16 priority polycyclic aromatic hydrocarbons (16 PAHs), remains limited. Zebrafish embryos and in vitro bioassays were utilized to investigate the embryotoxic, behavioral, and molecular effects of a soil sample from a former gasworks site in Sweden. Additionally, targeted chemical analysis was conducted to analyze 87 PACs in the soil, fish, water, and plate material. CALUX® assays were used to assess the activation of aryl hydrocarbon and estrogen receptors, as well as the inhibition of the androgen receptor. Larval behavior was measured by analyzing activity during light and darkness and in response to mechanical stimulation. Furthermore, qPCR analyses were performed on a subset of 36 genes associated with specific adverse outcomes, and the total lipid content in the larvae was measured. Exposure to the sample resulted in embryotoxic effects (LC50 = 0.480 mg dry matter soil/mL water). The mixture also induced hyperactivity in darkness and hypoactivity in light and in response to the mechanical stimulus. qPCR analysis revealed differential regulation of 15 genes, including downregulation of opn1sw1 (eye pigmentation) and upregulation of fpgs (heart failure). The sample caused significant responses in three bioassays (ERα-, DR-, and PAH-CALUX), and the exposed larvae exhibited elevated lipid levels. Chemical analysis identified benzo[a]pyrene as the predominant compound in the soil and approximately half of the total PAC concentration was attributed to the 16 PAHs. This study highlights the value of combining in vitro and in vivo methods with chemical analysis to assess toxic mechanisms at specific targets and to elucidate the possible interactions between various pathways in an organism. It also enhances our understanding of the risks associated with environmental mixtures of PACs and their distribution during toxicity testing.

Comparative Physiological and Transcriptomic Profiling Offers Insight into the Sexual Dimorphism of Hepatic Metabolism in Size-Dimorphic Spotted Scat (Scatophagus argus)

The spotted scat (Scatophagus argus) is an economically important cultured marine fish that exhibits a typical sexual size dimorphism (SSD). SSD has captivated considerable curiosity for farmed fish production; however, up till now the exact underlying mechanism remains largely unclear. As an important digestive and metabolic organ, the liver plays key roles in the regulation of fish growth. It is necessary to elucidate its significance as a downstream component of the hypothalamic-pituitary-liver axis in the formation of SSD. In this study, the liver physiological differences between the sexes were evaluated in S. argus, and the activity of several digestive and metabolic enzymes were affected by sex. Females had higher amylase, protease, and glucose-6-phosphate dehydrogenase activities, while males exhibited markedly higher hepatic lipase and antioxidant enzymes activities. A comparative transcriptomics was then performed to characterize the responsive genes. Illumina sequencing generated 272.6 million clean reads, which were assembled into 79,115 unigenes. A total of 259 differentially expressed genes were identified and a few growth-controlling genes such as igf1 and igfbp1 exhibited female-biased expression. Further analyses showed that several GO terms and pathways associated with metabolic process, particularly lipid and energy metabolisms, were significantly enriched. The male liver showed a more active mitochondrial energy metabolism, implicating an increased energy expenditure associated with reproduction. Collectively, the female-biased growth dimorphism of S. argus may be partially attributed to sexually dimorphic metabolism in the liver. These findings would facilitate further understanding of the nature of SSD in teleost fish.

Transcriptome responses in polar cod (Boreogadus saida) liver slice culture exposed to benzo[a]pyrene and ethynylestradiol: insights into anti-estrogenic effects

Polar cod (Boreogadus saida) is a key species in the arctic marine ecosystem vulnerable to effects of pollution, particularly from petroleum related activities. To facilitate studying the effects of those pollutants, we adapted a precision-cut liver slice culture protocol for this species. Using this system on board a research vessel, we studied gene expression in liver slice after exposure to the polycyclic aromatic hydrocarbon (PAH) benzo[a]pyrene (BaP), ethynylestradiol (EE2), and their mixtures, to map their molecular targets and examine possible anti-estrogenic effects of BaP. The exposure experiments were performed with BaP alone (0.1, 1, and 10 μM) or in combination with low concentrations of EE2 (5 nM) to mimic physiological estradiol levels in early vitellogenic female fish. Transcriptome analysis (RNA-seq) was performed after 72 h exposure in culture to map the genes and cellular pathways affected. The results provide a view of global transcriptome responses to BaP and EE2, which resulted in enrichment of many pathways such as the aryl hydrocarbon (Ahr) and estrogen receptor pathways. In the mixture exposure, BaP resulted in anti-estrogenic effects, shown by attenuation of EE2 activated transcription of many estrogen target genes. The results from this ex vivo experiment suggest that pollutants that activate the Ahr pathway such as the PAH compound BaP can result in anti-estrogenic effects that may lead to endocrine disruption in polar cod.

Investigation into Adaptation in Genes Associated with Response to Estrogenic Pollution in Populations of Roach (Rutilus rutilus) Living in English Rivers

Exposure of male fish to estrogenic substances from wastewater treatment works (WwTWs) results in feminization and reduced reproductive fitness. Nevertheless, self-sustaining populations of roach (Rutilus rutilus) inhabit river stretches polluted with estrogenic WwTW effluents. In this study, we examine whether such roach populations have evolved adaptations to tolerate estrogenic pollution by comparing frequency differences in single-nucleotide polymorphisms (SNPs) between populations sampled from rivers receiving either high- or low-level WwTW discharges. SNPs within 36 "candidate" genes, selected for their involvement in estrogenic responses, and 120 SNPs in reference genes were genotyped in 465 roaches. There was no evidence for selection in highly estrogen-dependent candidate genes, including those for the estrogen receptors, aromatases, and vitellogenins. The androgen receptor (ar) and cytochrome P450 1A genes were associated with large shifts in allele frequencies between catchments and in individual populations, but there is no clear link to estrogen pollution. Selection at ar in the effluent-dominated River Lee may have resulted from historical contamination with endocrine-disrupting pesticides. Critically, although our results suggest population-specific selection including at genes related to endocrine disruption, there was no strong evidence that the selection resulted from exposure to estrogen pollution.

Shedding new light on early sex determination in zebrafish

In contrast to established zebrafish gene annotations, the question of sex determination has still not been conclusively clarified for developing zebrafish, Danio rerio, larvae, 28 dpf or earlier. Recent studies indicate polygenic sex determination (PSD), with the genes being distributed throughout the genome. Early genetic markers of sex in zebrafish help unravel co-founding sex-related differences to apply to human health and environmental toxicity studies. A qPCR-based method was developed for six genes: cytochrome P450, family 17, subfamily A, polypeptide 1 (cyp17a1); cytochrome P450, family 19, subfamily A, polypeptide 1a (cyp19a1a); cytochrome P450, family 19, subfamily A, polypeptides 1b (cyp19a1b); vitellogenin 1 (vtg1); nuclear receptor subfamily 0, group B, member 1 (nr0b1), sry (sex-determining region Y)-box 9b (sox9b) and actin, beta 1 (actb1), the reference gene. Sry-box 9a (Sox9a), insulin-like growth factor 3 (igf3) and double sex and mab-3 related transcription factor 1 (dmrt1), which are also known to be associated with sex determination, were used in gene expression tests. Additionally, Next-Generation-Sequencing (NGS) sequenced the genome of two adult female and male and two juveniles. PCR analysis of adult zebrafish revealed sex-specific expression of cyp17a1, cyp19a1a, vtg1, igf3 and dmrt1, the first four strongly expressed in female zebrafish and the last one highly expressed in male conspecifics. From NGS, nine female and four male-fated genes were selected as novel for assessing zebrafish sex, 28 dpf. Differences in transcriptomes allowed allocation of sex-specific genes also expressed in juvenile zebrafish.

The Landscape of Immune Cells Infiltrating in Prostate Cancer

Background

This study was to explore the infiltration pattern of immune cells in the prostate cancer (PCa) microenvironment and evaluate the possibility of specific infiltrating immune cells as potential prognostic biomarkers in PCa.

Methods

Infiltrating percentage of 22 immune cells were extracted from 27 normalized datasets by CIBERSORT algorithm. Samples with CIBERSORT p-value < 0.05 were subsequently merged and divided into normal or tumor groups. The differences of 22 immune cells between normal and tumor tissues were analyzed along with potential infiltrating correlations among 22 immune cells and Gleason grades. SNV data from TCGA was used to calculate the TMB score. A univariate and multivariate regression were used to evaluate the prognostic effects of immune cells in PCa.

Results

Ten immune cells with significant differences were identified, including seven increased and three decreased infiltrating immune cells from 190 normal prostate tissues and 537 PCa tissues. Among them, the percentage of infiltration of resting NK cells increased the most, whereas the percentage of infiltration of resting mast cells decreased the most. In normal tissues, CD8+ T cells had the strongest infiltrating correlation with monocytes, while activated NK cells and naive B cells were the highest in PCa tissues. Moreover, the infiltration of five immune cells was significantly associated with TMB score and mutations of immune gene change the infiltration of immune cells. The Area Under Curve (AUC) of the multivariate regression model for the five- and 10-year survival prediction of PCa reached 0.796 and 0.862. The validation cohort proved that the model was reproducible.

Conclusions

This study demonstrated that different infiltrating immune cells in prostate cancer, especially higher infiltrating M1 macrophages and neutrophils in PCa tissue, are associated with patients’ prognosis, suggesting that these two immune cells might be potential targets for PCa diagnosis and prognosis of treatment.

Risk assessment of the lipid metabolism-disrupting effects of nitro-PAHs

Nitro-polycyclic aromatic hydrocarbons (NPAHs) are of increasing global concern due to their ubiquitous occurrence and long-range transport in the environment. However, their potential metabolism-disrupting effects, especially nuclear receptor-related lipid disorders, are still poorly understood. Targeting estrogen receptor α (ERα), this study for the first time evaluated the lipid metabolic effects of NPAHs using in vitro and in vivo models. The results indicated that four of the five NPAHs tested exhibited significant ERα agonistic activities, and induced increased secretion of 17β-estradiol (E₂) in HepG2 cells. Furthermore, lipidomic analysis showed that exposure to the candidate NPAH (3-nitrofluoranthene, 3-NFA) led to elevated hepatic levels of triacylglycerols (TAGs) and cholesteryl esters (CEs). Importantly, the lipid overload induced by 3-NFA was verified in the livers of zebrafish larvae using Oil Red O staining. Additionally, significant increases in E₂ production and the expression levels of associated genes (17βHSD and C/EBP-α) further supported the involvement of the ERα signaling pathway in the lipid metabolic perturbation induced by 3-NFA. These results provide novel insight into the lipid metabolism-disrupting effects induced by NPAHs and may offer a better understanding of the environmental risks of NPAHs.

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.

Does hepatotoxicity interfere with endocrine activity in zebrafish (Danio rerio)?

Vitellogenin (VTG), a well-established biomarker for the diagnosis of endocrine activity in fish, is used in multiple OECD test guidelines (TG) to identify activities of chemicals on hormonal pathways. However, the synthesis of VTG may not only be modified by typical endocrine-related pathways, but also through non-endocrine-mediated processes. In particular, hepatotoxicity, i.e. toxicant-induced impairment of liver structure and function, might influence VTG as a biomarker, since VTG is synthesized in hepatocytes. An intimate understanding of the interplay between endocrine-related and non-endocrine-related pathways influencing VTG production is crucial for the avoidance of erroneous diagnoses in hazard assessment for regulatory purposes of chemical compounds. In order to investigate whether hepatotoxicity may interfere with hepatic VTG synthesis, adult zebrafish (Danio rerio) were exposed to three well-known hepatotoxicants, acetaminophen, isoniazid and acetylsalicylic acid, according to OECD TG 230. Various hepatotoxicity- and endocrine system-related endpoints were recorded: mRNA expression of selected endocrine- and hepatotoxicity-related marker genes in the liver; VTG levels in head/tail homogenates; and liver histopathology. All three test compounds induced significant, but mild single cell necrosis of hepatocytes and transcriptional changes of hepatotoxicity-related marker genes, thus confirming hepatotoxic effects. A positive correlation between hepatotoxicity and reduced hepatic VTG synthesis was not observed, with the single exception of a weak increase in female zebrafish exposed to APAP. This suggests that - in studies conducted according to OECD TG 229 or 230 - it is unlikely that hepatotoxic chemicals will interfere with the hepatic capacity for VTG synthesis.

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.

Telencephalon transcriptome analysis of chronically stressed adult zebrafish

Chronic stress leads to disruptions in learning and memory processes. The effects of chronic stress experience on the adult zebrafish brain, particularly the memory associated telencephalon brain region, is unclear. The goal of this study was to identify gene expression changes in the adult zebrafish brain triggered by chronic unpredictable stress. Transcriptome analysis of the telencephalon revealed 155 differentially expressed genes. Of these genes, some are critical genes involved in learning and memory, such as cdk5 and chrna7, indicating effects of chronic unpredictable stress on zebrafish memory. Interestingly, several genes were annotated in the Orange domain, which is an amino acid sequence present in eukaryotic DNA-binding transcription repressors. Furthermore, we identified hsd11b2, a cortisol inactivating gene, as chronic stress-responsive in the whole zebrafish brain. Collectively, these findings suggest that memory associated gene expression changes in adult zebrafish telencephalon are affected by chronic stress experience.

Silencing of PPARαBb mRNA in brown trout primary hepatocytes: effects on molecular and morphological targets under the influence of an estrogen and a PPARα agonist

The crosstalk between peroxisome proliferator-activated receptor α (PPARα) and estrogenic pathways are shared from fish to humans. Salmonid fish had an additional genome duplication, and two PPARα isoforms (PPARαBa and PPARαBb) were previously identified. Since a negative regulation between estrogen signaling and PPARα was described, a post-transcriptional gene silencing for PPARαBb was designed in primary brown trout hepatocytes. The aims of the study were to: (i) decipher the effects of PPARαBb knock-down on peroxisome morphology and on mRNA expression of potential target genes, and (ii) to assess the cross-interferences caused by an estrogenic compound (17α-ethinylestradiol - EE2) and a PPARα agonist (Wy-14,643 - Wy) using the established knock-down model. A knock-down efficiency of 70% was achieved for PPARαBb and its silencing significantly reduced the volume density of peroxisomes, but did not alter mRNA levels of the studied genes. Exposure to Wy did not change peroxisome morphology or mRNA expression, but under silencing conditions Wy rescued the volume density of peroxisomes to control levels, and increased acyl-coenzyme A oxidase 1-3l (Acox1-3l) mRNA. Exposure to EE2 caused a reduction of peroxisome volume density, but under silencing conditions this effect was abolished and ApoA1 mRNA level was diminished. The morphological alterations of peroxisomes by WY and EE2 demonstrated that obtained results are PPARαBb dependent, and suggest the regulation of unknown downstream targets of PPARαBb. In summary, PPARαBb is involved in the control of peroxisome size and/or number, which opens future opportunities to explore its regulation and molecular targets.

Sex-steroids and hypolipidemic chemicals impacts on brown trout lipid and peroxisome signaling - Molecular, biochemical and morphological insights

Lipid metabolism involves complex pathways, which are regulated in a similar way across vertebrates. Hormonal and hypolipidemic deregulations cause lipid imbalance from fish to humans, but the underlying mechanisms are far from understood. This study explores the potential of using juvenile brown trout to evaluate the in vivo interferences caused by estrogenic (17α-ethinylestradiol - EE2), androgenic (testosterone - T), and hypolipidemic (clofibrate - CLF) compounds in lipidic and/or peroxisomal pathways. Studied endpoints were from blood/plasma biochemistry, plasma fatty acid profile, ultrastructure of hepatocytes and abundance of their peroxisomes to mRNA expression in the liver. Both T and CLF caused minimal effects when compared to EE2. Estrogenized fish had significantly higher hepatosomatic indexes, increased triglycerides and very-low density lipoproteins (VLDL) in plasma, compared with solvent control. Morphologically, EE2 fish showed increased lipid droplets in hepatocytes, and EE2 and T reduced volume density of peroxisomes in relation to the hepatic parenchyma. Polyunsaturated fatty acids (PUFA) in plasma, namely n-3 PUFA, increased with EE2. EE2 animals had increased mRNA levels of vitellogenin A (VtgA), estrogen receptor alpha (ERα), peroxisome proliferator-activated receptor alpha (PPARα), PPARαBa and acyl-CoA long chain synthetase 1 (Acsl1), while ERβ-1, acyl-CoA oxidase 1-3I (Acox1-3I), Acox3, PPARγ, catalase (Cat), urate oxidase (Uox), fatty acid binding protein 1 (Fabp1) and apolipoprotein AI (ApoAI) were down-regulated. In summary, in vivo EE2 exposure altered lipid metabolism and peroxisome dynamics in brown trout, namely by changing the mRNA levels of several genes. Our model can be used to study possible organism-level impacts, viz. in gonadogenesis.

Estrogenic effects following larval exposure to the putative anti-estrogen, fulvestrant, in the fathead minnow (Pimephales promelas)

The objective of the present study was to investigate the consequences of early-life exposure to fulvestrant on estrogenic gene expression in fathead minnow larvae. To address this objective, fathead minnow larvae were exposed to fulvestrant (ICI 182,780) during the window of sexual differentiation between 0 to 30 days post-hatch (dph). The four treatment groups in this study included: filtered water controls (never exposed), solvent controls (ethanol 0.01%), and nominally low (0.10μg/L) and high (10.0μg/L) doses of fulvestrant. Following 30 d exposure to their respective treatment, larvae were transferred to filtered water aquaria and assessed for alterations in endocrine-responsive gene expression (i.e., RT-qPCR), body size and survival. The remaining fish depurated in filtered water until reaching sexual maturity (180dph) for assessment of persistent effects on sex characteristics, reproductive performance and sex ratio. Following the 30-d early life exposure, larvae showed upregulations of the endocrine-responsive genes ar, erβ and vtg in response to both low and high doses of fulvestrant, but showed no differences in survival or body mass. Upon reaching sexual maturity under depuration conditions, male minnows previously exposed to fulvestrant as larvae showed reductions in gonad mass along with the feminization of secondary sex characteristics with no observed effects in females. Exposure to fulvestrant had no effects on gonadal histology, reproductive performance or final sex ratio as adults. Results from this study demonstrate that aqueous exposure to fulvestrant is estrogenic in fathead minnow larvae and is capable of feminizing male fish as adults following early life exposure.

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.

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.

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.

Nuclear receptor research in zebrafish

Nuclear receptors (NRs) form a superfamily of transcription factors that can be activated by ligands and are involved in a wide range of physiological processes. NRs are well conserved between vertebrate species. The zebrafish, an increasingly popular animal model system, contains a total of 73 NR genes, and orthologues of almost all human NRs are present. In this review article, an overview is presented of NR research in which the zebrafish has been used as a model. Research is described on the three most studied zebrafish NRs: the estrogen receptors (ERs), retinoic acid receptors (RARs) and peroxisome proliferator-activated receptors (PPARs). The studies on these receptors illustrate the versatility of the zebrafish as a model for ecotoxicological, developmental and biomedical research. Although the use of the zebrafish in NR research is still relatively limited, it is expected that in the next decade the full potential of this animal model will be exploited.

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.

Neuroendocrine disruption of organizational and activational hormone programming in poikilothermic vertebrates

In vertebrates, sexual differentiation of the reproductive system and brain is tightly orchestrated by organizational and activational effects of endogenous hormones. In mammals and birds, the organizational period is typified by a surge of sex hormones during differentiation of specific neural circuits; whereas activational effects are dependent upon later increases in these same hormones at sexual maturation. Depending on the reproductive organ or brain region, initial programming events may be modulated by androgens or require conversion of androgens to estrogens. The prevailing notion based upon findings in mammalian models is that male brain is sculpted to undergo masculinization and defeminization. In absence of these responses, the female brain develops. While timing of organizational and activational events vary across taxa, there are shared features. Further, exposure of different animal models to environmental chemicals such as xenoestrogens such as bisphenol A-BPA and ethinylestradiol-EE2, gestagens, and thyroid hormone disruptors, broadly classified as neuroendocrine disrupting chemicals (NED), during these critical periods may result in similar alterations in brain structure, function, and consequently, behaviors. Organizational effects of neuroendocrine systems in mammals and birds appear to be permanent, whereas teleost fish neuroendocrine systems exhibit plasticity. While there are fewer NED studies in amphibians and reptiles, data suggest that NED disrupt normal organizational-activational effects of endogenous hormones, although it remains to be determined if these disturbances are reversible. The aim of this review is to examine how various environmental chemicals may interrupt normal organizational and activational events in poikilothermic vertebrates. By altering such processes, these chemicals may affect reproductive health of an animal and result in compromised populations and ecosystem-level effects.

Gene Expression Profiling in Fish Toxicology: A Review

In this review, we present an overview of transcriptomic responses to chemical exposures in a variety of fish species. We have discussed the use of several molecular approaches such as northern blotting, differential display reverse transcription-polymerase chain reaction (DDRT-PCR), suppression subtractive hybridization (SSH), real time quantitative PCR (RT-qPCR), microarrays, and next-generation sequencing (NGS) for measuring gene expression. These techniques have been mainly used to measure the toxic effects of single compounds or simple mixtures in laboratory conditions. In addition, only few studies have been conducted to examine the biological significance of differentially expressed gene sets following chemical exposure. Therefore, future studies should focus more under field conditions using a multidisciplinary approach (genomics, proteomics and metabolomics) to understand the synergetic effects of multiple environmental stressors and to determine the functional significance of differentially expressed genes. Nevertheless, recent developments in NGS technologies and decreasing costs of sequencing holds the promise to uncover the complexity of anthropogenic impacts and biological effects in wild fish populations.

In silico prediction and in vivo analysis of antiestrogenic potential of 2-isopropylthioxanthone (2-ITX) in juvenile goldfish (Carassius auratus)

Previous studies have shown both anti-estrogenic and anti-androgenic activities of 2-isopropylthioxanthone (2-ITX), a well known food contaminant, in in vitro assays. However, no data are available on the anti-estrogenic potentials and risks of 2-ITX in aquatic organisms. This work evaluated the potential endocrine disrupting effects of 2-ITX at the level of estrogen receptor (ER) signaling cascade using juvenile goldfish (Carassius auratus) as model. Firstly, we investigated the ligand binding efficiency of 2-ITX to the ligand binding domains (LBD) of goldfish ER subtypes using a molecular docking approach. Secondly, we assessed the effects of 2-ITX on E2-induced hepatic expression of ERα1, ERβ1, ERβ2, and vitellogenin (VTG) in vivo. Crosstalk between ER-VTG and aryl hydrocarbon receptor 2 (AhR2)-cytochrome P4501A (CYP1A) was also investigated. Fish were injected with increasing doses of 2-ITX ranging from 2 to 10µg/g BW, and results were compared to the effect of tamoxifen, a well-known ER modulator. We observed that compared to ERβ, the interaction potentials of 2-ITX to goldfish ERα1 LBD was more stable in the inactive receptor conformation. The in silico docking simulation analysis also revealed that 2-ITX acted as agonist for the goldfish AhR2 LBDs suggesting the ability of this compound to activate the cross-talk between the ERα- and AhR-signaling pathways. In vivo experiments confirm in silico simulation predictions demonstrating that 2-ITX reduced the estrogenicity of E2 at both transcriptional and post-transcriptional levels, indicating a clear anti-estrogenic effect. Co-exposure of E2 and 2-ITX also resulted in a significant decrease of CYP1A gene expression with respect to 2-ITX alone. Results from these studies collectively revealed that the antiestrogenic property of 2-ITX can be ascribed to a combination of effects on multiple signaling pathways suggesting the potential for this environmental contaminant to affect the hormonal control of reproductive processes in fish.

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.

β-Estradiol and ethinyl-estradiol contamination in the rivers of the Carpathian Basin

17β-Estradiol (E2) and 17α-ethinyl estradiol (EE2), which are environmental estrogens, have been determined with LC-MS in freshwater. Their sensitive analysis needs derivatization and therefore is very hard to achieve in multiresidue screening. We analyzed samples from all the large and some small rivers (River Danube, Drava, Mur, Sava, Tisza, and Zala) of the Carpathian Basin and from Lake Balaton. Freshwater was extracted on solid phase and derivatized using dansyl chloride. Separation was performed on a Kinetex XB-C18 column. Detection was achieved with a benchtop orbitrap mass spectrometer using targeted MS analysis for quantification. Limits of quantification were 0.05 ng/L (MS1) and 0.1 ng/L (MS/MS) for E2, and 0.001 ng/L (MS1) and 0.2 ng/L (MS/MS) for EE2. River samples contained n.d.-5.2 ng/L E2 and n.d.-0.68 ng/L EE2. Average levels of E2 and EE2 were 0.61 and 0.084 ng/L, respectively, in rivers, water courses, and Lake Balaton together, but not counting city canal water. EE2 was less abundant, but it was still present in almost all of the samples. In beach water samples from Lake Balaton, we measured 0.076-0.233 E2 and n.d.-0.133 EE2. A relative high amount of EE2 was found in river Zala (0.68 ng/L) and in Hévíz-Páhoki canal (0.52 ng/L), which are both in the catchment area of Lake Balaton (Hungary).

Zebrafish neurobehavioral phenomics for aquatic neuropharmacology and toxicology research

Zebrafish (Danio rerio) are rapidly emerging as an important model organism for aquatic neuropharmacology and toxicology research. The behavioral/phenotypic complexity of zebrafish allows for thorough dissection of complex human brain disorders and drug-evoked pathological states. As numerous zebrafish models become available with a wide spectrum of behavioral, genetic, and environmental methods to test novel drugs, here we discuss recent zebrafish phenomics methods to facilitate drug discovery, particularly in the field of biological psychiatry. Additionally, behavioral, neurological, and endocrine endpoints are becoming increasingly well-characterized in zebrafish, making them an inexpensive, robust and effective model for toxicology research and pharmacological screening. We also discuss zebrafish behavioral phenotypes, experimental considerations, pharmacological candidates and relevance of zebrafish neurophenomics to other 'omics' (e.g., genomic, proteomic) approaches. Finally, we critically evaluate the limitations of utilizing this model organism, and outline future strategies of research in the field of zebrafish phenomics.

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.

Dose-dependent effects of metals on gene expression in the sydney rock oyster, Saccostrea glomerata

In the current study, we tested the effects of common environmental contaminants (the metals zinc and lead) on gene expression in Sydney rock oysters (Saccrostrea glomerata). Oysters were exposed to a range of metal concentrations under controlled laboratory conditions. The expression of 14 putative stress response genes was then measured using quantitative, real-time (q) PCR. The expression of all 14 genes was significantly affected (p < 0.05 vs. nonexposed controls) by at least one of the metals, and by at least one dose of metal. For 5 of the 14 target genes (actin, calmodulin, superoxide dismutase, topoisomerase I, and tubulin) the alteration of expression relative to controls was highest at intermediate (rather than high) doses of metals. Such responses may reflect adaptive (acclimation) reactions in gene expression at low to intermediate doses of contaminants, followed by a decline in expression resulting from exposure at higher doses. The data are discussed in terms of the intracellular pathways affected by metal contamination, and the relevance of such gene expression data to environmental biomonitoring.

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.

Evaluation of the sensitivity, responsiveness and reproducibility of primary rainbow trout hepatocyte vitellogenin expression as a screening assay for estrogen mimics

Vitellogenin (Vtg) induction in primary fish hepatocytes has been proposed as an in vitro screening assay for ER agonists and antagonists, but has not yet been used extensively as a high-throughput screening tool due to poor reproducibility, sensitivity and overall feasibility. The present work has evaluated the role of seasonal variation, normalization, optimal culture and assay conditions on the sensitivity, responsiveness and reproducibility of in vitro vtg gene mRNA and protein expression in rainbow trout (Oncorhynchus mykiss) primary hepatocytes using the xenoestrogen 17α-ethynylestradiol (EE2) as a test chemical. The results show that primary hepatocytes display a relatively high individual and seasonal variation in both Vtg mRNA and protein induction potential, although less variance was observed in assay sensitivity. Data normalization of assay response to maximum (3 nM EE2) and minimum (DMSO) Vtg production dramatically reduced this variance and led to improved assay reproducibility. A time-dependent response was observed both for mRNA and protein expression, reaching maximum Vtg induction after 96 h of exposure, although reproducible concentration response curves for both Vtg mRNA and protein could be obtained already after 48 h. A need for chemical re-exposure of the hepatocytes was identified to be important for sustaining exposure concentrations in extended studies (>48 h), whereas different plate formats (96, 24 or 6 wells) did not affect the bioassay performance. In conclusion, standardization of hepatocyte bioassay and test conditions as well as data-normalization procedures are proposed to be instrumental for more consistent and comparable results in future use of this in vitro assay.

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.

Nonmonotonic response of vitellogenin and estrogen receptor α gene expression after octylphenol exposure of Cichlasoma dimerus (Perciformes, Cichlidae)

In oviparous vertebrates, vitellogenin (VTG) is mainly produced by the liver in response to estrogen (E2) and its synthesis is traditionally coupled to estrogen receptor alpha induction. Even though VTG is a female-specific protein, chemicals that mimic natural estrogens, known as xenoestrogens, can activate its expression in males causing endocrine disruption to wildlife and humans. Alkylphenols such as nonylphenol (NP) and octylphenol (OP) are industrial additives used in the manufacture of a wide variety of plastics and detergents, and can disrupt endocrine functions in exposed animals. For more than a decade, the freshwater cichlid fish Cichlasoma dimerus has been used for ecotoxicological studies in our laboratory. We recently found an up-regulation of VTG gene expression in livers of male fish exposed to OP, from a silent state to values similar to those of E2-induced fish. To better understand the underlying mechanisms behind the action of xenoestrogens, the aim of this study was to analyze the dose-response relationship of C. dimerus VTG and estrogen receptors (ERs) gene expression after waterborne exposure to 0.15, 1.5, 15, and 150μg/L OP for up to 1 month (0, 1, 3, 7, 14, 21, and 28 days). At the end of the experiment, histological features of exposed fish included active hepatocytes with basophilic cytoplasm and high eosinophilic content in their vascular system due to augmented expression of VTG. In testis, high preponderance of sperm was found in fish exposed to 150μg/L OP. A classic dose-response down-regulation of the expression of Na(+)/K(+)-ATPase, a "non-gender specific gene" used for comparison, was found with increasing OP concentrations. No VTG and very low levels of ERα were detected in control male livers, but an up-regulation of both genes was found in males exposed to 0.15 or 150μg/L OP. Moreover, VTG transcripts were significant as early as day 3 or day 1 of exposure to these OP concentrations, respectively. Nearly no response was detected in 1.5 and 15μg/L OP exposed-fish. Data was curve-fitted evidencing a nonmonotonic dose-response curve. Interestingly, ERβ2 mRNA expression was augmented above baseline levels only when males were exposed to the lowest OP concentration. We speculate that genomic control of vitellogenesis is under control of multiple steroid receptors with different affinities for ligands. ERβ isoform, only up-regulated with very low concentrations of ligand, would act as a sensors of OP (or E2) to induce ERα and VTG. With high OP concentrations, the expression of ERα isoform is promptly augmented, with the concomitant VTG transactivation.

Short-term treatment of adult male zebrafish (Danio Rerio) with 17α-ethinyl estradiol affects the transcription of genes involved in development and male sex differentiation

The synthetic estrogen 17α-ethinyl estradiol (EE2) disturbs reproduction and causes gonadal malformation in fish. Effects on the transcription of genes involved in gonad development and function that could serve as sensitive biomarkers of reproductive effects in the field is, however, not well known. We have studied mRNA expression in testes and liver of adult zebrafish (Danio rerio) males treated with 0, 5 or 25 ng/L EE2for 14 days. qPCR analysis showed that the mRNA expression of four genes linked to zebrafish male sex determination and differentiation, Anti-Mullerian Hormone, Double sex and mab-related protein, Sry-related HMG box-9a and Nuclear receptor subfamily 5 group number 1b were significantly decreased by 25 ng/L, but not 5 ng/L EE2 compared with the levels in untreated fish. The decreased transcription was correlated with a previously shown spawning failure in these males (Reyhanian et al., 2011. Aquat Toxicol 105, 41-48), suggesting that decreased mRNA expression of genes regulating male sexual function could be involved in the functional sterility. The mRNA level of Cytochrome P-45019a, involved in female reproductive development, was unaffected by hormone treatment. The transcription of the female-specific Vitellogenin was significantly induced in testes. While testicular Androgen Receptor and the Estrogen Receptor-alpha mRNA levels were unchanged, Estrogen receptor-beta was significantly decreased by 25 ng/L EE2. Hepatic Estrogen Receptor-alpha mRNA was significantly increased by both exposure concentrations, while Estrogen Receptor-beta transcription was unaltered. The decreased transcription of male-predominant genes supports a demasculinization of testes by EE2 and might reflect reproductive disturbances in the environment.

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.

Endocrinology: advances through omics and related technologies

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

The antidepressant drug carbamazepine induces differential transcriptome expression in the brain of Atlantic salmon, Salmo salar

Concerns are being expressed recently over possible environmental effects of human pharmaceuticals. Although the likelihood of acute toxicity is low, the continuous discharge of pharmaceuticals into the aquatic environment means that sublethal effects on non-target organisms need to be seriously considered. One-year-old Atlantic salmon parr were exposed to 7.85±0.13μgL(-1) of the antidepressant drug Carbamazepine (CBZ) for five days to investigate changes of mRNA expression in the brain by means of a custom 17k Atlantic salmon cDNA microarray. The selected concentration is similar to upper levels that can be found in hospital and sewage treatment plant effluents. After treatment, 373 features were differently expressed with 26 showing up- or down-regulation of ≥2-fold (p≤0.05). Among the mRNAs showing the highest change were the pituitary hormones encoding features somatolactin, prolactin and somatotropin, or growth hormone. Functional enrichment and network analyses of up- and down-regulated genes showed that CBZ induced a highly different gene expression profile in comparison to untreated organisms. CBZ induced expression of essential genes of the focal adhesion and extracellular matrix - receptor interaction pathways most likely through integrin alpha-6 (itga6) activation. Negative regulation of apoptotic process, extracellular matrix organization and heme biosynthesis were the most enriched biological process related GO-terms, with the simultaneous enrichment of collagen and extracellular region related cellular component GO-terms, and extracellular matrix structural constituent, hormone activity and chromatin binding molecular function related GO-terms. These results show that relatively low doses of CBZ may affect brain physiology in exposed salmon parr, targeting similar processes as in human, indicating a high degree of conservation of targets of CBZ action. However, and since the mRNAs showing most changes in expression are critical for adaptation to different stressors and life history transitions in Atlantic salmon, more research should be undertaken to assess CBZ effects to avoid impairment of normal development and maintenance of natural populations.

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 .

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.

Disruption of zebrafish (Danio rerio) sexual development after full life-cycle exposure to environmental levels of triadimefon

In this study, zebrafish was exposed to environmental levels of triadimefon (0.125, 0.25, 0.5 μg/mL) from 24 h post fertilization to 120 days post fertilization. Several endpoints that related to reproductive function were evaluated. It was found that the body length, body weight and vitellogenin transcription were significantly reduced for fish exposed to triadimefon. Histological examination showed that the sex ratio of fish skewed to male and female exposed to 0.5 μg/mL triadimefon had immature ovary. The breeding success, as determined from data on egg production and spawning, was reduced in fish exposed to 0.25 μg/mL triadimefon. In the offspring, the reduced egg fertility, hatching rate and survival were observed in eggs exposed to 0.5 μg/mL triadimefon. These findings indicated that triadimefon had the potential to adversely affect the sexual development and breeding success through the multiple endocrine actions.

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.

Morpholino-mediated knockdown of ERα, ERβa, and ERβb mRNAs in zebrafish (Danio rerio) embryos reveals differential regulation of estrogen-inducible genes

Genetically distinct estrogen receptor (ER) subtypes (ERα and ERβ) play a major role in mediating estrogen actions in vertebrates, but their unique and overlapping functions are not entirely clear. Although mammals have 1 gene of each subtype (ESR1 and ESR2), teleost fish have a single esr1 (ERα) and 2 esr2 (ERβa and ERβb) genes. To determine the in vivo role of different ER isoforms in regulating estrogen-inducible transcription targets, zebrafish (Danio rerio) embryos were microinjected with esr-specific morpholino (MO) oligonucleotides to disrupt splicing of the exon III/intron III junction in the DNA-binding domain. Each MO knocked down its respective normal transcript and increased production of variants with a retained intron III (esr1 MO) or a deleted or mis-spliced exon III (esr2a and esr2b MOs). Both esr1 and esr2b MOs blocked estradiol induction of vitellogenin and ERα mRNAs, predominant hepatic genes, but esr2b was the only MO that blocked induction of cytochrome P450 aromatase B mRNA, a predominant brain gene. Knockdown of ERβa with the esr2a MO had no effect on estrogen induction of the 3 mRNAs but, when coinjected with esr1 MO, attenuated the effect of ERα knockdown. Results indicate that ERα and ERβb, acting separately or cooperatively on specific gene targets, are positive transcriptional regulators of estrogen action, but the role of ERβa, if any, is unclear. We conclude that MO technology in zebrafish embryos is an advantageous approach for investigating the interplay of ER subtypes in a true physiological context.

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.

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.

Estrogen-like effects in male goldfish co-exposed to fluoxetine and 17 alpha-ethinylestradiol

The antidepressant fluoxetine (FLX) and the synthetic estrogen, 17 alpha-ethinylestradiol (EE2), are present in municipal sewage discharges. To better understand possible interactions between them, male goldfish were exposed to an ethanol control or to nominal concentrations of FLX (0.54 μg/L) and EE2 (5 ng/L) alone and in combination for 14 days. Real-time reverse-transcription polymerase chain reaction was used to assess effects on hepatic gene expression and liquid chromatography tandem mass spectrometry to analyze the plasma proteome. The results showed an increase in estrogen receptor alpha (esr1) and vitellogenin (vtg) gene expression by 1.9-2.4-fold in the FLX and EE2 groups, but this did not reach statistical significance. In contrast, co-exposure up regulated esr1 and vtg gene expression by 5.5- and 5.3-fold, respectively. Fluoxetine and EE2 alone did not affect estrogen receptor beta (esr2), but the co-exposure down regulated esr2 expression by 50%. There was a significant increase in the number of plasma proteins that were related to endocrine system disorders in the FLX and FLX plus EE2 groups. The level of VTG protein was increased in the plasma from goldfish exposed to EE2, FLX, and FLX plus EE2. Our study demonstrates that low concentrations of FLX and EE2 in a simple mixture produce strong estrogen-like effects in the male goldfish.

A novel method for cross-species gene expression analysis

Background

Analysis of gene expression from different species is a powerful way to identify evolutionarily conserved transcriptional responses. However, due to evolutionary events such as gene duplication, there is no one-to-one correspondence between genes from different species which makes comparison of their expression profiles complex.

Results

In this paper we describe a new method for cross-species meta-analysis of gene expression. The method takes the homology structure between compared species into account and can therefore compare expression data from genes with any number of orthologs and paralogs. A simulation study shows that the proposed method results in a substantial increase in statistical power compared to previously suggested procedures. As a proof of concept, we analyzed microarray data from heat stress experiments performed in eight species and identified several well-known evolutionarily conserved transcriptional responses. The method was also applied to gene expression profiles from five studies of estrogen exposed fish and both known and potentially novel responses were identified.

Conclusions

The method described in this paper will further increase the potential and reliability of meta-analysis of gene expression profiles from evolutionarily distant species. The method has been implemented in R and is freely available at http://bioinformatics.math.chalmers.se/Xspecies/.