Inoculation of Freund's adjuvant in European eel (Anguilla anguilla) revealed key KEGG pathways and DEGs of host anti-Edwardsiella anguillarum infection

Abstract

Freund's complete (FCA) and incomplete adjuvants (FIA), generally applied in subunit fishery vaccine, have not been explored on the molecular mechanism of the nonspecific immune enhancement. In this study, we examined the RNA-seq in the spleen of European eel (Anguilla anguilla) inoculated with FCA and FIA (FCIA group) to elucidate the key KEGG pathways and differential expressed genes (DEGs) in the process of Edwardsiella anguillarum infection and A. anguilla anti-E. anguillarum infection using genome-wide transcriptome. After eels were challenged by E. anguillarum at 28 d post the first inoculation (dpi), compared to the control uninfected eels (Con group), the control infected eels (Con_inf group) showed severe pathological changes in the liver, kidney and spleen, although infected eels post the inoculation of FCIA (FCIA_inf group) also formed slight bleeding. Compared to the FCIA_inf group, there was more than 10 times colony forming unit (cfu) in the Con_inf group per 100 μg spleen, kidney or blood, and the relative percent survival (RPS) of eels was 44.4% in FCIA_inf vs Con_inf. Compared to the Con group, the SOD activity in the FCIA group increased significantly in the liver and spleen. Using high-throughput transcriptomics, DEGs were identified and 29 genes were verified using fluorescence real-time polymerase chain reaction (qRT-PCR). The result of DEGs clustering showed 9 samples in 3 groups of Con, FCIA and FCIA_inf were similar, contrast to distinct differences of 3 samples in the Con_inf group. We found 3795 up and 3548 down regulated DEGs in the compare of FCIA_inf vs Con_inf, of which 5 enriched KEGG pathways of “Lysosome”, “Autophagy”, “Apoptosis”, “C-type lectin receptor signaling” and “Insulin signaling” were ascertained, and 26 of 30 top GO terms in the compare were significantly enriched. Finally, protein-protein interactions between the DEGs of the 5 KEGG pathways and other DEGs were explored using Cytoscape 3.9.1. The compare of FCIA_inf vs Con_inf showed 110 DEGs from the 5 pathways and 718 DEGs from other pathways formed total of 9747° in a network, of which 9 hub DEGs play vital roles in anti-infection or apoptosis. Together, the interaction networks revealed that 9 DEGs involved in the 5 pathways underlies the key process of A. anguilla anti-E. anguillarum infection or host cell apoptosis.