Gene co-expression network analysis reveals mechanisms underlying ozone-induced carbamazepine toxicity in zebrafish (Danio rerio) embryos

Developmental toxicity of the emerging contaminant cyclophosphamide and the integrated biomarker response (IBRv2) in zebrafish

The safety of cyclophosphamide (CP) in the early developmental stages is not studied yet; it is important to study the responses at these stages because they might have relevance to CP-administered humans. We studied the developmental toxicity of CP by analysing physiological, morphological, and oxidative stress, neurotransmission enzymes, gene expression and histological endpoints in zebrafish embryos/larvae. The study lasted for 120 hpf at environmentally relevant concentrations of CP. No visible alterations were noticed in the control group. Delayed hatching, slow heart rate, yolk sac oedema, pericardial oedema, morphological deformities, the incompetence of oxidative stress biomarkers, excessive generation of ROS, apoptosis, inhibition of neurotransmitters and histopathological anomalies were observed in CP-treated groups. These alterations were found to be concentration- and duration-dependent effects for physiological and morphological endpoints, whereas concentration-dependent effects were antioxidants, ROS, apoptosis and histological endpoints. Biomarkers and gene expression were standardised using the integrated biomarker response-IBRv2 index. The IBRv2 index showed a concentration-dependent behaviour. A non-lethal developmental and teratogenic effect was observed in CP-treated zebrafish embryos/larvae at the studied concentrations. The studied biomarkers are sensitive, and the responses are interrelated; thus, their responses are useful to assess veiled and unseen hazards of pharmaceuticals.

Removal of micropollutants and ecotoxicity during combined biological activated carbon and ozone (BO3) treatment

Ozonation is a viable option to improve the removal of micropollutants (MPs) in wastewater treatment plants (WWTPs). Nevertheless, the application of ozonation is hindered by its high energy requirements and by the uncertainties regarding the formation of toxic transformation products in the process. Energy requirements of ozonation can be reduced with a pre-ozone treatment, such as a biological activated carbon (BAC) filter, that removes part of the effluent organic matter before ozonation. This study investigated a combination of BAC filtration followed by ozonation (the BO3 process) to remove MPs at low ozone doses and low energy input, and focused on the formation of toxic organic and inorganic products during ozonation. Effluent from a WWTP was collected, spiked with MPs (approximately 1 µg/L) and treated with the BO3 process. Different flowrates (0.25-4 L/h) and specific ozone doses (0.2-0.6 g O3/g TOC) were tested and MPs, ecotoxicity and bromate were analyzed. For ecotoxicity assessment, three in vivo (daphnia, algae and bacteria) and six in vitro CALUX assays (Era, GR, PAH, P53, PR, andNrf2 CALUX) were used. Results show that the combination of BAC filtration and ozonation has higher MP removal and higher ecotoxicity removal than only BAC filtration and only ozonation. The in vivo assays show a low ecotoxicity in the initial WWTP effluent samples and no clear trend with increasing ozone doses, while most of the in vitro assays show a decrease in ecotoxicity with increasing ozone dose. This suggests that for the tested bioassays, feed water and ozone doses, the overall ecotoxicity of the formed transformation products during ozonation was lower than the overall ecotoxicity of the parent compounds. In the experiments with bromide spiking, relevant formation of bromate was observed above specific ozone doses of approximately 0.4 O3/g TOC and more bromate was formed for the samples with BAC pre-treatment. This indirectly indicates the effectivity of the pre-treatment in removing organic matter and making ozone more available to react with other compounds (such as MPs, but also bromide), but also underlines the importance of controlling the ozone dose to be below the threshold to avoid formation of bromate. It was concluded that treatment of the tested WWTP effluent in the BO3 process at a specific ozone dose of 0.2 g O3/g TOC, results in high MP removal at limited energy input while no increase in ecotoxicity, nor formation of bromate was observed under this condition. This indicates that the hybrid BO3 process can be implemented to remove MPs and improve the ecological quality of this WWTP effluent with a lower energy demand than conventional MP removal processes such as standalone ozonation.

A pituitary gene network linking vgll3 to regulators of sexual maturation in male Atlantic salmon

Age at maturity is a key life history trait and a significant contributor to life history strategy variation. The maturation process is influenced by genetic and environmental factors, but specific causes of variation in maturation timing remain elusive. In many species, the increase in the regulatory gonadotropin-releasing hormone 1 (GnRH1) marks the onset of puberty. Atlantic salmon, however, lacks the gnrh1 gene, suggesting gnrh3 and/or other regulatory factors are involved in the maturation process. Earlier research in Atlantic salmon has found a strong association between alternative alleles of vgll3 and maturation timing. Recently we reported strong induction of gonadotropin genes (fshb and lhb) in the pituitary of Atlantic salmon homozygous for the early maturation allele (E) of vgll3. The induction of gonadotropins was accompanied by increased expression of their direct upstream regulators, c-jun and sf1 (nr5a1b) but the regulatory connection between vgll3 and these regulators has never been investigated in any organism. In this study, we investigated the potential regulatory connection between vgll3 genotypes and these regulators through a stepwise approach of identifying a gene regulatory network (GRN) containing c-jun and sf1, and transcription factor motif enrichment analysis. We found a GRN containing c-jun with predicted upstream regulators, e2f1, egr1, foxj1 and klf4, to be differentially expressed in the pituitary. Finally, we suggest a vgll3 and Hippo pathway -dependent model for transcriptional regulation of c-jun and sf1 in the pituitary, which may have broader implications across vertebrates.

Concentration-response gene expression analysis in zebrafish reveals phenotypically-anchored transcriptional responses to retene

Polycyclic aromatic hydrocarbons (PAHs) are ubiquitous environmental contaminants and are associated with human disease. Canonically, many PAHs induce toxicity via activation of the aryl hydrocarbon receptor (AHR) pathway. While the interaction between PAHs and the AHR is well-established, understanding which AHR-regulated transcriptional effects directly result in observable phenotypes and which are adaptive or benign is important to better understand PAH toxicity. Retene is a frequently detected PAH in environmental sampling and has been associated with AHR2-dependent developmental toxicity in zebrafish, though its mechanism of toxicity has not been fully elucidated. To interrogate transcriptional changes causally associated with retene toxicity, we conducted whole-animal RNA sequencing at 48 h post-fertilization after exposure to eight retene concentrations. We aimed to identify the most sensitive transcriptomic responses and to determine whether this approach could uncover gene sets uniquely differentially expressed at concentrations which induce a phenotype. We identified a concentration-response relationship for differential gene expression in both number of differentially expressed genes (DEGs) and magnitude of expression change. Elevated expression of cyp1a at retene concentrations below the threshold for teratogenicity suggested that while cyp1a expression is a sensitive biomarker of AHR activation, it may be too sensitive to serve as a biomarker of teratogenicity. Genes differentially expressed at only non-teratogenic concentrations were enriched for transforming growth factor-β (TGF-β) signaling pathway disruption while DEGs identified at only teratogenic concentrations were significantly enriched for response to xenobiotic stimulus and reduction-oxidation reaction activity. DEGs which spanned both non-teratogenic and teratogenic concentrations showed similar disrupted biological processes to those unique to teratogenic concentrations, indicating these processes were disrupted at low exposure concentrations. Gene co-expression network analysis identified several gene modules, including those associated with PAHs and AHR2 activation. One, Module 7, was strongly enriched for AHR2-associated genes and contained the strongest responses to retene. Benchmark concentration (BMC) of Module seven genes identified a median BMC of 7.5 µM, nearly the highest retene concentration with no associated teratogenicity, supporting the hypothesis that Module seven genes are largely responsible for retene toxicity.

Developmental Toxicity and Biotransformation of Two Anti-Epileptics in Zebrafish Embryos and Early Larvae

The zebrafish (Danio rerio) embryo is gaining interest as a bridging tool between in-vitro and in-vivo developmental toxicity studies. However, cytochrome P450 (CYP)-mediated drug metabolism in this model is still under debate. Therefore, we investigated the potential of zebrafish embryos and larvae to bioactivate two known anti-epileptics, carbamazepine (CBZ) and phenytoin (PHE), to carbamazepine-10,11-epoxide (E-CBZ) and 5-(4-hydroxyphenyl)-5-phenylhydantoin (HPPH), respectively. First, zebrafish were exposed to CBZ, PHE, E-CBZ and HPPH from 5¼- to 120-h post fertilization (hpf) and morphologically evaluated. Second, the formations of E-CBZ and HPPH were assessed in culture medium and in whole-embryo extracts at different time points by targeted LC-MS. Finally, E-CBZ and HPPH formation was also assessed in adult zebrafish liver microsomes and compared with those of human, rat, and rabbit. The present study showed teratogenic effects for CBZ and PHE, but not for E-CBZ and HPPH. No HPPH was detected during organogenesis and E-CBZ was only formed at the end of organogenesis. E-CBZ and HPPH formation was also very low-to-negligible in adult zebrafish compared with the mammalian species. As such, other metabolic pathways than those of mammals are involved in the bioactivation of CBZ and PHE, or, these anti-epileptics are teratogens and do not require bioactivation in the zebrafish.