Detrimental impacts and QSAR baseline toxicity assessment of Japanese medaka embryos exposed to methylparaben and its halogenated byproducts.

Despite the theoretical risk of forming halogenated methylparabens (halo-MePs) during water chlorination in the absence or presence of bromide ions, there remains a lack of in vivo toxicological assessments on vertebrate organisms for halo-MePs. This research addresses these gaps by investigating the lethal (assessed by embryo coagulation) or sub-lethal (assessed by hatching success/heartbeat rate) toxicity and teratogenicity (assessed by deformity rate) of MeP and its mono- and di-halogen derivatives (Cl- or Br-) using Japanese medaka embryos. In assessing selected apical endpoints to discern patterns in physiological or biochemical alterations, heightened toxic impacts were observed for halo-MePs compared to MeP. These include a higher incidence of embryo coagulation (4-36 fold), heartbeat rate decrement (11-36 fold), deformity rate increment (32-223 fold), hatching success decrement (11-59 fold), and an increase in Reactive Oxygen Species (ROS) level (1.2-7.4 fold)/Catalase (CAT) activity (1.7-2.8 fold). Experimentally determined LC50 values are correlated and predicted using a Quantitative Structure Activity Relationship (QSAR) based on the speciation-corrected liposome-water distribution ratio (Dlipw, pH 7.5). The QSAR baseline toxicity aligns well with (sub)lethal toxicity and teratogenicity, as evidenced by toxic ratio (TR) analysis showing TR < 10 for MeP exposure in all cases, while significant specific or reactive toxicity was found for halo-MeP exposure, with TR > 10 observed (excepting three values). Our extensive findings contribute novel insights into the intricate interplay of embryonic toxicity during the early-life-stage of Japanese medaka, with a specific focus on highlighting the potential hazards associated with halo-MePs compared to the parent compound MeP.

Vasa identifies germ cells in embryos and gonads of Oryzias celebensis.

Vasa is the most studied germ cell marker that is indispensable for germ cell development in teleost fishes. Here, a vasa full-length cDNA from Oryzias celebensis was isolated. Analysis of gene expression by reversed transcription polymerase chain reaction and in situ hybridization showed the vasa transcript was maternally inherited and specifically expressed in germ cells during embryogenesis and in adult gonads. During embryogenesis, vasa mRNA was widely distributed in the embryos until the somitogenesis stage and then specifically expressed in primordial germ cells (PGCs). In the testis, vasa expression was highest in spermatogonia and gradually decreased during spermatogenesis. In ovary, vasa expression was present predominantly in immature oocytes and persisted throughout oogenesis. Constructs containing green or red fluorescence proteins and vasa 3' UTR or dnd 3' UTR, confirmed stable vasa expression in the PGCs of O. celebensis and co-expression of the two genes. In summary, the conservation of vasa expression in embryonic and adult germ cells of both sexes compared to other vertebrates suggests its function is also widely conserved.

Effect of benzotriazole on oxidative stress response and transcriptional gene expression in Oryzias latipes and Danio rerio embryo.

Emerging contaminants (EC) such as benzotriazole are being released into the environment in various ways, therefore it is necessary to understand how organisms are affected by EC. In this study, we exposed medaka (Oryzias latipes) and zebrafish (Danio rerio) during their embryonic period (1 day after hatching) to benzotriazole to investigate its effects on oxidative stress (ROS, GSH, GST, SOD, CAT and MDA) and changes in gene expression patterns. In both medaka and zebrafish, the influence of oxidative stress was confirmed through an increased MDA level and changes in the ROS and GSH levels. Antioxidant enzymes such as GST, CAT, and SOD were affected by benzotriazole; however, medaka and zebrafish showed different patterns in the effects by benzotriazole. Results of oxidative stress genes expression showed that medaka had either no influence or had a decrease in the gene expression profile, whereas zebrafish had a statistically significant increase in the expression of some genes. The cyp1a gene expression was increased in both species. However, vtg gene expression was increased only in zebrafish but decreased in medaka, indicating no estrogenic effects in medaka. Apoptosis genes showed changes in expression in both the species but was these changes were not dose-dependent. However, zebrafish caspase-9 gene expression was increased in all of the exposed groups, suggesting the effects on the intrinsic pathway associated with caspase-9. In conclusion, the results indicate that the toxic effects of benzotriazole differ at various levels in the two small fish medaka and zebrafish embryos.

Local tissue interactions govern pLL patterning in medaka.

Vertebrate organs are arranged in a stereotypic, species-specific position along the animal body plan. Substantial morphological variation exists between related species, especially so in the vastly diversified teleost clade. It is still unclear how tissues, organs and systems can accommodate such diverse scaffolds. Here, we use the distinctive arrangement of neuromasts in the posterior lateral line (pLL) system of medaka fish to address the tissue-interactions defining a pattern. We show that patterning in this peripheral nervous system is established by autonomous organ precursors independent of neuronal wiring. In addition, we target the keratin 15 gene to generate stuck-in-the-midline (siml) mutants, which display epithelial lesions and a disrupted pLL patterning. By using siml/wt chimeras, we determine that the aberrant siml pLL pattern depends on the mutant epithelium, since a wild type epithelium can rescue the siml phenotype. Inducing epithelial lesions by 2-photon laser ablation during pLL morphogenesis phenocopies siml genetic mutants and reveals that epithelial integrity defines the final position of the embryonic pLL neuromasts. Our results using the medaka pLL disentangle intrinsic from extrinsic properties during the establishment of a sensory system. We speculate that intrinsic programs guarantee proper organ morphogenesis, while instructive interactions from surrounding tissues facilitates the accommodation of sensory organs to the diverse body plans found among teleosts.

The Embryonic Key Pluripotent Factor NANOG Mediates Glioblastoma Cell Migration via the SDF1/CXCR4 Pathway.

NANOG is a key transcription factor required for maintaining pluripotency of embryonic stem cells. Elevated NANOG expression levels have been reported in many types of human cancers, including lung, oral, prostate, stomach, breast, and brain. Several studies reported the correlation between NANOG expression and tumor metastasis, revealing itself as a powerful biomarker of poor prognosis. However, how NANOG regulates tumor progression is still not known. We previously showed in medaka fish that Nanog regulates primordial germ cell migration through Cxcr4b, a chemokine receptor known for its ability to promote migration and metastasis in human cancers. Therefore, we investigated the role of human NANOG in CXCR4-mediated cancer cell migration. Of note, we found that NANOG regulatory elements in the CXCR4 promoter are functionally conserved in medaka fish and humans, suggesting an evolutionary conserved regulatory axis. Moreover, CXCR4 expression requires NANOG in human glioblastoma cells. In addition, transwell assays demonstrated that NANOG regulates cancer cell migration through the SDF1/CXCR4 pathway. Altogether, our results uncover NANOG-CXCR4 as a novel pathway controlling cellular migration and support Nanog as a potential therapeutic target in the treatment of Nanog-dependent tumor progression.

Exposure of silver nanocolloids causes glycosylation disorders and embryonic deformities in medaka.

Silver nanomaterials such as silver nanocolloids (SNC) contribute to environmental pollution and have adverse ecological effects on aquatic organisms. In particular, chemical exposure of fish during embryogenesis leads to deformities and puts the population at risk. Although glycans and glycosylation are known to be important for proper morphology in embryogenesis, little glycobiology-based research has examined morphological disorders caused by environmental pollutants. This study addressed the glycobiological effects of SNC exposure on medaka embryogenesis. After exposure of medaka embryos to SNC, deformities such as small heads and deformed eyes were observed. The expression of five glycan-related genes (alg2, gnsb, b4galt2, b3gat1a, and b3gat2) was significantly altered, with changes depending on the embryonic stage at exposure, with more severe deformities with exposure at earlier stages. In situ hybridization analyses indicated that the five genes were expressed mainly in the head region; exposure of SNC suppressed alg2 and gnsb and enhanced b4galt2 and b3gat1a expression relative to controls on day 7. Loss (siRNA)- and gain (RNA overexpression)-of-function experiments confirmed that alg2, gnsb, and b4galt2 are essential for embryogenesis. The effects of SNC exposure on glycan synthesis were estimated by glycan structure analysis. In the medaka embryo, high mannose-type glycans were dominant, and SNC exposure altered glycan synthesis. The alteration was more significant when exposure occurred at an early stage of medaka embryogenesis. Thus, SNC exposure causes embryonic deformities in medaka embryos through disordered glycosylation.

Pre-validation of choriogenin H transgenic medaka eleutheroembryos as a quantitative estrogenic activity test method.

The choriogenin H - EGFP transgenic medaka (Oryzias melastigma) has been used to test estrogenic substances and quantify estrogenic activity into 17ß-estradiol (E2) equivalency (EEQ). The method uses 8 eleutheroembryos in 2 ml solution per well and 3 wells per treatment in 24-well plates at 26 ± 1 °C for 24 ± 2 h, with subsequent measurements of induced GFP signal intensity. EEQ measurements are calculated using a E2 probit regression model with a coefficient of determination (R2) > 0.90. The selectivity was confirmed evaluating 27 known estrogenic and 5 known non-estrogenic compounds. Limit of quantitation (LOQ), recovery rate and bias were calculated to be 1 ng/ml EEQ, 104% and 4% respectively. Robustness analysis revealed exposure temperature is a sensitive parameter that should be kept at 26 ± 1 °C. The repeatability of intra- and inter-laboratories achieved CV < 30% for most tested food and cosmetics samples. The lot-lot stability was confirmed by the stable EEQ qualitative control (QC, 1 ng/mL E2) and calibration curve results. The stability of standard reagents, samples and sample extracts was also investigated. These data demonstrated this method to be an accurate indicator of estrogenic activity for both chemicals and extracts.

Igf signaling couples retina growth with body growth by modulating progenitor cell division.

How the body and organs balance their relative growth is of key importance for coordinating size and function. This is of particular relevance in organisms, which continue to grow over their entire life span. We addressed this issue in the neuroretina of medaka fish (Oryzias latipes), a well-studied system with which to address vertebrate organ growth. We reveal that a central growth regulator, Igf1 receptor (Igf1r), is necessary and sufficient for proliferation control in the postembryonic retinal stem cell niche: the ciliary marginal zone (CMZ). Targeted activation of Igf1r signaling in the CMZ uncouples neuroretina growth from body size control, and we demonstrate that Igf1r operates on progenitor cells, stimulating their proliferation. Activation of Igf1r signaling increases retinal size while preserving its structural integrity, revealing a modular organization in which progenitor differentiation and neurogenesis are self-organized and highly regulated. Our findings position Igf signaling as a key module for controlling retinal size and composition, with important evolutionary implications.

Miniaturization: How many cells are needed to build a tooth?

BACKGROUND: Organs that develop early in life, and are replaced by a larger version as the animal grows, often represent a miniature version of the adult organ. Teeth constituting the first functional dentition in small-sized teleost fish, such as medaka (Oryzias latipes), are examples of such miniature organs. With a dentin cone as small as the size of one human cell, or even smaller, these teeth raise the question how many dentin-producing cells (odontoblasts) are required to build such a tooth, and whether this number can be as little as one. RESULTS: Based on detailed observations with transmission electron microscopy (TEM) and TEM-based 3D-reconstructions, we show that only one mesenchymal cell qualifies as a true odontoblast. A second mesenchymal cell potentially participates in dentin formation, but only at a late stage of tooth development. Moreover, the fate of these cells appears to be specified very early during tooth development. CONCLUSIONS: Our observations indicate that in this system, one single odontoblast fulfills roles normally exerted by a large and communicating cell population. First-generation teeth in medaka thus provide an exciting model to study integration of multiple functions into a single cell.

TMT-Opsins differentially modulate medaka brain function in a context-dependent manner.

Vertebrate behavior is strongly influenced by light. Light receptors, encoded by functional opsin proteins, are present inside the vertebrate brain and peripheral tissues. This expression feature is present from fishes to human and appears to be particularly prominent in diurnal vertebrates. Despite their conserved widespread occurrence, the nonvisual functions of opsins are still largely enigmatic. This is even more apparent when considering the high number of opsins. Teleosts possess around 40 opsin genes, present from young developmental stages to adulthood. Many of these opsins have been shown to function as light receptors. This raises the question of whether this large number might mainly reflect functional redundancy or rather maximally enables teleosts to optimally use the complex light information present under water. We focus on tmt-opsin1b and tmt-opsin2, c-opsins with ancestral-type sequence features, conserved across several vertebrate phyla, expressed with partly similar expression in non-rod, non-cone, non-retinal-ganglion-cell brain tissues and with a similar spectral sensitivity. The characterization of the single mutants revealed age- and light-dependent behavioral changes, as well as an impact on the levels of the preprohormone sst1b and the voltage-gated sodium channel subunit scn12aa. The amount of daytime rest is affected independently of the eyes, pineal organ, and circadian clock in tmt-opsin1b mutants. We further focused on daytime behavior and the molecular changes in tmt-opsin1b/2 double mutants, and found that-despite their similar expression and spectral features-these opsins interact in part nonadditively. Specifically, double mutants complement molecular and behavioral phenotypes observed in single mutants in a partly age-dependent fashion. Our work provides a starting point to disentangle the highly complex interactions of vertebrate nonvisual opsins, suggesting that tmt-opsin-expressing cells together with other visual and nonvisual opsins provide detailed light information to the organism for behavioral fine-tuning. This work also provides a stepping stone to unravel how vertebrate species with conserved opsins, but living in different ecological niches, respond to similar light cues and how human-generated artificial light might impact on behavioral processes in natural environments.

Both of marine fish species Oryzias melastigma and Pagrus major all failing in early localization at embryo stage by vasa RNA.

Germ cells are essential for gonadal development. As precursors of germ cells, primordial germ cells (PGCs) are particularly important for germline formation. However, the research on distribution patterns of PGCs in marine fish is very limited, especially for economic species. The vasa gene has been widely used as marker to identify PGCs origination and migration because of vasa RNA is a component of germ plasm. In this study, we isolated full-length vasa cDNA (Omvas and Pmvas) from marine medaka (Oryzias melastigma) and red seabream (Pagrus major), detected vasa transcripts in different tissues by RT-PCR and described vasa expression patterns during embryogenesis and gametogenesis by in situ hybridization. At the same time, we also explored the relationship between early distribution of germ plasm components and species evolution. The results demonstrated that deduced amino acid sequence of Omvas and Pmvas shared several conserved motifs of Vasa homologues and high identity with other teleost, and vasa transcripts were exclusively detected in early germ cells of gonad. During embryogenesis, vasa RNA of both fishes, like medaka (Oryzias latipes), failed to localize at cleavage furrows and distributed uniformly throughout each blastomere. This study firstly discovered that the marine economic fish, red seabream, lost vasa RNA early specific localization at cleavage furrows and distinctive distribution in germ cells. In addition, compared with other teleost, we found that early distribution of germ plasm might not relate to species evolution. This will improve our understanding of vasa localization modes in teleost, and facilitate fish germ cell manipulation.

Silver nanocolloid affects hindbrain vascular formation during medaka embryogenesis.

Angiogenesis is essential for the normal development of an embryo. Silver nanocolloid (SNC) is known to induce vascular malformation in the medaka embryo. We focused on the development of the central arteries (CtAs) in the hindbrain of Japanese medaka. The CtAs and the basilar artery from which they branch are essential for transporting the blood and nutrients necessary to support the hindbrain parenchyma and the development of the pons and cerebellum from the hindbrain. We exposed medaka embryos at developmental stage 21 (6 somite stage), to 0, 0.5, 5, or 10 mg/L SNC and evaluated hatching rate, number of thrombi per embryo, head size (length and width), body length, and angiogenesis. Although all SNC-exposed embryos hatched, their head size and body length were small in comparison to controls; in addition, the number of thrombi in the head increased and head size and body length decreased as the SNC concentration increased. To evaluate vasculogenic abnormalities, we performed whole-mount in situ hybridization using a vascular marker (eg, fl7) and visualized the CtAs in medaka embryos. In control embryos, CtAs started to sprout at stage 32 (somite completion stage) and their extension was complete by stage 35 (pectoral fin blood circulation stage). In contrast, CtAs failed to sprout in SNC-exposed embryos, and thrombi were present. Furthermore, qRT-PCR analysis showed that SNC significantly suppressed the egfl7 expression level at stage 35. Together, our findings suggest that SNC induced decreased developments of head and body in medaka embryos due to insufficient angiogenesis and hindbrain vascular formation.

Genetic developmental timing revealed by inter-species transplantations in fish.

The path from a fertilised egg to an embryo involves the coordinated formation of cell types, tissues and organs. Developmental modules comprise discrete units specified by self-sufficient genetic programs that can interact with each other during embryogenesis. Here, we have taken advantage of the different span of embryonic development between two distantly related teleosts, zebrafish (Danio rerio) and medaka (Oryzias latipes) (3 and 9 days, respectively), to explore modularity principles. We report that inter-species blastula transplantations result in the ectopic formation of a retina formed by donor cells - a module. We show that the time taken for the retina to develop follows a genetic program: an ectopic zebrafish retina in medaka develops with zebrafish dynamics. Heterologous transplantation results in a temporal decoupling between the donor retina and host organism, illustrated by two paradigms that require retina-host interactions: lens recruitment and retino-tectal projections. Our results uncover a new experimental system for addressing temporal decoupling along embryonic development, and highlight the presence of largely autonomous but interconnected developmental modules that orchestrate organogenesis.

Medaka (Oryzias latipes) Embryo as a Model for the Screening of Compounds That Counteract the Damage Induced by Ultraviolet and High-Energy Visible Light.

Continuous overexposure to sunlight increases its harmful effects on the skin. For this reason, there is a growing need to characterize economic models more representative of the negative effects and counteracting responses that irradiation causes on human skin. These models will serve for the screening of protective compounds against damage caused by ultraviolet (UV) and high energy visible light (HEV). Therefore, two common in vitro models employed for sunlight irradiation studies, namely human keratinocyte HaCat culture and reconstructed human epidermis (RHE), were compared with the medaka fish embryo model, traditionally used in other scientific disciplines. Using suberythemal doses of UVA and HEV to determine the level of Reactive Oxygen Species (ROS) generation and thymine dimers formed by UVB, we show that medaka embryo responds with a lower damage level, more comparable to human skin, than the other two models, probably due to the protective mechanisms that work in a complete organism. In the same way, the protective effects of antioxidant compounds have the greatest effect on medaka embryos. Taken together, these findings suggest that medaka embryos would be a good alternative in vitro model for sunlight effect studies, and for the screening of molecules with counteracting capacity against the damage caused by UV and HEV.

Distinct expression patterns of seven crucial microRNAs during early embryonic development in medaka (Oryzias latipes).

MicroRNAs (i.e. miRNAs) are small non-coding RNAs that play essential modulation roles in embryonic development in vertebrates. Paternal and maternal miRNAs contribute to the development of post-fertilization embryo and zygotic genome activation. The pattern of expression and their roles in embryonic development of medaka are not clearly understood. The present study, therefore, examined a temporal expression of seven miRNAs, ola-let-7a, ola-miR-202-3p, ola-miR-126-3p, ola-miR-122, ola-miR-92a, ola-miR-125a-3p and ola-miR-430a in sperm, oocytes, and embryos during early developmental stages. Three unique expression patterns of miRNAs were observed. ola-let7a, ola-miR-202-3p and ola-miR-126-3p showed both paternal and maternal expression, and ola-miR-122, ola-miR-92a, ola-miR-125a-3p showed maternal expression only. The expression of six out of seven miRNAs significantly decreased after maternal-zygotic transition (MZT), whereas ola-miR-430a expression initiated only after MZT. The temporal dynamic expression of these miRNAs suggests their potential roles in early embryogenesis and genome-zygotic activation in medaka.

Genetic and functional insights into the fractal structure of the heart.

The inner surfaces of the human heart are covered by a complex network of muscular strands that is thought to be a remnant of embryonic development1,2. The function of these trabeculae in adults and their genetic architecture are unknown. Here we performed a genome-wide association study to investigate image-derived phenotypes of trabeculae using the fractal analysis of trabecular morphology in 18,096 participants of the UK Biobank. We identified 16 significant loci that contain genes associated with haemodynamic phenotypes and regulation of cytoskeletal arborization3,4. Using biomechanical simulations and observational data from human participants, we demonstrate that trabecular morphology is an important determinant of cardiac performance. Through genetic association studies with cardiac disease phenotypes and Mendelian randomization, we find a causal relationship between trabecular morphology and risk of cardiovascular disease. These findings suggest a previously unknown role for myocardial trabeculae in the function of the adult heart, identify conserved pathways that regulate structural complexity and reveal the influence of the myocardial trabeculae on susceptibility to cardiovascular disease.

Developmental abnormalities and epigenetic alterations in medaka (Oryzias latipes) embryos induced by triclosan exposure.

Triclosan (TCS), an antibacterial and antifungal agent present in some consumer products, has been detected in the environment at varying concentrations. TCS exposure has been found to cause developmental abnormalities and endocrine disruption in various species of fish. It is not clearly understood whether TCS exposure causes epigenetic alterations in developing embryos and their germ cells. In the present study, we examined the effects of TCS exposure (0, 50, 100 and, 200 µg/L) on embryonic development and primordial germ cells (PGCs), which are precursors of sperm and eggs, in medaka (Oyzias latipes). Developmental TCS exposure from 8 h post-fertilization through 15 days post-fertilization (dpf) resulted in several developmental abnormalities, including enlarged yolk sac, decreased head trunk angle (HTA), and severe edema in the pericardial region. The male ratio increased in the 100 µg/L TCS exposure group, which was negatively correlated with the expression of cyp19ala (a gene encoding aromatase) and arα (androgen receptor alpha). Developmental 50 µg/L TCS exposure resulted in global hypomethylation in the whole body but not in the isolated PGCs. Expression of the gene encoding DNA methyltransferases (dnmt1 and dnmt3aa) was decreased by 50 µg/L TCS exposure both in the whole body and PGCs. TCS altered the expression of genes encoding enzymes involved in DNA methylation and demethylation in PGCs, suggesting epigenetic effects on germ cells. The present results demonstrate that the embryos exposed to the tested concentrations of TCS develop deformities during the early life stages and that the TCS within this range possesses endocrine disrupting properties potential enough to alter sex ratios of developing embryos.

Dynamic transcriptional and chromatin accessibility landscape of medaka embryogenesis.

Medaka (Oryzias latipes) has become an important vertebrate model widely used in genetics, developmental biology, environmental sciences, and many other fields. A high-quality genome sequence and a variety of genetic tools are available for this model organism. However, existing genome annotation is still rudimentary, as it was mainly based on computational prediction and short-read RNA-seq data. Here we report a dynamic transcriptome landscape of medaka embryogenesis profiled by long-read RNA-seq, short-read RNA-seq, and ATAC-seq. By integrating these data sets, we constructed a much-improved gene model set including about 17,000 novel isoforms and identified 1600 transcription factors, 1100 long noncoding RNAs, and 150,000 potential cis-regulatory elements as well. Time-series data sets provided another dimension of information. With the expression dynamics of genes and accessibility dynamics of cis-regulatory elements, we investigated isoform switching, as well as regulatory logic between accessible elements and genes, during embryogenesis. We built a user-friendly medaka omics data portal to present these data sets. This resource provides the first comprehensive omics data sets of medaka embryogenesis. Ultimately, we term these three assays as the minimum ENCODE toolbox and propose the use of it as the initial and essential profiling genomic assays for model organisms that have limited data available. This work will be of great value for the research community using medaka as the model organism and many others as well.

Transcriptional responses in newly-hatched Japanese medaka (Oryzias latipes) associated with developmental malformations following diluted bitumen exposure.

Japanese medaka embryos were exposed to water accommodated fractions (WAF) and chemically-enhanced WAF of two types of diluted bitumen (dilbit) at concentrations bracketing the EC50s for developmental malformations. Within these treatments, fish were grouped based on the presence or absence of developmental malformations (e.g., blue sac disease (BSD)), and analyzed for novel transcriptomic responses. Microarray analyses identified novel biomarkers and gene networks in dilbit-exposed malformed embryos that were not evident in dilbit-exposed fish without BSD or in controls without dilbit. The top differentially expressed genes (DEGs) included cytochrome P450 transcripts (cyp1) in fish from all dilbit treatments (malformed and non-malformed fish), as well as: fibroblast growth factor (fgf7), AHR repressor (ahrr), and squalene monooxygenase (sqle). In dilbit-exposed fish that did not develop BSD, the only reported individual DEG was eukaryotic translation initiation factor 3 subunit D (eif3d). However, a number of other pathways were enriched, including melatonin effects on circadian clock and the antioxidant response, estrogen and androgen metabolism as well as many receptor signaling pathways. Pathways associated with hedgehog, steroid biosynthesis, and Wnt signaling were significantly altered between low and high concentrations of dilbit exposure. An effect of the dispersant control on swim bladder development was observed at concentrations 10-fold higher than those used to disperse dilbit, and a number of gene targets unique to fish in this comparison were affected. This suggests that the toxic effects of dispersant may involve alternative mechanisms to dilbit, but cause similar phenotypic responses. This study identified novel biomarkers in fish exposed to dilbit, with or without visual malformations, that can be used to assess the risks of dilbit to aquatic ecosystem health.