Comparative transcriptomic analysis of hepatopancreas reveals that more genes are involved in the exposure response of Vibrio parahaemolyticus PirAvp compared to PirBvp

Vibrio parahaemolyticus (VP-AHPND) is regarded as one of the main pathogens that caused acute hepatopancreatic necrosis disease (AHPND) in the Pacific white shrimp Litopenaeus vannamei. PirAvp and PirBvp toxin proteins are the main pathogenic proteins of AHPND in shrimp. Knowledge about the mechanism of shrimp response to PirAvp or PirBvp toxin is very helpful for developing new prevention and control strategy of AHPND in shrimp. In this study, the pathological sections showed that after 4 h treatment, significant pathological changes were observed in the PirBvp treated group, and no obvious pathological changes was found in PirAvp treated group. In order to learn the mechanism of shrimp response to PirAvp and PirBvp, comparative transcriptome was applied to analyze the different expressions of genes in the hepatopancreas of shrimp after treatment with PirAvp or PirBvp. A total of 9978 differentially expressed genes (DEGs) were identified between PirAvp or PirBvp-treated and PBS control shrimp, including 6616 DEGs in the PirAvp treated group and 3362 DEGs in the PirBvp treated group. There were 2263 DEGs that were commonly expressed, 4353 DEGs were only expressed in PirAvp VS PBS group and 1099 DEGs were uniquely expressed in PirBvp VS PBS group. Among these DEGs, the anti-apoptosis related pathways and immune response related genes significantly expressed in the commonly expressed DEGs of PirAvp VS PBS group and PirBvp VS PBS group, and small GTPase-mediated signaling and DNA metabolic process might relate to the host special reaction towards PirAvp and PirBvp exposure. The data suggested that the differential expression of these immune and metabolic-related genes in hepatopancreas might contribute to the pathogenicity variations of shrimp to VP-AHPND. The identified genes in this study will be useful for clarifying the response mechanism of shrimp toward different toxins of VP-AHPND and will further provide molecular basis for understanding the pathogenic mechanism of VP-AHPND. •Significant pathological changes were found in the PirBvp treated group instead of PirAvp treated group.•More genes are involved in the exposure response of Vibrio parahaemolyticus PirAvp compared to PirBvp.•The anti-apoptosis related and immune related genes were significantly expressed in the commonly expressed DEGs.•Small GTPase signaling pathway and DNA metabolic process might related to the exposure of PirAvp and PirBvp respectively.