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

Transcriptome analysis of differentially expressed genes in the red swamp crayfish Procambarus clarkii challenged with Aeromonas hydrophila

As an important economic species in China, aquaculture of the crayfish Procambarus clarkii has suffered huge losses due to infection by pathogenic bacteria, mainly by Aeromonas hydrophila, which leads to high mortality and huge economic loss. To better understand the immune response of crayfish against bacterial infection, we compared and analyzed transcriptome data of hepatopancreatic tissue from P. clarkii that were either challenged with A. hydrophila or treated with PBS. After assembly and annotation of the data, 32,041 unigenes with an average length of 1512 base pairs were identified. Compared to control group, Differential gene expression (DEG) analysis revealed 608 DEGs were obtained, of which 274 unigenes were upregulated and 334 were downregulated in the A. hydrophila group. Furthermore, the expression levels of eight selected immune-related DEGs were validated by qRT-PCR, substantiating the reliability of RNA-seq results. This study not only provides effective data support for immune defense strategies of P. clarkii in response to bacterial infections, but also provides new information about the P. clarkii immune system and defense mechanisms, and a valuable basis for further studies to elucidate the molecular immune mechanisms of this species.

Comparative transcriptomic analysis of the different developmental stages of ovary in red swamp crayfish Procambarus clarkii

Background

The red swamp crayfish Procambarus clarkii is a freshwater species that possesses high adaptability, environmental tolerance, and fecundity. P. clarkii is artificially farmed on a large scale in China. However, the molecular mechanisms of ovarian development in P. clarkii remain largely unknown. In this study, we identified four stages of P. clarkii ovary development, the previtellogenic stage (stage I), early vitellogenic stage (stage II), middle vitellogenic stage (stage III), and mature stage (stage IV) and compared the transcriptomics among these four stages through next-generation sequencing (NGS).

Results

The total numbers of clean reads of the four stages ranged from 42,013,648 to 62,220,956. A total of 216,444 unigenes were obtained, and the GC content of most unigenes was slightly less than the AT content. Principal Component Analysis (PCA) and Anosim analysis demonstrated that the grouping of these four stages was feasible, and each stage could be distinguished from the others. In the expression pattern analysis, 2301 genes were continuously increase from stage I to stage IV, and 2660 genes were sharply decrease at stage IV compared to stages I-III. By comparing each of the stages at the same time, four clusters of differentially expressed genes (DEGs) were found to be uniquely highly expressed in stage I (136 genes), stage II (43 genes), stage III-IV (49 genes), and stage IV (22 genes), thus exhibiting developmental stage specificity. Moreover, in comparisons between adjacent stages, the number of DEGs between stage III and IV was the highest. GO enrichment analysis demonstrated that nutrient reservoir activity was highest at stage II and that this played a foreshadowing role in ovarian development, and the GO terms of cell, intracellular and organelle participated in the ovary maturation during later stages. In addition, KEGG pathway analysis revealed that the early development of the ovary was mainly associated with the PI3K-Akt signaling pathway and focal adhesion; the middle developmental period was related to apoptosis, lysine biosynthesis, and the NF-kappa B signaling pathway; the late developmental period was involved with the cell cycle and the p53 signaling pathway.

Conclusion

These transcriptomic data provide insights into the molecular mechanisms of ovarian development in P. clarkii. The results will be helpful for improving the reproduction and development of this aquatic species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07537-x.

Iso-seq analysis and functional annotation of the Santa Fe cave crayfish (Procambarus erythrops) transcriptome

The genus Procambarus represents a diverse genus of freshwater crayfish that includes epigean species, stygobitic species, and at least one parthenogenic species. Despite its evolutionary, ecological, and economic importance, most genomic and transcriptomic resources for this genus are limited to a couple of model species. We sequenced the transcriptome of a non-model species, P. erythrops, a geographically restricted stygobitic species from Florida. RNA isolated from gill, muscle and eye tissue was pooled to create a de novo transcriptome assembly using Single Molecule Real-Time sequencing (PacBio), resulting in 19,442 full-length isoforms. The assembly has been deposited in the NCBI (BioProject PRJNA657230). These data will make an important contribution to the comparative study of transcriptome evolution in crayfish and crustaceans.

Transcriptomic analysis of Procambarus clarkii affected by "Black May" disease

Each year from April to May, high mortality rates are reported in red swamp crayfish (Procambarus clarkii) cultured in Jiangsu and other regions, in China, and this phenomenon has come to be known as “Black May” disease (BMD). Therefore, in order to investigate the possible causes of this disease, this study gathered BMD-affected P. clarkii samples and performed transcriptome analysis on hepatopancreas, gill, and muscle tissues. A total of 19,995,164, 149,212,804, and 222,053,848 clean reads were respectively obtained from the gills, muscle, and hepatopancreas of BMD-affected P. clarkii, and 114,024 unigenes were identified. The number of differentially expressed genes (DEGs) in gill, muscle, and hepatopancreas was 1703, 964, and 476, respectively. GO and KEGG enrichment analyses of the DEGs were then conducted. Based on KEGG pathway enrichment analysis, the most significantly differentially expressed pathways were mainly those involved with metabolism, human disease, and cellular processes. Further analysis of the significantly DEGs revealed that they were mainly related to the mitochondrial-mediated apoptosis pathway and that the expression of these DEGs was mostly down-regulated. Moreover, the expression of genes related to immune and metabolism-related pathways was also significantly down-regulated, and these significantly-inhibited pathways were the likely causes of P. clarkii death. Therefore, our results provide a basis for the identification of BMD causes.

Differential transcriptomic analysis of crayfish (Procambarus clarkii) from a rice coculture system challenged by Vibrio parahaemolyticus

Rice-crayfish (Procambarus clarkii) coculture is an effective farming mode and has been promoted in various regions of China. However, infection in crayfish can be a significant economic drain. We found crayfish infected with Vibrio parahemolyticus (VP), and to understand the molecular mechanisms of the immune responses of crayfish to VP infection, Illumina sequencing was employed to identify changes in the mRNA of hepatopancreatic tissue. A total of 47.30 and 43.01million high-quality transcriptome reads were generated from the hepatopancreatic samples of the experimental group (EG) and control group (CG), respectively. We found 5559 genes were significantly differentially expressed, including 2521 up-regulated genes (45.35%) and 3038 down-regulated genes (54.65%). These genes were enriched in 126 GO terms and 76 KEGG pathways (P ≤ 0.05), including the MAPK and PI3K-Akt signaling pathways and cell adhesion molecules, with 23 up-regulated genes and 3 down-regulated genes related to immune responses in the EG relative to the CG. Histopathological analysis revealed that the epithelial cells of the hepatopancreatic tubules in the EG were severely atrophic, necrotic, and exfoliated, resulting in thin and collapsing hepatopancreatic tubules. The expression patterns of 8 differentially expressed genes involved in immune responses were validated by quantitative real-time RT-PCR. These results provide a valuable basis for the immune responses of crayfish to acute hepatopancreatic necrosis disease at transcriptome level.

Molecular characterization and expression analysis of tumor necrosis factor receptor-associated factor 6 (traf6) like gene involved in antibacterial innate immune of fresh water crayfish, Procambarus clarkii

In invertebrates, innate immunity was the crucial defending pattern against pathogenic microorganisms. For the past few years, Toll or Toll like receptors (TLRs) signaling pathway was studied extensively in crustaceans. Among the components of Toll or Toll like receptors (TLRs) signaling pathway, tumor necrosis factor receptor-associated factor 6 (TRAF6) acted as an important cytoplasmic adaptor, which was conserved from Drosophila to human. In this study, a new traf6 like gene was cloned from hepatopancreas of P. clarkii. After challenged respectively by S. aureus or E. ictaluri, the expression profiles were studied. And the results showed that the mRNA transcript of Pc-traf6 like gene was up-regulated significantly in the hemocytes, hepatopancreas, gills, and intestine of crayfish. After Pc-traf6 like gene was knocked down, the expression levels of transcription factor (Dorsal) and some crucial immunity effectors (ALF 3, Lysozyme 1, Lectin 1, and Crustin 2) in TLRs signaling pathway were dramatically suppressed. Simultaneously, the survival rate of crayfish challenged respectively by S. aureus or E. ictaluri was significantly decreased in RNAi assay. All these results indicated that Pc-traf6 like gene played an important role in regulating the expression of downstream effectors in the TLRs signaling pathway of crayfish.

The feeding of dietary Codonopsis pilosula polysaccharide enhances the immune responses, the expression of immune-related genes and the growth performance of red swamp crayfish (Procambarus clarkii)

Polysaccharides have many functions in aquatic animals and are widely used as immunopotentiators. However, despite the emergence of serious diseases, few studies have explored the effects of Codonopsis pilosula polysaccharide (CPP) on crustaceans. We studied the effects of CPP on the growth performance, nonspecific immunity, antioxidant activity and disease resistance of red swamp crayfish (Procambarus clarkii). Healthy crayfish (5.80 ± 0.1 g) were fed diets supplemented with 0% (control), 0.05%, 0.1%, 0.15%, 0.20%, and 0.30% CPP for 8 weeks. At the end of the 8-week feeding trial, the optimal final body weight (FBW), weight gain (WG), specific growth rate (SGR), and feed conversion ratio (FCR) were observed in the crayfish fed the diets with 0.15% and 0.20% CPP, followed by those fed the diet with 0.30% CPP and then those fed the diet with 0.10% CPP, whereas the values of these parameters were obtained with the control crayfish (P < 0.05). The crayfish fed the diets with 0.15% and 0.20% CPP exhibited a significantly higher total hemocyte count (THC) and significantly increased phenoloxidase (PO), lysozyme (LZM), hemocyte (Hc), acid phosphatase (ACP) and alkaline phosphatase (AKP) compared with those belonging to the other groups (P < 0.05). The crayfish fed the diets with 0.15% and 0.2% CPP exhibited significantly higher total superoxide dismutase (T-SOD) and glutathione peroxidase (GPx) activities, a significantly increased total antioxidant capacity (T-AOC) and a significantly lower malondialdehyde (MDA) content compared with the other groups (P < 0.05), which indicated that antioxidant capacity was significantly induced by the CPP-supplemented diets. Significantly upregulated expression of immune-related genes (anti-lipopolysaccharide factors (alf), peroxiredoxin (prx5), cathepsin B (ctsb), mitochondrial manganese superoxide dismutase (mtMnsod), cyclophilin A (cypa), glutathione peroxidase (gpx), Toll-like receptor 3 (tlr3), and heat shock protein 70 (hsp70)) was detected in the crayfish fed the diets supplemented with 0.15% and 0.20% CPP diet compared with the levels observed in the control crayfish. These results showed that dietary CPP supplementation greatly improved the growth, immunity and antioxidant capacities of crayfish, and according to the observed results, 0.15%-0.2% is the recommended optimal level of CPP dietary supplementation for crayfish.

Effects of hesperidin on the growth performance, antioxidant capacity, immune responses and disease resistance of red swamp crayfish (Procambarus clarkii)

We evaluated the effects of hesperidin on the nonspecific immunity, antioxidant capacity and growth performance of red swamp crayfish (Procambarus clarkii). A total of 900 healthy crayfish were randomly divided into six groups: the control group (fed the basal diet) and the HES25, HES50, HES75, HES100 and HES150 groups, which were fed the basal diet supplemented with 25, 50, 75, 100 and 150 mg kg^-1 hesperidin, respectively. The feeding experiment lasted 8 weeks. The results indicated that compared with the control group, the crayfish groups supplemented with 50-150 mg kg^-1 hesperidin had a decreased feed conversion ratio (FCR) and increased final body weight (FBW), specific growth rate (SGR) and weight gain (WG) (P < 0.05). The protein carbonyl content (PCC), reactive oxygen species (ROS) level and malondialdehyde (MDA) level in the hepatopancreas and hemocytes were significantly lower, while the total antioxidant capacity (T-AOC), glutathione peroxidase (GPx) activity, and superoxide dismutase (SOD) activity were significantly higher in the crayfish groups supplemented with 50-150 mg kg^-1 hesperidin than in the control group. Supplementation with 50-150 mg kg^-1 hesperidin significantly increased the activities of acid phosphatase (ACP), alkaline phosphatase (AKP), lysozyme (LZM), and phenoloxidase (PO) compared with the control group (P < 0.05); upregulated the mRNA expression of cyclophilin A (CypA), extracellular copper-zinc superoxide dismutase (ecCuZnSOD), GPxs, crustin, astacidin, Toll3 and heat shock protein 70 (HSP70) (P < 0.05); and decreased crayfish mortality following white spot syndrome virus (WSSV) infection. These findings indicate that dietary hesperidin supplementation at an optimum dose of 50-150 mg kg^-1 may effectively improve nonspecific immunity, antioxidant capacity and growth performance in crayfish.

A comparative synthesis of transcriptomic analyses reveals major differences between WSSV-susceptible Litopenaeus vannamei and WSSV-refractory Macrobrachium rosenbergii

Since the 1990s White Spot Syndrome Virus (WSSV) has severely affected shrimp aquaculture worldwide causing a global pandemic of White Spot Disease (WSD) in penaeid culture. However, not all decapod species that can be infected by WSSV show the same susceptibility to the virus, thus raising interesting questions regarding the potential genetic traits that might confer resistance to WSSV. In order to shed light into the genetic markers of WSSV resistance, we employed a dual approach: i) we initially analysed the transcriptomes derived from the hepatopancreas of two species, the susceptible white shrimp Litopenaeus vannamei and the refractory fresh water prawn Macrobrachium rosenbergii, both infected with WSSV. We found a large number of differentially expressed genes (DEGs) belonging to the immune system (mostly anti-microbial peptides and haemolymph clotting components) that were generally up-regulated in M. rosenbergii and down-regulated in L. vannamei. Further, in both species we identified many up-regulated DEGs that were related to metabolism (suggesting a metabolic shift during the infection) and, interestingly, in L. vannamei only, we found several DEGs that were related to moult and suggested an inhibition of the moult cycle in this species following WSSV infection. ii) we then identified a limited number of genetic markers putatively linked with WSD tolerance by employing an ecological genomics approach in which we compared published reports with our own RNA-seq datasets for different decapod species infected with WSSV. Using this second comparative approach, we found nine candidate genes which are consistently down-regulated in susceptible species and up-regulated in refractory species and which have a role in immune response. Together our data offer novel insights into gene expression differences that can be found in susceptible and refractory decapod species infected with WSSV and provide a valuable resource towards our understanding of the potential genetic basis of tolerance to WSSV.

Characterization and expression analysis of immune-related genes in the red swamp crayfish, Procambarus clarkii in response to lipopolysaccharide challenge

To learn more about red swamp crayfish related genes in response to bacterial infections, we investigated immune-related genes induced by lipopolysaccharide (LPS) in the hepatopancreas using high-throughput sequencing method. In present the study, a total of 55,107 unigenes were identified, with an average length of 678 bp. A total of 2215 differentially expressed genes (DEGs) were found, including 669 up-regulated genes and 1546 down-regulated genes. The result of Gene ontology (GO) analysis revealed that 3017 DEGs were enriched in 19 biological process subcategories, 17 cellular component subcategories and 15 molecular function subcategories. The top 20 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways showed that "ribosome" was the most abundant group, which had 34 DEGs. KEGG enrichment analysis identified several immune response pathways. Real-time quantitative reverse transcription-PCR (qRT-PCR) results exhibited that several immune responsive genes were greatly up-regulated following LPS stimulation as observed in the results of high-throughput sequencing. Overall, this study provides new insight into the immune defense mechanisms of P. clarkii against LPS infection.

New insight into the molecular basis of Fe (III) stress responses of Procambarus clarkii by transcriptome analysis

Iron in excess can have toxic effects on living organisms. In China, the freshwater crayfish Procambarus clarkii is a source of aquatic food with high-quality protein and has significant commercial value. P. clarkii shows oxidative stress on exposure to heavy metals, and antioxidant enzymes, such as ubiquitination enzymes and proteasomes, play important roles in oxidative stress. To understand the antioxidant defense system of P. clarkii, we analyzed the hepatopancreas transcriptomes of P. clarkii after stimulation with FeCl₃. In total, 5199 differentially expressed genes (DEGs) were identified (2747 upregulated and 2452 downregulated). GO analysis revealed that these DEGs belonged to 16 cellular component, 16 molecular function, and 19 biological process subcategories. A total of 1069 DEGs were classified into 25 categories by using COG. Some antioxidant defense pathways, such as "Ubiquitin mediated proteolysis" and "Glutathione metabolism," were identified using KEGG. In addition, quantitative real time-PCR (qRT-PCR) substantiated the up-regulation of a random selection of DEGs including antioxidant and immune defense genes. We obtained information for P. clarkii transcriptome databases and new insights into the responses of P. clarkii hepatopancreas to heavy metals.

Temporal changes in transcriptome profile provide insights of White Spot Syndrome Virus infection in Litopenaeus vannamei

Shrimp aquaculture is severely affected by WSSV. Despite an increasing effort to understand host/virus interaction by characterizing changes in gene expression (GE) following WSSV infection, the majority of published studies have focussed on a single time-point, providing limited insight on the development of host-pathogen interaction over the infection cycle. Using RNA-seq, we contrasted GE in gills of Litopenaeus vannamei at 1.5, 18 and 56 hours-post-infection (hpi), between WSSV-challenged and control shrimps. Time course analysis revealed 5097 differentially expressed genes: 63 DEGs were viral genes and their expression in WSSV group either peaked at 18 hpi (and decreased at 56 hpi) or increased linearly up to 56 hpi, suggesting a different role played by these genes during the course of infection. The remaining DEGs showed that WSSV altered the expression of metabolic, immune, apoptotic and cytoskeletal genes and was able to inhibit NF-κB and JAK/STAT pathways. Interestingly, GE changes were not consistent through the course of infection but were dynamic with time, suggesting the complexity of host-pathogen interaction. These data offer novel insights into the cellular functions that are affected during the course of infection and ultimately provide a valuable resource towards our understanding of the host-pathogen dynamics and its variation with time.

An E3 ubiquitin ligase TRIM9 is involved in WSSV infection via interaction with β-TrCP

The TRIpartite Motif (TRIM) proteins are known to play key roles in cell differentiation, apoptosis, development, autophagy and innate immunity. In the present study, a TRIM9 homolog (named LvTRIM9) was identified from the transcriptome of the Pacific whiteleg shrimp Litopenaeus vannamei. The deduced amino acid sequence of LvTRIM9 possessed typical features of TRIMs, consisting of a RING domain, two B-boxes, a coiled-coil domain, a FN3 domain, and a SPRY domain. The transcript of LvTRIM9 was detected in most tissues of the shrimp. Its expression level was obviously up-regulated at 3, 12 and 24 h post white spot syndrome virus (WSSV) infection. Knockdown of LvTRIM9 gene expression by double-strand RNA mediated interference could lead to a decrease of virus copy number in WSSV-infected shrimp. Yeast two-hybrid analysis showed that LvTRIM9 could directly interact with beta-transducin repeat-containing protein of shrimp (Lvβ-TrCP), an inhibitor of NF-κB pathway. Meanwhile, knockdown of LvTRIM9 could also up-regulate the expression levels of LvRelish and downstream production of antimicrobial peptides in the intestine of shrimp. These data indicated that WSSV might hijack the LvTRIM9 for its propagation through inhibition of NF-κB pathway and downstream antimicrobial peptides production via interaction of LvTRIM9 with Lvβ-TrCP in shrimp. The study improved our understanding about the impact of E3 ubiquitin ligases on the innate immune signaling pathway of shrimp and its role during WSSV infection.

High-throughput sequencing analysis of microRNAs in gills of red swamp crayfish, Procambarus clarkii infected with white spot syndrome virus

MicroRNAs (miRNAs) are important posttranscriptional regulators. They play an important role in the antiviral innate immunity of invertebrates. In the present study, high-throughput small RNAs Illumina sequencing systems were carried out to identify differentially expressed miRNAs (DEMs) in the gills of Procambarus clarkii, which was challenged with white spot syndrome virus (WSSV). Our results identified 11,617 known and 6 novel miRNAs in normal group (NG) and WSSV-challenged group (WG) small RNA libraries. Additionally, 27 DEMs were shown to participate in the antiviral innate immunity of P. clarkii and were significantly upregulated or downregulated. In addition, the results of the KEGG pathway prediction of the DEMs target genes showed that putative target genes of these 27 DEMs were related mainly to the RNA transport pathway, tight junction pathway, mRNA surveillance pathway, regulation actin cytoskeleton pathway, focal adhesion pathway, and MAPK signaling pathway. These results provide important information for future studies about the antiviral innate immunity of crustaceans.

Toxicity of Melaleuca alternifolia essential oil to the mitochondrion and NAD+/NADH dehydrogenase in Tribolium confusum

BACKGROUND

In our previous study, Melaleuca alternifolia essential oil (EO) was considered to have an insecticidal effect by acting on the mitochondrial respiratory chain in insects. However, the mode of action is not fully understood.

METHODS

In this study, we investigated the insecticidal efficacy of the M. alternifolia EO against another major stored-product pest, Tribolium confusum Jacquelin du Val. Rarefaction and vacuolization of the mitochondrial matrix were evident in oil-fumigated T. confusum adults.

RESULTS

Alterations to the mitochondria confirmed the insecticidal effect of the M. alternifolia EO. Furthermore, comparative transcriptome analysis of T. confusum using RNA-seq indicated that most of the differentially expressed genes were involved in insecticide detoxification and mitochondrial function. The biochemical analysis showed that the intracellular NAD+/NADH ratio is involved in the differential effect of the M. alternifolia EO.

DISCUSSION

These results led us to conclude that NAD+/NADH dehydrogenase may be the prime target site for the M. alternifolia EO in insects, leading to blocking of the mitochondrial respiratory chain.

Transcriptome analysis of Macrobrachium rosenbergii intestines under the white spot syndrome virus and poly (I:C) challenges

Intestine is a primary site of the white spot syndrome virus (WSSV) infection in most crustaceans. To date, little is known about its role in the anti-viral immune response in the freshwater prawn Macrobrachium rosenbergii. In this study, next-generation sequencing was employed to investigate the M. rosenbergii intestine transcriptomes following WSSV or poly I:C challenges. A total of 41.06 M, 39.58 M and 47.00 M clean reads were generated and assembled into 65,340, 71,241 and 70,614 transcripts from the negative control group (NG), WSSV challenge group (WG) and poly I:C treatment group (PG) respectively. Based on homology searches, functional annotation with 7 databases (NR, NT, GO, COG, KEGG, Swissprot and Interpro) for 88,412 transcripts was performed. After WSSV or poly (I:C) challenge, the numbers of up-regulated differentially expressed genes (DEGs) were greater than the down-regulated DEGs. Gene Ontology (GO) classification of the DEGs also distributed similarly, with the same top 10 annotations and were all assigned to the signaling pathways, including spliceosome, Rap1 signaling pathway, proteoglycans, PI3K-Akt signaling pathway, ECM receptor interaction. Results could contribute to a better understanding of the intestinal immune response to viral pathogens.

Characterization of Population Genetic Structure of red swamp crayfish, Procambarus clarkii, in China

The red swamp crayfish (Procambarus clarkii) is one of the most economically important farmed aquatic species in China. However, it is also a famous invasive species in the world. This invasive species was dispersed most via human activities including intentional or unintentional carry in China. Thus, P. clarkii naturally distributed in China provides us a desirable mode to investigate the genetic structure of an invasive species dispersed mainly by human-mediated factors. To reveal the impact of human-mediated dispersal on genetic structure of P. clarkii in China, a total of 22,043 genome-wide SNPs were obtained from approximately 7.4 billion raw reads using 2b-RAD technique in this study. An evident pattern of population genetic structure and the asymmetrical migrational rates between different regions were observed with 22 populations based on these SNPs. This study provide a better understanding of the population genetic structure and demographic history of P. clarkii populations in China, inferring that anthropogenic factors (aquaculture or by accident) and ecological factors (e.g., complicated topography and climatic environment), as well as its special biological traits could account for the current population structure pattern and dispersal history of P. clarkii.

Crayfish immunity - Recent findings

Freshwater crayfish is an important commodity as well as a successful model for studies on crustacean immunity. Due to the ease with which they are kept and the available methods for hemocyte separation and culture they have proven to be very useful. Here, recent progress regarding pattern recognition, immune effector production and antiviral mechanisms are discussed. Several cases of functional resemblance between vertebrate complement and the crayfish immune reactions are highlighted.

Transcriptome analysis of hepatopancraes of Procambarus clarkii challenged with polyriboinosinic polyribocytidylic acid (poly I:C)

Crustacean hepatopancreas regulates metabolic processes, biogenesis and innate immune processes, and the knowledge on its immune genes are crucial to understand antimicrobial mechanisms. In this study, we reported the transcriptomic profile of Procambarus clarkii hepatopancreas after poly I:C administration using high-throughput sequencing. Following de novo assembly 56,716 unigene sequences with an average length of 810 bp was obtained. The unigene sequences were annotated to three ontologies including cellular components, biological processes and molecular functions, further 56,716 unigene sequences were mapped to 25 COG categories. A total of 2497 differentially expressed genes (DEGs) were identified following the comparative analysis between poly I:C treated and control group, and then KEGG enrichment analysis were performed to detect immune related pathways. Quantitative real time polymerase chain reaction showed that the selected DEGs significantly up-regulated following poly I:C administration in comparison to control group. The transcriptomic sequence information will improve the knowledge of this economically important crustacean, and will shed light on its antiviral immune mechanisms.

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

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

Identification and characterization of intestine microRNAs and targets in red swamp crayfish, Procambarus clarkii infected with white spot syndrome virus

MicroRNAs (miRNAs) are small non-coding endogenous RNA molecules that play important roles in the innate immunity system of invertebrates, especially in the aspect of antivirus. In the present study, high-throughput small RNA Illumina sequencing systems were used to identify differentially expressed miRNAs (DEMs) from the intestines of Procambarus clarkii that were infected with white spot syndrome virus (WSSV). As a result, 39 known and 12 novel miRNAs were identified in both NG and WG small RNA libraries. Seven DEMs were determined to be involved in the antiviral innate immunity in the intestines of P. clarkii. The results of the target gene predictions of the DEMs showed that the putative target genes of these 7 DEMs are related to tight junctions, vascular smooth muscle contraction regulation of the actin cytoskeleton, focal adhesion, RNA transport, mRNA surveillance, viral carcinogenesis, and Salmonella infection. These results provide theoretical insights for future studies on the antiviral immunity of crustaceans.

Experimental inoculation of oriental river prawn Macrobrachium nipponense with white spot syndrome virus (WSSV)

The oriental river prawn Macrobrachium nipponense is an economically important species that is widely farmed in China. White spot syndrome virus (WSSV) is one of the most devastating pathogens of the cultured shrimp Litopenaeus vannamei, responsible for massive loss of its commercial products worldwide. We investigated the infectivity and pathogenicity of WSSV in adult M. nipponense using standardized conditions for L. vannamei. The median lethal dose of WSSV in adult M. nipponense was 103.84±0.06 copies g-1, which was about 1000-fold higher than in L. vannamei (100.59±0.22 copies g-1). WSSV was detected by 2-step PCR in the gills, hepatopancreas, muscle, stomach, heart, gut, nerve, integument, pereopod, eyestalk, testis, and ovary of experimentally infected dead M. nipponense. Lesions were observed histologically following WSSV injection, showing basophilic intranuclear inclusion bodies in the hepatopancreas and subsequently in the gills. The clearance of WSSV was observed in hepatopancreas and gills at 48 and 96 h post-inoculation, respectively. No histological lesions were detected in muscle from 0-96 h post-injection. The results show that the oriental river prawn M. nipponense can be infected by WSSV and the infections are self limiting over time; therefore, M. nipponense may serve as a useful model for studying resistance to WSSV.

Identification and characterization of lymph organ microRNAs in red swamp crayfish, Procambarus clarkii infected with white spot syndrome virus

MicroRNAs (miRNAs) were important post-transcriptional regulators and played vital roles in innate immunity system of invertebrates, especially in the aspect of antivirus. In this study, using high-throughput small RNAs Illumina sequencing system, differentially expressed miRNAs (DEMs) from lymph organs in red swamp crayfish, Procambarus clarkii, infected with white spot syndrome virus, were identified. As a result, 32 known miRNAs and 7 novel miRNAs were identified in crayfish lymph organ small RNAs library of NG and WG. Among them, 7 differentially expressed miRNAs (DEMs) were predicted to be involved in the lymph organ antiviral innate immunity of P. clarkii. Besides, the results showed that putative target genes of these DEMs were related with tight junction, RNA transport, regulation of actin cytoskeleton, focal adhesion, vascular smooth muscle contraction, mRNA surveillance pathway, NOD-like receptor signaling pathway, leukocyte transendothelial migration, and protein processing in endoplasmic reticulum. These results might provide the guiding theoretical foundation for future studies about crustaceans' antiviral innate immunity.

Transcriptome profiling of the Macrobrachium rosenbergii lymphoid organ under the white spot syndrome virus challenge

Macrobrachium rosenbergii is a crustacean with economic importance, and adult prawns are generally thought to be tolerant to white spot syndrome virus (WSSV) infection. Although certain genes are known to respond to WSSV infection and lymphoid tissue is an important immune organ, the response of lymphoid organ to WSSV infection is unclear. Next-generation sequencing was employed in this study to determine the transcriptome differences between WSSV infection and mock lymphoid organs. A total of 44,606,694 and 40,384,856 clean reads were generated and assembled into 73,658 and 72,374 unigenes from the control sample and the WSSV infection sample, respectively. Based on homology searches, KEGG, GO, and COG analysis, 21,323 unigenes were annotated. Among them, 4951 differential expression genes were identified and categorized into 244 metabolic pathways. Coagulation cascades, and pattern recognition receptor signaling pathways were used as examples to discuss the response of host to WSSV infection. We also identified 12,308 simple sequence repeats, which can be further used as functional markers. Results contribute to a better understanding of the immune response of prawn lymphoid organ to WSSV and provide information for identifying novel genes in the absence of the prawn genome.

Identification of ecdysteroid receptor-mediated signaling pathways in the hepatopancreas of the red swamp crayfish, Procambarus clarkii

The hepatopancreas of crustaceans plays an important role in lipid and carbohydrate metabolism, digestion of food, and biogenesis. In this study, the hepatopancreas transcriptome from the red crayfish Procambarus clarkii was characterized for the first time using high-throughput sequencing, producing approximately 41.4 million reads were obtained. After de novo assembly, 57,363 unigenes with an average length of 725bp were identified, Gene Ontology analysis categorized 22,580 as being involved in biological processes, among which metabolic process and cellular process groups were the most highly enriched. A total of 8034 unigenes were assigned to 223 metabolic pathways following mapping against the Kyoto encyclopedia of genes and genomes (KEGG) database. Ecdysteroid receptor (EcR)-mediated signaling pathways were investigated using digital gene expression (DGE) analysis following RNA interference targeting the EcR. A total of 529 differentially expressed genes (DEGs) were identified, including 322 downregulated and 207 upregulated unigenes. Of these, 445 (84.12%) were annotated successfully by alignment with known sequences, many of which were related to catalytic activity and binding functional categories. Using KEGG enrichment analysis, 183 DEGs were clustered into 78 pathways, and six significantly enriched pathways were predicted. The expression patterns of candidate genes identified by real-time PCR were consistent with the DGE results.

Insecticidal Activity of Melaleuca alternifolia Essential Oil and RNA-Seq Analysis of Sitophilus zeamais Transcriptome in Response to Oil Fumigation

BACKGROUND

The cereal weevil, Sitophilus zeamais is one of the most destructive pests of stored cereals worldwide. Frequent use of fumigants for managing stored-product insects has led to the development of resistance in insects. Essential oils from aromatic plants including the tea oil plant, Melaleuca alternifolia may provide environmentally friendly alternatives to currently used pest control agents. However, little is known about molecular events involved in stored-product insects in response to plant essential oil fumigation.

RESULTS

M. alternifolia essential oil was shown to possess the fumigant toxicity against S. zeamais. The constituent, terpinen-4-ol was the most effective compound for fumigant toxicity. M. alternifolia essential oil significantly inhibited the activity of three enzymes in S. zeamais, including two detoxifying enzymes, glutathione S-transferase (GST), and carboxylesterase (CarE), as well as a nerve conduction enzyme, acetylcholinesterase (AChE). Comparative transcriptome analysis of S. zeamais through RNA-Seq identified a total of 3,562 differentially expressed genes (DEGs), of which 2,836 and 726 were up-regulated and down-regulated in response to M. alternifolia essential oil fumigation, respectively. Based on gene ontology (GO) analysis, the majority of DEGs were involved in insecticide detoxification and mitochondrial function. Furthermore, an abundance of DEGs mapped into the metabolism pathway in the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway database were associated with respiration and metabolism of xenobiotics, including cytochrome P450s, CarEs, GSTs, and ATP-binding cassette transporters (ABC transporters). Some DEGs mapped into the proteasome and phagosome pathway were found to be significantly enriched. These results led us to propose a model of insecticide action that M. alternifolia essential oil likely directly affects the hydrogen carrier to block the electron flow and interfere energy synthesis in mitochondrial respiratory chain.

CONCLUSION

This is the first study to perform a comparative transcriptome analysis of S. zeamais in response to M. alternifolia essential oil fumigation. Our results provide new insights into the insecticidal mechanism of M. alternifolia essential oil fumigation against S. zeamais and eventually contribute to the management of this important agricultural pest.