Figures
Abstract
White spot syndrome virus (WSSV) is the most devastating virosis threatening the shrimp culture industry worldwide. Variations of environmental factors in shrimp culture ponds usually lead to the outbreak of white spot syndrome (WSS). In order to know the molecular mechanisms of WSS outbreak induced by temperature variation and the biological changes of the host at the initial stage of WSSV acute infection, RNA-Seq technology was used to analyze the differentially expressed genes (DEGs) in shrimp with a certain amount of WSSV cultured at 18°C and shrimp whose culture temperature were raised to 25°C. To analyze whether the expression changes of the DEGs were due to temperature rising or WSSV proliferation, the expression of selected DEGs was analyzed by real-time PCR with another shrimp group, namely Group T, as control. Group T didn’t suffer WSSV infection but was subjected to temperature rising in parallel. At the initial stage of WSSV acute infection, DEGs related to energy production were up-regulated, whereas most DEGs related to cell cycle and positive regulation of cell death and were down-regulated. Triose phosphate isomerase, enolase and alcohol dehydrogenase involved in glycosis were up-regulated, while pyruvate dehydrogenase, citrate synthase and isocitrate dehydrogenase with NAD as the coenzyme involved in TCA pathway were down-regulated. Also genes involved in host DNA replication, including DNA primase, DNA topoisomerase and DNA polymerase showed down-regulated expression. Several interesting genes including crustin genes, acting binding or inhibiting protein genes, a disintegrin and metalloproteinase domain-containing protein 9 (ADAM9) gene and a GRP 78 gene were also analyzed. Understanding the interactions between hosts and WSSV at the initial stage of acute infection will not only help to get a deep insight into the pathogenesis of WSSV but also provide clues for therapies.
Citation: Sun Y, Li F, Sun Z, Zhang X, Li S, Zhang C, et al. (2014) Transcriptome Analysis of the Initial Stage of Acute WSSV Infection Caused by Temperature Change. PLoS ONE 9(3): e90732. https://doi.org/10.1371/journal.pone.0090732
Editor: Irene Söderhäll, Uppsala University, Sweden
Received: November 16, 2013; Accepted: February 3, 2014; Published: March 4, 2014
Copyright: © 2014 Sun et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: This work was financially supported by the Major State Basic Research Development Program (2012CB114403), National High Tech Research and Development Program (2012AA10A404, 2012AA092205), National Natural Science Foundation Program (31072203). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
White spot syndrome (WSS) is the most devastating virosis threatening the shrimp culture industry worldwide [1]. The causative agent, white spot syndrome virus (WSSV), is an enveloped virus with a circular double strand DNA of about 300 kb, belonging to a new virus family Nimaviridae [2], [3]. The virus mainly infects animals of Decapoda and causes high lethality to shrimp [4]. Its outbreak depends on the interactions among the pathogen, host and environment. The interaction between WSSV and the hosts has been a research focus in recent years. Many host membrane proteins such as Rab7, β-integrin and chitin-binding protein showed interactions with WSSV envelope proteins in vitro and were proposed as possible WSSV receptors [5]–[7]. Shrimp STAT, Relish and Dorsal were annexed by WSSV to enhance viral immediate-early genes [8]–[11]. However, these studies are far from enough to illustrate the mechanisms of WSSV pathogenesis.
Environment has a profound impact on the outbreak of WSS. Water temperature is considered to be one of the most important environmental factors for shrimp. Field survey showed that in seasons with temperature lower than 20°C or higher than 30°C, the outbreak of WSS is abated; and temperatures between 22 to 30°C allow WSSV to replicate at a much higher rate [12], [13]. Studies in laboratory also verified the above phenomenon that temperature is crucial in determining WSSV proliferation [14]–[17]. Apoptosis might contribute to the increased survival of infected shrimp maintained at a hyperthermia state [18]. The main effect of hyperthermia on subcuticular epithelial cells was to reduce the expression of WSSV genes rather than to directly induce host genes that might contribute to control the infection [19]. A recent report showed that NAD-dependent aldehyde dehydrogenase (ALDH) and Hsp70 both play an important role in the inhibition of WSSV replication at high temperature [20]. Till present, there is no report about the initiation mechanisms of WSSV acute infection induced by temperature variation.
Exopalamon carincauda Holthuis has the potential to be an ideal experimental animal for crustacean with the virtues stated before [21]. In this study, shrimp E. carincauda, cultured at 18°C, were challenged with WSSV, followed by temperature rising to 25°C to mimic the state transition of the infected shrimp in culture ponds as demonstrated previously [12], [13]. It aims to get a better understanding of the biological changes of shrimp at the initial stage of the rapid propagation of WSSV and explore the initiation mechanisms of the acute infection caused by temperature rising.
Materials and Methods
1 Maintenance of Experimental Shrimp and Viral Inocula Preparation
Experimental shrimp (E. carincauda) were hatched in our lab. They were detected to be WSSV free through one-step PCR as described previously [22]. These shrimp were maintained in natural sea water with continuous aeration and were fed with commercial food pellet twice a day.
WSSV was prepared and quantified as previously described [21]. The WSSV solution, serially diluted to 700 copies/µl with PBS, was used as inocula.
2 Challenge Experiment and Sampling
Three groups, respectively named as Group T, Group W and Group WT, were set with 40 shrimp in each group. These shrimp, with an average length of 3.31±0.33 cm and an average weight of 0.51±0.17 g, were cultured at 18°C. Cephalothorax samples from 7 shrimp of the three groups were respectively collected as samples at 0 hour. For Group W, each shrimp was injected intramuscularly with WSSV inocula at a dose of 7000 copies, and kept at 18°C continuously. For Group WT, each shrimp was also injected with WSSV inocula at a dose of 7000 copies, and the culture temperature was raised to 25°C within 3 hours at 24 h post injection. They were kept at this temperature thereafter. For Group T, shrimp suffered no WSSV challenge. However the culture temperature of Group T was raised to 25°C within 3 hours at the same time with that of Group WT and then they were cultured at 25°C thereafter. At 6, 12 and 24 hours post temperature rising (hptr), 7 shrimp from each group were sacrificed and their cephalothoraxes were collected. Besides cephalothorax, pleopods were also collected from shrimp of Group W and Group WT. All the samples were immediately preserved in liquid nitrogen. The experimental scheme described above was shown in Fig. 1.
Three groups (Group T, Group W and Group WT) of shrimp were cultured at 18°C. Cephalothorax samples from 7 shrimp of the three groups were respectively collected as samples at 0 hour. For Group W, shrimp were challenged with WSSV inocula and kept at 18°C continuously. For Group WT, shrimp were also challenged with WSSV and the culture temperature was raised to 25°C at 24 h post challenge (hpc). For Group T, shrimp suffered no WSSV challenge. However the culture temperature of Group T was raised to 25°C in parallel with Group WT and then the two groups were maintained at 25°C thereafter. At 6, 12 and 24 hours post temperature rising (hptr), cephalothoraxes samples of 7 shrimp from each group were collected. Besides cephalothorax, pleopods were also collected from Group W and Group WT. Cephalothorax samples from Group W and Group WT at 12 hptr were subjected to Illumina sequencing.
3 DNA Extraction and Viral Load Quantification
DNA was extracted from pleopod samples using plant genomic DNA kit (Tiangen, China) following the instructions. Protease K was used additionally at a final concentration of 2 µg/µl for sample digestion. Extracted DNA was subjected to electrophoresis in 1% agarose gel and quantified by NanoDrop 1000 spectrophotometer (Thermo Fisher Scientific Inc., USA).
Four-fold serially diluted solutions of a plasmid harboring the gene encoding the extra-cellular part of WSSV envelope protein VP28 were used as standard to generate a standard curve in the quantitative real-time PCR, performed with primers QVP28F and QVP28R (Table 1) to determine the viral loads in extracted DNA using a mastercycler ep realplex (Eppendorf, Germany), just as previously described [22].
4 RNA Extraction and cDNA Synthesis
Total RNA was extracted from cephalothorax samples with RNAisol reagent (Takara, Japan), assessed by electrophoresis in 1% agarose gel, quantified by NanoDrop 1000 spectrophotometer and treated with RNase free DNase I (Promega, USA) to remove the contaminating DNA. The cDNA synthesis were proceeded at 37°C for 2 hours with 4 µg RNA, 40 µM random primers (Sangon, China), 50 mM dNTP, 160 U RNasin (Promega, USA), 20 µl 5×M-MLV buffer and 800 µM-MLV (Promega, USA) in a total volume of 100 µl.
5 Illumina Sequencing
RNA samples of Group W and Group WT at 12 hptr were subjected to Illumina sequencing using the paired-end RNA-seq method (Bentley et al., 2008). Briefly, poly (A) mRNA was isolated from total RNA with Sera-mag oligo(dT) beads and fragmented to a length of 200 bp to 700 bp. cDNA was synthesized with these fragments and random hexamer primers. After end reparation, addition of poly (A) and ligation with sequencing adapters, fragments with suitable length were recovered from agarose gel and amplified by PCR for 15 cycles. The products were purified with MinElute column (Qiagen), followed by paired-end sequencing using Illumina HiSeq™ 2000.
6 Reads Assembly and Sequence Annotation
Transcriptome de novo assembly is carried out by program Trinity and TIGR Gene Indices clustering tools (TGICL) [23]. The obtained unigenes were annotated by BLASTx with the sequences available in NR database and Swiss-Prot databases (E-value<10−5).
7 Analyzing the Biological Changes of the Host from the Transcriptomic Data
The expression of unigenes was calculated as RPKM (Reads per Kilobase of exon model per Million mapped reads), followed by multiple hypothesis testing and FDR (False Discovery Rate) control. DEGs were screened with the two following standards. Firstly, the FDR should be no larger than 0.001. Secondly, the expression difference between two groups should be larger than 3 times or the expression difference had been testified by real-time PCR to be larger than 2 times.
GO (Gene Ontology, http://www.geneontology.org/) functional classification annotation for DEGs was carried out and the GO ontology “biological process” was analyzed. The up and down-regulated DEGs for GO analysis were listed in Files S1 and S2 separately.
Based on the annotation information, genes involved in the glycometabolism pathways including glycosis and tricarboxylic acid cycle pathway (TCA) as well as DNA replication were manually screened and compared between Group W and Group WT at 12 hptr.
8 Quantitative Real-time PCR Analysis of Selected Genes from the Transcriptome
Some interesting genes were manually screened. Quantitative real-time PCR was performed with primers in Table 1 to analyze the expression levels of the interesting genes in Group T, Group W and Group WT at 6, 12 and 24 hptr as well as to verify the transcriptomic data. The 18S rRNA gene, which served as a stably expressed reference, was quantified with primers 18SF and 18SR (Table 1).
Results and Discussions
1 Mortality and the Viral Load Variations of Group W and Group WT
No mortality happened in Group T during the whole experiment course. In Group WT a few died shrimp began to emerge at about 24 hptr and then WSS broke out. No more death occurred in Group WT after 48 hptr with a final mortality of 94%. In Group W, no mortality appeared until 96 hptr, followed by the WSS outbreak beginning at 120 hptr. The situation in Group W stabilized at about 168 hptr and the final mortality was 96%. The data supported the previous reports that the temperature affected the WSSV infection and lower temperature could reduce WSSV pathogenicity [13], [15], [24], [25]. Lower temperature (18°C) could only delay the infection progression, rather than completely stopping the infection course.
WSSV loads in pleopod samples of Group W and Group WT were quantified by real-time PCR (Fig. 2). The WSSV amount in Group WT was 628.28±187.88 copies/ng DNA, 3.31×104±2.45×104 copies/ng DNA and 1.10×106±4.28×105 copies/ng DNA at 6, 12 and 24 hptr respectively, and t test analysis showed significant difference in the viral amount between samples at 24 hptr and that at 12 hptr (P<0.05). However it was only 201.86±182.08 copies/ng DNA, 395.06±226.99 copies/ng DNA and 364±434.82 copies/ng DNA in Group W at the corresponding time respectively and no significant difference was observed among the viral loads at the above time points by One-way ANOVA (P>0.05).
DNA was extracted from pleopods of Group W and Group WT sampled at 6, 12 and 24 hours post temperature rising. And viral loads in extracted DNA were determined by quantitative real-time PCR with serially diluted solutions of a plasmid harboring the gene encoding the extra-cellular part of WSSV envelope protein VP28 as the standard sample to generate a standard curve. Significant differences between the two groups were indicated by asterisks (P<0.05).
Previous study showed that the viral load in the pleopod at the light infection stage was at the level of no more than 103 copies per ng DNA [22]. When the viral load in the pleopod ranged from 103 to 105 copies per ng DNA, the infection progressed in the logarithmic phase of WSSV replication, during which the WSSV load experienced a sharp leap from 4.16×103 copies/ng DNA at 27 hours post infection (hpi) to 9.65×104 copies/ng DNA at 30 hpi, suggesting a turning point from chronic infection to acute infection [22]. In the present study, Group W was indicated to be at the stage of latent infection when the viral loads in pleopods were not more than 103 copies/ng DNA. The viral load in Group WT at 6 hptr showed no significant difference with that in Group W (P>0.05). However, the WSSV copy number in Group WT at 12 hptr was 3.31×104±2.45×104 copies/ng DNA, implying that Group WT was at the initial stage of the logarithmic phase of WSSV replication [22].
2 Reads Assembly and Sequence Annotation
Totally 54,156,824 clean reads in Group W and 54,258,768 clean reads in Group WT were obtained. After assembly, 61,886 unigenes were obtained (69,855 in Group W and 61,757 in Group WT) with an average length of 626 bp. The transcriptomic dataset was deposited in SRA database (No. SRR1105776) and the detailed information for sequence assembly was shown in Table 2.
Unigene annotation was performed by BLASTx to protein databases (E-value<10−5). A total of 25,399 unigenes (41.04% of all unigenes) were matched to known genes in NR database and 21,777 unigenes (35.19% of all unigenes) were matched to known genes in Swiss-Prot database. Other sequences have no homology to proteins in the above two databases. The phenomenon that in shrimp transcriptomic data, a large percentage of sequences with no matches in the protein database was also reported elsewhere and these unknown sequences, reflecting the limited knowledge that we have about the shrimp genomes [26], [27].
3 DEGs between Group WT and Group W
About 6,417 DEGs were obtained (File S3). 4687 DEGs were down-regulated in Group WT and the other 1,730 DEGs were up-regulated, of which 239 were WSSV genes (File S4), indicating the active transcription of WSSV genes at the initial stage of acute infection.
GO functional classification annotation for DEGs was carried out and 821 DEGs (643 down-regulated ones and 178 up-regulated ones in Group WT) had GO functional classification annotation. The DEGs were involved in various biological processes of the host, such as immune system process, metabolic process, energy production, organelle and cytoskeleton organization, cell cycle, signal transduction, death, cell adhesion, response to stimulus and so on, suggesting that the viral propagation after temperature rising had a comprehensive influence on the host.
3.1 DEGs involved in the energy production and glycometabolic process.
Go analysis showed that almost all the DEGs related to energy production were up-regulated in Group WT, suggesting an increased energy need and a boosted energy production (Table 3). This result was consistent with previous studies that proteins in electron transport and ATP synthesis such as cytochrome oxidase polypeptide, the alpha and beta subunits of mitochondrial ATP synthase, were up-regulated upon WSSV infection [28], [29]. The up-regulation of energy production could meet the energy requirements both of the host in the immune response to fight against viral infection and of the virus during the earlier phases of its rapid replication [29].
Genes encoding the enzymes involved in the glycometabolism were manually screened to analyze the pathway transition of glycometabolism (File S5). Triose phosphate isomerase, enolase and alcohol dehydrogenase which were involved in glycosis, were all up-regulated in Group WT, suggesting that the glycosis was enhanced (Table 4). Pyruvate dehydrogenase, citrate synthase and isocitrate dehydrogenase with NAD as the coenzyme were down-regulated in Group WT, suggesting that a reduced glycometabolism level through the tricarboxylic acid cycle pathway (TCA) (Table 4). However the isocitrate dehydrogenase with NADP as the coenzyme was up-regulated in Group WT, which may be the compensatory expression for the inhibited expression of isocitrate dehydrogenase with NAD as the coenzyme. The expression profiles of ATP-citrate synthase (CL4068.Contig2_All) and triosephosphate isomerase B (Unigene3409_All) in the three groups were analyzed by real-time PCR to check whether the changes in glycometabolism were caused by temperature variation or WSSV proliferation (Fig. 3). The expression of ATP-citrate synthase showed no significant difference between Group T and Group W at 12 hptr (P>0.05), but its expression in Group WT at the same time was significantly lower than those in Group T and Group W. As for triosephosphate isomerase B, its expression in Group WT was remarkably higher than those in Group W and Group T (P<0.05) at both 12 hptr and 24 hptr. These data suggested that the TCA pathway was possibly weakened and the glycosis was enhanced due to the rapid replication of WSSV. It was reported that the normal flux in TCA cycle was disrupted by the infectious hypodermal and hematopoietic necrosis virus (IHHNV) in shrimp at least temporally and it was reckoned that the interruption of the Krebs cycle may favor the inducement of a higher lipogenic and cholesterogenic activity, which may ultimately be essential for the virus replication [30]. Metabolic changes resembling the Warburg effect in shrimp hemocytes were detected in the early stage of WSSV infection and it was proposed that the enhanced aerobic glycosis and ATP energy production favored the virus during the earlier phases of its replication cycle [31].
Their expression in Group T, Group W and Group WT at 6, 12 and 24 hours post temperature rising was tested by real-time PCR. The values are shown as mean±SE. Significant differences among the three groups are indicated by different characters (P<0.05).
3.2 DEGs involved in the cell death process.
GO analysis showed that DEGs related to positive regulation of cell death were down-regulated and DEGs related to anti-apoptosis were up-regulated in Group WT, implying that the virus inhibited the host cell death to benefit its own propagation (Table 5). Apoptosis has been regarded as a nonspecific defense mechanism against virus infection in eukaryotic cells, aborting virus multiplication by a premature lysis of infected cells [32]–[34]. Viruses have evolved various strategies to escape apoptosis of the host through viral anti-apoptosis genes to guarantee the production of enough viral progeny [35], [36]. Apoptotic-related proteins showed a high expression level in the virus resistant shrimp, suggesting that apoptosis might play important roles in shrimp immunity against WSSV infection [37], [38]. The enhancement of apoptotic activity effectively inhibited the WSSV infection in shrimp, resulting in decreased shrimp mortality, while apoptosis inhibitor increased the mortality of shrimp infected by WSSV [39]. WSSV showed the ability to modulate apoptosis of shrimp cell. Two WSSV anti-apoptosis proteins AAP-1 and WSSV222 were identified [40]–[43]. AAP-1 is a direct caspase inhibitor and WSSV222 functions through ubiquitin-mediated degradation and blocks apoptosis induced by the tumor suppressor-like protein [36]. In the present study, AAP-1 (Unigene30118_All) and WSSV222 (Unigene30290_All) were all remarkably up-regulated in Group WT at 12 hptr suggesting that they might be responsible for the inhibited apoptosis at the initial stage of WSSV acute replication. Most DEGs with annotation to E3 ubiquitin-protein ligase (9/10, File S6) were down-regulated. Among them, two E3 ubiquitin-protein ligase genes (CL1901.Contig1_All and Unigene37574_All) were further confirmed by real-time PCR (Fig. 4). Ubiquitin (Ub) E3 ligases were regarded as key regulators of cell death, which controlled many pro- and anti-apoptotic molecules at the protein level by Ub-dependent degradation or non-degradative ubiquitylation events [44]. The present data showed that the malfunction to control cell death happened at the initial stage of acute infection.
Their expression in Group T, Group W and Group WT at 6, 12 and 24 hours post temperature rising was tested by real-time PCR. The values are shown as mean±SE. Significant differences among the three groups are indicated by different characters (P<0.05).
3.3 DEGs involved in the cell cycle process.
Go analysis displayed that almost all the DEGs related to cell cycle were down-regulated in Group WT, suggesting that the host cell proliferation was apparently inhibited (Table 6). The cell cycle consists of four phases, gap 1 (G1), synthesis (S), gap 2 (G2) and mitosis (M) in succession. In order to carry through this vital process orderly, a series of mechanisms have been developed and cyclins are a family of controller proteins. In the present study, 10 DEGs were annotated as cyclins, which were all down-regulated (File S7). The expression profiles of two cyclin genes (Unigene42832_All and CL2571.Contig1_All) in the three groups were further confirmed by real-time PCR (Fig. 5). The expression of both genes in Group WT at 12 hpi was nearly completely repressed, compared with that in Group T and Group W. These data further indicated that the inhibition of host cell proliferation was caused by WSSV; and the arrested state started at the initial stage of WSSV acute infection and lasted until 24 hptr at least when shrimp death began to occur. Some viruses were reported to be able to arrest host cells at certain cell phases with different mechanisms [45]–[49]. Autographa californica Nucleopolyhedrovirus (AcMNPV) arrested host cell at G2/M phase due to a virus-coded protein with cdc2-associated kinase activity [46]. Vpr protein of human immunodeficiency virus type 1 (HIV-1) activated cell cycle inhibitor p21/Waf1/Cip1, causing cellular arrest in G2/M phase [47]. Human cytomegalovirus (HCMV) induced a state of late G1 arrest and mouse hepatitis virus (MHV) infection resulted in an arrest at the G0/G1 phase [45], [48]. It was reported that alterations in cell cycle transit during virus infection might be very important for the virulence difference between Malignant fibroma virus (MV) and Shope fibroma virus (SFV) [50]. This may be also one of the causes of the high virulence of WSSV. The molecular mechanisms underlying the alteration in cell cycle transit needs to be further clarified to find a target for controlling the WSSV disease.
Their expression in Group T, Group W and Group WT at 6, 12 and 24 hours post temperature rising was tested by real-time PCR. The values are shown as mean±SE. Significant differences among the three groups are indicated by different characters (P<0.05).
3.4 DEGs involved in cellular DNA replication.
Proteins participated in DNA replication were manually screened. DNA primase, DNA topoisomerase and DNA polymerase showed down-regulated expression in Group WT (see in the appended File S8), whereas the WSSV ORFs wsv514 and wsv447, which are homologous to eukaryotic DNA polymerase and helicase respectively, were remarkably up-regulated (Table 7). The expression profiles of Unigene83_All (a DNA primase), Unigene20038_All (a DNA polymerase) and CL2500.Contig1_All (a DNA topoisomerase) in the three groups at 6, 12 and 24 hptr were further detected by real-time PCR (Fig. 6). The expression of these three genes in Group WT at 12 hptr was significantly lower than that in Group T and Group W at 12 hptr (P<0.05), indicating that their expression had been repressed due to the rapid replication of WSSV. And as for Unigene83_All and Unigene20038_All, the repressed state lasted till 24 hptr at least when shrimp death began to occur. All these data indicated that WSSV proliferation inhibited the cellular DNA replication process and WSSV may rely on its own polymerase and helicase for DNA replication. The adverse effect on the cellular DNA synthesis in the infection course of other viruses such as frog virus 3, herpes simplex virus, HCMV and AcMNPV were also reported and the repressed state of cellular DNA replication ensures that the virus will have unrivalled access to precursors of DNA synthesis [45], [51]–[53]. The decline of cellular DNA synthesis during the infection of frog virus 3 and herpes simplex virus might be caused by the infecting particles themselves and does not require the synthesis of new RNA or protein [51], [52]. The mechanisms of the repressed cellular DNA synthesis during WSSV infection still need further study.
Their expression in Group T, Group W and Group WT at 6, 12 and 24 hours post temperature rising was tested by real-time PCR. The values are shown as mean±SE. Significant differences among the three groups are indicated by different characters (P<0.05).
4 Analysis of Some Interested Genes from the Transcriptomic Data
Some interesting genes, including actin binding or inhibiting protein genes, ADAM9 (disintegrin and metalloproteinase domain-containing protein 9), crustin as well as HSP70, HSP90 and GRP78 (glucose-regulated protein 78) were manually screened and the expression levels of part of the above genes were analyzed by quantitative real-time PCR. The expression patterns of these genes in Group WT at 12 hptr compared to Group W were consistent with the transcriptomic data, suggesting that the transcriptomic data is credible and can reflect the actual gene expression pattern.
Actin and integrin were regarded as possible WSSV receptors [6], [54]. In the present study, DEGs annotated as actin binding or inhibiting protein and a DEG annotated as ADAM9 were all down-regulated at the initial stage of WSSV acute infection (File S9). ADMA9 were reported to be able to bind to integrin [55], [56]. Thus, these genes might compete with WSSV for the binding to actin or disintegrin. It was guessed that the down-regulation of these DEGs was possibly resulted from temperature rising and their down-regulation was responsible for the breakout of WSSV acute infection. Therefore, the expression profiles of two actin binding genes (Unigene7524_All, and Unigene42358_All) and the ADMA9 (Unigene37466_All) in the three groups at the sampling time were detected by real-time PCR to analyze whether the down-regulation of these genes was an inducement or a result of WSSV acute proliferation.
The expression of the two actin binding genes (Unigene7524_All and Unigene42358_All) in Group TW were significantly lower than that in Group T and Group W at both 12 hptr and 24 hptr (P<0.05) and no significant expression differences of the two genes were observed between Group T and Group W (P>0.05) (Fig. 7). These results indicated that their expression had been repressed as a result of the rapid replication of WSSV.
Their expression in Group T, Group W and Group WT at 6, 12 and 24 hours post temperature rising was tested by real-time PCR. The values are shown as mean±SE. Significant differences among the three groups are indicated by different characters (P<0.05).
The expression of ADMA9 was completely blocked in GroupWT at 12 hptr and 24 hptr and different with the expression of the above two actin binding genes, its expression in Group T at 12 hptr were also significantly lower than that in Group W (P<0.05) (Fig. 8), indicating that temperature rising resulted in the down-regulation of the expression of ADMA9 at 12 hptr. It was supposed that the down-regulation of ADMA9 after temperature rising resulted in the release of WSSV receptor integrin, bringing about the acute proliferation of WSSV and the expression of ADMA9 was further completely interdicted for the need of more integrin receptors in WSSV acute infection. So ADMA9 might be a possible inducer of WSSV acute infection.
Its expression in Group T, Group W and Group WT at 6, 12 and 24 hours post temperature rising was tested by real-time PCR. The values are shown as mean±SE. Significant differences among the three groups are indicated by different characters (P<0.05).
It was interesting that 6 out of 7 crustin genes in the DEGs were down regulated and one of them (Unigene37705_All) was completely repressed (File S10). When we detected the expression of Unigenes37705 by real-time PCR, we found that its expression was completely inhibited in Group WT at 12 hptr, the initial stage of WSSV acute infection. By contrast, it was significantly up-regulated in Group W at 6 hptr and 12 hptr, when the shrimp of this group were still at the stage of latent infection (Fig. 9). Crustin is a widely distributed family of antimicrobial peptide within decapod of crustaceans and was regarded as an important component of the innate immune system [57]. Most reports about crustins focused on their anti-bacterial function and little is known about their interaction with virus. A significant down-regulation of the crustin-like AMP was noted post WSSV challenge in Penaeus monodon. After administrating immunostimulants, such as the marine yeast Candida haemulonii, the expression of crustin-like AMP was enhanced, and the shrimp showed increased viability against WSSV infection [58]. In Macrobrachium rosenbergii, which has much stronger resistance to WSSV than penaeid shrimp, the expression of crustin remained up-regulation during WSSV infection [59]. Together with this study, it could be indicated that crustin may function as a WSSV-resistant gene and was possibly a responsible factor for the slowed infection pace of Group W. The virulence of viruses is partially related to their ability to suppress the host anti-viral defense [29], and the results demonstrated that WSSV is highly pathogenic.
Its expression in Group T, Group W and Group WT at 6, 12 and 24 hours post temperature rising was tested by real-time PCR. The values are shown as mean±SE. Significant differences among the three groups are indicated by different characters (P<0.05).
It is also very interesting to note that HSP70, HSP90 and glucose-regulated protein 78 (GRP78) genes were all up-regulated in Group WT at 12 hptr (File S11). Heat-shock cognate protein 70 (Hsc70) of shrimp was a binding partner of WSSV envelop protein VP28 during infection and their association was specific, ATP-dependent and Hsc70 concentration dependent [60]. Hsp70 could inhibit WSSV replication at high temperature [20]. Thus the roles of these heat shock proteins during WSSV infection deserve much attention. The expression profile of one GRP78 gene (Unigene9162_All) was detected in the three groups at the sampling time (Fig. 10). At 6 hptr, the expression of GRP78 in GroupWT was significantly higher than that in Group T (P<0.05) and its high expression level in Group WT remained thereafter. This result was consistent with previous report that GRP78 kept a high expression level in hepatopancreas all through the WSSV infection course [61]. At 24 hptr, its expression in Group WT and Group W were all significantly higher than that in Group T (P<0.05). These results implied that the up-regulated expression of this GRP78 gene was due to WSSV infection. GRP78 functions not only in unfolded protein response and oxidative stress, but also in anti-apoptosis process [62], [63]. The up-regulation of GRP78 was induced by the infection of many other viruses including hepatitis B virus (HBV), hepatitis C virus, HCMV, rotavirus and Dengue virus (DENV) [64]–[67]. However it exhibits distinct functions in different virus infections. GRP78 acts as an endogenous antiviral factor against HBV [65], whereas it could be annexed by other viruses to facilitate the infection. GRP78 may be utilized as a chaperone for viral-protein production by DENV, for virion release by HMCV and for infectious virion maturation by rotavirus [66]. The GRP78 was activated during WSSV infection but further study is needed to clarify whether it was involved in WSSV pathogenesis or in shrimp immunity against the virus.
Its expression in Group T, Group W and Group WT at 6, 12 and 24 hours post temperature rising was tested by real-time PCR. The values are shown as mean±SE. Significant differences among the three groups are indicated by different characters (P<0.05).
In summary, through the expression analyses of all the above genes as well as some other ones in Group W and Group WT at 12 hptr by real-time PCR, which were consistent with the transcriptomic data, it was testified that the transcriptomic data are creditable. At the initial stage of WSSV acute infection after temperature rising from 18°C to 25°C, WSSV exploited the cellular machinery for the seek of its own proliferation and the shrimp went through a lot of biological changes as this study demonstrated. For more information about the inducing mechanisms of WSSV acute infection after temperature change, sampling had better be carried out earlier when the infection of Group WT was at the end of the latent stage and the virus hasn’t undertaken rapid replication. Understanding the interactions between the host and WSSV at the initial stage of acute infection will not only help to get a deep insight into the pathogenesis of WSSV but also provide clues for therapies.
Supporting Information
File S1.
File submitted to the website for GO analysis of up-regulated DEGs.
https://doi.org/10.1371/journal.pone.0090732.s001
(OUT)
File S2.
File submitted to the website for GO analysis of down-regulated DEGs.
https://doi.org/10.1371/journal.pone.0090732.s002
(OUT)
File S5.
DEGs used for analysis of TCA and glycosis.
https://doi.org/10.1371/journal.pone.0090732.s005
(XLS)
File S6.
DEGs annotated as E3 ubiquin-protein.
https://doi.org/10.1371/journal.pone.0090732.s006
(XLS)
File S8.
DEGs used for analysis of DNA synthesis process.
https://doi.org/10.1371/journal.pone.0090732.s008
(XLS)
File S9.
DEGs annotated as actin binding or suppressor gene and a DEG annotated as ADAM gene.
https://doi.org/10.1371/journal.pone.0090732.s009
(XLS)
File S11.
DEGs annotated as GRP78, HSP70 or HSP90.
https://doi.org/10.1371/journal.pone.0090732.s011
(XLS)
Author Contributions
Conceived and designed the experiments: FL JX. Performed the experiments: YS SL. Analyzed the data: YS ZS. Contributed reagents/materials/analysis tools: XZ SL CZ. Wrote the paper: YS FL.
References
- 1. Chou H, Huang C, Wang C, Chiang H, Lo C (1995) Pathogenicity of a baculovirus infection causing white spot syndrome in cultured penaeid shrimp in Taiwan. Dis Aquat Organ 23: 165–173.
- 2. van Hulten MCW, Witteveldt J, Peters S, Kloosterboer N, Tarchini R, et al. (2001) The white spot syndrome virus DNA genome sequence. Virology 286: 7–22.
- 3. Yang F, He J, Lin X, Li Q, Pan D, et al. (2001) Complete genome sequence of the shrimp white spot bacilliform virus. J Virol 75: 11811–11820.
- 4.
Oidtmann B, Stentiford G (2011) White spot syndrome virus (WSSV) concentrations in crustacean tissues–a review of data relevant to assess the risk associated with commodity trade. Transbound Emerg Dis: 14.
- 5. Sritunyalucksana K, Wannapapho W, Lo CF, Flegel TW (2006) PmRab7 is a VP28-binding protein involved in white spot syndrome virus infection in shrimp. J Virol 80: 10734–10742.
- 6. Li DF, Zhang MC, Yang HJ, Zhu YB, Xu X (2007) [beta]-integrin mediates WSSV infection. Virology 368: 122–132.
- 7. Chen KY, Hsu TC, Huang PY, Kang ST, Lo CF, et al. (2009) Penaeus monodon chitin-binding protein (PmCBP) is involved in white spot syndrome virus (WSSV) infection. Fish Shellfish Immunol 27: 460–465.
- 8. Liu WJ, Chang YS, Wang AH, Kou GH, Lo CF (2007) White spot syndrome virus annexes a shrimp STAT to enhance expression of the immediate-early gene ie1. J Virol 81: 1461–1471.
- 9.
Huang XD, Zhao L, Zhang HQ, Xu XP, Jia XT, et al.. (2010) Shrimp NF-[kappa] B binds to the immediate-early gene ie1 promoter of white spot syndrome virus and upregulates its activity. Virology: 176–180.
- 10. Chen YH, Jia XT, Huang XD, Zhang S, Li M, et al. (2011) Two Litopenaeus vannamei HMGB proteins interact with transcription factors LvSTAT and LvDorsal to activate the promoter of white spot syndrome virus immediate-early gene ie1. Mol Immunol 48: 793–799.
- 11. Wang PH, Gu ZH, Wan DH, Zhang MY, Weng SP, et al. (2011) The shrimp NF-kappaB pathway is activated by white spot syndrome virus (WSSV) 449 to facilitate the expression of WSSV069 (ie1), WSSV303 and WSSV371. PLoS One 6: e24773.
- 12. Withyachumnarnkul B, Boonsaeng V, Chomsoong R, Flegel TW, Muangsin S, et al. (2003) Seasonal variation in white spot syndrome virus-positive samples in broodstock and post-larvae of Penaeus monodon in Thailand. Dis Aquat Organ 53: 167–171.
- 13. Moser JR, Álvarez DAG, Cano FM, Garcia TE, Molina DEC, et al. (2012) Water temperature influences viral load and detection of white spot syndrome virus (WSSV) in Litopenaeus vannamei and wild crustaceans. Aquaculture 326–329: 9–14.
- 14. Du HH, Li WF, Xu ZR, Kil ZS (2006) Effect of hyperthermia on the replication of white spot syndrome virus (WSSV) in Procambarus clarkii. Dis Aquat Organ 71: 175.
- 15. Du H, Dai W, Han X, Li W, Xu Y, et al. (2008) Effect of low water temperature on viral replication of white spot syndrome virus in Procambarus clarkii. Aquaculture 277: 149–151.
- 16. Gao H, Kong J, Li Z, Xiao G, Meng X (2011) Quantitative analysis of temperature, salinity and pH on WSSV proliferation in Chinese shrimp Fenneropenaeus chinensis by real-time PCR. Aquaculture 312: 26–31.
- 17. You X, Su Y, Mao Y, Liu M, Wang J, et al. (2010) Effect of high water temperature on mortality, immune response and viral replication of WSSV-infected Marsupenaeus japonicus juveniles and adults. Aquaculture 305: 133–137.
- 18. Granja CB, Aranguren LF, Vidal OM, Aragón L, Salazar M (2003) Does hyperthermia increase apoptosis in white spot syndrome virus (WSSV)-infected Litopenaeus vannamei? Dis Aquat Organ 54: 73–78.
- 19. Reyes A, Salazar M, Granja C (2007) Temperature modifies gene expression in subcuticular epithelial cells of white spot syndrome virus-infected Litopenaeus vannamei. Dev Comp Immunol 31: 23–29.
- 20. Lin YR, Hung HC, Leu JH, Wang HC, Kou GH, et al. (2011) The role of aldehyde dehydrogenase and Hsp70 in suppression of white spot syndrome virus replication at high temperature. J Virol 85: 3517–3525.
- 21. Sun Y, Li F, Chi Y, Xiang J (2013) Enhanced resistance of marine shrimp Exopalamon carincauda Holthuis to WSSV by injecting live VP28-recombinant bacteria. Acta Oceanologica Sinica 32: 52–58.
- 22. Sun Y, Li F, Xiang J (2012) Analysis on the dynamic changes of the amount of WSSV in Chinese shrimp Fenneropenaeus chinensis during infection. Aquaculture 376–379: 124–132.
- 23. Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, et al. (2011) Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol 29: 644–652.
- 24. Guan Y, Yu Z, Li C (2003) The effects of temperature on white spot syndrome infections in Marsupenaeus japonicus. J Invertebr Pathol 83: 257–260.
- 25. Jiravanichpaisal P, Söderhäll K, Söderhäll I (2004) Effect of water temperature on the immune response and infectivity pattern of white spot syndrome virus (WSSV) in freshwater crayfish. Fish Shellfish Immunol 17: 265–275.
- 26. Clavero-Salas A, Sotelo-Mundo RR, Gollas-Galvan T, Hernandez-Lopez J, Peregrino-Uriarte AB, et al. (2007) Transcriptome analysis of gills from the white shrimp Litopenaeus vannamei infected with white spot syndrome virus. Fish Shellfish Immunol 23: 459–472.
- 27. Robalino J, Carnegie RB, O’Leary N, Ouvry-Patat SA, de la Vega E, et al. (2009) Contributions of functional genomics and proteomics to the study of immune responses in the Pacific white leg shrimp Litopenaeus vannamei. Vet Immunol Immunopathol 128: 110–118.
- 28. Leu JH, Chang CC, Wu JL, Hsu CW, Hirono I, et al. (2007) Comparative analysis of differentially expressed genes in normal and white spot syndrome virus infected Penaeus monodon. BMC Genomics 8: 120.
- 29. Pongsomboon S, Tang S, Boonda S, Aoki T, Hirono I, et al. (2011) A cDNA microarray approach for analyzing transcriptional changes in Penaeus monodon after infection by pathogens. Fish Shellfish Immunol 30: 439–446.
- 30. Galvan-Alvarez D, Mendoza-Cano F, Hernandez-Lopez J, Sanchez-Paz A (2012) Experimental evidence of metabolic disturbance in the white shrimp Penaeus vannamei induced by the infectious hypodermal and hematopoietic necrosis Virus (IHHNV). J Invertebr Pathol 111: 60–67.
- 31. Chen IT, Aoki T, Huang YT, Hirono I, Chen TC, et al. (2011) White spot syndrome virus induces metabolic changes resembling the warburg effect in shrimp hemocytes in the early stage of infection. J Virol 85: 12919–12928.
- 32. Clem R, Miller L (1993) Apoptosis reduces both the in vitro replication and the in vivo infectivity of a baculovirus. J Virol 67: 3730–3738.
- 33. Everett H, McFadden G (1999) Apoptosis: an innate immune response to virus infection. Trends Microbiol 7: 160–165.
- 34. Koyama AH, Fukumori T, Fujita M, Irie H, Adachi A (2000) Physiological significance of apoptosis in animal virus infection. Microbes Infect 2: 1111–1117.
- 35. Hay S, Kannourakis G (2002) A time to kill: viral manipulation of the cell death program. J Gen Virol 83: 1547–1564.
- 36. Leu JH, Lin SJ, Huang JY, Chen TC, Lo CF (2012) A model for apoptotic interaction between white spot syndrome virus and shrimp. Fish Shellfish Immunol 34: 1011–1017.
- 37. He N, Qin Q, Xu X (2005) Differential profile of genes expressed in hemocytes of white spot syndrome virus-resistant shrimp (Penaeus japonicus) by combining suppression subtractive hybridization and differential hybridization. Antiviral Res 66: 39–45.
- 38. Pan D, He N, Yang Z, Liu H, Xu X (2005) Differential gene expression profile in hepatopancreas of WSSV-resistant shrimp (Penaeus japonicus) by suppression subtractive hybridization. Dev Comp Immunol 29: 103–112.
- 39. Zhi B, Tang W, Zhang X (2011) Enhancement of shrimp antiviral immune response through caspase-dependent apoptosis by small molecules. Mar Biotechnol 13: 575–583.
- 40. Wang Z, Hu L, Yi G, Xu H, Qi Y, et al. (2004) ORF390 of white spot syndrome virus genome is identified as a novel anti-apoptosis gene. Biochem Bioph Res Co 325: 899–907.
- 41. He F, Fenner BJ, Godwin AK, Kwang J (2006) White spot syndrome virus open reading frame 222 encodes a viral E3 ligase and mediates degradation of a host tumor suppressor via ubiquitination. J Virol 80: 3884–3892.
- 42. He F, Syed SM, Hameed AS, Kwang J (2009) Viral ubiquitin ligase WSSV222 is required for efficient white spot syndrome virus replication in shrimp. J Gen Virol 90: 1483–1490.
- 43. Leu JH, Chen LL, Lin YR, Kou GH, Lo CF (2010) Molecular mechanism of the interactions between white spot syndrome virus anti-apoptosis protein AAP-1 (WSSV449) and shrimp effector caspase. Dev Comp Immunol 34: 1068–1074.
- 44. Broemer M, Meier P (2009) Ubiquitin-mediated regulation of apoptosis. Trends Cell Biol 19: 130–140.
- 45. Bresnahan WA, Boldogh I, Thompson EA, Albrecht T (1996) Human cytomegalovirus inhibits cellular DNA synthesis and arrests productively infected cells in late G1. Virology 224: 150–160.
- 46. Braunagel SC, Parr R, Belyavskyi M, Summers MD (1998) Autographa californica nucleopolyhedrovirus infection results in Sf9 cell cycle arrest at G 2/M phase. Virology 244: 195–211.
- 47. Chowdhury IH, Wang XF, Landau NR, Robb ML, Polonis VR, et al. (2003) HIV-1 Vpr activates cell cycle inhibitor p21/Waf1/Cip1: a potential mechanism of G2/M cell cycle arrest. Virology 305: 371–377.
- 48. Chen CJ, Makino S (2004) Murine coronavirus replication induces cell cycle arrest in G0/G1 phase. J Virol 78: 5658–5669.
- 49. De Bolle L, Hatse S, Verbeken E, De Clercq E, Naesens L (2004) Human herpesvirus 6 infection arrests cord blood mononuclear cells in G2 phase of the cell cycle. FEBS Lett 560: 25–29.
- 50. Wali A, Strayer DS (1999) Comparative effects of virulent and avirulent poxviruses on cell cycle progression. Exp Mol Pathol 66: 31–38.
- 51. Guir J, Braunwald J, Kirn A (1971) Inhibition of host-specific DNA and RNA synthesis in KB cells following infection with frog virus 3. J Gen Virol 12: 293–301.
- 52. Fenwick M, Walker M (1978) Suppression of the synthesis of cellular macromolecules by herpes simplex virus. J Gen Virol 41: 37–51.
- 53. Braunagel SC, Parr R, Belyavskyi M, Summers MD (1998) Autographa californica Nucleopolyhedrovirus Infection Results in Sf9 Cell Cycle Arrest at G 2/M Phase. Virology 244: 195–211.
- 54. Xie X, Yang F (2005) Interaction of white spot syndrome virus VP26 protein with actin. Virology 336: 93–99.
- 55. Zhou M, Graham R, Russell G, Croucher PI (2001) MDC-9 (ADAM-9/Meltrin γ) functions as an adhesion molecule by binding the αvβ5 integrin. Biochem Biophys Res Commun 280: 574–580.
- 56. Eto K, Huet C, Tarui T, Kupriyanov S, Liu H-Z, et al. (2002) Functional classification of ADAMs based on a conserved motif for binding to integrin α9β1. J Biol Chem 277: 17804–17810.
- 57. Hauton C, Brockton V, Smith VJ (2006) Cloning of a crustin-like, single whey-acidic-domain, antibacterial peptide from the haemocytes of the European lobster, Homarus gammarus, and its response to infection with bacteria. Mol Immunol 43: 1490–1496.
- 58. Antony SP, Singh IS, Sudheer NS, Vrinda S, Priyaja P, et al. (2011) Molecular characterization of a crustin-like antimicrobial peptide in the giant tiger shrimp, Penaeus monodon, and its expression profile in response to various immunostimulants and challenge with WSSV. Immunobiology 216: 184–194.
- 59. Arockiaraj J, Gnanam AJ, Muthukrishnan D, Gudimella R, Milton J, et al. (2013) Crustin, a WAP domain containing antimicrobial peptide from freshwater prawn Macrobrachium rosenbergii: immune characterization. Fish Shellfish Immunol 34: 109–118.
- 60. Xu H, Yan F, Deng X, Wang J, Zou T, et al. (2009) The interaction of white spot syndrome virus envelope protein VP28 with shrimp Hsc70 is specific and ATP-dependent. Fish Shellfish Immunol 26: 414–421.
- 61. Luan W, Li F, Zhang J, Wang B, Xiang J (2009) Cloning and expression of glucose regulated protein 78 (GRP78) in Fenneropenaeus chinensis. Mol Biol Rep 36: 289–298.
- 62. Little E, Ramakrishnan M, Roy B, Gazit G, Lee AS (1994) The glucose-regulated proteins (GRP78 and GRP94): functions, gene regulation, and applications. Crit Rev Eukaryot Gene Expr 4: 1.
- 63. Rao RV, Peel A, Logvinova A, del Rio G, Hermel E, et al. (2002) Coupling endoplasmic reticulum stress to the cell death program: role of the ER chaperone GRP78. FEBS Lett 514: 122.
- 64. Liberman E, Fong YL, Selby MJ, Choo QL, Cousens L, et al. (1999) Activation of the grp78 andgrp94 Promoters by Hepatitis C Virus E2 Envelope Protein. J Virol 73: 3718–3722.
- 65. Ma Y, Yu J, Chan HL, Chen YC, Wang H, et al. (2009) Glucose-regulated protein 78 is an intracellular antiviral factor against hepatitis B virus. Mol Cell Proteomics 8: 2582–2594.
- 66. Wati S, Soo ML, Zilm P, Li P, Paton AW, et al. (2009) Dengue virus infection induces upregulation of GRP78, which acts to chaperone viral antigen production. J Virol 83: 12871–12880.
- 67. Ou DS C, Lee SB, Chu CS, Chang LH, Chung BC, et al. (2011) Transcriptional activation of endoplasmic reticulum chaperone GRP78 by HCMV IE1–72 protein. Cell Res 21: 642–653.