Transcriptome sequencing reveals improved ammonia nitrogen tolerance in Zebra II strain of the Manila clam Ruditapes philippinarum

Abstract

In this research, we identified genes associated with ammonia nitrogen (TAN) stress response and resistance in juveniles of the Zebra II strain and a wild population of the Manila clam Ruditapes philippinarum. Both groups were subjected to a 96 h acute toxicity test using TAN concentrations of 17.617 ± 0.634 and 16.670 ± 0.7 mg/l, respectively. We then collected samples, conducted transcriptome sequencing and screened the sequences for differentially expressed genes (DEGs) related to TAN stress response. We identified 2908 and 2861 DEGs in the Zebra II and wild clam groups, respectively, and the two groups had 626 DEGs in common. The verified DEGs had less of a detoxification effect in the wild population than that in the Zebra II group. Gene Ontology database analysis showed that Zebra II juveniles were mainly enriched in protein phosphorylation, purine nucleoside binding, and kinase activity, whereas the wild population juveniles were primarily enriched in oxidases activity, organic acid metabolic processes, and extracellular regions. Kyoto Encyclopedia of Genes and Genomes pathway analysis mainly highlighted aminoacyl tRNA biosynthesis in Zebra II juveniles and sphingolipid metabolism, FOXO signaling, biosynthesis of aminoacyl tRNA, and other pathways in the wild population. These results show that the toxic effect of TAN on the Manila clam is related to a variety of pathways, which are mainly related to immune response, inflammatory response, metabolic pathways, and nerve conduction. This study provides basic data and theoretical reference for revealing the molecular regulation mechanism of the improved TAN tolerance of Zebra II strain as compared with the wild population of Ruditapes philippinarum.