Skip to main content


  • Research article
  • Open Access

Transcriptome characterization and gene expression of Epinephelus spp in endoplasmic reticulum stress-related pathway during betanodavirus infection in vitro

BMC Genomics201213:651

  • Received: 7 February 2012
  • Accepted: 18 October 2012
  • Published:



Grouper (Epinephelus spp) is an economically important fish species worldwide. However, viral pathogens such as nervous necrosis virus (NNV) have been causing severe infections in the fish, resulting in great loss in the grouper aquaculture industry. Yet, the understanding of the molecular mechanisms underlying the pathogenicity of NNV is still inadequate, mainly due to insufficient genomic information of the host.


De novo assembly of grouper transcriptome in the grouper kidney (GK) cells was conducted by using short read sequencing technology of Solexa/Illumina. A sum of 66,582 unigenes with mean length of 603 bp were obtained, and were annotated according to Gene Ontology (GO) and Clusters of Orthologous Groups (COG). In addition, the tag-based digital gene expression (DGE) system was used to investigate the gene expression and pathways associated with NNV infection in GK cells. The analysis revealed endoplasmic reticulum (ER) stress response was prominently affected in NNV-infected GK cells. A further analysis revealed an interaction between the NNV capsid protein and the ER chaperone immunoglobulin heavy-chain binding protein (BiP). Furthermore, exogenous expression of NNV capsid protein was able to induce XBP-1 mRNA splicing in vivo, suggesting a role of the capsid protein in the NNV-induced ER stress.


Our data presents valuable genetic information for Epinephelus spp., which will benefit future study in this non-model but economically important species. The DGE profile of ER stress response in NNV-infected cells provides information of many important components associated with the protein processing in ER. Specifically, we showed that the viral capsid protein might play an important role in the ER stress response.


  • Endoplasmic Reticulum Stress
  • Bovine Viral Diarrhea Virus
  • Digital Gene Expression
  • Endoplasmic Reticulum Stress Response
  • African Swine Fever Virus


Because of the exquisite quality of flesh, grouper is one of the aquaculture species with high economic value. In past decades, the grouper industry has developed rapidly in many Asian countries. In Taiwan alone, the production volume of grouper was about 9,300 mt in 2006, reaching a total value of NT$ 1.73 billion, approximate US$60 million[1]. However, the sustainable development of grouper industry has been threatened by infection caused by nervous necrosis virus (NNV), which could cause mass mortality of up to 100% in the fish at larval and juvenile stages[25].

NNV, also known as viral encephalopathy and retinopathy (VER), belongs to the genus Betanodavirus of the family Nodaviridae[5]. NNV virion is non-enveloped, spherical-shaped particle with a diameter of 25–34 nm[2, 6, 7]. It consists of two single-stranded, positive-sense RNA genomes: RNA1 (3.1kb) encodes the viral RNA-dependent RNA polymerase (RdRp) and RNA2 encodes the viral capsid protein[8]. In addition, a subgenomic RNA3 from RNA1 contains two open reading frames that encode protein B1 (111 amino acids) and B2 (75 amino acids), which have been shown to be an anti-necrosis and a suppressor of host siRNA silencing factor, respectively[9, 10]. Protein B2 might also function as necrotic death factor[11]. The first NNV outbreak was reported in Japanese parrotfish, Oplegnathus fasciatus (Temminck & Schlegel) in 1990[2], and has been since widely reported in various teleost species, including grouper. In grouper, the infected fish exhibit signs of uncoordinated darting, spiral swimming, abnormalities of swimming bladder control, and vacuolation in brain and retina[2, 4, 1216].

NNV infection remains as now a major threat to the grouper industry, despite extensive efforts have been invested to develop therapeutic and prophylactic regiments against the viral infection. Extended understanding of the molecular mechanism of the viral pathogenicity can facilitate the development of antiviral strategy, yet such knowledge is currently inadequate. For examples, the infection route of NNV in the host is not fully illustrated, and the cellular response of the host toward NNV remains largely in puzzle. NNV is known to cause apoptotic death of infected cells. Recent studies have further indicated a connection between the viral-induced endoplasmic reticulum (ER) stress and the apoptotic death[1719]. Many studies have demonstrated that virus-induced ER stress could determine the fate of infected cell, survival or death[20, 21]. Cell which does not overcome the ER stress will succumb to apoptosis ultimately. In the stressed cell, the ER chaperone immunoglobulin heavy-chain binding protein (BiP), also known as glucose-regulated protein 78 (GRP78), functions as a master control to relieve the stress by initiating the unfolding protein response (UPR) pathway. BiP activates the three mediators of UPR: PKR-like ER kinase (PERK), activating transcription factor 6 (ATF6), and inositol-requiring 1 (IRE1). Activation of PERK could phosphrylate eIF-2α, resulting in attenuation of protein synthesis. On the other hand, ATF6 and IRE1 function at the later stage of UPR as compared to PERK. After being translocated into nucleus, the activated ATF6 stimulates the expression of several chaperones that are capable of refolding misfolded proteins. ATF6 also stimulates the expression of the X-box-binding protein (XBP1) mRNA, which would produce a transcription factor once one of its intron is removed by IRE1. The product of the spliced XBP1 mRNA stimulates the expression of proteins that could facilitate the degradation of misfolded proteins[17]. All together, these mediators alleviate ER stress by assisting protein folding, reducing translocation of newly synthesized proteins into ER, and facilitating protein degradation in ER lumen. Currently, the understanding of NNV-induced ER stress response is insufficient and requires further investigation.

One hurdle to the study of the molecular mechanism of NNV pathogenicity in grouper cells is the lack of sufficient genomic information of the host. The recent advance in the next generation sequencing (NGS) technology has exerted profound impact on the biological science research. The technique has been applied in decoding the genomes of several non-model organisms by transcriptome sequencing, providing valuable information in the understanding of gene function, cell responses and evolution. Transcriptome sequencing is cost effective compared to traditional Sanger sequencing method, and is a suitable alternative to whole genome sequencing as it can provide information of transcribed portion of genes at a lower cost[22].

In this study, grouper kidney cell (GK cell) was used to establish a transcriptome library by using the NGS technology. A total of 66,584 unigenes were obtained by SOAP de novo[23] transcriptome assembly software from 51 million reads of raw mRNA sequencing data. By matching the assembled unigenes with various databases, we have characterized these unigenes into different functional categories and pathways. The analysis resulted in the identification of 25 unigenes that are relevant to ER stress-associated genes, including BiP, PERK, ATF6, IRE1, and XBP-1. We further found that BiP could interact with the NNV capsid protein in the infected cells, suggesting a role of these two proteins in the ER stress response after viral infection. These data will benefit future study of NNV pathogenicity and development of therapeutic or prophylactic treatment against NNV.


De novo sequencing and reads assembly

The grouper kidney GK cell line was used in this study to generate a transcriptome database. After eliminating read with adaptors, read with unknown nucleotides larger than 5% and low quality read, 51,198,090 paired-end reads from total nucleotide of 4,607,828,100 were generated by using the SOAP de novo transcriptome assembly program. By overlapping reads into longer fragments, 204,517 contigs with mean length of 248bp were combined. The contigs were assembled into 79,678 scaffolds, which were further assembled into a total number of 66,582 unigenes with mean length of 603 bp (Table1). The sequence length with highest percentage was in the range of 100–1000 bp (84.5%), followed in succession by 1001–2000 bp (10.7%), 2001–3000 bp (3.2%), and longer than 3000 bp (1.6%).
Table 1

Summary of the grouper GK cell transcriptome

Total Reads


Total Nucleotides (nt)


Total Contig Number


Mean Length of Contig (bp)


Total Contig Length (nt)


Total Scaffold Number


Mean Length of Scaffold (bp)


Total Scaffold Length (nt)


Total Unigene Number


Mean Length of Unigene (bp)


Total Length of all Unigene (nt)


Gene ontology (GO) and clusters of orthologous groups (COG) annotation of unigenes

The GK transcriptome unigenes were annotated according to two functional annotation conventions, Gene Ontology (GO) and Clusters of Orthologous (COG). The unigenes were classified into the three categories of GO: molecular function, cellular component and biological process, and 51 sub-categories. Because multiple GO functions were assigned to a single unigene in some cases, 76,975 assignments were generated from the 66,582 unigenes. The biological process was the dominant category, which made up of 49.57%, followed by cellular component (35.45%), and molecular function (14.99%) (Additional file1: Figure S1). A total of 19,261 unigenes were grouped into the 25 COG categories. The largest group in COG is ‘general function prediction only’ (3673, 19.07%), followed by group ‘transcription’ (1708, 8.87%) and group ‘replication, recombination and repair’ (1680, 8.72%). The following groups: ‘defense mechanism’ (73, 0.38%), ‘extracellular structure’ (15, 0.08%) and ‘nuclear structure’ (8, 0.04%) represented the smallest groups (Additional file2: Figure S2).

Expression and functional annotation of unigenes associated with ER stress-related pathway

The tag-based DGE analysis was performed to characterize the gene expression profile associated with NNV infection in GK cells. The comparison of unigene expression was measured by the RPKM (Reads Per kb per Million reads) method, which can be directly used for comparing the difference of gene expression among samples[24]. The RPKM value was compared between NNV-infected and non-infected cells at 6 and 33 hr post-infection. As shown in Figure1, there was slight increase in the unigene expression in the NNV-infected sample at 6 hr post-infection, as indicated by the RPKM values scattering toward the X-axis (NNV-infected sample). At 33 hr post-infection, substantial numbers of unigenes were either up-regulated or down-regulated in the NNV-infected cells. To further characterize the biological functions that were affected by NNV, the 66,582 unigenes from the GK transcriptome database were subjected to analysis against Kyoto Encyclopedia of Genes and Genomes (KEGG) to assign their functions within known biological pathways. A total of 19505 unigenes were mapped into 218 KEGG pathways, among which 555 unigenes (2.85%) were related with the protein processing in ER.
Figure 1
Figure 1

Comparison of unigene expression among different DGE samples. The gene expression level was calculated by using RPKM method (Reads Per kb per Million reads). As the RPKM method can eliminate the influence of different gene length and sequencing discrepancy on the calculation of gene expression, the calculated gene expression can be directly used for comparing the difference of gene expression among samples. The expression of unigenes in response to NNV infection was analyzed by comparing the RPKM of unigenes in the infected versus non-infected GK cells at 6 (A) and 33 (B) hours post-infection. NNV: NNV-infected GK cells, Control: non-infected GK cells.

Analysis of the digital gene expression profile in NNV-infected GK cells revealed that the ER stress pathway was prominently affected by NNV infection (Figure2). Out of the 555 ER-associated unigenes, a total of 117 unigenes were affected in NNV-infected GK cells. Among the genes associated with the 117 unigenes, 16 (related to 25 unigenes) are main players of the UPR and ER-associated degradation (ERAD) processes. These unigenes have been deposited to the NCBI database and are summarized with the assigned accession numbers in Table2. At 33 hr post-infection, several major components of the UPR pathway were up-regulated, including BiP, PERK, ATF6, IRE1, GADD34, CHOP and XBP-1. Table3 lists the fold of up-regulation of these genes. Among the genes, CHOP and GADD34, both are components of the PERK pathway, were also up-regulated at 6 hr post-infection. To verify the expression profile of these major UPR genes in response to NNV infection, we performed a time-course study by using quantitative RT-PCR analysis. As shown in Figure3, GADD34 and CHOP were the two most affected genes at 6 hr post-infection, with 6-fold and 4.3-fold increase, respectively. Except for PERK, the relative expression of these genes at 33 hr post-infection detected by quantitative RT-PCR correlated well with the profile obtained by DGE analysis (Figure3 and Table3).
Figure 2
Figure 2

Identification of unigenes relevant to protein processing in endoplasmic reticulum. Analysis of the digital gene expression profile in NNV-infected GK cells at 33 hours post-infection identified a total of 117 unigenes that are assigned to genes known to be associated with protein processing in ER pathway. Red boxes refer to genes whose associated unigene(s) was up-regulated in the NNV-infected GK cells. Green boxes refer to genes whose associated unigene(s) was down-regulated. Boxes with both red and green colors indicate genes that might have isoforms.

Table 2

Summary of unigenes relevant to ER stress identified in NNV-infected GK cells


Function catagory


Accession Number

Name in NCBI





DNA damage-inducible transcript 3 protein-like







AMP-dependent transcription factor 6





Glucose-regulated protein 78 (GRP78/BiP)







x-box binding protein 1







Endoplasmic reticulum to nucleus signaling 2





Tumor necrosis factor receptor-associated factor 2







Mitogen-activaed protein kinase kinase kinase 5 (MAP3K5/ASK1)





Nuclear factor erythroid 2-related factor 2







Growth arrest and DNA damage-inducible protein/Protein phosphatase 1, regulatory subunit 15A


Protein regulation by luminal chaperones



Glucose-regulated protein 94





Heat shock protein 70





Heat shock protein 90





Endoplasmic oxidoreductase-1-like protein









Protein disulfide-isomerase A6









Translocon-associated protein subunit alpha





Apoptosis regulator BAX-like





Calpain-like protein

Table 3

Differential expression of ER-stress related unigenes in NNV infected GK cells at 6 and 33 hours post-infection


Expression Level

6 Hours Fold Change

33 Hours Fold Change

























Figure 3
Figure 3

Expression of CHOP, PERK, ATF6, BiP, XBP1 and GADD34 in NNV-infected GK cells. GK cells were infected with NNV at MOI of 10, and harvested at 6 and 33 hr post-infection. Triplicate samples were collected. The expression level of CHOP, PERK, ATF6, BiP, XBP1 and GADD34 was measured by quantitative RT-PCR, and the relative expression of each gene was normalized against the expression level of β-actin. The average relative expression of each set of sample is indicated by the numbers above bars. The significance of difference was analyzed by Student’s t-test (*0.01 ≤ p ≤ 0.05; **p ≤ 0.01). Vertical bars represent the standard deviation (n = 3).

Expression of BiP and the interaction between BiP and NNV capsid protein

BiP is a master regulator of UPR response. The up-regulation of BiP in NNV-infected cells inspired us to investigate how BiP might interact with viral components in the infected cells. From our de novo sequencing data, we have identified two fragments of grouper BiP cDNA. The first fragment (unigene44883) represents the first 109 residues of the open reading frame, while the second fragment (unigene52488) is of 155 a.a., with 3 residues overlapped with the C’-terminus of the first fragment (Figure4). A partial grouper BiP open reading frame obtained by compiling the two fragments showed greater than 93% and 87% in sequence similarity and identity with other vertebrate counterparts, respectively (Figure4).
Figure 4
Figure 4

Comparison of BiP amino sequence (partial) among grouper and other vertebrates. Two unigenes (unigene44883 and unigene52488) of grouper BiP were obtained. Unigene44883 encodes the first 109 residues of the Bip open reading frame, and unigene52488 encodes 155 a.a., with 3 residues overlapped with the C’-terminus of the fragment of the first 109 a.a. A partial grouper BiP open reading frame showed greater than 93% and 87% in sequence similarity and identity with other vertebrates, respectively.

The conservation in the amino sequence of BiP protein throughout evolution made it feasible for us to utilize commercially available BiP antibodies to further investigate the expression of BiP protein and its potential interaction with viral components. A series of functional assays were carried out with GF-1 cells, which has been popularly used as a model cell line for NNV infection[25]. As shown in Figure5, the expression of BiP protein co-localized with the NNV capsid protein in the infected GF-1 cells. In non-infected cells, BiP protein formed as aggregated speckles dispersed throughout the cell cytoplasm, while in the infected cells, NNV capsid protein displayed both nuclear localization and a distribution pattern of interconnected web apart from perinuclear region. Furthermore, BiP protein in infected cells also exhibited punctuated pattern and strictly localized to the cytoplasmic area where capsid protein expressed. These results indicate an interaction between BiP and the viral capsid protein.
Figure 5
Figure 5

Sub-cellular localization of NNV capsid protein and BiP protein in NNV-infected GF-1 cells. Grouper fin cell line GF-1 was infected with NNV at MOI of 10. At 24 hr post-infection, the cells were fixed and co-stained with rabbit anti-NNV capsid and mouse anti-BiP, followed by Texas Red conjugated anti-rabbit IgG and Fluorscein conjugated mouse IgG. The Nucleus was visualized by counter-staining with 4', 6-diamidino-2-phenylindole (DAPI). The NNV capsid protein (A) is shown to co-localize with BiP protein (B) in distinct area within the cytoplasm as indicated by the merged image (C). Uninfected GF1 cell was indicated by white arrow.

The interaction between BiP and NNV capsid protein was further confirmed by a co-immunoprecipitation assay, in which the capsid protein was co-precipitated with the BiP protein captured by anti-BiP antibodies (Figure6A). This result strongly supports the notion of an interaction between BiP and NNV capsid protein. In addition, it may suggest a potential role of the NNV capsid protein in the ER stress response. We thus injected the NNV capsid gene into zebrafish embryos to evaluate the effect of capsid protein on ER stress in vivo. The splicing pattern of XBP-1 mRNA was measured as an indicator of the ER stress response. As shown in Fig6B, exogenous expression of NNV capsid protein resulted in XBP-1 mRNA splicing as demonstrated by the production of XBP-1s. All together, these results support the notion that an interaction exists between NNV capsid protein and the UPR master regulator BiP, and that the interaction is associated with the ER stress response in the NNV-infected cells.
Figure 6
Figure 6

Role of NNV capsid protein in ER stress. (A) Co-immunoprecipitation of BiP protein and NNV capsid protein. GF-1 cells were infected with NNV at MOI of 5. The cell lysate was incubated with Protein A agarose and anti-BiP IgG. The immunocomplexes were separated on a 12% SDS-PAGE gel, and were further analyzed by western blotting with mouse anti-BiP and horseradish peroxidase-conjugated anti-mouse IgG. Chemiluminescence signal was captured after incubation with ECL substrate. For NNV capsid protein, the membrane was stripped, and re-probed with rabbit anti-NNV capsid antiserum, followed by horseradish peroxidase-conjugated anti-rabbit IgG. The result shows that BiP protein could interact with NNV capsid protein after NNV infection. (B) Induction of XBP1 mRNA splicing by NNV capsid protein in vivo. Zebrafish embryos at one-cell stage were mock-injected with 9.2 nl of sterile water or injected with 9.2 nl of pIR8 at the concentration of 150 ng/μl. The larvae were sacrificed at 20 hr post-fertilization and subjected to RNA extraction. RT- PCR was carried out to analyze the transcript of xbp-1 gene. Splicing of XBP1 mRNA, an indicator of ER stress, was detected in the embryos expressing exogenous NNV capsid protein.


Grouper is an economically important fish species, whose life cycle is severely threatened by many pathogens, particularly viruses. NNV causes devastative mortality in grouper at the larvae and juvenile stages, and yet the understanding of its pathogenecity is limited, mainly hampered by the lack of proper tools and fundamental genomic information of the host. The leap of the NGS technology, such as Illumina Genome Analyzer, is a valuable tool providing high-volume, parallel sequencing throughput, thus extremely beneficial to the study of non-model species of no reference genomic resources, such as grouper[26]. The technology has been widely applied in gene discovery, digital gene expression profiling and comparative genomics studies[2630]. In this study, we have employed the NGS technology to generate a transcriptome library of 66,582 unigenes from grouper GK cells. By using this library as a reference, we have identified in NNV-infected GK cells 117 unigenes representing genes that are known to be associated with protein processing in ER. The obtained sequence information will greatly assist in future studies on the molecular regulation of ER response in grouper cells under normal or pathological condition. Endoplasmic reticulum serves to maintain protein homeostasis, including the regulation of the concentration, conformation, folding, and trafficking of client protein in cells[31]. Many studies have revealed that the viral infection and replication of flavivirus family such as Japanese encephalitis virus, bovine viral diarrhea virus, and hepatitis C virus would induce ER dysfunction and stress[3235]. Some other viruses such as respiratory syncytial virus, African swine fever virus, herpes simplex virus and cytomegalovirus have also been shown to regulate ER stress[3639]. On the other hand, intracellular stress responses of UPR would induce apoptosis in cell if the ER stress is prolonged. Hence, the delicate modulation of UPR would determine the outcome of viral infection in host cells[40].

Recent studies have shown that ER stress plays a role in the apoptotic cell death caused by NNV infection[21]. However, the understanding of the molecular mechanism of NNV-induced ER stress response is limited. It has most recently been shown that BiP is up-regulated in grouper GF cells after NNV infection[21]. While the up-regulation of BiP after NNV infection is confirmed in this study, we have further shown that the three mediators (PERK, ATF6 and IRE1) controlled by BiP were up-regulated as well in the infected cells. Among the UPR-associated genes, CHOP and GADD34 were most prominently up-regulated at 6 hr post-infection. As both CHOP and GADD34 are components of the PERK pathway, the earlier up-regulation of these two genes suggests an earlier onset of the PERK pathway. This observation coincides with the notion that PERK functions at an earlier stage than ATF6 and IRE1 during ER stress response[17].

A surprising result of this study is that BiP protein may co-localize and interact with the capsid protein of NNV, as demonstrated by confocal microscopy and co-immunoprecipitation assays (Figures5 &6A). Furthermore, overexpression of the capsid protein would trigger ER stress response in vivo (Figure6B) and apoptotic cell death in the transfected cells (data not shown). Interestingly, it has been shown in a most recent study that BiP protein can interact with the RNA-dependent RNA polymerase (RdRp), and that the interaction of these two proteins may enhance mitochondria-mediated cell death[21]. Therefore, it appears that BiP protein may interact with at least two viral components (RdRp and capsid) in NNV-infected cells, and that the interaction between BiP and viral proteins is crucial to the outcome of NNV infection. Finally, we provide evidence that NNV capsid protein could enter the nucleus of infected cell during viral replication. This result indicates that NNV capsid protein might also play a role in the regulation of gene expression in the infected cells.


In summary, we have established a database of genes that are involved in the ER stress response in NNV-infected grouper cells. Additionally, we have discovered a potential role of NNV capsid protein in the ER stress response. Another interesting discovery is that NNV capsid protein might play a role in the regulation of gene expression in infected cells as suggested by the presence of capsid protein in the nucleus. These data provide new insights into the molecular mechanism of NNV-induced ER stress, and will be of great value to future studies of NNV pathogenicity.


Cells and virus

The grouper GK cells[41] and GF-1 cells[25], were grown at 28°C in Leibovitz L-15 medium supplemented with 5% and 10% fetal bovine serum (FBS) (Invitrogen), respectively. Malabar grouper nervous necrosis virus (MGNNV) was propagated and titrated in GF-1 cells as described previously[42]. For viral infection, confluent cells in 6-well plates (Corning) were rinsed with serum-free L15 and infected with the virus at multiplicity of infection (MOI) of 5 (for GF-1 cells) or 10 (for GK cells), unless mentioned otherwise.

RNA extraction

Total RNA was prepared by using Trizol® reagent (Invitrogen, USA) according to manufacturer’s instruction. Prior to Illumina sequencing, the quality of RNA was verified by Agilant 2100 Bioanalyzer (Agilent Technologies, USA) with a readout RIN value of 9.2, and the quantity of RNA was measured by using Nanodrop® ND-1000 spectrophotometer (LabTech, USA).

Transcriptomic sequencing and reads assembly

Solexa sequencing was performed as a commercial service in Beijing Genomic Institute (BGI). After collection of total RNA, beads coated with Oligo(dT) were used to isolate poly(A) mRNA from the samples. First-strand cDNA was synthesized with random hexamer-primers. The second-strand cDNA was synthesized and purified. After purification, end reparation and adding poly(A), the cDNA fragments were tailed with sequencing adaptors. These fragments were then subjected to PCR amplification, and the products were subsequently sequenced by using Illumina HiSeq 2000. After sequencing, the SOAP de novo software[23] was used to assemble the data. The software first combined reads and overlaps to form longer fragments, called contigs. The contigs were further assembled into scaffolds by joining the paired-end reads of contigs. Gap filling of scaffolds by paired-end-reads was applied to obtained sequences with least Ns and not extendable on either ends. These sequences are known as unigenes.

Unigenes annotation and functional classification

Homology searches of the assembled unigene sequences were conducted against public protein databases of Nr, Swiss-Prot, Kyoto Encyclopedia of Genes and Genomes (KEGG), and Cluster of Orthologous Groups (COG) by using blastx (the e-value was set to be < 0.00001). Highest sequence similarity of proteins were retrieved and assigned to specific unigene. Unigenes that cannot be aligned to any database were further scanned by ESTScan[43]. Functional annotation of unigenes was performed according to the KEGG Pathway, COG, and Gene Ontology (GO). GO was performed with non-redundant (nr) annotation by Blast2GO program[44].

Comparison of unigene expression among different DGE samples

The gene expression level was calculated by using RPKM method (Reads Per kb per Million reads)[24], according to the formula: RPKM(A) = 105C/(NL × 10-3). RPKM(A) stands for the expression of gene A, C for the number of reads that uniquely aligned to gene A, N for the total number of reads that uniquely aligned to all genes, and L for the number of bases on gene A. As the RPKM method can eliminate the influence of different gene length and sequencing discrepancy on the calculation of gene expression, the calculated gene expression can thus be directly used for comparing the difference of gene expression among samples. When there is more than one unigene for a gene, the longest one is used for the calculation of its expression level and coverage. The expression of unigenes in response to NNV infection was analyzed by comparing the RPKM of unigenes in the infected versus non-infected GK cells at 6 and 33 hours post-infection.

Identification and expression profile of ER-stress related genes in NNV-infected GK cells

Unigenes that are relevant to ER stress were identified by the analysis against the KEGG Pathway. Quantitative RT-PCR was carried out to quantitate the expression level of selective ER-stress related genes: CHOP, PERK, ATF6, BiP, XBP1 and GADD34. The RNA was reverse-transcribed into cDNA using Superscript III Reverse Transcriptase from Invitrogen by following the manufacturer’s recommendation. Real-time PCR was performed with the reagent Express SYBR GreenER qPCR Supermix Universal (Invitrogen) by following the manufacturer’s suggestions. The reaction was carried out with triplicates of each sample on a Rotor-Gene Q thermal cycler (Qiagen). Relative expression of each gene was normalized against the house-keeping gene β-actin. The PCR primers are listed in Additional file3: Table S1.

Sub-cellular localization of NNV capsid protein and BiP

GF-1 cells were infected with NNV at MOI of 10 for 24 hours. Infected cells were fixed with 4% formaldehyde at room temperature for 15 minutes. Fixed GF-1 cells were then permeabilized with PBS containing 0.2% triton X-100, followed by blocking with 3% bovine serum albumin (BSA) in PBS for 1 hour. Indirect immunofluorescence to stain NNV coat protein and BiP was performed with rabbit anti-NNV capsid antibody[45] at 1:1000 dilution and mouse anti-BiP antibody (Santa Cruz, sc-166490) at 1:50 dilution. The secondary antibodies against rabbit IgG conjugated with Texas Red (Invitrogen, T2767) or antibodies against mouse IgG conjugated with Fluorscein (KPL, Cat.02-18-06) were diluted at 1:500. All antibodies were diluted in PBS containing 3% BSA. 4',6-diamidino-2-phenylindole (DAPI) was further applied to stained cells for nuclear visualization. Stained cells were mounted onto glass slides with PBS containing 0.1g/ml 1,4-Diazabicyclo-[2]octane (DABCO) (Sigma, D-2522), 50% glycerol, and 0.1% Na3N. Confocal images were acquired with Olympus FV1000 inverted confocal microscope.

Co-immunoprecipitation assay

For co-immunoprecipitation (co-IP) assay, GF-1 cells at 90% confluency were infected with NNV (MOI = 5). The cells were harvested by trypsinization, washed three times with PBS, and lysed with cell lysis buffer (Cell Signaling, 9803S) at 0, 24, 48 and 72 hours post-infection. The cell lysate was collected after centrifugation at 14000 rpm for 20 min. The supernatant was precleared with Protein G-agarose beads (Sigma-Aldrich), and incubated overnight at 4°C with Protein A agarose (Invitrogen) and 2 μg anti-BiP IgG (Santa Cruz, sc-166490). The immunocomplexes were recovered by a brief centrifugation prior to separation on a 12% SDS–PAGE gel, and were subsequently transferred onto a nitrocellulose membrane. For BiP detection, the membranes were blocked with PBS-T containing 3% (w/v) non-fat milk (non-fat milk/PBS-T) for 1 h at room temperature. After washing, the membrane was incubated with mouse anti-BiP (1:1000, Santa Cruz) in 1% (w/v) non-fat milk/PBS-T overnight at 4°C. After blocking, the membrane was washed, and incubated with horseradish peroxidase-conjugated anti-mouse IgG (1:3000 in 0.1% PBS-T)(Zymed, 81–6520) for 1 h at room temperature. The membrane was incubated with ECL substrate (Pierce ECL Western Blotting Substrate, Thermo Scientific, 32209) for 1 minute, and visualized by a chemiluminescence image system (UVP, BioSpectrum Imaging System 500). For NNV capsid detection, the membrane was stripped with 20 ml of stripping buffer (Restore Western Blot Stripping Buffere, Thermo Scientific, 21059) at room temperature for 15 min, blocked with 3% (w/v) non-milk fat PBS-T at room temperature for 1 h, and re-probed with rabbit anti-NNV capsid antibody (1:3000 in 0.1% PBS-T)[46] overnight at 4°C. The membrane was incubated with secondary antibody horseradish peroxidase-conjugated anti-rabbit IgG (1:3000, PerkinElmar, NEF812001EA), and subsequently the ECL substrate.

Plasmid and transfection

The NNV capsid protein gene was cloned into a eukaryotic expression vector pIRES-e2EGFP (Clontech) at the Eco RI-Sal I sites of the vector to generate a plasmid termed pIR8. The plasmid was transformed into JM109 cells and the transformed cells were grown in large quantity to isolate the plasmid. Isolation and purification of the pIR8 plasmid was conducted by using QIAprep Spin Miniprep Kit (Qiagen). Transfection of pIR8 was carried out with Fugene® 6 Transfection Reagent (Roche) by following the manufacturer’s instruction.

Induction of XBP1 mRNA splicing by NNV capsid protein in vivo

Zebrafish embryos at one-cell stage were mock-injected or injected with 9.2 nl of pIR8 at the concentration of 150 ng/μl by using a microinjector (Leica EZ4 nanoliter 2000). The larvae were sacrificed at 20 hr post-fertilization and subjected to RNA extraction. The RNA was reversetranscribed to cDNA, followed by PCR analysis. The XBP-1 primers for the PCR are:5’-GTT CAG GTA CTG GAG TCC GC-3’ (forward primer) and 5’-GGA TGT CCA GAA TAC CAA GCA GG-3’ (reversed primer). The PCR product was analyzed by electrophoresis with a 3% agarose gel (UltraPure TM Agarose-1000, Invitrogen).



This study was partially supported by National Science Council (NSC98-2313-B019-006; NSC98-2324-B019-003-CC1) to MW Lu. This study is also partially supported by an Academia Sinica intramural funding to PP Chiou. We thank Dr. Thomas Chen at the University of Connecticut for reviewing this manuscript and discussion.

Authors’ Affiliations

Department of Aquaculture, National Taiwan Ocean University, Keelung, Taiwan
Center of Excellence for Marine Bioenvironment and Biotechnology, National Taiwan Ocean University, Keelung, Taiwan
Institute of Basic Medical Science, National Cheng Kung University, Tainan, Taiwan
Institute of Biotechnology, National Ilan University, Yilan, Taiwan
Marine Research Station, Institute of Cellular and Organismic Biology, Academia Sinica, Yilan, Taiwan


  1. Su M-S, Su W-C: The status and prospects of coastal aquaculture in Taiwan. Edited by: Institute FR. 2010, Keelung: Council of Aquaculture, 1-12.Google Scholar
  2. Yoshikoshi K, Inoue K: Viral nervous necrosis in hatchery-reared larvae and juveniles of Japanese parrotfish, Oplegnathus fasciatus (Temminck & Schlegel). J Fish Dis. 1990, 13: 69-77. 10.1111/j.1365-2761.1990.tb00758.x.View ArticleGoogle Scholar
  3. Mori K, Nakai T, Nagahara M, Muroga K, Mekuchi T, Kanno T: A viral disease in hatchery-reared larvae and juveniles of redspotted grouper. J Pathol. 1991, 26: 209-210.Google Scholar
  4. Mori K-i, Nakai T, Muroga K, Arimoto M, Mushiake K, Furusawa I: Properties of a new virus belonging to nodaviridae found in larval striped jack (Pseudocaranx dentex) with nervous necrosis. J Virol. 1992, 187: 368-371. 10.1016/0042-6822(92)90329-N.View ArticleGoogle Scholar
  5. Munday BL, Kwang J, Moody N: Betanodavirus infections of teleost fish: a review. J Fish Dis. 2002, 25: 127-142. 10.1046/j.1365-2761.2002.00350.x.View ArticleGoogle Scholar
  6. Nguyen HD, Nakai T, Muroga K: Progression of striped jack nervous necrosis virus (SJNNV) infection in naturally and experimentally infected stripted jack Pseudocaranx dentex larvae. Dis Aquat Organ. 1996, 24: 99-105.View ArticleGoogle Scholar
  7. Muroga K: Viral and bacterial diseases of marine fish and shellfish in Japansese hatcheries. J Aquaculture. 2001, 202: 23-44. 10.1016/S0044-8486(01)00597-X.View ArticleGoogle Scholar
  8. Schneemann A, Ball L, Delsert C, Johnson J, Nishizawa T: Family Nodaviridae. 2004, London: AcademicGoogle Scholar
  9. Chen L-J, Su Y-C, Hong J-R: Betanodavirus non-structural protein B1: A novel anti-necrotic death factor that modulates cell death in early replication cycle in fish cells. Virology. 2009, 385 (2): 444-454. 10.1016/j.virol.2008.11.048.View ArticlePubMedGoogle Scholar
  10. Fenner BJ, Thiagarajan R, Chua HK, Kwang J: Betanodavirus B2 Is an RNA Interference Antagonist That Facilitates Intracellular Viral RNA Accumulation. J Virol. 2005, 80 (1): 85-94.View ArticleGoogle Scholar
  11. Iwamoto T, Mise K, Takeda A, Okinaka Y, Mori K-I, Arimoto M, Okuno T, Nakai T: Characterization of Striped jack nervous necrosis virus subgenomic RNA3 and biological activities of its encoded protein B2. J Gen Virol. 2005, 86 (10): 2807-2816. 10.1099/vir.0.80902-0.View ArticlePubMedGoogle Scholar
  12. Glazebrook JS, Heasman MP, Beer SW: Picorna-like viral particles associated with mass mortalities in larval barramundi, Lates calcarifer Bloch. J Fish Dis. 1990, 13 (3): 245-249. 10.1111/j.1365-2761.1990.tb00780.x.View ArticleGoogle Scholar
  13. Bloch B, Gravningen K, Larsen JL: Encephalomyelitis among turbot associated with picornavirus-like agent. Dis Aquat Organ. 1991, 10: 65-70.View ArticleGoogle Scholar
  14. Breuil G, Bonami JR, Pepin JF, Pichot Y: Viral infection (picorna-like virus) associated with mass mortalities in hatchery-reared sea-bass (Dicentrarchus labrax) larvae and juveniles. Aquaculture. 1991, 97: 109-166. 10.1016/0044-8486(91)90258-9.View ArticleGoogle Scholar
  15. Boonyaratpalin S, Supamattaya K, Kasornchandra J, Hoffman RW: Picorna-like virus associated with mortality and a spongiosus encephalopathy in grouper Epinephelus malabaricus. Dis Aquat Organ. 1996, 26: 75-80.View ArticleGoogle Scholar
  16. Grotmol S, Totland GK, Kryvi H: Detection of a nodavirus-like agent in heart tissue from reared Atlantic salmon Salmo salar suffering from cardiac myopathy syndrome (CMS). Dis Aquat Organ. 1997, 29: 79-84.View ArticleGoogle Scholar
  17. He B: Viruses, endoplasmic reticulum stress, and interferon responses. Cell Death Differ. 2006, 13 (3): 393-403. 10.1038/sj.cdd.4401833.View ArticlePubMedGoogle Scholar
  18. Benali-Furet NL, Chami M, Houel L, De Giorgi F, Vernejoul F, Lagorce D, Buscail L, Bartenschlager R, Ichas F, Rizzuto R, et al: Hepatitis C virus core triggers apoptosis in liver cells by inducing ER stress and ER calcium depletion. Oncogene. 2005, 24 (31): 4921-4933. 10.1038/sj.onc.1208673.View ArticlePubMedGoogle Scholar
  19. Li X, Lankinen H, Putkuri N, Vapalahti O, Vaheri A: Tula hantavirus triggers pro-apoptotic signals of ER stress in Vero E6 cells. J Virol. 2005, 333: 180-189. 10.1016/j.virol.2005.01.002.View ArticleGoogle Scholar
  20. Su HL, Liao CL, Lin YL: Japanese Encephalitis Virus Infection Initiates Endoplasmic Reticulum Stress and an Unfolded Protein Response. J Virol. 2002, 76 (9): 4162-4171. 10.1128/JVI.76.9.4162-4171.2002.PubMed CentralView ArticlePubMedGoogle Scholar
  21. Su Y-C, Wu J-L, Hong J-R: Betanodavirus up-regulates chaperone GRP78 via ER stress: roles of GRP78 in viral replication and host mitochondria-mediated cell death. Apoptosis. 2011, 16 (3): 272-287. 10.1007/s10495-010-0565-x.View ArticlePubMedGoogle Scholar
  22. Gregory R, Darby AC, Irving H, Coulibaly MB, Hughes M, Koekemoer LL, Coetzee M, Ranson H, Hemingway J, et al: A De Novo Expression Profiling of Anopheles funestus, Malaria Vector in Africa, Using 454 Pyrosequencing. PLoS ONE. 2011, 6 (2): e17418-10.1371/journal.pone.0017418.PubMed CentralView ArticlePubMedGoogle Scholar
  23. Li R, Zhu H: De novo assembly of human genomes with massively parallel short read sequencing. Genome Res. 2010, 20 (2): 265-272. 10.1101/gr.097261.109.PubMed CentralView ArticlePubMedGoogle Scholar
  24. Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B: Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Meth. 2008, 5 (7): 621-628. 10.1038/nmeth.1226.View ArticleGoogle Scholar
  25. Chi S-C, Lin S-C, Su H-M, Hu W-W: Temperature effect on nervous necrosis virus infection in grouper cell line and in grouper larvae. Virus Res. 1999, 63 (1–2): 107-114.View ArticlePubMedGoogle Scholar
  26. Wang X-W, Luan J-B, Li J-M, Bao Y-Y, Zhang C-X, Liu S-S: De novo characterization of a whitefly transcriptome and analysis of its gene expression during development. BMC Genomics. 2010, 11 (1): 400-10.1186/1471-2164-11-400.PubMed CentralView ArticlePubMedGoogle Scholar
  27. Rosenkranz R, Borodina T, Lehrach H, Himmelbauer H: Characterizing the mouse ES cell transcriptome with Illumina sequencing. Genomics. 2008, 92 (4): 187-194. 10.1016/j.ygeno.2008.05.011.View ArticlePubMedGoogle Scholar
  28. Collins LJ, Biggs PJ, Voelckel C, Joly S: An approach to transcriptome analysis of non-model organism using short-read sequences. Genome Infomatics. 2008, 21: 3-14.Google Scholar
  29. Hegedus Z, Zakrzewska A, Agoston V, Ordas A, Racz P, Mink M, Spaink H, Meijer A: Deep sequencing of the zebrafish transcriptome response to mycobacterium infection. Mol Immunol. 2009, 46 (15): 2918-2930. 10.1016/j.molimm.2009.07.002.View ArticlePubMedGoogle Scholar
  30. Richardson PM, Nowrousian M, Stajich JE, Chu M, Engh I, Espagne E, Halliday K, Kamerewerd J, Kempken F, Knab B, et al: De novo Assembly of a 40 Mb Eukaryotic Genome from Short Sequence Reads: Sordaria macrospora, a Model Organism for Fungal Morphogenesis. PLoS Genet. 2010, 6 (4): e1000891-10.1371/journal.pgen.1000891.View ArticleGoogle Scholar
  31. Inagi R: Endoplasmic Reticulum Stress as a Target of Therapy Against Oxidative Stress and Hypoxia. Studies on Renal Disorders. Edited by: Miyata T, Eckardt K-U, Nangaku M. 2011, New York: Humana Press, 657-672. 672View ArticleGoogle Scholar
  32. Wu Y-P, Chang C-M, Hung C-Y, Tsai M-C, Schuyler SC, Wang R: Japanese encephalitis virus co-opts the ER-stress response protein GRP78 for viral infectivity. Virol J. 2011, 8 (1): 128-10.1186/1743-422X-8-128.PubMed CentralView ArticlePubMedGoogle Scholar
  33. Tardif K, Waris G, Siddiqui A: Hepatitis C virus, ER stress and oxidative stress. Trend Microbiol. 2005, 13: 159-163. 10.1016/j.tim.2005.02.004.View ArticleGoogle Scholar
  34. Jordan R, Wang L, Graczyk T, Block T, Romano P: Replication of a cytopathic strain of bovine viral diarrhea virus activates PERK and induces endoplasmic reticulum stress-mediated apoptosis od MDBK cells. J Virol. 2002, 76: 9588-9599. 10.1128/JVI.76.19.9588-9599.2002.PubMed CentralView ArticlePubMedGoogle Scholar
  35. Su H, Liao C, Lin Y: Japanese encephalitis virus infection initiates endoplasmic reticulum stress and an unfolded protein response. J Virol. 2002, 76: 4162-4171. 10.1128/JVI.76.9.4162-4171.2002.PubMed CentralView ArticlePubMedGoogle Scholar
  36. Netherton C, Parsley J, Wileman T: African swine fever virus inhibits induction of the stress-induced proapoptotic transcription factor CHOP/GADD153. J Virol. 2004, 78: 10825-10828. 10.1128/JVI.78.19.10825-10828.2004.PubMed CentralView ArticlePubMedGoogle Scholar
  37. Bitko V, Barik S: An endoplasmic reticulum-specific stress-activated caspase (caspase-12) is implicated in the apoptosis of A549 epithelial cells by respiratory syncytial virus. J Cell Biochem. 2001, 80: 441-454. 10.1002/1097-4644(20010301)80:3<441::AID-JCB170>3.0.CO;2-C.View ArticlePubMedGoogle Scholar
  38. Cheng G, Feng Z, He B: Herpes simplex virus 1 infection activates the endoplasmic reticulum resident kinase PERK and mediates elF-2α dephosphorylation by the γ134.5 protein. J Virol. 2005, 79: 1379-1388. 10.1128/JVI.79.3.1379-1388.2005.PubMed CentralView ArticlePubMedGoogle Scholar
  39. Isler J, Skalet A, Alwine J: Human cytomegalovirus infection activates and regulates the unfolded protein response. J Virol. 2005, 79: 6890-6899. 10.1128/JVI.79.11.6890-6899.2005.PubMed CentralView ArticlePubMedGoogle Scholar
  40. Umareddy I, Pluquet O, Wang Q, Vasudevan SG, Chevet E, Gu F: Dengue virus serotype infection specifies the activation of the unfolded protein response. Virol J. 2007, 4 (1): 91-10.1186/1743-422X-4-91.PubMed CentralView ArticlePubMedGoogle Scholar
  41. Lai YS, John JAC, Lin CH, Guo IC, Chen SC, Fang K, Chang CY: Establishment of cell lines from a tropical grouper, Epinephelus awoara (Temminck & Schlegel), and their susceptibility to grouper irido- and nodaviruses. J Fish Dis. 2003, 26 (1): 31-42. 10.1046/j.1365-2761.2003.00434.x.View ArticlePubMedGoogle Scholar
  42. Lu M-W, Chao Y-M, Guo T-C, Santi N, Evensen Ø, Kasani SK, Hong J-R, Wu J-L: The interferon response is involved in nervous necrosis virus acute and persistent infection in zebrafish infection model. Mol Immunol. 2008, 45 (4): 1146-1152. 10.1016/j.molimm.2007.07.018.View ArticlePubMedGoogle Scholar
  43. Iseli C, Jongeneel CV: ESTScan: a program for detecting, evaluating, and reconstructing potential coding regions in EST sequences. Proc Int Conf Intell Syst Mol Bio. 1999, 138-148.Google Scholar
  44. Conesa A, Gotz S: Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005, 21 (18): 3674-3676. 10.1093/bioinformatics/bti610.View ArticlePubMedGoogle Scholar
  45. Lu M-W, Liu W, Lin C-S: Infection competition against grouper nervous necrosis virus by virus-like particles produced in Escherichia coli. J Gen Virol. 2003, 84 (6): 1577-1582. 10.1099/vir.0.18649-0.View ArticlePubMedGoogle Scholar
  46. Lu MW: Infection competition against grouper nervous necrosis virus by virus-like particles produced in Escherichia coli. J Gen Virol. 2003, 84 (6): 1577-1582. 10.1099/vir.0.18649-0.View ArticlePubMedGoogle Scholar


© Lu et al.; licensee BioMed Central Ltd. 2012

This article is published under license to BioMed Central Ltd. 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 work is properly cited.