Transcriptomic characterization of adult zebrafish infected with Streptococcus agalactiae

Abstract

Streptococcus agalactiae is a major aquaculture pathogen infecting various saltwater and freshwater fish. To better understand the mechanism of the immune responses to S. agalactiae in wildtype zebrafish, the transcriptomic profiles of two organs containing mucosal-associated lymphoid tissues from S. agalactiae-infected and non-infected groups were obtained using RNA-seq techniques. In the intestines, 6735 and 12908 differently expressed genes (DEGs) were identified at 24 hpi and 48 hpi, respectively. Among 66 and 116 significantly enriched pathways, 15 and 21 pathways were involved in immune system or signal transduction at 24 hpi and 48 hpi, respectively. A number of genes involved in Toll-like receptor signaling pathway, RIG-I-like receptor signaling pathway, NOD-like receptor signaling pathway, T cell receptor signaling pathway, B cell receptor signaling pathway, Antigen processing and presentation, NF-kappa B signaling pathway and PI3K-Akt signaling pathway were significantly downregulated. In the skins, 3113 and 4467 DEGs were identified at 24 hpi and 48 hpi, respectively. Among 24 and 56 significantly enriched pathways, 4 and 13 pathways were involved in immune system or signal transduction at 24 hpi and 48 hpi, respectively. More immune-related signaling pathways including Leukocyte transendothelial migration, Cytokine-cytokine receptor interaction, PI3K-Akt signaling pathway, IL-17 signaling pathway, MAPK signaling pathway, TNF signaling pathway, Complement and coagulation cascades, Hematopoietic cell lineage and Jak-STAT signaling pathway were differently enriched for upregulated DEGs at 48 hpi, which were completely different from that in the intestines. Furthmore, comparative transcriptome analysis revealed that the downregulated 1618 genes and upregulated 1622 genes existed both at 24 hpi and 48 hpi for the intestine samples. In the skins, the downregulated 672 genes and upregulated 428 genes existed both at 24 hpi and 48 hpi. Three pathways related to immune processes were significantly enriched for downregulated DEGs both in the intestines and skins collected at 24 hpi and 48 hpi, which included Antigen processing and presentation, Intestinal immune network for IgA production and Hematopoietic cell lineage. Interaction network analysis of DEGs identified the main DEGs in the sub-network of complement and coagulation cascades both in the intestines and skins. Twenty of DEGs involved in complement and coagulation cascades were further validated by Real-time quantitative PCR. Altogether, the results obtained in this study will provide insight into the immune response of zebrafish against S. agalactiae XQ-1 infection in fatal conditions, and reveal the discrepant expression pattern of complement and coagulation cascades in the intestines and skins.