Developmental transcriptomics in Atlantic haddock: Illuminating pattern formation and organogenesis in non-model vertebrates

Structure-activity relationships for alkyl-phenanthrenes support two independent but interacting synergistic models for PAC mixture potency

Multiple lines of evidence at whole animal, cellular and molecular levels implicate polycyclic aromatic compounds (PACs) with three rings as drivers of crude oil toxicity to developing fish. Phenanthrene (P0) and its alkylated homologs (C1- through C4-phenanthrenes) comprise the most prominent subfraction of tricyclic PACs in crude oils. Among this family, P0 has been studied intensively, with more limited detail available for the C4-phenanthrene 1-methyl-7-isopropyl-phenanthrene (1-M,7-IP, or retene). While both compounds are cardiotoxic, P0 impacts embryonic cardiac function and development through direct blockade of K+ and Ca2+ currents that regulate cardiomyocyte contractions. In contrast, 1-M,7-IP dysregulates aryl hydrocarbon receptor (AHR) activation in developing ventricular cardiomyocytes. Although no other compounds have been assessed in detail across the larger family of alkylated phenanthrenes, increasing alkylation might be expected to shift phenanthrene family member activity from K+/Ca2+ ion current blockade to AHR activation. Using embryos of two distantly related fish species, zebrafish and Atlantic haddock, we tested 14 alkyl-phenanthrenes in both acute and latent developmental cardiotoxicity assays. All compounds were cardiotoxic, and effects were resolved into impacts on multiple, highly specific aspects of heart development or function. Craniofacial defects were clearly linked to developmental cardiotoxicity. Based on these findings, we suggest a novel framework to delineate the developmental toxicity of petrogenic PAC mixtures in fish, which incorporates multi-mechanistic pathways that produce interactive synergism at the organ level. In addition, relationships among measured embryo tissue concentrations, cytochrome P4501A mRNA induction, and cardiotoxic responses suggest a two-compartment toxicokinetic model that independently predicts high potency of PAC mixtures through classical metabolic synergism. These two modes of synergism, specific to the sub-fraction of phenanthrenes, are sufficient to explain the high embryotoxic potency of crude oils, independent of as-yet unmeasured compounds in these complex environmental mixtures.

Cardiac toxicity of phenanthrene depends on developmental stage in Atlantic cod (Gadus morhua)

Complex mixtures like crude oil, and single components such as Phenanthrene (Phe), induce cardiotoxicity by interfering with excitation-contraction coupling. However, recent work has demonstrated that the timing of pollutant exposure during embryogenesis greatly impacts the degree of cardiac dysfunction caused. Here, we aimed to clarify the temporal dependence of Phe toxicity and the downstream effects of cardiac dysfunction using Atlantic cod (Gadus morhua). Phe (nominal concentration, 1.12 μmol/L), or the L-type‑calcium channel blocker Nicardipine (Nic) (nominal concentration, 2 and 4 μmol/L), were individually applied to cod embryos either during cardiogenesis (early) or after the onset of cardiac function (late). Phe toxicity was highly dependent on the timing of exposure. Exposure after the onset of cardiac function (i.e. late) caused more severe cardiac and extracardiac abnormalities at 3 days post hatching (dph) than early exposure. Late Phe exposure resulted in a smaller ventricle, eliminated ventricular contraction, and reduced atrial contraction. In contrast, early Phe exposure did not have an effect on cardiac development and function. This temporal difference was not as evident in the Nic treatment. Early Nic exposure created similar morphological phenotypes to the late Phe exposure. The two treatments (early Nic and late Phe) also shared a cardiofunctional phenotype, comprised of eliminated ventricular, and reduced atrial, contraction. These data suggest that extracardiac abnormalities, such as the craniofacial deformities seen after late embryonic exposure to cardiotoxic oil components and mixtures, are mostly downstream effects of cardiac dysfunction.

Transcriptomic responses and apoptosis in larval red drum (Sciaenops ocellatus) co-exposed to crude oil and ultraviolet (UV) radiation

Ultraviolet (UV) radiation can significantly increase the toxicity of polycyclic aromatic hydrocarbons (PAHs) in crude oil to early life stage (ELS) fishes through photo-induced /photo-enhanced toxicity. However, little is known about the sub-lethal effects and mechanisms of photo-induced PAH toxicity in ELS fishes. The present study investigated apoptosis and global transcriptomic effects in larval red drum (Sciaenops ocellatus) (24-72 h post-fertilization) following co-exposure to oil (0.29-0.30 μg/L ∑PAH₅₀) and UV. Apoptosis was quantified using the TUNEL assay, and transcriptomic effects were assessed using RNA sequencing analysis. Apoptotic fluorescence was greatest in the eyes and skin following 24 and 48 h co-exposure to oil and UV, indicating photo-induced toxicity. Consistent with these phenotypic responses, pathways associated with phototransduction, eye development, and dermatological disease were among the top predicted pathways impacted. The present study is the first to provide global transcriptomic analysis of UV and oil co-exposure in an ELS fish.

Screening and Validation of p38 MAPK Involved in Ovarian Development of Brachymystax lenok

Brachymystax lenok (lenok) is a rare cold-water fish native to China that is of high meat quality. Its wild population has declined sharply in recent years, and therefore, exploring the molecular mechanisms underlying the development and reproduction of lenoks for the purposes of artificial breeding and genetic improvement is necessary. The lenok comparative transcriptome was analyzed by combining single molecule, real-time, and next generation sequencing (NGS) technology. Differentially expressed genes (DEGs) were identified in five tissues (head kidney, spleen, liver, muscle, and gonad) between immature [300 days post-hatching (dph)] and mature [three years post-hatching (ph)] lenoks. In total, 234,124 and 229,008 full-length non-chimeric reads were obtained from the immature and mature sequencing data, respectively. After NGS correction, 61,405 and 59,372 non-redundant transcripts were obtained for the expression level and pathway enrichment analyses, respectively. Compared with the mature group, 719 genes with significantly increased expression and 1,727 genes with significantly decreased expression in all five tissues were found in the immature group. Furthermore, DEGs and pathways involved in the endocrine system and gonadal development were identified, and p38 mitogen-activated protein kinases (MAPKs) were identified as potentially regulating gonadal development in lenok. Inhibiting the activity of p38 MAPKs resulted in abnormal levels of gonadotropin-releasing hormone, follicle-stimulating hormone, and estradiol, and affected follicular development. The full-length transcriptome data obtained in this study may provide a valuable reference for the study of gene function, gene expression, and evolutionary relationships in B. lenok and may illustrate the basic regulatory mechanism of ovarian development in teleosts.

Cardiac dysfunction affects eye development and vision by reducing supply of lipids in fish

Developing organisms are especially vulnerable to environmental stressors. Crude oil exposure in early life stages of fish result in multiple functional and developmental defects, including cardiac dysfunction and abnormal and smaller eyes. Phenanthrene (Phe) has a reversible impact on cardiac function, and under exposure Phe reduces cardiac contractility. Exposure to a known L-type channel blocker, nicardipine hydrochloride (Nic) also disrupts cardiac function and creates eye deformities. We aimed to investigate whether cardiac dysfunction was the major underlying mechanism of crude oil-, Phe- and Nic-induced eye malformations. We exposed Atlantic haddock (Melanogrammus aeglefinus) early embryos to Nic and crude oil (Oil) and late embryos/early larvae to Phe exposure. All three exposures resulted in cardiac abnormalities and lead to severe, eye, jaw and spinal deformities at early larval stages. At 3 days post hatching, larvae from the exposures and corresponding controls were dissected. Eyes, trunk, head and yolk sac were subjected to lipid profiling, and eyes were also subjected to transcriptomic profiling. Among most enriched pathways in the eye transcriptomes were fatty acid metabolism, calcium signaling and phototransduction. Changes in lipid profiles and the transcriptome suggested that the dysfunctional and abnormal eyes in our exposures were due to both disruption of signaling pathways and insufficient supply of essential fatty acids and other nutrients form the yolk.

Low-level embryonic crude oil exposure disrupts ventricular ballooning and subsequent trabeculation in Pacific herring

There is a growing awareness that transient, sublethal embryonic exposure to crude oils cause subtle but important forms of delayed toxicity in fish. While the precise mechanisms for this loss of individual fitness are not well understood, they involve the disruption of early cardiogenesis and a subsequent pathological remodeling of the heart much later in juveniles. This developmental cardiotoxicity is attributable, in turn, to the inhibitory actions of crude oil-derived mixtures of polycyclic aromatic compounds (PACs) on specific ion channels and other proteins that collectively drive the rhythmic contractions of heart muscle cells via excitation-contraction coupling. Here we exposed Pacific herring (Clupea pallasi) embryos to oiled gravel effluent yielding ΣPAC concentrations as low as ~ 1 μg/L (64 ng/g in tissues). Upon hatching in clean seawater, and following the depuration of tissue PACs (as evidenced by basal levels of cyp1a gene expression), the ventricles of larval herring hearts showed a concentration-dependent reduction in posterior growth (ballooning). This was followed weeks later in feeding larvae by abnormal trabeculation, or formation of the finger-like projections of interior spongy myocardium, and months later with hypertrophy (overgrowth) of the spongy myocardium in early juveniles. Given that heart muscle cell differentiation and migration are driven by Ca²⁺-dependent intracellular signaling, the observed disruption of ventricular morphogenesis was likely a secondary (downstream) consequence of reduced calcium cycling and contractility in embryonic cardiomyocytes. We propose defective trabeculation as a promising phenotypic anchor for novel morphometric indicators of latent cardiac injury in oil-exposed herring, including an abnormal persistence of cardiac jelly in the ventricle wall and cardiomyocyte hyperproliferation. At a corresponding molecular level, quantitative expression assays in the present study also support biomarker roles for genes known to be involved in muscle contractility (atp2a2, myl7, myh7), cardiomyocyte precursor fate (nkx2.5) and ventricular trabeculation (nrg2, and hbegfa). Overall, our findings reinforce both proximal and indirect roles for dysregulated intracellular calcium cycling in the canonical fish early life stage crude oil toxicity syndrome. More work on Ca²⁺-mediated cellular dynamics and transcription in developing cardiomyocytes is needed. Nevertheless, the highly specific actions of ΣPAC mixtures on the heart at low, parts-per-billion tissue concentrations directly contravene classical assumptions of baseline (i.e., non-specific) crude oil toxicity.

Untangling mechanisms of crude oil toxicity: Linking gene expression, morphology and PAHs at two developmental stages in a cold-water fish

Early life stages of fish are highly sensitive to crude oil exposure and thus, short term exposures during critical developmental periods could have detrimental consequences for juvenile survival. Here we administered crude oil to Atlantic haddock (Melanogrammus aeglefinus) in short term (3-day) exposures at two developmental time periods: before first heartbeat, from gastrulation to cardiac cone stage (early), and from first heartbeat to one day before hatching (late). A frequent sampling regime enabled us to determine immediate PAH uptake, metabolite formation and gene expression changes. In general, the embryotoxic consequences of an oil exposure were more severe in the early exposure animals. Oil droplets on the eggshell resulted in severe cardiac and craniofacial abnormalities in the highest treatments. Gene expression changes of Cytochrome 1 a, b, c and d (cyp1a, b, c, d), Bone morphogenetic protein 10 (bmp10), ABC transporter b1 (abcb1) and Rh-associated G-protein (rhag) were linked to PAH uptake, occurrence of metabolites of phenanthrene and developmental and functional abnormalities. We detected circulation-independent, oil-induced gene expression changes and separated phenotypes linked to proliferation, growth and disruption of formation events at early and late developmental stages. Changes in bmp10 expression suggest a direct oil-induced effect on calcium homeostasis. Localized expression of rhag propose an impact on osmoregulation. Severe eye abnormalities were linked to possible inappropriate overexpression of cyp1b in the eyes. This study gives an increased knowledge about developmentally dependent effects of crude oil toxicity. Thus, our findings provide more knowledge and detail to new and several existing adverse outcome pathways of crude oil toxicity.

Environmental effects of offshore produced water discharges: A review focused on the Norwegian continental shelf

Produced water (PW), a large byproduct of offshore oil and gas extraction, is reinjected to formations or discharged to the sea after treatment. The discharges contain dispersed crude oil, polycyclic aromatic hydrocarbons (PAHs), alkylphenols (APs), metals, and many other constituents of environmental relevance. Risk-based regulation, greener offshore chemicals and improved cleaning systems have reduced environmental risks of PW discharges, but PW is still the largest operational source of oil pollution to the sea from the offshore petroleum industry. Monitoring surveys find detectable exposures in caged mussel and fish several km downstream from PW outfalls, but biomarkers indicate only mild acute effects in these sentinels. On the other hand, increased concentrations of DNA adducts are found repeatedly in benthic fish populations, especially in haddock. It is uncertain whether increased adducts could be a long-term effect of sediment contamination due to ongoing PW discharges, or earlier discharges of oil-containing drilling waste. Another concern is uncertainty regarding the possible effect of PW discharges in the sub-Arctic Southern Barents Sea. So far, research suggests that sub-arctic species are largely comparable to temperate species in their sensitivity to PW exposure. Larval deformities and cardiac toxicity in fish early life stages are among the biomarkers and adverse outcome pathways that currently receive much attention in PW effect research. Herein, we summarize the accumulated ecotoxicological knowledge of offshore PW discharges and highlight some key remaining knowledge needs.

Differential impacts of carp and salmon pituitary extracts on induced oogenesis, egg quality, molecular ontogeny and embryonic developmental competence in European eel

Low egg quality and embryonic survival are critical challenges in aquaculture, where assisted reproduction procedures and other factors may impact egg quality. This includes European eel (Anguilla anguilla), where pituitary extract from carp (CPE) or salmon (SPE) is applied to override a dopaminergic inhibition of the neuroendocrine system, preventing gonadotropin secretion and gonadal development. The present study used either CPE or SPE to induce vitellogenesis in female European eel and compared impacts on egg quality and offspring developmental competence with emphasis on the maternal-to-zygotic transition (MZT). Females treated with SPE produced significantly higher proportions of floating eggs with fewer cleavage abnormalities and higher embryonic survival. These findings related successful embryogenesis to higher abundance of mRNA transcripts of genes involved in cell adhesion, activation of MZT, and immune response (dcbld1, epcam, oct4, igm) throughout embryonic development. The abundance of mRNA transcripts of cldnd, foxr1, cea, ccna1, ccnb1, ccnb2, zar1, oct4, and npm2 was relatively stable during the first eight hours, followed by a drop during MZT and low levels thereafter, indicating transfer and subsequent clearance of maternal mRNA. mRNA abundance of zar1, epcam, and dicer1 was associated with cleavage abnormalities, while mRNA abundance of zar1, sox2, foxr1, cldnd, phb2, neurod4, and neurog1 (before MZT) was associated with subsequent embryonic survival. In a second pattern, low initial mRNA abundance with an increase during MZT and higher levels persisting thereafter indicating the activation of zygotic transcription. mRNA abundance of ccna1, npm2, oct4, neurod4, and neurog1 during later embryonic development was associated with hatch success. A deviating pattern was observed for dcbld1, which mRNA levels followed the maternal-effect gene pattern but only for embryos from SPE treated females. Together, the differences in offspring production and performance reported in this study show that PE composition impacts egg quality and embryogenesis and in particular, the transition from initial maternal transcripts to zygotic transcription.

mRNA-miRNA-Seq Reveals Neuro-Cardio Mechanisms of Crude Oil Toxicity in Red Drum ( Sciaenops ocellatus)

Polycyclic aromatic hydrocarbons (PAHs) present in crude oil can cause global gene dysregulation and developmental impairment in fish. However, the mechanisms that alter gene regulation are not well understood. In this study, larval red drum ( Sciaenops ocellatus) were exposed to water accommodated fractions of source oil (6.8, 13.7, and 35.9 μg/L total PAHs) and weathered slick oil (4.7, 8.1, and 18.0 μg/L total PAHs) from the Deepwater Horizon (DWH) oil spill. The global mRNA-microRNA functional networks associated with the toxicity of DWH oil were explored by next-generation sequencing and in-depth bioinformatics analyses. Both source and slick oil significantly altered the expression of miR-18a, miR-27b, and miR-203a across all exposure concentrations. Consistent with the observed concentration-dependent morphological changes, the target mRNAs of these microRNAs were predominantly involved in neuro-cardio system development processes and associated key signaling pathways such as axonal guidance signaling, cAMP-response-element-binding protein signaling in neurons, calcium signaling, and nuclear-factor-of-activated T cells signaling in cardiac hypertrophy. The results indicated that the developmental toxicity of crude oil may result from the abnormal expression of microRNAs and associated target genes, especially for the nervous system. Moreover, we provide a case study for systematic toxicity evaluation leveraging mRNA-microRNA-seq data using nonmodel species.

Larval Red Drum (Sciaenops ocellatus) Sublethal Exposure to Weathered Deepwater Horizon Crude Oil: Developmental and Transcriptomic Consequences

The Deepwater Horizon (DWH) incident resulted in extensive oiling of the pelagic zone and shoreline habitats of many commercially important fish species. Exposure to the water-accommodated fraction (WAF) of oil from the spill causes developmental toxicity through cardiac defects in pelagic fish species. However, few studies have evaluated the effects of the oil on near-shore estuarine fish species such as red drum (Sciaenops ocellatus). Following exposure to a certified weathered slick oil (4.74 μg/L ∑PAH₅₀) from the DWH event, significant sublethal impacts were observed ranging from impaired nervous system development [average 17 and 22% reductions in brain and eye area at 48 h postfertilization (hpf), respectively] to abnormal cardiac morphology (100% incidence at 24, 48, and 72 hpf) in red drum larvae. Consistent with the phenotypic responses, significantly differentially expressed transcripts, enriched gene ontology, and altered functions and canonical pathways predicted adverse outcomes in nervous and cardiovascular systems, with more pronounced changes at later larval stages. Our study demonstrated that the WAF of weathered slick oil of DWH caused morphological abnormalities predicted by a suite of advanced bioinformatic tools in early developing red drum and also provided the basis for a better understanding of molecular mechanisms of crude oil toxicity in fish.

Developmental transcriptomic analyses for mechanistic insights into critical pathways involved in embryogenesis of pelagic mahi-mahi (Coryphaena hippurus)

Mahi-mahi (Coryphaena hippurus) is a commercially and ecologically important species of fish occurring in tropical and temperate waters worldwide. Understanding early life events is crucial for predicting effects of environmental stress, which is largely restricted by a lack of genetic resources regarding expression of early developmental genes and regulation of pathways. The need for anchoring developmental stages to transcriptional activities is highlighted by increasing evidence on the impacts of recurrent worldwide oil spills in this sensitive species during early development. By means of high throughput sequencing, we characterized the developmental transcriptome of mahi-mahi at three critical developmental stages, from pharyngula embryonic stage (24 hpf) to 48 hpf yolk-sac larva (transition 1), and to 96 hpf free-swimming larva (transition 2). With comparative analysis by multiple bioinformatic tools, a larger number of significantly altered genes and more diverse gene ontology terms were observed during transition 2 than transition 1. Cellular and tissue development terms were more significantly enriched in transition 1, while metabolism related terms were more enriched in transition 2, indicating a switch progressing from general embryonic development to metabolism during the two transitions. Special focus was given on the most significant common canonical pathways (e.g. calcium signaling, glutamate receptor signaling, cAMP response element-binding protein signaling, cardiac β-adrenergic signaling, etc.) and expression of developmental genes (e.g. collagens, myosin, notch, glutamate metabotropic receptor etc.), which were associated with morphological changes of nervous, muscular, and cardiovascular system. These data will provide an important basis for understanding embryonic development and identifying molecular mechanisms of abnormal development in fish species.

Novel adverse outcome pathways revealed by chemical genetics in a developing marine fish

Crude oil spills are a worldwide ocean conservation threat. Fish are particularly vulnerable to the oiling of spawning habitats, and crude oil causes severe abnormalities in embryos and larvae. However, the underlying mechanisms for these developmental defects are not well understood. Here, we explore the transcriptional basis for four discrete crude oil injury phenotypes in the early life stages of the commercially important Atlantic haddock (Melanogrammus aeglefinus). These include defects in (1) cardiac form and function, (2) craniofacial development, (3) ionoregulation and fluid balance, and (4) cholesterol synthesis and homeostasis. Our findings suggest a key role for intracellular calcium cycling and excitation-transcription coupling in the dysregulation of heart and jaw morphogenesis. Moreover, the disruption of ionoregulatory pathways sheds new light on buoyancy control in marine fish embryos. Overall, our chemical-genetic approach identifies initiating events for distinct adverse outcome pathways and novel roles for individual genes in fundamental developmental processes.