Hepatopancreas transcriptome analysis of Chinese mitten crab (Eriocheir sinensis) with white hepatopancreas syndrome

Metabolomic and transcriptomic profiling reveals the effect of dietary protein and lipid levels on growth performance in loach (Paramisgurnus dabryanus)

The subject of this study was to explore the optimum requirements of loach (Paramisgurnus dabryanus) regarding dietary proteins and lipids and discuss the underlying mechanism. We designed nine diets to determine the effects of different levels of dietary crude protein (CP: 30%, 35%, and 40%) and ether extract (EE: 6%, 10%, and 14%) on the growth performance and metabolism of P. dabryanus. In total, 2160 healthy P. dabryanus (5.19 ± 0.01 g) were divided into nine groups with four replications at 60 fish per barrel stocking density. The trial lasted for eight weeks. Serum and liver samples were gathered for metabolomic and transcriptomic analyses. The results showed that the specific growth rate of P. dabryanus in the CP40EE10 group was the fastest and notably higher than that in other groups (P< 0.05). Analysis of the metabolome results found that the mTOR signaling pathway, glycerophospholipid metabolism, D-arginine and D-ornithine metabolism were significantly enriched pathways in the CP40EE10 group compared with the other groups (P< 0.05). Moreover, the transcriptomic analysis of differentially expressed genes (DEGs) showed that the expression of ARG (arginase) involved in protein synthesis was significantly upregulated in the CP40EE10 group compared to the slowest growing group (P< 0.05). Additionally, the expression of SPLA2 (secretory phospholipase A2) involved in lipid metabolism and FBP (fructose-1,6-bisphosphatase) involved in glucose metabolism were all significantly downregulated in the CP30EE6 group compared with the CP40EE10 group (P< 0.05). Furthermore, the analysis of differentially expressed metabolites (DEMs) and DEGs co-enriched in the KEGG pathway revealed that the significantly enriched pathways were arginine and proline metabolism, glycerophospholipid metabolism, and glycolysis/gluconeogenesis in CP40EE10 compared with other groups (P< 0.05). We conclude that including 40% CP and 10% EE in the P. dabryanus diet could result in a better growth rate. We hypothesized from metabolomic and transcriptomic analyses that the CP40EE10 diet might promote the growth of P. dabryanus by promoting protein synthesis, lipid metabolism, and energy production.

Comparative transcriptome analysis reveals the effects of different feeding times on the hepatopancreas of Chinese mitten crabs

Feeding rhythms affect a range of physiological functions in crustaceans. To investigate their effect on the physiological functions of Eriocheir sinensis, herein, we analyzed the influence of different feeding times on the hepatopancreas transcriptome via high-throughput sequencing. We harvested the hepatopancreas of crabs at 12:00 on day 11 of the experiment. We weighted the crabs before and after the experiment and found that those in the 06:00 group had the highest weight gain rate. In addition, 512 differentially expressed genes (DEGs) were grouped into nine distinct clusters. Functional enrichment analysis of DEGs showed that E. sinensis metabolic and immune processes were affected by the feeding time. Furthermore, we mapped the DEGs involved in retinol metabolism and the lysosome pathway. To our knowledge, this is the first comparative transcriptomic analysis of the hepatopancreas of E. sinensis based on different feeding times, which provides multi-level information to reveal the mechanism underlying the regulation of feeding rhythms in E. sinensis.

Alternative splicing derived invertebrate variable lymphocyte receptor displays diversity and specificity in immune system of crab Eriocheir sinensis

Variable lymphocyte receptors (VLRs) play vital roles in adaptive immune system of agnathan vertebrate. In the present study, we first discover a novel VLR gene, VLR2, from an invertebrate, the Chinese mitten crab, Eriocheir sinensis. VLR2 has ten different isoforms formed via alternative splicing, which is different from that in agnathan vertebrate with the assembly of LRR modules. The longest isoform, VLR2-L, responds to Gram-positive bacteria Staphylococcus aureus challenge specifically, while shows no response to Gram-negative bacteria Vibrio parahaemolyticus challenge, confirmed by recombinant expression and bacterial binding experiments. Interestingly, VLR2s with short LRRs regions (VLR2-S8 and VLR2-S9) tend to bind to Gram-negative bacteria rather than Gram-positive bacteria. Antibacterial activity assay proves six isoforms of VLR2 have pluralistic antibacterial effects on bacteria which were never reported in invertebrate. These results suggest that the diversity and specificity of VLR2 resulted from alternative splicing and the length of the LRRs region. This pathogen-binding receptor diversity will lay the foundation for the study of immune priming. Furthermore, studying the immune function of VLR2 will provide a new insight into the disease control strategy of crustacean culture.

New insight of variable lymphocyte receptor-likes in anti-bacteria activity from Eriocheir sinensis

The Chinese mitten crab, Eriocheir sinensis, is a vital freshwater aquaculture species in China, however, is also facing various crab disease threats. In the present study, we identify three novel variable lymphocyte receptor-like (VLR-like) genes-VLR-like1, VLR-like3 and VLR-like4-from E. sinensis, which play vital roles in adaptive immune system of agnathan vertebrates. The bacterial challenge, bacterial binding and antibacterial-activity experiments were applied to study immune functions of VLR-likes, and the transcriptomic data from previous E. sinensis bacterial challenge experiments were analyzed to speculate the possible signaling pathway. VLR-like1 and VLR-like4 can respond to Staphylococcus aureus challenge and inhibit S. aureus specifically. VLR-like1 and VLR-like4 possess broad-spectrum bacteria-binding ability whereas VLR-like3 do not. VLR-likes in E. sinensis could associate with the Toll-like receptor (TLR) signaling pathway. The above results suggest that VLR-likes defend against bacteria invasion though exerting anti-bacteria activity, and probably connect with the TLR signaling pathway. Furthermore, studying the immune functions of these VLR-likes will provide a new insight into the disease control strategy of crustacean culture.

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.

Hallmarks of crustacean immune hemocytes at single-cell resolution

In invertebrates, hemocytes are the key factors in innate immunity. However, the types of invertebrate immune hemocytes are unclassified due to the limitation of morphological classification. To determine the immune hemocytes of crustaceans, the heterogeneity of hemocytes of shrimp Marsupenaeus japonicus and crayfish Procambarus clarkii, two representative crustacean species, were characterized in this study. The results of single-cell RNA sequencing indicated that shrimp and crayfish contained 11 and 12 types of hemocytes, respectively. Each of different types of hemocytes specifically expressed the potential marker genes. Based on the responses of shrimp and crayfish to the infection of white spot syndrome virus (WSSV) and the challenge of lipopolysaccharide (LPS), four types of immune hemocytes of crustaceans were classified, including semi-granular hemocytes involved in antimicrobial peptide production, granular hemocytes responsible for the production of antimicrobial peptides, hemocytes related to cell proliferation and hemocytes in immunity-activated state. Therefore, our study provided the first classification of crustacean hemocytes as well as of immune hemocytes of crustaceans at the single-cell resolution, which would be helpful to understand the innate immunity of invertebrates.

Transcriptome analysis reveals the molecular mechanism of long-term exposure of Eriocheir sinensis to low concentration of trichlorfon

Trichlorphon, a common organophosphorus pesticide (OPs), is widely used in aquaculture to prevent aquatic insects from infecting cultured objects as well as to control the excessive proliferation of plankton in water bodies. However, its repeated use time can contaminate water bodies and impart direct/indirect toxicity to beneficial aquatic species. However, the underlying mechanism regarding toxicity and cellular metabolism remains unclear. Understanding the mechanism would enable the standardized use and management of OPs and their use in the aquatic environment. Here, low concentration of trichlorphon (5 × 10^-5 g/L) was used to construct a hepatopancreatic transcriptional library 30 d, 60 d and 90 d after exposure using RNA-Seq. We detected 649, 148, and 2949 DEGs in the hepatopancreas of E. sinensis for the Tri01 vs. Ctr01, Tri02 vs. Ctr02 and Tri03 vs. Ctr03 library, respectively. The results of KEGG pathway enrichment analysis showed that DEGs were mainly enriched in signal transduction, carbohydrate metabolism, transport and catabolism, endocrine system, and digestive system. Also, under trichlorfon stress, DEGs of E. sinensis were enriched in thyroid hormone signaling pathways, protein digestion and absorption, cancer pathways, etc. The significant DEGs were mainly related to metabolism and the apoptosis and autophagy pathways. This study lays a foundation for further revealing the effects of long-term trichlorfon stress on E. sinensis as well as the potential physiological toxicity. The relevant transcriptome data could provide a reference for the molecular toxicological evaluation of trichlorfon in aquaculture.

Comparison on nutritional quality of adult female Chinese mitten crab (Eriocheir sinensis) with different colored hepatopancreases

White hepatopancreas syndrome has recently emerged in Chinese mitten crab (Eriocheir sinensis) aquaculture, causing considerable economic loss. The hepatopancreas color of diseased crabs becomes gradually lighter, turning from yellow to yellow-white to white. Therefore, this study was conducted to investigate the changes in nutrient composition in three edible parts (hepatopancreas, ovaries, and muscle) of adult females with different colored hepatopancreases. Three groups were assessed in this study, including a yellow hepatopancreas group (control, L * = 63.92, a * = 22.14, b * = 60.95), a yellow-white hepatopancreas group (YWHG, L * = 65.06, a * = 22.35, b * = 57.80), and a white hepatopancreas group (WHG, L * = 65.72, a * = 10.70, b * = 30.52). No statistically significant differences in average weight, tissue indices, and total edible yield were observed among the three crab groups (P >0.05). The moisture content of the hepatopancreases and ovaries in WHG was 56.12% and 9.23% higher than the control values (P <0.05), whereas hepatopancreas crude fat and ovary crude protein levels were 62.23% and 11.45% lower than the control values (P < 0.05). The total carbohydrate levels of the three edible tissues were significantly higher and the crude protein content of ovaries was significantly lower in YWHG (P < 0.05). Most amino acid levels in the WHG muscle and ovaries were significantly lower than the control (P < 0.05). Moreover, the hepatopancreas levels of total polyunsaturated fatty acids (PUFA) and n-6PUFA in WHG were 24.88% and 31.83% lower than in control group, whereas the hepatopancreas levels of total PUFA and n-6PUFA in YWHG were also 21.88% and 23.20% lower compared to the controls (P < 0.05). Overall, the growth and the edible parts were not affected in YWHG and WHG. Moreover, YWHG crabs exhibited few effects on nutritional value; however, the fatty acid composition of crabs was significantly changed. In contrast, WHG crabs exhibited poor nutritional quality. Nonetheless, the consumption of crabs with yellow-white or white hepatopancreases is not recommended since the animal also referred to as diseased crabs.

Reconstruction of Eriocheir sinensis Y-organ Genome-Scale Metabolic Network and Differential Analysis After Eyestalk Ablation

Genome-scale metabolic network (GSMN) has been proven to be a useful tool for the system analysis of organism metabolism and applied to deeply explore the metabolic functions or mechanisms in many organisms, including model or non-model organisms. However, the systematic studies on the metabolisms of aquatic animals are seldom reported, especially the aquatic crustaceans. In this work, we reconstructed an Eriocheir sinensis Y-organ GSMN based on the transcriptome sequencing of Y-organ, which includes 1,645 reactions, 1,885 unigenes, and 1,524 metabolites distributed in 100 pathways and 11 subsystems. Functional module and centrality analysis of the GSMN show the main metabolic functions of Y-organ. Further analysis of the differentially expressed unigenes in Y-organ after eyestalk ablation reveals that 191 genes in the network were up-regulated and 283 were down-regulated. The unigenes associated with the ecdysone synthetic pathway were all up-regulated, which is consistent with the report on the increasing secretion of ecdysone after eyestalk ablation. Besides, we compared the Y-organ GSMN with that of E. sinensis eyestalk and hepatopancreas, and we analyzed the specific metabolisms in each organ. The specific metabolisms and pathways of these three networks are closely related to their corresponding metabolic functions. The GSMN reconstructed in this work provides a new method and many novel clues for further understanding the physiological function of Y-organ. It also supplies a new platform for the investigation of the interactions among different organs in the growth process of E. sinensis.

Transcriptome analysis of hepatopancreas of Eriocheir sinensis with hepatopancreatic necrosis disease (HPND)

Hepatopancreatic necrosis disease (HPND) is a newly emerging disease in the Chinese mitten crab, Eriocheir sinensis, which has resulted in large economic losses. However, the underlying cause of this disease remains unclear. To better understand the pathogenesis and pathogenic mechanism of HPND, we compared the transcriptome differences of the hepatopancreas of E. sinensis with and without HPND. The analysis yielded > 30 million reads for each sample of three test (with HPND) and three control groups (without HPND). We observed 978 downregulated genes and 644 upregulated genes. Among the gene ontology categories "biological process," "cellular component," and "molecular function", the subcategories cellular process, single-organism process, biological regulation, metabolic process, cell part, organelle, organelle part, binding, and catalytic were enriched. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis showed that "metabolism of xenobiotics by cytochrome P450," "drug metabolism-cytochrome P450," "chemical carcinogenesis," and "material metabolism" were the "five" most significantly enriched pathways in the hepatopancreas of E. sinensis with HPND. The results revealed that material metabolic abnormalities and drug effects from the external environment might be associated with HPND in the Chinese mitten crab. Considering the wide use of pyrethroids for pond cleaning in Xinghua city, we speculated that pyrethroids might cause HPND in the Chinese mitten crab. Our study provided useful information about the cause and pathogenetic mechanisms of HPND and could help to prevent this disease in production practice.

The Single-molecule long-read sequencing of Scylla paramamosain

Scylla paramamosain is an important aquaculture crab, which has great economical and nutritional value. To the best of our knowledge, few full-length crab transcriptomes are available. In this study, a library composed of 12 different tissues including gill, hepatopancreas, muscle, cerebral ganglion, eyestalk, thoracic ganglia, intestine, heart, testis, ovary, sperm reservoir, and hemocyte was constructed and sequenced using Pacific Biosciences single-molecule real-time (SMRT) long-read sequencing technology. A total of 284803 full-length non-chimeric reads were obtained, from which 79005 high-quality unique transcripts were obtained after error correction and sequence clustering and redundant. Additionally, a total of 52544 transcripts were annotated against protein database (NCBI nonredundant, Swiss-Prot, KOG, and KEGG database). A total of 23644 long non-coding RNAs (lncRNAs) and 131561 simple sequence repeats (SSRs) were identified. Meanwhile, the isoforms of many genes were also identified in this study. Our study provides a rich set of full-length cDNA sequences for S. paramamosain, which will greatly facilitate S. paramamosain research.

Distribution, fate and risk assessment of PAHs in water and sediments from an aquaculture- and shipping-impacted subtropical lake, China

The spatial-temporal distribution of polycyclic aromatic hydrocarbons (PAHs), their source, and potential health risks were determined in overlying water and surface sediments from Chinese Lake Guchenghu, adjacent commercial mitten crab ponds and the connected Wushen Canal to assess the contamination profile of the area. The total PAHs concentrations in sediment and water were 86.7-1790 ng g^-1 dry weight (dw) and 184-365 ng L^-1 in summer and 184-3140 ng g^-1 dw and 410-1160 ng L^-1 in winter. Two- and 3-ring PAHs were the predominant compounds in water, while PAHs with 4-6 rings dominated in the sediment at both upstream and downstream sites. PAHs concentrations in water and sediment correlated significantly. Diagnostic ratios and positive matrix factorization (PMF) analyses indicated a strong influence of pyrogenic sources, principally biomass combustion and vehicle emission, on the concentrations of PAHs. The distribution, source identification, and mean effects range median quotients (mERMQ) analyses suggested that the most contaminated area was located downstream and upstream of the Wushen Canal, followed by Lake Guchenghu and a commercial crab pond area. From an ecological point of view, PAHs posed a potential risk to drinking water sources as the concentrations exceeded the guideline value of 0.05 μg L^-1. The risk posed by sediment PAHs appeared to be low except for the downstream sites, which showed a low to medium ecotoxicological risk. The total incremental lifetime cancer risks ranged between 10^-7 and 10^-5, indicating a potential health risk for the local population when exposed to sediment from the area.

Transcriptomic analysis of Portunus trituberculatus reveals a critical role for WNT4 and WNT signalling in limb regeneration

The swimming crab (Portunus trituberculatus) is among the most economically important seawater crustacean species in Asia. Despite its commercial importance and being well-studied status, genomic and transcriptomic data are scarce for this crab species. In the present study, limb bud tissue was collected at different developmental stages post amputation for transcriptomic analysis. Illumina RNA-sequencing was applied to characterise the limb regeneration transcriptome and identify the most characteristic genes. A total of 289,018 transcripts were obtained by clustering and assembly of clean reads, producing 150,869 unigenes with an average length of 956 bp. Subsequent analysis revealed WNT signalling as the key pathway involved in limb regeneration, with WNT4 a key mediator. Overall, limb regeneration appears to be regulated by multiple signalling pathways, with numerous cell differentiation, muscle growth, moult, metabolism, and immune-related genes upregulated, including WNT4, LAMA, FIP2, FSTL5, TNC, HUS1, SWI5, NCGL, SLC22, PLA2, Tdc2, SMOX, GDH, and SMPD4. This is the first experimental study done on regenerating claws of P. trituberculatus. These findings expand existing sequence resources for crab species, and will likely accelerate research into regeneration and development in crustaceans, particularly functional studies on genes involved in limb regeneration.