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

To determine and compare the toxic effects of Iranian heavy crude oil (IHCO) on the embryonic development of two fish species, we examined transcriptome profiles using RNA-seq. The assembled contigs were 66,070 unigenes in olive flounder embryos and 76,498 unigenes in spotted seabass embryos. In the differential gene expression (DEG) profiles, olive flounder embryos showed different up- and down-regulated patterns than spotted seabass embryos in response to fresh IHCO (FIHCO) and weathered IHCO (WIHCO). In this work, we categorized DEG profiles into six pathways: ribosome, oxidative phosphorylation, Parkinson's disease, Alzheimer's disease, Huntington's disease, and cardiac muscle contraction, validating the expression patterns of 13 DEGs using real-time quantitative RT-PCR. The expression of the CYP1A, CYP1B1, and CYP1C1 genes in spotted seabass embryos was higher than in olive flounder embryos, whereas genes related to cell processing, development, and the immune system showed the opposite trend. Orthologous gene cluster analysis showed that olive flounder embryos were sensitive (fold change of genes with cutoff P < 0.05) to both FIHCO and WIHCO, but spotted seabass embryos exhibited higher sensitivity to WIHCO than FIHCO, indicating that species-specific differences are likely to be reflected in population levels after oil spills. Overall, our study provides new insight on the different embryonic susceptibilities of two marine fish species to FIHCO and WIHCO and a better understanding of the underlying molecular mechanisms via RNA-seq and DEGs.