Co-Expression Network Analysis of Spleen Transcriptome in Rock Bream (Oplegnathus fasciatus) Naturally Infected with Rock Bream Iridovirus (RBIV)

Rock bream iridovirus (RBIV) is a notorious agent that causes high mortality in aquaculture of rock bream (Oplegnathus fasciatus). Despite severity of this virus, no transcriptomic studies on RBIV-infected rock bream that can provide fundamental information on protective mechanism against the virus have been reported so far. This study aimed to investigate physiological mechanisms between host and RBIV through transcriptomic changes in the spleen based on RNA-seq. Depending on infection intensity and sampling time point, fish were divided into five groups: uninfected healthy fish at week 0 as control (0C), heavy infected fish at week 0 (0H), heavy mixed RBIV and bacterial infected fish at week 0 (0MH), uninfected healthy fish at week 3 (3C), and light infected fish at week 3 (3L). We explored clusters from 35,861 genes with Fragments Per Kilo-base of exon per Million mapped fragments (FPKM) values of 0.01 or more through signed co-expression network analysis using WGCNA package. Nine of 22 modules were highly correlated with viral infection (|gene significance (GS) vs. module membership (MM) |> 0.5, p-value < 0.05). Expression patterns in selected modules were divided into two: heavy infected (0H and 0MH) and control and light-infected groups (0C, 3C, and 3L). In functional analysis, genes in two positive modules (5448 unigenes) were enriched in cell cycle, DNA replication, transcription, and translation, and increased glycolysis activity. Seven negative modules (3517 unigenes) built in this study showed significant decreases in the expression of genes in lymphocyte-mediated immune system, antigen presentation, and platelet activation, whereas there was significant increased expression of endogenous apoptosis-related genes. These changes lead to RBIV proliferation and failure of host defense, and suggests the importance of blood cells such as thrombocytes and B cells in rock bream in RBIV infection. Interestingly, a hub gene, pre-mRNA processing factor 19 (PRPF19) showing high connectivity (kME), and expression of this gene using qRT-PCR was increased in rock bream blood cells shortly after RBIV was added. It might be a potential biomarker for diagnosis and vaccine studies in rock bream against RBIV. This transcriptome approach and our findings provide new insight into the understanding of global rock bream-RBIV interactions including immune and pathogenesis mechanisms.