Figures
Abstract
Pacific white shrimp (Litopenaeus vannamei) is the most extensively farmed crustacean species in the world. White spot syndrome virus (WSSV) is one of the major pathogens in the cultured shrimp. However, the molecular mechanisms of the host-virus interaction remain largely unknown. In this study, the impact of WSSV infection on host gene expression in the hepatopancreas of L. vannamei was investigated through the use of 454 pyrosequencing-based RNA-Seq of cDNA libraries developed from WSSV-challenged shrimp or normal controls. By comparing the two cDNA libraries, we show that 767 host genes are significantly up-regulated and 729 genes are significantly down-regulated by WSSV infection. KEGG analysis of the differentially expressed genes indicated that the distribution of gene pathways between the up- and down-regulated genes is quite different. Among the differentially expressed genes, several are found to be involved in various processes of animal defense against pathogens such as apoptosis, mitogen-activated protein kinase (MAPK) signaling, toll-like receptor (TLR) signaling, Wnt signaling and antigen processing and presentation pathways. The present study provides valuable information on differential expression of L. vannamei genes following WSSV infection and improves our current understanding of this host-virus interaction. In addition, the large number of transcripts obtained in this study provides a strong basis for future genomic research on shrimp.
Citation: Chen X, Zeng D, Chen X, Xie D, Zhao Y, Yang C, et al. (2013) Transcriptome Analysis of Litopenaeus vannamei in Response to White Spot Syndrome Virus Infection. PLoS ONE 8(8): e73218. https://doi.org/10.1371/journal.pone.0073218
Editor: Shoba Ranganathan, Macquarie University, Australia
Received: December 9, 2012; Accepted: July 18, 2013; Published: August 26, 2013
Copyright: © 2013 Chen 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 supported by the National Natural Science Foundation of China (31160531) and the “Twelfth Five-Year” National Science and Technology Project in Rural Areas (#2012AA10A404-3). 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
Pacific white shrimp (Litopenaeus vannamei) has become the most extensively cultured crustacean species in the world due to its high-yield and low-demand for concentration of salt [1]. Over the last two decades, viral diseases have seriously threatened the shrimp aquaculture industry [2]–[5]. Among the viral diseases, white spot syndrome virus (WSSV) is the most important viral pathogens of shrimp farming, causing high mortality and large economic loss [6]. WSSV is an enveloped virus with a large, double stranded, circular genome (∼300 kb) [7], [8] and has a wide host range, including many aquatic crustaceans and even aquatic insects [9], [10]. In cultured shrimp, WSSV causes a cumulative mortality can reach 100% within 3–10 days [4]. To date, no effective cure for this disease has yet been found. As a result, the response of shrimp to viral pathogens has become an increasingly important subject for study [11]. An understanding of the molecular responses and defense mechanism in shrimp against pathogens will promote the development of effective techniques to prevent viral diseases of shrimp.
In recent years, some attempts have been made to investigate the effects of viral infection on shrimp transcriptome using the suppression subtractive hybridization (SSH) or the cDNA microarray technology [12]–[14]. Several immune-related proteins in shrimp, such as hemocyanin [15], [16], lectins [17], [18], antimicrobial peptide [19] and superoxide dismutase [20] have been identified. However, such studies have been hampered by a relative lack of molecular genetic information and tools suitable for high throughput assessment of gene expression [11].
Recently, several next-generation sequencing technologies, such as the 454 Life Sciences (Roche) pyrosequencing platform, the Illumina Genome Analyzer and the Applied Biosystems Solid platform, provide a rapid, cost-effective and high-throughput tool for sequencing of the transcriptome (RNA-Seq) [21], [22]. The RNA-Seq permits deep, robust assessments of transcript abundances and transcript structure [23]. When gene expression profiling is applied to the interaction between WSSV and shrimp, insights into the mechanisms this virus use to suppress and subvert the host defense response can be made.
In this study, we used 454 pyrosequencing to conduct a comparative transcriptome profiling analysis between the hepatopancreas of WSSV-infected L. vannamei shrimp and uninfected controls, with the aim of investigating candidate immune-related genes in shrimp, improving the current understanding of the host-virus interaction, and providing a substantial dataset that increases the publicly available DNA sequence resources for this crustacean species.
Materials and Methods
Experimental shrimp and WSSV challenge
The experiment was conducted using a specific pathogen-free (SPF) shrimp L. vannamei strain (National and Guangxi Shrimp Genetic Breeding Center, Guangxi Province, China). The shrimp (11–12 g body weight) were maintained in the environmentally controlled 1000-liter glass salt water tanks (32-ppt salinity, 25 to 26°C) and fed an artificial pellet feed. In the challenge experiment, a batch of shrimp was divided into a challenge group and a normal control group. The challenge group contained 40 shrimp, and the control group contained 20 shrimp. Shrimp in the challenge group were fed once a day for 3 consecutive days with minced WSSV-infected tail tissue at 5% of their body weight (a dose of about 1×105 WSSV copies/g sufficient to cause 100% mortality in 5–7 days). The WSSV dose was determined in a previous experiment, in which RT-PCR method and reference samples containing gradient concentrations of WSSV were used. The WSSV strain was isolated by our laboratory in 2008 from a L. vannamei shrimp in China. In the challenge experiment, the mortality rate of challenge group was 0, 5% and 32.5% at 24 h, 48 h and 72 h, respectively, after WSSV challenge. In parallel, shrimp in the control group were fed once a day for 3 consecutive days with minced PCR-confirmed [24] healthy tail tissue at 5% of their body weight [25]. The hepatopancreas tissues of shrimp (20 shrimp each group) were taken at 72 h post challenge and stored in liquid nitrogen (−196°C) until RNA isolation.
RNA extraction, cDNA library construction and RNA-seq
Total RNA was extracted using TRIzol reagent (Qiagen) following the manufacturer's instructions. RNA concentrations were measured using a spectrophotometer, and integrity was ensured through the analysis on a 1.5% (w/v) agarose gel. After RNA extraction, mRNAs were purified using the PolyATtract® mRNA Isolation Systems (Promega) and concentrated using the Reasy RNA Cleaning Kit (Qiagen). Equal amounts of the high-quality mRNA samples from each group were then pooled for cDNA synthesis and sequencing. The pooled mRNAs were fragmented into small pieces using RNA fragment reagent (Qiagen), and the pieces were collected using the Reasy RNA Cleaning Kit (Qiagen). Subsequently, cleaved RNA fragments were copied into first strand cDNA using MMLV reverse transcriptase and random primers. This was followed by second strand cDNA synthesis using DNA polymerase I and RNase H. Library construction and a 1/2 plate pyrosequencing run was performed by Beijing Autolab Biotechnology Co., Ltd.(China) on a 454 GS FLX system (Roche) according to the manufacturer's instructions.
De novo assembly and gene annotation
Before the assembly, the raw sequencing reads were qualities trimmed and adaptor sequences removed using the SeqClean program (http://compbio.dfci.harvard.edu/tgi/software). The screened high-quality sequences (cleaned reads) were de novo assembled using the iAssembler program (http://bioinfo.bti.cornell.edu/tool/iAssembler) with default parameters [26]. The overall assembly was performed using the combined sequence data from both the WSSV-infected and the control samples. The contigs and singletons were generally referred to as unigenes. Open reading frames (ORF) of each unigene were identified by getorf (http://emboss.bioinformatics.nl/cgi-bin/emboss/getorf). For functional annotation analysis, all unigenes were compared against sequences in NCBI non-redundant (nr) protein and UniProtKB/Swiss-Prot (UniProt release 2013_07 - Jun 26, 2013) [27] database using the BLASTX programs (E-value <10−5) [28]. Genes were tentatively identified according to the best hits against known sequences. Functional annotation by gene ontology terms (GO) was analyzed using a BLAST2GO program (http://www.BLAST2go.org/). The COG and KEGG pathway annotations of unigenes was performed using the BLASTX software against the COG database and the KEGG database, respectively [29], [30].
Identification of differentially expressed genes
For differential gene expression analyses, the transcript levels were measured as RPKM (Reads Per Kilobase of exon model per Million mapped reads) values to determine relative transcript abundance. Statistical comparison between two different libraries was conducted using a web tool IDEG6 (http://telethon.bio.unipd.it/bioinfo/IDEG6_form/) [31]. FDR (false discovery rate) <0.001 was used as the threshold of p-value in multiple test to judge the significance of gene expression difference [32]. Genes were considered differentially expressed in a given library when the p-value was less than 0.001, and a greater than two-fold change in expression across libraries was observed.
Quantitative RT-PCR analysis
To validate our 454 sequencing data, eight genes including hemocyanin (RL-all_rep_c8318), glycine N-methyltransferase (RL-all_rep_c13364), nucleoside diphosphate kinase (RL-all_rep_c11491), lectin (RL-all_rep_c12149), cathepsin L (RL-all_c1576), trypsin (RL-all_c1114), zinc proteinase (RL-all_c2336), and acidic mammalian chitinase (RL-all_rep_c12797) were selected randomly for quantitative RT-PCR analysis, using the same RNA samples as for the 454 sequencing. First strand cDNA was synthesized from 1 µg of RNA using M-MuLV reverse transcriptase (Qiagen). The qPCR reaction mixture (20 µL) consisted of 2× Power SYBR Green PCR Master mix, 0.9 M each of the forward and reverse primers, and 1 µL of template cDNA. PCR amplification was performed under the following conditions: 50°C for 2 min and 95°C for 30 s, followed by 40 cycles of 95°C for 15 s and 62°C for 1 min, and a final extension at 72°C for 5 min. Primer sets were designed using the Primer Premier 5.0 software (Table S1), and the parallel reactions using beta-actin primers [33] were performed to normalize the amount of template cDNA added to each reaction.
Results and Discussion
Sequencing and de novo assembly
As described in the Methods, cDNA libraries from pooled mRNAs extracted from the hepatopancreas of twenty WSSV-challenged shrimp or of twenty control shrimp were subjected to a 1/2-plate production run on the 454 GS-FLX sequencing instrument, resulting in 159,742 and 101,725 raw reads, respectively. Files containing these data were deposited in the Short Read Archive (SRA) database of the National Center for Biotechnology Information (NCBI) with accession numbers of SRR554363 (WSSV- challenged) and SRR554364 (control). After removal of repetitive, low-complexity, and low-quality reads, a total of 151023 cleaned reads, ranging from 70 bp to 611 bp were harvested in the WSSV-challenged library, with an average length of 374 bp and N50 length of 451 bp. In the control library, a total of 95372 cleaned reads, ranging from 70 bp to 620 bp were obtained, with an average length of 370 bp and N50 length of 444 bp (Table 1). The overall de novo assembly was performed using the combined reads from the two libraries, generating a total of 14538 unigenes. The length distribution of these unigenes is shown in Table 2. The lengths of the unigenes range from 100 to 8750 bp, with the average length of 574 bp, a greater number than that obtained by Illumina sequencing in L. vannamei [34]. The relatively long read lengths allow for greater statistical confidence in the de novo assembly [35], [36].
Similarity searches
For functional annotation, all unigenes were subjected to BLASTX similarity search against the NCBI non-redundant (NR) protein databases with a cutoff E-value <10−5. Among the 14538 unigenes from both the WSSV-infected library and the control library, 4605 (31.68%) showed significant matches in the UniProtKB/Swiss-Prot database. An additional 10625 (73.08%) unigenes showed significant matches in the NR database. In total, 10648 (73.24%) unigenes were annotated in UniProtKB/Swiss-Prot or NR database. The species distribution of the best match result for each sequence is shown in Figure 1. The L. vannamei unigenes showed 13.68% matches with Pediculus humanus corporis sequences followed by Bos taurus (9.93%) and Nematostella vectensis (9.40%). Unexpectedly, only 5.19% unigenes were found to be similar to the proteome of L. vannamei. This might be attributed to the reason that the amount of the L. vannamei protein sequences available in the NR protein databases is relatively limited, and the number of matched sequences does not completely reflect the phylogenetic relationships.
The figure shows the species distribution of unigene BLASTX results against the NCBI non-redundant protein databases with a cutoff E value <10−5. Different colors represent different species. Species with proportions of more than 1% are shown.
Transcriptome comparison between WSSV-infected and uninfected shrimp
GO is an international standardized gene functional classification system which offers a dynamic-updated controlled vocabulary and a strictly defined concept to comprehensively describe properties of genes and their products in any organis [37]. To obtain an overview of the functional categories represented by the L. vannamei transcriptome, we used the GO assignment to classify the functions of the unigenes.Comparison of the transcriptome pattern of the hepatopancreas of WSSV-infected shrimp for GO terms with those of the hepatopancreas of uninfected shrimp showed high similarity in the distribution of unigenes across GO categories. In total, 7325 (45.4%) and 6527 (45.5%) unigenes map to biological processes, 5672 (35.2%) and 5108 (35.6%) unigenes map to cellular components, and 3129 (19.4%) and 2725 (20.0%) unigenes map to molecular functions in the WSSV-infected and uninfected shrimp, respectively (Figure 2). In both infected and uninfected shrimp, most of the corresponding biological process genes were involved in cellular processes, metabolic processes and localization; most of the cellular component genes encode proteins associated with parts of cells and cell organelles; most of the molecular function genes were associated with catalytic activity, binding and molecular transducer activity. The result showed that the functional distribution of the genes of the two shrimp groups was similar, indicating that the number of genes expressed in each GO category was not significantly affected by WSSV-infection.
The x axis shows subgroups of molecular functions from GO classification and the y axis shows the number of the matched unigenes.
Identification of differentially expressed genes
The significant differences in the expression of different unigenes were determined by comparing the relative transcript abundance of each unigene between the WSSV-infected shrimp and uninfected shrimp. We identified 767 host genes that are significantly up-regulated and 729 genes that are significantly down-regulated by infection with WSSV (Table S2). For functional annotation, all the differentially expressed unigenes were searched against NR, Nt, UniProtKB/Swiss-Prot, GO, COG, and KEGG databases by BLASTX with a cut-off E-value of 10−5 (Table S2). By this method, 1208 (80.7%) returned an above cut-off BLAST result. Of them, 1206 unigenes were annotated by NR (80.6%), and 819 (54.7%), 840 (56.1%), 534 (35.7%) and 758 (50.7%) unigenes by Nt, UniProtKB/Swiss-Prot, GO and COG, respectively. All differentially expressed unigenes with their Nr, Nt, UniProtKB/Swiss-Prot, COG, GO and ORF analysis are presented in additional Table S2. It is noteworthy that of the 1496 differentially expressed unigenes (up- and down-regulated), 288 (19.3%) unigenes have low sequence homology to known sequences in public databases, suggesting that they may represent non-coding RNA, misassembled contigs or unknown genes of L. vannamei involved in response to WSSV-infection.
To validate our RNA-seq results, eight unigenes were selected randomly for quantitative real time-PCR (qRT-PCR) analysis. The qRT-PCR results confirmed the data obtained from 454 sequencing analysis showing similar trends in up- or down-regulation of host genes (Figure 3). For example, based on 454 sequencing analysis, hemocyanin, glycine N-methyltransferase, nucleoside diphosphate kinase and lectin were up-regulated 1.94, 2.74, 2.71 and 2.87 log2-fold, respectively, and showed 3.50, 4.12, 2.16 and 3.52 log2-fold changes, respectively in qRT-PCR analyses (Figure 3). Likewise, in 454 sequencing analysis, cathepsin L, trypsin, zinc proteinase and acidic mammalian chitinase were down-regulated −1.04, −1.69, −1.18 and −1.20 log2-fold, respectively, and showed −1.85, −1.35, −2.36 and −1.05 log2-fold changes, respectively in qRT-PCR analyses. Although the results of qRT-PCR do not match the 454 sequencing analysis perfectly, which might be caused by the sequencing biases [6], [23], the qRT-PCR analysis confirmed the change direction detected by the 454 sequencing analysis.
Target gene abbreviations are as follows: HCY - hemocyanin, GNMT - glycine N-methyltransferase, NDK - nucleoside diphosphate kinase, LT - lectin, CTSL - cathepsin L, TPS - trypsin, ZP - zinc proteinase, AMC - acidic mammalian chitinase.
KEGG analysis of the differentially expressed genes
KEGG is a pathway-based categorization of orthologous genes that provides useful information for predicting functional profiles of genes [38]. To analyze the interplay between the up- and down-regulated genes, all the differentially expressed genes were mapped to the referential canonical pathways in the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. The top 30 pathways are presented in Figure 4. The distribution of gene pathways between the up- and down-regulated genes is quite different. Genes related to ‘Metabolic pathways’, ‘Phagosome’, ‘Focal adhesion’, ‘Tight junction’, ‘Adherens junction’ and ‘Biosynthesis of secondary metabolites’ were concentrated in the up-regulated gene cohort. Among these, more than 10 fold as many genes related to ‘Adherens junction’ and ‘Biosynthesis of secondary metabolites’ were present in the up-regulated gene group as compared to the group of down-regulated genes. While genes related to ‘Metabolic pathways’, ‘Phagosome’, ‘Focal adhesion’, ‘Tight junction’, ‘ABC transporters’, ‘Betalain biosynthesis’, ‘Tyrosine metabolism’, ‘Riboflavin metabolism’, ‘Isoquinoline alkaloid biosynthesis’, and ‘beta-Alanine metabolism’ were more suppressed in the presence of the viruses. Among these, the number of genes related to ‘ABC transporters’, ‘Betalain biosynthesis’ and ‘beta-Alanine metabolism’ were much larger (9–103 fold) in the down-regulated gene group as compared to the up-regulated genes group. Interestingly, analysis of the distribution of genes by pathways defined by the KEGG revealed clustering of the differentially expressed genes in cell junction pathways. Cell junctions include ‘focal adhesions’, ‘tight junctions’, ‘gap junctions’, ‘adherens junctions’, ‘desmosomes’, and ‘hemidesmosomes’, and each one plays a specific role in adhesion, communication, and transduction of mechanical force. Cell junctions regulate many normal biological processes, such as cell and metabolite transport across epithelial layers, development and differentiation programs, excitable cell contraction, immune responses, and neural activity [39]. Our results suggest that cell junctions might play an important role in the immune responses of shrimp against viral infection.
The top 30 pathways are shown. The x axis shows pathways from KEGG classification and the y axis shows the frequency of the matched unigenes.
Candidate genes involved in L. vannamei immune response
Among the genes that were found to be differentially expressed in the WSSV-infected shrimp compared to the uninfected controls, several are involved in various processes of animal defense against pathogens such as apoptosis, MAPK signaling, toll-like receptor, Wnt signaling and antigen processing and presentation (Table 3). Apoptosis pathway is known to be related to the cell hypersensitivity response, blocking pathogen progression and systemic resistance [40], [41]. Apoptosis has been reported to play an essential role in animal defense against viral pathogens. In crustaceans, the occurrence of apoptosis upon viral infection has long been observed [42]. Upon WSSV infection, apoptosis has been detected in several viral target tissues of shrimp [43]. Several studies have shown the changes in the level of apoptosis-related gene expression in shrimp against viral-infection [44]–[47]. In this study, we found six differentially expressed genes, including beclin, catalase, oncoprotein nm23, cathepsin B, deoxyribonuclease I, prohibitin and caspase, are involved in apoptosis. These genes might respond to the WSSV infection through controlling the extent of the apoptosis in the defense response, but the detailed mechanisms still need further investigation. Mitogen-activated protein kinase (MAPK) signaling is another noteworthy pathway. MAPK is an important process for animal defense against pathogens, and was activated during virus infection and contributed to virus replication in animal cells [48]. Among the differentially expressed genes, we found heat shock protein 70, max protein and P38 mapk shared homology to signaling molecules of the MAPK signaling pathway. These genes were over-expressed in the WSSV-infected shrimp, with the fold changes (log2 ratio) ranging from 1.15 to 4.03, indicating they might play an important role in response to WSSV infection. However, the underlying molecular mechanisms also remain unclear. To our knowledge, this is the first report showing the changes in the level of MAPK-related gene expression in WSSV-infected shrimp. Several differentially expressed genes seem to be involved in toll-like receptor (TLR) signaling pathway. TLRs are a class of proteins that have been shown to play a critical role in the activation of innate immunity by recognizing distinct pathogen-associated molecular patterns (PAMPs) derived from viruses, bacteria and fungi [49]. However, there are few reports of interactions between viruses and the TLR signaling pathway in crustacean. Here we found that polyubiquitin-C and caspase, belonging to TLR signaling pathway, were both up-regulated after WSSV infection. Several other genes seem to belong to other immune-relevant pathways, such as Wnt signaling pathway (including ring box protein, transmembrane receptor, and caprin-1) and antigen processing and presentation pathway (including cathepsin L, cathepsin B, autoantigen, legumain, and gamma-interferon-inducible lysosomal thiol reductase). These pathways have also been reported to play critical role in animal defense against viral pathogens [50], [51].
Conclusion
Understanding the molecular interactions between L. vannamei and WSSV is useful, not only for studies on the molecular responses and defense mechanism in shrimp against pathogens, but also for agricultural practice that aims to develop effective strategies to prevent viral disease. However, these studies have been hampered by the lack of genomic resources of L. vannamei. In this study, we investigated the transcriptome profile of WSSV-infected L. vannamei using the 454 RNA-seq technology. In conclusion, this study allowed us to (1) Obtain over 14538 unigenes from L. vannamei, (2) Compare the transcriptomes between WSSV-infected and uninfected shrimp, and (3) Identify potential pathways and candidate genes involved in L. vannamei response to WSSV infection. Further studies on the functionality of these genes will provide valuable information on the antiviral mechanism in shrimp. Furthermore, the large number of transcripts obtained in this study provides a strong basis for future genomic research on shrimp.
Supporting Information
Table S1.
Primers used in qRT-PCR for validation of differentially expressed genes.
https://doi.org/10.1371/journal.pone.0073218.s001
(XLS)
Table S2.
List of differentially expressed genes with Nr, Nt, UniProtKB/Swiss-Prot, GO, COG and ORF analysis.
https://doi.org/10.1371/journal.pone.0073218.s002
(XLS)
Author Contributions
Conceived and designed the experiments: XHC DGZ. Performed the experiments: XLC DXX YML NM YZZ CLY ML QY ZPL HW. Analyzed the data: DGZ. Wrote the paper: DGZ.
References
- 1. Zhou J, Fang W, Yang X, Zhou S, Hu L, et al. (2012) A nonluminescent and highly virulent Vibrio harveyi strain is associated with “bacterial white tail disease” of Litopenaeus vannamei shrimp. PLoS One 7: e29961.
- 2. Naylor RL, Goldburg RJ, Primavera JH, Kautsky N, Beveridge MC, et al. (2000) Effect of aquaculture on world fish supplies. Nature 405: 1017–1024.
- 3. Valles-Jimenez R, Cruz P, Perez-Enriquez R (2004) Population genetic structure of Pacific white shrimp (Litopenaeus vannamei) from Mexico to Panama: microsatellite DNA variation. Mar Biotechnol (NY) 6: 475–484.
- 4. Lightner DV, Redman RM, Poulos BT, Nunan LM, Mari JL, et al. (1997) Risk of spread of penaeid shrimp viruses in the Americas by the international movement of live and frozen shrimp. Rev Sci Tech 16: 146–160.
- 5. Escobedo-Bonilla CM, Alday-Sanz V, Wille M, Sorgeloos P, Pensaert MB, et al. (2008) A review on the morphology, molecular characterization, morphogenesis and pathogenesis of white spot syndrome virus. J Fish Dis 31: 1–18.
- 6. t Hoen PA, Ariyurek Y, Thygesen HH, Vreugdenhil E, Vossen RH, et al. (2008) Deep sequencing-based expression analysis shows major advances in robustness, resolution and inter-lab portability over five microarray platforms. Nucleic Acids Res 36: e141.
- 7. 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.
- 8. van Hulten MC, Witteveldt J, Peters S, Kloosterboer N, Tarchini R, et al. (2001) The white spot syndrome virus DNA genome sequence. Virology 286: 7–22.
- 9. Wang YG, Lee KL, Najiah M, Shariff M, Hassan MD (2000) A new bacterial white spot syndrome (BWSS) in cultured tiger shrimp Penaeus monodon and its comparison with white spot syndrome (WSS) caused by virus. Dis Aquat Organ 41: 9–18.
- 10. Tan YW, Shi ZL (2011) Genotyping of white spot syndrome virus in Chinese cultured shrimp during 1998-1999. Virol Sin 26: 123–130.
- 11. Robalino J, Almeida JS, McKillen D, Colglazier J, Trent HF 3rd, et al. (2007) Insights into the immune transcriptome of the shrimp Litopenaeus vannamei: tissue-specific expression profiles and transcriptomic responses to immune challenge. Physiol Genomics 29: 44–56.
- 12. Flegel TW, Sritunyalucksana K (2011) Shrimp molecular responses to viral pathogens. Mar Biotechnol (NY) 13: 587–607.
- 13. James R, Thampuran N, Lalitha KV, Rajan LA, Joseph TC (2010) Differential gene expression profile of the hepatopancreas of white spot syndrome virus infected Fenneropenaeus indicus by suppression subtractive hybridization. Fish Shellfish Immunol 29: 884–889.
- 14. Nayak S, Singh SK, Ramaiah N, Sreepada RA (2010) Identification of upregulated immune-related genes in Vibrio harveyi challenged Penaeus monodon postlarvae. Fish Shellfish Immunol 29: 544–549.
- 15. Scherbaum S, Ertas B, Gebauer W, Burmester T (2010) Characterization of hemocyanin from the peacock mantis shrimp Odontodactylus scyllarus (Malacostraca: Hoplocarida). J Comp Physiol B 180: 1235–1245.
- 16. Linzen B (1989) [Blue blood: structure and evolution of hemocyanin]. Naturwissenschaften 76: 206–211.
- 17. Song KK, Li DF, Zhang MC, Yang HJ, Ruan LW, et al. (2010) Cloning and characterization of three novel WSSV recognizing lectins from shrimp Marsupenaeus japonicus. Fish Shellfish Immunol 28: 596–603.
- 18. Liu YC, Li FH, Dong B, Wang B, Luan W, et al. (2007) Molecular cloning, characterization and expression analysis of a putative C-type lectin (Fclectin) gene in Chinese shrimp Fenneropenaeus chinensis. Mol Immunol 44: 598–607.
- 19. Tassanakajon A, Amparyup P, Somboonwiwat K, Supungul P (2011) Cationic antimicrobial peptides in penaeid shrimp. Mar Biotechnol (NY) 13: 639–657.
- 20. Gomez-Anduro GA, Ascencio-Valle F, Peregrino-Uriarte AB, Campa-Cordova A, Yepiz-Plascencia G (2012) Cytosolic manganese superoxide dismutase genes from the white shrimp Litopenaeus vannamei are differentially expressed in response to lipopolysaccharides, white spot virus and during ontogeny. Comp Biochem Physiol B Biochem Mol Biol 162: 120–125.
- 21. Margulies M, Egholm M, Altman WE, Attiya S, Bader JS, et al. (2005) Genome sequencing in microfabricated high-density picolitre reactors. Nature 437: 376–380.
- 22. Liu L, Li Y, Li S, Hu N, He Y, et al. (2012) Comparison of next-generation sequencing systems. J Biomed Biotechnol 2012: 251364.
- 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. Kiatpathomchai W, Boonsaeng V, Tassanakajon A, Wongteerasupaya C, Jitrapakdee S, et al. (2001) A non-stop, single-tube, semi-nested PCR technique for grading the severity of white spot syndrome virus infections in Penaeus monodon. Dis Aquat Organ 47: 235–239.
- 25. Bonnichon V, Lightner DV, Bonami JR (2006) Viral interference between infectious hypodermal and hematopoietic necrosis virus and white spot syndrome virus in Litopenaeus vannamei. Dis Aquat Organ 72: 179–184.
- 26. Zheng Y, Zhao L, Gao J, Fei Z (2011) iAssembler: a package for de novo assembly of Roche-454/Sanger transcriptome sequences. BMC Bioinformatics 12: 453.
- 27. Apweiler R, Martin MJ, O'Donovan C, Magrane M, Alam-Faruque Y, et al. (2013) Update on activities at the Universal Protein Resource (UniProt) in 2013. Nucleic Acids Res 41: D43–D47.
- 28.
Mount DW (2007) Using the Basic Local Alignment Search Tool (BLAST). CSH Protoc 2007 : pdb top17.
- 29. Tatusov RL, Galperin MY, Natale DA, Koonin EV (2000) The COG database: a tool for genome-scale analysis of protein functions and evolution. Nucleic Acids Res 28: 33–36.
- 30. Wixon J, Kell D (2000) The Kyoto encyclopedia of genes and genomes–KEGG. Yeast 17: 48–55.
- 31. Romualdi C, Bortoluzzi S, D'Alessi F, Danieli GA (2003) IDEG6: a web tool for detection of differentially expressed genes in multiple tag sampling experiments. Physiol Genomics 12: 159–162.
- 32. Storey J (2002) A direct approach to false discovery rates. Journal of the Royal Statistical Society 64: 19.
- 33. Cheng W, Liu CH, Tsai CH, Chen JC (2005) Molecular cloning and characterisation of a pattern recognition molecule, lipopolysaccharide- and beta-1,3-glucan binding protein (LGBP) from the white shrimp Litopenaeus vannamei. Fish Shellfish Immunol 18: 297–310.
- 34. Chaozheng Li SW, Yonggui Chen , Xiaoqiang Yu , Ling Lu , Haiqing Zhang , et al. (2012) Analysis of Litopenaeus vannamei Transcriptome Using the Next-Generation DNA Sequencing Technique. 7: e47442.
- 35. You FM, Huo N, Deal KR, Gu YQ, Luo MC, et al. (2011) Annotation-based genome-wide SNP discovery in the large and complex Aegilops tauschii genome using next-generation sequencing without a reference genome sequence. BMC Genomics 12: 59.
- 36. Franssen SU, Shrestha RP, Brautigam A, Bornberg-Bauer E, Weber AP (2011) Comprehensive transcriptome analysis of the highly complex Pisum sativum genome using next generation sequencing. BMC Genomics 12: 227.
- 37. Ye J, Fang L, Zheng H, Zhang Y, Chen J, et al. (2006) WEGO: a web tool for plotting GO annotations. Nucleic Acids Res 34: W293–297.
- 38. Kanehisa M, Goto S (2000) KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res 28: 27–30.
- 39.
Bruce Alberts AJ, Julian Lewis, Martin Raff, Keith Roberts, and Peter Walter (2002) Molecular Biology of the Cell. New York: Garland Science.
- 40. Barakat A, DiLoreto DS, Zhang Y, Smith C, Baier K, et al. (2009) Comparison of the transcriptomes of American chestnut (Castanea dentata) and Chinese chestnut (Castanea mollissima) in response to the chestnut blight infection. BMC Plant Biol 9: 51.
- 41. Cohen JJ (1993) Apoptosis: the physiologic pathway of cell death. Hosp Pract (Off Ed) 28: 35–43.
- 42. Tey SK, Khanna R (2012) Autophagy mediates transporter associated with antigen processing-independent presentation of viral epitopes through MHC class I pathway. Blood 120: 994–1004.
- 43. Tey SK, Khanna R (2012) Host immune system strikes back: autophagy-mediated antigen presentation bypasses viral blockade of the classic MHC class I processing pathway. Autophagy 8: 1839–1841.
- 44. Compeer EB, Flinsenberg TW, van der Grein SG, Boes M (2012) Antigen processing and remodeling of the endosomal pathway: requirements for antigen cross-presentation. Front Immunol 3: 37.
- 45. Purdie AC, Plain KM, Begg DJ, de Silva K, Whittington RJ (2012) Expression of genes associated with the antigen presentation and processing pathway are consistently regulated in early Mycobacterium avium subsp. paratuberculosis infection. Comp Immunol Microbiol Infect Dis 35: 151–162.
- 46. Tiwari N, Garbi N, Reinheckel T, Moldenhauer G, Hammerling GJ, et al. (2007) A transporter associated with antigen-processing independent vacuolar pathway for the MHC class I-mediated presentation of endogenous transmembrane proteins. J Immunol 178: 7932–7942.
- 47. Lee JA, Sinkovits RS, Mock D, Rab EL, Cai J, et al. (2006) Components of the antigen processing and presentation pathway revealed by gene expression microarray analysis following B cell antigen receptor (BCR) stimulation. BMC Bioinformatics 7: 237.
- 48. Andrade AA, Silva PN, Pereira AC, De Sousa LP, Ferreira PC, et al. (2004) The vaccinia virus-stimulated mitogen-activated protein kinase (MAPK) pathway is required for virus multiplication. Biochem J 381: 437–446.
- 49. Dobrovolskaia MA, Medvedev AE, Thomas KE, Cuesta N, Toshchakov V, et al. (2003) Induction of in vitro reprogramming by toll-like receptor (TLR)2 and TLR4 agonists in murine macrophages: Effects of TLR “homotolerance” versus “heterotolerance” on NF-kappa B signaling pathway components. Journal of Immunology 170: 508–519.
- 50. Gan DD, Reiss K, Carrill T, Del Valle L, Croul S, et al. (2001) Involvement of Wnt signaling pathway in murine medulloblastoma induced by human neurotropic JC virus. Oncogene 20: 4864–4870.
- 51. Brodsky FM (1992) Antigen processing and presentation: close encounters in the endocytic pathway. Trends Cell Biol 2: 109–115.