Comparative Transcriptome Analysis in the Hepatopancreas Tissue of Pacific White Shrimp Litopenaeus vannamei Fed Different Lipid Sources at Low Salinity

Transcriptomic variations associated with salinity stress in Penaeus indicus

BACKGROUND

The Indian white shrimp, Penaeus indicus a native species of India, is important brackishwater aquaculture species. Shrimps are euryhaline in nature and they regulate osmotic and ionic concentrations by osmoregulatory process. However, variations in abiotic factors such as salinity result in stress to the shrimps during culture period affecting their growth and immunity.

METHODS AND RESULTS

To understand the adaptive mechanism to stress in low salinity conditions, RNA-seq was used to compare the transcriptomic response of P. indicus upto 3 weeks. De novo assembly using Trinity assembler generated a total of 173,582 transcripts. The assembly had a mean length of 854 bp, N50 value of 1243 bp and GC content of 42.33%. Differential gene expression analysis, resulted in identification of 2130, 3090, and 5351 DEGs in 7 days, 14 days and 21 days respectively of salinity stress period. The pathway prediction of the assembled trinity transcripts using KEGG database showed total number of 329 pathways linking 12,430 transcripts. KEGG pathway enrichment analyses led to the identification of several enriched pathways related to lipid metabolism, amino acid metabolism, glycolysis, signalling pathways etc. Selected genes involved in osmoregulatory process and immune response in shrimps were validated and analysed for the gene expression levels by quantitative real-time PCR (qPCR).

CONCLUSION

This study on the adaptive transcriptomic response of P. indicus to low salinity, will further help in our understanding of the molecular mechanisms underlying osmoregulation mechanism in shrimps.

Single Nucleotide Polymorphism Markers with Applications in Conservation and Exploitation of Aquatic Natural Populations

Simple Summary

In recent decades, societies, states and local authorities have become increasingly aware that for effective long-term management and protection of aquatic ecosystems and populations, it is necessary to take into account the genetic changes occurring in these populations. One type of high-resolution molecular marker suitable for studying the neutral and adaptive genetic diversity of populations is single nucleotide polymorphism (SNP). This review is an attempt to show the benefits of using SNPs to recognize natural populations of aquatic animals and detect the threats to them from accidentally or intentionally released farm animals, fishery and global climate changes. It is postulated that conservation actions should protect not only pristine natural populations that are endangered or overfished, but also protect populations of non-threatened species from unnecessarily released semi-domesticated animals. The enhancement of natural populations with farmed material usually reduces their genetic diversity. Experimental size-selective catches of artificially created populations have caused evolutionary changes in the life cycles of fishes. However, fishery-induced evolution in natural populations is difficult to observe. The negative measurable effects on populations can be expected when the number of breeding individuals is reduced below 100, which occurs very rarely in the sea and more often in fragmented freshwater streams, ponds and seasonal rivers.

Abstract

An increasing number of aquatic species have been studied for genetic polymorphism, which extends the knowledge on their natural populations. One type of high-resolution molecular marker suitable for studying the genetic diversity of large numbers of individuals is single nucleotide polymorphism (SNP). This review is an attempt to show the range of applications of SNPs in studies of natural populations of aquatic animals. In recent years, SNPs have been used in the genetic analysis of wild and enhanced fish and invertebrate populations in natural habitats, exploited migratory species in the oceans, migratory anadromous and freshwater fish and demersal species. SNPs have been used for the identification of species and their hybrids in natural environments, to study the genetic consequences of restocking for conservation purposes and the negative effects on natural populations of fish accidentally escaping from culture. SNPs are very useful for identifying genomic regions correlated with phenotypic variants relevant for wildlife protection, management and aquaculture. Experimental size-selective catches of populations created in tanks have caused evolutionary changes in life cycles of fishes. The research results have been discussed to clarify whether the fish populations in natural conditions can undergo changes due to selective harvesting targeting the fastest-growing fishes.

Multi-Omics Analysis of Lipid Metabolism for a Marine Probiotic Meyerozyma guilliermondii GXDK6 Under High NaCl Stress

Investigating microbial lipid regulation contributes to understanding the lipid-dependent signal transduction process of cells and helps to improve the sensitivity of microorganisms to environmental factors by interfering with lipid metabolism, thus beneficial for constructing advanced cell factories of novel molecular drugs. Integrated omics technology was used to systematically reveal the lipid metabolism mechanism of a marine Meyerozyma guilliermondii GXDK6 under high NaCl stress and test the sensitivity of GXDK6 to antibiotics when its lipid metabolism transformed. The omics data showed that when GXDK6 perceived 10% NaCl stress, the expression of AYR1 and NADPH-dependent 1-acyldihydroxyacetone phosphate reductase was inhibited, which weaken the budding and proliferation of cell membranes. This finding was further validated by decreased 64.39% of OD₆₀₀ under 10% NaCl stress when compared with salt-free stress. In addition, salt stress promoted a large intracellular accumulation of glycerol, which was also verified by exogenous addition of glycerol. Moreover, NaCl stress remarkably inhibited the expression of drug target proteins (such as lanosterol 14-alpha demethylase), thereby increasing sensitivity to fluconazole. This study provided new insights into the molecular mechanism involved in the regulation of lipid metabolism in Meyerozyma guilliermondii strain and contributed to developing new methods to improve the effectiveness of killing fungi with lower antibiotics.

Multi-omics analysis to examine microbiota, host gene expression and metabolites in the intestine of black tiger shrimp (Penaeus monodon) with different growth performance

Understanding the correlation between shrimp growth and their intestinal bacteria would be necessary to optimize animal's growth performance. Here, we compared the bacterial profiles along with the shrimp's gene expression responses and metabolites in the intestines between the Top and the Bottom weight groups. Black tiger shrimp (Penaeus monodon) were collected from the same population and rearing environments. The two weight groups, the Top-weight group with an average weight of 36.82 ± 0.41 g and the Bottom-weight group with an average weight of 17.80 ± 11.81 g, were selected. Intestines were aseptically collected and subjected to microbiota, transcriptomic and metabolomic profile analyses. The weighted-principal coordinates analysis (PCoA) based on UniFrac distances showed similar bacterial profiles between the two groups, suggesting similar relative composition of the overall bacterial community structures. This observed similarity was likely due to the fact that shrimp were from the same genetic background and reared under the same habitat and diets. On the other hand, the unweighted-distance matrix revealed that the bacterial profiles associated in intestines of the Top-weight group were clustered distinctly from those of the Bottom-weight shrimp, suggesting that some unique non-dominant bacterial genera were found associated with either group. The key bacterial members associated to the Top-weight shrimp were mostly from Firmicutes (Brevibacillus and Fusibacter) and Bacteroidetes (Spongiimonas), both of which were found in significantly higher abundance than those of the Bottom-weight shrimp. Transcriptomic profile of shrimp intestines found significant upregulation of genes mostly involved in nutrient metabolisms and energy storage in the Top-weight shrimp. In addition to significantly expressed metabolic-related genes, the Bottom-weight shrimp also showed significant upregulation of stress and immune-related genes, suggesting that these pathways might contribute to different degrees of shrimp growth performance. A non-targeted metabolome analysis from shrimp intestines revealed different metabolic responsive patterns, in which the Top-weight shrimp contained significantly higher levels of short chain fatty acids, lipids and organic compounds than the Bottom-weight shrimp. The identified metabolites included those that were known to be produced by intestinal bacteria such as butyric acid, 4-indolecarbaldehyde and L-3-phenyllactic acid as well as those produced by shrimp such as acyl-carnitines and lysophosphatidylcholine. The functions of these metabolites were related to nutrient absorption and metabolisms. Our findings provide the first report utilizing multi-omics integration approach to investigate microbiota, metabolic and transcriptomics profiles of the host shrimp and their potential roles and relationship to shrimp growth performance.

Transcriptome analysis of metabolisms related to fruit cracking during ripening of a cracking-susceptible grape berry cv. Xiangfei (Vitis vinifera L.)

BACKGROUND

Grape is an economically valuable fruit around the world. However, some cultivars are prone to fruit cracking during ripening, leading to severe losses.

OBJECTIVE

We aimed to find important metabolisms related to fruit cracking during ripening process.

METHODS

RNA-Sequence and analysis was applied to the pericarp of cracking-susceptible 'Xiang Fei' at 1 (W1), 2 (W2) and 3 weeks (W3) after veraison on Illumina HiSeq xten; RESULTS: Compared with W1, the berry cracking rate increased significantly in W2 and W3. Through transcriptomic analysis, a total of 22,609 genes were expressed in the grape pericarp, among which 805 and 2758 genes were significantly differentially regulated in W1-vs.-W2 and W1-vs.-W3 comparison, respectively. Besides, 304 and 354 genes were up- and down-regulated in both comparisons. The significantly enriched GO terms of both W1-W2 and W1-W3 are related to cell wall and wax biosynthesis. And lipid metabolism, which are involved in the top 20 enriched KEGG pathways of both comparisons, was related to wax biosynthesis. Further, GO enrichment analysis of differentially expressed genes (DEGs) with same regulatory changes also indicated that the continuously up-regulated DEGs are significantly enriched in cell wall component biosynthesis and hydrolase.

CONCLUSION

These findings suggested that genes related to cell wall metabolism and cuticle biosynthesis may play important roles in regulating grape berry cracking. Our results provide a reference for further studies on the molecular mechanism underlying fruit cracking.

Transcriptome Analysis Reveals the Tolerance Mechanism of Mantis Shrimp (Oratosquilla oratoria) under a Lipopolysaccharide Challenge

Lipopolysaccharide (LPS), a major cell wall component of Gram-negative bacteria, is considered to lead to some disease development in commercial crustaceans. However, mantis shrimps Oratosquilla oratoria (Crustacea: Stomatopoda) have a strong vitality and ability to resist disease. To study the tolerance mechanism of mantis shrimp, transcriptome analyses were conducted in hepatopancreas of O. oratoria under LPS challenge investigation. Totally, 84 547 044 clean reads were obtained from transcriptomes (43 159 230 in OP (control), 41 387 814 in OL (treatment), respectively). Unigenes, the longest transcript of each gene, with a total length of 68 318 880 bp and the total number of 100 978 were obtained. 8369 (8.28%) of unigenes were successfully annotated in all databases and 54 888 (54.35%) were annotated in at least one database. Finally, 1012 differentially expressed genes (DEGs) including 439 and 573 showed significantly upregulated and downregulated were determined between OL and OP, respectively. Moreover, those DEGs only expressed in OL or OP accounted for 8.99%. The functional classification based on GO and KEGG indicated that the common enrichment categories for the DEGs are "amino sugar metabolic" and "cellular homeostasis" and that the progress of nutrient metabolic and homeostasis in cells is important in facing variable environmental conditions. Protein-protein interaction analysis elucidated proteins, β-actin (ACTB_G1), T-complex protein subunits (TCPs), heat shock proteins (HSPs), hydroxysteroid dehydrogenase-like protein 2 (HSDL2), kinesin family member 5 (KIF5), methylglutaconyl-CoA hydratase (AUH), and myosin heavy chain (MYH) may play key roles in response to an LPS challenge. This study laid a foundation to further investigate the possible adaptation way that O. oratoria survives in a bacterial challenge.

Regulation of Dietary Lipid Sources on Tissue Lipid Classes and Mitochondrial Energy Metabolism of Juvenile Swimming Crab, Portunus trituberculatus

An 8-weeks feeding trial with swimming crab, Portunus trituberculatus, was conducted to investigate the effects of different dietary lipid sources on the lipid classes, lipid metabolism, and mitochondrial energy metabolism relevant genes expression. Six isonitrogenous and isolipidic experimental diets were formulated to contain fish oil (FO), krill oil (KO), palm oil (PO), rapeseed oil (RO), soybean oil (SO), and linseed oil (LO), respectively. A total of 270 swimming crab juveniles (initial weight 5.43 ± 0.03 g) were randomly divided into six diets with three replications, each consisted of 45 juvenile crabs. The results revealed that crabs fed KO had highest lipid content in hepatopancreas and free fatty acids in serum among all diets. The anabolic pathway relevant genes: fas and acc were up-regulated in KO diet. The catabolic pathway relevant genes, hsl, was up-regulated in LO diet, while cpt1 was up-regulated in KO diet. Whereas, the genes involved in the transport and uptake of fatty acids such as fabp1 and fatp4 were down-regulated in crab fed PO and RO diets. Furthermore, the gene expression levels of transcription factors: srebp-1 and hnf4α in KO and SO diets were the highest among all diets. FO and KO diets had significantly higher unsaturation index of mitochondrial membrane than others. The genes related to mitochondrial energy metabolism, such as Atpase6, sirt1, and sirt3 were significantly up-regulated in KO and SO diets. In summary, dietary KO and SO supplementation could improve the lipid metabolism, promote energy production for juvenile swimming crab and improve physiological process and function including molting. These findings could contribute to deepen the understanding of the physiological metabolism of dietary fatty acids for swimming crab.

Comparative transcriptome analysis of Eriocheir japonica sinensis response to environmental salinity

Chinese mitten crabs (Eriocheir japonica sinensis) are catadromous, spending most of their lives in fresh water, but moving to a mixed salt-fresh water environment for reproduction. The characteristics of this life history might imply a rapidly evolutionary transition model for adaptation to marine from freshwater habitats. In this study, transcriptome-wide identification and differential expression on Chinese mitten crab groups were analysed. Results showed: clean reads that were obtained totalled 93,833,096 (47,440,998 in Group EF, the reference, and 46,392,098 in Group ES, the experimental) and 14.08G (7.12G in Group EF 6.96G in Group ES); there were 11,667 unigenes (15.29%) annotated, and they were located to 230 known KEGG pathways in five major categories; in differential expression analysis, most of the top 20 up-regulated pathways were connected to the immune system, disease, and signal transduction, while most of the top 20 down-regulated pathways were related to the metabolism system; meanwhile, 8 representative osmoregulation-related genes (14-3-3 epsilon, Cu2+ transport ATPase, Na+/K+ ATPase, Ca2+ transporting ATPase, V-ATPase subunit A, Putative arsenite-translocating ATPase, and Cation transport ATPase, Na+/K+ symporter) showed up-regulation, and 1 osmoregulation-related gene (V-ATPase subunit H) showed down-regulation. V-ATPase subunit H was very sensitive to the transition of habitats. These results were consistent with the tests of qRT-PCR. The present study has provided a foundation to further understand the molecular mechanism in response to salinity changing in water.

Transcriptome Analysis of Hepatopancreas from the Cr (VI)-Stimulated Mantis Shrimp ( Oratosquilla oratoria) by Illumina Paired-End Sequencing: Assembly, Annotation, and Expression Analysis

Cr (VI), the pathogenicity factor, is widely known to cause toxic effects in living organisms. Given the economic importance of the mantis shrimp ( Oratosquilla oratoria), the understanding of impacts by Cr (VI) is considered important. In this study, transcriptome of mantis shrimp was characterized by a comparison between control and Cr (VI)-treated samples using RNA-seq approach. Totally, 88 234 826 bp and 13.24G clean reads were obtained. The total length and number of unigenes were 68 411 206 bp and 100 918, respectively. The maximal and average length of unigenes was 24 906 bp and 678 bp, respectively (N50, 798 bp). 7115 of these unigenes accounted for 7.05% of the total that were annotated in all databases. After annotation of assembled unigenes, 35 619 of them were assigned into 3 functional categories and 56 subcategories using Gene Ontology; 18 580 of them were assigned into 26 functional categories using Clusters of Orthologous Groups of proteins; 16 864 of them were assigned into 5 major categories and 32 subclasses using KEGG. Finally, 1730 genes were differentially expressed (DGEs), 9 up-regulated pathways (protein digestion and absorption, neuroactive ligand-receptor interaction, pancreatic secretion, tyrosine metabolism, amoebiasis, ECM-receptor interaction, riboflavin metabolism, amino sugar and nucleotide sugar metabolism and AGE-RAGE signaling pathway in diabetic complications) were significantly enriched ( q < 0.05), and one down-regulated pathway ( Staphylococcus aureus infection) was significantly enriched ( q < 0.05). Up-regulation of genes in pathways of protein digestion/absorption ( PepT1/SLC15A and ATP1B) and environment information processing ( COL1AS, COL4A; LAMA3_5, LAMB3; FN1 and TN) may imply the potentially positive toxicity resistance mechanisms.

Phenotypic plasticity in gene expression and physiological response in red drum Sciaenops ocellatus exposed to a long-term freshwater environment

Red drum (Sciaenops ocellatus) is a euryhaline fish commonly found in the Gulf of Mexico and along the Atlantic coast of North America. Because of high commercial demand and its euryhaline characteristics, aquaculture of this species has diversified from marine to low-salinity aquaculture systems. In recent years, interest in the feasibility of producing red drum in inland freshwater systems has grown and this prompted us to investigate its osmoregulatory capacity after rearing for 8 months in a freshwater aquaculture system. We compared the activities of several genes and enzymes involved in the osmoregulatory process in freshwater-acclimatized (FW) and seawater (SW) red drum. The gene expression profiles were variable: the expression of genes encoding Na⁺/K⁺-ATPase (NKA) and the cystic fibrosis transmembrane regulator (CFTR) was slightly higher in SW than FW fish, while phosphoenolpyruvate carboxykinase (PEPCK) and the glucocorticoid receptor messenger RNA (mRNA) levels were higher in FW red drum. The total plasma K concentration was 60.3% lower, and gill NKA activity was 63.5% lower in FW than in SW fish. PEPCK activity was twofold higher in FW than in SW red drum. Similarly, liver glycogen was 60% higher in FW fish. In summary, both gene expression and the enzyme activity data support the phenotypic plasticity of red drum and suggest that the limited capacity for ion homeostasis observed, in particular the low plasma K concentration, was due to the composition of freshwater and does not necessarily reflect a physiological inability to osmoregulate.

Effects of dietary lipids on the hepatopancreas transcriptome of Chinese mitten crab (Eriocheir sinensis)

Fish oil supplies worldwide have declined sharply over the years. To reduce the use of fish oil in aquaculture, many studies have explored the effects of fish oil substitutions on aquatic animals. To illustrate the effects of dietary lipids on Chinese mitten crab and to improve the use of vegetable oils in the diet of the crabs, 60 male juvenile Chinese mitten crabs were fed one of five diets for 116 days: fish oil (FO), soybean oil (SO), linseed oil (LO), FO + SO (1:1, FSO), and FO + LO (1:1, FLO). Changes in the crab hepatopancreas transcriptome were analyzed using RNA sequencing. There were a total 55,167 unigenes obtained from the transcriptome, of which the expression of 3030 was significantly altered in the FLO vs. FO groups, but the expression of only 412 unigenes was altered in the FSO vs. FO groups. The diets significantly altered the expression of many enzymes involved in lipid metabolism, such as pancreatic lipase, long-chain acyl-CoA synthetases, carnitine palmitoyltransferase I, acetyl-CoA carboxylase, fatty acid synthase, and fatty acyl Δ9-desaturase. The dietary lipids also affected the Toll-like receptor and Janus activated kinase-signal transducers and activators of transcription signaling pathways. Our results indicate that substituting fish oil with vegetable oils in the diet of Chinese mitten crabs might decrease the digestion and absorption of dietary lipids, fatty acids biosynthesis, and immunologic viral defense, and increase β-oxidation by altering the expression of the relevant genes. Our results lay the foundation for further understanding of lipid nutrition in Chinese mitten crab.

Comparative transcriptome analysis of Rimicaris sp. reveals novel molecular features associated with survival in deep-sea hydrothermal vent

Shrimp of the family Alvinocarididae are the predominant megafauna of deep-sea hydrothermal vents. However, genome information on this family is currently unavailable. In the present study, by employing Illumina sequencing, we performed the first de novo transcriptome analysis of the gills of the shrimp Rimicaris sp. from the hydrothermal vent in Desmos, Manus Basin. The analysis was conducted in a comparative manner with the shrimp taken directly from the vent (GR samples) and the shrimp that had been maintained for ten days under normal laboratory condition (mGR samples). Among the 128,938 unigenes identified, a large number of differentially expressed genes (DEGs) between the GR and mGR samples were detected, including 2365 and 1607 genes significantly upregulated and downregulated, respectively, in GR. The DEGs covered diverse functional categories. Most of the DEGs associated with immunity were downregulated in GR, while most of the DEGs associated with sulfur metabolism and detoxification were upregulated in GR. These results provide the first comprehensive transcriptomic resource for hydrothermal vent Rimicaris and revealed varied categories of genes likely involved in deep-sea survival.

Purification and characterization of bioactive peptides RYamide and CCHamide in the kuruma shrimp Marsupenaeus japonicus

To understand the regulation systems of appetite, bioactive peptides from the kuruma shrimp Marsupenaeus japonicus (Mj) were isolated and purified by reverse pharmacological assays using CHO cells expressing the Drosophila melanogaster G-protein-coupled receptors (GPCRs) CG5811 (a RYamide receptor) or CG14593 (a CCHamide-2 receptor). Four peptides having binding activity to GPCRs were obtained and named Mj RYamide-1, Mj RYamide-2, Mj RYamide-3, and Mj CCHamide. Genes encoding the prepropeptides of these peptides were identified using kuruma shrimp transcriptome databases. The Mj prepro-RYamide gene encodes a 130-amino acid polypeptide containing Mj RYamide-1, Mj RYamide-2, and Mj RYamide-3, whereas the Mj prepro-CCHamide gene encodes a 119-amino acid polypeptide containing a single Mj CCHamide peptide. The expression of these genes was confirmed in various neuronal organs including the brain and ventral nerve cord. In addition, prepro-RYamide gene expression is significantly reduced in the brain after starvation. RYamides may thus be associated with regulation of feeding or digestion. Changes in kayak (the c-fos ortholog in invertebrates) gene expression after administration of synthetic peptides were also investigated. Mj kayak expression levels are upregulated in hepatopancreas after treatment with Mj RYamide-3 or CCHamide. Thus, the peptides isolated in this study may have some regulatory effect on cellular metabolism in aquacultured invertebrates.

Comparative transcriptome and potential antiviral signaling pathways analysis of the gills in the red swamp crayfish, Procambarus clarkii infected with White Spot Syndrome Virus (WSSV)

Red swamp crayfish is an important model organism for research of the invertebrate innate immunity mechanism. Its excellent disease resistance against bacteria, fungi, and viruses is well-known. However, the antiviral mechanisms of crayfish remain unclear. In this study, we obtained high-quality sequence reads from normal and white spot syndrome virus (WSSV)-challenged crayfish gills. For group normal (GN), 39,390,280 high-quality clean reads were randomly assembled to produce 172,591 contigs; whereas, 34,011,488 high-quality clean reads were randomly assembled to produce 182,176 contigs for group WSSV-challenged (GW). After GO annotations analysis, a total of 35,539 (90.01%), 14,931 (37.82%), 28,221 (71.48%), 25,290 (64.05%), 15,595 (39.50%), and 13,848 (35.07%) unigenes had significant matches with sequences in the Nr, Nt, Swiss-Prot, KEGG, COG and GO databases, respectively. Through the comparative analysis between GN and GW, 12,868 genes were identified as differentially up-regulated DEGs, and 9,194 genes were identified as differentially down-regulated DEGs. Ultimately, these DEGs were mapped into different signaling pathways, including three important signaling pathways related to innate immunity responses. These results could provide new insights into crayfish antiviral immunity mechanism.

Transcriptome Analysis of the Tadpole Shrimp (Triops longicaudatus) by Illumina Paired-End Sequencing: Assembly, Annotation, and Marker Discovery

The tadpole shrimp (Triops longicaudatus) is an aquatic crustacean that helps control pest populations. It inhabits freshwater ponds and pools and has been described as a living fossil. T. longicaudatus was officially declared an endangered species South Korea in 2005; however, through subsequent protection and conservation management, it was removed from the endangered species list in 2012. The limited number of available genetic resources on T. longicaudatus makes it difficult to obtain valuable genetic information for marker-aided selection programs. In this study, whole-transcriptome sequencing of T. longicaudatus generated 39.74 GB of clean data and a total of 269,822 contigs using the Illumina HiSeq 2500 platform. After clustering, a total of 208,813 unigenes with an N₅₀ length of 1089 bp were generated. A total of 95,105 unigenes were successfully annotated against Protostome (PANM), Unigene, Eukaryotic Orthologous Groups (KOG), Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases using BLASTX with a cut-off of 1E−5. A total of 57,731 unigenes were assigned to GO terms, and 7247 unigenes were mapped to 129 KEGG pathways. Furthermore, 1595 simple sequence repeats (SSRs) were detected from the unigenes with 1387 potential SSR markers. This is the first report of high-throughput transcriptome analysis of T. longicaudatus, and it provides valuable insights for genetic research and molecular-assisted breeding of this important species.

In-depth comparative transcriptome analysis of intestines of red swamp crayfish, Procambarus clarkii, infected with WSSV

Crayfish has become one of the most important farmed aquatic species in China due to its excellent disease resistance against bacteria and viruses. However, the antiviral mechanism of crayfish is still not very clear. In the present study, many high-quality sequence reads from crayfish intestine were obtained using Illumina-based transcriptome sequencing. For the normal group (GN), 44,600,142 high-quality clean reads were randomly assembled to produce 125,394 contigs. For the WSSV-challenged group (GW), 47,790,746 high-quality clean reads were randomly assembled to produce 148,983 contigs. After GO annotation, 39,482 unigenes were annotated into three ontologies: biological processes, cellular components, and molecular functions. In addition, 15,959 unigenes were mapped to 25 different COG categories. Moreover, 7,000 DEGs were screened out after a comparative analysis between the GN and GW samples, which were mapped into 250 KEGG pathways. Among these pathways, 36 were obviously changed (P-values < 0.05) and 28 pathways were extremely significantly changed (P-values < 0.01). Finally, five key DEGs involved in the JAK-STAT signaling pathway were chosen for qRT-PCR. The results showed that these five DEGs were obviously up-regulated at 36 h post WSSV infection in crayfish intestine. These results provide new insight into crayfish antiviral immunity mechanisms.