Camilla A. Santos1*
Sónia C. S. Andrade2
Ana K. Teixeira3Flávio Farias3Karin Kurkjian4Ana C. Guerrelhas4João L. Rocha3
Pedro M. Galetti Jr.1
Patrícia D. Freitas1- 1Departamento de Genética e Evolução, Universidade Federal de São Carlos, São Carlos, Brazil
- 2Departamento de Genética e Biologia Evolutiva, Instituto de Biociências, Universidade de São Paulo, São Paulo, Brazil
- 3Genearch Aquacultura, Rio do Fogo, Brazil
- 4Aquatec Larvicultura de Camarão Marinho, Canguaretama, Brazil
Introduction
The Pacific white shrimp Litopenaeus vannamei represents one of the main marine shrimp species commercialized worldwide (Gucic et al., 2013; FAO, 2014; Senapati et al., 2017). Despite its importance on aquaculture and fishery activities, the exposure of wild and captive stocks to pathogens, such as fungi, bacteria and mainly virus, has commonly leading to concerns on healthy and population viability, causing damages related to quality of the shrimp used for local and international trading and human feeding (Wang and Wang, 2013; Liu et al., 2014; Rao et al., 2015). The White Spot Syndrome Virus (WSSV), for instance, may spread very quickly in the open sea, with even more severe harms for the aquaculture stocks, sometimes leading to 100% mortality (FAO, 2017). In addition to the animal health, growth is also an important value for the shrimp aquaculture (Liu et al., 2005; Jung et al., 2013; Lv et al., 2014). This scenario causes large economic losses on the global shrimp farming (Chen et al., 2015) and also on the annual fishing catches (Xue et al., 2013).
Despite the remarkable importance of penaeid shrimp, few data about determinant genes related to growth and immunity are already described for crustacean (Jindra et al., 2013; Jung et al., 2013; Lv et al., 2014), added to the fact that molecular mechanisms involved in growth and immunity are not yet known in detail. In this manner, it becomes evident the need of expanding the knowledge about genes and transcripts involved in cellular responses and physiology of this group. The lack of a shrimp reference genome highlights the importance of studies addressing next generation sequencing (NGS) as a relevant source of information (Guo et al., 2013; Santos et al., 2014; Powell et al., 2015).
In this study, we used the RNA-seq approach for assembling a robust L. vannamei transcriptome obtained from muscle tissue of captive individuals evaluated for growth and high survival performances, and also from hepatopancreas and muscle of animals exposed to White Spot Virus. Sequencing was carried out using Illumina technology (San Diego, California, USA), and the transcripts blastX annotation against Swissprot database returned more than 20 thousand unigenes with known gene products, identifying a higher amount than those previously reported for other marine shrimp transcriptomes (Nguyen et al., 2016; Rao et al., 2016). These findings certainly constitute an important source of genetic information for further studies, considering gene differential expression analysis and single nucleotide polymorphism (SNPs) identification in crustacean species.
Value of the Data
- L. vannamei has enormous commercial value among penaeid shrimp, which account for one-third of the annual global crustacean fisheries.
- The transcriptome reported herein was produced from hepatopancreas and muscle tissue samples of specific pathogen free (SPF) shrimp submitted to selective breeding programs for high rates of survival and fast growth; and to exposure to White Spot Virus.
- The percentage of unigenes presenting SwissProt hits described in this study is higher than those observed for other shrimp Illumina transcriptomes (e.g., Nguyen et al., 2016; Rao et al., 2016).
- These data certainly constitute an important source of genomic information for L. vannamei and other penaeid species in addition to provide resources for the aquaculture studies aiming higher growth gains and WSSV-resistant shrimp.
Materials and Methods
Experimental Design and Sampling
The L. vannamei samples which were used in this work are SPF shrimp submitted to selective breeding programs for high rates of survival and fast growth in two shrimp commercial companies (Genearch and Aquatec), both located in Rio Grande do Norte (RN) state, Brazil. All samples were kept in a sanitary environment monitored for the presence of many pathogens, including WSSV and were collected in the first semester of 2015. The genetic improvement programs have been conducted based on performance of closed familiar lines, according to quantitative analyses and genetic gain parameters evaluated by the companies (data not available). For the sampling of individuals evaluated for growth and survival performance, 20 closed families in triplicates were assessed, and four families with the higher and four with the lower growth rates were selected. We sampled muscle tissue (pleopods) of one animal from each family, totalizing eight animals evaluated for growth performance. The individuals were about 2 months aged and weighed around 1.5 g. For the WSSV exposed shrimp sampling, we used SPF individuals of a hybrid line from Aquatec, which was originated from the crossing of two different SPF lines, one from Aquatec and another from Genearch. These hybrid shrimps were exposed to WSSV in a farm pound with a high incidence of the disease. Before and after the exposure to the virus, shrimp hemolymph samples were submitted to quantitative real-time PCR (qPCR) analyzes in order to confirm the presence or absence of the WSSV, following a procedure proposed by Silva et al. (2011). The efficiency of qPCR amplification was determined by the method described by Pfaffl et al. (2002) and the specificity of the amplification was confirmed by the analysis of the dissociation curve (Ririe et al., 1997). After the contact with the WSSV and qPCR tests, we selected one positive animal showing WSSV symptoms, and a negative one with no symptoms, both weighting about 6 g. Later, hepatopancreas and muscle tissues from each individual were sampled, totalizing four samples from these two shrimps evaluated to WSSV. Hepatopancreas and muscle tissues were chosen given their relevance to immune response in crustacean (Fan et al., 2016), and as a target for pathogen infections, including WSSV (Yu et al., 2017), respectively. The tissue samples were conditioned in RNA later (Thermo Fisher Scientific, Waltham, MA, USA) and kept in biofreezer (−80°C) for the isolation of total RNA.
RNA Isolation and cDNA Library Preparation
In total, 12 cDNA libraries were established. Eight libraries were constructed using muscle tissue of individuals evaluated for growth performance; the four remaining ones were constructed using muscle and also hepatopancreas from the two individuals exposed and evaluated to WSSV by qPCR. Total RNA isolation was performed using the Trizol®/chloroform protocol (Chomczynski and Mackey, 1995), modified for the centrifugation time in ethanol washing (12 min) and added a purification stage using lithium chloride (7.5 M) for elimination of phenol traces. The quality was assessed in a Qμbit fluorometer (Thermo Fisher Scientific), quantity was evaluated in a NanoDrop spectrophotometer (Thermo Fisher Scientific) and samples with ratios between 1.8 and 2.2 were considered pure. The integrity of the samples was confirmed in a BioAnalyser equipment (Agilent Technologies Inc., Santa Clara, CA, USA) and those with RNA Integrity Number (RIN) ≥ 6 were stored at −80°C. A cDNA library was produced from each sample using a TruSeq RNA Library Preparation V2 kit (Illumina Inc., San Diego, California, USA).
Sequencing, de Novo Assembly and BlastX Annotation
Sample quantification was carried out in qPCR with the Universal Library Quantification Kit (KAPA Biosystems, Wilmington, MA, USA), prior to loading into the sequencer. The libraries were grouped and ran on the Illumina HiSeq 2500 platform with 2 × 100 bp paired-end, using a TruSeq SBS V3 kit (Illumina Inc., Thermo Fisher Scientific), at the Laboratório de Biotecnologia Animal, Escola Superior de Agricultura “Luiz de Queiroz”—Universidade de São Paulo (ESALQ-USP). The quality of the raw data generated after sequencing was visualized in the FastQC software (version 0.10.1) (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/). All reads were filtered for Phred quality (QS) 23 (sequence inside) and 30 (sequence edges) and minimum length of 65 bp using the SeqyClean software (v.1.9.9) (https://github.com/ibest/seqyclean). The same package was also used to remove contaminant sequences (primers and vectors) using the Univec database (https://www.ncbi.nlm.nih.gov/tools/vecscreen/univec/). After filtering, the reference de novo assembly was performed in the Trinity software (Grabherr et al., 2011). The sequences were normalized using Trinity's insilico_read_normalization.pl script, with 50X adopted as the minimum coverage value. For the assembly, 300 bp was defined as the minimum size of the assembled contigs. The CD-Hit package (Li and Godzik, 2006) and DustMasker (Morgulis et al., 2006) were used to remove redundant contigs with more than 95% similarity and low complexity sequences, respectively. The TransDecoder package (http://transdecoder.sourceforge.net/) was used to identify the contigs candidate coding regions. Meanwhile, Trinotate pipeline (https://trinotate.github.io/) was employed for annotation of the sequences through the following databases: Uniprot (uniref90 + SwissProt) with cut-off value of 1e10−5, Gene Ontology (GO) (Ashburner et al., 2000) for the GO terms Biological Process, Molecular Function and Cellular Component and KEGG (Kyoto Encyclopedia of Genes and Genomes) (Kanehisa et al., 2012), with the identification of the main metabolic pathways found. The paired reads fastq files from all the samples used in transcriptome analysis (BioProject: Litopenaeus vannamei hepatopancreas and muscle Transcriptome) are available in the NCBI Sequence Read Archive (SRA) under the accession number SRP128934. This Transcriptome Shotgun Assembly project has been deposited at DDBJ/EMBL/GenBank under the accession GGKO00000000. The version described in this paper is the first version, GGKO01000000.
Results and Discussion
The 488,681,310 transcripts from the 12 libraries were filtered on SeqyClean for contaminants, rRNA and low quality sequences. From these, about 30 million paired reads (29,376,967) remained after normalization in Trinity package, resulting in 28,654,475 (97.5%) fragments used in the reference transcriptome after filtering. We obtained 63,105 transcripts with 2,511 bp mean length and N50 of 3,464 bp. From this total, 20,865 (33%) showed blastX hits on SwissProt database (Supplementary Material 1; Table 1). The contigs filtered to <95% of similarity were considered unigenes after the reference assembly. The contig average size found here was greater than those previously reported in L. vannamei (Ghaffari et al., 2014), Penaeus monodon (Nguyen et al., 2016) and Macrobrachium rosenberguii (Rao et al., 2015, 2016) Illumina transcriptomes. The longest contig size observed in our study was equal to 20,061 bp, being superior to the value of 13,578 bp reported for P. monodon transcriptome (Nguyen et al., 2016). After elimination of the redundant contigs, low complexity sequences and isoforms, candidate-coding regions were identified, remaining 14,124 unigenes. The de novo assembly and the automatic annotation results are shown in Table 1.
Table 1. Statistics of the reference assembly after normalization and an overview of the functional annotation of transcripts.
Following the sequences alignment against SwissProt database, the focus was given to genes that showed hits for arthropod species and fitness-related functions, such as growth and disease resistance. Among the most frequent proteins exclusively identified in the animals evaluated for growth performance are those related to myosin (MYSA) (41%), tropomyosin and myosin heavy chain (9%), chitinase (6%) and molt-inhibiting hormone-like (2%) genes (Figure 1). Myosin acts on the formation of the muscular myofibrils and, along with the paramyosin, is responsible for the stabilization of the sarcomere (Liu et al., 2005). There is evidence that increased myosin production and reorganization of myofibrils is associated with growth during and after shrimp molt (Cesar and Yang, 2007; Li et al., 2014). Chitin constituent proteins are among those responsible for the digestion of arthropod chitin exoskeleton, enabling molt and animal growth (Li et al., 2015). Among the most frequent genes in the WSSV-exposed shrimp samples returning blastX hits were mainly those related to crustacyanin 1 (CRA1) and 2 (CRA2) (13%), lectins (8%) and hemocyanin B (HCYB), and C (HCYC) heavy chain (6%) (Figure 1). Hemocyanin and crustacyanin proteins have a wide range of action against viruses, especially WSSV, acting on animal stress and survival responses (Fan et al., 2016), agglutination of the pathogen on hemolymph and cell lysis (Shockey et al., 2009; Zheng et al., 2016). On the other hand, lectins act on the pathogen specific recognition through molecular structures to a given molecule (Song et al., 2010).
Figure 1. SwissProt blastX hits found for the top 10 proteins related to growth and immunity in arthropod species. The protein name and blastX hit in arthropod species are shown. Most of the blastX returned hits for myosin heavy chain (constituent of muscle) (41%) and crustacyanin (13%) (acting against virus in crustacean) proteins.
The de novo transcriptome assembly from muscle and hepatopancreas of non-infected and infected shrimp, enabled the identification of important genes with potential fitness-related functions involved in growth and immunity. Such assembly strategy, gathering two different tissues with diverse functions, was performed in the effort to reach a wider range and representativeness of transcripts, producing the most robust transcriptome dataset for L. vannamei reported to date. It becomes of particular interest for the species lacking a reference genome, such as the majority of crustacean, including L. vannamei. Considering the high accuracy and wide coverage of RNA-seq produced here, this dataset can be useful for a gamma of further approaches, such as differential gene expression studies, analyses of isoforms originated from alternative splicing (Rao et al., 2015; Sun et al., 2015), SNPs calling (Cui et al., 2014; Santos et al., 2014; Yu et al., 2014), with potential to high-density chip and linkage map (Baranski et al., 2014), apart from genome wide association (GWAS) studies (Yu et al., 2017).
Author Contributions
CS has contributed to the transcriptome experiments, data analysis, and article writing and discussion; SA has conducted the bioinformatics analyses and participated in the article writing and discussion; KK has performed the qPCR analyses and the experiments related to the WSSV exposure; JR, FF, AT, and AG have conducted the selective breeding programs in the Genearch and Aquatec companies; PG has collaborated in the article writing and discussion; PF has lead the work, participating in the data analyses, article writing and discussion.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
This research was supported by a grant to PF from the Brazilian agency Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP–2012/17322-8 and 2012/13069-6). PG thanks Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq–308385/2014-4) and CS thanks Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) and FAPESP. The authors thank Professor Luiz L. Coutinho and the Centro de Genômica Funcional do Laboratório de Biotecnologia Animal (ESALQ-USP), Piracicaba-SP, Brazil.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2018.00120/full#supplementary-material
Supplementary Material 1. All the blastX hits found for 20,865 transcripts in the SwissProt database obtained for Litopenaeus vannamei.
Abbreviations
PLs, post-larvae; SPF, Specific Pathogen Free; WSSV, White Spot Syndrome Virus.
References
Ashburner, M., Ball, C. A., Blake, J. A., Botstein, D., Butler, H., Cherry, J. M., et al. (2000). Gene ontology: tool for the identification of biology. Nat. Genet. 25, 25–29. doi: 10.1038/75556
Baranski, M., Gopikrishna, G., Robinson, N. A., Katneni, V. K., Shekhar, M. S., Shanmugakarthik, J., et al. (2014). The development of a high density linkage map for black tiger shrimp (Penaeus monodon) based on cSNPs. PloS ONE 9:e85413. doi: 10.1371/journal.pone.0085413
Cesar, J. R. O., and Yang, J. (2007). Expression patterns of ubiquitin, heat shock protein 70, ɑ-actin and β-actin over the molt cycle in the abdominal muscle of marine shrimp Litopenaeus vannamei. Mol. Reprod. Dev. 74, 554–559. doi: 10.1002/mrd.20605
Chen, K., Li, E., Li, T., Xu, C., Wang, X., Lin, H., et al. (2015). Transcriptome and molecular pathway analysis of the hepatopancreas in the Pacific white shrimp Litopenaeus vannamei under chronic low-salinity stress. PLoS ONE 10:e0131503. doi: 10.1371/journal.pone.0131503
Chomczynski, P., and Mackey, K. (1995). Short technical reports. Modification of the TRI reagent procedure for isolation of RNA from polysaccharide- and proteoglycan-rich sources. BioTechniques 19, 942–945.
Cui, J., Wang, H., Liu, S., Qiu, X., Jang, Z., and Wang, X. (2014). Transcriptome analysis of the gill of Takifugu rubripes using illumina sequencing for discovery of SNPs. Comp. Biochem. Physiol. D Genomics Proteomics 10, 44–51. doi: 10.1016/j.cbd.2014.03.001
Fan, L., Wang, A., Miao, Y., Liao, S., Ye, C., and Lin, Q. (2016). Comparative proteomic identification of the hepatopancreas response to cold stress in white shrimp, Litopenaeus vannamei. Aquaculture 454, 27–34. doi: 10.1016/j.aquaculture.2015.10.016
Food Agriculture Organization of the United Nations. (FAO) (2014). Cultured Aquatic Species Information Programme – Litopenaeus vannamei. Available online at: http://www.fao.org/fishery/culturedspecies/Penaeus_vannamei/en (Accessed September 10, 2017).
Food Agriculture Organization of the United Nations (FAO) (2017). Threats and Risks of Introducing Alien Shrimp Species. Available online at: http://www.fao.org/docrep/009/a0086e/a0086e10.htm (Accessed September 10, 2017).
Ghaffari, N., Sanchez-Flores, A., Doan, R., Garcia-Orozco, K. D., Chen, P. L., Ochoa-Leyva, A., et al. (2014). Novel transcriptome assembly and improved annotation of the whiteleg shrimp (Litopenaeus vannamei), a dominant crustacean in global seafood mariculture. Sci. Rep. 4:7081. doi: 10.1038/srep07081
Grabherr, M. G., Haas, B. J., Yassour, M., Levin, J. Z., Thompson, D. A., Amit, I., et al. (2011). Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat. Biotechnol. 29, 644–652. doi: 10.1038/nbt.1883
Gucic, M., Cortés-Jacinto, E., Civera-Cerecedo, R., Ricque-Marie, D., and Martínez-Córdova, L. R. (2013). Apparent carbohydrate and lipid digestibility of feeds for whiteleg shrimp, Litopenaeus vannamei (Decapoda: Penaeidae), cultivated at different salinities. Rev. Biol. Trop. 61, 1201–1213.
Guo, H., Ye, C.-X., Wang, A.-L., Xian, J.-A., Liao, S.-A., Miao, Y.-T., et al. (2013). Transcriptome analysis of the Pacific white shrimp Litopenaeus vannamei exposed to nitrite by RNA-seq. Fish Shellfish Immunol. 35, 2008–2016. doi: 10.1016/j.fsi.2013.09.019
Jindra, M., Palli, S. R., and Riddiford, L. M. (2013). The juvenile hormone signaling pathway in insect development. Annual Rev. Entomol. 58, 181–204. doi: 10.1146/annurev-ento-120811-153700
Jung, H., Lyons, R. E., Hurwood, D. A., and Mather, P. B. (2013). Genes and growth performance in crustacean species: a review of relevant genomic studies in crustaceans and other taxa. Rev. Aquaculture 5, 77–110. doi: 10.1111/raq.12005
Kanehisa, M., Goto, S., Sato, Y., Furumichi, M., and Tanabe, M. (2012). KEGG for integration and interpretation of large-scale molecular data sets. Nucleic Acids Res. 40, D109–D114. doi: 10.1093/nar/gkr988
Li, E., Wang, S., Li, C., Wang, X., Chen, K., and Chen, L. (2014). Transcriptome sequencing revealed the genes and pathways involved in salinity stress of Chinese mitten crab, Eriocheir sinensis. Physiol. Genomics 46, 177–190. doi: 10.1152/physiolgenomics.00191.2013
Li, W., and Godzik, A. (2006). Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics 22, 1658–1659. doi: 10.1093/bioinformatics/btl158
Li, X., Xu, Z., Zhou, G., Lin, H., Zhou, J., Zeng, Q., et al. (2015). Molecular characterization and expression analysis of five chitinases associated with molting in the Chinese mitten crab, Eriocheir sinensis. Comp. Biochem. Physiol. Part B Biochem. Mol. Biol. 187, 110–120. doi: 10.1016/j.cbpb.2015.05.007
Liu, H., Miller, M. S., Swank, D. M., Kronert, W. A., Maughan, D. W., and Bernstein, S. I. (2005). Paramyosin phosphorylation site disruption affects indirect flight muscle stiffness and power generation in Drosophila melanogaster. Proc. Natl. Acad. Sci. U.S.A. 102, 10522–10527. doi: 10.1073/pnas.0500945102
Liu, Y., Hou, F., He, S., Qian, Z., Wang, X., Mao, A., et al. (2014). Identification, characterization and functional analysis of a serine protease inhibitor (Lvserpin) from the Pacific white shrimp, Litopenaeus vannamei. Dev. Comp. Immunol. 43, 35–46. doi: 10.1016/j.dci.2013.10.012
Lv, J., Liu, P., Gao, B., Wang, Y., Wang, Z., Chen, P., et al. (2014). Transcriptome analysis of the Portunus trituberculatus: De novo assembly, growth-related gene identification and marker discovery. PLoS ONE 9:e94055. doi: 10.1371/journal.pone.0094055
Morgulis, A., Gertz, E. M., Schäffer, A. A., and Agarwala, R. (2006). A fast and symmetric DUST implementation to mask low-complexity DNA sequences. J. Comput. Biol. 13, 1028–1040. doi: 10.1089/cmb.2006.13.1028
Nguyen, C., Nguyen, T. G., Nguyen, L. V., Pham, H. Q., Nguyen, T. H., Pham, H. T., et al. (2016). De novo assembly and transcriptome characterization of major growth-related genes in various tissues of Penaeus monodon. Aquaculture 464, 545–553. doi: 10.1016/j.aquaculture.2016.08.003
Pfaffl, M. W., Horgan, G. W., and Dempfle, L. (2002). Relative expression software tool (REST©) for group-wise comparison and statistical analysis of relative expression results in real-time PCR. Nucleic Acids Res. 30:e36. doi: 10.1093/nar/30.9.e36
Powell, D., Knibb, W., Remilton, C., and Elizur, A. (2015). De-novo transcriptome analysis of the banana shrimp (Fenneropenaeus merguiensis) and identification of genes associated with reproduction and development. Mar. Genomics 22, 71–78. doi: 10.1016/j.margen.2015.04.006
Rao, R., Bhassu, S., Bing, R. Z. Y., Alinejad, T., Hassan, S. S., and Wang, J. (2016). A transcriptome study on Macrobrachium rosenbergii hepatopancreas experimentally challenged with white spot syndrome virus (WSSV). J. Invert. Pathol. 136, 10–22. doi: 10.1016/j.jip.2016.01.002
Rao, R., Bing Zhu, Y., Alinejad, T., Tiruvayipati, S., Lin Thong, K., Wang, J., et al. (2015). RNA-seq analysis of Macrobrachium rosenbergii hepatopancreas in response to Vibrio parahaemolyticus infection. Gut. Pathog. 7:6. doi: 10.1186/s13099-015-0052-6
Ririe, K. M., Rasmussen, R. P., and Wittwer, C. T. (1997). Product differentiation by analysis of DNA melting curves during the polymerase chain reaction. Anal. Biochem. 245, 154–160. doi: 10.1006/abio.1996.9916
Santos, C. A., Blanck, D. V., and de Freitas, P. D. (2014). RNA-seq as a powerful tool for penaeid shrimp genetic progress. Front. Genet. 5:298. doi: 10.3389/fgene.2014.00298
Senapati, S. R., Kumar, G. P., Singh, C. B., Xavier, K. M., Chouksey, M. K., Nayak, B. B., et al. (2017). Melanosis and quality attributes of chill stored farm raised whiteleg shrimp (Litopenaeus vannamei). J. Appl. Nat. Sci. 9, 626–631.
Shockey, J. E., O'Leary, N. A., de la Vega, E., Browdy, C. L., Baatz, J. E., and Gross, P. S. (2009). The role of crustins in Litopenaeus vannamei in response to infection with shrimp pathogens: an in vivo approach. Dev. Comp. Immunol. 33, 668–673. doi: 10.1016/j.dci.2008.11.010
Silva, S. M. B. C., Pinheiro, A. C. A. S., and Coimbra, M. R. M. (2011). Quantitation of infectious myonecrosis virus indifferent tissues of naturally infected Pacific white shrimp, Litopenaeus vannamei, using real-time PCR with SYBR Green chemistry. J. Virol. Methods 177, 197–201. doi: 10.1016/j.jviromet.2011.08.001
Song, K. K., Li, D. F., Zhang, M. C., Yang, H. J., Ruan, L. W., and Xu, X. (2010). Cloning and characterization of three novel WSSV recognizing lectins from shrimp Marsupenaeus japonicus. Fish Shellfish Immunol. 28, 596–603. doi: 10.1016/j.fsi.2009.12.015
Sun, S., Xuan, F., Fu, H., Zhu, J., Ge, X., and Gu, Z. (2015). Transciptomic and histological analysis of hepatopancreas, muscle and gill tissues of oriental river prawn (Macrobrachium nipponense) in response to chronic hypoxia. BMC Genomics 16:491. doi: 10.1186/s12864-015-1701-3
Wang, X. W., and Wang, J. X. (2013). Pattern recognition receptors acting in innate immune system of shrimp against pathogen infections. Fish Shellfish Immunol. 34, 981–989. doi: 10.1016/j.fsi.2012.08.008
Xue, S., Liu, Y., Zhang, Y., Sun, Y., Geng, X., and Sun, J. (2013). Sequencing and de novo analysis of the hemocytes transcriptome in Litopenaeus vannamei response to white spot syndrome virus infection. PLoS ONE 8:e76718. doi: 10.1371/journal.pone.0076718
Yu, Y., Liu, J., Li, F., Zhang, X., Zhang, C., and Xiang, J. (2017). Gene set based association analyses for the WSSV resistance of Pacific white shrimp Litopenaeus vannamei. Sci. Rep. 7:40549. doi: 10.1038/srep40549
Yu, Y., Wei, J., Zhang, X., Liu, J., Liu, C., Li, F., et al. (2014). SNP discovery in the transcriptome of white Pacific shrimp Litopenaeus vannamei by next generation sequencing. PLoS ONE 9:e87218. doi: 10.1371/journal.pone.0087218
Keywords: RNA-seq, growth, immune genes, WSSV, Litopenaeus vannamei
Citation: Santos CA, Andrade SCS, Teixeira AK, Farias F, Kurkjian K, Guerrelhas AC, Rocha JL, Galetti PM Jr and Freitas PD (2018) Litopenaeus vannamei Transcriptome Profile of Populations Evaluated for Growth Performance and Exposed to White Spot Syndrome Virus (WSSV). Front. Genet. 9:120. doi: 10.3389/fgene.2018.00120
Received: 16 January 2018; Accepted: 26 March 2018;
Published: 10 April 2018.
Edited by:
Youri I. Pavlov, University of Nebraska Medical Center, United StatesReviewed by:
Pierre Boudry, Institut Français de Recherche pour l'Exploitation de la Mer (IFREMER), FranceChu-Fang Lo, Center for Shrimp Disease Control and Genetic Improvement, Taiwan
Copyright © 2018 Santos, Andrade, Teixeira, Farias, Kurkjian, Guerrelhas, Rocha, Galetti and Freitas. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Camilla A. Santos, Y2FtaWxsYS5hbHZlc0B5YWhvby5jb20uYnI=; cGF0ZGZAdWZzY2FyLmJy