Sequencing and de novo analysis of the hemocytes transcriptome in Litopenaeus vannamei response to white spot syndrome virus infection

Correlative analysis of transcriptome and proteome in Penaeus vannamei reveals key signaling pathways are involved in IFN-like antiviral regulation mediated by interferon regulatory factor (PvIRF)

Interferon regulatory factors (IRFs) are crucial transcription factors that regulate interferon (IFN) induction in response to pathogen invasion. The regulatory mechanism of IRF has been well studied in vertebrates, but little has been known in arthropods. Therefore, in order to obtain new insights into the potential molecular mechanism of Peneaus vannamei IRF (PvIRF) in response to viral infection, comprehensive comparative analysis of the transcriptome and proteome profiles in shrimp infected with WSSV after knocking down PvIRF was conducted by using RNA sequencing (RNA-seq) and isobaric tags for relative and absolute quantification (iTRAQ). The sequence characterization, molecular functional evolution and 3D spatial structure of PvIRF were analyzed by using bioinformatics methods. PvIRF share the higher homology with different species in N-terminal end (containing DNA binding domain (DBD) including DNA sequence recognition sites and metal binding site) than that in C-terminal end. Within 4 IRF subfamilies of vertebrates, PvIRF had closer relationship with IRF1 subfamily. The DBD of PvIRF and C. gigas IRF1a were composed of α-helices and β-folds which was similar with the DBD structure of M. musculus IRF2. Interestingly, different from the five Tryptophan repeats highly homologous in the DBD of vertebrate IRF, the first and fifth tryptophans of PvIRF mutate to Phenylalanine and Leucine respectively, while the mutations were conserved among shrimp IRFs. RNAi knockdown of PvIRF gene by double-strand RNA could obviously promote the in vivo propagation of WSSV in shrimp and increase the mortality of WSSV-infected shrimp. It suggested that PvIRF was involved in inhibiting the replication of WSSV in shrimp. A total of 8787 transcripts and 2846 proteins were identified with significantly differential abundances in WSSV-infected shrimp after PvIRF knockdown, among which several immune-related members were identified and categorized into 10 groups according to their possible functions. Furthermore, the variation of expression profile from members of key signaling pathways involving JAK/STAT and Toll signaling pathway implied that they might participate IRF-mediated IFN-like regulation in shrimp. Correlative analyses indicated that 722 differentially expressed proteins (DEPs) shared the same expression profiles with their corresponding transcripts, including recognition-related proteins (CTLs and ITGs), chitin-binding proteins (peritrophin), and effectors (ALFs and SWD), while 401 DEPs with the opposite expression profiles across the two levels emphasized the critical role of post-transcriptional and post-translational modification. The results provide candidate signaling pathway including pivotal genes and proteins involved in the regulatory mechanism of interferon mediated by IRF on shrimp antiviral response. This is the first report in crustacean to explore the IFN-like antiviral regulation pathway mediated by IRF on the basis of transcriptome and proteomics correlative analysis, and will provide new ideas for further research on innate immune and defense mechanisms of crustacean.

Transcriptomic analysis reveals the role of Glycolysis pathway in Litopenaeus vannamei during DIV1 infection

In recent years, shrimp farming has experienced significant losses due to the emergence of DIV1 (Decapod iridescent virus 1), an infectious virus with a high fatality rate among shrimp. In this study, we conducted transcriptomic analyses on shrimp Litopenaeus vannamei hemocytes following DIV1 infection and focused on the function of genes in the Glycolysis pathway during DIV1 infection. A total of 2197 differentially expressed genes (DEGs) were identified, comprising 1506 up-regulated genes and 691 down-regulated genes. These genes were primarily associated with Phagosome, ECM-Receptor Interaction, Drug Metabolism-Other Enzymes, and the AGE-RAGE signaling pathway in diabetic complications. KEGG pathway enrichment analysis of the DEGs revealed a noteworthy correlation with metabolic pathways, with a specific focus on glucose metabolism. Specifically, the Glycolysis/Gluconeogenesis pathway exhibited significant upregulation following DIV1 infection. In line with this, we observed an augmented accumulation of glycolytic-related metabolites in the hemolymph following DIV1 challenge along with upregulation of the relative mRNA expression of several glycolytic-related genes. Moreover, we found that the inhibition of lactate dehydrogenase (LDH) activity through RNAi or the use of an inhibitor resulted in reduced lactate production, effectively safeguarding shrimp from DIV1 infection. These findings not only provide a comprehensive dataset for further investigation into DIV1 pathogenesis but also offer valuable insights into the immunometabolism mechanisms that govern shrimp responses to DIV1 infection.

The pyruvate kinase of the whiteleg shrimp Litopenaeus vannamei: Gene structure and responses to short term hypoxia

The shrimp Litopenaeus vannamei is the main farmed crustaceans worldwide. This crustacean suffers environmental changes in oxygen availability that affect its energy metabolism. Pyruvate kinase (PK) catalyzes the last reaction of glycolysis and is key for the regulation of glycolysis and gluconeogenesis. There is ample knowledge about mammalian PK, but in crustaceans, the information is very scarce. In this study, we analyzed in silico the structures of the PK gene and protein. Also, the effects of hypoxia on gene expression, enzymatic activity, glucose, and lactate in hepatopancreas and muscle were analyzed. The PK gene is 15,103 bp and contains 11 exons and 10 introns, producing four mRNA variants by alternative splicing and named PK1, PK2, PK3 and PK4, and two proteins with longer C-terminus and two with a 12 bp insertion. The promoter contains putative binding sites for transcription factors (TF) that are typically involved in stress responses. The deduced amino acid sequences contain the classic domains, binding sites for allosteric effectors and potential reversible phosphorylation residues. Protein modeling indicates a homotetramer with highly conserved structure. The effect of hypoxia for 6 and 12 h showed tissue-specific patterns, with higher expression, enzyme activity and lactate in muscle, but higher glucose in hepatopancreas. Changes in response to hypoxia were detected at 12 h in expression with induction in muscle and reduction in hepatopancreas, while enzyme activity was maintained, and glucose and lactate decreased. These results show rapid changes in expression and metabolites, while enzyme activity was maintained to cope with short-term hypoxia.

scRNA-seq Analysis of Hemocytes of Penaeid Shrimp Under Virus Infection

The classification of cells in non-model organisms has lagged behind the classification of cells in model organisms that have established cluster of differentiation marker sets. To reduce fish diseases, research is needed to better understand immune-related cells, or hemocytes, in non-model organisms like shrimp and other marine invertebrates. In this study, we used Drop-seq to examine how virus infection affected the populations of hemocytes in kuruma shrimp, Penaeus japonicus, which had been artificially infected with a virus. The findings demonstrated that virus infection reduced particular cell populations in circulating hemolymph and inhibited the expression of antimicrobial peptides. We also identified the gene sets that are likely to be responsible for this reduction. Additionally, we identified functionally unknown genes as novel antimicrobial peptides, and we supported this assumption by the fact that these genes were expressed in the population of hemocytes that expressed other antimicrobial peptides. In addition, we aimed to improve the operability of the experiment by conducting Drop-seq with fixed cells as a source and discussed the impact of methanol fixation on Drop-seq data in comparison to previous results obtained without fixation. These results not only deepen our understanding of the immune system of crustaceans but also demonstrate that single-cell analysis can accelerate research on non-model organisms.

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.

Chronic toxicity of shrimp feed added with silver nanoparticles (Argovit-4®) in Litopenaeus vannamei and immune response to white spot syndrome virus infection

In recent years, the application of silver nanoparticles (AgNPs) as antibacterial compounds has been widely used in human and veterinary medicine. In this work, we investigated the effects of AgNPs (Argovit-4®) as feed additives (feed-AgNPs) on shrimp (Litopenaeus vannamei) using three different methods: 1) chronic toxicity after 28 days of feeding, 2) Effects against white spot syndrome virus (WSSV) challenged by oral route, and 3) transcriptional responses of immune-related genes (PAP, ProPO, CTL-3, Crustin, PEN3, and PEN4) following WSSV infection. The results showed that the feed-AgNPs did not interfere with the growth and survival of shrimp. Also, mild lesions in the hepatopancreas were recorded, proportional to the frequency of the feed-AgNP supply. Challenge test versus WSSV showed that feeding every 7 days with feed-AgNPs reduced mortality, reaching a survival rate of 53%, compared to the survival rates observed in groups fed every 4 days, daily and control groups of feed-AgNPs for the 30%, 10%, and 7% groups, respectively. Feed-AgNPs negatively regulated the expression of PAP, ProPO, and Crustin genes after 28 days of treatment and altered the transcriptional responses of PAP, ProPO, CTL-3, and Crustin after WSSV exposure. The results showed that weekly feeding-AgNPs could partially prevent WSSV infection in shrimp culture. However, whether or not transcriptional responses against pathogens are advantageous remains to be elucidated.

LncRNA-mRNA integrated profiling analysis in response to white spot syndrome virus in hepatopancreas in Penaeus japonicus

Penaeus japonicas is an important shrimp species, which is exposed to stressors including a variety of epidemic diseases. To date, little is known about the mechanisms involved in the response to white spot syndrome virus (WSSV) mediated by long non-coding RNAs (lncRNAs). A total of 6544 putative lncRNAs were identified in the hepatopancreas in P. japonicas, which provides a useful lncRNA reference resource for use in future studies. In addition, a total of 444 differentially expressed mRNAs and 457 differentially expressed lncRNAs were identified at 6, 12, and 24 h after WSSV infection in the hepatopancreas of P. japonicas. Functional enrichment analysis showed that the differentially expressed mRNAs were enriched in terms related to immune response and viral infectivity such as defense response, aminopeptidase activity, whereas the differentially expressed lncRNA partner genes were enriched in ubiquitin-dependent protein catabolic process, lipoprotein metabolic process, and antigen processing and presentation. Moreover, several lncRNAs were induced by WSSV infection, indicating these lncRNAs might participate in regulating many immune processes referring to their partner genes. Co-expression analysis of the lncRNAs and their partner genes identified some high lncRNA-mRNA correlations. These results suggest that WSSV stimulates the immune response in the hepatopancreas potentially through an important coding and non-coding gene network, thereby providing valuable information regarding non-coding responses to WSSV in Penaeus species.

Lipid metabolism, immune and apoptosis transcriptomic responses of the hepatopancreas of Chinese mitten crab to the exposure to microcystin-LR

Global warming is favouring the incidence, intensity and duration of harmful cyanobacterial blooms. Microcystin-LR (MC-LR), a hepatotoxic agent, is produced during cyanobacterial blooms. To understand the molecular mechanisms of acute hepatotoxic effect of low doses of MC-LR in crab, we examined differentially expressed genes in samples of the hepatopancreas of Chinese mitten crab (Eriocheir sinensis) collected in 48 h after injections of MC-LR at doses of 0, 25, 50, and 75 µg/kg. The results revealed that MC-LR induced changes in corresponding gene led to the accumulation of triglycerides. MC-LR exposure affected sterol metabolism. Apoptosis-related genes such as Fas-L, Bcl-XL, Cytc, AiF, p53, PERK, calpain, CASP2, CASP7, α-tubulin, PARP, GF, G12, and PKC were upregulated. Conversely, expression levels of CASP10 and ASK1 were downregulated. Genes related to the regulation of actin cytoskeleton (Rho, ROCK, MLCP, MLC, PAK, and PFN) were upregulated. Further, expression levels of genes encoding fatty acid elongation-related enzymes were upregulated, but the expression of genes related to fatty acid synthesis was slightly down regulated. Taken together, these results demonstrated the hepatic toxicity and molecular mechanisms of changes in lipid metabolism, immune and apoptosis in Chinese mitten crab under the MC-LR-induced stress, which is the first report on crabs and performs a comprehensive analysis and a new insight of the molecular toxicological responses in crabs.

Protein phosphorylation in hemocytes of Fenneropenaeus chinensis in response to white spot syndrome virus infection

Protein phosphorylation and dephosphorylation are the most common and important regulatory mechanisms in signal transduction, which play a vital role in immune defense response. Our previous study has found the level of tyrosine phosphorylation was significantly changed in the hemocytes of Fenneropenaeus chinensis upon white spot syndrome virus (WSSV) infection. In order to explore the relationship between protein phosphorylation and WSSV infection, the quantitative phosphoproteomics was employed to identify differential phosphorylated proteins in hemocytes of F. chinensis before and after WSSV infection, and elucidate the role of key differential phosphorylated proteins in WSSV infection process. The results showed that a total of 147 differential phosphorylated proteins were identified in the hemocytes, including 64 phosphorylated proteins and 83 dephosphorylated proteins, which were mostly enriched in pyruvate metabolism, TCA cycle, glycolysis, and ribosomal biosynthesis. Functional analysis of differential phosphorylated proteins showed that they were involved in cell apoptosis, cell phagocytosis, cell metabolism and antiviral infection. A total of 236 differential phosphorylation sites were found, including 91 modified sites in the phosphorylation proteins and 145 modified sites in the dephosphorylation proteins. Motif analysis showed that these phosphorylation sites could activate mitogen-activated protein kinase, P70 S6 kinase and other kinases in hemocytes. Moveover, the phosphorylation levels of eukaryotic protein initiation factor 4E binding proteins and histone H3 were further determined by ELISA and Western blotting, which both exhibited a significant increase post WSSV infection and reach their peak levels at 6 and 12 h, respectively. Moreover, we found that lactate, a metabolite closely related to pyruvate metabolism, TCA cycle and glycolysis, was significantly increased in the hemocytes after WSSV infection. This study revealed the protein phosphorylation response in hemocytes of F. chinensis to WSSV infection, which help to clarify the response characteristics and virus resistance mechanism of hemocytes in F. chinensis, and also facilitate further understanding of the interaction between WSSV and shrimp hemocytes.

Transcriptome profiling reveals the novel immunometabolism-related genes against WSSV infection from Fenneropenaeus merguiensis

The white spot syndrome virus (WSSV) has been considered a serious threat to shrimp aquaculture. Besides, the activation of cell metabolism as an immune reaction to the virus is now recognized as a piece of the pivotal puzzle of the antiviral responses. Hence, this study explores the relationship between metabolic gene expression and antiviral responses in shrimp using transcriptome analysis. The RNA-seq libraries of Fenneropenaeus merguensis hemocytes after WSSV challenge at early (6 hpi) and late (24 hpi) stages of infection were analyzed to identify differentially expressed genes (DEGs) that the WSSV subverted the expression. One-hundred-thirty-three DEGs that were expressed in response to WSSV infection at both stages were identified. Based on the GO annotation, they were related to innate immunity and metabolic pathway. The expression correlation between "full term" (NGS) and qRT-PCR of 16 representative DEGs is shown. Noticeably, the expression profiles of all the selected metabolic genes involved in glucose metabolism, lipid metabolism, amino acid metabolism, and nucleotide metabolism showed a specific correlation between NGS and qRT-PCR upon WSSV infection. Of these, we further characterized the function related to the WSSV response of glutamine: fructose-6-phosphate aminotransferase (FmGFAT), the rate-limiting enzyme of the hexosamine biosynthesis pathway, which was found to be up-regulated at the late stage of WSSV infection. Suppression of FmGFAT by RNA interference resulted in postponing the death of WSSV-infected shrimp and reduction of viral copy number. These results suggested that the FmGFAT is linked between metabolic change and WSSV responses in shrimp, where the virus-induced metabolic rewiring hijack biological compounds and/or energy sources to benefit the viral replication process.

Comparative transcriptome analysis of Chinese grass shrimp (Palaemonetes sinensis) hepatopancreas under ectoparasitic isopod (Tachaea chinensis) infection

Tachaea chinensis, a parasitic isopod, negatively affects the production of several commercially important shrimp species. To better understand the interaction between shrimp immunity and isopod infection, we performed a transcriptome analysis of the hepatopancreas of Palaemonetes sinensis challenged with T. chinensis. After assembly and annotation, 75,980 high-quality unigenes were obtained using RNA-seq data. Differential gene expression analysis revealed 896 significantly differently expressed genes (DEGs) after infection, with 452 and 444 upregulated and downregulated genes, respectively. Specifically, expression levels of genes involved in detoxification, such as the interferon regulatory factor, venom carboxylesterase-6, serine proteinase inhibitor, and cytochrome P450, were upregulated. Furthermore, expression levels of genes corresponding to retinol dehydrogenase, triosephosphate isomerase, variant ionotropic glutamate receptor, and phosphoenolpyruvate carboxykinase were significantly upregulated after isopod parasitization, indicating that the shrimp's visual system was influenced by isopod parasitization. Moreover, quantitative real-time PCR of 10 DEGs helped validate the RNA-seq findings. These results provide a valuable basis for future studies on the elucidation of immune responses of P. sinensis to T. chinensis infection.

Global mRNA and miRNA Analysis Reveal Key Processes in the Initial Response to Infection with WSSV in the Pacific Whiteleg Shrimp

White Spot Disease (WSD) presents a major barrier to penaeid shrimp production. Mechanisms underlying White Spot Syndrome Virus (WSSV) susceptibility in penaeids are poorly understood due to limited information related to early infection. We investigated mRNA and miRNA transcription in Penaeus vannamei over 36 h following infection. Over this time course, 6192 transcripts and 27 miRNAs were differentially expressed—with limited differential expression from 3–12 h post injection (hpi) and a more significant transcriptional response associated with the onset of disease symptoms (24 hpi). During early infection, regulated processes included cytoskeletal remodelling and alterations in phagocytic activity that may assist WSSV entry and translocation, novel miRNA-induced metabolic shifts, and the downregulation of ATP-dependent proton transporter subunits that may impair cellular recycling. During later infection, uncoupling of the electron transport chain may drive cellular dysfunction and lead to high mortalities in infected penaeids. We propose that post-transcriptional silencing of the immune priming gene Dscam (downregulated following infections) by a novel shrimp miRNA (Pva-pmiR-78; upregulated) as a potential mechanism preventing future recognition of WSSV that may be suppressed in surviving shrimp. Our findings improve our understanding of WSD pathogenesis in P. vannamei and provide potential avenues for future development of prophylactics and treatments.

Ubiquitination as an Important Host-Immune Response Strategy in Penaeid Shrimp: Inferences From Other Species

Shrimp aquaculture is an essential economic venture globally, but the industry faces numerous challenges, especially pathogenic infections. As invertebrates, shrimp rely mainly on their innate immune system for protection. An increasing number of studies have shown that ubiquitination plays a vital role in the innate immune response to microbial pathogens. As an important form of posttranslational modification (PTM), both hosts and pathogens have exploited ubiquitination and the ubiquitin system as an immune response strategy to outwit the other. This short review brings together recent findings on ubiquitination and how this PTM plays a critical role in immune modulation in penaeid shrimps. Key findings inferred from other species would help guide further studies on ubiquitination as an immune response strategy in shrimp-pathogen interactions.

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.

Research into the hemocyte immune response of Fenneropenaeus merguiensis under decapod iridescent virus 1 (DIV1) challenge using transcriptome analysis

The banana shrimp (Fenneropenaeus merguiensis) is a common cultural species worldwide. With the development of the shrimp farming industry, increasing number of diseases have emerged and cause huge impacts. Decapod iridescent virus 1 (DIV1) is a new virus of the family Iridoviridae isolated in China that causes very high mortality in shrimp. In this study, DIV1 and PBS were injected into two groups of shrimp, and hemocytes were collected for comparative transcriptomic analysis. We confirmed that F. merguiensis was the new host of DIV1 by nested PCR. A total of 100,759 unigenes were assembled from the control group and the DIV1 infected group, with an average length of 733.06 bp and N50 of 1136 bp. Significant hits were found in 21,465 unigenes compared to known sequences in major databases including COG (33.30%), GO (42.17%), KEGG (46.76%), KOG (61.37%), Pfam (66.90%), Swissprot (54.21%) and Nr (93.86%). A total of 1003 differentially expressed genes (DEGs) were identified, including 929 up-regulated genes and 74 down-regulated genes. Several known immune-related genes, including caspase, C-type lectin, Wnt5 and integrin, were among the differentially expressed transcripts. A total of 14,459 simple sequence repeats, including 8128 monomers, 3276 dimers, 1693 trimers, 150 quadmers, 4 pentamers and 16 hexamers, were found in the transcriptomic dataset. Our study is the first comprehensive investigation of the transcriptomic response to DIV1 infection in F. merguiensis. Collectively, these results not only provide valuable information for characterizing the immune mechanisms of the shrimp responses to DIV1 infection, they open new ways for the study of the molecular mechanisms of DIV1 infection in F. merguiensis.

Comparative Transcriptome Analysis of Litopenaeus vannamei Reveals That Triosephosphate Isomerase-Like Genes Play an Important Role During Decapod Iridescent Virus 1 Infection

Decapod iridescent virus 1 (DIV1) results in severe economic losses in shrimp aquaculture. However, little is known about the physiological effect of DIV1 infection on the host. In this study, we found that the lethal dose 50 of DIV1-infected Litopenaeus vannamei after 48, 72, 96, and 156 h were 4.86 × 10⁶, 5.07 × 10⁵, 2.13 × 10⁵, and 2.38 × 10⁴ copies/μg DNA, respectively. In order to investigate the mechanisms of DIV1 infection, a comparative transcriptome analysis of hemocytes from L. vannamei, infected or not with DIV1, was conducted. The BUSCO analysis showed that the transcriptome was with high completeness (complete single-copy BUSCOs: 57.3%, complete duplicated BUSCOs: 41.1%, fragmentation: 0.8%, missing: 0.8%). A total of 168,854 unigenes were assembled, with an average length of 601 bp. Based on homology searches, Kyoto Encyclopedia of Genes and Genomes (KEGG), gene ontology (GO), and cluster of orthologous groups of proteins (KOG) analysis, 62,270 (36.88%) unigenes were annotated. Among them, 1,112 differentially expressed genes (DEGs) were identified, of which 889 genes were up-regulated and 223 genes were down-regulated after DIV1 infection. These genes were mainly annotated to the major metabolic processes such as fructose and mannose metabolism, carbon metabolism, and inositol phosphate metabolism. Among these metabolic pathways, the triosephosphate isomerase (TPI) family was the most eye-catching DEG as it participates in several metabolic processes. Three types of TPI, LvTPI-like, LvTPI-Blike, and LvTPI-Blike1, were obtained for gene silencing by RNA interference. The results showed that LvTPI-like and LvTPI-Blike1 silencing caused a high mortality rate among L. vannamei. However, LvTPI-like and LvTPI-Blike silencing reduced DIV1 replication in DIV1-infected L. vannamei. All the results indicated that TPI-like genes play an important role during DIV1 infection, which provides valuable insight into the infection mechanism of DIV1 in shrimp and may aid in preventing viral diseases in shrimp culture.

LuxS/AI-2 Quorum Sensing System in Edwardsiella piscicida Promotes Biofilm Formation and Pathogenicity

LuxS/AI-2 is an important quorum sensing system which affects the growth, biofilm formation, virulence, and metabolism of bacteria. LuxS is encoded by the luxS gene, but how this gene is associated with a diverse array of physiological activities in Edwardsiella piscicida (E. piscicida) is not known. Here, we constructed an luxS gene mutant strain, the △luxS strain, to identify how LuxS/AI-2 affects pathogenicity. The results showed that LuxS was not found in the luxS gene mutant strain, and this gene deletion decreased E. piscicida growth compared to that of the wild-type strain. Meanwhile, the wild-type strain significantly increased penetration and motility in mucin compared to levels with the △luxS strain. The 50% lethal dose (LD₅₀) of the E. piscicida △luxS strain for zebrafish was significantly higher than that of the wild-type strain, which suggested that the luxS gene deletion could attenuate the strain's virulence. The AI-2 activities of EIB202 were 56-fold higher than those in the △luxS strain, suggesting that the luxS gene promotes AI-2 production. Transcriptome results demonstrated that between cells infected with the △luxS strain and those infected with the wild-type strain 46 genes were significantly differentially regulated, which included 34 upregulated genes and 12 downregulated genes. Among these genes, the largest number were closely related to cell immunity and signaling systems. In addition, the biofilm formation ability of EIB202 was significantly higher than that of the △luxS strain. The supernatant of EIB202 increased the biofilm formation ability of the △luxS strain, which suggested that the luxS gene and its product LuxS enhanced biofilm formation in E. piscicida All results indicate that the LuxS/AI-2 quorum sensing system in E. piscicida promotes its pathogenicity through increasing a diverse array of physiological activities.

A new insight to biomarkers related to resistance in survived-white spot syndrome virus challenged giant tiger shrimp, Penaeus monodon

The emergence of diseases such as white spot disease has become a threat to Penaeus monodon cultivation. Although there have been a few studies utilizing RNA-Seq, the cellular processes of host-virus interaction in this species remain mostly anonymous. In the present study, P. monodon was challenged with WSSV by intramuscular injection and survived for 12 days. The effect of the host gene expression by WSSV infection in the haemocytes, hepatopancreas and muscle of P. monodon was studied using Illumina HiSeq 2000. The RNA-Seq of cDNA libraries was developed from surviving WSSV-challenged shrimp as well as from normal healthy shrimp as control. A comparison of the transcriptome data of the two groups showed 2,644 host genes to be significantly up-regulated and 2,194 genes significantly down-regulated as a result of the infection with WSSV. Among the differentially expressed genes, our study discovered HMGB, TNFSF and c-Jun in P. monodon as new potential candidate genes for further investigation for the development of potential disease resistance markers. Our study also provided significant data on the differential expression of genes in the survived WSSV infected P. monodon that will help to improve understanding of host-virus interactions in this species.

Comparative transcriptomic analysis of surf clams (Paphia undulate) infected with two strains of Vibrio spp. reveals the identity of key immune genes involved in host defense

BACKGROUND

Vibrio spp. is the major infection-producing marine bacteria in commercially important bivalve Paphia undulata. The host resistance is the major determining factor for the development of pathogenesis. To explore defense mechanisms, researchers have focused primarily on the study of differential expression of individual or specific groups of host immune genes during pathogen-challenge.

RESULTS

We compared the expression profile in the surf clams infected with avirulent V. alginolyticus and virulent V. parahaemolyticus to mark the possible molecular mechanisms of pathogenesis. Comparison of the differentially expressed genes between the two groups of Vibrio-infected clams revealed that the number of down-regulate genes in V. parahaemolyticus injected clams (1433) were significantly higher than the other group (169). Based on Gene Ontology classification, a large proportion of these down-regulate genes were found to be associated with cellular and molecular mechanisms for pathogen recognition, and immunity development thereby explaining the low survival rate for the V. parahaemolyticus-treated clams and suggesting a higher virulence of this bacterium towards the surf clams. Quantitative real-time PCR of 24 candidate genes related to immunity involving the JAK-STAT signaling pathway, complementary cascade, cytokine signaling pathway, oxidative stress, phagocytosis and apoptosis down regulated under V. parahaemolyticus infection, indicating compromised host defense. Furthermore, we could demonstrate a central role of JAK-STAT pathway in bacterial clearance. dsRNA mediated depletion of a clam STAT homolog gene results in dramatic increase in the infection by V. alginolyticus, a mildly pathogenic strain under control conditions.

CONCLUSIONS

The difference in gene expression profiles in surf clams treated with two Vibrio species with a differential pathogenicity to P. undulate and downstream molecular analysis could enlighten on the probable molecular mechanisms of the Vibrio pathogenesis and the virulence of V. parahaemolyticus in surf clams, which also benefits to develop new strategies for disease control in surf calm aquaculture.

A novel complement C3 like gene (Lv-C3L) from Litopenaeus vannamei with bacteriolytic and hemolytic activities and its role in antiviral immune response

As a core component of the complement system, complement component 3 (C3) plays a central role in the opsonization of pathogens, immune defense and immune regulation in the mammalian for its activation is required to trigger classical as well as alternative complement pathways. However, the molecular mechanism underlying C3 activation in invertebrates remains unknown. Several C3 genes have been characterized in invertebrates but very few in crustacean. To understand the molecular characterization and immunological functions of shrimp C3, we characterized a novel complement C3 like gene (designated Lv-C3L) with full-length cDNA sequence identified from pacific white shrimp Litopenaeus vannamei in the present study. The full length cDNA of Lv-C3L sequence was 4769 bp (GenBank accession number: MH638255) containing a 4077 bp open reading frame (ORF), which encodes 1358 amino acids contained a putative signal peptide of 17 amino acids. Six model motifs of C3 were found in Lv-C3L including typical A2M domain, a highly conserved thioester region (GCGEQ) and proteolytic cleavage site of ANATO. In addition to typical conservative domains, Lv-C3L also contains a particular GLN-rich region which might be involved in the protein interaction and transcriptional activation. The transcripts of Lv-C3L were mainly detected in hemocytes and gill which might be involved in defense response. At 36 h post V.parahaemolyticus and B.thuringensis infection, the expression level of Lv-C3L gene in hemocytes were significantly upregulated. At 48 h and 72 h post WSSV infection, the expression level of Lv-C3L gene in hemocytes and gill were significantly upregulated. These results indicated that Lv-C3L gene play a pivotal role in innate immune responses to the WSSV and G⁺/G^- bacterial infection. The obvious immune function of Lv-C3L was described as an effective membrane rupture in bacteriolytic and hemolytic activities on V.parahaemolyticus, V.anguillarum and rabbit erythrocytes. Combining with WSSV copy number, WSSV-VP28 gene expression profile and shrimp cumulative mortality analysis, RNAi knockdown of Lv-C3L gene could obviously promote the in vivo propagation of WSSV in shrimp. This is the first report in crustaceans that Lv-C3L, as a key complement like components, is involved in shrimp antiviral immune response. It is speculated that complicated complement response cascade may exist in shrimp. These results collectively indicated that the complement pathway in shrimp might play an important protective role against pathogenic infection and activation of complement pathway including C3 could restrict the propagation of WSSV.

Transcriptome analysis reveals the activation of neuroendocrine-immune system in shrimp hemocytes at the early stage of WSSV infection

Background

Functional communications between nervous, endocrine and immune systems are well established in both vertebrates and invertebrates. Circulating hemocytes act as fundamental players in this crosstalk, whose functions are conserved during the evolution of the main groups of metazoans. However, the roles of the neuroendocrine-immune (NEI) system in shrimp hemocytes during pathogen infection remain largely unknown.

Results

In this study, we sequenced six cDNA libraries prepared with hemocytes from Litopenaeus vannamei which were injected by WSSV (white spot syndrome virus) or PBS for 6 h using Illumina Hiseq 4000 platform. As a result, 3444 differentially expressed genes (DEGs), including 3240 up-regulated genes and 204 down-regulated genes, were identified from hemocytes after WSSV infection. Among these genes, 349 DEGs were correlated with innate immunity and categorized into seven groups based on their predictive function. Interestingly, 18 genes encoded putative neuropeptide precursors were induced significantly by WSSV infection. Furthermore, some genes were mapped to several typical processes in the NEI system, including proteolytic processing of prohormones, amino acid neurotransmitter pathways, biogenic amine biosynthesis and acetylcholine signaling pathway.

Conclusions

The data suggested that WSSV infection triggers the activation of NEI in shrimp, which throws a light on the pivotal roles of NEI system mediated by hemocytes in shrimp antiviral immunity.

Electronic supplementary material

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

Transcriptome of white shrimp Litopenaeus vannamei induced with rapamycin reveals the role of autophagy in shrimp immunity

Autophagy plays a vital role in innate and adaptive immunity against invading microorganisms, such as virus and bacteria. However, the mechanism underlying autophagy in shrimp is still limited. In our study, we challenged white shrimp L. vannamei with rapamycin to induce autophagy and employed Solexa/Illumina high-throughput RNA-seq method to examine the differences of transcriptome from gills of shrimps treated with or without rapamycin. More than 22.64 Gb raw data were produced, which were assembled into 62, 503 unigenes, with 14,126 unigenes over 1 kb in length. We then performed differential expression analysis and identified a total of 3050 differentially expressed genes (DEGs). Among them, 1456 were upregulated and 1594 were downregulated. We further annotated DEGs by matching against non-redundant protein sequence (Nr), Swiss-Prot, Kyoto Encyclopedia of Genes and Genomes (KEGG), Clusters of Orthologous Groups of proteins (COG), euKaryotic Orthologous Groups (KOG), Gene ontology (GO), and Pfam databases. The assembled and annotated DEGs will facilitate our understanding of the molecular mechanism underlying autophagy and promote the studies on the role of autophagy in innate immunity of L. vannamei and other crustaceans.

The Protein-Protein Interaction Network of Litopenaeus vannamei Haemocytes

Protein–protein interaction networks (PINs) have been constructed in various organisms and utilized to conduct evolutionary analyses and functional predictions. Litopenaeus vannamei is a high-valued commercial aquaculture species with an uncharacterized interactome. With the development of RNA-seq techniques and systems biology, it is possible to obtain genome-wide transcriptional information for L. vannamei and construct a systematic network based on these data. In this work, based on the RNA-seq of haemocytes we constructed the first L. vannamei PIN including 4,858 proteins and 104,187 interactions. The PIN constructed here is the first large-scale PIN for shrimp. The confidence scores of interactions in the PIN were evaluated on the basis of sequence homology and genetic relationships. The immune-specific sub-network was extracted from global PIN, and more than a third of proteins were found in signaling pathways in the sub-network, which indicates an inseparable relationship between signaling processes and immune mechanisms. Six selected signaling pathways were constructed at different age groups based on evolutionary analyses. Furthermore, we showed that the functions of the pathways’ proteins were associated with their evolutionary history based on the evolutionary analyses combining with protein functional analyses. In addition, the functions of 1,955 unclassified proteins which were associated with 3,191 unigenes were assigned using the PIN, which account for approximately 70.3 and 44.9% of the previously unclassified proteins and unigenes in the network, respectively. The annotation of unclassified proteins and unigenes based on the PIN provides new candidates for further functional studies. The immune-specific sub-network and the pathways extracted from the PIN provide a novel information source for studying of immune mechanisms and disease resistances in shrimp.

Effects of environmental stress on shrimp innate immunity and white spot syndrome virus infection

The shrimp aquaculture industry is plagued by disease. Due to the lack of deep understanding of the relationship between innate immune mechanism and environmental adaptation mechanism, it is difficult to prevent and control the diseases of shrimp. The shrimp innate immune system has received much recent attention, and the functions of the humoral immune response and the cellular immune response have been preliminarily characterized. The role of environmental stress in shrimp disease has also been investigated recently, attempting to clarify the interactions among the innate immune response, the environmental stress response, and disease. Both the innate immune response and the environmental stress response have a complex relationship with shrimp diseases. Although these systems are important safeguards, allowing shrimp to adapt to adverse environments and resist infection, some pathogens, such as white spot syndrome virus, hijack these host systems. As shrimp lack an adaptive immune system, immunization therapy cannot be used to prevent and control shrimp disease. However, shrimp diseases can be controlled using ecological techniques. These techniques, which are based on the innate immune response and the environmental stress response, significantly reduce the impact of shrimp diseases. The object of this review is to summarize the recent research on shrimp environmental adaptation mechanisms, innate immune response mechanisms, and the relationship between these systems. We also suggest some directions for future research.

Single-molecule long-read sequencing facilitates shrimp transcriptome research

Although shrimp are of great economic importance, few full-length shrimp transcriptomes are available. Here, we used Pacific Biosciences single-molecule real-time (SMRT) long-read sequencing technology to generate transcripts from the Pacific white shrimp (Litopenaeus vannamei). We obtained 322,600 full-length non-chimeric reads, from which we generated 51,367 high-quality unique full-length transcripts. We corrected errors in the SMRT sequences by comparison with Illumina-produced short reads. We successfully annotated 81.72% of all unique SMRT transcripts against the NCBI non-redundant database, 58.63% against Swiss-Prot, 45.38% against Gene Ontology, 32.57% against Clusters of Orthologous Groups of proteins (COG), and 47.83% against Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. Across all transcripts, we identified 3,958 long non-coding RNAs (lncRNAs) and 80,650 simple sequence repeats (SSRs). Our study provides a rich set of full-length cDNA sequences for L. vannamei, which will greatly facilitate shrimp transcriptome research.

RNAi screening identifies a new Toll from shrimp Litopenaeus vannamei that restricts WSSV infection through activating Dorsal to induce antimicrobial peptides

The function of Toll pathway defense against bacterial infection has been well established in shrimp, however how this pathway responds to viral infection is still largely unknown. In this study, we report the Toll4-Dorsal-AMPs cascade restricts the white spot syndrome virus (WSSV) infection of shrimp. A total of nine Tolls from Litopenaeus vannamei namely Toll1-9 are identified, and RNAi screening in vivo reveals the Toll4 is important for shrimp to oppose WSSV infection. Knockdown of Toll4 results in elevated viral loads and renders shrimp more susceptible to WSSV. Furthermore, Toll4 could be a one of upstream pattern recognition receptor (PRR) to detect WSSV, and thereby leading to nuclear translocation and phosphorylation of Dorsal, the known NF-κB transcription factor of the canonical Toll pathway. More importantly, silencing of Toll4 and Dorsal contributes to impaired expression of a specific set of antimicrobial peptides (AMPs) such as anti-LPS-factor (ALF) and lysozyme (LYZ) family, which exert potent anti-WSSV activity. Two AMPs of ALF1 and LYZ1 as representatives are demonstrated to have the ability to interact with several WSSV structural proteins to inhibit viral infection. Taken together, we therefore identify that the Toll4-Dorsal pathway mediates strong resistance to WSSV infection by inducing some specific AMPs.

The TLR pathway mediated antiviral immune response is well identified in mammals, yet, Toll pathway governing this protection in invertebrates remains unknown. In the present study, we uncover that a shrimp Toll4 from a total of nine Tolls in L. vannamei confers resistance to WSSV thought inducing the NF-κB transcription factor Dorsal to inspire the production of some antimicrobial peptides (AMPs) with antiviral activity. The anti-LPS-factor (ALF) and lysozyme (LYZ) family are identified as the Toll4-Dorsal pathway targeted genes with the ability to interact with viral structural proteins in response to WSSV infection. These results suggest that the Toll receptor induces the expression of AMPs with antiviral activity could be a general antiviral mechanism in invertebrates and Toll pathway established antiviral defense could be conserved during evolution.

De novo Assembly and Transcriptome Characterization of Opisthopappus (Asteraceae) for Population Differentiation and Adaption

Opisthopappus Shih (Asteraceae), an endangered genus endemic to the Taihang Mountains of China, is a high-value ornamental and medicinal plant consisting of two species, Opisthopappus longilobus shih and Opisthopappus taihangensis (Ling) Shih. However, the evolutionary relationships and the taxonomic characteristics between the two species remain unknown. In this study, high-throughput transcriptome sequencing was used to analyze the differential metabolic activity and gene expression and screened special molecular markers for exploring the genetic variation and species differentiation in Opisthopappus Shih. The results showed that 33,974 unigenes with an average size of 801 bp were obtained with optimization of de novo assembly. The comprehensive functional annotation based on Gene Ontology (GO), Cluster of Orthologous Group (COG) and Kyoto Encyclopedia of Genes and Genomes pathway database (KEGG) revealed that these unigenes were mainly related to many physiological, metabolic, and molecular processes. Furthermore, the comparative transcriptome analysis indicated that 3,410 differentially expressed genes were mainly involved in lipid, carbohydrate and amino acid metabolism, xenobiotics biodegradation and metabolism as well as environment adaptation via KEGG. Such as the CYP710A, GST, HSP90A and so on, could be the potential candidate genes for further investigating the molecular mechanism of physiological variations between O. taihangensis and O. longilobus. In addition, the potential 71,804 high quality single nucleotide polymorphisms (SNPs) and 1,444 simple sequence repeats (SSRs) were estimated. Based on the predicted SNP, we have developed eight SNP markers for population genetic analysis in Opisthopappus Shih. A significantly high level of genetic differentiation between the populations of O. longilobus and O. taihangensis were found, and they were clearly grouped into two distinct genetic clusters. These results conformed to the record of Flora Reipublicae Popularis Sinicae (FRPS) and unsupported the taxonomic status in the Flora of China. The transcriptome analysis of Opisthopappus Shih can contribute to in-depth exploring of internal mechanisms in species variation and differentiation based on molecular evidence. With the rich and valuable data resources, the more novel structural, functional, and comparative genomic studies will provide comprehensive insights into the evolutionary relationships between O. taihangensis and O. longilobus.

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.

White spot syndrome virus (WSSV) suppresses penaeidin expression in Marsupenaeus japonicus hemocytes

Penaeidins are a unique family of antimicrobial peptides specific to penaeid shrimp and have been reported mainly function as anti-bacterial and anti-fungal. In order to investigate whether penaeidins could also respond to virus or not, we examined the effect of WSSV on MjPen-II (penaeidin in kuruma shrimp, Marsupenaeus japonicus) expression. In the control group, MjPen-II transcript level can be detected in almost all test tissues but was expressed most strongly in hemocytes. After WSSV infection, MjPen-II transcript level was significantly downregulated in hemocytes. Moreover, the proportion of MjPen-II+ hemocytes was not significantly different between non-infected and WSSV-infected shrimp, but the number of MjPen-II+ highly expressing hemocytes decreased after infection. In addition, MjPen-II was observed in the cytoplasm of granule-containing hemocytes. These results suggest that WSSV suppresses MjPen-II expression in hemocytes.

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.

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.

Deep transcriptome profiling sheds light on key players in nucleus implantation induced immune response in the pearl oyster Pinctada martensii

Immunological rejection of the pearl oysters following nucleus implantation is a major issue limiting the successful rate of cultured pearls. To date, the molecular mechanism of immune tolerance during pearl formation in the pearl oysters is still largely unknown. Through the RNA sequencing platform and comparative transcriptomic analysis, we investigated the chronic gene expression changes at seven time points (0, 5, 10, 15, 20, 30, 60 days post implantation or dpi) over a period of 60 days following nucleus implantation in the pearl oyster Pinctada martensii. A total of 81,390 unique transcripts (or unigenes) with a combined length of 96.8 million bp and a N50 value of 2227 bp were obtained. When compared with sequences in the nr, nt, Swiss-Prot, KEGG, COG and GO databases, 36,380 unigenes can find homologous genes. Pairwise comparison of gene expression among all the samples showed that the largest number (or 6846) of differentially expressed genes was observed at 10 dpi. The number then decreased to below 5000 at 15, 20 and 30 dpi and increased again to 6679 at 60 dpi. PCA analysis further showed that the seven time points can be roughly divided into four groups. Comparative transcriptomic analysis between the four groups identified a variety of genes showing differential expression at different time points, including many immune-related genes such as those encoding for toll-like receptor, lectin, scavenger receptor, and peroxidase. In addition, GO and KEGG enrichment analysis revealed that these differentially expressed genes were mainly associated with metabolism, ribosome function, immune response, signaling transduction, and cytoskeleton organization. Notably, two KEGG pathways, namely "cell adhesion molecules" and "primary immunodeficiency" were significantly enriched during the whole process. This finding indicates that genes in these pathways are likely to play critical roles in the immune tolerance of the pearl oysters. To conclude, the data obtained contribute to a better understanding of the molecular mechanisms of nucleus implantation induced immune response in the pearl oysters, and will facilitate the development of effective measures to improve the performance of pearl culture.

Integrated analysis of mRNA-seq in the haemocytes of Eriocheir sinensis in response to Spiroplasma eriocheiris infection

The Chinese mitten crab Eriocheir sinensis is an important economic crustacean that has been exposed to various diseases. Spiroplasma eriocheiris, isolated from tremor-diseased E. sinensis, was first identified as a lethal pathogen of freshwater crustaceans. To understand the pathogenesis of S. eriocheiris to E. sinensis, the transcriptomic profiles of haemocytes in the experimental and control groups at 1 d and 7 d post-injection were obtained using Illumina HiSeq 2500. These results showed that 40,358,724, 44,462,112, 45,516,576 and 37,713,728 paired-end clean reads were obtained from the cDNA libraries of DZ1 (the control group at 1 d), DZ7 (the control group at 7 d), SY1 (the experimental group at 1 d) and SY7 (the experimental group at 7 d), respectively. In total, 106,641 unique transcript fragments (unigenes) were assembled, with an average length of 710 bp. On the first day of stimulation, 33,084 up-regulated transcripts and 19,208 down-regulated transcripts were found in the experimental group compared with those in the control group. On the seventh day of stimulation, 40,198 up-regulated transcripts and 12,032 down-regulated transcripts were found in the experimental group compared with those in the control group. Some canonical immune-related pathways were identified via KEGG pathway analysis, including complement and coagulation cascades, the VEGF signalling pathway, the Wnt signalling pathway, natural killer cell-mediated cytotoxicity, the MAPK signalling pathway, neuroactive ligand-receptor interactions, and the Lysosome pathway. We found important immune-related genes (GNPTAB, MASP2, F7, F5, NFATC, TRAF6, MAP3K5, and TRa) in the KEGG pathway, and those genes were confirmed by qRT-PCR analysis. In addition, the significantly enriched neuroactive ligand-receptor interaction pathway was associated with intense paroxysmal tremors of infected crabs. Our results provide valuable information for the further analysis of the mechanisms of E. sinensis defence against S. eriocheiris invasion.

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.

Insights into Sexual Precocity of Female Oriental River Prawn Macrobrachium nipponense through Transcriptome Analysis

Background

The oriental river prawn (Macrobrachium nipponense) is the most prevalent aquaculture species in China. The sexual precocity in this species has received considerable attention in recent years because more and more individuals matured at a small size, which devalues the commercial production. In this study, we developed deep-coverage transcriptomic sequencing data for the ovaries of sexually precocious and normal sexually mature M. nipponense using next-generation RNA sequencing technology and attempted to provide the first insight into the molecular regulatory mechanism of sexual precocity in this species.

Results

A total of 63,336 unigenes were produced from the ovarian cDNA libraries of sexually precocious and normal sexually mature M. nipponense using Illumina HiSeq 2500 platform. Through BLASTX searches against the NR, STRING, Pfam, Swissprot and KEGG databases, 15,134 unigenes were annotated, accounting for 23.89% of the total unigenes. 5,195 and 3,227 matched unigenes were categorized by GO and COG analysis respectively. 15,908 unigenes were consequently mapped into 332 KEGG pathways, and many reproduction-related pathways and genes were identified. Moreover, 26,008 SSRs were identified from 18,133 unigenes. 80,529 and 80,516 SNPs were yielded from ovarian libraries of sexually precocious and normal sexually mature prawn, respectively, and 29,851 potential SNPs between these two groups were also predicted. After comparing the ovarian libraries of sexually precocious and normal sexually mature prawn, 549 differentially expressed genes (DEGs) and 9 key DEGs that may be related to sexual precocity of M. nipponense were identified. 20 DEGs were selected for validation by quantitative real-time PCR (QPCR) and 19 DEGs show consistent expression between QPCR and RNAseq-based differential expression analysis datasets.

Conclusion

This is the first report on the large-scale RNA sequencing of ovaries of sexually precocious and normal sexually mature M. nipponense. The annotated transcriptome data will provide fundamental support for future research into the reproduction biology of M. nipponense. The large number of candidate SNPs and SSRs detected in this study could be used as genetic markers for population genetics and functional genomics in this species. More importantly, many DEGs, especially nine key DEGs between sexually precocious and normal sexually mature prawns were identified, which will dramatically improve understanding of molecular regulatory mechanism of sexual precocity of this species.

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.

A transcriptome study on Macrobrachium rosenbergii hepatopancreas experimentally challenged with white spot syndrome virus (WSSV)

The world production of shrimp such as the Malaysian giant freshwater prawn, Macrobrachium rosenbergii is seriously affected by the white spot syndrome virus (WSSV). There is an urgent need to understand the host pathogen interaction between M. rosenbergii and WSSV which will be able to provide a solution in controlling the spread of this infectious disease and lastly save the aquaculture industry. Now, using Next Generation Sequencing (NGS), we will be able to capture the response of the M. rosenbergii to the pathogen and have a better understanding of the host defence mechanism. Two cDNA libraries, one of WSSV-challenged M. rosenbergii and a normal control one, were sequenced using the Illumina HiSeq™ 2000 platform. After de novo assembly and clustering of the unigenes from both libraries, 63,584 standard unigenes were generated with a mean size of 698bp and an N50 of 1137bp. We successfully annotated 35.31% of all unigenes by using BLASTX program (E-value <10-5) against NCBI non-redundant (Nr), Swiss-Prot, Kyoto Encyclopedia of Genes and Genome pathway (KEGG) and Orthologous Groups of proteins (COG) databases. Gene Ontology (GO) assessment was conducted using BLAST2GO software. Differentially expressed genes (DEGs) by using the FPKM method showed 8443 host genes were significantly up-regulated whereas 5973 genes were significantly down-regulated. The differentially expressed immune related genes were grouped into 15 animal immune functions. The present study showed that WSSV infection has a significant impact on the transcriptome profile of M. rosenbergii's hepatopancreas, and further enhanced the knowledge of this host-virus interaction. Furthermore, the high number of transcripts generated in this study will provide a platform for future genomic research on freshwater prawns.

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.

Transcriptome and Molecular Pathway Analysis of the Hepatopancreas in the Pacific White Shrimp Litopenaeus vannamei under Chronic Low-Salinity Stress

The Pacific white shrimp Litopenaeus vannamei is a euryhaline penaeid species that shows ontogenetic adaptations to salinity, with its larvae inhabiting oceanic environments and postlarvae and juveniles inhabiting estuaries and lagoons. Ontogenetic adaptations to salinity manifest in L. vannamei through strong hyper-osmoregulatory and hypo-osmoregulatory patterns and an ability to tolerate extremely low salinity levels. To understand this adaptive mechanism to salinity stress, RNA-seq was used to compare the transcriptomic response of L. vannamei to changes in salinity from 30 (control) to 3 practical salinity units (psu) for 8 weeks. In total, 26,034 genes were obtained from the hepatopancreas tissue of L. vannamei using the Illumina HiSeq 2000 system, and 855 genes showed significant changes in expression under salinity stress. Eighteen top Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were significantly involved in physiological responses, particularly in lipid metabolism, including fatty-acid biosynthesis, arachidonic acid metabolism and glycosphingolipid and glycosaminoglycan metabolism. Lipids or fatty acids can reduce osmotic stress in L. vannamei by providing additional energy or changing the membrane structure to allow osmoregulation in relevant organs, such as the gills. Steroid hormone biosynthesis and the phosphonate and phosphinate metabolism pathways were also involved in the adaptation of L. vannamei to low salinity, and the differential expression patterns of 20 randomly selected genes were validated by quantitative real-time PCR (qPCR). This study is the first report on the long-term adaptive transcriptomic response of L. vannamei to low salinity, and the results will further our understanding of the mechanisms underlying osmoregulation in euryhaline crustaceans.

Transcriptome analysis of the Capra hircus ovary

BACKGROUND

Capra hircus is an important economic livestock animal, and therefore, it is necessary to discover transcriptome information about their reproductive performance. In this study, we performed de novo transcriptome sequencing to produce the first transcriptome dataset for the goat ovary using high-throughput sequencing technologies. The result will contribute to research on goat reproductive performance.

METHOD AND RESULTS

RNA-seq analysis generated more than 38.8 million clean paired end (PE) reads, which were assembled into 80,069 unigenes (mean size = 619 bp). Based on sequence similarity searches, 64,824 (60.6%) genes were identified, among which 29,444 and 11,271 unigenes were assigned to Gene Ontology (GO) categories and Clusters of Orthologous Groups (COG), respectively. Searches in the Kyoto Encyclopedia of Genes and Genomes pathway database (KEGG) showed that 27,766 (63.4%) unigenes were mapped to 258 KEGG pathways. Furthermore, we investigated the transcriptome differences of goat ovaries at two different ages using a tag-based digital gene expression system. We obtained a sequencing depth of over 5.6 million and 5.8 million tags for the two ages and identified a large number of genes associated with reproductive hormones, ovulatory cycle and follicle. Moreover, many antisense transcripts and novel transcripts were found; clusters with similar differential expression patterns, enriched GO terms and metabolic pathways were revealed for the first time with regard to the differentially expressed genes.

CONCLUSIONS

The transcriptome provides invaluable new data for a functional genomic resource and future biological research in Capra hircus, and it is essential for the in-depth study of candidate genes in breeding programs.

RNA-seq analysis of Macrobrachium rosenbergii hepatopancreas in response to Vibrio parahaemolyticus infection

Background

The Malaysian giant freshwater prawn, Macrobrachium rosenbergii, is an economically important crustacean worldwide. However, production of this prawn is facing a serious threat from Vibriosis disease caused by Vibrio species such as Vibrio parahaemolyticus. Unfortunately, the mechanisms involved in the immune response of this species to bacterial infection are not fully understood. We therefore used a high-throughput deep sequencing technology to investigate the transcriptome and comparative expression profiles of the hepatopancreas from this freshwater prawn infected with V. parahaemolyticus to gain an increased understanding of the molecular mechanisms underlying the species’ immune response to this pathogenic bacteria.

Result

A total of 59,122,940 raw reads were obtained from the control group, and 58,385,094 reads from the Vibrio-infected group. Via de novo assembly by Trinity assembler, 59,050 control unigenes and 73,946 Vibrio-infected group unigenes were obtained. By clustering unigenes from both libraries, a total of 64,411 standard unigenes were produced. The standard unigenes were annotated against the NCBI non-redundant, Swiss-Prot, Kyoto Encyclopaedia of Genes and Genome pathway (KEGG) and Orthologous Groups of Proteins (COG) databases, with 19,799 (30.73%), 16,832 (26.13%), 14,706 (22.83%) and 7,856 (12.19%) hits respectively, giving a final total of 22,455 significant hits (34.86% of all unigenes). A Gene Ontology (GO) analysis search using the Blast2GO program resulted in 6,007 unigenes (9.32%) being categorized into 55 functional groups. A differential gene expression analysis produced a total of 14,569 unigenes aberrantly expressed, with 11,446 unigenes significantly up-regulated and 3,103 unigenes significantly down-regulated. The differentially expressed immune genes fall under various processes of the animal immune system.

Conclusion

This study provided an insight into the antibacterial mechanism in M. rosenbergii and the role of differentially expressed immune genes in response to V. parahaemolyticus infection. Furthermore, this study has generated an abundant list of transcript from M.rosenbergii which will provide a fundamental basis for future genomics research in this field.

Electronic supplementary material

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

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.