Isobaric Tags for Relative and Absolute Quantitation-Based Proteomics Analysis Provides a New Perspective Into Unsynchronized Growth in Kuruma Shrimp (Marsupenaeus japonicus)

Unsynchronized growth is a common phenomenon in farmed crustaceans. The underlying molecular mechanism of unsynchronized growth of crustaceans is unclear. In this study, a comparative proteomic analysis focusing on growth differences was performed using kuruma shrimp Marsupenaeus japonicus, an economic crustacean species, as the model. The study analyzed kuruma shrimp at fast growth stage and steady growth stage from both fast growth group and slow growth group by an Isobaric tags for relative and absolute quantitation (iTRAQ)-based quantitative proteomic analysis method. A total of 1,720 proteins, including 12,291 peptides, were identified. Fifty-two and 70 differentially expressed proteins (DEPs) were identified in the fast growth stage and steady growth stage, respectively. Interestingly, 10 DEPs, including 14-3-3-epsilon-like, GPI, GPD1, MHC-1a, and MHC-1b, were presented in both growth stages. In addition, all these 10 DEPs shared the same expression tendency at these two growth stages. The results indicated that these 10 DEPs are potential growth biomarkers of M. japonicus. Proteins associated with faster growth of M. japonicus may promote cell growth and inhibit cell apoptosis through the Hippo signaling pathway. The fast growth group of M. japonicus may also achieve growth superiority by activating multiple related metabolic pathways, including glycolysis, glycerophospholipid metabolism and Citrate cycle. The present study provides a new perspective to explore the molecular mechanism of unsynchronized growth in crustacean species.

Comparative analysis of transcriptomic data shows the effects of multiple evolutionary selection processes on codon usage in Marsupenaeus japonicus and Marsupenaeus pulchricaudatus

Background Kuruma shrimp, a major commercial shrimp species in the world, has two cryptic or sibling species, Marsupenaeus japonicus and Marsupenaeus pulchricaudatus. Codon usage analysis would contribute to our understanding of the genetic and evolutionary characteristics of the two Marsupenaeus species. In this study, we analyzed codon usage and related indices using coding sequences (CDSs) from RNA-seq data. Results Using CodonW 1.4.2 software, we performed the codon bias analysis of transcriptomes obtained from hepatopancreas tissues, which indicated weak codon bias. Almost all parameters had similar correlations for both species. The gene expression level (FPKM) was negatively correlated with A/T3s. We determined 12 and 14 optimal codons for M. japonicus and M. pulchricaudatus, respectively, and all optimal codons have a C/G-ending. The two Marsupenaeus species had different usage frequencies of codon pairs, which contributed to further analysis of transcriptional differences between them. Orthologous genes that underwent positive selection (ω > 1) had a higher correlation coefficient than that of experienced purifying selection (ω < 1). Parity Rule 2 (PR2) and effective number of codons (ENc) plot analysis showed that the codon usage patterns of both species were influenced by both mutations and selection. Moreover, the average observed ENc value was lower than the expected value for both species, suggesting that factors other than GC may play roles in these phenomena. The results of multispecies clustering based on codon preference were consistent with traditional classification. Conclusions This study provides a relatively comprehensive understanding of the correlations among codon usage bias, gene expression, and selection pressures of CDSs for M. japonicus and M. pulchricaudatus. The genetic evolution was driven by mutations and selection pressure. Moreover, the results point out new insights into the specificities and evolutionary characteristics of the two Marsupenaeus species.

Effects of artificial substrates on the growth and immunology of postlarvae of Marsupenaeus japonicus (Spence Bate, 1888) (Decapoda: Dendrobranchiata: Penaeidae) reared in biofloc

We investigated the effects of biofloc and artificial substrates (net and brush) on the growth and immunology of the postlarvae of Marsupenaeus japonicus (Spence Bate, 1888) reared for 30 days. The mRNA expressions of innate immune-related genes (prophenoloxidase, masquerade-like serine proteinase, and lysozyme) in the postlarvae were analyzed using quantitative reverse transcription PCR (qRT-PCR). The highest specific growth rate of larvae (3.66 ± 0.02% d–1, 3.75 ± 0.02% d–1) were shown in the net and control groups. The highest survival rate (90 ± 7.1%) was observed in the brush group (P &lt; 0.05). The net yield was significantly higher in the artificial substrate groups (brush 36.60 ± 6.19 g m–2, net 36.46 ± 2.36 g m–2) than in the control (33.79 ± 0.16 g m–2) (P &lt; 0.05). Total suspended, volatile suspended, and suspended solids were significantly lower in the brush group than in the other groups. The immune-related genes showed significantly higher expressions in the artificial substrate groups than in the control. These results support the benefits of artificial substrate in M. japonicus nursery culture to mitigate space competition and suppress cannibalism, with a positive effect on survival. The additional surface area provided by the biofloc, which the larvae can graze on, enhanced the expression of immune-related genes in individuals.

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 LC50 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.

Genome and transcriptome assemblies of the kuruma shrimp, Marsupenaeus japonicus

The kuruma shrimp Marsupenaeus japonicus (order Decapoda, family Penaeidae) is an economically important crustacean that occurs in shallow, warm seas across the Indo-Pacific. Here, using a combination of Illumina and Oxford Nanopore Technologies platforms, we produced a draft genome assembly of M. japonicus (1.70 Gbp; 18,210 scaffolds; scaffold N50=234.9 kbp; 34.38% GC, 93.4% BUSCO completeness) and a complete mitochondrial genome sequence (15,969 bp). As with other penaeid shrimp genomes, the M. japonicus genome is extremely rich in simple repeats, which occupies 27.4% of the assembly. A total of 26,381 protein-coding gene models (94.7% BUSCO completeness) were predicted, of which 18,005 genes (68.2%) were assigned functional description by at least one method. We also produced an Illumina-based transcriptome shotgun assembly (40,991 entries; 93.0% BUSCO completeness) and a PacBio Iso-Seq transcriptome assembly (25,415 entries; 67.5% BUSCO completeness). We envision that the M. japonicus genome and transcriptome assemblies will serve as useful resources for the basic research, fisheries management, and breeding programs of M. japonicus.