RNA-Seq-based transcriptome profiling and expression of 16 cytochrome P450 genes in the benzo[α]pyrene-exposed estuarine copepod Eurytemora affinis

The calanoid copepod Eurytemora affinis is one of the most abundant estuarine species and is considered to be an ideal candidate species for ecotoxicological research. An RNA-Seq-based transcriptome was developed from whole bodies of this species. Among 142,442 contigs of the de novo assembly by Trinity, 48,480 open reading frame (ORF) contigs were found using TransDecoder. A total of 17,762 genes were identified by BLAST analysis, which covers about 75% of the annotated genes in the E. affinis genome. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis showed that most annotated sequences were related to metabolism pathways, including xenobiotic biodegradation. Using transcriptome data, we identified putative transcripts related to xenobiotic processing genes including phase I enzymes, phase II enzymes, transporters, and transcription factors. To understand the CYP-mediated detoxification metabolism of xenobiotics, we measured the transcriptional levels of 16 CYPs (within full sequences) of E. affinis in response to benzo[α]pyrene (B[α]P). Most Ea-CYP genes were significantly down- and/or up-regulated (P < 0.05) in response to B[α]P, suggesting that Ea-CYP genes are likely involved in detoxification (mainly in biotransformation of xenobiotics) with particular genes, demonstrating significant upregulation or downregulation compared to others, as shown in other copepod model species (e.g. Tigriopus japonicus and Paracyclopina nana). This study will provide insight into the potential role of E. affinis in response to various toxic or xenobiotic chemicals in the marine environment.