RNA seq- and DEG-based comparison of developmental toxicity in fish embryos of two species exposed to Iranian heavy crude oil (2023)

Marine sponges play important roles in benthic ecosystems. More than providing shelter and food to other species, they help maintain water quality by regulating nitrogen and ammonium levels in the water, and bioaccumulate heavy metals. This system, however, is particularly sensitive to sudden environmental changes including catastrophic pollution event such as oil spills. Hundreds of oil platforms are currently actively extracting oil and gas in the Gulf of Mexico. To test the vulnerability of the benthic ecosystems to oil spills, we utilized the Caribbean reef sponge, Cinachyrella alloclada, as a novel experimental indicator. We have exposed organisms to crude oil and oil dispersant for up to 24h and measured resultant gene expression changes. Our findings indicate that 1-hour exposure to water accommodated fractions (WAF) was enough to elicit massive shifts in gene expression in sponges and host bacterial communities (8052 differentially expressed transcripts) with the up-regulation of stress related pathways, cancer related pathways, and cell integrity pathways. Genes that were upregulated included heat shock proteins, apoptosis, oncogenes (Rab/Ras, Src, CMYC), and several E3 ubiquitin ligases. 24-hour exposure of chemically enhanced WAF (CE-WAF) had the greatest impact to benthic communities, resulting in mostly downregulation of gene expression (4248 differentially expressed transcripts). Gene deregulation from 1-hour treatments follow this decreasing trend of toxicity: WAF>CE-WAF>Dispersant, while the 24-hour treatment showed a shift to CE-WAF>Dispersant>WAF in our experiments. Thus, this study supports the development of Cinachyrella alloclada as a research model organism and bioindicator species for Florida reefs and underscores the importance of developing more efficient and safer ways to remove oil in the event of a spill catastrophe.

A new generation of booster biocides that include metal pyrithiones (PTs) such as copper pyrithione (CuPT) and zinc pyrithione (ZnPT) are being used as tributyltin alternatives. In the marine environment, ZnPT can easily transchelate Cu to form CuPT, and the environmental fate and persistence of these two metal pyrithiones are closely related. Although some data on the toxicity of biocides on marine fish are available, little is known about their toxicity and toxic pathway. We thus compared the toxic effects of CuPT and ZnPT on embryonic olive flounder (Paralichthys olivaceus) by investigating their adverse effects based on developmental morphogenesis and transcriptional variation. In our study, the toxic potency of CuPT was greater with respect to developmental malformation and mortality than ZnPT. Consistent with the developmental effects, the expression of genes related to tail fin malformation (including plod2, furin, and wnt3a) was higher in embryonic flounder exposed to CuPT than in those exposed to ZnPT. Genes related to muscle and nervous system development exhibited significant changes on differential gene expression profiles using RNA sequencing (cutoff value P< 0.05). Gene ontology analysis of embryos exposed to CuPT revealed affected cellular respiration and kidney development, whereas genes associated with cell development, nervous system development and heart development showed significant variation in embryonic flounder exposed to ZnPT. Overall, our study clarifies the common and unique developmental toxic effects of CuPT and ZnPT through transcriptomic analyses in embryonic flounder.

Aquatic sediments act as a storage for diverse mixtures of organic and inorganic contaminants. Nevertheless, most evaluations of contaminated sediments have been limited to the assessment of concentrations of target compounds and lethal effects on some test species. To identify the organic contaminants causing sub-lethal effects of contaminated sediment, this study combined chemical and biological (non)target analysis involving comprehensive two-dimensional gas chromatography coupled with a time-of-flight Mass Spectrometer (GCxGC/ToF-MS) analysis, embryonic malformation and high-throughput sequencing (RNA-seq) analysis on developing flounder. Polycyclic aromatic hydrocarbons were more abundant in the sediment extract of Yeongil Bay (SEY), while Jinhae Bay (SEJ) was contaminated with a large amount of unidentified chemicals. The unidentified chemicals of SEJ included branched alkanes, oxygenated cycloalkanes, heterogeneous hydrocarbons, and other unknown compounds. Percentage of pericardial edema was the highest in embryonic flounder exposed to SEY. Consistent with the morphogenesis results, the expression level of genes related to heart formation including the nkx2.5 and robo1 was greater in embryonic flounder exposed to SEY. In the analyses of differential gene expression profiles (cutoff P < 0.05), by RNA-seq, embryos exposed to SEJ showed changes related to cell differentiation, cell part morphogenesis, neurogenesis, and neuron development. Genes related to neurogenesis and positive regulation of molecular functions variated significantly in embryos exposed to SEY. These results demonstrated the advantages of combining target and non-target analysis to accurately evaluate the major chemical groups causing sediment toxicity. Therefore, this work provided a useful approach to tracking and revealing the causes of toxic effects and identifying potential toxic mechanisms.

To further understand the underlying mechanisms involved in the developmental toxicity of crude oil and chemically dispersed crude oil on fish early-life stages (ELS), zebrafish (Danio rerio) embryos were exposed to GM-2 chemical dispersant (DISP), low-energy water-accommodated fractions (LEWAF), and chemically enhanced WAF (CEWAF) of Merey crude oil at sublethal concentrations for 120h. We employed the General Morphology Score (GMS) and General Teratogenic Score (GTS) systems in conjunction with high-throughput RNA-Seq analysis to evaluate the phenotypic and transcriptomic responses in zebrafish ELS. Results showed that ΣPAHs concentrations in LEWAF and CEWAF solutions were 507.63±80.95ng·L⁻¹ and 4039.51±241.26ng·L⁻¹, respectively. The GMS and GTS values indicated that CEWAF exposure caused more severe developmental delay and higher frequencies of teratogenic effects than LEWAF exposure. Moreover, no significant change in heart rate was observed in LEWAF treatment, while CEWAF exposure caused a significant reduction in heart rate. LEWAF and CEWAF exposure exhibited an overt change in eye area, with a reduction of 4.0% and 25.3% (relative to the control), respectively. Additionally, no obvious impact on phenotypic development was observed in zebrafish embryo-larvae following DISP exposure. Significant changes in gene expression were detected in LEWAF and CEWAF treatments, with a total of 957 and 2062 differentially expressed genes (DEGs), respectively, while DISP exposure altered only 91 DEGs. Functional enrichment analysis revealed that LEWAF and CEWAF exposure caused significant perturbations in the pathways associated with phototransduction, retinol metabolism, metabolism of xenobiotics by cytochrome P450, and immune response-related pathways. Our results provide more valid evidence to corroborate the previous suggestion that ocular impairment is an equal or possibly more sensitive biomarker than cardiotoxicity in fish ELS exposed to oil-derived PAHs. All these findings could gain further mechanistic insights into the effects of crude oil and chemical dispersant on fish ELS.

Chemical contaminants can be discharged by vessel hull cleaning processes, such as scraping, jet spraying, and painting, all of which produce readily transportable contaminants into the marine environment, where they are referred to as ‘hotspots’ of contamination in coastal areas. However, many countries have not yet established effective evaluation methods for disposal of waste mixtures or management guidelines for areas of hull cleaning. To define the toxic effects of wastewater from vessel hull cleaning in dry docks on resident non-target organisms, we investigated the chemical concentrations and developmental toxicity on embryonic flounder, which is an organism sensitive to chemical contamination. In this study, the dominant inorganic metal discharged was zinc when cleaning Ship A (300 tons) and copper for Ship B (5,000 tons). The wastewater from high-pressure water blasting (WHPB) of Ship A (300 tons) and Ship B (5,000 tons) produced a largely overlapping suite of developmental malformations including pericardial edema, spinal curvature, and tail fin defects. Forty-eight hours after exposure, the frequency percentage of malformation began to increase in embryos exposed to a 500-fold dilution of WHPB from Ships A and B. We performed transcriptome sequencing to characterize the toxicological developmental effects of WHPB exposure at the molecular level. The results of the analysis revealed significantly altered expression of genes associated with muscle cell differentiation, actin-mediated cell contraction, and nervous system development (cutoff P < 0.01) in embryonic flounder exposed to high-pressure cleaning effluent from Ship A. Genes associated with chromatin remodeling, cell cycling, and insulin receptor signaling pathways were significantly altered in embryonic flounder exposed to WHPB of Ship B (cutoff P < 0.01). These findings provide a greater understanding of the developmental toxicity and potential effects of WHPB effluent on coastal embryonic fish. Furthermore, our results could inform WHPB effluent management practices to reduce impacts on non-target coastal organisms.

Ecotoxicology and Environmental Safety, Volume 146, 2017, pp. 76-82

During the Deepwater Horizon Natural Resource Damage Assessment, gross morphologic cardiac abnormalities, including softer, more distensible musculature, were noted upon gross necropsy in hearts from laughing gulls and double-crested cormorants exposed to weathered MC252 crude oil. A species specific, echocardiographic technique was developed for antemortem evaluation of function that was used to evaluate and better characterize cardiac dysfunction. Control (n=12) and treated (n=13) cormorant groups of similar sex-ratio and ages were dermally treated with approximately 13ml of water or weathered MC252 crude oil, respectively, every 3 days for 6 dosages. This resulted in a low to moderate external exposure. Upon visualization and clinical assessment of the hearts of all test subjects, comprehensive diagnostic cardiographic measurements were taken twice, prior to oil application and after a 21day dermal oil exposure. Oil-treated birds showed a decrease in cardiac systolic function, as characterized by an increased left ventricular internal dimension-systole and left ventricular stroke volume as well as concurrent decreased left ventricular ejection fraction and left ventricular fractional shortening when compared to both control birds’ and the treated birds’ time zero values. These changes are indicative of a possible dilative cardiomyopathy induced by oil exposure, although further elucidation of possible collagen damage is recommended. Arrhythmias including tachycardia in two treated birds and bradycardia in all treated birds were documented, indicating further clinically significant abnormalities induced by MC252 oil that warrant further investigation. A statistically significant increase in free calcium concentration, important to muscular and neurologic function in treated birds was also noted. This study documents that weathered MC252 oil caused clinically significant cardiac dysfunction that could result in mortality and decrease recruitment.

Fish & Shellfish Immunology, Volume 94, 2019, pp. 99-112

The lack of practical control measures for pharyngeal myxobolosis is becoming an important limiting factor for the sustainable development of the gibel carp (Carassius auratus gibelio) culture industry in China. Myxobolus honghuensis has been identified as the causative agent of this pandemic disease, which exclusively infects the pharynx of gibel carp, a potential important mucosal lymphoid-associated tissue (MLAT). Myxozoa generally initiate invasion through the mucosal tissues of fish, where some of them also complete their sporogonial stages. However, the pharynx-associated immune responses of teleost against myxosporeans infection remain unknown. Here, a de novo transcriptome assembly of the pharynx of gibel carp naturally infected with M. honghuensis was performed for the first time, using RNA-seq. Comparative analysis of severely infected and mildly infected pharyngeal tissues (SI group and MI group) from the same fish individuals and control pharyngeal tissues (C group) from the uninfected fish was carried out to investigate the potential mucosal immune function of the fish pharynx, and characterize the panoramic picture of pharynx local mucosal immune responses of gibel carp against the M. honghuensis infection. A total of 242,341 unigenes were obtained and pairwise comparison resulted in 13,009 differentially-expressed genes (DEGs) in the SI/C group comparison, 6014 DEGs in the MI/C group comparison, and 9031 DEGs in the SI/MI group comparison. Comprehensive analysis showed that M. honghuensis infection elicited a significant parasite load-dependent alteration of the expression of numerous innate and adaptive immune-related genes in the local lesion tissue. Innate immune molecules, including mucins, toll-like receptors, C-type lectin, serum amyloid A, cathepsins and complement components were significantly up-regulated in the SI group compared with the C group. Up-regulation of genes involved in apoptosis signaling pathway and the IFN-mediated immune system were found in the SI group, suggesting these two pathways played a crucial role in innate immune response to M. honghuensis infection. Up-regulation of chemokines and chemokine receptors and the induction of the leukocyte trans-endothelial migration pathways in the severely and mildly infected pharynx suggested that many leucocytes were recruited to the local infected sites to mount a strong mucosal immune responses against the myxosporean infection. Up-regulation of CD3D, CD22, CD276, IL4/13A, GATA3, arginase 2, IgM, IgT and pIgR transcripts provided strong evidences for the presence of T/B cells and specific mucosal immune responses at local sites with M. honghuensis infection. Our results firstly demonstrated the mucosal function of the teleost pharynx and provided evidences of intensive local immune defense responses against this mucosa-infecting myxosporean in the gibel carp pharynx. Pharyngeal myxobolosis was shaped by a prevailing anti-inflammatory response pattern during the advanced infection stages. Further understanding of the functional roles of fish immune molecules involved in the initial invasion and/or final sporogony site may facilitate future development of control strategies for this myxobolosis.

Fish & Shellfish Immunology, Volume 79, 2018, pp. 244-255

Aquatic Toxicology, Volume 201, 2018, pp. 174-186

Polycyclic aromatic hydrocarbons such as benzo[a]pyrene (BaP) that activate the aryl hydrocarbon receptor (Ahr) pathway, and endocrine disruptors acting through the estrogen receptor pathway are among environmental pollutants of major concern. In this work, we exposed Atlantic cod (Gadus morhua) precision-cut liver slices (PCLS) to BaP (10 nM and 1000 nM), ethynylestradiol (EE2) (10 nM and 1000 nM), and equimolar mixtures of BaP and EE2 (10 nM and 1000 nM) for 48 h, and performed RNA-Seq based transcriptome mapping followed by systematic bioinformatics analyses. Our gene expression analysis showed that several genes were differentially expressed in response to BaP and EE2 treatments in PCLS. Strong up-regulation of genes coding for the cytochrome P450 1a (Cyp1a) enzyme and the Ahr repressor (Ahrrb) was observed in BaP treated PCLS. EE2 treatment of liver slices strongly up-regulated genes coding for precursors of vitellogenin (Vtg) and eggshell zona pellucida (Zp) proteins. As expected, pathway enrichment and network analysis showed that the Ahr and estrogen receptor pathways are among the top affected by BaP and EE2 treatments, respectively. Interestingly, two genes coding for fibroblast growth factor 3 (Fgf3) and fibroblast growth factor 4 (Fgf4) were up-regulated by EE2 in this study. To our knowledge, the fgf3 and fgf4 genes have not previously been described in relation to estrogen signaling in fish liver, and these results suggest the modulation of the FGF signaling pathway by estrogens in fish. The signature expression profiles of top differentially expressed genes in response to the single compound (BaP or EE2) treatment were generally maintained in the expression responses to the equimolar binary mixtures. However, in the mixture-treated groups, BaP appeared to have anti-estrogenic effects as observed by lower number of differentially expressed putative EE2 responsive genes. Our in-depth quantitative analysis of changes in liver transcriptome in response to BaP and EE2, using PCLS tissue culture provides further mechanistic insights into effects of the compounds. Moreover, the analyses demonstrate the usefulness of PCLS in cod for omics experiments.

Fish & Shellfish Immunology, Volume 52, 2016, pp. 57-64

Lateolabrax japonicus is one of the main marine aquatic fish species, and is popularly cultured in East Asia due to its highly commercial value. In recent years, because of large-scale and intensive farming and seawater pollution, fish diseases keep breaking out. However, systematic study on L.japonicus immunogenetics is limited due to the deficiency of deep sequencing technologies and genome backgrounds. In this study, the widely analysis at the transcriptome level for L.japonicus that infected with Vibrio anguillarum was performed. In total, 334,388,688 high quality reads were obtained in six libraries (HK-VA, HK-PBS, LI-VA, LI-PBS, SP-VA and SP-PBS) and de novo assembled into 101,860 Unigenes with an average unigene length of 879 bp. Based on sequence similarity 30,142 unigenes (29.59%) were annotated in the public databases. Comparative analysis revealed, 1,202, 3034 and 3519 differentially expressed genes (DEGs) were identified in three comparisons (HK-PBS VS HK-VA, LI-PBS VS LI-VA and SP-PBS VS SP-VA). Enrichment and pathway analysis of the DEGs was also carried out to excavate the candidate genes related to immunity. In conclusion, this study identifies and evaluates dozen of potential immune related pathways and candidate genes, which are indispensable for padding genomic resources of L.japonicus, and would lay the foundation for further studying and illuminating the mechanism of host-pathogen interactions.

Ecotoxicology and Environmental Safety, Volume 182, 2019, Article 109464

Chlorpyrifos (CPF) is an environmental pollutant due to its high toxicity to aquatic animals. Because CPF was detected in aquatic environments in many countries, it has been widely concerned by researchers. Although the immunotoxicity of CPF to fish had been reported, the immunotoxicity mechanism is still not clear. Recently, transcriptome analysis has become a major method to study the toxic mechanism of pollutants in environmental toxicology. However, the immunotoxicity identification of CPF on fish had not been reported by transcriptome analysis. In the present study, we examined the effects of CPF on organismal system in the spleen of common carp by transcriptome analysis. We have successfully constructed a database of transcriptome analysis of carp spleens under exposure to CPF and found 773 differentially expressed genes (DEGs) (including 498 up-regulated DEGs and 275 down-regulated DEGs) and 4 branches (containing 33 known KEGG pathways). Some genes associated with the 4 pathways (Complement and coagulation cascades, PPAR signaling pathway, Fat digestion and absorption, and Collecting duct acid secretion) contained in organismal system were validated by quantitative real-time PCR and showed significant improvement compared with the control group. Our results indicated that exposure to CPF caused a change in the signal pathways of organismal system in carp spleens. The present study provides new insights into the immunotoxicity mechanism and risk assessment of CPF, as well as references for comparative medicine.