Testicular transcriptome alterations in zebrafish (Danio rerio) exposure to 17β-estradiol

Full sexual maturity-cycle exposure to environmentally relevant concentrations of 17β-estradiol decreases reproductive capacity of zebrafish

17β-Estradiol (E2) has been widely detected in natural water and treatment with E2 induces potential endocrine disrupting effects in fish. However, effects on fish fecundity and steroid system after treatment with environmentally relevant concentrations of E2 for the full sexual maturation cycle remain unclear. In this study, zebrafish were treatment with 0, 10 or 100 ng/L E2 from embryo to adult stage, and effects on gonadal development and differentiation, steroid hormone levels, transcription of genes associated with the hypothalamic-pituitary-gonadal-liver (HPG) axis in adults and fertilization rate of offspring were assessed. The results showed that treatment with E2 lead to increased number of feminization in zebrafish. In females, E2 decreased cumulative amount of spawning and inhibited the maturation of oocyte. In males, E2 inhibited the maturation and motility of sperm, as well as decreased the movement speed of sperm. These adverse effects on sperm might be responsible for the reduced fertilization observed in offspring. In addition, treatment with E2 changed the levels of steroid hormones in zebrafish gonad and altered the transcriptional levels of genes associated with HPG axis, which is responsible for the regulation of germ cells maturation and gonadal development in zebrafish. Overall, these results suggested that treatment with environmentally relevant concentrations of E2 for the full sexual maturity cycle resulted in adverse effects on reproduction in zebrafish.

17β-Estradiol Induced Sex Reversal and Gonadal Transcriptome Analysis in the Oriental River Prawn (Macrobrachium nipponense): Mechanisms, Pathways, and Potential Harm

Sex reversal induced by 17β-estradiol (E₂) has shown the potential possibility for monoculture technology development. The present study aimed to determine whether dietary supplementation with different concentrations of E₂ could induce sex reversal in M. nipponense, and select the sex-related genes by performing the gonadal transcriptome analysis of normal male (M), normal female (FM), sex-reversed male prawns (RM), and unreversed male prawns (NRM). Histology, transcriptome analysis, and qPCR were performed to compare differences in gonad development, key metabolic pathways, and genes. Compared with the control, after 40 days, feeding E₂ with 200 mg/kg at PL25 (PL: post-larvae developmental stage) resulted in the highest sex ratio (female: male) of 2.22:1. Histological observations demonstrated the co-existence of testis and ovaries in the same prawn. Male prawns from the NRM group exhibited slower testis development without mature sperm. RNA sequencing revealed 3702 differentially expressed genes (DEGs) between M vs. FM, 3111 between M vs. RM, and 4978 between FM vs. NRM. Retinol metabolism and nucleotide excision repair pathways were identified as the key pathways for sex reversal and sperm maturation, respectively. Sperm gelatinase (SG) was not screened in M vs. NRM, corroborating the results of the slice D. In M vs. RM, reproduction-related genes such as cathepsin C (CatC), heat shock protein cognate (HSP), double-sex (Dsx), and gonadotropin-releasing hormone receptor (GnRH) were expressed differently from the other two groups, indicating that these are involved in the process of sex reversal. Exogenous E₂ can induce sex reversal, providing valuable evidence for the establishment of monoculture in this species.

Developmental toxicity of TCBPA on the nervous and cardiovascular systems of zebrafish (Danio rerio): A combination of transcriptomic and metabolomics

Tetrachlorobisphenol A (TCBPA), a widely used halogenated flame retardant, is frequently detected in environmental compartments and human samples. However, unknown developmental toxicity and mechanisms limit the entire understanding of its effects. In this study, zebrafish (Danio rerio) embryos were exposed to various concentrations of TCBPA while a combination of transcriptomics, behavioral and biochemical analyzes as well as metabolomics were applied to decipher its toxic effects and the potential mechanisms. We found that TCBPA could interfere with nervous and cardiovascular development through focal adhesion and extracellular matrix-receptor (ECM-receptor) interaction pathways through transcriptomic analysis. Behavioral and biochemical analysis results indicated abnormal swimming behavior of zebrafish larvae. Morphological observations revealed that TCBPA could cause the loss of head blood vessels. Metabolomic analysis showed that arginine-related metabolic pathways were one of the main pathways leading to TCBPA developmental toxicity. Our study demonstrated that by using omics, TCBPA was shown to have neurological and cardiovascular developmental toxicity and the underlying mechanisms were uncovered and major pathways identified.

Transcriptomes of Testes at Different Developmental Stages in the Opsariichthys bidens Predict Key Genes for Testis Development and Spermatogenesis

Testis development is a complex process involving multiple genes, and the molecular mechanisms underlying testis development in Opsariichthys bidens remain unclear. We performed transcriptome sequencing analysis on a total of 12 samples of testes from stages II, III, IV, and V of O. bidens and obtained a total of 79.52 Gb clean data, as well as 288,573 transcripts and 116,215 unigenes. Differential expression analysis showed that 22,857 differentially expressed genes (DEGs) were screened in six comparison groups (III vs. II, IV vs. II, V vs. II, IV vs. III, V vs. III, and V vs. IV). Kyoto Encyclopedia of Genes and Genomes enrichment analysis of DEGs showed that six comparison groups were significantly enriched for a total of 20 significantly up- or down-regulated pathways, including six pathways related to signal transduction, three pathways related to energy metabolism, five pathways related to disease, and two pathways related to ribosomes. Furthermore, our investigation revealed that DEGs were enriched in several important functional pathways, such as Huntington's disease signaling pathway, TGF-β signaling pathway, and ribosome signaling pathway. Protein-protein interaction network analysis of DEGs identified 63 up-regulated hub genes, including 9 kinesin genes and 2 cytoplasmic dynein genes, and 39 down-regulated hub genes, including 13 ribosomal protein genes. This result contributes to the knowledge of spermatogenesis and testis development in O. bidens.

Histology study and transcriptome analysis of the testis of Loach(Misgurnus anguillicaudatus) in response to phenanthrene exposure

Phenanthrene (PHE) is one of the most abundant polycyclic aromatic hydrocarbon compounds (PAHs) in the aquatic environment. The loaches were exposed at concentrations of 0.30、1.00、3.00 mg L^-1 for 60 days. The effects of PHE on the testis development were evaluated by calculating the survival rate, observing the structure of testis and analyzing transcriptome. Firstly, PHE markedly decreased the survival rate in a dose-dependent manner. Then, the number and density of spermatogonia, primary spermatocytes, secondary spermatocytes and spermatids were substantially reduced under PHE exposure. The space in the seminiferous tubule obviously increased in the high PHE concentration group. Meanwhile, transcriptome comparative analysis identified 5329 differentially expressed genes (DEGs) including 2928 up-regulated and 2401 down-regulated in the testis of loach exposed PHE for 60 days. Meiotic cell cycle, arganelle fission, ATPase activity and adenylate nucleotide binding were significantly differences by GO (Gene Ontology) enrichment. KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway analysis revealed that TNF (Tumor Necrosis Factor) signaling pathway, CAMs (Cell Adhesion Molecules), cytochrome P450 and lipid metabolism were markedly regulated. In addition, eight DEGs were randomly selected from the testis transcriptomics results for qPCR verification, the results were consistent with RNA-Seq. Finally, related genes (piwil2, dmc1, vasa, ubr2, dnd, rnf17, plcb2, c-fos, gpx4) of testis development were further confirmed and they were differentially regulated after PHE exposure. In summary, a survey of the mechanism of loach testis response to PHE was performed, and a large number of gene expression levels regarding metabolism, spermatogenesis and immunity genes were acquired from RNA-seq. This study provide informations for elucidating the molecular mechanism of PHE affected the testis development of loach.

Feminization of channel catfish with 17β-oestradiol involves methylation and expression of a specific set of genes independent of the sex determination region

Exogenous oestrogen 17β-oestradiol (E2) has been shown to effectively induce feminization in teleosts. However, the molecular mechanisms underlying the process remain unclear. Here, we determined global DNA methylation and gene expression profiles of channel catfish (Ictalurus punctatus) during early sex differentiation after E2 treatment. Overall, the levels of global DNA methylation after E2 treatment were not significantly different from those of controls. However, a specific set of genes were differentially methylated, which included many sex differentiation-related pathways, such as MARK signalling, adrenergic signalling, Wnt signalling, GnRH signalling, ErbB signalling, and ECM-receptor interactions. Many genes involved in these pathways were also differentially expressed after E2 treatment. Specifically, E2 treatments resulted in upregulation of female-related genes and downregulation of male-related genes in genetic males during sex reversal. However, E2-induced sex reversal did not cause sex-specific changes in methylation profiles or gene expression within the sex determination region (SDR) on chromosome 4, suggesting that E2-induced sex reversal was a downstream process independent of the sex determination process that was regulated by sex-specific methylation within the SDR.

Transcriptomic and targeted metabolomic analysis revealed the toxic effects of prochloraz on larval zebrafish

Prochloraz (PCZ), an imidazole fungicide, has been extensively used in horticulture and agriculture to protect against pests and diseases. To investigate the potential toxicity of PCZ on aquatic organisms, larval zebrafish, as a model, were exposed to a series of concentrations (0, 20, 100, and 500 μg/L) of PCZ for 7 days. With transcriptomic analysis, we found that exposure to high dose PCZ could produce 76 downregulated and 345 upregulated differential expression genes (DEGs). Bioinformatics analysis revealed that most of the DEGs were characterized in the pathways of glycolipid metabolism, amino acid metabolism and oxidative stress in larval zebrafish. Targeted metabolomic analysis was conducted to verify the effects of PCZ on the levels of acyl-carnitines and some amino acids in larval zebrafish. In addition, biochemical indicators related to glycolipid metabolism were affected obviously, manifested as elevated triglyceride (TG) levels and decreased glucose (Glu) levels in whole larvae. The expression levels of genes associated with glycolipid metabolism were affected in larvae after exposure to PCZ (PK, GK, PEPckc, SREBP, ACO). Interestingly, we further confirmed that PCZ could induce oxidative stress by the changing enzyme activities (T-GSH, GSSG) and upregulating several related genes levels in larval zebrafish. Generally, our results revealed that the endpoints related to glycolipid metabolism, amino acid metabolism and oxidative stress were influenced by PCZ in larval zebrafish.

Insights into the effects of difenoconazole on the livers in male mice at the biochemical and transcriptomic levels

Difenoconazole (DFZ) is a broad-spectrum triazole fungicide, that is extensively used in agriculture. Studies have shown that residues of DFZ and other fungicides have toxic effects on nontarget organisms. However, its hepatoxicity in mammals remains unclear. Here, we characterized the toxic hepatic effects in male C57BL/6 mice exposed to 30 and 100 mg/kg bw DFZ for 14 and 56 days, respectively. The results revealed that DFZ could increase the relative liver weights, however, the relative fat and spleen weights decreased. More importantly, DFZ exposure changed the hepatic morphology and induced hepatic oxidative stress. Gene expression analysis suggested that DFZ could induce a glycolipid metabolism disorder. Moreover, hepatic transcriptomic analysis revealed the effects of DFZ exposure on the transcriptional levels of various genes, and enrichment analysis of differentially expressed genes (DEGs) showed that energy metabolism and immune-associated pathways were mainly affected. We validated the results from transcriptomic analysis and found that some key genes related to energy metabolism were affected. In addition, flow cytometry showed that the CD3+/CD4+ and CD3+ /CD8+ levels declined in the spleen of mice. Taken together, these findings combined with transcriptome analysis highlighted that DFZ caused different endpoints in the liver, which could provide more evidence for investigating the toxic effects of DFZ in mammals.

Transcriptomic analysis of pituitary in female and male spotted scat (Scatophagus argus) after 17β-estradiol injection

Spotted scat (Scatophagus argus) is a popular species of marine fish cultured in China. It shows normal sexual growth dimorphism. Female spotted scat grows quicker and bigger than males. Growth and reproduction are the most important traits in aquaculture. In vertebrates, the pituitary gland occupies an important position in the growth and reproduction axis. Estrogen is involved in regulating growth and reproduction in the pituitary gland in an endocrine fashion. Transcriptome sequencing of the pituitary was performed in female and male fish at 6 h after 17β-estradiol injection (4.0 μg E₂/g body weight, BW). Compared with the pituitary of female and male groups, 144 and 64 genes [|log₂^{(fold change)}| ≥ 1.0 and false discovery rate (FDR) < 0.05] were significantly differentially expressed in E₂-injected females and males, respectively (p < 0.05). Of these, 59 and 48 were up-regulated, and 85 and 16 were down-regulated. According to the Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) pathway analyses, DEGs were involved in signal pathways, such as growth, reproduction, oocyte meiosis and steroid biosynthesis. Of these, estrogen affected the expression of some sex steroid synthesis and receptor genes in the pituitary gland through feedback, such as hsd17b7, pgr and cyp19a1b, regulating the reproductive activities. Besides, some growth-related genes, such as gap43, junbb, mstn2 and insm1a responded to estrogen. E₂ might affect the expression level of gh mRNA by regulating the expression levels of growth-related genes. Our results provide a theoretical basis for studying the molecular mechanism of growth and reproduction regulation at the pituitary level of spotted scat responded to E₂.

In vivo and in vitro short-term bisphenol A exposures disrupt testicular energy metabolism and negatively impact spermatogenesis in zebrafish

This study investigated the in vitro and short-term in vivo effects of Bisphenol A (BPA) on testicular energy metabolism and morphology in the zebrafish (Danio rerio). Testes were incubated in vitro for 1 h or fish were exposed in vivo to BPA in the tank water for 12 h. Testicular lactate, glycogen and cholesterol were measured and ¹⁴C-deoxy-d-glucose uptake and activity of lactate dehydrogenase (LDH), alanine aminotransferase (ALT) and aspartate aminotransferase (AST) were determined. In addition, testis samples from the in vivo exposures were subject to digital analysis of testicular cells using Ilastik software and the Pixel Classification module and estimation of apoptosis by Terminal deoxynucleotidyl transferase (TdT) dUTP Nick-End Labeling (TUNEL) immunohistochemical analysis. Our results from in vitro studies showed that BPA at 10 pM and 10 μM decreased testicular lactate content, glycogen content and LDH activity, but increased testicular AST activity. In addition, only BPA at 10 pM significantly decreased testicular ALT activity and cholesterol content. However, ¹⁴C-deoxy-d-glucose uptake was not changed. Furthermore, our results from in vivo studies showed that 10 pM BPA but not 10 μM BPA reduced testicular content of lactate and glycogen. In addition, both BPA concentrations decreased AST activity, whereas only BPA at 10 μM reduced ALT activity. However, LDH activity was not changed. Additionally, both concentrations of BPA induced spermatocyte apoptosis and a decrease in the proportion of the surface area of spermatids and spermatozoa. Collectively these data suggest that short-term BPA exposure affects energy metabolism and spermatogenesis in male zebrafish.

Developing a transcriptomic framework for testing testosterone-mediated handicap hypotheses

Sexually selected traits are hypothesized to be honest signals of individual quality due to the costs associated with their maintenance, development, and/or production. Testosterone, a sex steroid associated with the development and/or production of sexually selected traits, has been proposed to enforce the honesty of sexually selected traits via its immunosuppressive effects (i.e., the Immunocompetence Handicap Hypothesis) and/or by influencing an individual's exposure/susceptibility to oxidative stress (i.e., the Oxidation Handicap Hypothesis). Previous work testing these hypotheses has primarily focused on physiological measurements of immunity or oxidative stress, but little is known about the molecular pathways by which testosterone could influence immunity and/or oxidative stress pathways. To further understand the transcriptomic consequences of experimentally elevated testosterone in the context of handicap hypotheses, we used previously published RNA-seq data from studies that measured the transcriptome of individuals treated with either a testosterone-filled or an empty (i.e., control) implant. Two studies encompassing three species of bird and three tissue types fit our selection criteria and we reanalyzed the data using weighted gene co-expression network analysis. Testosterone-treated individuals exhibited signatures of immunosuppression and our results describe the molecular pathways underlying this effect. We also provide some evidence to suggest that the transcriptomic signature of immunosuppression is evolutionarily conserved between the three species. While our results provide no evidence to suggest testosterone mediates handicaps via pathways associated with oxidative stress, they do support the hypothesis that testosterone enforces the honesty of sexually-selected traits by influencing an individual's immunocompetence. Overall, this study develops a framework for testing testosterone-mediated handicap hypotheses and provides guidelines for future integrative and comparative studies focused on the proximate mechanisms mediating sexually selected traits.

Vitellogenins and choriogenins are biomarkers for monitoring Oryzias curvinotus juveniles exposed to 17 β - estradiol

The effect of estrogens on Oryzias curvinotus juveniles were investigated by sequencing the transcriptome of O. curvinotus juveniles exposed to 17 β - estradiol for 24 h. A total of 69,071,524 and 71,210,528 raw reads were obtained for the control group (NC) and 17 β - estradiol exposure group (E2), respectively. After de novo assembly, total 133,210 unigenes were identified, and 85,837 unigenes (64.44% of 133,210) were annotated. Analysis of the transcriptome showed that exposure to 2 μg/L 17 β - estradiol led to the up-regulation of 19 genes and down-regulation of 18 genes. The eef1b and rps4x was most suitable as controls for quantitative real-time PCR (qPCR) using Reffinder. Different expression genes enrichment analysis found that exposed to 2 μg/L 17 β - estradiol affected various physiological processes, including spliceosome, phototransduction, amino sugar and nuclear sugar metabolism, hypotaurine metabolism, and renin-angiotensin system, etc. Exposing O. curvinotus juveniles to increasing concentrations of 17 β - estradiol (2 ng/L, 20 ng/L, 200 ng/L and 2 μg/L) led to significant up-regulation of vitellogenins (vtgs) and choriogenins (chgs) mRNA expression. The present study is the first high-throughput transcriptome sequencing of O. curvinotus juveniles, which will be useful for future functional analysis of genes related to environmental estrogen exposed, and development of biomarkers.

Transcriptome alterations in female Daphnia (Daphnia magna) exposed to 17β-estradiol

The molecular mechanism of evaluating 17β-estradiol (E₂)-induced toxicity in female Daphnia magna has not been determined. In this study, the transcriptome of D. magna was analyzed after exposure to three different concentrations (0, 10, and 100 ng L^-1) of E₂ at 3, 6, and 12 h. The results showed 351-17,221 significantly up-regulated and 505-10,282 significantly down-regulated genes (P < 0.05). Overall, the selected largest 10,282 (10 ng L^-1vs control at 12 h) down-regulated and 17,221 (100 vs 10 ng L^-1) up-regulated genes were identified; following annotation, pathways in cancer and RNA transport were found to be enriched according to the interaction network. Among all completed comparisons, KEGG pathways related to the immune system, cancer, disease infection, and active compound metabolism were identified by short time series expression miner analysis. A different set of genes fluctuated in a "U"-shaped pattern over time and at different concentrations of E₂, whereas some genes associated with disintoxication showed a reverse "U"-shaped response as E₂ administration was increased. These results suggest that E₂ exposure caused transcriptional changes in the immune system, disintoxication, disease prevention, and the protein degradation pathway.

Polystyrene microplastic exposure disturbs hepatic glycolipid metabolism at the physiological, biochemical, and transcriptomic levels in adult zebrafish

Microplastics (MPs), which are new types of environmental pollutants, have recently received widespread attention worldwide. MPs can accumulate in the bodies of animals and in plants, and they can also enter the human body through the food chain. However, knowledge of the effects of MPs on the health of animals is still limited. In this experiment, adult male zebrafish were exposed to 20 or 100 μg/L of 5 μm polystyrene MP for 21 days in an attempt to determine the hepatic effects related to glycolipid metabolism at the biochemical and transcriptomic levels. It was found that body weight and condition factor decreased significantly in zebrafish after exposure to 20 and 100 μg/L polystyrene MP for 21 days. The transcription levels of major genes related to glycolipid metabolism decreased significantly in the liver. Correspondingly, the levels of major biochemical parameters, including Glu, pyruvic acid, α-ketoglutaric acid and IDH, were also decreased in the livers of exposed zebrafish, especially those in the 100 μg/L polystyrene MP-treated group. Moreover, the data on the hepatic transcriptome also confirmed that some genes related to fatty acid metabolism, amino acid metabolism and carbon metabolism tended to be decreased in the livers of exposed zebrafish. Taken together, our data confirmed that polystyrene PS-MP can induce hepatic glycolipid metabolism disorder at the physiological, biochemical, and transcriptomic levels in adult zebrafish after 21 days of exposure.

Cosmetics with hormonal composition for bioindicators Artemia salina L. and Allium cepa L. toxic potential

The emerging pollutants cover a wide range of synthetic chemicals that are indispensable to modern society but with little known effects for aquatic animals and for people who consume polluted waters with these products. Gels manipulated with hormones are widely used for hormone replacement, muscle growth, among other purposes. However, only a small part of these hormones are absorbed into the skin, and so these can be transferred to the domestic sewage during the washing of exposed body regions. Thus, the aim of this study was to evaluate the toxicity levels for the bioindicator Artemia salina L., and cytotoxicity and mutagenicity for the bioindicator Allium cepa L. of gels handled with 1% testosterone and 0.1% estradiol. Data from immobile/dead A. salina numbers after 24 h of exposure showed that the highest concentrations of testosterone (10 and 25 μg/mL) and the highest concentration of estradiol (15 μg/mL) were toxic to this bioindicator. For the bioindicator A. cepa, mitotic indices and chromosomal aberrations did not indicate statistical differences between the groups treated with the testosterone gels (1, 10, 50 μg/mL) and estradiol (0.03, 0.30, 1.5 μg/mL) and the control group. However, all concentrations of the testosterone-containing gel decreased the percentage of cell division in relation to the time 0 h of each treatment and to the time 24 h of the negative control. Therefore, it is concluded that rivers or aquatic environments can be polluted if wastewater with the toxic concentrations found of these hormonal gels is discarded without previous treatment, compromising the life of organisms that live there. And, it encourages the development of techniques for treating sewage and water to reduce/eliminate the hormones present in them.

Next generation sequencing data for use in risk assessment

Next generation sequencing (NGS) represents several powerful platforms that have revolutionized RNA and DNA analysis. The parallel sequencing of millions of DNA molecules can provide mechanistic insights into toxicology and provide new avenues for biomarker discovery with growing relevance for risk assessment. The evolution of NGS technologies has improved over the last decade with increased sensitivity and accuracy to foster new biomarker assays from tissue, blood and other biofluids. NGS sequencing technologies can identify transcriptional changes and genomic targets with base pair precision in response to chemical exposure. Further, there are several exciting movements within the toxicology community that incorporate NGS platforms into new strategies for more rapid toxicological characterizations. These include the Tox21 in vitro high throughput transcriptomic screening program, development of organotypic spheroids, alternative animal models, mining archival tissues, liquid biopsy and epigenomics. This review will describe NGS-based technologies, demonstrate how they can be used as tools for target discovery in tissue and blood, and suggest how they might be applied for risk assessment.

Molecular insights into ovary degeneration induced by environmental factors in female oriental river prawns Macrobrachium nipponense

The oriental river prawn, Macrobrachium nipponense, is an important breeding species in China. The ovary development of this prawn is regulated by the genetic factors and external environmental factors and has obvious seasonal regularity. However, the molecular mechanism of regulating ovary degradation in M. nipponense remains unclear. To address this issue, we performed transcriptome sequencing and gene expression analyses of eyestalks, cerebral ganglia (CG) and thoracic ganglia (TG) of female M. nipponense between the full ovary stage and degenerate ovary stage. Differentially expressed genes enrichment analysis results identified several important pathways such as "phototransduction-fly," "circadian rhythm-fly" and "steroid hormone biosynthesis secretion." In the period of ovarian degeneration, the expressions of Tim, Per2 and red pigment concentration hormone (RPCH) were significantly decreased in the eyestalk, CG and TG. And expression of 7 genes in the steroid synthesis pathway, including steryl-sulfatase, cytochrome P450 family 1 subfamily A polypeptide 1, estradiol 17β-dehydrogenase 2, glucuronosyltransferase, 3-oxo-5-alpha-steroid 4-dehydrogenase 1, estradiol 17-dehydrogenase 1 and estrone sulfotransferase was significantly decreased in the CG. Food and light signals affect the expression of clock genes and thereby decrease the expression of RPCH and the estradiol synthesis-related genes in the nervous system, which may be the main cause of ovarian degeneration in M. nipponense. The results will contribute to a better understanding of the molecular mechanisms of ovarian development regulation in crustaceans.