RNA-Seq analysis reveals genes associated with resistance to Taura syndrome virus (TSV) in the Pacific white shrimp Litopenaeus vannamei

Temporal Transcriptomic Profiling Reveals Dynamic Changes in Gene Expression of Giant Freshwater Prawn upon Acute Saline-Alkaline Stresses

Bicarbonate and sulfate are among two primary ion constituents of saline-alkaline water, with excessive levels potentially causing metabolic disorders in crustaceans, affecting their molting and interrupting development. As an economically important crustacean species, the molecular adaptive mechanism of giant freshwater prawn Macrobrachium rosenbergii in response to the stress of bicarbonate and sulfate remains unexplored. To investigate the mechanism underlying NaHCO3, Na2SO4, and mixed NaHCO3, Na2SO4 stresses, M. rosenbergii larvae were exposed to the above three stress conditions, followed by total RNA extraction and high-throughput sequencing at eight distinct time points (0, 4, 8, 12, 24, 48, 72, and 96 h). Subsequent analysis revealed 13, 16, and 13 consistently identified differentially expressed genes (DEGs) across eight time points under three stress conditions. These consistently identified DEGs were significantly involved in the Gene Ontology (GO) terms of chitin-based cuticle development, protein-carbohydrate complex, structural constituent of cuticle, carnitine biosynthetic process, extracellular matrix, and polysaccharide catabolic process, indicating that alkaline stresses might potentially impact the energy metabolism, growth, and molting of M. rosenbergii larvae. Particularly, the transcriptome data revealed that DEGs associated with energy metabolism, immunity, and amino acid metabolism were enriched across multiple time points under three stress conditions. These DEGs are linked to Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, including glycolysis/glucogenesis, amino sugar and nucleotide sugar metabolism, and lysine degradation. Consistent enrichment findings across the three stress conditions support conclusions above. Together, these insights are instrumental in enhancing our understanding of the molecular mechanisms underlying the alkaline response in M. rosenbergii larvae. Additionally, they offer valuable perspectives on the regulatory mechanisms of freshwater crustaceans amid saline-alkaline water development.

MicroRNA sequencing analysis reveals immune responses in hepatopancreas of Fenneropenaeus penicillatus under white spot syndrome virus infection

White Spot Disease is one of the most harmful diseases of the red tail shrimp, which can cause devastating economic losses due to the highest mortality up to 100% within a few days. MicroRNAs (miRNAs) are large class of small noncoding RNAs with the ability to post-transcriptionally repress the translation of target mRNAs. MiRNAs are considered to have a significant role in the innate immune response of crustaceans, particularly in relation to antiviral defense mechanisms. Numerous crustacean miRNAs have been verified to be required in host immune defense against viral infection, however, till present, the miRNAs functions of F. penicillatus defense WSSV infection have not been studied yet. Here in this study, for the first time, miRNAs involved in the F. penicillatus immune defense against WSSV infection were identified using high-throughput sequencing platform. A total of 432 miRNAs were obtained including 402 conserved miRNAs and 30 novel predicted miRNAs. Comparative analysis between the WSSV-challenged group and the control group revealed differential expression of 159 microRNAs in response to WSSV infection. Among these, 48 were up-regulated and 111 were down-regulated. Ten candidate MicroRNAs associated with immune activities were randomly selected for qRT-PCR analysis, which confirming the expression profiling observed in the MicroRNA sequencing data. As a result, most differentially expressed miRNAs were down-regulated lead to increase the expression of various target genes that mediated immune reaction defense WSSV infection, including genes related to signal transduction, Complement and coagulation cascade, Phagocytosis, and Apoptosis. Furthermore, the genes expression of the key members in Toll and Imd signaling pathways and apoptotic signaling were mediated by microRNAs to activate host immune responses including apoptosis against WSSV infection. These results will help to understand molecular defense mechanism against WSSV infection in F. penicillatus and to develop an effective WSSV defensive strategy in shrimp farming.

Applying genetic technologies to combat infectious diseases in aquaculture

Disease and parasitism cause major welfare, environmental and economic concerns for global aquaculture. In this review, we examine the status and potential of technologies that exploit genetic variation in host resistance to tackle this problem. We argue that there is an urgent need to improve understanding of the genetic mechanisms involved, leading to the development of tools that can be applied to boost host resistance and reduce the disease burden. We draw on two pressing global disease problems as case studies-sea lice infestations in salmonids and white spot syndrome in shrimp. We review how the latest genetic technologies can be capitalised upon to determine the mechanisms underlying inter- and intra-species variation in pathogen/parasite resistance, and how the derived knowledge could be applied to boost disease resistance using selective breeding, gene editing and/or with targeted feed treatments and vaccines. Gene editing brings novel opportunities, but also implementation and dissemination challenges, and necessitates new protocols to integrate the technology into aquaculture breeding programmes. There is also an ongoing need to minimise risks of disease agents evolving to overcome genetic improvements to host resistance, and insights from epidemiological and evolutionary models of pathogen infestation in wild and cultured host populations are explored. Ethical issues around the different approaches for achieving genetic resistance are discussed. Application of genetic technologies and approaches has potential to improve fundamental knowledge of mechanisms affecting genetic resistance and provide effective pathways for implementation that could lead to more resistant aquaculture stocks, transforming global aquaculture.

New insights into the regulation mechanism of red claw crayfish (Cherax quadricarinatus) hepatopancreas under air exposure using transcriptome analysis

Red claw crayfish (Cherax quadricarinatus) is an important freshwater shrimp species worldwide with enormous economic value. Waterless transportation is an inherent feature of red claw crayfish transportation. However, the high mortality of red claw crayfish is a severe problem in the aquaculture of crayfish after waterless transportation. In this study, we investigated the responses of the hepatopancreas from the red claw crayfish undergoing air exposure stress and normal conditions on transcriptome levels. We used Illumina-based RNA sequencing (RNA-Seq) to perform a transcriptome analysis from the hepatopancreas of red claw crayfish challenged by air exposure. An average of 57,148,800 clean reads per library was obtained, and 33,567 unigenes could be predicted and classified according to their homology with matches in the National Center for Biotechnology Information (NCBI) non-redundant protein sequences (Nr), Gene Ontology (GO), a manually annotated and reviewed protein sequence database (Swiss-Prot), protein families (Pfam), Clusters of Orthologous Groups (COG) of proteins, and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. 690 and 3407 differentially expressed genes (DEGs) were identified between the two stress stages of the red claw crayfish. More DEGs were identified in 12 h, indicating that gene expressions were largely changed at 12 h. Some immune-related pathways and genes were identified according to KEGG and GO enrichment analysis. A total of 12 DEGs involved in immune response and trehalose mechanism were verified by quantitative real-time-polymerase chain reaction (qRT-PCR). The results indicated that the red claw crayfish might counteract the stress of air exposure at the transcriptomic level by increasing expression levels of antioxidant-, immune-, and trehalose metabolism-related genes. These transcriptome results from the hepatopancreas provide significant insights into the influence mechanism of air exposure to the trehalose mechanism and immune response in the red claw crayfish.

RNA-seq analysis revealing the immune response of Neocaridina denticulata sinensis gill to Vibrio parahaemolyticus infection

Vibrio parahaemolyticus causes serious economic losses to the shrimp farming industry. There is still a lack of adequate understanding of the changes in the overall response of N. denticulata sinensis caused by V. parahaemolyticus, particularly with respect to gill tissue, which is severely damaged by the pathogen. In this study, a total of 1358 differentially expressed genes were identified between the PBS control and Vibrio stimulation groups using transcriptome sequencing techniques. After further screening and analysis, many immune-related genes were obtained, involving lysosome pathway, metabolic process, chitin-binding protein, and serine protease family members. In addition, we randomly selected six DEGs in the lysosome pathway for qRT-PCR verification, and the results showed that their expression patterns were consistent with the RNA-seq. The results demonstrate the molecular regulation of the gill tissue response to V. parahaemolyticus infection in N. denticulata sinensis, contributing to the understand of the complex and efficient innate immune system and defense mechanisms in crustaceans.

The Molecular Mechanism of Hemocyte Immune Response in Marsupenaeus japonicus Infected With Decapod Iridescent Virus 1

As a new type of shrimp lethal virus, decapod iridescent virus 1 (DIV1) has caused huge economic losses to shrimp farmers in China. Up to now, DIV1 has been detected in a variety of shrimps, but there is no report in Marsupenaeus japonicus. In the current study, we calculated the LC₅₀ to evaluate the toxicity of DIV1 to M. japonicus and determined through nested PCR that M. japonicus can be the host of DIV1. Through enzyme activity study, it was found that DIV1 can inhibit the activities of superoxide dismutase, catalase, lysozyme, and phenoloxidase, which could be a way for DIV1 to achieve immune evasion. In a comprehensive study on the transcriptomic changes of M. japonicus in response to DIV1 infection, a total of 52,287 unigenes were de novo assembled, and 20,342 SSR markers associated with these unigenes were obtained. Through a comparative transcriptomic analysis, 6,900 differentially expressed genes were identified, including 3,882 upregulated genes and 3,018 downregulated genes. The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis showed that some GO terms related to virus invasion, replication, and host antiviral infection were promoted under DIV1 infection, such as carbohydrate binding, chitin binding, chitin metabolic process, and DNA replication initiation, and some KEGG pathways related to immune response were significantly influenced by DIV1 infection, including Toll and IMD signaling pathway, JAK-STAT signaling pathway, IL-17 signaling pathway, C-type lectin receptor signaling pathway, complement and coagulation cascades, antigen processing and presentation, necroptosis, apoptosis, NOD-like receptor signaling pathway, apoptosis-multiple species, and TNF signaling pathway. Further analysis showed that STAT, Dorsal, Relish, heat shock protein 70 (HSP70), C-type lectins, and caspase play an important role in DIV1 infection. This is the first detailed study of DIV1 infection in M. japonicus, which initially reveals the molecular mechanism of DIV1 infection in M. japonicus by using the transcriptome analysis of hemocytes combined with enzyme activity study.

A novel Dicistro-like virus discovered in Procambarus clarkii with "Black May" disease

The peak period of morbidity and death in cultured Procambarus clarkii is around May each year and is called the "Black May" disease. The pathogen causing "Black May" disease is believed to be a white spot syndrome virus (WSSV). In 2018, a significant number of P. clarkii died in the pond culture of Xinglong Township, Xuyi County. Two sampling tests on the affected pond showed that, in addition to WSSV, a novel Dicistro-like virus (PcDV) was present. Genomic sequence analysis indicated that this new virus belongs to the Dicistroviridae family, Picornaviridaes order. A high number of spherical particles were detected in gill tissues of P. clarkii with "Black May" disease by electron microscopy, a finding consistent with viruses from the Picornaviridaes order. From October 2018 to September 2019, we took monthly samples from Hubei, Jiangsu and Anhui provinces, and tested for the presence of PcDV and WSSV in P. clarkii. The detection rates of PcDV in P. clarkii peaked from April to June, consistent with the onset of the "Black May" disease. In conclusion, we believe that the discovery of PcDV will provide new research directions for investigating the pathogens causing "Black May" disease in P. clarkii.

Analysis of transcriptome difference between rapid-growing and slow-growing in Penaeus vannamei

Penaeus vannamei is the principle cultured shrimp species in China. However, with the increase of culture density, the growth difference between individuals is also expanding. Here, we make use of RNA-seq to study the growth mechanisms of P. vannamei. After 120 days, we examined the transcriptomes of rapid-growing individuals (RG) and slow-growing individuals (SG). A total of 2116 and 176 differentially expressed genes (DEGs) were found in SG and RG, respectively. Moreover, the main DEGs are opsin, heat shock protein (HSP), actin, myosin, superoxide dismutase (SOD), cuticle protein, and chitinase. GO analysis further revealed that the DEGs were enriched in biological processes significantly, such as "sensory perception," "sensory perception of light stimulus," "response to stimulus," and "response to stress." Additionally, KEGG enrichment analysis showed that the DEGs were mainly enriched in "pentose and glucuronate interconversions," "amino sugar and nucleotide sugar metabolism," "glycophospholipid biosynthesis," and "glutathione metabolism." Interestingly, the upstream genes in the ecdysone signaling pathway, including molting inhibition hormone (MIH) and crustacean hyperglycemic hormone (CHH), did not differ significantly between RG and SG, which suggests that the cause for the inconsistent growth performance is due to the stress levels rather than the ecdysone signal pathway. In summary, this work provides data that will be useful for future studies on shrimp growth and development.

Comparative transcriptome reveals the response of oriental river prawn (Macrobrachium nipponense) to sulfide toxicity at molecular level

Aquatic environmental pollutants have various impacts on aquaculture. Specifically, sulfide has been established as being toxic to aquatic animals including the oriental river prawn Macrobrachium nipponense. In response, the hepatopancreas has been broadly studied, as it plays a pivotal role in arthropod nutrient digestion and absorption, energy supply, and organ development as well as in crustacean immunity. However, the underlying molecular mechanisms of hepatopancreas's response to sulfide toxicity are still poorly understand. Herein, we used Nova-seq 6000 platform to conduct a comparative transcriptome analysis of gene expression profiles in the hepatopancreas of M. nipponense, while it was under the influence of a semi-lethal sulfide concentration (3.20 mg/L at 48 h). A total of 139 million raw reads were obtained, in which 67,602 transcripts were clustered into 37,041 unigenes for further analysis. After constant sulfide exposure for 48 h, 235 differentially expressed genes, i.e., DEGs (151 up-regulated and 84 down-regulated) were identified in the sulfide treatment group (TGHP) compared with the control group (CGHP). We used GO and KEGG databases to annotate all the DEGs into 1180 functions and 123 pathways, respectively. The metabolic pathways included proximal tubule bicarbonate reclamation, sulfur metabolism, glycolysis and gluconeogenesis, and the TCA cycle; while immune-related pathways contained Ras, Rap1, focal adhesion and platelet activation. Additionally, apoptosis-involved pathways e.g., lysosome, also exhibited remarkable alteration in the presence of sulfide stress. Notably, responses to external stimuli and detoxification genes- such as GSKIP, CRT2, APOD, TRET1, CYP4C3 and HR39- were significantly altered by the sulfide stress, indicating that significant toxicity was transferred through energy metabolism, growth, osmoregulatory processes and immunity. Finally, we demonstrated that in the present of sulfide stress, M. nipponense altered the expression of detoxification- and extracellular stimulation-related genes, and displayed positive resistance via tight junction activation and lysosome pathways. The results of these novel experiments shed light on the hepatopancreas's molecular response to sulfide stress resistance and the corresponding adaptation mechanism; and enable us to identify several potential biomarkers for further studies.

A proteomic approach reveals possible molecular mechanisms and roles for endosymbiotic bacteria in begomovirus transmission by whiteflies

BACKGROUND

Many plant viruses are vector-borne and depend on arthropods for transmission between host plants. Begomoviruses, the largest, most damaging and emerging group of plant viruses, infect hundreds of plant species, and new virus species of the group are discovered each year. Begomoviruses are transmitted by members of the whitefly Bemisia tabaci species complex in a persistent-circulative manner. Tomato yellow leaf curl virus (TYLCV) is one of the most devastating begomoviruses worldwide and causes major losses in tomato crops, as well as in many agriculturally important plant species. Different B. tabaci populations vary in their virus transmission abilities; however, the causes for these variations are attributed among others to genetic differences among vector populations, as well as to differences in the bacterial symbionts housed within B. tabaci.

RESULTS

Here, we performed discovery proteomic analyses in 9 whitefly populations from both Middle East Asia Minor I (MEAM1, formerly known as B biotype) and Mediterranean (MED, formerly known as Q biotype) species. We analysed our proteomic results on the basis of the different TYLCV transmission abilities of the various populations included in the study. The results provide the first comprehensive list of candidate insect and bacterial symbiont (mainly Rickettsia) proteins associated with virus transmission.

CONCLUSIONS

Our data demonstrate that the proteomic signatures of better vector populations differ considerably when compared with less efficient vector populations in the 2 whitefly species tested in this study. While MEAM1 efficient vector populations have a more lenient immune system, the Q efficient vector populations have higher abundance of proteins possibly implicated in virus passage through cells. Both species show a strong link of the facultative symbiont Rickettsia to virus transmission.

Shedding the Light on Litopenaeus vannamei Differential Muscle and Hepatopancreas Immune Responses in White Spot Syndrome Virus (WSSV) Exposure

White Spot Syndrome Virus (WSSV) is one of the main threats to farming Litopenaeus vannamei, the most important crustacean commercialized in aquaculture worldwide. Here, we performed RNA-seq analyses in hepatopancreas and muscle from WSSV-negative (healthy) and WSSV-positive (unhealthy) L. vannamei, previously exposed to the virus, to obtain new insights about the molecular basis of resistance to WSSV. We detected 71% of our reads mapped against the recently described L. vannamei genome. This is the first report mapping RNA-seq transcripts from shrimps exposed to WSSV against the species reference genome. Differentially expressed gene (DEG) analyses were performed for four independent comparisons, and 13,338 DEGs were identified. When the redundancies and isoforms were disregarded, we observed 8351 and 6514 DEGs, respectively. Interestingly, after crossing the data, we detected a common set of DEGs for hepatopancreas and healthy shrimps, as well as another one for muscle and unhealthy shrimps. Our findings indicate that genes related to apoptosis, melanization, and the Imd pathway are likely to be involved in response to WSSV, offering knowledge about WSSV defense in shrimps exposed to the virus but not infected. These data present potential to be applied in further genetic studies in penaeids and other farmed shrimp species.

Transcriptomic profiles reveal that inactivated iridovirus and rhabdovirus bivalent vaccine elicits robust adaptive immune responses against lethal challenge in marbled sleepy goby

Oxyeleotris marmoratus iridovirus (OMIV) and Oxyeleotris marmoratus rhabdovirus (OMRV) are the two major causative agents of disease leading to massive mortality and severe economic losses in marbled sleepy goby (Oxyeleotris marmoratus) industry. It's urgent to develop an effective vaccine against these fatal diseases. In this study, we developed bivalent inactivated vaccine against OMIV and OMRV and evaluated its protective effect in Oxyeleotris marmoratus. The intraperitoneally vaccinated fish were protected against challenge with OMIV and OMRV with both relative percent survival (RPS) of 100%. In addition, deep RNA sequencing was used to analyze the transcriptomic profiles of the spleen tissues at progressive time points post-vaccination with bivalent inactivated vaccine and challenge with OMIV and OMRV infection. Results showed that adaptive immune response was induced in Oxyeleotris marmoratus injected with bivalent inactivated vaccine. Furthermore, robust adaptive immune responses were also detected in vaccinated fish at 7 d and 2 d post-challenge with OMIV and OMRV. Taken together, these results indicated that bivalent inactivated vaccine activated adaptive immune responses in Oxyeleotris marmoratus, and provided protection against OMIV and OMRV lethal challenge.

Analysis of the transcriptome data in Litopenaeus vannamei reveals the immune basis and predicts the hub regulation-genes in response to high-pH stress

Soil salinization erodes the farmlands and poses a serious threat to human life, reuse of the saline-alkali lands as cultivated resources becomes increasingly prominent. Pacific white shrimp (Litopenaeus vannamei) is an important farmed aquatic species for the development and utilization of the saline-alkali areas. However, little is known about the adaptation mechanism of this species in terms of high-pH stress. In the present study, a transcriptome analysis on the gill tissues of L. vannamei in response to high-pH stress (pH 9.3 ± 0.1) was conducted. After analyzing, the cyclic nucleotide gated channel-Ca2+ (CNGC-Ca2+) and patched 1 (Ptc1) were detected as the majority annotated components in the cAMP signaling pathway (KO04024), indicating that the CNGC-Ca2+ and Ptc1 might be the candidate components for transducing and maintaining the high-pH stress signals, respectively. The immunoglobulin superfamily (IgSF), heat shock protein (HSP), glutathione s-transferase (GST), prophenoloxidase/phenoloxidase (proPO/PO), superoxide dismutase (SOD), anti-lipopolysaccharide factor (ALF) and lipoprotein were discovered as the major transcribed immune factors in response to high-pH stress. To further detect hub regulation-genes, protein-protein interaction (PPI) networks were constructed; the genes/proteins "Polymerase (RNA) II (DNA directed) polypeptide A" (POLR2A), "Histone acetyltransferase p300" (EP300) and "Heat shock 70kDa protein 8" (HSPA8) were suggested as the top three hub regulation-genes in response to acute high-pH stress; the genes/proteins "Heat shock 70kDa protein 4" (HSPA4), "FBJ murine osteosarcoma viral oncogene homolog" (FOS) and "Nucleoporin 54kDa" (NUP54) were proposed as the top three hub regulation-genes involved in adapting endurance high-pH stress; the protein-interactions of "EP300-HSPA8" and "HSPA4-NUP54" were detected as the most important biological interactions in response to the high-pH stress; and the HSP70 family genes might play essential roles in the adaptation of the high-pH stress environment in L. vannamei. These findings provide the first insight into the molecular and immune basis of L. vannamei in terms of high-pH environments, and the construction of a PPI network might improve our understanding in revealing the hub regulation-genes in response to abiotic stress in shrimp species and might be beneficial for further studies.

De novo assembly, characterization, functional annotation and expression patterns of the black tiger shrimp (Penaeus monodon) transcriptome

The black tiger shrimp (Penaeus monodon) remains the second most widely cultured shrimp species globally; however, issues with disease and domestication have seen production levels stagnate over the past two decades. To help identify innovative solutions needed to resolve bottlenecks hampering the culture of this species, it is important to generate genetic and genomic resources. Towards this aim, we have produced the most complete publicly available P. monodon transcriptome database to date based on nine adult tissues and eight early life-history stages (BUSCO - Complete: 98.2% [Duplicated: 51.3%], Fragmented: 0.8%, Missing: 1.0%). The assembly resulted in 236,388 contigs, which were then further segregated into 99,203 adult tissue specific and 58,678 early life-history stage specific clusters. While annotation rates were low (approximately 30%), as is typical for a non-model organisms, annotated transcript clusters were successfully mapped to several hundred functional KEGG pathways. Transcripts were clustered into groups within tissues and early life-history stages, providing initial evidence for their roles in specific tissue functions, or developmental transitions. We expect the transcriptome to provide an essential resource to investigate the molecular basis of commercially relevant-significant traits in P. monodon and other shrimp species.

The State of "Omics" Research for Farmed Penaeids: Advances in Research and Impediments to Industry Utilization

Elucidating the underlying genetic drivers of production traits in agricultural and aquaculture species is critical to efforts to maximize farming efficiency. "Omics" based methods (i.e., transcriptomics, genomics, proteomics, and metabolomics) are increasingly being applied to gain unprecedented insight into the biology of many aquaculture species. While the culture of penaeid shrimp has increased markedly, the industry continues to be impeded in many regards by disease, reproductive dysfunction, and a poor understanding of production traits. Extensive effort has been, and continues to be, applied to develop critical genomic resources for many commercially important penaeids. However, the industry application of these genomic resources, and the translation of the knowledge derived from "omics" studies has not yet been completely realized. Integration between the multiple "omics" resources now available (i.e., genome assemblies, transcriptomes, linkage maps, optical maps, and proteomes) will prove critical to unlocking the full utility of these otherwise independently developed and isolated resources. Furthermore, emerging "omics" based techniques are now available to address longstanding issues with completing keystone genome assemblies (e.g., through long-read sequencing), and can provide cost-effective industrial scale genotyping tools (e.g., through low density SNP chips and genotype-by-sequencing) to undertake advanced selective breeding programs (i.e., genomic selection) and powerful genome-wide association studies. In particular, this review highlights the status, utility and suggested path forward for continued development, and improved use of "omics" resources in penaeid aquaculture.

Transcriptome analysis of 'Huanghai No. 2' Fenneropenaeus chinensis response to WSSV using RNA-seq

White spot syndrome (WSS) is one of the most damaging phenomena in the culturing of shrimp. To characterize the mechanisms of the molecular responses to WSSV infection in 'Huanghai No. 2'' Fenneropenaeus chinensis, we used next-generation sequencing to observe the transcriptome after oral infection. A total of 108.6 million clean reads were obtained and assembled into 64,103 final unigenes with an average length of 845 bp (N50 = 1534 bp). The assembled unigenes contained 14,263 significant unigenes after BLASTX against the Nr database (E-value cut-off of 10^-5). After comparison of digital gene expression data between challenged and control shrimp, a total of 896 DEGs after WSSV infection were identified. Gene pathway analysis indicated that 92, 131 and 142 metabolic pathways were affected at early, peak and late phases respectively. Some pathways were related to the immune response, such as the phagosome, complement and coagulation cascades, the antigen processing and presentation pathway and so on. Many immune-related genes were also identified after pathway analysis. Interestingly, some growth-related genes, such as cathepsin L, myosin regulatory light chain 2 smooth muscle, and alpha-amylase were also differentially expressed after WSSV infection, and the correlation between growth trait and WSSV-resistance trait need further research. The expression patterns of eight DEGs were confirmed by quantitative real-time reverse transcription polymerase chain reaction, and there was good agreement between RNA-seq and qRT-PCR. These data will provide valuable information for characterizing the immune mechanism of the response of shrimp's to WSSV.

Transcriptomic analysis reveals differentially expressed genes and a unique apoptosis pathway in channel catfish ovary cells after infection with the channel catfish virus

The channel catfish virus (CCV) can cause lethal hemorrhagic infection in juvenile channel catfish, thereby resulting in a huge economic loss to the fish industry. The genome of the CCV has been fully sequenced, and its prevalence is well documented. However, less is known about the molecular mechanisms and pathogenesis of the CCV. Herein, the channel catfish ovary cells (CCO) were infected with CCV and their transcriptomic sketches were analyzed using an RNA sequencing technique. In total, 72,686,438 clean reads were obtained from 73,231,128 sequence reads, which were further grouped into 747,168 contigs. These contigs were assembled into 49,119 unigenes, of which 20,912 and 18,333 unigenes were found in Nr and SwissProt databases and matched 15,911 and 14,625 distinctive proteins, respectively. From these, 3641 differentially expressed genes (DEGs), comprising 260 up-regulated and 3381 down-regulated genes, were found compared with the control (non-infected) cells. For verification, 16 DEGs were analyzed using qRT-PCR. The analysis of the DEGs and their related cellular signaling pathways revealed a substantial number of DEGs that were involved in the apoptosis pathway induced by CCV infection. The apoptosis pathways were further elucidated using standard apoptosis assays. The results showed that CCV could induce extrinsic apoptosis pathway (instead of a mitochondrial intrinsic apoptosis pathway) in CCO cells. This study helps our understanding of the pathogenesis of CCV and contributes to the prevention of CCV infection in channel catfish.

Novel Insights into Antiviral Gene Regulation of Red Swamp Crayfish, Procambarus clarkii, Infected with White Spot Syndrome Virus

White spot syndrome virus (WSSV), one of the major pathogens of Procambarus clarkii, has caused severe disruption to the aquaculture industry of P. clarkii in China. To reveal the gene regulatory mechanisms underlying WSSV infection, a comparative transcriptome analysis was performed among WSSV-infected susceptible individuals (GS), viral resistant individuals (GR), and a non-infected control group (GC). A total of 61,349 unigenes were assembled from nine libraries. Subsequently, 515 and 1033 unigenes exhibited significant differential expression in sensitive and resistant crayfish individuals compared to the control group (GC). Many differentially expressed genes (e.g., C-type lectin 4, Peroxinectin, Prophenoloxidase, and Serine/threonine-protein kinase) observed in GR and GS play critical roles in pathogen recognition and viral defense reactions after WSSV infection. Importantly, the glycosaminoglycan biosynthesis-chondroitin sulfate/dermatan sulfate pathway was identified to play critical roles in defense to WSSV infection for resistant crayfish individuals by upregulating the chondroitin sulfate related genes for the synthesis of WSSV-sensitive, functional chondroitin sulfate chains containing E units. Numerous genes and the key pathways identified between resistant and susceptible P. clarkii individuals provide valuable insights regarding antiviral response mechanisms of decapoda species and may help to improve the selective breeding of P. clarkii WSSV-resistance.

Transcriptome analysis of Kuruma shrimp (Marsupenaeus japonicus) hepatopancreas in response to white spot syndrome virus (WSSV) under experimental infection

Kuruma shrimp (Marsupenaeus japonicus) is one of the most valuable crustacean species in capture fisheries and mariculture in the Indo-West Pacific. White spot syndrome virus (WSSV) is a highly virulent pathogen which has seriously threatened Kuruma shrimp aquaculture sector. However, little information is available in relation to underlying mechanisms of host-virus interaction in Kuruma shrimp. In this study, we performed a transcriptome analysis from the hepatopancreas of Kuruma shrimp challenged by WSSV, using Illumina-based RNA-Seq. A total of 39,084,942 pair end (PE) reads, including 19,566,190 reads from WSSV-infected group and 19,518,752 reads from non-infected (control) group, were obtained and assembled into 33,215 unigenes with an average length of 503.7 bp and N50 of 601 bp. Approximately 17,000 unigenes were predicted and classified based on homology search, gene ontology, clusters of orthologous groups of proteins, and biological pathway mapping. Differentially expressed genes (DEGs), including 2150 up-regulated and 1931 down-regulated, were found. Among those, 805 DEGs were identified and categorized into 14 groups based on their possible functions. Many genes associated with JAK-STAT signaling pathways, Integrin-mediated signal transduction, Ras signaling pathways, apoptosis and phagocytosis were positively modified after WSSV challenge. The proteolytic cascades including Complement-like activation and Hemolymph coagulations likely participated in antiviral immune response. The transcriptome data from hepatopancreas of Kuruma shrimp under WSSV challenge provided comprehensive information for identifying novel immune related genes in this valuable crustacean species despite the absence of the genome database of crustaceans.

Transcriptomic information from Pacific white shrimp (Litopenaeus vannamei) ovary and eyestalk, and expression patterns for genes putatively involved in the reproductive process

The increased use of massive sequencing technologies has enabled the identification of several genes known to be involved in different mechanisms associated with reproduction that so far have only been studied in vertebrates and other model invertebrate species. In order to further investigate the genes involved in Litopenaeus vannamei reproduction, cDNA and SSH libraries derived from female eyestalk and gonad were produced, allowing the identification of expressed sequences tags (ESTs) that potentially have a role in the regulation of gonadal maturation. In the present study, different transcripts involved in reproduction were identified and a number of them were characterized as full-length. These transcripts were evaluated in males and females in order to establish their tissue expression profiles during developmental stages (juvenile, subadult and adult), and in the case of females, their possible association with gonad maturation was assessed through expression analysis of vitellogenin. The results indicated that the expression of vitellogenin receptor (vtgr) and minichromosome maintenance (mcm) family members in the female gonad suggest an important role during previtellogenesis. Additionally, the expression profiles of genes such as famet, igfbp and gpcr in brain tissues suggest an interaction between the insulin/insulin-like growth factor signaling pathway (IIS) and methyl farnesoate (MF) biosynthesis for control of reproduction. Furthermore, the specific expression pattern of farnesoic acid O-methyltransferase suggests that final synthesis of MF is carried out in different target tissues, where it is regulated by esterase enzymes under a tissue-specific hormonal control. Finally, the presence of a vertebrate type steroid receptor in hepatopancreas and intestine besides being highly expressed in female gonads, suggest a role of that receptor during sexual maturation.

Transcriptome profiles of the clam Meretrix petechialis hepatopancreas in response to Vibrio infection

Microbial diseases have received much attention due to their enormous destruction of aquaculture, and Vibrio parahaemolyticus is one of the main pathogens that cause bacterial disease in the clam Meretrix petechialis. To better understand the molecular mechanisms of the immune response to Vibrio in M. petechialis, RNA-Seq was applied to explore global expression changes of hepatopancreas from this clam after Vibrio challenge. There were 199,318,966 clean reads obtained by Illumina sequencing, which were further assembled into 214,577 transcripts, and then 147,255 unigenes with an N50 of 1393 bp were identified. Gene ontology (GO) analysis revealed 21 biological process subcategories, 15 cellular component subcategories and 12 molecular function subcategories. A total of 8358 unigenes were mapped onto 267 biological signaling pathways by KEGG, among which there were 16 pathways related to the immune system. In total, 206 differentially expressed genes (DEGs) were identified, including 113 up-regulated unigenes and 93 down-regulated unigenes. In these DEGs, 96 DEGs were annotated in at least one database, accounting for 46.60% of all significant DEGs. To validate the transcriptome dataset, 15 DEGs were selected for real-time qPCR confirmation and the results showed that expression patterns of 13 genes (86.7%) agreed well with the RNA-Seq analysis. Fourteen of the 206 DEGs were annotated to be immune-related genes, and we examined the expression patterns of four immune-related DEGs using clams post immersion challenge. This study enriched the M. petechialis transcriptome database and provided insight into the immune response of M. petechialis against Vibrio infection.

Advances in the application of high-throughput sequencing in invertebrate virology

Over the last decade, advances in high-throughput sequencing technologies have revolutionised biological research, making it possible for DNA/RNA sequencing of any organism of interest to be undertaken. Sequencing approaches are now routinely used in the detection and characterisation of (novel) viruses, investigation of host-pathogen interactions, and effective development of disease treatment strategies. For the sequencing and identification of viruses of interest, metagenomics approaches using infected host tissue are frequently used, as it is not always possible to culture and isolate these pathogens. High-throughput sequencing can also be used to investigate host-pathogen interactions by investigating (temporal) transcriptomic responses of both the host and virus, potentially leading to the discovery of novel opportunities for treatment and drug targets. In addition, viruses in environmental samples (e.g. water or soil samples) can be identified using eDNA/metagenomics approaches. The promise that recent developments in sequencing brings to the field of invertebrate virology are not devoid of technical challenges, including the need for better laboratory and bioinformatics strategies to sequence and assemble virus genomes within complex tissue or environmental samples, and the difficulties associated with the annotation of the large number of novel viruses being discovered.

Aquaculture genomics, genetics and breeding in the United States: current status, challenges, and priorities for future research

Advancing the production efficiency and profitability of aquaculture is dependent upon the ability to utilize a diverse array of genetic resources. The ultimate goals of aquaculture genomics, genetics and breeding research are to enhance aquaculture production efficiency, sustainability, product quality, and profitability in support of the commercial sector and for the benefit of consumers. In order to achieve these goals, it is important to understand the genomic structure and organization of aquaculture species, and their genomic and phenomic variations, as well as the genetic basis of traits and their interrelationships. In addition, it is also important to understand the mechanisms of regulation and evolutionary conservation at the levels of genome, transcriptome, proteome, epigenome, and systems biology. With genomic information and information between the genomes and phenomes, technologies for marker/causal mutation-assisted selection, genome selection, and genome editing can be developed for applications in aquaculture. A set of genomic tools and resources must be made available including reference genome sequences and their annotations (including coding and non-coding regulatory elements), genome-wide polymorphic markers, efficient genotyping platforms, high-density and high-resolution linkage maps, and transcriptome resources including non-coding transcripts. Genomic and genetic control of important performance and production traits, such as disease resistance, feed conversion efficiency, growth rate, processing yield, behaviour, reproductive characteristics, and tolerance to environmental stressors like low dissolved oxygen, high or low water temperature and salinity, must be understood. QTL need to be identified, validated across strains, lines and populations, and their mechanisms of control understood. Causal gene(s) need to be identified. Genetic and epigenetic regulation of important aquaculture traits need to be determined, and technologies for marker-assisted selection, causal gene/mutation-assisted selection, genome selection, and genome editing using CRISPR and other technologies must be developed, demonstrated with applicability, and application to aquaculture industries.Major progress has been made in aquaculture genomics for dozens of fish and shellfish species including the development of genetic linkage maps, physical maps, microarrays, single nucleotide polymorphism (SNP) arrays, transcriptome databases and various stages of genome reference sequences. This paper provides a general review of the current status, challenges and future research needs of aquaculture genomics, genetics, and breeding, with a focus on major aquaculture species in the United States: catfish, rainbow trout, Atlantic salmon, tilapia, striped bass, oysters, and shrimp. While the overall research priorities and the practical goals are similar across various aquaculture species, the current status in each species should dictate the next priority areas within the species. This paper is an output of the USDA Workshop for Aquaculture Genomics, Genetics, and Breeding held in late March 2016 in Auburn, Alabama, with participants from all parts of the United States.

Two hemocyte sub-populations of kuruma shrimp Marsupenaeus japonicus

Hemocytes in the circulating hemolymph play important roles for immune responses in shrimp. Previous studies on immune responses by hemocytes in penaeid shrimp were based on gene expression analyses of the entire population of hemocytes and thus may have missed different immune responses of different hemocyte sub-populations. In this study, we separated hemocytes into two sub-populations by Percoll gradient centrifugation, morphological characteristics of each population were then analyzed by May-Giemsa staining, flow cytometry, and FACSCalibur. Results showed hemocytes were divided into an upper layer basophilic, and lower layer of eosinophilic hemocytes. Basophilic hemocytes were larger in size compared to eosinophilic hemocytes, which were more granulated than the basophilic hemocytes. Transcriptome analysis was then conducted through RNA-seq analysis by Miseq, which revealed 16 differentially-expressed transcripts between the two sub-populations. In the upper-layer, the highly expressed transcripts that were homologous to immune-related genes that suggest hemocytes from this layer may play as the regulator of immune system and control the action of other cells to eliminate pathogen. On the other hand, transcripts that were highly expressed in the lower-layer were homologous to the antimicrobial peptide (AMP) crustin, which supports that hemocytes on this layer have granules as crustins are normally secreted from hemocyte granules. The high expression of crustin in the lower-layer also provides insight on the mechanism of the anti-microbial function, where hemocytes produce and store AMPs in its granules. These differentially expressed genes are potential hemocyte molecular markers, and among them we identified one of the highly expressed genes in the hemocytes from the upper-layer (c11736_g1) to be a promising candidate molecular marker predicted to be a surface molecule, which is a common characteristic for molecular markers.

Differentially expressed transcripts in stomach of Penaeus monodon in response to AHPND infection

Acute Hepatopancreatic Necrosis Disease (AHPND) is an emerging disease in aquacultured shrimp caused by a pathogenic strain of Vibrio parahaemolyticus. As with several pathogenic bacteria, colonization of the stomach appeared to be the initial step of the infection for AHPND-causing Vibrio. To understand the immune responses in the stomach of black tiger shrimp (Penaeus monodon), differentially expressed transcripts (DETs) in the stomach during V. parahaemolyticus strain 3HP (VP3HP) infection was examined using Ion Torrent sequencing. From the total 42,998 contigs obtained, 1585 contigs representing 1513 unigenes were significantly differentially expressed with 1122 and 391 unigenes up- and down-regulated, respectively. Among the DETs, there were 141 immune-related unigenes in 10 functional categories: antimicrobial peptide, signal transduction pathway, proPO system, oxidative stress, proteinases/proteinase inhibitors, apoptotic tumor-related protein, pathogen recognition immune regulator, blood clotting system, adhesive protein and heat shock protein. Expression profiles of 20 of 22 genes inferred from RNA sequencing were confirmed with the results from qRT-PCR. Additionally, a novel isoform of anti-lipopolysaccharide factor, PmALF7 whose transcript was induced in the stomach after challenge with VP3HP was discovered. This study provided a fundamental information on the molecular response in the shrimp stomach during the AHPND infection that would be beneficial for future research.

Close ecological relationship among species facilitated horizontal transfer of retrotransposons

Background

Horizontal transfer (HT) of genetic materials is increasingly being found in both animals and plants and mainly concerns transposable elements (TEs). Many crustaceans have big genome sizes and are thus likely to harbor high TE contents. Their habitat might offer them ample opportunities to exchange genetic materials with organisms that are ecologically close but taxonomically distant to them.

Results

In this study, we analyzed the transcriptome of Pacific white shrimp (Litopenaeus vannamei), an important economic crustacean, to explore traces of HT events. From a collection of newly assembled transcripts, we identified 395 high reliable TE transcripts, most of which were retrotransposon transcripts. One hundred fifty-seven of those transcripts showed highest similarity to sequences from non-arthropod organisms, including ray-finned fishes, mollusks and putative parasites. In total, 16 already known L. vannamei TE families are likely to be involved in horizontal transfer events. Phylogenetic analyses of 10 L. vannamei TE families and their homologues (protein sequences) revealed that L. vannamei TE families were generally more close to sequences from aquatic species. Furthermore, TEs from other aquatic species also tend to group together, although they are often distantly related in taxonomy. Sequences from parasites and microorganisms were also widely present, indicating their possible important roles in HT events. Expression profile analyses of transcripts in two NCBI BioProjects revealed that transcripts involved in HT events are likely to play important roles in antiviral immunity. More specifically, those transcripts might act as inhibitors of antiviral immunity.

Conclusions

Close ecological relationship, especially predation, might greatly facilitate HT events among aquatic species. This could be achieved through exchange of parasites and microorganisms, or through direct DNA flow. The occurrence of HT events may be largely incidental, but the effects could be beneficial for recipients.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0767-0) contains supplementary material, which is available to authorized users.

Transcriptome analysis of the unfolded protein response in hemocytes of Litopenaeus vannamei

In this study, Litopenaeus vannamei was injected with double-stranded RNA (dsRNA) against L. vannamei immunoglobulin heavy chain binding protein (LvBip) to activating UPR in the hemocytes, shirmps injected dsRNA against enhanced green fluorescence protein (eGFP) as control group. And genes expression in hemocytes of then were analyzed using Illumina Hiseq 2500 (PE100). By comparing the analyzed results, 1418 unigenes were significantly upregulated, and 596 unigenes were significantly down-regulated upon UPR. Analysis of the differentially expressed genes against known databases indicated that the distribution of gene pathways between the upregulated and down-regulated genes were substantially different. A total of 208 genes of UPR system were obtained, and 69 of them were differentially expressed between the two groups. Results also showed that L. vannamei UPR was involved in various metabolic processes, such as glycometabolism, lipid metabolism, amino acid metabolism, and nucleic acid metabolism. In addition, UPR was emgaged in immune-assicoated signaling pathways, such as NF-κB signaling pathway, NOD-like receptor signaling pathway, Hippo signaling pathway, p38 MAPK signaling pathway and Wnt signaling pathway in L. vannamei. These results improved our current understanding of the L. vannamei UPR, and highlighted its importance in cell homeostasis upon environmental stress.

Exploring developmental gene toolkit and associated pathways in a potential new model crustacean using transcriptomic analysis

The crustaceans are one of the largest, most diverse, and most successful groups of invertebrates. The diversity among the crustaceans is also reflected in embryonic development models. However, the molecular genetics that regulates embryonic development is not known in those crustaceans that have a short germ-band development with superficial cleavage, such as Macrobrachium olfersi. This species is a freshwater decapod and has great potential to become a model for developmental biology, as well as for evolutionary and environmental studies. To obtain sequence data of M. olfersi from an embryonic developmental perspective, we performed de novo assembly and annotation of the embryonic transcriptome. Using a pooling strategy of total RNA, paired-end Illumina sequencing, and assembly with multiple k-mers, a total of 25,636,097 pair reads were generated. In total, 99,751 unigenes were identified, and 20,893 of these returned a Blastx hit. KEGG pathway analysis mapped a total of 6866 unigenes related to 129 metabolic pathways. In general, 21,845 unigenes were assigned to gene ontology (GO) categories: molecular function (19,604), cellular components (10,254), and biological processes (13,841). Of these, 2142 unigenes were assigned to the developmental process category. More specifically, a total of 35 homologs of embryonic development toolkit genes were identified, which included maternal effect (one gene), gap (six), pair-rule (six), segment polarity (seven), Hox (four), Wnt (eight), and dorsoventral patterning genes (three). In addition, genes of developmental pathways were found, including TGF-β, Wnt, Notch, MAPK, Hedgehog, Jak-STAT, VEGF, and ecdysteroid-inducible nuclear receptors. RT-PCR analysis of eight genes related to embryonic development from gastrulation to late morphogenesis/organogenesis confirmed the applicability of the transcriptome analysis.

Transcriptional Profiling of Banana Shrimp Fenneropenaeus merguiensis with Differing Levels of Viral Load

Viral pathogens are of serious concern to the culture of penaeid shrimps worldwide. However, little is known about the molecular response of shrimp to viral infection. Selective breeding has been suggested as an effective long-term strategy to manage viral disease, though more information on gene function is needed to help inform breeding programs. The study of cultured banana shrimp (Fenneropenaeus merguiensis) infected with hepatopancreatic parvo-like virus (HPV) provides a unique opportunity to explore the host response to viral infection independent of challenge testing. To gain insight into the genetic mechanisms underlying resistance to high levels of HPV, we examined hepatopancreas tissue from six full-sib groups of banana shrimp with differing levels of HPV infection for differences in gene expression. A total of 404 differentially expressed genes were identified with 180 being over-expressed and 224 under-expressed among high-HPV full-sib groups. Based on homology analysis, a large proportion of these genes were associated with processes reported to be involved in the immune response of crustaceans, including pattern recognition proteins, antimicrobial peptides, components of the prophenoloxidase system, and antiviral activity. The results indicate shrimp from high-HPV full-sib groups appear to have a lower presence of important immune response elements, yet possess upregulated putative antiviral pathways. Within the differentially expressed genes, over 4000 sequence variants were identified to be exclusive to either the high- or low-HPV full-sib groups. To our knowledge, this is the first report of differential expression analysis using RNA-Seq to explore differences in viral load among high- and low-HPV full-sib groups of cultured shrimp. This research has provided additional insight into our understanding of the mechanisms involved in the response of this shrimp species to a naturally occurring viral pathogen. Sequence variants identified in this study offer an exceptional resource for mining gene-associated markers of HPV resistance.

Gonadal transcriptomic analysis and differentially expressed genes in the testis and ovary of the Pacific white shrimp (Litopenaeus vannamei)

Background

The Pacific white shrimp (Litopenaeus vannamei) is the world’s most prevalent cultured crustacean species. However, the supply of high-quality broodstocks is limited and baseline information related to its reproductive activity and molecular issues related to gonad development are scarce. In this study, we performed transcriptome sequencing on the gonads of adult male and female L. vannamei to identify sex-related genes.

Results

A total of 25.16 gigabases (Gb) of sequences were generated from four L. vannamei gonadal tissue libraries. After quality control, 24.11 Gb of clean reads were selected from the gonadal libraries. De-novo assembly of all the clean reads generated a total of 65,218 unigenes with a mean size of 1021 bp and a N50 of 2000 bp. A search of all-unigene against Nr, SwissProt, KEGG, COG and NT databases resulted in 26,482, 23,062, 20,659, 11,935 and 14,626 annotations, respectively, providing a total of 30,304 annotated unigenes. Among annotated unigenes, 12,320 unigenes were assigned to gene ontology categories and 20,659 unigenes were mapped to 258 KEGG pathways. By comparing the ovary and testis libraries, 19,279 testicular up-regulated and 3,529 ovarian up-regulated unigenes were identified. Enrichment analysis of differentially expressed unigenes resulted in 1060 significantly enriched GO terms and 34 significantly enriched KEGG pathways. Nine ovary-specific, 6 testis-specific, 45 testicular up-regulated and 39 ovarian up-regulated unigenes were then confirmed by semi-quantitative PCR and quantitative real-time PCR. In addition, using all-unigenes as a reference, a total of 13,233 simple sequence repeats (SSRs) were identified in 10,411 unigene sequences.

Conclusions

The present study depicts the first large-scale RNA sequencing of shrimp gonads. We have identified many important sex-related functional genes, GO terms and pathways, all of which will facilitate future research into the reproductive biology of shrimp. We expect that the SSRs detected in this study can then be used as genetic markers for germplasm evaluation of breeding and imported populations.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2219-4) contains supplementary material, which is available to authorized users.

Transcriptomic analysis of Mandarin fish brain cells infected with infectious spleen and kidney necrosis virus with an emphasis on retinoic acid-inducible gene 1-like receptors and apoptosis pathways

Infectious spleen and kidney necrosis virus (ISKNV) has caused significant economic losses in the cultured Mandarin fish (Siniperca chuatsi) industry. The molecular mechanisms that underlie the pathogenesis of the viral infection remain poorly understood. In this study, deep RNA sequencing technique was used to analyze the transcriptomic profiles of Mandarin fish brain cells (CPB) at progressive time points after ISKNV infection. A total of 96,206,040 clean data from 98,235,240 sequence reads were obtained. These raw data were assembled into 66,787 unigenes. Among these unigenes, 33,225 and 29,210 had significant hit the Nr and SwissProt databases where they matched 27,537and 19,638 unique protein accessions, respectively. In the samples harvested at 24 or 72 h post of the infection, a total of 10,834 or 7584 genes were differentially expressed in infected CPB cells compared to non-infected cells, including 5445 or 3766 up-regulated genes and 5389 or 3818 down-regulated genes, respectively. In addition, 12 differentially expressed genes (DEGs) were validated by quantitative PCR. These DEGs were involved in many pathways of viral pathogenesis. Further analysis of the major DEGs genes involved in the RLRs and apoptosis pathways revealed some interesting findings. In the RLRs pathway, ISKNV infection inhibited the activation of NF-κB via over expression of the IKKB-α and IKKB-β and lessened expression of interleukin-1 receptor-associated kinase 4 (IRAK4). In the apoptosis pathway, ISKNV infection could induce apoptosis mainly via tumor necrosis factor (TNF) mediated extrinsic pathway. The cellular apoptosis induced by ISKNV infection was confirmed using annexinV-FITC/PI and DAPI staining methods.

De novo Transcriptome Assemblies of Rana (Lithobates) catesbeiana and Xenopus laevis Tadpole Livers for Comparative Genomics without Reference Genomes

In this work we studied the liver transcriptomes of two frog species, the American bullfrog (Rana (Lithobates) catesbeiana) and the African clawed frog (Xenopus laevis). We used high throughput RNA sequencing (RNA-seq) data to assemble and annotate these transcriptomes, and compared how their baseline expression profiles change when tadpoles of the two species are exposed to thyroid hormone. We generated more than 1.5 billion RNA-seq reads in total for the two species under two conditions as treatment/control pairs. We de novo assembled these reads using Trans-ABySS to reconstruct reference transcriptomes, obtaining over 350,000 and 130,000 putative transcripts for R. catesbeiana and X. laevis, respectively. Using available genomics resources for X. laevis, we annotated over 97% of our X. laevis transcriptome contigs, demonstrating the utility and efficacy of our methodology. Leveraging this validated analysis pipeline, we also annotated the assembled R. catesbeiana transcriptome. We used the expression profiles of the annotated genes of the two species to examine the similarities and differences between the tadpole liver transcriptomes. We also compared the gene ontology terms of expressed genes to measure how the animals react to a challenge by thyroid hormone. Our study reports three main conclusions. First, de novo assembly of RNA-seq data is a powerful method for annotating and establishing transcriptomes of non-model organisms. Second, the liver transcriptomes of the two frog species, R. catesbeiana and X. laevis, show many common features, and the distribution of their gene ontology profiles are statistically indistinguishable. Third, although they broadly respond the same way to the presence of thyroid hormone in their environment, their receptor/signal transduction pathways display marked differences.

Digital gene expression analysis in hemocytes of the white shrimp Litopenaeus vannamei in response to low salinity stress

The white shrimp Litopenaeus vannamei has been greatly impacted by low salinity stress. To gain knowledge on the immune response in L. vannamei under such stress, we investigated digital gene expression (DEG) in L. vannamei hemocytes using the deep-sequencing platform Illumina HiSeq 2000. In total, 38,155 high quality unigenes with average length 770 bp were generated; 145 and 79 genes were identified up- or down-regulated, respectively. Functional categorization and pathways of the differentially expressed genes revealed that immune signaling pathways, cellular immunity, humoral immunity, apoptosis, cellular protein synthesis, lipid transport and energy metabolism were the differentially regulated processes occurring during low salinity stress. These results will provide a resource for subsequent gene expression studies regarding environmental stress and a valuable gene information for a better understanding of immune mechanisms of L. vannamei under low salinity stress.

Novel transcriptome assembly and improved annotation of the whiteleg shrimp (Litopenaeus vannamei), a dominant crustacean in global seafood mariculture

We present a new transcriptome assembly of the Pacific whiteleg shrimp (Litopenaeus vannamei), the species most farmed for human consumption. Its functional annotation, a substantial improvement over previous ones, is provided freely. RNA-Seq with Illumina HiSeq technology was used to analyze samples extracted from shrimp abdominal muscle, hepatopancreas, gills and pleopods. We used the Trinity and Trinotate software suites for transcriptome assembly and annotation, respectively. The quality of this assembly and the affiliated targeted homology searches greatly enrich the curated transcripts currently available in public databases for this species. Comparison with the model arthropod Daphnia allows some insights into defining characteristics of decapod crustaceans. This large-scale gene discovery gives the broadest depth yet to the annotated transcriptome of this important species and should be of value to ongoing genomics and immunogenetic resistance studies in this shrimp of paramount global economic importance.

In-depth transcriptome analysis of the red swamp crayfish Procambarus clarkii

The red swamp crayfish Procambarus clarkii is a highly adaptable, tolerant, and fecund freshwater crayfish that inhabits a wide range of aquatic environments. It is an important crustacean model organism that is used in many research fields, including animal behavior, environmental stress and toxicity, and studies of viral infection. Despite its widespread use, knowledge of the crayfish genome is very limited and insufficient for meaningful research. This is the use of next-generation sequencing techniques to analyze the crayfish transcriptome. A total of 324.97 million raw reads of 100 base pairs were generated, and a total of 88,463 transcripts were assembled de novo using Trinity software, producing 55,278 non-redundant transcripts. Comparison of digital gene expression between four different tissues revealed differentially expressed genes, in which more overexpressed genes were found in the hepatopancreas than in other tissues, and more underexpressed genes were found in the testis and the ovary than in other tissues. Gene ontology (GO) and KEGG enrichment analysis of differentially expressed genes revealed that metabolite- and immune-related pathway genes were enriched in the hepatopancreas, and DNA replication-related pathway genes were enriched in the ovary and the testis, which is consistent with the important role of the hepatopancreas in metabolism, immunity, and the stress response, and with that of the ovary and the testis in reproduction. It was also found that 14 vitellogenin transcripts were highly expressed specifically in the hepatopancreas, and 6 transcripts were highly expressed specifically in the ovary, but no vitellogenin transcripts were highly expressed in both the hepatopancreas and the ovary. These results provide new insight into the role of vitellogenin in crustaceans. In addition, 243,764 SNP sites and 43,205 microsatellite sequences were identified in the sequencing data. We believe that our results provide an important genome resource for the crayfish.

De novo transcriptome sequencing analysis and comparison of differentially expressed genes (DEGs) in Macrobrachium rosenbergii in China

Giant freshwater prawn (GFP; Macrobrachium rosenbergii) is an exotic species that was introduced into China in 1976 and thereafter it became a major species in freshwater aquaculture. However the gene discovery in this species has been limited to small-scale data collection in China. We used the next generation sequencing technology for the experiment; the transcriptome was sequenced of samples of hepatopancreas organ in individuals from 4 GFP groups (A1, A2, B1 and B2). De novo transcriptome sequencing generated 66,953 isogenes. Using BLASTX to search the Non-redundant (NR), Search Tool for the Retrieval of Interacting Genes (STRING), and Kyoto Encyclopedia of Genes and Genome (KEGG) databases; 21,224 unigenes were annotated, 9,552 matched unigenes with the Gene Ontology (GO) classification; 5,782 matched unigenes in 25 categories of Clusters of Orthologous Groups of proteins (COG) and 20,859 unigenes were consequently assigned to 312 KEGG pathways. Between the A and B groups 147 differentially expressed genes (DEGs) were identified; between the A1 and A2 groups 6,860 DEGs were identified and between the B1 and B2 groups 5,229 DEGs were identified. After enrichment, the A and B groups identified 38 DEGs, but none of them were significantly enriched. The A1 and A2 groups identified 21,856 DEGs in three main categories based on functional groups: biological process, cellular_component and molecular function and the KEGG pathway defined 2,459 genes had a KEGG Ortholog-ID (KO-ID) and could be categorized into 251 pathways, of those, 9 pathways were significantly enriched. The B1 and B2 groups identified 5,940 DEGs in three main categories based on functional groups: biological process, cellular_component and molecular function, and the KEGG pathway defined 1,543 genes had a KO-ID and could be categorized into 240 pathways, of those, 2 pathways were significantly enriched. We investigated 99 queries (GO) which related to growth of GFP in 4 groups. After enrichment we identified 23 DEGs and 1 KEGG PATHWAY 'ko04711' relation with GFP growth.

Comparative transcriptomic characterization of the early development in Pacific white shrimp Litopenaeus vannamei

Penaeid shrimp has a distinctive metamorphosis stage during early development. Although morphological and biochemical studies about this ontogeny have been developed for decades, researches on gene expression level are still scarce. In this study, we have investigated the transcriptomes of five continuous developmental stages in Pacific white shrimp (Litopenaeus vannamei) with high throughput Illumina sequencing technology. The reads were assembled and clustered into 66,815 unigenes, of which 32,398 have putative homologues in nr database, 14,981 have been classified into diverse functional categories by Gene Ontology (GO) annotation and 26,257 have been associated with 255 pathways by KEGG pathway mapping. Meanwhile, the differentially expressed genes (DEGs) between adjacent developmental stages were identified and gene expression patterns were clustered. By GO term enrichment analysis, KEGG pathway enrichment analysis and functional gene profiling, the physiological changes during shrimp metamorphosis could be better understood, especially histogenesis, diet transition, muscle development and exoskeleton reconstruction. In conclusion, this is the first study that characterized the integrated transcriptomic profiles during early development of penaeid shrimp, and these findings will serve as significant references for shrimp developmental biology and aquaculture research.

RNA-seq as a powerful tool for penaeid shrimp genetic progress

The sequences of all different RNA transcripts present in a cell or tissue that are related to the gene expression and its functional control represent what it is called a transcriptome. The transcripts vary between cells, tissues, ontogenetic and environmental conditions, and the knowledge that can be gained through them is of a solid relevance for genetic applications in aquaculture. Some of the techniques used in transcriptome studies, such as microarrays, are being replaced for next-generation sequencing approaches. RNA-seq emerges as a new possibility for the transcriptome complexity analysis as well as for the candidate genes and polymorphisms identification of penaeid species. Thus, it may also help to understand the determination of complex traits mechanisms and genetic improvement of stocks. In this review, it is first introduced an overview of transcriptome analysis by RNA-seq, followed by a discussion of how this approach may be applied in genetic progress within penaeid stocks.

The promotion of cytoskeleton integration and redox in the haemocyte of shrimp Litopenaeus vannamei after the successive stimulation of recombinant VP28

VP28 protein has been reported to work as a "vaccine" to protect the host from white spot syndrome, but the detailed mechanism of vaccination with VP28 protein in shrimp is still far from well understood. In the present study, whole transcriptomes of shrimp haemocytes were sequenced using the SOLiD4 platform after the successive VP28 stimulation. Eight single-end fragment libraries were constructed and sequenced in the four groups including the VP28-VP28, PBS-VP28, PBS-PBS and BLANK group, and there were 243,949,667 single-end reads with length of 50bp obtained totally, with 14,800 genes further identified. After reads mapping and transcript assembling, 1027, 1539, 1158, 1091 and 1300 genes in five differentially expressed gene lists were obtained in the comparison of VP28-VP28 versus PBS-VP28, VP28-VP28 versus PBS-PBS, VP28-VP28 versus BLANK, PBS-VP28 versus PBS-PBS and PBS-VP28 versus BLANK, respectively. There were 555 differentially expressed genes responsive to the single VP28 stimulation after grouping the PBS-VP28_BLANK and PBS-VP28_PBS-PBS gene lists, and 269 ones responsive to the successive VP28 stimulation after grouping the VP28-VP28_BLANK, VP28-VP28_PBS-PBS and VP28-VP28_PBS-VP28 gene lists. In the GO enrichment analysis of the genes responsive to the single VP28 stimulation, five immune-related GO terms were observed among 14 increased terms, which included defense response to bacterium, response to stimulus, disruption of cells of other organism, killing of cells of other organism and response to bacterium. It was worth noting that the GO terms, response to stimulus and response to stress, were the most common annotation ones which accounted 28.7% and 18.8% of the total differently expressed genes, respectively. For the genes responsive to the successive VP28 stimulation, terms including actin filament-based movement and myosin heavy chain binding were mostly enriched in the Biological Process and Molecular Function category, respectively. In the Cellular Component category, the enriched GO terms were myosin VII complex, myosin V complex, myosin VI complex and myosin II complex. Furthermore, the most abundant GO term was oxidation-reduction process, followed by single-organism transport, neurogenesis and translation for 214 genes only responsive to successive VP28 stimulation. These results collectively indicated that the successive VP28 stimulation could modulate cytoskeleton integration and redox to promote the phagocytosis activity of shrimp haemocytes, which might protect effectively for shrimp against WSSV infection.

Antiviral defense in shrimp: from innate immunity to viral infection

The culture of penaeid shrimp is rapidly developing as a major business endeavor worldwide. However, viral diseases have caused huge economic loss in penaeid shrimp culture industries. Knowledge of shrimp innate immunity and antiviral responses has made important progress in recent years, allowing the design of better strategies for the prevention and control of shrimp diseases. In this study, we have updated information on shrimp antiviral immunity and interactions between shrimp hosts and viral pathogens. Current knowledge and recent progress in immune signaling pathways (e.g., Toll/IMD-NF-κB and JAK-STAT signaling pathways), RNAi, phagocytosis, and apoptosis in shrimp antiviral immunity are discussed. The mechanism of viral infection in shrimp hosts and the interactions between viruses and shrimp innate immune systems are also analyzed.