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

In silico identification and functional study of long non-coding RNA involved in acute hepatopancreatic necrosis disease caused by Vibrio parahaemolyticus infection in white shrimp, Litopenaeus vannamei

Acute hepatopancreatic necrosis disease (AHPND) caused by toxin-producing Vibrio parahaemolyticus (VpAHPND) has severely affected shrimp production. Long non-coding RNA (lncRNA), a regulatory non-coding RNA, which can play important function in shrimp disease responses. This study aimed to identify and investigate the role of lncRNA involved in VpAHPND infection in Pacific white shrimp, Litopenaeus vannamei. From a total of 368,736 de novo assembled transcripts, 67,559 were identified as putative lncRNAs, and only 72 putative lncRNAs showed differential expression between VpAHPND-infected and normal shrimp. The six candidate lncRNAs were validated for their expression profiles during VpAHPND infection and tissue distribution using RT-qPCR. The role of lnc2088 in response to VpAHPND infection was investigated through RNA interference. The result indicated that the suppression of lnc2088 expression led to an increase in shrimp mortality after VpAHPND infection. To explore the set of genes involved in lnc2088 knockdown, RNA sequencing was performed. A total of 275 differentially expressed transcripts were identified in the hepatopancreas of lnc2088 knockdown shrimp. The expression profiles of five candidate metabolic and immune-related genes were validated in lnc2088 knockdown and VpAHPND-infected shrimp. The result showed that the expression of ChiNAG was significantly increased, while that of NCBP1, WIPF2, and NFKB1 was significantly downregulated in ds2088-injected shrimp. Additionally, the expression of NFKB1, NCBP1 and WIPF2 was significantly increased, whereas that of ChiNAG and CUL5 were significantly decreased after infection with VpAHPND. Our work identified putative lncRNA profiles in L. vannamei in response to VpAHPND infection and investigated the role of lncRNA in shrimp immunity.

Metagenomic dissection of the intestinal microbiome in the giant river prawn Macrobrachium rosenbergii infected with Decapod iridescent virus 1

Microbiome in the intestines of aquatic invertebrates plays pivotal roles in maintaining intestinal homeostasis, especially when the host is exposed to pathogen invasion. Decapod iridescent virus 1 (DIV1) is a devastating virus seriously affecting the productivity and success of crustacean aquaculture. In this study, a metagenomic analysis was conducted to investigate the genomic sequences, community structure and functional characteristics of the intestinal microbiome in the giant river prawn Macrobrachiumrosenbergii infected with DIV1. The results showed that DIV1 infection could significantly reduce the diversity and richness of intestinal microbiome. Proteobacteria represented the largest taxon at the phylum level, and at the species level, the abundance of Gonapodya prolifera and Solemya velum gill symbiont increased significantly following DIV1 infection. In the infected prawns, four metabolic pathways related to purine metabolism, pyrimidine metabolism, glycerophospholipid metabolism, and pentose phosphate pathway, and five pathways related to nucleotide excision repair, homologous recombination, mismatch repair, base excision repair, and DNA replication were significantly enriched. Moreover, several immune response related pathways, such as shigellosis, bacterial invasion of epithelial cells, Salmonella infection, and Vibrio cholerae infection were repressed, indicating that secondary infection in M. rosenbergii may be inhibited via the suppression of these immune related pathways. DIV1 infection led to the induction of microbial carbohydrate enzymes such as the glycoside hydrolases (GHs), and reduced the abundance and number of antibiotic-resistant ontologies (AROs). A variety of AROs were identified from the microbiota, and mdtF and lrfA appeared as the dominant genes in the detected AROs. In addition, antibiotic efflux, antibiotic inactivation, and antibiotic target alteration were the main antibiotic resistance mechanisms. Collectively, the data would enable a deeper understanding of the molecular response of intestinal microbiota to DIV1, and offer more insights into its roles in prawn resistance to DIVI infection.

The Supplementation of Berberine in High-Carbohydrate Diets Improves Glucose Metabolism of Tilapia (Oreochromis niloticus) via Transcriptome, Bile Acid Synthesis Gene Expression and Intestinal Flora

Berberine is an alkaloid used to treat diabetes. This experiment aimed to investigate the effects of berberine supplementation in high-carbohydrate diets on the growth performance, glucose metabolism, bile acid synthesis, liver transcriptome, and intestinal flora of Nile tilapia. The six dietary groups were the C group with 29% carbohydrate, the H group with 44% carbohydrate, and the HB1-HB4 groups supplemented with 25, 50, 75, and 100 mg/kg of berberine in group H. The results of the 8-week trial showed that compared to group C, the abundance of Bacteroidetes was increased in group HB2 (p < 0.05). The cholesterol-7α-hydroxylase (CYP7A1) and sterol-27-hydroxylase (CYP27A1) activities were decreased and the expression of FXR was increased in group HB4 (p < 0.05). The pyruvate carboxylase (PC) and phosphoenolpyruvate carboxykinase (PEPCK) activities was decreased in group HB4 (p < 0.05). The liver transcriptome suggests that berberine affects carbohydrate metabolic pathways and primary bile acid synthesis pathways. In summary, berberine affects the glucose metabolism in tilapia by altering the intestinal flora structure, enriching differentially expressed genes (DEGs) in the bile acid pathway to stimulate bile acid production so that it promotes glycolysis and inhibits gluconeogenesis. Therefore, 100 mg/kg of berberine supplementation in high-carbohydrate diets is beneficial to tilapia.

HIF-1-mediated regulation of LDH gene unravels key insights into MCDV-1 pathogenesis in mud crabs Scylla paramamosain

Hypoxia-inducible factors -1 (HIF-1) is a crucial transcription factor that regulates the expression of glycolytic genes. Our previous study proved that the Mud crab dicistrovirus-1 (MCDV-1) can induce aerobic glycolysis that favors viral replication in mud crab Scylla paramamosain. However, the role of HIF-1 on key glycolytic genes during the MCDV-1 infection has not been examined. In this study, the intricate interplay between HIF-1 and the key glycolysis enzyme, lactate dehydrogenase (LDH), was investigated after MCDV-1 infection. The expression of LDH was significant increased after MCDV-1 infection. Additionally, the expression of HIF-1α was upregulated following MCDV-1 infection, potentially attributed to the downregulation of prolyl hydroxylase domains 2 expression. Subsequent examination of the SpLDH promoter identified the presence of hypoxia response elements (HREs), serving as binding sites for HIF-1α. Intriguingly, experimental evidence demonstrated that SpHIF-1α actively promotes SpLDH transcription through these HREs. To further elucidate the functional significance of SpHIF-1α, targeted silencing was employed, resulting in a substantial reduction in SpLDH expression, activity, and lactate concentrations in MCDV-1-infected mud crabs. Notably, SpHIF-1α-silenced mud crabs exhibited higher survival rates and lower viral loads in hepatopancreas tissues following MCDV-1 infection. These results highlight the critical role of SpHIF-1α in MCDV-1 pathogenesis by regulating LDH gene dynamics, providing valuable insights into the molecular mechanisms underlying the virus-host interaction.

Transcriptome analysis of the hepatopancreas from the Litopenaeus vannamei infected with different flagellum types of Vibrio alginolyticus strains

Vibrio alginolyticus, one of the prevalently harmful Vibrio species found in the ocean, causes significant economic damage in the shrimp farming industry. Its flagellum serves as a crucial virulence factor in the invasion of host organisms. However, the processes of bacteria flagella recognition and activation of the downstream immune system in shrimp remain unclear. To enhance comprehension of this, a ΔflhG strain was created by in-frame deletion of the flhG gene in V. alginolyticus strain HN08155. Then we utilized the transcriptome analysis to examine the different immune responses in Litopenaeus vannamei hepatopancreas after being infected with the wild type and the mutant strains. The results showed that the ΔflhG strain, unlike the wild type, lost its ability to regulate flagella numbers negatively and displayed multiple flagella. When infected with the hyperflagella-type strain, the RNA-seq revealed the upregulation of several immune-related genes in the shrimp hepatopancreas. Notably, two C-type lectins (CTLs), namely galactose-specific lectin nattectin and macrophage mannose receptor 1, and the TNF receptor-associated factor (TRAF) 6 gene were upregulated significantly. These findings suggested that C-type lectins were potentially involved in flagella recognition in shrimp and the immune system was activated through the TRAF6 pathway after flagella detection by CTLs.

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

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

A simple sequence repeats marker of disease resistance in shrimp Litopenaeus vannamei and its application in selective breeding

The polymorphism of the simple sequence repeat (SSR) in the 5' untranslated coding region (5'-UTR) of the antiviral gene IRF (LvIRF) has been shown to be implicated in the resistance to viral pathogens in shrimp Litopenaeus vannamei (L. vannamei). In this study, we explored the potential of this (CT)n-SSR marker in disease resistance breeding and the hereditary property of disease resistance traits in offspring. From 2018 to 2021, eight populations were generated through crossbreeding by selecting individuals according to microsatellite genotyping. Our results demonstrated that shrimp with the shorter (CT)n repeat exhibited higher resistance to white spot syndrome virus (WSSV) or Decapod iridescent virus 1 (DIV1); meanwhile, these resistance traits could be inherited in offspring. Interestingly, we observed that the longer (CT)n repeats were associated with bacterial resistance traits. Accordingly, shrimp with longer (CT)n repeats exhibited higher tolerance to Vibrio parahaemolyticus infection. Taken together, these results indicate that the single (CT)n-SSR marker could be used to selective breeding for both resistance to virus and bacteria in shrimps.

Immune responses and histopathological analyses of giant river prawn (Macrobrachium rosenbergii, De Man 1879) challenged with a sub-lethal dose of decapod iridescent virus 1 (DIV1) and chemical control investigation

Decapod iridescent virus 1 (DIV1) is a lethal virus that has a significant influence on the shrimp and prawn culture industries. The mechanism through which infected prawns respond to the DIV1 virus is currently unknown. Here, we examined in detail the clinical signs, histopathology, and humoral, cellular, and immune-related gene responses after a sub-lethal dose of DIV1 during the acute infection period of 0-120 hours post infection (hpi). Interestingly, at the end of the experiment, DIV1-infected prawns had black lesions on several external regions. The DIV1-infected prawns also exhibited few karyopyknotic nuclei in the gills and intestine tissues and exhibited increasing immunological responses, as revealed by significant increases in all examined parameters, including total hemocytes, phagocytosis, lysozyme, and overall bactericidal activity, from 6 to 48 hpi. In addition, between 72 and 120 hpi, all immune response activities of DIV1-infected prawn were impaired compared with those of normal prawns, indicating negative impacts on immunological parameters. A viral load analysis of various tissues by qPCR indicated that hemocytes were the dominant initial viral target tissues, followed by the gills and hepatopancreas. An expression analysis of crucial immune-related genes by qRT‒PCR revealed various expression patterns in response to DIV1 infection; in particular, fold changes in the relative expression of anti-lipopolysaccharide factors (ALFs), prophenoloxidase (proPO), lipopolysaccharide and β-1,3-glucan binding protein (LGBP) were observed. Additionally, five common chemicals, calcium hypochlorite [Ca(OCl)₂] at 16.25-130 ppm, hydrogen peroxide (H₂O₂) at 8.75-70 ppm, povidone iodine (PVP-I) at 3-24 ppm, benzalkonium chloride (BKC) at 20-160 ppm, and formalin at 25-200 ppm, had a significant effect on the killing of DIV1 particles in vitro within 24 h after exposure. These data will be helpful for determining the health status and immune defense mechanisms of giant river prawns during DIV1 infection periods. The study performed the first application of very common disinfectants, and the obtained information will be useful for implementing effective strategies to prevent and control DIV1 infection in both hatchery and grow-out ponds.

Host-microbiota interactions and responses of Metapenaeus ensis infected with decapod iridescent virus 1

Introduction

Decapod iridescent virus 1 (DIV1) has caused severe economic losses in shrimp aquaculture. So far, Researchs on DIV1-infected shrimp have mainly focused on the hemocytes immune response, while studies on the host-intestine microbiota interactions during DIV1 infection have been scarce.

Methods

This study determined the lethal concentration 50 (LC₅₀) of DIV1 to Metapenaeus ensis, preliminarily determining that M. ensis could serve as a susceptible object for DIV1. The interactions and responses between the immune and intestine microbiota of shrimp under DIV1 infection were also investigated.

Results and Discussion

DIV1 infection decreases intestine bacterial diversity and alters the composition of intestine microbiota. Specifically, DIV1 infection decreases the abundance of potentially beneficial bacteria (Bacteroidetes, Firmicutes, and Actinobacteria), and significantly increases the abundance of pathogenic bacteria such as Vibrio and Photobacterium, thereby increasing the risk of secondary bacterial infections. The results of PICRUSt functional prediction showed that altered intestine microbiota induces host metabolism disorders, which could be attributed to the bioenergetic and biosynthetic requirements for DIV1 replication in shrimp. The comparative transcriptomic analysis showed that some metabolic pathways related to host immunity were significantly activated following DIV1 infection, including ncRNA processing and metabolic process, Ascorbate and aldarate metabolism, and Arachidonic acid metabolism. M. ensis may against DIV1 infection by enhancing the expression of some immune-related genes, such as Wnt16, heat shock protein 90 (Hsp90) and C-type lectin 3 (Ctl3). Notably, correlation analysis of intestinal microbial variation with host immunity showed that expansion of pathogenic bacteria (Vibrio and Photobacterium) in DIV1 infection could increased the expression of NF-κB inhibitors cactus-like and Toll interacting protein (Tollip), which may limit the TLR-mediated immune response and ultimately lead to further DIV1 infection.

Significance and Impact of the Study

This study enhances our understanding of the interactions between shrimp immunity and intestinal microbiota. The ultimate goal is to develop novel immune enhancers for shrimp and formulate a safe and effective DIV1 defense strategy.

Identification and characterization of an envelope protein 168L in Cherax quadricarinatus iridovirus (CQIV)

Cherax quadricarinatus iridovirus (CQIV), a new member of family Iridoviridae, mainly infects the shrimps and crayfish with a high mortality rate. Previous gel-based LC-MS/MS study on CQIV has identified 30 structural proteins. In this study, one of the structural proteins, CQIV-168L, was selected for further analysis. RT-PCR and Western-blotting (WB) detection revealed that the transcript and the protein appeared late during infection of C. quadricarinatus cells and that the transcript was blocked by viral DNA replication inhibitor, indicating that CQIV-168L is a late expression gene. The specific antiserum against CQIV-168L was raised and immunofluorescence analysis showed that CQIV-168L was localized in the cytoplasm and associated with virus factories. Western-blotting (WB) assay suggested that CQIV-168L antiserum bound specifically to a 57-kDa protein in both the intact virions and the envelope fraction. As revealed by immunogold labeling, CQIV-168L was a component of the viral envelope. Findings in this work help to further understand the structure and entry mechanism of CQIV.

Changes in the gene expression and gut microbiome to the infection of decapod iridescent virus 1 in Cherax quadricarinatus

As a new emerging viral pathogen, Decapod iridescent virus 1 (DIV1) seriously threatens crustacean farming in recent years. However, limited research progresses have been made on the immune mechanism between host and viral factors in response to DIV1 infection. In the current study, a natural occurrence of DIV1 infection with obvious clinical signs was found in farmed redclaw crayfish Cherax quadricarinatus, and confirmed by nested PCR detection and histopathological examination. Besides, gene expression profiles were analyzed after being challenged with DIV1, and results showed that 27 immune related genes were upregulated compared with the control group. Moreover, the gut microbiota from healthy and DIV1-infected crayfish were investigated by 16S rDNA high-throughput sequencing. Results showed that significant differences in the microbial composition and function were observed after DIV1 challenge. Furthermore, we discovered that changes in gene expression profiles were correlated with microbiota alterations under DIV1 challenge. Taken together, our findings will provide new insights into the immune response mechanism of DIV1 infection in crustaceans.

Transcriptome analysis provides novel insights into the immune mechanisms of Macrobrachium nipponense during molting

Molting is a basic physiological behavior of the Oriental river prawn (Macrobrachium nipponense), however, the gene expression patterns and immune mechanisms during the molting process of Oriental river prawn are unclear. In the current study, the gene expression levels of the hepatopancreas of the Oriental river prawn at different molting stages (pre-molting, Prm; mid-molting, Mm; and post-molting, Pom) were detected by mRNA sequencing. A total of 1721, 551, and 1054 differentially expressed genes (DEGs) were identified between the Prm hepatopancreas (PrmHe) and Mm hepatopancreas (MmHe), MmHe and Pom hepatopancreas (PomHe) and PrmHe and PomHe, respectively. The results showed that a total of 1151 DEGs were annotated into 316 signaling pathways, and the significantly enriched immune-related pathways were "Lysosome", "Hippo signaling pathway", "Apoptosis", "Autophagy-animal", and "Endocytosis". The qRT-PCR verification results of 30 randomly selected DEGs were consistent with RNA-seq. The expression patterns of eight immune related genes in different molting stages of the Oriental river prawn were analyzed by qRT-PCR. The function of Caspase-1 (CASP1) was further investigated by bioinformatics, qRT-PCR, and RNAi analysis. CASP1 has two identical conserved domains: histidine active site and pentapeptide motif, and the expression of CASP1 is the highest in ovary. The expression levels of triosephosphate isomerase (TPI), Cathepsin B (CTSB) and Hexokinase (HXK) were evaluated after knockdown of CASP1. This research provides a valuable basis to improve our understanding the immune mechanisms of Oriental river prawns at different molting stages. The identification of immune-related genes is of great significance for enhancing the immunity of the Oriental river prawn, or other crustaceans, by transgenic methods in the future.

Transcriptomic Signature of Horseshoe Crab Carcinoscorpius rotundicauda Hemocytes' Response to Lipopolysaccharides

Carcinoscorpius rotundicauda (C. rotundicauda) is one of the four species of horseshoe crabs (HSCs). The HSC hemocytes store defense molecules that are released upon encountering invading pathogens. The HSCs rely on this innate immunity to continue its existence as a living fossil for more than 480 million years. To gain insight into the innate mechanisms involved, transcriptomic analysis was performed on isolated C. rotundicauda hemocytes challenged with lipopolysaccharides (LPS), the main components of the outer cell membrane of gram-negative bacteria. RNA-sequencing with Illumina HiSeq platform resulted in 232,628,086 and 245,448,176 raw reads corresponding to 190,326,253 and 201,180,020 high-quality mappable reads from control and LPS-stimulated hemocytes, respectively. Following LPS-stimulation, 79 genes were significantly upregulated and 265 genes were downregulated. The differentially expressed genes (DEGs) were related to multiple immune functional categories and pathways such as those of the cytoskeleton, Toll and Imd, apoptosis, MAP kinase (MAPK), inositol phosphate metabolism, phagosome, leucocyte endothelial migration, and gram-negative bacterial infection, among others. This study provides important information about the mechanisms of response to LPS, which is relevant for the understanding the HSCs' immune response.

Integrated analysis of intestinal microbiota and metabolomic reveals that decapod iridescent virus 1 (DIV1) infection induces secondary bacterial infection and metabolic reprogramming in Marsupenaeus japonicus

In recent years, with global warming and increasing marine pollution, some novel marine viruses have become widespread in the aquaculture industry, causing huge losses to the aquaculture industry. Decapod iridescent virus 1 (DIV1) is one of the newly discovered marine viruses that has been reported to be detected in a variety of farmed crustacean and wild populations. Several previous studies have found that DIV1 can induce Warburg effect-related gene expression. In this study, the effects of DIV1 infection on intestinal health of shrimp were further explored from the aspects of histological, enzymatic activities, microorganisms and metabolites using Marsupenaeus japonicus as the object of study. The results showed that obvious injury in the intestinal mucosa was observed after DIV1 infection, the oxidative and antioxidant capacity of the shrimp intestine was unbalanced, the activity of lysozyme was decreased, and the activities of digestive enzymes were disordered, and secondary bacterial infection was caused. Furthermore, the increased abundance of harmful bacteria, such as Photobacterium and Vibrio, may synergized with DIV1 to promote the Warburg effect and induce metabolic reprogramming, thereby providing material and energy for DIV1 replication. This study is the first to report the changes of intestinal microbiota and metabolites of M. japonicus under DIV1 infection, demonstrating that DIV1 can induce secondary bacterial infection and metabolic reprogramming. Several bacteria and metabolites highly associated with DIV1 infection were screened, which may be leveraged for diagnosis of pathogenic infections or incorporated as exogenous metabolites to enhance immune response.

Transcriptome analysis of IHHNV infection in Penaeus vannamei at different developmental stages

Infectious hypodermal and haematopoietic necrosis virus (IHHNV) is the smallest known virus in shrimp, which causes runt-deformity syndrome (RDS) and leads to huge economic loss every year in penaeid shrimp farming. Previous studies have shown that the juvenile Penaeus vannamei is more susceptible to IHHNV infection than the adults, but the mechanism is still unclear. In order to investigate the mechanism of pathogenic differences in IHHNV infection of P. vannamei at different developmental stages, the juvenile and adult P. vannamei were studied by transcriptome high-throughput sequencing to analyze their response to IHHNV infection. GO and KEGG enrichment were analyzed to search for differentially expressed genes (DEGs) related to immunity, growth and metabolism. The results showed that many immune-related genes of the juvenile and adult P. vannamei responded differently to IHHNV infection. For the adult P. vannamei, the expression of most immune-related genes was significantly up-regulated, which means that a cellular defense response was triggered after IHHNV infection. However, most immune-related genes in juvenile P. vannamei were inhibited, indicating that the immune system of juvenile the P. vannamei is imperfect and makes it to be more susceptible to IHHNV. Similarly, the growth-related genes of P. vannamei were changed during IHHNV infection. For the juvenile P. vannamei, the growth-related genes were significantly down-regulated, which resulted in a growth hormone disorder and prevented the juvenile P. vannamei from growth. In the adult P. vannamei, most molting-related genes were significantly up-regulated, indicating that IHHNV infection leads the adult P. vannamei to early molting to eliminate pathogen in the body. Metabolic process data showed that energy metabolism pathway was affected when P. vannamei infected with IHHNV. The adult P. vannamei infected with IHHNV can cause energetically costly and lead to the disturbance of the metabolism, activate complex immune systems to resist the invasion of pathogens. The results of this study clarified the response mechanism of P. vannamei at different developmental stages to IHHNV infection, which can provide new insights to IHHNV effective control and a reference for the study of sensitive period of different shrimp virus to host infection.

Integrated Analysis of mRNA-Seq and MiRNA-Seq Reveals the Molecular Mechanism of the Intestinal Immune Response in Marsupenaeus japonicus Under Decapod Iridescent Virus 1 Infection

The intestine is not only an important digestive organ but also an important immune organ for shrimp; it plays a key role in maintaining homeostasis. Decapod iridescent virus 1 (DIV1) is a new type of shrimp-lethal virus that has received extensive attention in recent years. To date, most studies of the shrimp intestinal immune response under viral infections have relied on single omics analyses; there is a lack of systematic multi-omics research. In the current study, intestinal mRNA-seq and microRNA (miRNA)-seq analyses of Marsupenaeus japonicus under DIV1 infection were performed. A total of 1,976 differentially expressed genes (DEGs) and 32 differentially expressed miRNAs (DEMs) were identified. Among them, 21 DEMs were negatively correlated with 194 DEGs from a total of 223 correlations. Functional annotation analysis revealed that M. japonicus can regulate glycosaminoglycan biosynthesis (chondroitin sulfate, dermatan sulfate, and keratan sulfate), vitamin metabolism (retinol metabolism and ascorbate and aldarate metabolism), immune pathway activation (Toll and IMD signaling pathways, Wnt signaling pathway, IL-17 signaling pathway, and Hippo signaling pathway), immunity enzyme activity promotion (triose-phosphate isomerase), antimicrobial peptide (AMP) expression, reactive oxygen species (ROS) production, and cell apoptosis through miRNAs to participate in the host’s antiviral immune response, while DIV1 can influence Warburg effect-related pathways (pyruvate metabolism, glycolysis/gluconeogenesis, and citrate cycle), glycosphingolipid biosynthesis-related pathways (glycosphingolipid biosynthesis—globo and isoglobo series and glycosphingolipid biosynthesis—lacto and neolacto series), and the tight junction and adhesion junction of the intestinal mucosal epithelium through the host’s miRNAs and mRNA to promote its own invasion and replication. These results indicate that intestinal miRNAs play important roles in the shrimp immune response against DIV1 infection. This study provides a basis for further study of the shrimp intestinal antiviral immune response and for the formulation of effective new strategies for the prevention and treatment of DIV1 infection.

Interaction between 038R and 125R of Cherax quadricarinatus iridovirus (CQIV) and their effects on virus replication

Iridovirids are a group icosahedral viruses containing linear double-stranded DNA, and mainly infect invertebrates and poikilothermic vertebrates. Cherax quadricarinatus iridovirus (CQIV) is a new species of the family Iridoviridae and can cause high mortality in shrimps. In CQIV genome, there are 25 conserved genes and the putative products are involved in several viral processes. In this study, three core protein including CQIV-032R, CQIV-125R and CQIV-160L were identified to interact with CQIV-038R by yeast two-hybrid (Y2H), and the interaction between CQIV-038R and CQIV-125R was further confirmed by co-immunoprecipitation (Co-IP) assays. In the expression system, EGFP-038R and mCherry-125R were colocalized in the cytoplasm when co-expressed in Sf9 cells. Moreover, silencing the expression of 038R, 125R or both of these two proteins respectively in C. quadricarinatus cells by small interfering RNAs showed significantly inhibit CQIV replication. Collectively, we identified the interaction between 038R and 125R, and demonstrated they are essential for CQIV replication.

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

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

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

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

Growth trait gene analysis of kuruma shrimp (Marsupenaeus japonicus) by transcriptome study

Growth traits are a vital standard for the animal culture industry. The molecular mechanism of growth traits remains poorly understood, especially in aquaculture, which hinders the development of the selective breeding industry. Genomic resources discovered by next-generation sequencing (NGS) have been widely applied in certain species. However, accurate assembly and downstream analysis by NGS are still major challenges for species without reference genomes. In this study, a comparative transcriptome analysis of an economic crustacean species (Marsupenaeus japonicus) between a fast growth group and slow growth group at different stages was performed by SMRT (single molecule real time) and NGS. A high-quality full-length transcriptome (e.g., mean length of unigenes was longer than those unigenes assembled by Illumina clean reads from previous reports, and annotation rate was higher than Illumina sequencing in the same studies) was generated and analyzed. Several differentially expressed genes (DEGs) related to growth were identified and validated by quantitative real-time PCR (qPCR). The results showed that compared with the late stage, more DEGs were identified at the early stage, indicating that the growth-related physiological activity differences between different individuals at the early stage were higher than at the late stage. Moreover, 215 DEGs were shared between the early stage and late stage, and 109 had divergent functions during development. These 109 genes may play an important role in regulating the specific growth rate (SGR) of kuruma shrimp. In addition, twelve growth-related pathways were shared between the two comparative groups. Among these pathways, the fly Hippo signaling pathway and its key gene Mj14-3-3-like were identified for the first time to be involved in growth traits in crustaceans. Further analysis showed that Mj14-3-3-like was significantly downregulated in the fast growth group at the early stage and late stage; its expression level was reduced to its lowest level at the intermolt stage (C), the most important growth stage in shrimp, suggesting that Mj14-3-3-like may inhibit the growth of kuruma shrimp. Our study helps to elucidate the genes involved in the molecular mechanisms governing growth traits in kuruma shrimp, which is valuable for future shrimp developmental research.

The T cell factor, pangolin, from Litopenaeus vannamei play a positive role in the immune responses against white spot syndrome virus infection

As a downstream interactor of β-catenin, Pangolin which is the homologous protein of the T cell factor/lymphoid enhancer factor (TCF/LEF) in vertebrates is less understood in the research field of immunity. In this study, two isoforms of Litopenaeus vannamei Pangolin (LvPangolin1 and LvPangolin2) were identified. Phylogenetic tree analysis revealed that all of the Pangolin proteins from invertebrates were represented the same lineage. The mRNA expression profiles of the LvPangolin1 and LvPangolin2 genes differed across different tissues. The expression of LvPangolin1 and the amount of LvPangolin1and LvPangolin2 combined (LvPangolinComb) were significantly increased in the haemocyte, intestine and gill but reduced in the hepatopancreas after white spot syndrome virus (WSSV) challenge. The inhibition of LvPangolin1 but not LvPangolinComb significantly reduced the survival rates of L. vannamei after WSSV infection, while significantly higher WSSV viral loads in both LvPangolin1-inhibited and LvPangolinComb-inhibited L. vannamei were observed. Knockdown of LvPangolin by RNAi could distinctly decrease the expression of antimicrobial peptide (AMP) genes and their related transcription factors. All of these results indicate that LvPangolin plays a positive role in the response to WSSV infection and that this may be mediated through regulating the immune signalling pathways which control the expression of AMPs with antiviral abilities.

GSK3β Plays a Negative Role During White Spot Syndrome Virus (WSSV) Infection by Regulating NF-κB Activity in Shrimp Litopenaeus vannamei

Glycogen synthase kinase-3 (GSK3), a cytoplasmic serine/threonine-protein kinase involved in a large number of key cellular processes, is a little-known signaling molecule in virus study. In this study, a GSK3 protein which was highly similar to GSK3β homologs from other species in Litopenaeus vannamei (designated as LvGSK3β) was obtained. LvGSK3β was expressed constitutively in the healthy L. vannamei, at the highest level in the intestine and the lowest level in the eyestalk. White spot syndrome virus (WSSV) reduced LvGSK3β expression was in immune tissues including the hemocyte, intestine, gill and hepatopancreas. The inhibition of LvGSK3β resulted in significantly higher survival rates of L. vannamei during WSSV infection than the control group, and significantly lower WSSV viral loads in LvGSK3β-inhibited L. vannamei were observed. Knockdown of LvGSK3β by RNAi resulted in increases in the expression of LvDorsal and several NF-κB driven antimicrobial peptide (AMP) genes (including ALF, PEN and crustin), but a decrease in LvCactus expression. Accordingly, overexpression of LvGSK3β could reduce the promoter activity of LvDorsal and several AMPs, while the promoter activity of LvCactus was increased. Electrophoretic mobility shift assays (EMSA) showed that LvDorsal could bind to the promoter of LvGSK3β. The interaction between LvGSK3β and LvDorsal or LvCactus was confirmed using co-immunoprecipitation (Co-IP) assays. In addition, the expression of LvGSK3β was dramatically reduced by knockdown of LvDorsal. In summary, the results presented in this study indicated that LvGSK3β had a negative effect on L. vannamei by mediating a feedback regulation of the NF-κB pathway when it is infected by WSSV.

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.