Transcriptome analysis of Procambarus clarkii infected with infectious hypodermal and haematopoietic necrosis virus

Parvoviruses of Aquatic Animals

Family Parvoviridae consists of small, non-enveloped viruses with linear, single-stranded DNA genomes of approximately 4-6 kilobases, subdivided into three subfamilies, Parvovirinae, Densovirinae, and Hamaparvovirinae, and unassigned genus Metalloincertoparvovirus. Parvoviruses of aquatic animals infect crustaceans, mollusks, and finfish. This review describes these parvoviruses, which are highly host-specific and associated with mass morbidity and mortality in both farmed and wild aquatic animals. They include Cherax quadricarinatus densovirus (CqDV) in freshwater crayfish in Queensland, Australia; sea star-associated densovirus (SSaDV) in sunflower sea star on the Northeastern Pacific Coast; Clinch densovirus 1 in freshwater mussels in the Clinch River, Virginia, and Tennessee, USA, in subfamily Densovirinae; hepatopancreatic parvovirus (HPV) and infectious hypodermal and hematopoietic necrosis virus (IHHNV) in farmed shrimp worldwide; Syngnathid ichthamaparvovirus 1 in gulf pipefish in the Gulf of Mexico and parts of South America; tilapia parvovirus (TiPV) in farmed tilapia in China, Thailand, and India, in the subfamily Hamaparvovirinae; and Penaeus monodon metallodensovirus (PmMDV) in Vietnamese P. monodon, in unassigned genus Metalloincertoparvovirus. Also included in the family Parvoviridae are novel parvoviruses detected in both diseased and healthy animals using metagenomic sequencing, such as zander parvovirus from zander in Hungary and salmon parvovirus from sockeye salmon smolts in British Columbia, Canada.

Integration of transcriptome and whole-genome re-sequencing analyses reveal growth-related candidate genes in Procambarus clarkii

Growth is a crucial economic trait of all aquaculture species. It is important to explore the molecular regulation on growth, which could help improve the growth rate of species. Mining the growth-related genes is the foundation for revealing its molecular regulation on growth. Presently, the molecular regulation of growth in Procambarus clarkii is not clear, and the study on exploring growth-related genes is limited. In this study, RNA-Seq was used to compare gene expression profiles of the individuals with different growth rates involved in four groups including Big Male (BM), Big Female (BF), Small male (SM), and Small Female (SF) from one P. clarkii family, and the analyses were performed in combination with sex. Meanwhile, whole-genome resequencing data was used to get growth-specific SNP (Single Nucleotide Polymorphism)/InDel (Insertion/Deletion) sites information. Totally, we identified 16,127 genes, of which 9065 were successfully annotated in the GO database. Among these, 1328 DEGs were identified in BM vs. SM, with 357 up-regulated and 971 down-regulated. Additionally, 3507 DEGs were identified in BF vs. SF, with 241 up-regulated and 3266 down-regulated. 96 DEGs were up-regulated and 820 DEGs were down-regulated in Growth-related Group. The expression levels of nine DEGs were validated by RT-qPCR to verify the analysis results of sequencing. 684,040 growth-related SNPs and 182,050 growth-related InDels were obtained after screened. These findings provide candidate growth-related genes and growth-specific SNP/InDel sites for regulation of growth traits in P. clarkii, and new insight into the molecular regulation of P. clarkii growth.

Salvia miltiorrhiza polysaccharide promotes the health of crayfish (Procambarus clarkii) by promoting hemocyte phagocytosis, protecting hepatopancreas and enhancing intestinal barrier function

Plant polysaccharides as immunomodulators are considered one of the effective measures to reduce antibiotic therapy in aquaculture. The immunomodulatory function of Salvia miltiorrhiza polysaccharides (SMP) has been demonstrated and begun to be applied in vertebrates, but its potential effect on crustaceans is unclear. In this study, crayfish (Procambarus clarkii) was fed with 0 %, 0.3 %, 0.7 %, 1.1 %, and 1.5 % SMP for 4 weeks to investigate the effects of SMP on hemocytes phagocytosis, hepatopancreatic function, and intestinal barrier function. The results revealed that hemocyte phagocytic activity was increased in all SMP groups. During the process of hemocytes phagocytic recognition and formation of phagosomes and phagolysosomes, the mRNA expression levels of mas, hem, rab3, ctsb, and lamp-1 were up-regulated mainly in the 0.3 % SMP group. During the clearance phase of phagocytosis, respiratory burst activity, ROS level, T-SOD, CAT, GST, and LZM activities were mainly increased in the 1.5 % SMP group. Hepatopancreas AKP and GOT activity were no significant change in all SMP groups. ACP activity was significantly enhanced in the 1.1 % SMP group. The GPT activity of 0.3-0.7 % SMP group was significantly decreased. The 0.7 % SMP group had the highest intestinal fold height. The highest index values of OTUs, Ace, Chao, and Shannon were in the 0.3 % SMP group. The dietary addition of 0.3 % SMP led to a tendency of increased relative abundance of Firmicutes and Bacteroidota at the phylum level, while the relative abundance of Proteobacteria at the phylum level decreased. In conclusion, dietary SMP could promote crayfish health by enhancing phagocytosis, protecting hepatopancreas and enhancing intestinal barrier function. This study contributes to the theoretical foundation for exploring the potential application of plant polysaccharides in crustaceans.

Effects of Aminomethylphosphonic Acid on the Transcriptome and Metabolome of Red Swamp Crayfish, Procambarus clarkii

Red swamp crayfish, Procambarus clarkii (P. clarkii), is an important model crustacean organism used in many types of research. However, the effects of different doses of aminomethylphosphonic acid (AMAP) on the transcriptome and metabolites of P. clarkii have not been explored. Thus, this study investigated the molecular and metabolic mechanisms activated at the different exposure dosages of AMAP in P. clarkii to provide new insights into the strategies of P. clarkii in response to the high concentrations of AMAP in the environment. In the present study, the P. clarkii were divided into three groups (control group; low-dosage AMAP exposure; high-dosage AMAP exposure), and hepatopancreatic tissue samples were dependently taken from the three groups. The response mechanisms at the different dosages of AMAP were investigated based on the transcriptome and metabolome data of P. clarkii. Differentially expressed genes and differentially abundant metabolites were identified in the distinct AMAP dosage exposure groups. The genes related to ribosome cell components were significantly up-regulated, suggesting that ribosomes play an essential role in responding to AMAP stress. The metabolite taurine, involved in the taurine and hypotaurine metabolism pathway, was significantly down-regulated. P. clarkii may provide feedback to counteract different dosages of AMAP via the upregulation of ribosome-related genes and multiple metabolic pathways. These key genes and metabolites play an important role in the response to AMAP stress to better prepare for survival in high AMAP concentrations.

Analysis of mRNA and MicroRNA Expression Profiles of Nervous Tissues and Reproductive Tissues in Male Procambarus clarkii After Silencing IAG

The insulin-like androgenic gland hormone gene (IAG), primarily expressed in the androgenic gland (AG), plays a crucial role in controlling male sex differentiation and maintaining male secondary sexual characteristics in decapods. In this study, we investigated the mRNA and microRNA expression profiles of male Procambarus clarkii to understand the transcriptomic regulatory mechanism of IAG after the injection of an efficient siRNA (GsiRNA) designed based on IAG. The results revealed that several differentially expressed genes were enriched in reproduction-related pathways, such as the wnt signaling pathway, MAPK signaling pathway, and GnRH signaling pathway. In the testis (Te), the injection of GsiRNA led to the up-regulation of many ovary-related genes and down-regulation of testis-related genes. Moreover, the brain (Br) and abdominal nerve cord (AN) appeared to be involved in the regulation of IAG, with numerous differentially expressed genes found in Br and AN. Notably, the expression of five neuropeptide genes, Crustacean hyperglycemic hormone, pigment-dispersing hormone, red pigment concentrating hormone precursor, corazonin, and gonadotropin-releasing hormone II receptor isoform X1 in Br/AN, was significantly changed. Additionally, three ovary-related miRNAs (miR-263a, miR-263b, miR-133) highly expressed in Te/AG showed significant up-regulation after GsiRNA injection. Furthermore, the long-term interference of GsiRNA was found to inhibit the development of male external sexual characteristics during the juvenile stage and delay it during the adult stage. This research provides valuable insights into the molecular regulatory mechanism and function of IAG in P. clarkii.

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 analysis of the immunomodulatory effects of Salvia miltiorrhiza polysaccharide on hemocyte immune response in Procambarus clarkii

Natural plant polysaccharide as immune modulator is considered an effective strategy for healthy aquaculture to reduce medicine treatment. Salvia miltiorrhiza polysaccharides (SMP) had applications to regulate immune activity and enhance antioxidant in vertebrates, but the potential function has been rarely reported in crustaceans. In this study, the immunological effects of SMP on hemocytes of Procambarus clarkii were analyzed. Results showed that total superoxide dismutase (T-SOD), phenoloxidase (PO) activity and respiratory burst were up-regulated after SMP treatment. After high-throughput sequencing, 2170 differentially expressed genes (DEGs) including 1294 up-regulated and 876 down-regulated genes were identified. KEGG function enrichment analysis indicated that DEGs are involved in crustaceans cellular immune-related signaling pathways, including lysosome, phagosome and endocytosis. Transcriptome mining and qRT-PCR showed that SMP up-regulated humoral immunity factors gene expression. Diets supplemented with 0.8% SMP significantly up-regulated the total number of hemocytes (THC), T-SOD and PO activity, improved the survival of crayfish after Citrobacter freundii infection. This study suggested that SMP could improve the cellular and humoral immunity of P. clarkii. Furthermore, this finding supplied a molecular foundation for further comprehending the immunopotentiator effects of plant polysaccharides in crustaceans.

Transcriptome and 16S rRNA analysis revealed the response of largemouth bass (Micropterus salmoides) to Rhabdovirus infection

To better understand the response of largemouth bass (Micropterus salmoides) to Micropterus salmoides rhabdovirus (MSRV) infection, we investigated the intestinal bacterial flora and transcriptome profile of fish at 72 hours post-infection (hpi). Total of 1574 differentially expressed genes (DEGs) were identified in largemouth bass spleen following MSRV infection, including 573 upregulated and 1001 downregulated genes. KEGG and GO enrichment analysis revealed that upregulated genes were enriched in certain antiviral related signaling pathway, including NOD-like receptor (NLR), RIG-I like receptors (RLR) and regulation of the interferon (IFN)-γ-mediated signaling pathway, whereas some immune-related DEGs enriched in focal adhesion (FA) and ECM-receptor interaction(ECM-RI) were downregulated, as well as genes associated with metabolic processes, such as peroxisome proliferator-activated receptors (PPAR), adipocytokine signaling pathway, Glycerolipid and Retinol metabolism. Furthermore, the principal component analysis (PCA) and phylogenetic analysis revealed that MSRV infection significantly affected the microbiota of largemouth bass intestine; the LEfSe analysis showed that relative abundances of Streptococcus were significantly increased, while the content of Akkermansia, Enterococcus and Lactobacillus were remarkably decreased in the fish intestine following MSRV infection. Additionally, a high correlation was determined between the expressions of interferon-related upregulated genes and the relative abundance of Streptococcus by redundancy analysis (RDA). These results collectively illustrated that intestinal microbiota composition might be associated with the immune-related gene expression in largemouth bass in response to MSRV infection.

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.

The Role of Intestinal Microbial Metabolites in the Immunity of Equine Animals Infected With Horse Botflies

The microbiota and its metabolites play an important role in regulating the host metabolism and immunity. However, the underlying mechanism is still not well studied. Thus, we conducted the LC-MS/MS analysis and RNA-seq analysis on Equus przewalskii with and without horse botfly infestation to determine the metabolites produced by intestinal microbiota in feces and differentially expressed genes (DEGs) related to the immune response in blood and attempted to link them together. The results showed that parasite infection could change the composition of microbial metabolites. These identified metabolites could be divided into six categories, including compounds with biological roles, bioactive peptides, endocrine-disrupting compounds, pesticides, phytochemical compounds, and lipids. The three pathways involving most metabolites were lipid metabolism, amino acid metabolism, and biosynthesis of other secondary metabolites. The significant differences between the host with and without parasites were shown in 31 metabolites with known functions, which were related to physiological activities of the host. For the gene analysis, we found that parasite infection could alarm the host immune response. The gene of “cathepsin W” involved in innate and adaptive immune responses was upregulated. The two genes of the following functions were downregulated: “protein S100-A8” and “protein S100-A9-like isoform X2” involved in chemokine and cytokine production, the toll-like receptor signaling pathway, and immune and inflammatory responses. GO and KEGG analyses showed that immune-related functions of defense response and Th17 cell differentiation had significant differences between the host with and without parasites, respectively. Last, the relationship between metabolites and genes was determined in this study. The purine metabolism and pyrimidine metabolism contained the most altered metabolites and DEGs, which mainly influenced the conversion of ATP, ADP, AMP, GTP, GMP, GDP, UTP, UDP, UMP, dTTP, dTDP, dTMP, and RNA. Thus, it could be concluded that parasitic infection can change the intestinal microbial metabolic activity and enhance immune response of the host through the pathway of purine and pyrimidine metabolism. This results will be a valuable contribution to understanding the bidirectional association of the parasite, intestinal microbiota, and host.

Integrative omics analysis highlights the immunomodulatory effects of the parasitic dinoflagellate hhematodinium on crustacean hemocytes

Parasitic dinoflagellates in genus Hematodinium have caused substantial economic losses to multiple commercially valuable marine crustaceans around the world. Recent efforts to better understand the life cycle and biology of the parasite have improved our understanding of the disease ecology. However, studies on the host-parasite interaction, especially how Hematodinium parasites evade the host immune response are lacking. To address this shortfall, we used the comprehensive omics approaches (miRNA transcriptomics, iTRAQ-based proteomics) to get insights into the host-parasite interaction between hemocytes from Portunus trituberculatus and Hematodinium perezi in the present study. The parasitic dinoflagellate H. perezi remodeled the miRNome and proteome of hemocytes from challenged hosts, modulated the host immune response at both post-transcriptional and translational levels and caused post-transcriptional regulation to the host immune response. Multiple important cellular and humoral immune-related pathways (ex. Apoptosis, Endocytosis, ECM-receptor interaction, proPO activation pathway, Toll-like signaling pathway, Jak-STAT signaling pathway) were significantly affected by Hematodinium parasites. Through modulation of the host miRNome, the host immune responses of nodulation, proPO activation and antimicrobial peptides were significantly suppressed. Cellular homeostasis was imbalanced via post-transcriptional dysregulation of the phagosome and peroxisome pathways. Cellular structure and communication was seriously impacted by post-transcriptional downregulation of ECM-receptor interaction and focal adhesion pathways. In conclusion, H. perezi parasites could trigger striking changes in the miRNome and proteome of crustacean hemocytes, and this parasite exhibited multifaceted immunomodulatory effects and potential immune-suppressive mechanisms in crustacean hosts.

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.

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.

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.