Abstract
Benzo[a]pyrene (BaP) is one of the most studied targets among polycyclic aromatic hydrocarbons (PAHs). Because of the complexity of the toxicity mechanism in BaP, little is known about the molecular mechanism at the level of transcription of BaP in marine fishes. The primary objective of this study was to investigate the molecular basis of the effects of BaP on marine fish, using Mugilogobius chulae (Smith 1932) as the model. A closed colony of M. chulae was used for the BaP toxicity test. Two fish liver samples per replicate from each group were excised and blended into one sample by pooling an equal amount of liver tissue. Total RNA of all samples was extracted separately. Equal quantities of total RNA from the three replicates of the two groups were pooled for sequencing. The sequencing cDNA libraries were sequenced using Illumina HiSeq 2000 system. Differentially expressed genes were detected with the DEGSeq R package. In total, 52,364,032 and 53,771,748 clean nucleotide reads were obtained in the control and BaP-exposed libraries, respectively, with N50 lengths of 1277 and 1288 bp, respectively. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analyses revealed a significant enrichment of genes related to detoxification, transportation, and lipid metabolism. We also identified, for the first time, an association between endoplasmic reticulum dysfunction and lipid metabolism resulting from BaP exposure. Using quantitative real-time PCR, some effective molecular biomarkers for monitoring of BaP-polluted seawater were identified. The results demonstrate that BaP enhanced the expression of genes involved in detoxification in M. chulae and inhibited that of genes related to lipid metabolism, possibly by suppressing the expression of numerous ER-related genes involved in fat digestion and absorption.
Similar content being viewed by others
References
Ahn J, Lee H, Jung CH, Ha T (2012) Lycopene inhibits hepatic steatosis via microRNA-21-induced downregulation of fatty acid-binding protein 7 in mice fed a high- fat diet. Mol Nutr Food Res 56(11):1665–1674
Andersen Ø, Frantzen M, Rosland M, Timmerhaus G, Skugor A, Krasnov A (2015) Effects of crude oil exposure and elevated temperature on the liver transcriptome of polar cod (Boreogadus saida). Aquat Toxicol 165:9–18
Anzenbacher P, Anzenbacherová E (2001) Cytochromes P450 and metabolism of xenobiotics. Cell Mol Life Sci 58(5):737–747
Barat A, Sahoo PK, Kumar R, Goel C, Singh AK (2016) Transcriptional response to heat shock in liver of snow trout (Schizothorax richardsonii)-a vulnerable Himalayan Cyprinid fish. Funct Integr Genomics 16(2):203–213
Brindley DN, Pilquil C, Sariahmetoglu M, Reue K (2009) Phosphatidate degradation: phosphatidate phosphatases (lipins) and lipid phosphate phosphatases. Biochim Biophys Acta 1791(9), 956–961
Cai L, Chen XQ, Zheng WQ, Li JJ (2015) Isolation and characterization of polymorphic microsatellitemarkers in Mugilogobius chulae. Acta Lab Anim Sci Sin 23(1):57–62
Cai L, Huang R, Yu LJ, Li JJ (2016) Complete mitochondrial genome of Mugilogobius chulae (Perciformes: Gobiidae). Mitochondrial DNA A 27(6):4054–4055
Cai L, Li JJ, Yu LJ, Wei YZ, Miao ZY, Chen ML, Huang R (2018) De novo transcriptome assembly of the new marine fish model of goby, Mugilogobius chulae. Mar Genomics. https://doi.org/10.1016/j.margen.2018.02.001
Carman GM, Han GS (2006) Roles of phosphatidate phosphatase enzymes in lipid metabolism. Trends Biochem Sci 31(12):694–699
Carvalho RN, Bopp SK, Lettieri T (2011) Transcriptomics responses in marine diatom Thalassiosira pseudonana exposed to the polycyclic aromatic hydrocarbon benzo[a]pyrene. PLoS ONE 6(11):e26985
Chen QL, Luo Z, Huang C, Pan YX, Wu K (2016) De novo characterization of the liver transcriptome of javelin goby Synechogobius hasta and analysis of its transcriptomic profile following waterborne copper exposure. Fish Physiol Biochem 42(3):979–994. https://doi.org/10.1007/s10695-015-0190-2
Denissenko MF, Pao A, Tang MS, Pfeifer GP (1996) Preferential formation of benzo [a] pyrene adducts at lung cancer mutational hotspots in P53. Science 274:430–432
Eastmond PJ, Quettier AL, Kroon JT, Craddock C, Adams N, Slabas AR (2010) PHOSPHATIDIC ACID PHOSPHOHYDROLASE1 and 2 regulate phospholipid synthesis at the endoplasmic reticulum in Arabidopsis. Plant Cell 22(8):2796–2811
Eisler R (1987) Polycyclic aromatic hydrocarbon hazards to fish, wildlife, and invertebrates: a synoptic review (No. PB-87-189825/XAB; BIOLOGICAL-85 (1.11)). Patuxent Wildlife Research Center, Laurel. https://www.pwrc.usgs.gov/eisler/CHR_11_PAHs.pdf. Accessed 14 June 2018
Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, Adiconis X, Fan L, Raychowdhury R, Zeng Q, Chen Z (2011) Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol 29(7):644–652
Guo Z, Gao N, Wu Y, Zhang L (2017) The simultaneous uptake of dietary and waterborne Cd in gastrointestinal tracts of marine yellowstripe goby Mugilogobius chulae. Environ Pollut 223:31–41
Harney E, Dubief B, Boudry P, Basuyaux O, Schilhabel MB, Huchette S, Paillard C, Nunes FL (2016) De novo assembly and annotation of the European abalone Haliotis tuberculata transcriptome. Mar Genomics 28:11–16
Hassanin AA, Kaminishi Y, Funahashi A, Itakura T (2012) Cytochrome P450 1C1 complementary DNA cloning, sequence analysis and constitutive expression induced by benzo-a-pyrene in Nile tilapia (Oreochromis niloticus). Aquat Toxicol 109:17–24
Higgins CF (2007) Multiple molecular mechanisms for multidrug resistance transporters. Nature 446(12):749–757
Hu F, Pan L, Cai Y, Liu T, Jin Q (2015) Deep sequencing of the scallop Chlamys farreri transcriptome response to tetrabromobisphenol A (TBBPA) stress. Mar Genomics 19:31–38
Huang Q, Dong S, Fang C, Wu X, Ye T, Lin Y (2012) Deep sequencing-based transcriptome profiling analysis of Oryzias melastigma exposed to PFOS. Aquat Toxicol 120:54–58
Huang L, Zuo Z, Zhang Y, Wu M, Lin JJ, Wang C (2014) Use of toxicogenomics to predict the potential toxic effect of benzo(a)pyrene on zebrafish embryos: ocular developmental toxicity. Chemosphere 108:55–61
International Agency for Research on Cancer (2017) Agents classified by the IARC monographs. http://monographs.iarc.fr/ENG/Classification/latest_classif.php/. Accessed 24 Feb 2018
Iseli C, Jongeneel CV, Bucher P (1999) ESTScan: a program for detecting, evaluating, and reconstructing potential coding regions in EST sequences. In ISMB 99:138–148
Jakoby WB (1978) The glutathione S-transferases: a group of multifunctional detoxification proteins. Adv Enzymol Relat Areas Mol Biol 46:383–414
Kim KH, Jahan SA, Kabir E, Brown RJ (2013) A review of airborne polycyclic aromatic hydrocarbons (PAHs) and their human health effects. Environ Int 60:71–80
Labib S, Guo CH, Williams A, Yauk CL, White PA, Halappanavar S (2013) Toxicogenomic outcomes predictive of forestomach carcinogenesis following exposure to benzo(a)pyrene: relevance to human cancer risk. Toxicol Appl Pharmacol 273(2):269–280
Larson HK (2001) A revision of the gobiid fish genus Mugilogobius (Teleostei: Gobioidei), and its systematic placement. Western Australian Museum, Perth
Leslie EM, Deeley RG, Cole SP (2001) Toxicological relevance of the multidrug resistance protein 1, MRP1 (ABCC1) and related transporters. Toxicology 167(1):3–23
Li B, Dewey CN (2011) RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinform 12:323
Li Q, Zhang X, Yan C (2010) Polycyclic aromatic hydrocarbon contamination of recent sediments and marine organisms from Xiamen Bay, China. Arch Environ Contam Toxicol 58(3):711–721
Li JJ, Lin ZT, Chen XQ, Huang R (2014) Single and joint toxicity of four heavy metal ions on Mugilogobius chulae. Mar Environ Sci 33(2):236–241
Li C, Ling Q, Ge C, Ye Z, Han X (2015) Transcriptome characterization and SSR discovery in large-scale loach Paramisgurnus dabryanus (Cobitidae, Cypriniformes). Gene 557(2):201–208
Liu T, Pan L, Jin Q, Cai Y (2015) Differential gene expression analysis of benzo(a) pyrene toxicity in the clam, Ruditapes philippinarum. Ecotoxicol Environ Saf 115:126–136
Livak KJ, Schmittgen TD (2001) Analysis of relative gene expression data using real-time quantitative PCR and the 2−∆∆CT method. Methods 25(4):402–408
Long AS, Lemieux CL, Arlt VM, White PA (2016) Tissue-specific in vivo genetic toxicity of nine polycyclic aromatic hydrocarbons assessed using the Muta™ Mouse transgenic rodent assay. Toxicol Appl Pharmacol 290:31–42
Lu TT, Lin Q, Ke CL, Sun RX (2012) Polycyclic aromatic hydrocarbons and risk assessment in the surface sediments from Lingdingyang, Pearl River Estuary. J Fish Sci China 19(2):336–347
Marioni JC, Mason CE, Mane SM, Stephens M, Gilad Y (2008) RNA-seq: an assessment of technical reproducibility and comparison with gene expression arrays. Genome Res 18(9):1509–1517
Men B, He M, Tan L, Lin C, Quan X (2009) Distributions of polycyclic aromatic hydrocarbons in the Daliao River Estuary of Liaodong Bay, Bohai Sea (China). Mar Pollut Bull 58(6):818–826
Miao LH, Lin Y, Pan WJ, Huang X, Ge XP, Zhou QL, Liu B, Ren MC, Zhang WX, Liang HL, Yu H (2018) Comparative transcriptome analysis reveals the gene expression profiling in bighead carp (Aristichthys nobilis) in response to acute nitrite toxicity. Fish Shellfish Immunol 79:244–255
Miyazaki H, Sawada T, Kiyohira M, Yu Z, Nakamura K, Yasumoto Y, Kagawa Y, Ebrahimi M, Islam A, Sharifi K, Kawamura S (2014) Fatty acid binding protein 7 regulates phagocytosis and cytokine production in kupffer cells during liver injury. Am J Pathol 184(9):2505–2515
Pascual F, Carman GM (2013) Phosphatidate phosphatase, a key regulator of lipid homeostasis. Biochim Biophys Acta 1831(3):514–522
Qian B, Xue L (2016) Liver transcriptome sequencing and de novo annotation of the large yellow croaker (Larimichthy crocea) under heat and cold stress. Mar Genomics 25:95–102
Regnault C, Worms IA, Oger-Desfeux C, MelodeLima C, Veyrenc S, Bayle ML, Combourieu B, Bonin A, Renaud J, Raveton M, Reynaud S (2014) Impaired liver function in Xenopus tropicalis exposed to benzo[a]pyrene: transcriptomic and metabolic evidence. BMC Genomics 15:666
Rickborn AJ, Buston PM (2015) Life-history transitions of the coral-reef fish Elacatinus lori. J Fish Biol 86(2):637–650
Seemann F, Peterson DR, Witten PE, Guo BS, Shanthanagouda AH, Rui RY, Zhang G, Au DW (2015) Insight into the transgenerational effect of benzo[a]pyrene on bone formation in a teleost fish (Oryzias latipes). Comp Biochem Physiol C 178:60–67
Sinaei M, Mashinchian A (2014) Polycyclic aromatic hydrocarbons in the coastal sea water, the surface sediment and Mudskipper Boleophthalmus dussumieri from coastal areas of the Persian Gulf: source investigation, composition pattern and spatial distribution. J Environ Health Sci Eng 12:59
Squadrone S, Favaro L, Abete MC, Vivaldi B, Prearo M (2014) Polycyclic aromatic hydrocarbon levels in European catfish from the upper Po River basin. Environ Monit Assess 186(4):2313–2320
Su X, Tan QS, Parikh BH, Tan A, Mehta MN, Wey YS, Tun SB, Li LJ, Han XY, Wong TY, Hunziker W (2016) Characterization of fatty acid binding protein 7 (FABP7) in the murine retina characterization of FABP7 in the murine retina. Invest Ophthalmol Vis Sci 57(7):3397–3408
Takaoka N, Takayama T, Teratani T, Sugiyama T, Mugiya S, Ozono S (2011) Analysis of the regulation of fatty acid binding protein 7 expression in human renal carcinoma cell lines. BMC Mol Biol 12:31
Takaoka N, Takayama T, Ozono S (2017) Functional analysis of fatty acid binding protein 7 and its effect on fatty acid of renal cell carcinoma cell lines. BMC Cancer 17:192
Tian Y, Zheng T, Wang X, Luo Y (2003) Contamination characteristics of polycyclic aromatic hydrocarbons in Maluan Bay mariculture area of Xiamen. Mar Environ Sci 22(1):29–33
Tian S, Pan L, Sun X (2013) An investigation of endocrine disrupting effects and toxic mechanisms modulated by benzo[a]pyrene in female scallop Chlamys farreri. Aquat Toxicol 144:162–171
Tong C, Zhang C, Zhang R, Zhao K (2015) Transcriptome profiling analysis of naked carp (Gymnocypris przewalskii) provides insights into the immune-related genes in highland fish. Fish Shellfish Immunol 46(2):366–377
Uno S, Dalton TP, Dragin N, Curran CP, Derkenne S, Miller ML, Shertzer HG, Gonzalez FJ, Nebert DW (2006) Oral benzo[a]pyrene in Cyp1 knockout mouse lines: CYP1A1 important in detoxication, CYP1B1 metabolism required for immune damage independent of total-body burden and clearance rate. Mol Pharmacol 69(4):1103–1114
Wakabayashi K, Tamura A, Saito H, Onishi Y, Ishikawa T (2006) Human ABC transporter ABCG2 in xenobiotic protection and redox biology. Drug Metab Rev 38(3):371–391
Wan Y, Jin X, Hu J, Jin F (2007) Trophic dilution of polycyclic aromatic hydrocarbons (PAHs) in a marine food web from Bohai Bay, North China. Environ Sci Technol 41(9):3109–3114
Wang L, Feng Z, Wang X, Wang X, Zhang X (2009) DEGseq: an R package for identifying differentially expressed genes from RNA-seq data. Bioinformatics 26(1):136–138
Wang Z, Liu Z, Xu K, Mayer LM, Zhang Z, Kolker AS, Wu W (2014) Concentrations and sources of polycyclic aromatic hydrocarbons in surface coastal sediments of the northern Gulf of Mexico. Geochem Trans 15:2
Weisiger RA (2002) Cytosolic fatty acid binding proteins catalyze two distinct steps in intracellular transport of their ligands. Mol Cell Biochem 239(1):35–43
Wen J, Pan L (2016) Short-term exposure to benzo[a]pyrene causes oxidative damage and affects haemolymph steroid levels in female crab Portunus trituberculatus. Environ Pollut 208:486–494
Wu G, Fang YZ, Yang S, Lupton JR, Turner ND (2004) Glutathione metabolism and its implications for health. J Nutr 134(3):489–492
Xu C, Li CYT, Kong ANT (2005) Induction of phase I, II and III drug metabolism/transport by xenobiotics. Arch Pharm Res 28(3):249–268
Xu EG, Mager EM, Grosell M, Hazard ES, Hardiman G, Schlenk D (2017) Novel transcriptome assembly and comparative toxicity pathway analysis in mahi-mahi (Coryphaena hippurus) embryos and larvae exposed to Deepwater Horizon oil. Sci Rep 7:44546
Yang YH, Miranda CL, Henderson MC, Wang-Buhler JL, Buhler DR (2000) Heterologous expression of Cyp2k1 and identification of the expressed protein (bv-Cyp2k1) as lauric acid (ω-1)-hydroxylase and Aflatoxin B1exo-epoxidase. Drug Metab Dispos 28(11):1279–1283
Ye R, Holland WL, Gordillo R, Wang M, Wang QA, Shao M, Morley TS, Gupta RK, Stahl A, Scherer PE (2014) Adiponectin is essential for lipid homeostasis and survival under insulin deficiency and promotes β-cell regeneration. Elife 3:e03851
Zhang L, Qiu L, Wu H, Liu X, You L, Pei D, Chen L, Wang Q, Zhao J (2012) Expression profiles of seven glutathione S-transferase (GST) genes from Venerupis philippinarum exposed to heavy metals and benzo[a]pyrene. Comp Biochem Physiol C 155(3):517–527
Acknowledgements
This work was supported by the National Key Technologies R & D Program of China (Grant No. 2015BAI09B05).
Author information
Authors and Affiliations
Corresponding author
Ethics declarations
Conflict of interest
The authors would like to declare no conflict of interest.
Rights and permissions
About this article
Cite this article
Cai, L., Li, J., Yu, L. et al. Characterization of transcriptional responses mediated by benzo[a]pyrene stress in a new marine fish model of goby, Mugilogobius chulae. Genes Genom 41, 113–123 (2019). https://doi.org/10.1007/s13258-018-0743-8
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s13258-018-0743-8