Aquatic organisms have to produce proteins or factors that help maintain a stable relationship with microbiota and prevent colonization by pathogenic microorganisms. In crustaceans and other aquatic invertebrates, relatively few of these host factors have been characterized. In this study, we show that the respiratory glycoprotein hemocyanin is a crucial host factor that modulates microbial composition and diversity in the hepatopancreas of penaeid shrimp. Diseased penaeid shrimp (Penaeus vannamei), had an empty gastrointestinal tract with atrophied hepatopancreas, expressed low hemocyanin, and high total bacterial abundance, with Vibrio as the dominant bacteria. Similarly, shrimp depleted of hemocyanin had mitochondrial depolarization, increased reactive oxygen species (ROS) levels, and dysregulation of several energy metabolism-related genes. Hemocyanin silencing together with ROS scavenger (N-acetylcysteine) treatment improved microbial diversity and decreased Vibrio dominance in the hepatopancreas. However, fecal microbiota transplantation after hemocyanin knockdown could not restore the microbial composition in the hepatopancreas. Collectively, our data provide, to our knowledge, new insight into the pivotal role of hemocyanin in modulating microbial composition in penaeid shrimp hepatopancreas via its effect on mitochondrial integrity, energy metabolism, and ROS production.
Crustacean molting is highly related to energy and lipid metabolism. This study was conducted to detect the changes of total lipids (TL), triacylglyceride (TAG), phospholipid (PL) and lipid droplets in hepatopancreas, and then to investigate the gene expression patterns related to hepatopancreatic lipid metabolism during the molting cycle of Chinese mitten crab Eriocheir sinensis. Hepatopancreatic TL and TAG increased significantly from post-molt stage to pre-molt stage, then decreased significantly from pre-molt stage to ecdysis stage, which is consistent to the changes of neutral lipid-rich adipocytes in hepatopancreas. By transcriptomic analysis, 65,325 transcripts were sequenced and assembled, and 28,033 transcripts were annotated. Most genes were related to energy metabolism, and the enriched genes were involved in carbohydrate and lipid metabolism and biosynthesis, especially in de novo synthesis of fatty acids and TAG, and ketone body production. Compared to the inter-molt stages, acetyl-CoA carboxylase, fatty acid synthase and other genes related to the synthesis of fatty acids were upregulated in the pre-molt stage. TAG synthesis related genes, including Glycerol-3-phosphate acyltransferase and 1-acylglycerol-3-phosphate acyltransferases, were upregulated in the post-molt stage compared to the inter-molt stage. The expression of ketone body-related genes had no significant changes during the molting cycle. Compared to the TAG synthetic pathway, ketone body biosynthesis may contribute less/secondarily to fatty acid metabolic processes, which could be involved in the other physiological processes or metabolism. In conclusion, these results showed that TAG is the major lipid deposition during inter- and pre-molt stages, and the most genes are related to the fatty acids and TAG metabolism in the hepatopancreas during the molting cycle of E. sinensis.
Plastics are widely produced for industrial and domestic applications due to their unique properties, and studies on the toxic effects of nanoplastics (NPs) on aquatic animals are essential. In this study, we investigated the transcriptomic patterns of Litopenaeus vannamei after NPs exposure. We found that the lysosome pathway was activated when after NPs exposure, with up-regulated DEGs, including glucocerebrosidase (GBA), hexosaminidase A (HEXA), sphingomyelin phosphodiesterase-1 (SMPD1), and solute carrier family 17 member 5 (SLC17A5). In addition, the PI3K-Akt signaling pathway was strongly affected by NPs, and the upstream genes of PI3K-Akt, including epidermal growth factor receptor (EGFR), integrin subunit beta 1 (ITGB1) and heat shock protein 90 (HSP90) were up-regulation. Other genes involved in lipogenesis, such as sterol regulatory element binding transcription factor 1 (SREBP-1c), fatty acid synthase (FASN) and stearoyl-CoA desaturase (SCD-1), were down-regulated. However, the contents of triglycerides (TG) and total cholesterol (TCH) in L. vanname hepatopancreas were reduced, which indicated that the ingestion of NPs led to the disturbance of hepatic lipid metabolism. What more, NPs treatment of L. vannamei also caused oxidative stress. In addition, NPs can damage part of the tissue structure and affect the physiological function of shrimps. The results of this study provide valuable ecotoxicological data to improve the understanding of the biological fate and effects of nanoplastics in L. vannamei.
MicroRNAs (miRNAs) are ~20-22 nucleotides, non protein-coding RNA regulatory genes that post-transcriptionally regulate many protein-coding genes, influencing critical biological and metabolic processes. While the number of known microRNA is increasing, there is currently no published data for miRNA from giant freshwater prawns, Macrobrachium rosenbergii (M. rosenbergii), a commercially cultured and economically important food species. In this study, we identified novel miRNAs in the gill and hepatopancreas of M. rosenbergii. Through a deep parallel sequencing analysis and an in silico data analysis approach, 327 miRNA families were identified from small RNA libraries with reference to both the de novo transcriptome of M. rosenbergii obtained from RNA-Seq and to miRBase (Release 18.0, November 2012). Based on the identified mature miRNA and recovered precursor sequences that form appropriate hairpin structures, three conserved miRNA (miR125, miR750, miR993) and 27 novel miRNA candidates encoding messenger-like non-coding RNA were identified. miR-125, miR-750, G-m0002/H-m0009, G-m0005, G-m0008/H-m0016, G-m0011/H-m0027 and G-m0015 were selected for experimental validation with stem-loop quantitative RT-PCR and were found to be coherent with the expression profile of deep sequencing data as evaluated with Pearson's correlation coefficient (r = 0.835178 for miRNA in gill, r = 0.724131 for miRNA in hepatopancreas). Using a combinatorial approach of pathway enrichment analysis and inverse expression relationship of miRNA and mRNA, four co-expressed novel miRNA candidates (G-m0005, G-m0008/H-m0016, G-m0011/H-m0027, and G-m0015) were found to be associated with energy metabolism. In addition, the expression of the three novel miRNA candidates (G-m0005, G-m0008/H-m0016, and G-m0011/H-m0027) were also found to be significantly reduced at 9 and 24 h post infection in M. rosenbergii challenged with infectious hypodermal and hematopoietic necrosis virus, suggesting a functional role of these miRNAs in crustacean immune defense.
Gene discovery in the Malaysian giant freshwater prawn (Macrobrachium rosenbergii) has been limited to small scale data collection, despite great interest in various research fields related to the commercial significance of this species. Next generation sequencing technologies that have been developed recently and enabled whole transcriptome sequencing (RNA-seq), have allowed generation of large scale functional genomics data sets in a shorter time than was previously possible. Using this technology, transcriptome sequencing of three tissue types: hepatopancreas, gill and muscle, has been undertaken to generate functional genomics data for M. rosenbergii at a massive scale. De novo assembly of 75-bp paired end Ilumina reads has generated 102,230 unigenes. Sequence homology search and in silico prediction have identified known and novel protein coding candidate genes (∼24%), non-coding RNA, and repetitive elements in the transcriptome. Potential markers consisting of simple sequence repeats associated with known protein coding genes have been successfully identified. Using KEGG pathway enrichment, differentially expressed genes in different tissues were systematically represented. The functions of gill and hepatopancreas in the context of neuroactive regulation, metabolism, reproduction, environmental stress and disease responses are described and support relevant experimental studies conducted previously in M. rosenbergii and other crustaceans. This large scale gene discovery represents the most extensive transcriptome data for freshwater prawn. Comparison with model organisms has paved the path to address the possible conserved biological entities shared between vertebrates and crustaceans. The functional genomics resources generated from this study provide the basis for constructing hypotheses for future molecular research in the freshwater shrimp.
Considering the importance of heat shock proteins (HSPs) in the innate immune system of prawn, a comparative molecular approach was proposed to study the crustacean large HSPs 60, 70 and 90. Three different large HSPs were identified from freshwater prawn Macrobrachium rosenbergii (Mr) cDNA library during screening. The structural and functional characteristic features of HSPs were studied using various bioinformatics tools. Also, their gene expression and mRNA regulation upon various pathogenic infections was studied by relative quantification using 2(-ΔΔCT) method. MrHSP60 contains a long chaperonin 60 domain at 46-547 which carries a chaperonin 60 signature motif between 427 and 438, whereas MrHSP70 contains a long HSP70 domain at 21-624 and MrHSP90 carries a HSP90 domain at 188-719. The two dimensional analysis showed that MrHSP60 contains more amino acids (52%) in helices, whereas MrHSP70 (40.6%) and MrHSP90 (51.8%) carried more residues in coils. Gene expression results showed significant (P
Apoptosis is genetically programmed cellular killing processes that execute unnecessary or infected cells. It plays an important role in embryogenesis, homeostasis, insect metamorphosis and immunity. Apoptosis inhibitor (MrIAP) was sequenced from the freshwater giant prawn Macrobrachium rosenbergii using Illumina Solexa Genome Analyzer Technique. MrIAP consisted of 1753 base pair nucleotides encoded 535 polypeptide with an estimated molecular mass of 60 kDa. MrIAP amino acid sequence contains IAP superfamily domain between 5 and 490. The deduced amino acid sequences of the MrIAP were aligned with the other IAP family members. The highest sequence similarity was observed in IAP-5 from ant Camponotus floridanus (67%) followed by IAP from body louse Pediculus humanus corporis (66%) and the lowest (62%) in IAP-5 isoform-5 from common chimpanzee Pan troglodytes and IAP-5 from Aedes aegypti. The IAP phylogenetic tree showed that MrIAP closely related to other arthropod blacklegged tick Ixodes scapularis, formed a sister group with IAP from a hemichordate acorn worm Saccoglossus kowalevskii and finally clustered together with IAPs from fish groups. The quantitative real time PCR analysis revealed that significantly (P < 0.05) highest expression was noticed in hepatopancreas and significantly (P < 0.05) lowest expression in pleopods. Based on the results of gene expression analysis, MrIAP mRNA transcription in M. rosenbergii challenged to infectious hypodermal and hematopoietic necrosis virus (IHHNV) was highly induced in hepatopancreas. The collective results of this study indicate that the MrIAP is an essential immune gene and influences the immune response against IHHNV infection in M. rosenbergii.
Understanding Macrobrachium rosenbergii ovarian maturation control at the genome level is an important aspect for increasing larvae production. In this study, an ovarian maturation related gene, M. rosenbergii vWD domain and three Kazal-type domains of a gene (MrvWD-Kazal) have been studied. The MrvWD-Kazal gene was isolated using a rapid amplification of cDNA end (RACE) method and the relative abundances of MrvWD-Kazal mRNA in the ovary, hepatopancreas, stomach, intestine and gill were determined by using the quantitative PCR technique. The MrvWD-Kazal gene is composed of 2194 bp with an open reading frame (ORF) of 1998 bp encoding 665 amino acids and has great similarity to the M. nipponense vWD-Kazal gene (91%). The qPCR analyses indicated the relative abundance of MrvWD-Kazal mRNA transcript varied among different stages of ovarian function (P < 0.05), but there were no differences abundance in hepatopancreas, stomach, intestine and gill (P> 0.05). In the ovary, relative abundance of MrvWD-Kazal mRNA transcript gradually increased with ovarian maturation from Stages 1 (Spent; 1.00-fold), to 2 (Proliferative; 3.47-fold) to 3 (Premature; 6.18-fold) and decreased at Stage 4 (Mature; 1.31-fold). Differential relative abundances of MrvWD-Kazal mRNA transcript in the ovary indicate the MrvWD-Kazal protein may have an important function in ovarian maturation of M. rosenbergii. The results of this study also indicate the MrvWD-Kazal is not involved in regulation of the reproductive related function of the hepatopancreas, digestive system (stomach and intestine) and respiratory system (gill).