Skip to main content
Log in

Characterization of transcriptional responses mediated by benzo[a]pyrene stress in a new marine fish model of goby, Mugilogobius chulae

  • Research Article
  • Published:
Genes & Genomics Aims and scope Submit manuscript

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.

This is a preview of subscription content, log in via an institution to check access.

Access this article

Price excludes VAT (USA)
Tax calculation will be finalised during checkout.

Instant access to the full article PDF.

Fig. 1
Fig. 2
Fig. 3
Fig. 4

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

    Article  CAS  PubMed  Google Scholar 

  • 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

    Article  CAS  PubMed  Google Scholar 

  • Anzenbacher P, Anzenbacherová E (2001) Cytochromes P450 and metabolism of xenobiotics. Cell Mol Life Sci 58(5):737–747

    Article  CAS  PubMed  Google Scholar 

  • 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

    Article  CAS  PubMed  Google Scholar 

  • 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

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  • 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

    CAS  Google Scholar 

  • Cai L, Huang R, Yu LJ, Li JJ (2016) Complete mitochondrial genome of Mugilogobius chulae (Perciformes: Gobiidae). Mitochondrial DNA A 27(6):4054–4055

    Article  CAS  Google Scholar 

  • 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

    Article  PubMed  Google Scholar 

  • Carman GM, Han GS (2006) Roles of phosphatidate phosphatase enzymes in lipid metabolism. Trends Biochem Sci 31(12):694–699

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  • 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

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  • 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

    Article  CAS  PubMed  Google Scholar 

  • 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

    Article  CAS  PubMed  Google Scholar 

  • 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

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  • 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

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  • 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

    Article  CAS  PubMed  Google Scholar 

  • 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

    Article  PubMed  Google Scholar 

  • 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

    Article  CAS  PubMed  Google Scholar 

  • Higgins CF (2007) Multiple molecular mechanisms for multidrug resistance transporters. Nature 446(12):749–757

    Article  CAS  PubMed  Google Scholar 

  • 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

    Article  PubMed  Google Scholar 

  • 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

    Article  PubMed  CAS  Google Scholar 

  • 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

    Article  CAS  PubMed  Google Scholar 

  • 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

    Google Scholar 

  • Jakoby WB (1978) The glutathione S-transferases: a group of multifunctional detoxification proteins. Adv Enzymol Relat Areas Mol Biol 46:383–414

    CAS  PubMed  Google Scholar 

  • 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

    Article  CAS  PubMed  Google Scholar 

  • 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

    Article  CAS  PubMed  Google Scholar 

  • Larson HK (2001) A revision of the gobiid fish genus Mugilogobius (Teleostei: Gobioidei), and its systematic placement. Western Australian Museum, Perth

    Google Scholar 

  • 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

    Article  CAS  PubMed  Google Scholar 

  • Li B, Dewey CN (2011) RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinform 12:323

    Article  CAS  Google Scholar 

  • 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

    Article  CAS  PubMed  Google Scholar 

  • 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

    Google Scholar 

  • 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

    Article  CAS  PubMed  Google Scholar 

  • 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

    Article  CAS  PubMed  Google Scholar 

  • 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

    Article  CAS  PubMed  Google Scholar 

  • 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

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  • 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

    CAS  Google Scholar 

  • 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

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  • 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

    Article  CAS  PubMed  Google Scholar 

  • 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

    Article  CAS  PubMed  Google Scholar 

  • 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

    Article  CAS  PubMed  Google Scholar 

  • Pascual F, Carman GM (2013) Phosphatidate phosphatase, a key regulator of lipid homeostasis. Biochim Biophys Acta 1831(3):514–522

    Article  CAS  PubMed  Google Scholar 

  • 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

    Article  PubMed  Google Scholar 

  • 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

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  • Rickborn AJ, Buston PM (2015) Life-history transitions of the coral-reef fish Elacatinus lori. J Fish Biol 86(2):637–650

    Article  CAS  PubMed  Google Scholar 

  • 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

    CAS  Google Scholar 

  • 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

    Article  PubMed  PubMed Central  Google Scholar 

  • 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

    Article  CAS  PubMed  Google Scholar 

  • 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

    Article  CAS  PubMed  Google Scholar 

  • 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

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  • 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

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  • 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

    CAS  Google Scholar 

  • 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

    Article  PubMed  CAS  Google Scholar 

  • 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

    Article  CAS  PubMed  Google Scholar 

  • 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

    Article  CAS  PubMed  Google Scholar 

  • 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

    Article  CAS  PubMed  Google Scholar 

  • 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

    Article  CAS  PubMed  Google Scholar 

  • 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

    Article  CAS  PubMed  Google Scholar 

  • 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

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  • 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

    Article  CAS  PubMed  Google Scholar 

  • 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

    Article  CAS  PubMed  Google Scholar 

  • Wu G, Fang YZ, Yang S, Lupton JR, Turner ND (2004) Glutathione metabolism and its implications for health. J Nutr 134(3):489–492

    Article  CAS  PubMed  Google Scholar 

  • 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

    Article  CAS  PubMed  Google Scholar 

  • 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

    Article  PubMed  PubMed Central  Google Scholar 

  • 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

    CAS  PubMed  Google Scholar 

  • 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

    Article  PubMed Central  CAS  Google Scholar 

  • 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

    CAS  Google Scholar 

Download references

Acknowledgements

This work was supported by the National Key Technologies R & D Program of China (Grant No. 2015BAI09B05).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Ren Huang.

Ethics declarations

Conflict of interest

The authors would like to declare no conflict of interest.

Rights and permissions

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

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

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1007/s13258-018-0743-8

Keywords

Navigation