Transcriptome analysis reveals molecular mechanisms of lymphocystis formation caused by lymphocystis disease virus infection in flounder (Paralichthys olivaceus)

Lymphocystis disease is frequently prevalent and transmissible in various teleost species worldwide due to lymphocystis disease virus (LCDV) infection, causing unsightly growths of benign lymphocystis nodules in fish and resulting in huge economic losses to aquaculture industry. However, the molecular mechanism of lymphocystis formation is unclear. In this study, LCDV was firstly detected in naturally infected flounder (Paralichthys olivaceus) by PCR, histopathological, and immunological techniques. To further understand lymphocystis formation, transcriptome sequencing of skin nodule tissue was performed by using healthy flounder skin as a control. In total, RNA-seq produced 99.36%-99.71% clean reads of raw reads, of which 91.11%-92.89% reads were successfully matched to the flounder genome. The transcriptome data showed good reproducibility between samples, with 3781 up-regulated and 2280 down-regulated differentially expressed genes. GSEA analysis revealed activation of Wnt signaling pathway, Hedgehog signaling pathway, Cell cycle, and Basal cell carcinoma associated with nodule formation. These pathways were analyzed to interact with multiple viral infection and tumor formation pathways. Heat map and protein interaction analysis revealed that these pathways regulated the expression of cell cycle-related genes such as ccnd1 and ccnd2 through key genes including ctnnb1, lef1, tcf3, gli2, and gli3 to promote cell proliferation. Additionally, cGMP-PKG signaling pathway, Calcium signaling pathway, ECM-receptor interaction, and Cytokine-cytokine receptor interaction associated with nodule formation were significantly down-regulated. Among these pathways, tnfsf12, tnfrsf1a, and tnfrsf19, associated with pro-apoptosis, and vdac2, which promotes viral replication by inhibiting apoptosis, were significantly up-regulated. Visual analysis revealed significant down-regulation of cytc, which expresses the pro-apoptotic protein cytochrome C, as well as phb and phb2, which have anti-tumor activity, however, casp3 was significantly up-regulated. Moreover, bcl9, bcl11a, and bcl-xl, which promote cell proliferation and inhibit apoptosis, were significantly upregulated, as were fgfr1, fgfr2, and fgfr3, which are related to tumor formation. Furthermore, RNA-seq data were validated by qRT-PCR, and LCDV copy numbers and expression patterns of focused genes in various tissues were also investigated. These results clarified the pathways and differentially expressed genes associated with lymphocystis nodule development caused by LCDV infection in flounder for the first time, providing a new breakthrough in molecular mechanisms of lymphocystis formation in fish.

New Insights into Lymphocystis Disease Virus Genome Diversity

Lymphocystis disease viruses (LCDVs) are viruses that infect bony fish which has been found in different locations across the globe. Four virus species have been classified by the International Committee on Taxonomy of Viruses (ICTV), despite remarkable discrepancies in genome size. Whole genome sequencing and phylogenetic analysis of LCDVs from wild fish from the North Sea and partial sequences from gilthead sea bream of an aquafarm located in the Aegean Sea in Turkey confirm that the LCDV1 genome at 100 kb is approximately half the size of the genomes of LCDV2-4. Since the fish species, of which LCDV1 was isolated, differ taxonomically at the order level, co-speciation can be excluded as the driver of the adaptation of the genome of this nucleocytoplasmic large DNA virus, but may represent an adaptation to the lifestyle of this demersal fish in the northeast Atlantic.

Initial Evidence That Gilthead Seabream (Sparus aurata L.) Is a Host for Lymphocystis Disease Virus Genotype I

Simple Summary

Nodular lesions were observed on the skin and fins of 95% of one and a half million juvenile gilthead seabreams cultured in Egypt, shortly after importation from Europe. We undertook a study to describe the clinical disease course, identify the causative agent, and investigate the origin of the causative agent. Preliminary diagnosis based on gross lesions and postmortem examination suggested lymphocystis disease caused by lymphocystis disease virus, Lymphocystivirus, Iridoviridae. Histopathological and ultrastructural pictures were typical of lymphocystis disease virus infections. Polymerase Chain Reaction followed by sequencing and phylogenetic analysis of the major capsid protein gene demonstrated the presence of lymphocystis disease virus genotype I, originally associated with lymphocystis disease in Northern European countries, with 99.7% and 100% nucleotide and deduced amino acid identity values, respectively. Lymphocystis disease virus genotype I has never been reported in this species or in the region. Regardless of whether it has maintained a previously undetected state of endemicity in Egypt or was introduced through importation or contamination of ship ballast water, the findings of this study add to existing knowledge about the lymphocystis disease’s ecology, and lymphocystis disease virus genotypes and their host range.

Abstract

Marine and brackish water aquacultures are rapidly expanding in the Mediterranean basin. In this context, Egypt recently received a shipment of a 1.5 million juvenile gilthead seabream (Sparus aurata L.) from European Mediterranean facility. Within a few weeks of their arrival, 95% of the imported fish developed nodules on their skin and fins that lasted for several months. This study was undertaken to describe the clinical disease course, to identify the causative agent, and to investigate its origin. Preliminary diagnosis based on gross lesions and postmortem examination suggested lymphocystis disease (LCD), caused by the lymphocystis disease virus (LCDV; genus Lymphocystivirus, family Iridoviridae). Histopathological and ultrastructural features were typical of LCDV infections. PCR followed by sequencing and phylogenetic analysis of a 306-bp fragment of the major capsid protein (MCP) gene demonstrated the presence of LCDV genotype I, originally associated with LCD in Northern European countries, with 99.7% and 100% nucleotide and deduced amino acid identity values, respectively. LCDV genotype I has neither been reported in this species nor in the region. Regardless of the source of infection, findings of this study add to existing knowledge about the ecology of LCDV genotype I and its host range.

Immune gene expression in gilthead seabream (Sparus aurata) after Lymphocystis disease virus (LCDV-Sa) challenge resulting in asymptomatic infection

AIMS

To determine the immune gene expression response of gilthead seabream (Sparus aurata) that is experimentally infected with the lymphocystivirus LCDV-Sa.

METHODS AND RESULTS

Viral DNA and transcripts were detected by qPCR in all samples from fish injected with LCDV-Sa, demonstrating that the virus establish a systemic and asymptomatic infection. The expression of 23 immune-related genes was also analysed by RT-qPCR in the head kidney (HK) and intestine at several times post-infection (dpi). In HK, the expression of five type I interferon (IFN)-related genes (ifn, irf3, mx2, mx3 and isg15), il10 and ck10 was upregulated at 1-3 dpi, while genes related to the inflammation process (tnfα, il1ß, il6, casp1) were not differentially expressed or even downregulated. The expression profile in the intestine was different regarding type I INF-related genes. An upregulated c3 and ighm expression was observed in both HK and intestine at 3-8 dpi. Finally, the transcription of nccrp1 and mhcIIα was induced in HK, whereas tcrβ expression was downregulated in both organs.

CONCLUSIONS

LCDV-Sa seems to trigger an immune response in gilthead seabream characterized by a partial activation of type I IFN system and a lack of systemic inflammatory response which may be related to viral persistence.

SIGNIFICANCE AND IMPACT OF THE STUDY

The immune response observed in gilthead seabream infected by LCDV-Sa could be implicated in the establishment of an asymptomatic persistent infection.

Voltage-Dependent Anion Channel Protein 2 (VDAC2) and Receptor of Activated Protein C Kinase 1 (RACK1) Act as Functional Receptors for Lymphocystis Disease Virus Infection

In previous research, a 27.8-kDa protein in flounder Paralichthys olivaceus gill (FG) cells was identified as a putative cellular receptor (27.8R), which mediated lymphocystis disease virus (LCDV) infection via interaction with a 32-kDa viral attachment protein (VAP) of LCDV, and monoclonal antibodies (MAbs) against 27.8R and 32-kDa VAP were developed. In this study, the 27.8R was identified as voltage-dependent anion channel protein 2 (VDAC2) and receptor of activated protein C kinase 1 (RACK1) of flounder. Recombinant VDAC2 (rVDAC2) and RACK1 (rRACK1) were obtained by prokaryotic expression, and rabbit anti-VDAC2/RACK1 polyclonal antibodies were prepared. The rVDAC2, rRACK1, and 27.8-kDa proteins in FG cells were recognized by anti-27.8R MAbs and anti-VDAC2/RACK1 polyclonal antibodies simultaneously. Preincubation of FG cells with anti-VDAC2/RACK1 polyclonal antibodies significantly decreased the percentages of LCDV-infected cells and LCDV copy numbers, blocked virus infection, and delayed the development of cytopathic effect. The mRNA expressions of VDAC2 and RACK1 in FG cells were upregulated to maximum levels 12 h and 48 h after LCDV infection, respectively. VDAC2/RACK1 knockdown through short interfering RNA (siRNA) significantly reduced VDAC2/RACK1 expression and LCDV copy numbers in FG cells compared with negative controls, while VDAC2/RACK1 expression on LCDV-nonpermissive epithelial papillosum cells (EPCs) conferred susceptibility to LCDV infection, indicating the VDAC2 and RACK1 were sufficient to allow LCDV entry and infection. All these results collectively showed that VDAC2 and RACK1 function as receptors for LCDV entry and infection.IMPORTANCE Lymphocystis disease virus (LCDV) is the causative agent of lymphocystis disease in fish, which has caused huge economic losses to the aquaculture industry worldwide, but the molecular mechanism underlying the LCDV-host interaction remains unclear. Here, the 27.8-kDa putative cellular receptor for LCDV was identified as voltage-dependent anion channel protein 2 (VDAC2) and receptor of activated protein C kinase 1 (RACK1), and our results revealed that VDAC2 and RACK1 expression was sufficient to allow LCDV entry and that they are functional receptors that initiate LCDV infection for the first time, which leads to a better understanding of the molecular mechanism underlying LCDV infection and virus-host interactions.

Development and Characterization of Monoclonal Antibodies to the 32 kDa Viral Attachment Protein of Lymphocystis Disease Virus and Their Neutralizing Ability in Vitro

In previous research, a 32 kDa protein in lymphocystis disease virus (LCDV) was identified as viral attachment protein (VAP) that specifically interacted with the 27.8 kDa cellular receptor from flounder Paralichthys olivaceus gill (FG) cells, and the recombinant VAP (rVAP) was expressed in Escherichia coli strain BL21 (DE3). In this study, monoclonal antibodies (MAbs) against 32 kDa VAP are produced by immunization of BALB/c mice with the rVAP. Seven hybridoma secreting MAbs were screened by enzyme-linked immunosorbent assay, five of which designated as 1C6, 1C8, 3B5, 3D11 and 3H10 are cloned by the limiting dilution method, depending on the strongly positive results of ELISA. Western blotting analysis shows that the five MAbs can specifically react with the 32 kDa protein of LCDV and the purified 50 kDa rVAP, and the subtype of the MAbs is identified as IgG. Immunofluorescence results demonstrate that the specific fluorescence signals for LCDV appear in the cytoplasm of FG cells at 24 h post LCDV infection. Neutralization assay results indicate that pre-incubations of LCDV with the five MAbs can significantly decrease the LCDV copy numbers and delay the development of the cytopathic effect in FG cells, revealing that the five MAbs can neutralize the LCDV particles and block viral infection in vitro. The neutralizing MAbs against 32 kDa VAP would be useful for the study on the LCDV⁻host interaction and might be promising inhibitors of LCDV infection in fish.

Transcriptome Analysis of Flounder (Paralichthys olivaceus) Gill in Response to Lymphocystis Disease Virus (LCDV) Infection: Novel Insights into Fish Defense Mechanisms

Lymphocystis disease virus (LCDV) infection may induce a variety of host gene expression changes associated with disease development; however, our understanding of the molecular mechanisms underlying host-virus interactions is limited. In this study, RNA sequencing (RNA-seq) was employed to investigate differentially expressed genes (DEGs) in the gill of the flounder (Paralichthys olivaceus) at one week post LCDV infection. Transcriptome sequencing of the gill with and without LCDV infection was performed using the Illumina HiSeq 2500 platform. In total, RNA-seq analysis generated 193,225,170 clean reads aligned with 106,293 unigenes. Among them, 1812 genes were up-regulated and 1626 genes were down-regulated after LCDV infection. The DEGs related to cellular process and metabolism occupied the dominant position involved in the LCDV infection. A further function analysis demonstrated that the genes related to inflammation, the ubiquitin-proteasome pathway, cell proliferation, apoptosis, tumor formation, and anti-viral defense showed a differential expression. Several DEGs including β actin, toll-like receptors, cytokine-related genes, antiviral related genes, and apoptosis related genes were involved in LCDV entry and immune response. In addition, RNA-seq data was validated by quantitative real-time PCR. For the first time, the comprehensive gene expression study provided valuable insights into the host-pathogen interaction between flounder and LCDV.

Characterization of the gilthead seabream (Sparus aurata L.) immune response under a natural lymphocystis disease virus outbreak

Lymphocystis or lymphocystis disease virus (LCDV) is distributed worldwide and affects many fresh and marine water fish species. LCDV is commonly found in aquaria fish species but also in farmed fish species, among them the gilthead seabream (Sparus aurata L.). The immune status of gilthead seabream (S. aurata) specimens under a natural outbreak of LCDV was studied. The replication of the virus was demonstrated in infected fish, but not in control fish. The results showed decreased total serum IgM levels and increased innate cellular immune response (peroxidase and respiratory burst activities) of head kidney leucocytes in LCDV-infected fish, compared to the values obtained in uninfected specimens. In addition, transcription of antiviral genes (ifn and irf3) was down-regulated in the skin of LCDV-positive fish as well as genes involved in cellular immunity (csf1r, mhc2a, tcra and ighm) that were down-regulated in skin and head kidney of infected fish. By contrast, the transcription of nccrp1 was up-regulated in head kidney after LCDV infection. These present results show that head kidney leucocytes are activated to encounter the virus at the sites of replication.

Tissue Localization of Lymphocystis Disease Virus (LCDV) Receptor-27.8 kDa and Its Expression Kinetics Induced by the Viral Infection in Turbot (Scophthalmus maximus)

The 27.8 kDa membrane protein expressed in flounder (Paralichthys olivaceus) gill cells was proved to be a receptor mediating lymphocystis disease virus (LCDV) infection. In this study, SDS-PAGE and Western blotting demonstrated that 27.8 kDa receptor (27.8R) was shared by flounder and turbot (Scophthalmus maximus). Indirect immunofluorescence assay (IIFA) and immunohistochemistry showed that 27.8R was widely expressed in tested tissues of healthy turbot. The indirect enzyme-linked immunosorbent assay indicated that 27.8R expression was relatively higher in stomach, gill, heart, and intestine, followed by skin, head kidney, spleen, blood cells, kidney and liver, and lower in ovary and brain in healthy turbot, and it was significantly up-regulated after LCDV infection. Meanwhile, real-time quantitative PCR demonstrated that LCDV was detected in heart, peripheral blood cells, and head kidney at 3 h post infection (p.i.), and then in other tested tissues at 12 h p.i. LCDV copies increased in a time-dependent manner, and were generally higher in the tissues with higher 27.8R expression. Additionally, IIFA showed that 27.8R and LCDV were detected at 3 h p.i. in some leukocytes. These results suggested that 27.8R also served as a receptor in turbot, and LCDV can infect some leukocytes which might result in LCDV spreading to different tissues in turbot.

Identification and characterization of a novel lymphocystis disease virus isolate from cultured grouper in China

Grouper Epinephelus spp. is one of the most important mariculture fish species in China and South-East Asian countries. The emerging viral diseases, evoked by iridovirus which belongs to genus Megalocytivirus and Ranavirus, have been well characterized in recent years. To date, few data on lymphocystis disease in grouper which caused by lymphocystis disease virus (LCDV) were described. Here, a novel LCDV isolate was identified and characterized. Based on the sequence of LCDV major capsid protein (MCP) and DNA polymerase gene, we found that the causative agents from different species of diseased groupers were the same one and herein were uniformly defined as grouper LCDV (GLCDV). Furthermore, H&E staining revealed that the nodules on the skin were composed of giant cells that contained inclusion bodies in the cytoplasm. Numerous virus particles with >210 nm in diameter and with hexagonal profiles were observed in the cytoplasm. In addition, phylogenetic analysis based on four iridovirus core genes, MCP, DNA polymerase, myristoylated membrane protein (MMP) and ribonucleotide reductase (RNR), consistently showed that GLCDV was mostly related to LCDV-C, followed by LCDV-1. Taken together, our data firstly provided the molecular evidence that GLCDV was a novel emerging iridovirus pathogen in grouper culture.