Single-cell RNA sequencing reveals tissue-specific transcriptomic changes induced by perfluorooctanesulfonic acid (PFOS) in larval zebrafish (Danio rerio)

Perfluorooctanesulfonic acid (PFOS) elicits adverse effects on numerous organs and developmental processes but the mechanisms underlying these effects are not well understood. Here, we use single-cell RNA-sequencing to assess tissue-specific transcriptomic changes in zebrafish (Danio rerio) larvae exposed to 16 µM PFOS or dimethylsulfoxide (0.01 %) from 3-72 h post fertilization (hpf). Data analysis was multi-pronged and included pseudo-bulk, untargeted clustering, informed pathway queries, and a cluster curated for hepatocyte biomarkers (fabp10a, and apoa2). Overall, 8.63 % (2390/27698) genes were significantly differentially expressed. Results from untargeted analysis revealed 22 distinct clusters that were manually annotated to specific tissues using a weight-of-evidence approach. The clusters with the highest number of significant differentially expressed genes (DEGs) were digestive organs, muscle, and otolith. Additionally, we assessed the distribution of pathway-specific genes known to be involved in PFOS toxicity: the PPAR pathway, β-oxidation of fatty acids, the Nfe2l2 pathway, and epigenetic modifications by DNA methylation, across clusters and identified the blood-related tissue to be the most sensitive. The curated hepatocyte cluster showed 220 significant DEGs and was enriched for the Notch signaling pathway. These findings provide insights into both established and novel sensitive target tissues and molecular mechanisms of developmental toxicity of PFOS.

High-Resolution Magic-Angle Spinning Nuclear Magnetic Resonance Identifies Impairment of Metabolism by T-2 Toxin, in Relation to Toxicity, in Zebrafish Embryo Model

Among the widespread trichothecene mycotoxins, T-2 toxin is considered the most toxic congener. In the present study, we utilized high-resolution magic-angle spinning nuclear magnetic resonance (HRMAS NMR), coupled to the zebrafish (Danio rerio) embryo model, as a toxicometabolomics approach to elucidate the cellular, molecular and biochemical pathways associated with T-2 toxicity. Aligned with previous studies in the zebrafish embryo model, exposure to T-2 toxin was lethal in the high parts-per-billion (ppb) range, with a median lethal concentration (LC50) of 105 ppb. Exposure to the toxins was, furthermore, associated with system-specific alterations in the production of reactive oxygen species (ROS), including decreased ROS production in the liver and increased ROS in the brain region, in the exposed embryos. Moreover, metabolic profiling based on HRMAS NMR revealed the modulation of numerous, interrelated metabolites, specifically including those associated with (1) phase I and II detoxification, and antioxidant pathways; (2) disruption of the phosphocholine lipids of cell membranes; (3) mitochondrial energy metabolism, including apparent disruption of the tricarboxylic acid (TCA) cycle, and the electron transport chain of oxidative phosphorylation, as well as "upstream" effects on carbohydrate, i.e., glucose metabolism; and (4) several compensatory catabolic pathways. Taken together, these observations enabled development of an integrated, system-level model of T-2 toxicity in relation to human and animal health.

Triphenyl Phosphate Alters Methyltransferase Expression and Induces Genome-Wide Aberrant DNA Methylation in Zebrafish Larvae

Emerging environmental contaminants, organophosphate flame retardants (OPFRs), pose significant threats to ecosystems and human health. Despite numerous studies reporting the toxic effects of OPFRs, research on their epigenetic alterations remains limited. In this study, we investigated the effects of exposure to 2-ethylhexyl diphenyl phosphate (EHDPP), tricresyl phosphate (TMPP), and triphenyl phosphate (TPHP) on DNA methylation patterns during zebrafish embryonic development. We assessed general toxicity and morphological changes, measured global DNA methylation and hydroxymethylation levels, and evaluated DNA methyltransferase (DNMT) enzyme activity, as well as mRNA expression of DNMTs and ten-eleven translocation (TET) methylcytosine dioxygenase genes. Additionally, we analyzed genome-wide methylation patterns in zebrafish larvae using reduced-representation bisulfite sequencing. Our morphological assessment revealed no general toxicity, but a statistically significant yet subtle decrease in body length following exposure to TMPP and EHDPP, along with a reduction in head height after TPHP exposure, was observed. Eye diameter and head width were unaffected by any of the OPFRs. There were no significant changes in global DNA methylation levels in any exposure group, and TMPP showed no clear effect on DNMT expression. However, EHDPP significantly decreased only DNMT1 expression, while TPHP exposure reduced the expression of several DNMT orthologues and TETs in zebrafish larvae, leading to genome-wide aberrant DNA methylation. Differential methylation occurred primarily in introns (43%) and intergenic regions (37%), with 9% and 10% occurring in exons and promoter regions, respectively. Pathway enrichment analysis of differentially methylated region-associated genes indicated that TPHP exposure enhanced several biological and molecular functions corresponding to metabolism and neurological development. KEGG enrichment analysis further revealed TPHP-mediated potential effects on several signaling pathways including TGFβ, cytokine, and insulin signaling. This study identifies specific changes in DNA methylation in zebrafish larvae after TPHP exposure and brings novel insights into the epigenetic mode of action of TPHP.

A review on aquatic toxins - Do we really know it all regarding the environmental risk posed by phytoplankton neurotoxins?

Aquatic toxins are potent natural toxins produced by certain cyanobacteria and marine algae species during harmful cyanobacterial and algal blooms (CyanoHABs and HABs, respectively). These harmful bloom events and the toxins produced during these events are a human and environmental health concern worldwide, with occurrence, frequency and severity of CyanoHABs and HABs being predicted to keep increasing due to ongoing climate change scenarios. These contexts, as well as human health consequences of some toxins produced during bloom events have been thoroughly reviewed before. Conversely, the wider picture that includes the non-human biota in the assessment of noxious effects of toxins is much less covered in the literature and barely covered by review works. Despite direct human exposure to aquatic toxins and related deleterious effects being responsible for the majority of the public attention to the blooms' problematic, it constitutes a very limited fraction of the real environmental risk posed by these toxins. The disruption of ecological and trophic interactions caused by these toxins in the aquatic biota building on deleterious effects they may induce in different species is paramount as a modulator of the overall magnitude of the environmental risk potentially involved, thus necessarily constraining the quality and efficiency of the management strategies that should be placed. In this way, this review aims at updating and consolidating current knowledge regarding the adverse effects of aquatic toxins, attempting to going beyond their main toxicity pathways in human and related models' health, i.e., also focusing on ecologically relevant model organisms. For conciseness and considering the severity in terms of documented human health risks as a reference, we restricted the detailed revision work to neurotoxic cyanotoxins and marine toxins. This comprehensive revision of the systemic effects of aquatic neurotoxins provides a broad overview of the exposure and the hazard that these compounds pose to human and environmental health. Regulatory approaches they are given worldwide, as well as (eco)toxicity data available were hence thoroughly reviewed. Critical research gaps were identified particularly regarding (i) the toxic effects other than those typical of the recognized disease/disorder each toxin causes following acute exposure in humans and also in other biota; and (ii) alternative detection tools capable of being early-warning signals for aquatic toxins occurrence and therefore provide better human and environmental safety insurance. Future directions on aquatic toxins research are discussed in face of the existent knowledge, with particular emphasis on the much-needed development and implementation of effective alternative (eco)toxicological biomarkers for these toxins. The wide-spanning approach followed herein will hopefully stimulate future research more broadly addressing the environmental hazardous potential of aquatic toxins.

From multiallele fish to nonstandard environments, how ZFIN assigns phenotypes, human disease models, and gene expression annotations to genes

Danio rerio is a model organism used to investigate vertebrate development. Manipulation of the zebrafish genome and resultant gene products by mutation or targeted knockdown has made the zebrafish a good system for investigating gene function, providing a resource to investigate genetic contributors to phenotype and human disease. Phenotypic outcomes can be the result of gene mutation, targeted knockdown of gene products, manipulation of experimental conditions, or any combination thereof. Zebrafish have been used in various genetic and chemical screens to identify genetic and environmental contributors to phenotype and disease outcomes. The Zebrafish Information Network (ZFIN, zfin.org) is the central repository for genetic, genomic, and phenotypic data that result from research using D. rerio. Here we describe how ZFIN annotates phenotype, expression, and disease model data across various experimental designs, how we computationally determine wild-type gene expression, the phenotypic gene, and how these results allow us to propagate gene expression, phenotype, and disease model data to the correct gene, or gene related entity.

Enrichment of nucleobase adducts from genomic DNA in the cytoplasm by solid-phase extraction

The exact levels of some DNA adducts, like N7-deoxyguanosine (N7-dG), can be under-calculated since these adducts may depurinate due to their chemical instability, leading to corresponding nucleobase adducts being released into the cytoplasm. To accurately quantify the levels of DNA adducts, it is necessary to consider those modified nucleobases. However, high levels and diversity of cytoplasmic small molecule metabolites (SMMs) can strongly interfere with the detection of adducts, and it is almost impossible to remove them with nucleobase adducts being well retained. Therefore, we aimed to establish an optimized enrichment method based on solid-phase extraction (SPE) to reduce the co-elution of SMMs with target analytes. In this vein, we employed three bisphenols (BPA, BPF, and BPAF) as examples, prepared corresponding N7-guanine (N7-Gua) adducts, loaded on an Oasis hydrophilic-lipophilic balance ® (HLB) cartridge, used a series of mobile phases containing different percentage of methanol for elution, and evaluated the levels of these adducts in each eluent. First, we found that neutral samples led to the best retention for all three adducts compared with acidified or basified ones. We next employed normal distribution fitting model to characterize the elution of analytes from H₂O/methanol with different pHs and observed that neutral mobile phases resulted in more hydrophobic elution for all three adducts. Besides, N7-BPA-Gua and N7-BPF-Gua obtained narrow fitted peaks at neutral pH, while N7-BPAF-Gua had minimized elution windows at low pH. After optimization, we exposed 293T cells to the aforementioned bisphenols and quantified the N7-Gua adducts in the cytoplasm and the corresponding N7-dG adducts in genomic DNA. The results revealed that with the same levels of BPs exposure, BPAF led to the highest levels of adducts in both cytoplasm and genomic DNA samples, followed by BPA and BPF in order. In summary, our research established an appropriate model to describe the elution of DNA adducts in the SPE, applied it to optimize the loading and elution conditions, and discussed the genotoxicity of bisphenols by accurate quantification of both cleaved and uncleaved N7-dG adducts.

Roles of H19/miR-29a-3p/COL1A1 axis in COE-induced lung cancer

Occupational lung cancer caused by coke oven emissions (COE) has attracted increasing attention, but the mechanism is not clear. Many evidences show ceRNA (competing endogenous RNA) networks play important regulatory roles in cancers. In this study, we aimed to construct and verify the ceRNA regulatory network in the occurrence of COE-induced lung squamous cell carcinoma (LUSC). We performed RNA sequencing with lung bronchial epithelial cell (16HBE) and COE induced malignant transformed cell (Rf). Furthermore, we analyzed RNA sequencing data of LUSC and adjacent tissues in the cancer genome atlas (TCGA) database. Combined our data and TCGA data to determine the differentially expressed lncRNAs, miRNAs, mRNAs. lncBASE, miRDB and miRTarBase were used to predict the binding relationship between lncRNA and miRNA, miRNA and mRNA. Based on these, we construct the ceRNA network. FREMSA, dual-luciferase reporter assay, quantitative real-time PCR (qRT-PCR), western-blot were used to verify the regulatory axis. CCK8 assay, phalloidin staining, p53 detection were used to explore the roles of this axis in the COE induced malignant transformation. Results showed 7 lncRNAs, 7 miRNAs and 146 mRNAs were identified. Among these, we constructed a ceRNA network including 1 lncRNA, 2 miRNAs and 9 mRNAs. Further verification confirmed the trend of lncRNA H19, miR-29a-3p and COL1A1 were consistent with sequencing results. H19 and COL1A1 were significantly higher in Rf than in 16HBE and miR-29a-3p was reverse. Regulatory investigation revealed H19 increased COL1A1 expression by sponging miR-29a-3p. Knockdown of H19, COL1A1 or overexpression of miR-29a-3p in Rf cells could inhibit cell proliferation, increased cell adhesion and p53 level. However, knockdown of H19 while suppressing the miR-29a-3p partially rescue the malignant phenotype of Rf caused by H19. In conclusion, all these indicated H19 functioned as a ceRNA to increase COL1A1 by sponging miR-29a-3p and promoted COE-induced cell malignant transformation.

H2O2/DEM-Promoted Maft Promoter Demethylation Drives Nrf2/ARE Activation in Zebrafish

The Nrf2/ARE signaling pathway is a cell survival response pathway in response to environmental stresses. The Nrf2/ARE signaling pathway can be activated by stimulating cysteine residues at different positions in the Keap1. However, the epigenetic mechanisms of the Nrf2/ARE pathway under different stimuli are still poorly understood. In this study, we found that both hydrogen peroxide (H₂O₂) and Diethyl Maleate (DEM) activated the Nrf2/ARE signaling pathway at 120 hpf in zebrafish. H₂O₂ regulated the demethylation of the maft promoter by inhibiting the expression of methyltransferase. This promotes the mRNA expression of the Nrf2 binding factor maft, thereby promoting the downstream antioxidant genes. The methylation of the Nrf2/ARE signaling pathway was not significantly regulated by DEM. However, under oxidative stress, the methyltransferase inhibitors (decitabine and azacitidine) demethylated the promoter region of maft. It activated the expression of the maft, further improving the Nrf2/ARE signal pathway. At last, antioxidant target genes were activated. It was shown that H₂O₂ and DEM cooperated with methyltransferase inhibitors, providing an important reference for the treatment of oxidative stress-related diseases and breaking new ground for the study of the mechanism of methyltransferase inhibitors in the process of tumor chemotherapy.

Effect of aerobic exercise as a treatment on type 2 diabetes mellitus with depression-like behavior zebrafish

BACKGROUND

Depression is the most known complication of type 2 diabetes mellitus (T2DM). Aerobic exercise improves glycemic control in T2DM, although the underlying mechanisms of comorbid depression-like behaviors in T2DM have not yet been fully elucidated.

METHODS

120 zebrafish were randomly assigned to four groups: Control, T2DM, T2DM + metformin, and T2DM + aerobic exercise. Then, all animals except the control group were fed with high glucose fairy shrimp (~40 g/kg/day) and exposed reserpine (40 μg/ml for 20 min) for 10 days. Here, behavioral tests were used for model verification. Following the verification, all groups were treated as before. Additionally, the T2DM + metformin group received metformin (~10.6 mg/kg/day) at the same time, while the T2DM + aerobic exercise group received aerobic exercise 30 min/day. Finally, blood glucose and behavioral tests, as well as protein and molecular levels were determined at Day 11 and 12.

RESULTS

Aerobic exercise alleviated depressive-like behavior and enhanced the levels of antidepressant biomarkers (NE, 5-HIAA) in zebrafish after 10 consecutive days of exercise. Additionally, 10 consecutive days of aerobic exercise decreased the levels of inflammatory biomarkers (IFN-γ, IL-1, IL-4) and depressive biomarkers (cortisol). Meanwhile, it also aided in the reduction of CD11b, IL-6, IL-6R, and caspase-3 expression to combat the neuroinflammation induced by T2DM, mediated the BDNF-TrkB pathway, and increased Bcl-2/Bax levels.

CONCLUSION

Given the remarkable similarity in neurochemistry between humans and zebrafish, this study supports the effectiveness of aerobic exercise as clinical guidance in preventing and treating T2DM complicated with depression.

CircRNA-IGLL1/miR-15a/RNF43 axis mediates ammonia-induced autophagy in broilers jejunum via Wnt/β-catenin pathway

With the continued increase of global ammonia emission, the damage to human or animal caused by ammonia pollution has attracted wide attention. The noncoding RNAs have been reported to regulate a variety of biological processes under different environmental stimulation via ceRNA (competing endogenous RNA) networks. Autophagy is a hallmark of tissue damage from air pollution. However, the specific role of circular RNAs (circRNAs) in the injury of intestinal tissue caused by autophagy remains unclear. Here, we established 42-days old ammonia-exposed broiler models and observed that autophagy flux in broiler jejunum was activated under ammonia exposure. Meanwhile, a total of eight significantly dysregulated expressed circRNAs were obtained and a circRNAs-miRNAs-genes interaction networks were constructed by bioinformatics analysis. Furthermore, an axis named circRNA-IGLL1/miR-15a/RNF43 was predicted to participate in the excessive autophagy by targeting RNF43. The target relationship was proved by dual-luciferase reporter assay in vitro. Mechanistically, downregulated circRNA-IGLL1 could suppress the expression of RNF43 in ammonia-exposed jejunum and the Wnt/β-catenin pathway was activated. Inhibition of miR-15a reversed autophagy caused by downregulated circRNA-IGLL1. CircRNA-IGLL1 could competitively bind miR-15a to regulate RNF43 expression, thus modulating the occurrence of autophagy. Taken together, our results showed that circRNA-IGLL1/miR-15a/RNF43 axis is involved in ammonia-induced intestinal autophagy in broilers.

Lethal and Sub-Lethal Effects and Modulation of Gene Expression Induced by T Kinase Inhibitors in Zebrafish (Danio Rerio) Embryos

Tyrosine kinase inhibitors (TKIs) are designed for targeted cancer therapy. The consumption of these drugs during the last 20 years has been constantly rising. In the zebrafish (Danio rerio) embryo toxicity test, we assessed the toxicity of six TKIs: imatinib mesylate, erlotinib, nilotinib, dasatinib, sorafenib and regorafenib. Imatinib mesylate and dasatinib induced lethal effects, while regorafenib, sorfenib and dasatinib caused a significant increase of sub-lethal effects, predominantly oedema, no blood circulation and formation of blood aggregates. The analyses of the changes in the expression of selected genes associated with the hormone system after the exposure to imatinib mesylate, dasatinib and regorafenib demonstrated that all three tested TKIs deregulated the expression of oestrogen receptor esr1, cytochrome P450 aromatase (cypa19b) and hydroxysteroid-dehydrogenase (hsd3b), regorafenib, and also thyroglobulin (tg). The expression of genes involved in the DNA damage response (gadd45 and mcm6) and apoptosis (bcl2) was deregulated only by exposure to regorafenib. The data indicate that common mechanisms, namely antiangiogenic activity and interference with steroidogenesis are involved in the TKI induced sub-lethal effects and potential hormone disrupting activity, respectively. The residues of TKIs may represent an environmental hazard; therefore, further ecotoxicological studies focusing also on the effects of their mixtures are warranted.

Advantages of omics technology for evaluating cadmium toxicity in zebrafish

In the last decade, several advancements have been made in omics technologies and they have been applied extensively in diverse research areas. Especially in toxicological research, omics technology can efficiently and accurately generate relevant data on the molecular dynamics associated with adverse outcomes. Toxicomics is defined as the combination of toxicology and omics technologies and encompasses toxicogenomics, toxicoproteomics, and toxicometabolomics. This paper reviews the trend of applying omics technologies to evaluate cadmium (Cd) toxicity in zebrafish (D. rerio). Cd is a toxic heavy metal posing several environmental concerns; however, it is being used widely in everyday life. Zebrafish embryos and larvae are employed as standard models for many toxicity tests because they share 71.4% genetic homology with humans. This study summarizes the toxicity of Cd on the nerves, liver, heart, skeleton, etc. of zebrafish and introduces detailed omics techniques to understand the results of the toxicomic studies. Finally, the trend of toxicity evaluation in the zebrafish model of Cd based on omics technology is presented.

An investigation into the biological effects of indirect potable reuse water using zebrafish embryos

Advanced treatment technologies are being assessed as a proactive measure to assist with the transformation of treated wastewater into a source of water for potable water production. We investigated the biological effects along an advanced water treatment pilot plant, using zebrafish embryos throughout early development. The study compared phenotypic observations with global transcriptome responses, enabling us to keep an open mind about the chemicals that might influence the biological activity. There was no evidence of acute toxicity at any treatment stage, but skeletal, cardiovascular and pigmentation changes occurred in a small proportion of embryos along the treatment process, and in a tap water; not detected in the aquarium water control. Reverse osmosis (RO) reduced the concentration of measured chemical contaminants in the water the most, while eliminating the occurrence of abnormalities detected in fish embryos. Conversely, advanced oxidation reversed the benefits of RO treatment by increasing the frequency of teratogenic and sub-lethal abnormalities seen. Using the molecular responses of zebrafish embryos to different IPR water, we report the bioactivity within the water at different stages of advanced treatment and associate these to perturbed biological functions. Transcriptomic analysis revealed alterations to the retinoid system, which was consistent with the observed teratogenic effects. Changes to tryptophan metabolism (associated with the production of melatonin required for the control of normal circadian rhythms) and somatolactin-beta (associated with normal pigmentation in fish) were also found. We show that underexplored forms of biological activity occur in treated wastewater effluent, and/or may be created depending on the type of advanced treatment process used. By integrating the available analytical chemistry we highlight chemical groups associated to this response. Our study shows that more detailed and in-depth characterisation of chemicals and biological pathways associated with advanced treatment water systems are needed to mitigate possible risks to downstream organisms.

Gene co-expression network analysis in zebrafish reveals chemical class specific modules

BACKGROUND

Zebrafish is a popular animal model used for high-throughput screening of chemical hazards, however, investigations of transcriptomic mechanisms of toxicity are still needed. Here, our goal was to identify genes and biological pathways that Aryl Hydrocarbon Receptor 2 (AHR2) Activators and flame retardant chemicals (FRCs) alter in developing zebrafish. Taking advantage of a compendium of phenotypically-anchored RNA sequencing data collected from 48-h post fertilization (hpf) zebrafish, we inferred a co-expression network that grouped genes based on their transcriptional response.

RESULTS

Genes responding to the FRCs and AHR2 Activators localized to distinct regions of the network, with FRCs inducing a broader response related to neurobehavior. AHR2 Activators centered in one region related to chemical stress responses. We also discovered several highly co-expressed genes in this module, including cyp1a, and we subsequently show that these genes are definitively within the AHR2 signaling pathway. Systematic removal of the two chemical types from the data, and analysis of network changes identified neurogenesis associated with FRCs, and regulation of vascular development associated with both chemical classes. We also identified highly connected genes responding specifically to each class that are potential biomarkers of exposure.

CONCLUSIONS

Overall, we created the first zebrafish chemical-specific gene co-expression network illuminating how chemicals alter the transcriptome relative to each other. In addition to our conclusions regarding FRCs and AHR2 Activators, our network can be leveraged by other studies investigating chemical mechanisms of toxicity.

Stem cells of aquatic invertebrates as an advanced tool for assessing ecotoxicological impacts

Environmental stressors are assessed through methods that quantify their impacts on a wide range of metrics including species density, growth rates, reproduction, behaviour and physiology, as on host-pathogen interactions and immunocompetence. Environmental stress may induce additional sublethal effects, like mutations and epigenetic signatures affecting offspring via germline mediated transgenerational inheritance, shaping phenotypic plasticity, increasing disease susceptibility, tissue pathologies, changes in social behaviour and biological invasions. The growing diversity of pollutants released into aquatic environments requires the development of a reliable, standardised and 3R (replacement, reduction and refinement of animals in research) compliant in vitro toolbox. The tools have to be in line with REACH regulation 1907/2006/EC, aiming to improve strategies for potential ecotoxicological risks assessment and monitoring of chemicals threatening human health and aquatic environments. Aquatic invertebrates' adult stem cells (ASCs) are numerous and can be pluripotent, as illustrated by high regeneration ability documented in many of these taxa. This is of further importance as in many aquatic invertebrate taxa, ASCs are able to differentiate into germ cells. Here we propose that ASCs from key aquatic invertebrates may be harnessed for applicable and standardised new tests in ecotoxicology. As part of this approach, a battery of modern techniques and endpoints are proposed to be tested for their ability to correctly identify environmental stresses posed by emerging contaminants in aquatic environments. Consequently, we briefly describe the current status of the available toxicity testing and biota-based monitoring strategies in aquatic environmental ecotoxicology and highlight some of the associated open issues such as replicability, consistency and reliability in the outcomes, for understanding and assessing the impacts of various chemicals on organisms and on the entire aquatic environment. Following this, we describe the benefits of aquatic invertebrate ASC-based tools for better addressing ecotoxicological questions, along with the current obstacles and possible overhaul approaches.

Multiple toxicity of propineb in developing zebrafish embryos: Neurotoxicity, vascular toxicity, and notochord defects in normal vertebrate development

A dithiocarbamate (DTC) fungicide, propineb, affects thyroid function and exerts immunotoxicity, cytotoxicity, and neurotoxicity in humans. Long-term exposure to propineb is associated with carcinogenicity, teratogenicity, malfunction of the reproductive system, and abnormalities in vital signs during organ development. However, there is no evidence of acute toxicity attributable to propineb in zebrafish. Therefore, in the present study, we assessed the toxicity of propineb in zebrafish by studying its adverse effects on embryo development, angiogenesis, and notochord development. Embryos with propineb exposure developed morphological and physiological defects and in larvae, apoptosis and notochord defects were induced in the early development stage. Transgenic fli1:eGFP zebrafish exposed to propineb showed abnormal larval development with defects in angiogenesis and deformed vasculature. Propineb induced irreversible damage to the neural development of embryos and neurogenic defects in developing zebrafish in transgenic olig2:dsRED zebrafish. These results show that exposure to propineb triggers abnormalities in different organ systems of zebrafish and suggests the physiological complexity of the response to propineb.

DNA methylome signatures as epigenetic biomarkers of hexanal associated with lung toxicity

Background

Numerous studies have investigated the relationship of environmental exposure, epigenetic effects, and human diseases. These linkages may contribute to the potential toxicity mechanisms of environmental chemicals. Here, we investigated the epigenetic pulmonary response of hexanal, a major indoor irritant, following inhalation exposure in F-344 rats.

Methods

Based on DNA methylation profiling in gene promoter regions, we identified hexanal-characterized methylated sites and target genes using an unpaired t-test with a fold-change cutoff of ≥ 3.0 and a p-value < 0.05. We also conducted an integrated analysis of DNA methylation and mRNA expression data to identify core anti-correlated target genes of hexanal exposure. To further investigate the potential key biological processes and pathways of core DNA methylated target genes, Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis were performed.

Results

Thirty-six dose-dependent methylated genes and anti-correlated target genes of DNA methylation and mRNA in lung tissue of hexanal exposed F-344 rats were identified. These genes were involved in diverse biological processes such as neuroactive ligand-receptor interaction, protein kinase cascade, and intracellular signaling cascade associated with pulmonary toxicity. These results suggest that novel DNA methylation-based epigenetic biomarkers of exposure to hexanal and elucidate the potential pulmonary toxicological mechanisms of action of hexanal.

Deregulation of whole-transcriptome gene expression in zebrafish (Danio rerio) after chronic exposure to low doses of imatinib mesylate in a complete life cycle study

Imatinib mesylate (IM) is an anticancer drug that belongs to tyrosine kinase inhibitors. We report the results of the first investigation of the chronic exposure of zebrafish (Danio rerio) to IM. The exposure to IM (0.01, 1 and 100 μg/L) was initiated in adult fish and continued through hatching and the offspring generation for seven months. In addition to standard toxicological endpoints, induction of genotoxic effects and whole-genome transcriptome of liver samples of offspring generation of zebrafish were analysed. Exposure to IM did not affect the survival and growth of zebrafish, did not cause any histopathological changes, but it induced a marginal increase in the chromosomal damage in blood cells. The whole-genome transcriptome analyses demonstrated dose-dependent increase in the number of differentially expressed genes with a significantly higher number of deregulated genes in female fish compared to male. Differentially expressed genes included genes involved in response to DNA damage, cell cycle control and regulation of circadian rhythm. Based on the low genotoxic activity and the pattern of the changes in DNA damage responsive genes we consider that at current environmental exposure levels, IM represents low risk for genotoxic effects in aquatic organisms. Exposure to IM also induced deregulation of the expression of genes associated with steroidogenesis and hormone metabolism and function, which indicates hormone-disrupting activity of IM that has not been studied so far. The study provide new information on the potential consequences of chronic exposure to the residues of tyrosine kinase inhibitors, which remain to be further explored.

Exposure to low doses of inorganic arsenic induces transgenerational changes on behavioral and epigenetic markers in zebrafish (Danio rerio)

The ability of environmental pollutants to alter the epigenome with resultant development of behavioral alterations has received more attention in recent years. These alterations can be transmitted and affect later generations that have not been directly in contact with the contaminant. Arsenic (As) is a neurotoxicant and potent epigenetic disruptor that is widespread in the environment; however, the precise potential of As to produce transgenerational effects is unknown. Our study focused on the possible transgenerational effects on behavior by ancestral exposure to doses relevant to the environment of As, and the epigenetic mechanisms that could be involved. Embryos of F0 (ancestral generation) were directly exposed to 50 or 500 ppb of As for 150 days. F0 adults were raised to produce the F1 generation (intergeneration) and subsequently the F2 generation (transgeneration). We evaluated motor and cognitive behavior, neurodevelopment-related genes, and epigenetic markers on the F0 and F2 generation. As proposed in our hypothesis, ancestral arsenic exposure altered motor activity through the development and increased anxiety-like behaviors which were transmitted to the F2 generation. Additionally, we found a reduction in brain-derived neurotrophic factor expression between the F0 and F2 generation, and an increase in methylation on histone H3K4me3 in the nervous system.

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

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

Ocean acidification promotes broad transcriptomic responses in marine metazoans: a literature survey

For nearly a decade, the metazoan-focused research community has explored the impacts of ocean acidification (OA) on marine animals, noting that changes in ocean chemistry can impact calcification, metabolism, acid-base regulation, stress response and behavior in organisms that hold high ecological and economic value. Because OA interacts with several key physiological processes in marine organisms, transcriptomics has become a widely-used method to characterize whole organism responses on a molecular level as well as inform mechanisms that explain changes in phenotypes observed in response to OA. In the past decade, there has been a notable rise in studies that examine transcriptomic responses to OA in marine metazoans, and here we attempt to summarize key findings across these studies. We find that organisms vary dramatically in their transcriptomic responses to pH although common patterns are often observed, including shifts in acid-base ion regulation, metabolic processes, calcification and stress response mechanisms. We also see a rise in transcriptomic studies examining organismal response to OA in a multi-stressor context, often reporting synergistic effects of OA and temperature. In addition, there is an increase in studies that use transcriptomics to examine the evolutionary potential of organisms to adapt to OA conditions in the future through population and transgenerational experiments. Overall, the literature reveals complex organismal responses to OA, in which some organisms will face more dramatic consequences than others. This will have wide-reaching impacts on ocean communities and ecosystems as a whole.

Genotoxic activity of bisphenol A and its analogues bisphenol S, bisphenol F and bisphenol AF and their mixtures in human hepatocellular carcinoma (HepG2) cells

The use of bisphenol A (BPA) in manufacturing of plastics is being gradually replaced by presumably safer analogues such as bisphenol S (BPS), bisphenol F (BPF) and bisphenol AF (BPAF). Despite their widespread occurrence in the environment, there is a knowledge gap in their toxicological profiles. We investigated cytotoxic/genotoxic effects as well as changes in the expression of selected genes involved in the xenobiotic metabolism, response to oxidative stress and DNA damage upon exposure to BPs and their mixtures in human hepatocellular carcinoma HepG2 cells. BPS and BPF slightly decreased the viability of HepG2 cells, while BPAF was the most cytotoxic compound tested. BPA, BPF and BPAF induced the formation of DNA double strand breaks determined with γH2AX assay, while BPS was inactive (5-20 μg/mL). All four BPs up-regulated the expression of CYP1A1 and UGT1A1, while BPS up-regulated and BPAF down-regulated also the expression of GST1A. Only BPA up-regulated oxidative stress responsive gene GCLC, while BPAF up-regulated the expression of CDKN1A and GADD45a. At concentrations relevant for human exposure (ng/mL range) BPA and its analogues as individual compounds and in mixtures did not exert genotoxic activity, whereas BPA and BPAF as well as the mixtures up-regulated the expressions of CYP1A1 and UGT1A1.

Common responses of fish embryos to metals: an integrated analysis of transcriptomes and methylomes in zebrafish embryos under the stress of copper ions or silver nanoparticles

This study demonstrated the common responses of differentially expressed genes (DEGs) and differentially methylated regions (DMRs) under Cu²⁺ or AgNPs stresses in zebrafish, and verified the correlation of the gene transcription and the methylation status of some common DMGs.

Comparative transcriptomics analysis of the river pufferfish (Takifugu obscurus) by tributyltin exposure: Clues for revealing its toxic injury mechanism

TBT residual in water had become a noticeable ecological problem for aquatic ecosystems. The river pufferfish (Takifugu obscurus) is a kind of an anadromous fish species and widely distributed in the East China Sea and the Yellow Sea. Because of the water contamination, the pufferfish wild resource had a sudden decline in recent years. Therefore, the study on the response of pufferfish to the TBT exposure may contribute to reveal toxic injury mechanism of T. obscurus under TBT exposure. In this study, the transcriptional library of T. obscurus liver and gill was constructed and sequenced by an improved Illumina HiseqX10 high-throughput sequencing platform under different concentrations of TBT acute stress. The blood cell numbers distinctly decreased after TBT exposure, showing the adverse effects of TBT invasion and self-adjusting ability of the pufferfish. The production of reactive oxygen species increased, demonstrating the oxidation resistance of T. obscurus when exposed to TBT. The obtained data were compared with the genome data of Takifugu rubripes and transcriptional resource database. On this basis, gene function annotation, analysis and classification were carried out by bioinformatics method, and differential genes related to toxic injury function were screened out. Meanwhile, new toxic related genes and related signal pathways were sought to provide new theoretical guidance for the pathogenesis of T. obscurus exposed to TBT. This study not only enriched the transcriptome data of T. obscurus under environmental stress, but also provided a new research method for the response mechanism of T. obscurus under the stimulation of environmental factors.

Toxicity and non-harmful effects of the soya isoflavones, genistein and daidzein, in embryos of the zebrafish, Danio rerio

Based on the assumed oestrogenic and apoptotic properties of soya isoflavones (genistein, daidzein), and following the current OECD test-guidelines and principle of 3Rs, we have studied the potential toxicity of phytochemicals on the zebrafish embryos test (ZFET). For this purpose, zebrafish embryos at 2-3 h post-fertilisation (hpf) were exposed to both soya isoflavones (from 1.25 mg/L to 20 mg/L) and assayed until 96 hpf. Lethal and sub-lethal endpoints (mortality, hatching rates and malformations) were estimated in the ZFET, which was expanded to potential gene expression markers, determining the lowest observed effect (and transcriptional) concentrations (LOEC, LOTEC), and the no-observable effect (and transcriptional) concentrations (NOEC, NOTEC). The results revealed that genistein is more toxic (LC50-96 hpf: 4.41 mg/L) than daidzein (over 65.15 mg/L). Both isoflavones up-regulated the oestrogen (esrrb) and death receptors (fas) and cyp1a transcript levels. Most thyroid transcript signals were up-regulated by genistein (except for thyroid peroxidase/tpo), and the hatching enzyme (he1a1) was exclusively up-regulated by daidzein (from 1.25 mg/L onwards). The ZFET proved suitable for assessing toxicant effects of both isoflavones and potential disruptions (i.e. oestrogenic, apoptotic, thyroid, enzymatic) during the embryogenesis and the endotrophic larval period.

Toxicology tailored low density oligonucleotide microarray for the thicklip grey mullets (Chelon labrosus): Biomarker gene transcription profile after caging in a polluted harbour

In aquatic organisms inhabiting polluted waters genes are activated to build an adaptive/compensatory defence against the possible effects of pollutants. Such responses can be used as biomarkers of exposure to chemical compounds, outlining the molecular mechanisms activated under specific pollution scenarios. With the aim of exploiting such approach in environmental health assessment, toxicologically relevant gene fragments were sequenced in the thicklip grey mullet (Chelon labrosus) and a toxicologically tailored low-density (160 genes) oligonucleotide microarray was customised. The tool was validated comparing organ/sex specific gene expression profiles and characterising responses under laboratory exposure to model chemicals. Finally, juvenile mullets were caged in a polluted harbour and hepatic gene expression profiles analysed after 5 and 21 days of deployment. Cages were deployed in the inner (IH) and outer (OH) Pasaia harbour, Bay of Biscay. Mussels (Mytilus galloprovincialis) were also caged as biological matrix for chemical bioaccumulation analysis and stress biomarkers measurements. Slightly higher concentrations of chemicals (metals, tributyltin, PAHs, phthalates) were quantified in IH than in OH, fish bile metabolites also revealing higher availability of PAHs in IH. Lysosome membrane stability in mussels was reduced, indicating stress condition in both sites. The developed microarray discriminated mullets showing distinctive expression profiles depending on site and deployment time. Genes related to immune and hypoxia responses were regulated comparing IH and OH at day 5. Phase I and II biotransformation genes, such as cyp2, cyp3 and ugt, were up-regulated in IH, together with the aryl hydrocarbon receptor 2 (ahr2) and the ahr repressor. Similarly, TBT-binding proteins and genes involved in lipid metabolism (pparγ, cyp7) were up-regulated with deployment time. Even if nowadays higher throughput approaches for gene expression analyses are available, the developed mullet tool constitutes a comprehensive tool to assess molecular responses of mullets exposed to pollutants, although it remains to be explored whether it can be applied to assess pollutant exposure in active pollution monitorings and in environmental health assessment.

Use of 5-azacytidine in a proof-of-concept study to evaluate the impact of pre-natal and post-natal exposures, as well as within generation persistent DNA methylation changes in Daphnia

Short-term exposures at critical stages of development can lead to delayed adverse effects long after the initial stressor has been removed, a concept referred to as developmental origin of adult disease. This indicates that organisms' phenotypes may epigenetically reflect their past exposure history as well as reflecting chemicals currently present in their environment. This concept has significant implications for environmental monitoring. However, there is as yet little or no implementation of epigenetics in environmental risk assessment. In a proof-of-principle study we exposed Daphnia magna to 5-azacytidine, a known DNA de-methylating agent. Exposures covered combinations of prenatal and postnatal exposures as well as different exposure durations and recovery stages. Growth, the transcription of genes and levels of metabolites involved in regulating DNA methylation, and methylation levels of several genes were measured. Our data shows that prenatal exposures caused significant changes in the methylome of target genes, indicating that prenatal stages of Daphnia are also susceptible to same level of change as post-natal stages of Daphnia. While the combination of pre- and postnatal exposures caused the most extreme reduction in DNA methylation compared to the control group. Furthermore, some of the changes in the methylation patterns were persistent even after the initial stressor was removed. Our results suggest that epigenetic biomarkers have the potential to be used as indicators of past chemical exposure history of organisms and provide strong support for implementing changes to the current regimes for chemical risk assessment to mimic realistic environmental scenarios.

The differential expression of Chironomus spp genes as useful tools in the search for pollution biomarkers in freshwater ecosystems

Insects of the Chironomidae family are characterized by a wide ecological diversity in freshwater ecosystems. The larvae have the physiological potential to tolerate environmental stress even when there is a low concentration of oxygen, the presence of toxic substances or when there are changes in temperature and salinity. On the other hand, it is important to consider that at a cellular level, when individual insects are exposed to environmental changes, it induces responses of groups of genes that govern the molecular mechanisms related to such tolerance. In this review, using fourth instar larvae of Chironomus spp. in natural conditions and of Chironomus columbiensis under controlled conditions, we will discuss the genetic expression of a group of genes that respond to detoxification and also the biological functions involved and impacted on by mining stressors. The study of macroinvertebrate bioindicator species and their gene expression as a result of mining activity opens a window on the search for genetic biomarkers that could be used in environmental pollution assessments in freshwater ecosystems.

The Transcriptome of the Zebrafish Embryo After Chemical Exposure: A Meta-Analysis

Numerous studies have been published in the past years investigating the transcriptome of the zebrafish embryo (ZFE) upon being subjected to chemical stress. Aiming at a more mechanistic understanding of the results of such studies, knowledge about commonalities of transcript regulation in response to chemical stress is needed. Thus, our goal in this study was to identify and interpret genes and gene sets constituting a general response to chemical exposure. Therefore, we aggregated and reanalyzed published toxicogenomics data obtained with the ZFE. We found that overlap of differentially transcribed genes in response to chemical stress across independent studies is generally low and the most commonly differentially transcribed genes appear in less than 50% of all treatments across studies. However, effect size analysis revealed several genes showing a common trend of differential expression, among which genes related to calcium homeostasis emerged as key, especially in exposure settings up to 24 h post-fertilization. Additionally, we found that these and other downregulated genes are often linked to anatomical regions developing during the respective exposure period. Genes showing a trend of increased expression were, among others, linked to signaling pathways (e.g., Wnt, Fgf) as well as lysosomal structures and apoptosis. The findings of this study help to increase the understanding of chemical stress responses in the developing zebrafish embryo and provide a starting point to improve experimental designs for this model system. In future, improved time- and concentration-resolved experiments should offer better understanding of stress response patterns and access to mechanistic information.

Epigenetics in teleost fish: From molecular mechanisms to physiological phenotypes

While the field of epigenetics is increasingly recognized to contribute to the emergence of phenotypes in mammalian research models across different developmental and generational timescales, the comparative biology of epigenetics in the large and physiologically diverse vertebrate infraclass of teleost fish remains comparatively understudied. The cypriniform zebrafish and the salmoniform rainbow trout and Atlantic salmon represent two especially important teleost orders, because they offer the unique possibility to comparatively investigate the role of epigenetic regulation in 3R and 4R duplicated genomes. In addition to their sequenced genomes, these teleost species are well-characterized model species for development and physiology, and therefore allow for an investigation of the role of epigenetic modifications in the emergence of physiological phenotypes during an organism's lifespan and in subsequent generations. This review aims firstly to describe the evolution of the repertoire of genes involved in key molecular epigenetic pathways including histone modifications, DNA methylation and microRNAs in zebrafish, rainbow trout, and Atlantic salmon, and secondly, to discuss recent advances in research highlighting a role for molecular epigenetics in shaping physiological phenotypes in these and other teleost models. Finally, by discussing themes and current limitations of the emerging field of teleost epigenetics from both theoretical and technical points of view, we will highlight future research needs and discuss how epigenetics will not only help address basic research questions in comparative teleost physiology, but also inform translational research including aquaculture, aquatic toxicology, and human disease.

Omics in Zebrafish Teratogenesis

The genome revolution represents a complete change on our view of biological systems. The quantitative determination of changes in all major molecular components of the living cells, the "omics" approach, opened whole new fields for all health sciences. Genomics, transcriptomics, proteomics, metabolomics, and others, together with appropriate prediction and modeling tools, will mark the future of developmental toxicity assessment both for wildlife and humans. This is especially true for disciplines, like teratology, which rely on studies in model organisms, as studies at lower levels of organization are difficult to implement. Rodents and frogs have been the favorite models for studying human reproductive and developmental disorders for decades. Recently, the study of the development of zebrafish embryos (ZE) is becoming a major alternative tool to adult animal testing. ZE intrinsic characteristics makes this model a unique system to analyze in vivo developmental alterations that only can be studied applying in toto approaches. Moreover, under actual legislations, ZE is considered as a replacement model (and therefore, excluded from animal welfare regulations) during the first 5 days after fertilization. Here we review the most important components of the zebrafish toolbox available for analyzing early stages of embryotoxic events that could eventually lead to teratogenesis.

Integration of metabolomics and transcriptomics in nanotoxicity studies

Biomedical research involving nanoparticles has produced useful products with medical applications. However, the potential toxicity of nanoparticles in biofluids, cells, tissues, and organisms is a major challenge. The '-omics' analyses provide molecular profiles of multifactorial biological systems instead of focusing on a single molecule. The 'omics' approaches are necessary to evaluate nanotoxicity because classical methods for the detection of nanotoxicity have limited ability in detecting miniscule variations within a cell and do not accurately reflect the actual levels of nanotoxicity. In addition, the 'omics' approaches allow analyses of in-depth changes and compensate for the differences associated with high-throughput technologies between actual nanotoxicity and results from traditional cytotoxic evaluations. However, compared with a single omics approach, integrated omics provides precise and sensitive information by integrating complex biological conditions. Thus, these technologies contribute to extended safety evaluations of nanotoxicity and allow the accurate diagnoses of diseases far earlier than was once possible in the nanotechnology era. Here, we review a novel approach for evaluating nanotoxicity by integrating metabolomics with metabolomic profiling and transcriptomics, which is termed "metabotranscriptomics". [BMB Reports 2018; 51(1): 14-20].

Toxicological analysis of triadimefon on endocrine disruption and oxidative stress during rare minnow (Gobiocypris rarus) larvae development

Triadimefon (TDF) is a systemic wide-spectrum antifungal compound that is widely used in agriculture to inhibit fungal growth on various crops. Since previous studies focused on the embryo and adult life stages in the investigation of ecological impact, here we investigated the long-term effects of TDF (1, 10, 100 μg/L) on rare minnow during its larvae development. TDF caused an anti-estrogenic effect by decreasing vitellogenin (VTG) and CYP19a mRNA level, and inhibiting the aromatase activity and VTG levels after a 3, 6, 10, or 14-day exposure in rare minnow larvae. TDF also disturbed the endocrine disruption by regulating the transcription of estrogen receptors ERα, ERβ1 and ERβ2, CYP1a, CYP11, CYP17, steroidogenic acute regulator (STAR), doublesex and mab-3 related transcription factor (DMRT1), gonadotropin-releasing hormone (GnRH2), GnRH3, GnRHR1A, and GnRHR1B. Furthermore, TDF induced the accumulation of reactive oxygen species (ROS) and the activity of antioxidant proteins glutathione peroxidase (GPX), superoxide dismutase (SOD), and catalase (CAT), and significantly increased the transcriptions of stress response genes P53, growth arrest and DNA damage-inducible 45 alpha (Gadd45α), and COX1, suggested that TDF might cause oxidative stress during larvae development. The changes in transcript and biological levels represented the potential adaptive or compensatory responses to impaired oxidative stress and endocrine system after TDF exposure in rare minnow during its larvae development.

Transcriptomics in toxicology

Xenobiotics, of which many are toxic, may enter the human body through multiple routes. Excessive human exposure to xenobiotics may exceed the body's capacity to defend against the xenobiotic-induced toxicity and result in potentially fatal adverse health effects. Prevention of the adverse health effects, potentially associated with human exposure to the xenobiotics, may be achieved by detecting the toxic effects at an early, reversible and, therefore, preventable stage. Additionally, an understanding of the molecular mechanisms underlying the toxicity may be helpful in preventing and/or managing the ensuing adverse health effects. Human exposures to a large number of xenobiotics are associated with hepatotoxicity or pulmonary toxicity. Global gene expression changes taking place in biological systems, in response to exposure to xenobiotics, may represent the early and mechanistically relevant cellular events contributing to the onset and progression of xenobiotic-induced adverse health outcomes. Hepatotoxicity and pulmonary toxicity resulting from exposure to xenobiotics are discussed as specific examples to demonstrate the potential application of transcriptomics or global gene expression analysis in the prevention of adverse health effects associated with exposure to xenobiotics.

Integrated assessment of toxic effects of maghemite (γ-Fe2O3) nanoparticles in zebrafish

The increasing use of nanotechnology in the last decade has raised concerns about the impact of nanoparticles in the environment. In particular, the potential harmful effects of iron oxide nanoparticles (IONPs) in aquatic organisms have been poorly addressed. We analyze here the toxic effects induced by IONPs in zebrafish using a combination of classical (genotoxicity, oxidative stress) and molecular (transcriptomic) methodologies. Adult animals were exposed for 96h to five sub-lethal IONP concentrations, ranging from 4.7 to 74.4mg/L. Comet and micronucleus assays revealed a significant number of DNA lesions induced by IONPs at all concentrations tested. Conversely, the thiobarbituric acid reactive substances (TBARS) test detected only a mild oxidative damage in liver cells (∼1.5-fold increase of malondialdehyde concentrations) and only at the two higher IONP concentrations tested. Microarray analysis of liver samples identified 953 transcripts (927 unique genes) differentially expressed between controls and IONP-exposed samples. Subsequent functional analysis identified genes related to cation/metal ion binding, membrane formation, and morphogenesis among the transcripts overrepresented upon IONP treatments, whereas mRNAs encompassing genes associated with RNA biogenesis, translation, ribosomes, and several metabolic processes became underrepresented in treated samples. Taken together, these results indicate considerable genotoxic effects of IONPs combined with general negative effect on cell growth and on the ability of the cell produce new proteins. On the contrary, IONPs showed only a limited capacity to induce oxidative stress. To our knowledge, this is the first study on IONPs toxicity using such an integrative approach in an aquatic organism.

RNA-sequencing and pathway analysis reveal alteration of hepatic steroid biosynthesis and retinol metabolism by tributyltin exposure in male rare minnow (Gobiocypris rarus)

Tributyltin (TBT) is widely spread in aquatic ecosystems. Although adverse effects of TBT on reproduction and lipogenesis are observed in fishes, the underlying mechanisms, especially in livers, are still scarce and inconclusive. Thus, RNA-sequencing runs were performed on the hepatic libraries of adult male rare minnow (Gobiocypris rarus) after TBT exposure for 60d. After differentially expressed genes were identified, enrichment analysis and validation by quantitative real-time PCR were conducted. The results showed that TBT up-regulated the profile of hepatic genes in the steroid biosynthesis pathway and down-regulated the profile of hepatic genes in the retinol metabolism pathway. In the hepatic steroid biosynthesis pathway, TBT might induce biosynthesis of cholesterol, which could affect the bioavailability of steroid hormones. More important, 3beta-hydroxysteroid 3-dehydrogenase, a key enzyme in the biosynthesis of all active steroid hormones, was up-regulated by TBT exposure. In the hepatic retinol metabolism pathway, TBT impaired retinoic acid homeostasis which plays essential roles in both reproduction and lipogenesis. The results of two pathways offered new mechanisms underlying the toxicology of TBT and represented a starting point from which detailed mechanistic links should be explored.

A multi-omic approach to elucidate low-dose effects of xenobiotics in zebrafish (Danio rerio) larvae

Regulatory-approved toxicity assays such as the OECD Fish Embryo Toxicity Assay (TG236) allow correlation of chemical exposure to adverse morphological phenotypes. However, these assays are ineffective in assessing sub-lethal (i.e. low-dose) effects, or differentiating between similar phenotypes induced by different chemicals. Inclusion of multi-omic analyses in studies investigating xenobiotic action provides improved characterization of biological response, thereby enhancing prediction of toxicological outcomes in whole animals in the absence of morphological effects. In the current study, we assessed perturbations in both the metabolome and transcriptome of zebrafish (Danio rerio; ZF) larvae exposed from 96 to 120h post fertilization to environmental concentrations of acetaminophen (APAP), diphenhydramine (DH), carbamazepine (CBZ), and fluoxetine (FLX); common pharmaceuticals with known mechanisms of action. Multi-omic responses were evaluated independently and integrated to identify molecular interactions and biological relevance of the responses. Results indicated chemical- and dose-specific changes suggesting differences in the time scale of transcript abundance and metabolite production. Increased impact on the metabolome relative to the transcriptome in FLX-treated animals suggests a stronger post-translational effect of the treatment. In contrast, the transcriptome showed higher sensitivity to perturbation in DH-exposed animals. Integration of 'omic' responses using multivariate approaches provided additional insights not obtained by independent 'omic' analyses and demonstrated that the most distinct overall response profiles were induced following low-dose exposure for all 4 pharmaceuticals. Importantly, changes in transcript abundance corroborated with predictions from metabolomic enrichment analyses and the identified perturbed biological pathways aligned with known xenobiotic mechanisms of action. This work demonstrates that a multi-omic toxicological approach, coupled with a sensitive animal model such as ZF larvae, can help characterize the toxicological relevance of acute low-dose chemical exposures.

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.

Toxicogenomics to Evaluate Endocrine Disrupting Effects of Environmental Chemicals Using the Zebrafish Model

The extent of our knowledge on the number of chemical compounds related to anthropogenic activities that can cause damage to the environment and to organisms is increasing. Endocrine disrupting chemicals (EDCs) are one group of potentially hazardous substances that include natural and synthetic chemicals and have the ability to mimic endogenous hormones, interfering with their biosynthesis, metabolism, and normal functions. Adverse effects associated with EDC exposure have been documented in aquatic biota and there is widespread interest in the characterization and understanding of their modes of action. Fish are considered one of the primary risk organisms for EDCs. Zebrafish (Danio rerio) are increasingly used as an animal model to study the effects of endocrine disruptors, due to their advantages compared to other model organisms. One approach to assess the toxicity of a compound is to identify those patterns of gene expression found in a tissue or organ exposed to particular classes of chemicals, through new technologies in genomics (toxicogenomics), such as microarrays or whole-genome sequencing. Application of these technologies permit the quantitative analysis of thousands of gene expression changes simultaneously in a single experiment and offer the opportunity to use transcript profiling as a tool to predict toxic outcomes of exposure to particular compounds. The application of toxicogenomic tools for identification of chemicals with endocrine disrupting capacity using the zebrafish model system is reviewed.

DNA methyltransferases and stress-related genes expression in zebrafish larvae after exposure to heat and copper during reprogramming of DNA methylation

DNA methylation, a well-studied epigenetic mark, is important for gene regulation in adulthood and for development. Using genetic and epigenetic approaches, the present study aimed at evaluating the effects of heat stress and copper exposure during zebrafish early embryogenesis when patterns of DNA methylation are being established, a process called reprogramming. Embryos were exposed to 325 μg Cu/L from fertilization (<1 h post fertilization - hpf) to 4 hpf at either 26.5 °C or 34 °C, followed by incubation in clean water at 26.5 °C till 96 hpf. Significant increased mortality rates and delayed hatching were observed following exposure to combined high temperature and Cu. Secondly, both stressors, alone or in combination, significantly upregulated the expression of de novo DNA methyltransferase genes (dnmt3) along with no differences in global cytosine methylation level. Finally, Cu exposure significantly increased the expression of metallothionein (mt2) and heat shock protein (hsp70), the latter being also increased following exposure to high temperature. These results highlighted the sensitivity of early embryogenesis and more precisely of the reprogramming period to environmental challenges, in a realistic situation of combined stressors.

Advancing toxicology research using in vivo high throughput toxicology with small fish models

Small freshwater fish models, especially zebrafish, offer advantages over traditional rodent models, including low maintenance and husbandry costs, high fecundity, genetic diversity, physiology similar to that of traditional biomedical models, and reduced animal welfare concerns. The Collaborative Workshop on Aquatic Models and 21st Century Toxicology was held at North Carolina State University on May 5-6, 2014, in Raleigh, North Carolina, USA. Participants discussed the ways in which small fish are being used as models to screen toxicants and understand mechanisms of toxicity. Workshop participants agreed that the lack of standardized protocols is an impediment to broader acceptance of these models, whereas development of standardized protocols, validation, and subsequent regulatory acceptance would facilitate greater usage. Given the advantages and increasing application of small fish models, there was widespread interest in follow-up workshops to review and discuss developments in their use. In this article, we summarize the recommendations formulated by workshop participants to enhance the utility of small fish species in toxicology studies, as well as many of the advances in the field of toxicology that resulted from using small fish species, including advances in developmental toxicology, cardiovascular toxicology, neurotoxicology, and immunotoxicology. We alsoreview many emerging issues that will benefit from using small fish species, especially zebrafish, and new technologies that will enable using these organisms to yield results unprecedented in their information content to better understand how toxicants affect development and health.

Environmentally induced epigenetic toxicity: potential public health concerns

Throughout our lives, epigenetic processes shape our development and enable us to adapt to a constantly changing environment. Identifying and understanding environmentally induced epigenetic change(s) that may lead to adverse outcomes is vital for protecting public health. This review, therefore, examines the present understanding of epigenetic mechanisms involved in the mammalian life cycle, evaluates the current evidence for environmentally induced epigenetic toxicity in human cohorts and rodent models and highlights the research considerations and implications of this emerging knowledge for public health and regulatory toxicology. Many hundreds of studies have investigated such toxicity, yet relatively few have demonstrated a mechanistic association among specific environmental exposures, epigenetic changes and adverse health outcomes in human epidemiological cohorts and/or rodent models. While this small body of evidence is largely composed of exploratory in vivo high-dose range studies, it does set a precedent for the existence of environmentally induced epigenetic toxicity. Consequently, there is worldwide recognition of this phenomenon, and discussion on how to both guide further scientific research towards a greater mechanistic understanding of environmentally induced epigenetic toxicity in humans, and translate relevant research outcomes into appropriate regulatory policies for effective public health protection.

Developmental profiles and expression of the DNA methyltransferase genes in the fathead minnow (Pimephales promelas) following exposure to di-2-ethylhexyl phthalate

DNA methylation is an epigenetic regulator of gene expression, and this process has been shown to be disrupted by environmental contaminants. Di-2-(ethylhexyl) phthalate (DEHP) and related phthalate esters have been shown to affect development in early life stages of fish and can alter genomic methylation patterns in vertebrates. The objectives of this study were the following: (1) Describe the expression patterns of the DNA methyltransferase (dnmt) genes during early fathead minnow (FHM) development. These genes are critical for methylation and imprinting during development. (2) Determine the effects of DEHP on the development of FHM larvae [1 and 14 days post-hatch (dph)]. (3) Determine the effect of DEHP on dnmt expression and global methylation status in larval FHM. FHMs were first collected over a developmental time course [1, 3, 5, 6, and 14 days post-fertilization (dpf)] to investigate the expression patterns of five dnmt isoforms. The expression of dnmt1 and dnmt7 was relatively high in embryos at 1 dpf but was variable in expression, and these transcripts were later expressed at a lower level (>3 dpf); dnmt3 was significantly higher in embryos at 1 dpf compared to those at 3 dpf. Dnmt6 showed more of a constitutive pattern of expression during the first 2 weeks of development, and the mRNA levels of dnmt8 were higher in embryos at 5 and 6 dpf compared to those at 1 and 3 dpf, corresponding to the hatching period of the embryos. A waterborne exposure to three concentrations of DEHP (1, 10 and 100 µg/L) was conducted on 1-day FHM embryos for 24 h and on larval fish for 2 weeks, ending at 14 dpf. DEHP did not negatively affect survival, hatch rate, or the expression of dnmt isoforms in FHMs. There were no differences in global cytosine methylation following DEHP treatments in 14 dpf larvae, suggesting that environmentally relevant levels of DEHP may not affect global methylation at this stage of FHM development. However, additional targeted methylome studies are required to determine whether specific gene promoters are differently methylated following exposure to DEHP. This study offers new insight into the roles of the dnmt enzymes during FHM development.

Modulation of DNA methylation machineries in Japanese rice fish (Oryzias latipes) embryogenesis by ethanol and 5-azacytidine

As a sequel of our investigations on the impact of epigenome in inducing fetal alcohol spectrum disorder (FASD) phenotypes in Japanese rice fish, we have investigated on several DNA methylation machinery genes including DNA methyl transferase 3ba (dnmt3ba) and methyl binding proteins (MBPs), namely, mbd1b, mbd3a, mbd3b, and mecp2 at the transcription level. Studies were made during normal development, from 0day post fertilization (dpf) to hatching, and also exposing the fertilized eggs to ethanol or a DNMT inhibitor, 5-azacytidine (5-azaC). We observed that during development, all these genes followed distinct expression patterns, generally high mRNA copies in early phases (0-1dpf) and significantly low mRNA copies prior to or after hatching. Ethanol (100-500mM, 0-2dpf) was unable to alter any of these mRNAs in 2dpf; additional four day (2-6dpf) maintenance of these embryos in ethanol-free environment, on 6dpf, was also unable to establish any significant difference in these mRNA levels in comparison with the corresponding controls. However, continuous exposure of fertilized eggs in 300mM ethanol, 0-6dpf, showed significantly high mRNA copies only in MBPs (mbd1b, mbd3a, mbd3b, mecp2). 5-azaC (2mM) on 2dpf was able to enhance only mbd3b mRNA. Removal of 5-azaC and maintenance of these embryos in clean medium, 2-6dpf, showed significantly enhanced mbd3b and mecp2 mRNAs compared to corresponding controls on 6dpf. Our studies showed that in Japanese rice fish embryogenesis both ethanol and 5-azaC have the potential to specifically modulate the developmental rhythm of DNA methylation machineries.

Functional genomics to assess biological responses to marine pollution at physiological and evolutionary timescales: toward a vision of predictive ecotoxicology

Marine pollution is ubiquitous, and is one of the key factors influencing contemporary marine biodiversity worldwide. To protect marine biodiversity, how do we surveil, document and predict the short- and long-term impacts of pollutants on at-risk species? Modern genomics tools offer high-throughput, information-rich and increasingly cost-effective approaches for characterizing biological responses to environmental stress, and are important tools within an increasing sophisticated kit for surveiling and assessing impacts of pollutants on marine species. Through the lens of recent research in marine killifish, we illustrate how genomics tools may be useful for screening chemicals and pollutants for biological activity and to reveal specific mechanisms of action. The high dimensionality of transcriptomic responses enables their usage as highly specific fingerprints of exposure, and these fingerprints can be used to diagnose environmental problems. We also emphasize that molecular pathways recruited to respond at physiological timescales are the same pathways that may be targets for natural selection during chronic exposure to pollutants. Gene complement and sequence variation in those pathways can be related to variation in sensitivity to environmental pollutants within and among species. Furthermore, allelic variation associated with evolved tolerance in those pathways could be tracked to estimate the pace of environmental health decline and recovery. We finish by integrating these paradigms into a vision of how genomics approaches could anchor a modernized framework for advancing the predictive capacity of environmental and ecotoxicological science.

Zebrafish as a model to study the role of DNA methylation in environmental toxicology

Environmental epigenetics is a rapidly growing field which studies the effects of environmental factors such as nutrition, stress, and exposure to compounds on epigenetic gene regulation. Recent studies have shown that exposure to toxicants in vertebrates is associated with changes in DNA methylation, a major epigenetic mechanism affecting gene transcription. Zebra fish, a well-known model in toxicology and developmental biology, are emerging as a model species in environmental epigenetics despite their evolutionary distance to rodents and humans. In this review, recent insights in DNA methylation during zebra fish development are discussed and compared to mammalian models in order to evaluate zebra fish as a model to study the role of DNA methylation in environmental toxicology. Differences exist in DNA methylation reprogramming during early development, whereas in later developmental stages, tissue distribution of both 5-methylcytosine and 5-hydroxymethylcytosine seems more conserved between species, as well as basic DNA (de)methylation mechanisms. All DNA methyl transferases identified so far in mammals are present in zebra fish, as well as a number of major demethylation pathways. However, zebra fish appear to lack some methylation pathways present in mammals, such as parental imprinting. Several studies report effects on DNA methylation in zebra fish following exposure to environmental contaminants, such as arsenic, benzo[a]pyrene, and tris(1,3-dichloro-2-propyl)phosphate. Though more research is needed to examine heritable effects of contaminant exposure on DNA methylation, recent data suggests the usefulness of the zebra fish as a model in environmental epigenetics.

Environmental epigenetics: A promising venue for developing next-generation pollution biomonitoring tools in marine invertebrates

Environmental epigenetics investigates the cause-effect relationships between specific environmental factors and the subsequent epigenetic modifications triggering adaptive responses in the cell. Given the dynamic and potentially reversible nature of the different types of epigenetic marks, environmental epigenetics constitutes a promising venue for developing fast and sensible biomonitoring programs. Indeed, several epigenetic biomarkers have been successfully developed and applied in traditional model organisms (e.g., human and mouse). Nevertheless, the lack of epigenetic knowledge in other ecologically and environmentally relevant organisms has hampered the application of these tools in a broader range of ecosystems, most notably in the marine environment. Fortunately, that scenario is now changing thanks to the growing availability of complete reference genome sequences along with the development of high-throughput DNA sequencing and bioinformatic methods. Altogether, these resources make the epigenetic study of marine organisms (and more specifically marine invertebrates) a reality. By building on this knowledge, the present work provides a timely perspective highlighting the extraordinary potential of environmental epigenetic analyses as a promising source of rapid and sensible tools for pollution biomonitoring, using marine invertebrates as sentinel organisms. This strategy represents an innovative, groundbreaking approach, improving the conservation and management of natural resources in the oceans.