Transcriptome analysis and discovery of genes involved in immune pathways from hepatopancreas of microbial challenged mitten crab Eriocheir sinensis

A Review on the Involvement of Heat Shock Proteins (Extrinsic Chaperones) in Response to Stress Conditions in Aquatic Organisms

Heat shock proteins (HSPs) encompass both extrinsic chaperones and stress proteins. These proteins, with molecular weights ranging from 14 to 120 kDa, are conserved across all living organisms and are expressed in response to stress. The upregulation of specific genes triggers the synthesis of HSPs, facilitated by the interaction between heat shock factors and gene promoter regions. Notably, HSPs function as chaperones or helper molecules in various cellular processes involving lipids and proteins, and their upregulation is not limited to heat-induced stress but also occurs in response to anoxia, acidosis, hypoxia, toxins, ischemia, protein breakdown, and microbial infection. HSPs play a vital role in regulating protein synthesis in cells. They assist in the folding and assembly of other cellular proteins, primarily through HSP families such as HSP70 and HSP90. Additionally, the process of the folding, translocation, and aggregation of proteins is governed by the dynamic partitioning facilitated by HSPs throughout the cell. Beyond their involvement in protein metabolism, HSPs also exert a significant influence on apoptosis, the immune system, and various characteristics of inflammation. The immunity of aquatic organisms, including shrimp, fish, and shellfish, relies heavily on the development of inflammation, as well as non-specific and specific immune responses to viral and bacterial infections. Recent advancements in aquatic research have demonstrated that the HSP levels in populations of fish, shrimp, and shellfish can be increased through non-traumatic means such as water or oral administration of HSP stimulants, exogenous HSPs, and heat induction. These methods have proven useful in reducing physical stress and trauma, while also facilitating sustainable husbandry practices such as vaccination and transportation, thereby offering health benefits. Hence, the present review discusses the importance of HSPs in different tissues in aquatic organisms (fish, shrimp), and their expression levels during pathogen invasion; this gives new insights into the significance of HSPs in invertebrates.

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.

Transcriptome, intestinal microbiome and histomorphology profiling of differences in the response of Chinese sea bass (Lateolabrax maculatus) to Aeromonas hydrophila infection

The Chinese sea bass (Lateolabrax maculatus) is an important aquaculture fish, but diseases caused by Aeromonas hydrophila have led to severe economic losses to the aquaculture industry in recent years. To date, only a few studies have focused on the relationship between the intestinal immune response and changes in intestinal microbes by A. hydrophila infection. Here, we report the transcriptome and intestinal changes in infected sea bass. Histopathological results showed that severe steatosis and vacuolation occurred in the liver and that the intestinal villi and mesentery were seriously affected after infection. By extracting total RNA from intestinal tissue and studying the transcriptome profile, 1,678 genes (1,013 upregulated and 665 downregulated) were identified as significantly differentially expressed genes (DEGs). These genes are involved in many immune-related signalling pathways, such as the NOD-like receptor, C-type lectin receptor, and Toll-like receptor signalling pathways. Moreover, the intestinal microbes of sea bass changed significantly after infection. Interestingly, at the genus level, there was an increase in Serratia, Candida arthromitus and Faecalibacterium as well as a decrease in Akkermansia and Parabacteroides after infection. The results also indicated that some of the DEGs involved in the immune response were related to the genus level of intestinal microbiota. Finally, there was a relationship between gene expression patterns and the bacterial structure in the host intestine. Our study provides a reference for the study of the immune response and particular functions of intestinal microbes of sea bass after pathogen infection.

Full-length transcriptome sequencing and comparative transcriptome analysis of Eriocheir sinensis in response to infection by the microsporidian Hepatospora eriocheir

As a new generation of high-throughput sequencing technology, PacBio Iso-Seq technology (Iso-Seq) provides a better alternative sequencing method for the acquisition of full-length unigenes. In this study, a total of 22.27 gigabyte (Gb) subread bases and 128,614 non-redundant unigenes (mean length: 2,324 bp) were obtained from six main tissues of Eriocheir sinensis including the heart, nerve, intestine, muscle, gills and hepatopancreas. In addition, 74,732 unigenes were mapped to at least one of the following databases: Non-Redundant Protein Sequence Database (NR), Gene Ontology (GO), Kyoto Encyclopaedia of Genes and Genomes (KEGG), KEGG Orthology (KO) and Protein family (Pfam). In addition, 6696 transcription factors (TFs), 28,458 long non-coding RNAs (lncRNAs) and 94,230 mRNA-miRNA pairs were identified. Hepatospora eriocheir is the primary pathogen of E. sinensis and can cause hepatopancreatic necrosis disease (HPND); the intestine is the main target tissue. Here, we attempted to identify the key genes related to H. eriocheir infection in the intestines of E. sinensis. By combining Iso-Seq and Illumina RNA-seq analysis, we identified a total of 12,708 differentially expressed unigenes (DEUs; 6,696 upregulated and 6,012 downregulated) in the crab intestine following infection with H. eriocheir. Based on the biological analysis of these DEUs, several key processes were identified, including energy metabolism-related pathways, cell apoptosis and innate immune-related pathways. Twelve selected genes from these DEUs were subsequently verified by quantitative real-time PCR (qRT-PCR) analysis. Our findings enhance our understanding of the E. sinensis transcriptome and the specific association between E. sinensis and H. eriocheir infection.

Garlic Powder Supplementation Improves Growth, Nonspecific Immunity, Antioxidant Capacity, and Intestinal Flora of Chinese Mitten Crabs (Eriocheir sinensis)

This study was conducted to survey the effects of garlic powder on growth performance, nonspecific immunity, antioxidant capacity, and intestinal flora structure of Chinese mitten crabs. Altogether, 216 crabs which originally weigh 20.71 ± 0.13 g were randomly allocated into three treatment groups with 6 replicates of 12 crabs per replicate. The control group (CN) was fed a basal diet, while the other two groups were fed the basal diet supplemented with 1000 mg/kg (GP1000) and 2000 mg/kg (GP2000) garlic powder, respectively. This trial lasted 8 weeks. The results showed that the supplementation of garlic powder improved the final body weight, weight gain rate, and specific growth rate of the crabs (P < 0.05). Meanwhile, in serum, better nonspecific immune was confirmed by the enhancement of phenoloxidase and lysozyme levels, with the improvement of phosphatase activities in GP1000 and GP2000 (P < 0.05). On the other hand, the levels of total antioxidant capacity, glutathione peroxidases, and total superoxide dismutase in serum and hepatopancreas were increased (P < 0.05) while malondialdehyde content declined (P < 0.05) as the garlic powder was added to the basal diet. And, catalase in serum also shows an increase (P < 0.05). In both GP1000 and GP2000, genes related to antioxidant and immunity, for instance, Toll-like receptor 1, glutathione peroxidase, catalase, myeloid differentiation factor 88, TuBe, Dif, relish, crustins, antilipopolysaccharide factor, lysozyme, and prophenoloxidase mRNA expression levels, were increased (P < 0.05). The abundance of Rhizobium and Rhodobacter was reduced by adding garlic powder (P < 0.05). This study indicated that dietary addition of garlic powder promoted growth, enhanced nonspecific immunity and antioxidant capacity, activated Toll pathway, IMD pathway, and proPO system, increased antimicrobial peptide expression, while simultaneously improving the intestinal flora of Chinese mitten crabs.

Comparison of immune defense and antioxidant capacity between broodstock and hybrid offspring of juvenile shrimp (Macrobrachium nipponense): Response to acute ammonia stress

Ammonia is a major environmental pollutant in the aquatic system that poses a great threat to the health of shrimp. Macrobrachium nipponense, as one of the large-yield farmed shrimp, is facing germplasm degradation. Genetic improvement through hybridization is one of the effective methods to solve this problem. However, there are few studies on the effects of ammonia nitrogen on the germplasm resources of M. nipponense. In this study, the broodstock populations (Dianshan, DS) and hybrid offspring (DS ♀ × CD [Changjiang, CJ ♂ × Dongting, DT ♀], SCD) were exposed to 0, 5, or 20 mg/L of ammonia for 96 h. The survival rate of the SCD group was greater than the DS group, although there were no significant differences in weight gain rate and length gain rate (p > 0.05). The number of positive cells and apoptosis rates in the DS group were significantly greater than in the SCD group after ammonia exposure (p < 0.05). As the ammonia concentration increased, the antioxidant enzyme activities in the SCD group were significantly higher than DS group, while the hepatotoxicity enzyme activities in the SCD group were significantly lower than DS group (p < 0.05). The trends in the expression of antioxidant- and immune-related genes were generally consistent with the activities of antioxidant enzymes. Our study found that the hybrid population had stronger stress resistance than their parent populations at the same ammonia concentration. This study confirms our speculation that hybrid population has a greater advantage in antioxidant immunity, which also provides reference for the follow-up study of chronic ammonia toxicity.

Insights into the species evolution of Calanus copepods in the northern seas revealed by de novo transcriptome sequencing

Copepods of the zooplankton genus Calanus play a key role in marine ecosystems in the northern seas. Although being among the most studied organisms on Earth, due to their ecological importance, genomic resources for Calanus spp. remain scarce, mostly due to their large genome size (from 6 to 12 Gbps). As an alternative to whole-genome sequencing in Calanus spp., we sequenced and de novo assembled transcriptomes of five Calanus species: Calanus glacialis, C. hyperboreus, C. marshallae, C. pacificus, and C. helgolandicus. Functional assignment of protein families based on clusters of orthologous genes (COG) and gene ontology (GO) annotations showed analogous patterns of protein functions across species. Phylogenetic analyses using maximum likelihood (ML) of 191 protein-coding genes mined from RNA-seq data fully resolved evolutionary relationships among seven Calanus species investigated (five species sequenced for this study and two species with published datasets), with gene and site concordance factors showing that 109 out of 191 protein-coding genes support a separation between three groups: the C. finmarchicus group (including C. finmarchicus, C. glacialis, and C. marshallae), the C. helgolandicus group (including C. helgolandicus, C. sinicus, and C. pacificus) and the monophyletic C. hyperboreus group. The tree topology obtained in ML analyses was similar to a previously proposed phylogeny based on morphological criteria and cleared certain ambiguities from past studies on evolutionary relationships among Calanus species.

The Innate Immune Response to Infection by Polyascus gregaria in the Male Chinese Mitten Crab (Eriocheir sinensis), Revealed by Proteomic Analysis

The Chinese mitten crab (Eriocheir sinensis) is a representative catadromous invertebrate of the Yangtze River and a commercial species widely cultivated in China. Both cultivated and wild crabs suffer from a variety of parasites and pathogens, which can result in catastrophic economic losses in aquaculture revenue. Polyascus gregaria, a parasitic barnacle with a highly derived morphology, is specialized in invading these crabs. This study examines the immunological mechanism in E. sinensis infected with P. gregaria. Tandem mass tags (TMT), a specialized method of mass-spectrometry, was used to analyze the infection by P. gregaria resistance at the protein level. In the hepatopancreas of infected crabs, 598 proteins differentially expressed relating to physiological change, of which, 352 were upregulated and 246 were downregulated. Based on this differential protein expression, 104 GO terms and 13 KEGG pathways were significantly enriched. Differentially expressed proteins, such as ATG, cathepsin, serpin, iron-related protein, Rab family, integrin, and lectin, are associated with the lysosome GO term and the autophagy-animal KEGG pathways, both of which likely relate to the immune response to the parasitic P. gregaria infection. These results show the benefit of taking a detailed, protein-level approach to understanding the innate immune response of aquatic invertebrates to macroparasite infection.

iTRAQ-based quantitative proteomic analysis of differentially expressed proteins in the hepatopancreas of Litopenaeus vannamei after WSSV infection

White spot syndrome virus (WSSV) is the most destructive virus among invertebrates. In this study, we analyzed the immune response after WSSV infection in Pacific white shrimp Litopenaeus vannamei using isobaric tags for relative and absolute quantitation (iTRAQ). We identified 325 differentially expressed proteins (DEPs) in the hepatopancreas of L. vannamei. Among them, 212 were up-regulated proteins, and several of them might be related to immunity (e.g. arginine kinase and peroxiredoxin). Of the 113 down-regulated proteins, some were related to immunity (e.g. cathepsin C and cathepsin L) and others to the antioxidant defense process (e.g. glutathione peroxidase and catalase). One down-regulated DEP (C7M84_014268) and 3 up-regulated DEPs (C7M84_003456, C7M84_020702, and C7M84_007135) were randomly selected and analyzed using parallel reaction monitoring. This study is an important step for a comprehensive understanding of the immune relationship between L. vannamei and WSSV and provides valuable information for the prevention of viral diseases in the crustacean aquaculture industry.

A newly identified Hippo homologue from the oriental river prawn Macrobrachium nipponense is involved in the antimicrobial immune response

The Hippo signalling pathway plays a vital role in organ size control, cell proliferation, apoptosis, and immune regulation. In this study, a Hippo homologue with three isoforms (named MnHippo-a, MnHippo-b, and MnHippo-c) was isolated and characterized for the first time from the freshwater prawn Macrobrachium nipponense. The deduced amino acid sequences of MnHippo-a (698 aa), MnHippo-b (688 aa), and MnHippo-c (656 aa) were highly similar, and they all contained an N-terminal S_TKc (serine/threonine protein kinase catalytic) domain and a C-terminal Mst1_SARAH (Sav/Rassf/Hpo) domain. MnHippo-a and MnHippo-c were derived from alternative splicing. Phylogenetic analysis was performed, and the results revealed that MnHippo was a member of the clade containing STPK4 and Hippo of Penaeus vannamei. The expression distribution showed that MnHippo was constitutively expressed in various tissues of uninfected prawns and highly expressed in the hepatopancreas and intestine. In prawns challenged with Vibrio parahaemolyticus and Staphylococcus aureus, the expression of MnHippo in haemocytes was significantly upregulated. Furthermore, in MnHippo-knockdown prawns injected with V. parahaemolyticus or S. aureus, the transcription levels of five antimicrobial peptides were downregulated. MnHippo silencing weakened the clearance of V. parahaemolyticus and S. aureus in prawns. The survival rate of the MnHippo-dsRNA group was obviously decreased from 2 to 6 days post-injection with V. parahaemolyticus or S. aureus. Hence, MnHippo might be involved in the antibacterial immune defence of M. nipponense.

Toxic effects of metal copper stress on immunity, metabolism and pathologic changes in Chinese mitten crab (Eriocheir japonica sinensis)

Copper (Cu²⁺), which represents a major physiological challenge for crab culture, is ubiquitous in the aquatic culture environment, and gills are the first organs that come into direct contact with the environment. However, the molecular basis of the response of crabs to Cu²⁺ stress remains unclear. Here, we conducted a transcriptome and differential expression analysis on the gills from Chinese mitten crab unexposed and exposed to Cu²⁺ for 24 h. The comparative transcriptome analysis identified 2486 differentially expressed genes (DEGs). GO functional analysis and KEGG pathway analysis revealed some DEGs, which were mostly related to immunity, metabolism, osmotic regulation, Cu²⁺ homeostasis regulation, antioxidant activity, and detoxification process. Some pathways related to humoral and cellular immunity, such as phagosome, peroxisome, lysosome, mTOR signaling pathway, PI3K-Akt signaling pathway, Toll-like receptor signaling pathway, and T cell receptor signaling pathway were enhanced under Cu²⁺ stress. In addition, Cu²⁺ stress altered the expression patterns of key phagocytosis and apoptosis genes (lectin, cathepsin L, Rab7, and HSP70), confirming that Cu²⁺ can induce oxidative stress and eventually even apoptosis. Histological analysis revealed that the copper can induce damage at the cellular level. This comparative transcriptome analysis provides valuable molecular information to aid future study of the immune mechanism of Chinese mitten crab in response to Cu²⁺ stress and provides a foundation for further understanding of the effects of metal toxicity.

Effects of dietary icariin supplementation on the ovary development-related transcriptome of Chinese mitten crab (Eriocheir sinensis)

The Chinese mitten crab (Eriocheir sinensis) is an economically important aquaculture species in China, with distinct differences in ovarian maturation status between crabs fed with natural diets and artificial diets during the listing period, thus, leading to selling price differentiation. Our previous study showed that dietary supplementation with 100 mg/kg icariin can effectively promote ovarian development of E. sinensis. However, the internal molecular mechanism has not yet been elucidated because of a lack of comprehensive genome sequence information. We compared the ovary transcriptomes of E. sinensis fed with two diets containing 0 and 100 mg/kg ICA using the BGISEQ-500 platform. This yielded 12.54 Gb clean bases and 54,794 unigenes, 13,832 of which were found to be differentially expressed after icariin exposure. Twenty pathways closely related to gonadal development were selected through KEGG analysis. Seven differentially expressed genes relevant to vitellogenesis and oocyte maturation (serine/threonine-protein kinase mos-like, Eg2, cytoplasmic polyadenylation element-binding protein, cyclin B, vitellogenin 1, cathepsin D, and juvenile hormone esterase-like carboxylesterase 1) were validated by qRT-PCR, and four proteins (MEK1/2, ERK1/2, Cyclin B and Cdc2) associated with the progesterone mediated oocyte maturation pathway (i.e., MAPK/MPF pathway) were analyzed by western-blot. The results showed that icariin could promote the synthesis, processing and deposition of vitellogenin in oocytes, and that it also has the potential to promote oocyte maturation (resumption of Meiosis I) by altering the expression of the relevant genes and proteins.

Hepatopancreas immune response during molt cycle in the mud crab, Scylla paramamosain

Molt is a critical developmental process in crustaceans. Recent studies have shown that the hepatopancreas is an important source of innate immune molecules, yet hepatopancreatic patterns of gene expression during the molt cycle which may underlie changes in immune mechanism are unknown. In this study, we performed Illumina sequencing for the hepatopancreas of the mud crab, Scylla paramamosain during molt cycle (pre-molt stage, post-molt stage, and inter-molt stage). A total of 44.55 Gb high-quality reads were obtained from the normalized cDNA of hepatopancreas. A total of 70,591 transcripts were assembled; 55,167 unigenes were identified. Transcriptomic comparison revealed 948 differentially expressed genes (DEGs) in the hepatopancreas from the three molt stages. We found that genes associated with immune response patterns changed in expression during the molt cycle. Antimicrobial peptide genes, inflammatory response genes, Toll signaling pathway factors, the phenoloxidase system, antioxidant enzymes, metal-binding proteins and other immune related genes are significantly up-regulated at the post-molt stage and inter-molt stage compared with the pre-molt stage, respectively. These genes are either not expressed or are expressed at low levels at the pre-molt stage. To our knowledge, this is the first systematic transcriptome analysis of genes capable of mobilizing a hepatopancreas immune response during the molt cycle in crustaceans, and this study will contribute to a better understanding of the hepatopancreas immune system and mud crab prophylactic immune mechanisms at the post-molt stage.

Proteome and transcriptome reveal the involvement of heat shock proteins and antioxidant system in thermotolerance of Clematis florida

Clematis florida Thun (CfT) is an ornamental and medicinal plant. It is a cold resistant but heat sensitive species and deserves to be further investigated to improve its adaptability to heat stress. Exploring the molecular mechanism potential via an omic-analysis constitutes a promising approach towards improving heat tolerance of CfT. Two CfT lines, heat resistance (HR) and heat sensitive (HS), with differential thermotolerance capacities were used for the integrative analyses of proteomics and transcriptomes. Transcriptomes analysis showed that various pathways were significantly enriched including plant hormone signal transduction and carbon fixation pathways in prokaryotes. Proteomics study revealed the enrichment of some other pathways comprising antioxidant activity and carbohydrates metabolism. Based on combined transcriptomes and proteomics analyses and following heat stress treatment, a total of 1724 annotated genes were overlapped between both CfT lines. Particularly, 84 differential expressed genes (DEGs) were overlapped in both CfT lines. Fifteen out of these 84 genes were up-regulated solely for HR line (PS) but not for HS one (SG). This strongly suggests a potential prominent role for these genes in the thermotolerance process in PS line. We corroborate that two Hsps (Hsp18 and Hsp70) out of 20 detected proteins with higher expression levels in PS than in SG based on either global transcripts or proteins levels. According to the transcriptomes and proteomics analyses, 6 proteins and their corresponding genes were found to be significantly abundant in HR line (PS). Data are available via ProteomeXchange with identifier PXD018192. The expressions levels of these 6 genes were checked also for both CfT lines to evaluate their potential contributions in the heat tolerance process. Thus, their expression levels were approximately 2~4 times higher in HR than in HS line. We provided as well a representative schematic model to highlight the key genes involved in ROS scavenging and photorespiratory pathway in CfT. This model could be helpful also in understanding the mechanism of heat tolerance in CfT.

Characterization of calcineurin A and B genes in the abalone, Haliotis diversicolor, and their immune response role during bacterial infection

Calcineurin (CN) is known to be involved in many biological processes, particularly, the immune response mechanism in many invertebrates. In this study, we characterized both HcCNA and HcCNB genes in Haliotis diversicolor, documented their expression in many tissues, and discerned their function as immune responsive genes against Vibrio parahaemolyticus infection. Similar to other mollusk CNs, the HcCNA gene lacked a proline-rich domain and comprised only one isoform of its catalytic unit, in contrast to CNs found in mammals. HcCNB was highly conserved in both sequence and domain architecture. Quantitative PCR and in situ hybridization revealed that the genes were broadly expressed and were not restricted to tissues traditionally associated with immune function. Upon infection of H. diversicolor with V. parahaemolyticus (a bacteria that causes serious disease in crustaceans and mollusks), both HcCNA and HcCNB genes were highly up-regulated at the early phase of bacterial infection. HcCNB was expressed significantly higher than HcCNA in response to bacterial challenge, suggesting its independent or more rapid response to bacterial infection. Together, the two CN genes are unique in their gene structure (particular HcCNA) and distribution in mollusk species and likely function as immune responsive genes along with many other genes that are enhanced in the early phase of V. parahaemolyticus infection in abalone.

Dynamic immune and metabolism response of clam Meretrix petechialis to Vibrio challenge revealed by a time series of transcriptome analysis

Meretrix petechialis is an important commercial aquaculture species in China. During the clam culture period, mass mortality events often occurred due to the Vibrio infection. In this paper, M. petechialis were challenged with Vibrio parahaemolyticus immersion to simulate a natural infection, and the infection process were divided into four phases including latency, prodrome, onset and recovery phases based on the clam mortality data. Then, the dynamic response of clams to Vibrio infection at different infection phases were investigated by transcriptome analysis. A total of 38,067 differentially expressed genes (DEGs) were identified at different infection phases. DEG annotations showed that immune-related and metabolism-related signaling pathways were enriched, indicating that immune defense and metabolism process play key roles during bacterial infection. Three kinds of expression pattern were classified by cluster analysis, including U-shape, L-shape and inverted V-shape. Anabolism and cellular growth proliferation related signaling pathways were repressed (long-lasting or transient) during bacterial infection. However, the immune related signaling pathways with different immune functions showed induction expression or repression expression against bacterial infection, which indicated that immune system take different strategies against bacterial infection. Furthermore, some signaling pathways such as PI3K-Akt signaling pathway both involved in immune defense and cell metabolism. This study provides a sight that the dynamic immunity and metabolic responses may be integrated to improve the host survival and shift more energy for immune defense.

Transcriptome analysis of response mechanism to ammonia stress in Asian clam (Corbicula fluminea)

Corbicula fluminea is highly sensitive to ammonia, and its response mechanism to ammonia stress is unclear. In this study, C. fluminea was exposed to different levels of ammonia (control group, 10 mg/L, and 25 mg/L) for 24 h and 48 h. A comparative analysis of transcriptome sequencing (RNA-seq) of C. fluminea digestive gland showed that the expression of 6742 genes (11.54%) was significantly affected by ammonia stress. The TLR, NF-κB, FOXO, and apoptotic signaling pathways were involved in the regulation. The differential expression of 14 genes was confirmed by real-time PCR. In summary, the response mechanism of C. fluminea digestive gland under ammonia stress may be different from that of oxidative stress in marine vertebrates. Also, the NMDAR-mediated pathway may not be the main mechanism in the response to ammonia stress in C. fluminea. The present study is a preliminary study for further investigation into ammonia toxicity in shellfish.

Comparative transcriptome analysis explores maternal to zygotic transition during Eriocheir sinensis early embryogenesis

The maternal genome directs almost all aspects of early animal development. As development proceeds, the elimination of maternal gene products and zygotic genome activation (ZGA) occur during the maternal to zygotic transition (MZT). To study the molecular mechanisms regulating this developmental event in Eriocheir sinensis, RNA-Seq technology was applied to generate comprehensive information on transcriptome dynamics during early embryonic stages. In total, 32,088 annotated unigenes were obtained from the transcriptomes of fertilized eggs and embryos at the cleavage (2-4 cell) and blastula stage. A total of 566 maternal genes and 1165 zygotic genes were isolated, among which 103 and 266 genes were predicted conserved maternal transcripts (COMATs) and conserved zygotic transcripts (COZYTs), respectively. The COMATs performed housekeeping gene functions and may be essential for initiating early embryogenesis of the Bilateria. Furthermore, 87, 76 and 117 differentially expressed genes associated with the MZT, morphogenesis and immunity were identified when compared the three transcriptomic datasets. We also unmask that the MZT takes place around the cleavage stage, when the genes involved in the clearance of maternal gene products and the ZGA were significantly up-regulated. Taken together, these datasets provide a valuable resource for understanding the mechanisms of early developmental events in E. sinensis, and facilitate further studies on molecular mechanisms of asynchronous development in crabs.

Transcriptome analysis of Aeromonas hydrophila infected hybrid sturgeon (Huso dauricus×Acipenser schrenckii)

The hybrid sturgeon (Huso dauricus × Acipenser schrenckii) is an economically important species in China. With the increasing aquaculture of hybrid sturgeon, the bacterial diseases are a great concern of the industry. In this study, de novo sequencing was used to compare the difference in transcriptome in spleen of the infected and mock infected sturgeon with Aeromonas hydrophila. Among 187,244 unigenes obtained, 87,887 unigenes were annotated and 1,147 unigenes were associated with immune responses genes. Comparative expression analysis indicated that 2,723 differently expressed genes between the infected and mock-infected group were identified, including 1,420 up-regulated and 1,303 down-regulated genes. 283 differently expressed anti-bacterial immune related genes were scrutinized, including 168 up-regulated and 115 down-regulated genes. Ten of the differently expressed genes were further validated by qRT-PCR. In this study, toll like receptors (TLRs) pathway, NF-kappa B pathway, class A scavenger receptor pathway, phagocytosis pathway, mannose receptor pathway and complement pathway were shown to be up-regulated in Aeromonas hydrophila infected hybrid sturgeon. Additionally, 65,040 potential SSRs and 2,133,505 candidate SNPs were identified from the hybrid sturgeon spleen transcriptome. This study could provide an insight of host immune genes associated with bacterial infection in hybrid sturgeon.

Transcriptomic analysis of immune-related genes in the lipopolysaccharide-stimulated hepatopancreas of the mudflat crab Helice tientsinensis

The mudflat crab Helice tientsinensis is one of the most economically important aquaculture species in China. Nevertheless, it is susceptible to various diseases caused by viruses, bacteria and rickettsia-like organisms. A better understanding of the immune system and genes related to the responses to bacterial and viral infection is required. Herein, the hepatopancreas transcriptome of H. tientsinensis was analyzed by comparing control and lipopolysaccharide (LPS)-stimulated RNA-Seq data, yielding 91,885,038 bp and 13.78 Gb of clean reads. Following assembly and annotation, 93,207 unigenes with an average length of 883 bp were identified, of which 31,674 and 13,700 were annotated in Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases, respectively. Following LPS, 4845 differentially expressed genes (DEGs) were identified, of which 2491 and 2354 were up- and down-regulated, respectively. To further investigate immune-related DEGs, KEGG enrichment analysis identified immune response pathways, most notably the peroxisome and Toll-like receptor signaling pathways. Quantitative real time-PCR (qRT-PCR) confirmed the up-regulation of a random selection of DEGs. This systematic transcriptomic analysis of the innate immune pathway in H. tientsinensis expands our understanding of the immune system in crabs.

Transcriptome analysis and discovery of genes involved in immune pathways from coelomocytes of Onchidium struma after bacterial challenge

Onchidium struma widely distributes in subtidal and low-tidal zones, which is considered to be an economical species with rich nutrition, a valuable biomonitor for heavy metal pollution and a representative species for evolution from ocean to land. However, there is limited genetic information available for O. struma development. This study compared transcriptomic profiles of coelomocytes from normal and bacteria infected O. struma by Illumina-based paired-end sequencing to explore the molecular immune mechanism of O. struma against bacterial infection. After assembly, a total of 92,450 unigenes with an average length of 1019 bp were obtained. Approximately 34,964 (37.82%) unigenes were annotated in the Nr NCBI database and 40.1% of unigenes were similar with that of Aplysia californica. Among them, 7609 unigenes were classified into three Gene Ontology (GO) categories: biological process (3250 unigenes, 42.7%), cellular component (2,281, 30.0%) and molecular function (2078 unigenes, 27.3%). A total of 22,776 unigenes were aligned to the Clusters of Orthologous Groups (COG) of proteins and classified into 25 functional categories. Following bacterial infection, 10,623 differently expressed unigenes (DEGs) were identified, including 7644 up-regulated and 2979 down-regulated unigenes. Further KEGG analysis annotated 11,681 DEGs to 42 pathways, and 11 pathways were identified to be related with diseases and immune system. To our knowledge, it was first time to analyze transcriptome profiles of O. struma. Results of the present study will provide valuable theoretical resources for future genetic and genomic research on O. struma. The research results will be helpful for improving the efficiency and quality of artificial breeding, establishing genetic linkage map, and enhancing health management for this species.

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.

Transcriptional responses in the hepatopancreas of Eriocheir sinensis exposed to deltamethrin

Deltamethrin is an important pesticide widely used against ectoparasites. Deltamethrin contamination has resulted in a threat to the healthy breeding of the Chinese mitten crab, Eriocheir sinensis. In this study, we investigated transcriptional responses in the hepatopancreas of E. sinensis exposed to deltamethrin. We obtained 99,087,448, 89,086,478, and 100,117,958 raw sequence reads from control 1, control 2, and control 3 groups, and 92,094,972, 92,883,894, and 92,500,828 raw sequence reads from test 1, test 2, and test 3 groups, respectively. After filtering and quality checking of the raw sequence reads, our analysis yielded 79,228,354, 72,336,470, 81,859,826, 77,649,400, 77,194,276, and 75,697,016 clean reads with a mean length of 150 bp from the control and test groups. After deltamethrin treatment, a total of 160 and 167 genes were significantly upregulated and downregulated, respectively. Gene ontology terms "biological process," "cellular component," and "molecular function" were enriched with respect to cell killing, cellular process, other organism part, cell part, binding, and catalytic. Pathway analysis using the Kyoto Encyclopedia of Genes and Genomes showed that the metabolic pathways were significantly enriched. We found that the CYP450 enzyme system, carboxylesterase, glutathione-S-transferase, and material (including carbohydrate, lipid, protein, and other substances) metabolism played important roles in the metabolism of deltamethrin in the hepatopancreas of E. sinensis. This study revealed differentially expressed genes related to insecticide metabolism and detoxification in E. sinensis for the first time and will help in understanding the toxicity and molecular metabolic mechanisms of deltamethrin in E. sinensis.

Transcriptomic analysis of porcine PBMCs in response to FMDV infection

Foot-and-mouth disease (FMD) is a significant zoonotic infectious disease. It has an important economic impact throughout the world. As well, it is a considerable threat to food security. At present, the molecular mechanism of FMDV infection is not clear to a large extent. Innate immune response is the first line of defense against infectious diseases. The systematic analysis of the host immune response to infection has an important role in understanding the pathogenesis of infection. However, there are few reports about effect of immune regulation on virus replication in the interaction of virus and host cellular. High-throughput RNA-seq technology as a powerful and efficient means for transcript analysis provides a new insight into FMDV study. In this study, RNA extracted from pig PBMCs infected with O subtype FMDV at 4 dpi. A total of 29942658 and 31452917 Illumina read pairs were obtained from the non-infected (NI) group and infected (I) group, respectively. The clean bases for all samples are 3.61G (NI group) and 3.79G (I group), respectively. The clean reads of the NI and I group that mapped to pig genome data were 47195073 (81.82%) and 46556714 (76.85%), respectively. Most of the clean reads were distributed in the exon region, followed by intron region and intergenic region. Differently expressed (DE) genes were analyzed using edgeR software. 451 genes were differentially expressed between the infected and the non-infected groups. According to the comparison analysis, more genes were down-regulated in the non-infected samples than in those infected with FMDV.66 out of 451 genes were down-regulated, 385 out of 451 genes were up-regulated following FMDV infection. For function classification and pathway analysis, among 17741 assembled unigenes, there are 349 genes which are different genes of GO notes. Moreover, 49 genes were down-regulated, 300 genes were up-regulated associate with GO term. 1621 were successfully annotated by GO assignments, belonging to one or more of the three categories: biological process, cellular component, and molecular function. According to KEGG analysis,the main pathway was represented including protein processing in endoplasmic reticulum, phagosome, cell cycle and cytokine-cytokine receptor interaction. Some key DE genes related to immune process and signaling pathways were analyzed and quantified by RT-PCR. This is the first systematical transcriptome analysis of pig PBMCs infected by FMDV. These findings will help us better understand the host Cell-FMDV interaction and its relationship to pathogenesis, as well as contribute to the prevention and control of FMDV.

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.

Biomolecular changes that occur in the antennal gland of the giant freshwater prawn (Machrobrachium rosenbergii)

In decapod crustaceans, the antennal gland (AnG) is a major primary source of externally secreted biomolecules, and some may act as pheromones that play a major role in aquatic animal communication. In aquatic crustaceans, sex pheromones regulate reproductive behaviours, yet they remain largely unidentified besides the N-acetylglucosamine-1,5-lactone (NAGL) that stimulates male to female attraction. In this study, we used an AnG transcriptome of the female giant freshwater prawn (Macrobrachium rosenbergii) to predict the secretion of 226 proteins, including the most abundantly expressed transcripts encoding the Spaetzle protein, a serine protease inhibitor, and an arthropodial cuticle protein AMP 8.1. A quantitative proteome analysis of the female AnG at intermolt, premolt and postmolt, identified numerous proteins of different abundances, such as the hemocyanin subunit 1 that is most abundant at intermolt. We also show that hemocyanin subunit 1 is present within water surrounding females. Of those metabolites identified, we demonstrate that the NAGL and N-acetylglucosamine (NAG) can bind with high affinity to hemocyanin subunit 1. In summary, this study has revealed components of the female giant freshwater prawn AnG that are released and contribute to further research towards understanding crustacean conspecific signalling.

A second generation SNP and SSR integrated linkage map and QTL mapping for the Chinese mitten crab Eriocheir sinensis

The Chinese mitten crab Eriocheir sinensis is the most economically important cultivated crab species in China, and its genome has a high number of chromosomes (2n = 146). To obtain sufficient markers for construction of a dense genetic map for this species, we employed the recently developed specific-locus amplified fragment sequencing (SLAF-seq) method for large-scale SNPs screening and genotyping in a F1 full-sib family of 149 individuals. SLAF-seq generated 127,677 polymorphic SNP markers, of which 20,803 valid markers were assigned into five segregation types and were used together with previous SSR markers for linkage map construction. The final integrated genetic map included 17,680 SNP and 629 SSR markers on the 73 linkage groups (LG), and spanned 14,894.9 cM with an average marker interval of 0.81 cM. QTL mapping localized three significant growth-related QTL to a 1.2 cM region in LG53 as well as 146 sex-linked markers in LG48. Genome-wide QTL-association analysis further identified four growth-related QTL genes named LNX2, PAK2, FMRFamide and octopamine receptors. These genes are involved in a variety of different signaling pathways including cell proliferation and growth. The map and SNP markers described here will be a valuable resource for the E. sinensis genome project and selective breeding programs.

In vivo activity and the transcriptional regulatory mechanism of the antimicrobial peptide SpHyastatin in Scylla paramamosain

A new gene homologous to the reported antimicrobial peptide (AMP) hyastatin from Hyas araneus was screened in the SSH library constructed from the hemocytes of Scylla paramamosain, and named SpHyastatin. In vivo study showed that SpHyastatin was predominantly expressed in hemocytes of S. paramamosain. With the challenge of either Vibrio parahaemolyticus or lipopolysaccharide (LPS), SpHyastatin showed a positive response, meaning that it was probably involved in the immune reaction against bacterial infection in vivo. A distinctive feature of SpHyastatin in comparison with six other known AMPs tested was that SpHyastatin could maintain a higher transcription level from megalopas to the adult crab, indicating a potential consistent resistance against pathogens conferred by this peptide existing in the blood circulation of crabs. RNA interference assay was performed to inhibit SpHyastatin transcription in vivo and the result demonstrated that silencing SpHyastatin mRNA transcripts could decrease the survival rate of crabs challenged with V. parahaemolyticus. To further understand the molecular mechanisms that regulate SpHyastatin expression, a 576 bp 5'-flanking sequence of SpHyastatin was obtained using genome walking. Here, we focused our experiments on investigating the roles of the putative NF-κB binding site in LPS-mediated transcriptional regulation of the SpHyastatin gene using endothelial progenitor cells and Hela cells. Luciferase reporter analyses demonstrated that the putative NF-κB element acted as a positive regulatory element and was essential for the induction of SpHyastatin promoter by LPS. These results should shed light on the in vivo functional property and the molecular mechanism of regulation for the crab AMP SpHyastatin.

Diversity of insulin-like peptide signaling system proteins in Calanus finmarchicus (Crustacea; Copepoda) - Possible contributors to seasonal pre-adult diapause

Calanus finmarchicus, an abundant calanoid copepod in the North Atlantic Ocean, is both a major grazer on phytoplankton and an important forage species for invertebrate and vertebrate predators. One component of the life history of C. finmarchicus is the overwintering dormancy of sub-adults, a feature key for the annual recruitment of this species in early spring. While little is known about the control of dormancy in C. finmarchicus, one hypothesis is that it is an insect-like diapause, where the endocrine system is a key regulator. One group of hormones implicated in the control of insect diapause is the insulin-like peptides (ILPs). Here, C. finmarchicus transcriptomic data were used to predict ILP signaling pathway proteins. Four ILP precursors were identified, each possessing a distinct A- and B-chain peptide; these peptides are predicted to form bioactive heterodimers via inter-chain disulfide bridging. Two ILP receptors, which likely represent splice variants of a common gene, were identified. Three insulin-degrading enzymes were also discovered, as were proteins encoding the transcription factor FOXO, a downstream target of ILP that has been implicated in the regulation of insect diapause, and insulin receptor substrate, a protein putatively linking the ILP receptor and FOXO. RNA-Seq data suggest that some C. finmarchicus insulin pathway transcripts are differentially expressed across development. As in insects, the ILP signaling system may be involved in controlling C. finmarchicus' organism-environment interactions (e.g., regulation of seasonal sub-adult diapause), a hypothesis that can now be investigated using these data.

Transcriptome Analysis Revealed Changes of Multiple Genes Involved in Haliotis discus hannai Innate Immunity during Vibrio parahemolyticus Infection

Abalone (Haliotis discus hannai) is one of the most valuable marine aquatic species in Korea, Japan and China. Tremendous exposure to bacterial infection is common in aquaculture environment, especially by Vibrio sp. infections. It’s therefore necessary and urgent to understand the mechanism of H. discus hannai host defense against Vibrio parahemolyticus infection. However studies on its immune system are hindered by the lack of genomic resources. In the present study, we sequenced the transcriptome of control and bacterial challenged H. discus hannai tissues. Totally, 138 MB of reference transcriptome were obtained from de novo assembly of 34 GB clean bases from ten different libraries and annotated with the biological terms (GO and KEGG). A total of 10,575 transcripts exhibiting the differentially expression at least one pair of comparison and the functional annotations highlight genes related to immune response, cell adhesion, immune regulators, redox molecules and mitochondrial coding genes. Mostly, these groups of genes were dominated in hemocytes compared to other tissues. This work is a prerequisite for the identification of those physiological traits controlling H. discus hannai ability to survive against Vibrio infection.

Identification of Genes Relevant to Pesticides and Biology from Global Transcriptome Data of Monochamus alternatus Hope (Coleoptera: Cerambycidae) Larvae

Monochamus alternatus Hope is the main vector in China of the Pine Wilt Disease caused by the pine wood nematode Bursaphelenchus xylophilus. Although chemical control is traditionally used to prevent pine wilt disease, new strategies based in biological control are promising ways for the management of the disease. However, there is no deep sequence analysis of Monochamus alternatus Hope that describes the transcriptome and no information is available about gene function of this insect vector. We used next generation sequencing technology to sequence the whole fourth instar larva transcriptome of Monochamus alternatus Hope and successfully built a Monochamus alternatus Hope transcriptome database. In total, 105,612 unigenes were assigned for Gene Ontology (GO) terms, information for 16,730 classified unigenes was obtained in the Clusters of Orthologous Groups (COGs) database, and 13,024 unigenes matched with 224 predicted pathways in the Kyoto Encyclopedia of Genes and Genome (KEGG). In addition, genes related to putative insecticide resistance-related genes, RNAi, the Bt receptor, intestinal digestive enzymes, possible future insect control targets and immune-related molecules are described. This study provides valuable basic information that can be used as a gateway to develop new molecular tools for Monochamus alternatus Hope control strategies.

Transcriptome assembly and expression profiling of molecular responses to cadmium toxicity in hepatopancreas of the freshwater crab Sinopotamon henanense

Cadmium (Cd) pollution is a serious global problem, which causes irreversible toxic effects on animals. Freshwater crab, Sinopotamon henanense, is a useful environmental indicator since it is widely distributed in benthic habitats whereby it tends to accumulate Cd and other toxicants. However, its molecular responses to Cd toxicity remain unclear. In this study, we performed transcriptome sequencing and gene expression analyses of its hepatopancreas with and without Cd treatments. A total of 7.78 G clean reads were obtained from the pooled samples, and 68,648 unigenes with an average size of 622 bp were assembled, in which 5,436 were metabolism-associated and 2,728 were stimulus response-associated that include 380 immunity-related unigenes. Expression profile analysis demonstrated that most genes involved in macromolecular metabolism, oxidative phosphorylation, detoxification and anti-oxidant defense were up-regulated by Cd exposure, whereas immunity-related genes were down-regulated, except the genes involved in phagocytosis were up-regulated. The current data indicate that Cd exposure alters gene expressions in a concentration-dependent manner. Therefore, our results provide the first comprehensive S.henanense transcriptome dataset, which is useful for biological and ecotoxicological studies on this crab and its related species at molecular level, and some key Cd-responsive genes may provide candidate biomarkers for monitoring aquatic pollution by heavy metals.

Identification and characterization of nascent polypeptide-associated complex alpha from Chinese mitten crab (Eriocheir sinensis): A novel stress and immune response gene in crustaceans

Disease in aquatic animals is tightly linked to environmental challenges and their immune responses are greatly modified by their external environment. The chaperone protein nascent polypeptide-associated complex alpha (NACA) has been suggested to play important roles in the cellular response to stress and immune challenges, while the related biological functions remain largely unknown in invertebrates. In the present study we identified a NACA gene (termed EsNACA) from Chinese mitten crab Eriocheir sinensis and analyzed its expression changes in response to ambient (salinity and pH) stresses and immune challenges. The EsNACA protein is comprised of 209 amino acid residues with a conserved DNA binding domain, a C-Jun binding domain, a NAC domain and an ubiquitin-associated domain and shows the highest sequence identity (87%) with its counterpart in shrimp Penaeus monodon. EsNACA mRNA transcripts are presented in all tested normal tissues with predominant expression in hepatopancreas and lower expression in hemocytes. In addition, EsNACA expression was significantly altered in response to the ambient salinity (15‰ and 30‰ salinities) and pH (pH 6 and 8.5) stresses in gill, hepatopancreas, muscle, hemocytes and intestine tissues. Furthermore, EsNACA gene expression was substantially induced upon LPS and Poly(I:C) immune stimulations in E. sinensis hemocytes in vitro. Finally, EsNACA expression was up-regulated in E. sinensis hemocytes, gill, hepatopancreas, intestine and muscle tissues in response to Vibrio anguillarum challenges in vivo. Taken together, our findings for the first time show that EsNACA is an inducible gene involved in stress and immune response in crustaceans.

Large-Scale Transcriptome Analysis of Cucumber and Botrytis cinerea during Infection

Cucumber gray mold caused by Botrytis cinerea is considered one of the most serious cucumber diseases. With the advent of Hi-seq technology, it is possible to study the plant-pathogen interaction at the transcriptome level. To the best of our knowledge, this is the first application of RNA-seq to identify cucumber and B. cinerea differentially expressed genes (DEGs) before and after the plant-pathogen interaction. In total, 248,908,688 raw reads were generated; after removing low-quality reads and those containing adapter and poly-N, 238,341,648 clean reads remained to map the reference genome. There were 3,512 cucumber DEGs and 1,735 B. cinerea DEGs. GO enrichment and KEGG enrichment analysis were performed on these DEGs to study the interaction between cucumber and B. cinerea. To verify the reliability and accuracy of our transcriptome data, 5 cucumber DEGs and 5 B. cinerea DEGs were chosen for RT-PCR verification. This is the first systematic transcriptome analysis of components related to the B. cinerea-cucumber interaction. Functional genes and putative pathways identified herein will increase our understanding of the mechanism of the pathogen-host interaction.

Transcriptome analysis of genes involved in defence response in Polyporus umbellatus with Armillaria mellea infection

Polyporus umbellatus, a species symbiotic with Armillaria mellea and it also exhibits substantial defence response to Armillaria mellea infection. There are no genomics resources databases for understanding the molecular mechanism underlying the infection stress of P. umbellatus. Therefore, we performed a large-scale transcriptome sequencing of this fungus with A. mellea infection using Illumina sequencing technology. The assembly of the clean reads resulted in 120,576 transcripts, including 38,444 unigenes. Additionally, we performed a gene expression profiling analysis upon infection treatment. The results indicated significant differences in the gene expression profiles between the control and the infection group. In total, 10933 genes were identified between the two groups. Based on the differentially expressed genes, a Gene Ontology annotation analysis showed many defence-relevant categories. Meanwhile, the Kyoto Encyclopedia of Genes and Genomes pathway analysis uncovered some important pathways. Furthermore, the expression patterns of 13 putative genes that are involved in defence response resulting from quantitative real-time PCR were consistent with their transcript abundance changes as identified by RNA-seq. The sequenced genes covered a considerable proportion of the P. umbellatus transcriptome, and the expression results may be useful to strengthen the knowledge on the defence response of this fungus defend against Armillaria mellea invasion.

Ovarian Transcriptome Analysis of Portunus trituberculatus Provides Insights into Genes Expressed during Phase III and IV Development

Enhancing the production of aquatic animals is crucial for fishery management and aquaculture applications. Ovaries are specialized tissues that play critical roles in producing oocytes and hormones. Significant biochemical changes take place during the sexual maturation of Portunus trituberculatus, but the genetics of this process has not been extensively studied. Transcriptome sequencing can be used to determine gene expression changes within specific periods. In the current study, we used transcriptome sequencing to produce a comprehensive transcript dataset for the ovarian development of P. trituberculatus. Approximately 100 million sequencing reads were generated, and 126,075 transcripts were assembled. Functional annotation of the obtained transcripts revealed important pathways in ovarian development, such as those involving the vitellogenin gene. Also, we performed deep sequencing of ovaries in phases III and IV of sexual maturation in P. trituberculatus. Differential analysis of gene expression identified 506 significantly differentially expressed genes, which belong to 20 pathway, transporters, development, transcription factors, metabolism of other amino acids, carbohydrate and lipid, solute carrier family members, and enzymes. Taken together, our study provides the first comprehensive transcriptomic resource for P. trituberculatus ovaries, which will strengthen understanding of the molecular mechanisms underlying the sexual maturation process and advance molecular nutritional studies of P. trituberculatus.

Transcriptome Analysis and Discovery of Genes Involved in Immune Pathways from Coelomocytes of Sea Cucumber (Apostichopus japonicus) after Vibrio splendidus Challenge

Vibrio splendidus is identified as one of the major pathogenic factors for the skin ulceration syndrome in sea cucumber (Apostichopus japonicus), which has vastly limited the development of the sea cucumber culture industry. In order to screen the immune genes involving Vibrio splendidus challenge in sea cucumber and explore the molecular mechanism of this process, the related transcriptome and gene expression profiling of resistant and susceptible biotypes of sea cucumber with Vibrio splendidus challenge were collected for analysis. A total of 319,455,942 trimmed reads were obtained, which were assembled into 186,658 contigs. After that, 89,891 representative contigs (without isoform) were clustered. The analysis of the gene expression profiling identified 358 differentially expression genes (DEGs) in the bacterial-resistant group, and 102 DEGs in the bacterial-susceptible group, compared with that in control group. According to the reported references and annotation information from BLAST, GO and KEGG, 30 putative bacterial-resistant genes and 19 putative bacterial-susceptible genes were identified from DEGs. The qRT-PCR results were consistent with the RNA-Seq results. Furthermore, many DGEs were involved in immune signaling related pathways, such as Endocytosis, Lysosome, MAPK, Chemokine and the ERBB signaling pathway.

Molecular characterization and expression profile of three Fem-1 genes in Eriocheir sinensis provide a new insight into crab sex-determining mechanism

The FEM-1 protein of Caenorhabditis elegans plays a crucial role in the nematode sex-determination pathway. Here, we reported the characterization of three members of Fem-1 gene family in Eriocheir sinensis (designated EsFem-1a, EsFem-1b, and EsFem-1c), which were homologs of the nematode FEM-1 protein. The amino acid sequences of EsFem-1a, EsFem-1b, and EsFem-1c contained eight, nine, and eight ankyrin repeats, respectively. None of the ankyrin repeats had its own specific signature, and the evolution of ankyrin repeat was not completely independent. The predicted three-dimensional structure of EsFem-1 proteins exhibited highly similar superhelical conformation, especially the N-terminal six contiguous ankyrin repeats, which provided a binding surface for the protein-protein interaction. Phylogenetic tree based on the amino acid sequences revealed that EsFem-1a, EsFem-1b, and EsFem-1c were divided into three obvious separated clades. EsFem-1 genes were highly expressed in fertilized egg, 2-4 cell and blastula stage comparing with larval stage (P<0.01), which suggested they might be maternal genes. They also showed a certain degree of sexually dimorphic expression in some tissues. Notably, the highest expression of EsFem-1a was in the hepatopancreas, with EsFem-1b in testes and EsFem-1c in muscle (P<0.05), which indicated their potential role in a broad array of tissues. In addition, the genes initially involved in sex differentiation were not limited to those specifically expressed in the developing gonad. Taken together, these results suggested that EsFem-1 might function in crab early sex determination and late gonad development. The identification of Fem-1 gene family in E. sinensis provides a new insight into crab sex-determination mechanism.

Transcriptome analysis and microsatellite discovery in the blunt snout bream (Megalobrama amblycephala) after challenge with Aeromonas hydrophila

The blunt snout bream, Megalobrama amblycephala, is a herbivorous freshwater fish species native to China and a major aquaculture species in Chinese freshwater polyculture systems. In recent years, the bacterium Aeromonas hydrophila has been reported to be its pathogen causing great losses of farmed fish. To understand the immune response of the blunt snout bream to A. hydrophila infection, we used the Solexa/Illumina technology to analyze the transcriptomic profile after artificial bacterial infection. Two nonnormalized cDNA libraries were synthesized from tissues collected from control blunt snout bream or those injected with A. hydrophila. After assembly, 155,052 unigenes (average length 692.8 bp) were isolated. All unigenes were annotated using BLASTX relative to several public databases: the National Center for Biotechnology Information nonreduntant (Nr) database, SwissProt, Eukaryotic Orthologous Groups of proteins (KOG), Kyoto Encyclopedia of Genes and Genomes (KEGG), and Gene Ontology (GO). The sequence similarity (86%) of the assembled unigenes was to zebrafish based on the Nr database. A number of unigenes (n = 30,482) were assigned to three GO categories: biological processes (25,242 unigenes), molecular functions (26,096 unigenes), and cellular components (22,778 unigenes). 20,909 unigenes were classified into 25 KOG categories and 28,744 unigenes were assigned into 315 specific signaling pathways. In total, 238 significantly differentially expressed unigenes (mapped to 125 genes) were identified: 101 upregulated genes and 24 downregulated genes. Another 303 unigenes were mapped to unknown or novel genes. Among the known expressed genes identified, 53 were immune-related genes and were distributed in 71 signaling pathways. The expression patterns of selected up- and downregulated genes from the control and injected groups were determined with reverse transcription-quantitative PCR (RT-qPCR). Microsatellites (n = 10,877), including di-to pentanucleotide repeat motifs, were also identified in the blunt snout bream transcriptome profiles. This study extends our understanding of the immune defense mechanisms of the blunt snout bream against A. hydrophila and provides useful data for further studies of the immunogenetics of this species.

De novo assembly of the Carcinus maenas transcriptome and characterization of innate immune system pathways

Background

The European shore crab, Carcinus maenas, is used widely in biomonitoring, ecotoxicology and for studies into host-pathogen interactions. It is also an important invasive species in numerous global locations. However, the genomic resources for this organism are still sparse, limiting research progress in these fields. To address this resource shortfall we produced a C. maenas transcriptome, enabled by the progress in next-generation sequencing technologies, and applied this to assemble information on the innate immune system in this species.

Results

We isolated and pooled RNA for twelve different tissues and organs from C. maenas individuals and sequenced the RNA using next generation sequencing on an Illumina HiSeq 2500 platform. After de novo assembly a transcriptome was generated encompassing 212,427 transcripts (153,699 loci). The transcripts were filtered, annotated and characterised using a variety of tools (including BLAST, MEGAN and RSEM) and databases (including NCBI, Gene Ontology and KEGG). There were differential patterns of expression for between 1,223 and 2,741 transcripts across tissues and organs with over-represented Gene Ontology terms relating to their specific function. Based on sequence homology to immune system components in other organisms, we show both the presence of transcripts for a series of known pathogen recognition receptors and response proteins that form part of the innate immune system, and transcripts representing the RNAi, Toll-like receptor signalling, IMD and JAK/STAT pathways.

Conclusions

We have produced an assembled transcriptome for C. maenas that provides a significant molecular resource for wide ranging studies in this species. Analysis of the transcriptome has revealed the presence of a series of known targets and functional pathways that form part of their innate immune system and illustrate tissue specific differences in their expression patterns.

Electronic supplementary material

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

High-density linkage mapping aided by transcriptomics documents ZW sex determination system in the Chinese mitten crab Eriocheir sinensis

The sex determination system in crabs is believed to be XY-XX from karyotypy, but centromeres could not be identified in some chromosomes and their morphology is not completely clear. Using quantitative trait locus mapping of the gender phenotype, we revealed a ZW-ZZ sex determination system in Eriocheir sinensis and presented a high-density linkage map covering ~98.5% of the genome, with 73 linkage groups corresponding to the haploid chromosome number. All sex-linked markers in the family we used were located on a single linkage group, LG60, and sex linkage was confirmed by genome-wide association studies (GWAS). Forty-six markers detected by GWAS were heterozygous and segregated only in the female parent. The female LG60 was thus the putative W chromosome, with the homologous male LG60 as the Z chromosome. The putative Z and W sex chromosomes were identical in size and carried many homologous loci. Sex ratio (5:1) skewing towards females in induced triploids using unrelated animals also supported a ZW-ZZ system. Transcriptome data were used to search for candidate sex-determining loci, but only one LG60 gene was identified as an ankyrin-2 gene. Double sex- and mab3-related transcription factor 1 (Dmrt1), a Z-linked gene in birds, was located on a putative autosome. With complete genome sequencing and transcriptomic data, more genes on putative sex chromosomes will be characterised, thus leading towards a comprehensive understanding of the sex determination and differentiation mechanisms of E. sinensis, and decapod crustaceans in general.

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

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

The invertebrate midintestinal gland ("hepatopancreas") is an evolutionary forerunner in the integration of immunity and metabolism

The immune system has an impact on the metabolic performance in vertebrates, thus the metabolic effects of immune cells are receiving intense attention today in the biomedical field. However, the evolutionary origin of the immunity-metabolism interaction is still uncertain. In this review, I show that mollusks and crustaceans integrate immune functions to a metabolic organ, the midintestinal gland ("hepatopancreas"). In these animals, the epithelial cells of the midintestinal gland are major sources of immune molecules, such as lectins, hemocyanin, ferritin, antibacterial and antiviral proteins, proteolytic enzymes and nitric oxide. There is crosstalk between midintestinal gland cells and phagocytes, which aids the initiation of the immune response and the clearance of pathogens. The midintestinal gland is thereby an integrated organ of immunity and metabolism. It is likely that immunity was the primary function of the midintestinal gland cells and that their role in the intermediate metabolism has evolved during the course of their further specialization.

The protein-protein interaction network of eyestalk, Y-organ and hepatopancreas in Chinese mitten crab Eriocheir sinensis

Background

The protein-protein interaction network (PIN) is an effective information tool for understanding the complex biological processes inside the cell and solving many biological problems such as signaling pathway identification and prediction of protein functions. Eriocheir sinensis is a highly-commercial aquaculture species with an unclear proteome background which hinders the construction and development of PIN for E. sinensis. However, in recent years, the development of next-generation deep-sequencing techniques makes it possible to get high throughput data of E. sinensis tanscriptome and subsequently obtain a systematic overview of the protein-protein interaction system.

Results

In this work we sequenced the transcriptional RNA of eyestalk, Y-organ and hepatopancreas in E. sinensis and generated a PIN of E. sinensis which included 3,223 proteins and 35,787 interactions. Each protein-protein interaction in the network was scored according to the homology and genetic relationship. The signaling sub-network, representing the signal transduction pathways in E. sinensis, was extracted from the global network, which depicted a global view of the signaling systems in E. sinensis. Seven basic signal transduction pathways were identified in E. sinensis. By investigating the evolution paths of the seven pathways, we found that these pathways got mature in different evolutionary stages. Moreover, the functions of unclassified proteins and unigenes in the PIN of E. sinensis were predicted. Specifically, the functions of 549 unclassified proteins related to 864 unclassified unigenes were assigned, which respectively covered 76% and 73% of all the unclassified proteins and unigenes in the network.

Conclusions

The PIN generated in this work is the first large-scale PIN of aquatic crustacean, thereby providing a paradigmatic blueprint of the aquatic crustacean interactome. Signaling sub-network extracted from the global PIN depicts the interaction of different signaling proteins and the evolutionary paths of the identified signal transduction pathways. Furthermore, the function assignment of unclassified proteins based on the PIN offers a new reference in protein function exploration. More importantly, the construction of the E. sinensis PIN provides necessary experience for the exploration of PINs in other aquatic crustacean species.

Transcriptome profiling analysis on whole bodies of microbial challenged Eriocheir sinensis larvae for immune gene identification and SNP development

To study crab immunogenetics of individuals, newly hatched Eriocheir sinensis larvae were stimulated with a mixture of three pathogen strains (Gram-positive bacteria Micrococcus luteus, Gram-negative bacteria Vibrio alginolyticus and fungi Pichia pastoris; 10⁸ cfu·mL^-1). A total of 44,767,566 Illumina clean reads corresponding to 4.52 Gb nucleotides were generated and assembled into 100,252 unigenes (average length: 1,042 bp; range: 201-19,357 bp). 17,097 (26.09%) of 65,535 non-redundant unigenes were annotated in NCBI non-redundant protein (Nr) database. Moreover, 23,188 (35.38%) unigenes were assigned to three Gene Ontology (GO) categories, 15,071 (23.00%) to twenty-six Clusters of orthologous Groups (COG) and 8,574 (13.08%) to six Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, respectively. Numerous genes were further identified to be associated with multiple immune pathways, including Toll, immune deficiency (IMD), janus kinase (JAK)-signal transducers and activators of transcription (STAT) and mitogen-activated protein kinase (MAPK) pathways. Some of them, such as tumor necrosis factor receptor associated factor 6 (TRAF6), fibroblast growth factor (FGF), protein-tyrosine phosphatase (PTP), JNK-interacting protein 1 (JIP1), were first identified in E. sinensis. TRAF6 was even first discovered in crabs. Additionally, 49,555 single nucleotide polymorphisms (SNPs) were developed from over 13,309 unigenes. This is the first transcriptome report of whole bodies of E. sinensis larvae after immune challenge. Data generated here not only provide detail information to identify novel genes in genome reference-free E. sinensis, but also facilitate our understanding on host immunity and defense mechanism of the crab at whole transcriptome level.