Characterization, expression, and function analysis of gonad-inhibiting hormone in Oriental River prawn, Macrobrachium nipponense and its induced expression by temperature

Whole transcriptome RNA sequencing provides novel insights into the molecular dynamics of ovarian development in mud crab, Scylla paramamosain after mating

Ovarian development in animals is a complicated biological process, requiring the simultaneous coordination among various genes and pathways. To understand the dynamic changes and molecular regulatory mechanisms of ovarian development in mud crab (Scylla paramamosain), both histological observation and whole transcriptome sequencing of ovarian tissues at different mating stages were implemented in this study. The histological results revealed that ovarian development was delayed in unmated females (60 days after courtship behavior but not mating), who exhibited an oocyte diameter of 56.38 ± 15.17 μm. Conversely, mated females exhibited accelerated the ovarian maturation process, with females reaching ovarian stage III (proliferative stage) 23 days after mating and attained an average oocyte diameter of 132.19 ± 15.07 μm. Thus, mating process is essential in promoting the rapid ovarian development in mud crab. Based on the whole transcriptome sequencing analysis, a total of 518 mRNAs, 1502 lncRNAs, 18 circRNAs and 151 miRNAs were identified to be differentially expressed between ovarian tissues at different mating stages. Notably, six differentially expressed genes (DEGs) associated with ovarian development were identified, including ovary development-related protein, red pigment concentrating hormone receptor, G2/mitotic-specific cyclin-B3-like, lutropin-chorio gonadotropic hormone receptor, renin receptor, and SoxB2. More importantly, both DEGs and targets of differentially expressed non-coding RNAs (DEncRNAs) were enriched in renin-angiotensin system, TGF-β signaling, cell adhesion molecules, MAPK signaling pathway, and ECM-receptor interaction, suggesting that these pathways may play significant roles in the ovarian development of mud crabs. Moreover, competition endogenous RNA (ceRNA) networks were constructed while mRNAs were differentially expressed between mating stages were involved in Gene Ontology (GO) biological processes such as developmental process, reproduction, and growth. These findings could provide solid foundations for the future development of female mud crab maturation enhancement strategy, and improve the understanding of the ovarian maturation process in crustaceans.

A new type II CHH neuropeptide involves ovarian development in the peppermint shrimp, Lysmata vittata

Type II crustacean hyperglycemic hormone (CHH) neuropeptides play diverse roles in crustaceans. In the hermaphrodite shrimp Lysmata vittata, two transcripts of type II CHHs (molt-inhibiting hormone/gonad-inhibiting hormone, MIH/GIH1 and MIH/GIH2) were identified by transcriptome sequencing, and MIH/GIH1 was later named Lvit-GIH1 for its inhibitory effect on ovarian development. Based on the high similarity of MIH/GIH2 to Lvit-GIH1, we named tentatively MIH/GIH2 as Lvit-GIH2 and explored the role of Lvit-GIH2 in ovarian development. The open reading frame (ORF) of Lvit-GIH2 was 333 bp in length, encoding a precursor consisted of a 32-aa signal peptide and a 78-aa mature peptide, which shared high sequence similarity with the type II subfamily peptides in crustaceans. Notably, Lvit-GIH2 was widely expressed in multiple tissues. The qRT-PCR findings indicated a rising trend in the expression of Lvit-GIH2 from the male phase to the euhermaphrodite phase. Both RNA interference and addition of GIH2 recombinant proteins (rGIH2) experiments showed that Lvit-GIH2 suppressed Lvit-Vg expression in hepatopancreas and Lvit-VgR expression in ovary. To further investigate the role of Lvit-GIH2 in ovarian development, the RNA-sequence analysis was performed to examine the changes in ovary after addition of rGIH2. The results showed that the pathways (Cysteine and methionine metabolism, Apoptosis-multiple species, etc.) and the genes (17bHSD8, IGFR, CHH, etc.) related to ovarian development were negatively regulated by rGIH2. In brief, Lvit-GIH2 might inhibit the ovarian development in L. vittata.

NPC Intracellular Cholesterol Transporter 1 Regulates Ovarian Maturation and Molting in Female Macrobrachium nipponense

NPC intracellular cholesterol transporter 1 (NPC1) plays an important role in sterol metabolism and transport processes and has been studied in many vertebrates and some insects, but rarely in crustaceans. In this study, we characterized NPC1 from Macrobrachium nipponense (Mn-NPC1) and evaluated its functions. Its total cDNA length was 4283 bp, encoding for 1344 amino acids. It contained three conserved domains typical of the NPC family (NPC1_N, SSD, and PTC). In contrast to its role in insects, Mn-NPC1 was mainly expressed in the adult female hepatopancreas, with moderate expression in the ovary and heart. No expression was found in the embryo (stages CS-ZS) and only weak expression in the larval stages from hatching to the post-larval stage (L1-PL15). Mn-NPC1 expression was positively correlated with ovarian maturation. In situ hybridization showed that it was mainly located in the cytoplasmic membrane and nucleus of oocytes. A 25-day RNA interference experiment was employed to illustrate the Mn-NPC1 function in ovary maturation. Experimental knockdown of Mn-NPC1 using dsRNA resulted in a marked reduction in the gonadosomatic index and ecdysone content of M. nipponense females. The experimental group showed a significant delay in ovarian maturation and a reduction in the frequency of molting. These results expand our understanding of NPC1 in crustaceans and of the regulatory mechanism of ovarian maturation in M. nipponense.

Insulin-like Androgenic Gland Hormone Induced Sex Reversal and Molecular Pathways in Macrobrachium nipponense: Insights into Reproduction, Growth, and Sex Differentiation

This study investigated the potential to use double-stranded RNA insulin-like androgenic gland hormone (dsIAG) to induce sex reversal in Macrobrachium nipponense and identified the molecular mechanisms underlying crustacean reproduction and sex differentiation. The study aimed to determine whether dsIAG could induce sex reversal in PL30-male M. nipponense during a critical period. The sex-related genes were selected by performing the gonadal transcriptome analysis of normal male (dsM), normal female (dsFM), neo-female sex-reversed individuals (dsRM), and unreversed males (dsNRM). After six injections, the experiment finally resulted in a 20% production of dsRM. Histologically, dsRM ovaries developed slower than dsFM, but dsNRM spermathecae developed normally. A total of 1718, 1069, and 255 differentially expressed genes were identified through transcriptome sequencing of the gonads in three comparison groups, revealing crucial genes related to reproduction and sex differentiation, such as GnRHR, VGR, SG, and LWS. Principal Component Analysis (PCA) also distinguished dsM and dsRM very well. In addition, this study predicted that the eyestalks and the "phototransduction-fly" photoperiodic pathways of M. nipponense could play an important role in sex reversal. The enrichment of related pathways and growth traits in dsNRM were combined to establish that IAG played a significant role in reproduction, growth regulation, and metabolism. Finally, complete sex reversal may depend on specific stimuli at critical periods. Overall, this study provides valuable findings for the IAG regulation of sex differentiation, reproduction, and growth of M. nipponense in establishing a monoculture.

Dual roles of CYP302A1 in regulating ovarian maturation and molting in Macrobrachium nipponense

CYP302A1 is a member of the Halloween genes in the cytochrome P450 supergene family, which play an important regulatory role in the synthesis of 20-hydroxyecdysone (20E) in crustaceans and insects. In this study, we found that the Mn-CYP302A1sequence included typical CYP450 conserved domains. Phylogenic showed that it is closely related to crustaceans and insects. q-PCR analysis indicated that Mn-CYP302A1 was highly expressed in the ovaries and peaked before ovarian maturation. Mn-CYP302A1 expression was higher at the post-larval stage of day 15 than at other stages of embryogenesis. In situ hybridization indicated that Mn-CYP302A1 was mainly distributed in the nucleus, yolk granules, cell membrane and cytoplasm To further establish the function of CYP302A1, a 21-day RNA interference experiment was conducted. On day 16, the Gonad Somatic Index of the control group and the experimental group showed significant differences, with GSI of 11.72% in the control group and 3.21% in the experimental group. The cumulative proportion of the second entry into stage O-Ⅲ was 100% in the control group, while it was 41.67% in the experimental group on day 21. The ecdysone content was 8.91nmol/L in the control group and 6.11nmol/L in the experimental group on day 9. A significant difference in the molting proportion between the control group and the experimental group was also observed (49% in the control group and 34% in the experimental group) on day 16. Statistical results showed that the average molting cycle of the control group was 14.5 days, while that of the experimental group was 16.5 days. However, the morphological structure of ovarian tissue did not abnormal change. Therefore, the results of this study suggest that Mn-CYP302A1 can promote ovarian maturation and molting in female M. nipponense.

RNAi Analysis of Potential Functions of Cyclin B3 in Reproduction of Male Oriental River Prawns (Macrobrachium nipponense)

Cyclin B3 (CycB3) is involved in the metabolic pathway of the cell cycle, playing essential roles in the regulation of cell proliferation and mitosis. CycB3 is also predicted to be involved in the reproduction of male oriental river prawns (Macrobrachium nipponense). In this study, the potential functions of CycB3 in M. nipponense were investigated using quantitative real-time PCR, RNA interference, and histological observations. The full-length DNA sequence of CycB3 in M. nipponense was 2147 base pairs (bp) long. An open reading frame of 1500 bp was found, encoding 499 amino acids. A highly conserved destruction box and two conserved cyclin motifs were found in the protein sequence of Mn-CycB3. Phylogenetic tree analysis revealed that this protein sequence was evolutionarily close to that of CycB3s of crustacean species. Quantitative real-time PCR analysis results suggested that CycB3 was involved in the process of spermiogenesis, oogenesis, and embryogenesis in M. nipponense. RNA interference analysis showed that CycB3 had a positive regulatory relationship with insulin-like androgenic gland hormone (IAG) in M. nipponense. In addition, sperm were rarely observed in the testis of double-stranded CycB3-injected prawns after 14 days of treatment, and sperm abundance was dramatically lower than that in the double-stranded GFP-injected prawns on the same day. This result indicated that CycB3 can regulate the testis reproduction in M. nipponense through inhibiting the IAG expressions. Overall, these results indicated that CycB3 plays essential roles in the regulation of male reproduction in M. nipponense, which may promote the studies of male reproduction in other crustacean species.

Function Analysis of Cholesterol 7-Desaturase in Ovarian Maturation and Molting in Macrobrachium nipponense: Providing Evidence for Reproductive Molting Progress

The Cholesterol 7-desaturase gene plays an important role in insect ecdysone synthesis, but its role in ovarian development has not been reported. In this study, characteristics and the phylogenetic relationship of Cholesterol 7-desaturase were identified by bioinformatics. qPCR showed that the Mn-CH7D gene was highly expressed in the ovary, which was much higher than that in other tissues, and the expression level of Mn-CH7D reached the highest level at the third stage of the ovarian development stage (O-III). During embryonic development, the Mn-CH7D gene expression was highest in the zoea stage. The function of the Mn-CH7D gene was explored by RNA interference. The experimental group was injected with Mn-CH7D dsRNA through the pericardial cavity of M. nipponense, while the control group was injected with the same volume of dsGFP. Statistical analysis of gonadal development and GSI calculation showed that the silencing of Mn-CH7D resulted in the suppression of gonadal development. In addition, the molting frequency of the experimental group was significantly lower than that of the control group during the second molting cycle after silencing Mn-CH7D. On the seventh day after silencing, ecdysone content in the experimental group was significantly reduced. These results demonstrated that the Mn-CH7D gene played a dual role in ovarian maturation and molting of M. nipponense.

Hepatopancreas Proteomic Analysis Reveals Key Proteins and Pathways in Regulatory of Ovary Maturation of Macrobrachium nipponense

A TMT-based (Tandem Mass Tag) liquid chromatography-tandem mass spectrometry (LC-MS/MS) proteomics approach was employed to explore differentially expressed proteins (DEPs) and KEGG pathways in hepatopancreas of 5 ovary stages. In total, 17,999 peptides were detected, among which 3395 proteins were identified. Further analysis revealed 26, 24, 37, and 308 DEPs in HE-I versus HE-II, HE-II versus HE-Ⅲ, HE-Ⅲ versus HE-Ⅳ, and HE-Ⅳ versus HE-Ⅴ, respectively (HE-I, HE-II, HE-III, HE-IV, and HE-V means hepatopancreas sampled from ovary stage I to V.). Gene ontology (GO) analysis indicated that DEPs were significantly enriched in "catalytic activity", "metabolic process", and "cell" of 4 comparison groups in turn. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment results showed that in hepatopancreas, as the ovaries developed to maturation, carbohydrate metabolism, lipid metabolism, amino acid metabolism, and lysosome played important roles in turn. The mRNA expression of 15 selected DEPs were consistent with proteome results by qPCR analysis. Further mRNA expression investigation results suggested 4 proteins (fatty acid-binding protein, NPC intracellular cholesterol transporter 1, Serine hydroxymethyltransferase, and Crustapin) were involved in ovary maturation. These results enhance the understanding of the regulatory role of hepatopancreas in M. nipponense ovary maturation and provide new insights for understanding the crustacean regulation mechanisms.

CRISPR/Cas9 establishment-mediated targeted mutagenesis in Macrobrachium nipponense

Introduction: CRISPR/Cas9 is a gene-editing technology which could specifically cleave dsDNA and induce target gene mutation. CRISPR/Cas9 has been widely used in gene functional studies in many fields, such as medicine, biology, and agriculture due to its simple design, low cost, and high efficiency. Although it has been well developed in model fish and freshwater fish for gene function analysis, it is still novel in the studies dealing with economic crustacean species. Methods: In this study, we established a CRISPR/Cas9 system based on microinjection for M. nipponense, an important economic crustacean aquaculture species. The vitellogenin (Vg) gene and the eyeless (Ey) gene were selected as the targeted genes for mutation. Two sgRNAs were designed for Mn-Vg and Mn-Ey gene editing, respectively. Results and Discussion: For sg-Vg-1, the gastrula survival ratio was 8.69%, and the final hatching ratio was 4.83%. The blastula mutant ratio was 10%, and the hatching individual mutant ratio was 30%. For sg-Vg-2, the gastrula survival ratio was 5.85%, and the final hatching ratio was 3.89%. The blastula mutant ratio was 16.67%, and no mutant sequences were detected in hatching individuals. For sg-Ey-1, the gastrula survival ratio was 6.25%, and the final hatching ratio was 2.34%. The blastula mutant ratio was 10.00%, and the hatching individual mutant ratio was 66.67%. For sg-Ey-2, the gastrula survival ratio was 6.00%, and the final hatching ratio was 2.67%. No mutant sequence was detected in both blastula stage and hatching individuals. There were no significant morphological changes observed in the Mn-Vg group. Two deformed types were detected in sg-Ey-1-injected embryos. An evident developmental delay of the compound eye was detected in Ey-sg1-H1 in the zoea stage. The compound eyes of the Ey-sg1-H2 embryo could not form well-defined spheres, and the whole compound eye appeared to diffuse at the end of the late zoea stage. The establishment of a gene-editing platform based on CRISPR/Cas9 will not only provide an efficient and convenient method for gene function analysis but also provide a powerful tool for molecular-assisted breeding of Macrobrachium nipponense.

Characterization and functional analysis of pax3 in body color transition of polychromatic Midas cichlids (Amphilophus citrinellus)

As the representative genetic and economic trait of ornamental fish, skin color has a strong impact on speciation and adaptation. However, the genetic basis of skin color pigmentation, differentiation and change is still not understood. The Midas cichlid fish with three typical body color transition stages of "black-gray‑gold" is an ideal model system for investigating the formation and change of fish body color. In this study, to investigate the regulatory role of the pair box 3 (pax3) gene in the early body color fading process of Midas cichlids, the complete cDNA sequence (3513 bp) of pax3 was successfully isolated from Midas cichlids (Amphilophus Citrinellus), and found to encode polypeptides of 491 amino acids. Expression patterns of the pax3 gene in tissues of Midas cichlids during different periods, including embryonic development and body color fading stages were detected by quantitative real-time PCR. The qRT-PCR analysis showed that pax3 was expressed in all tissues of adult fish, with a higher expression level in muscle and skin. The highest expression level in muscle tissue was significantly higher than that in other tissues (P < 0.05). During embryonic development, the expression tendency of pax3 was first increased and then decreased. In the three typical stages of early skin color fading from black to gold, pax3 expression in skin, caudal fin and scales all showed a downward trend. The expression level in the black stage was significantly higher than that in other stages (P < 0.05). Positive signal of pax3 protein was detected in the three typical skin color conversion stages, and the highest positive signal intensity was detected in the black stage, which was consistent with qRT-PCR results. After pax3 RNA interference, pax3 and the downstream genes mitf and tyr all decreased, while dct mRNA expression increased in the skin of fish. Western blotting also showed a decrease in pax3 protein concentration. Those results suggest that pax3 plays an important role in skin color formation, distribution and change in Midas cichlids through the melanogenesis pathway.

RNA interference analysis of potential functions of cyclin A in the reproductive development of male oriental river prawns (Macrobrachium nipponense)

Cyclin A (CycA) plays essential roles in regulating multiple steps of the cell cycle, and it affects gonad development in mammals and invertebrates. Previous RNA interference (RNAi) analysis revealed that knocking-down the expression of CycA in female oriental river prawns (Macrobrachium nipponense) inhibited ovarian development. CycA was also predicted to have regulatory roles in reproductive development of male M. nipponense based on significant changes of Mn-CycA expression after eyestalk ablation. The goal of this study was to investigate the potential functions of CycA in the reproductive development of male M. nipponense using RNAi and histological observations. Quantitative real-time PCR analysis revealed that both single-side and double-side eyestalk ablation stimulated the expressions of Mn-CycA, and the expression was higher in prawns with double-side eyestalk ablation (p < 0.05). Mn-CycA expression was significantly higher in the testis and androgenic gland during the reproductive season than during the non-reproductive season (p < 0.05). In the RNAi analysis, Mn-CycA expression significantly decreased after prawns were injected with dsCycA, and the expression of insulin-like androgenic gland hormone (Mn-IAG) also decreased as Mn-CycA expression decreased. This result indicated that CycA positively regulated the expression of IAG in M. nipponense. Histological observations revealed that the number of sperm decreased dramatically to <5% of the total cells in the testis of the dsCycA-treated group compared to that of control group on day 14, indicating that knockdown of Mn-CycA expression inhibited testis development by affecting the expression of Mn-IAG in M. nipponense. These results highlighted the functions of CycA in male reproductive development of M. nipponense, which can be applied to future studies of male reproduction in other crustacean species.

RNA Interference Analysis of the Functions of Cyclin B in Male Reproductive Development of the Oriental River Prawn (Macrobrachium nipponense)

Cyclin B (CycB) plays essential roles in cell proliferation and promotes gonad development in many crustaceans. The goal of this study was to investigate the regulatory roles of this gene in the reproductive development of male oriental river prawns (Macrobrachium nipponense). A phylo-genetic tree analysis revealed that the protein sequence of Mn-CycB was most closely related to those of freshwater prawns, whereas the evolutionary distance from crabs was much longer. A quantitative PCR analysis showed that the expression of Mn-CycB was highest in the gonad of both male and female prawns compared to that in other tissues (p < 0.05), indicating that this gene may play essential roles in the regulation of both testis and ovary development in M. nipponense. In males, Mn-CycB expression in the testis and androgenic gland was higher during the reproductive season than during the non-reproductive season (p < 0.05), implying that CycB plays essential roles in the reproductive development of male M. nipponense. An RNA interference analysis revealed that the Mn-insulin-like androgenic gland hormone expression decreased as the Mn-CycB expression decreased, and that few sperm were detected 14 days after the dsCycB treatment, indicating that CycB positively affects testis development in M. nipponense. The results of this study highlight the functions of CycB in M. nipponense, and they can be applied to studies of male reproductive development in other crustacean species.

Effects of Low Temperature on Antioxidant and Heat Shock Protein Expression Profiles and Transcriptomic Responses in Crayfish (Cherax destructor)

Low temperature is a critical factor restricting the growth and survival of aquatic animals, but research on the mechanism of response to low temperature in Cherax destructor is limited. C. destructor is one of the most important freshwater crustaceans with strong adaptability in Australia, and it has been commercialized gradually in recent years. Here, growth indicators, antioxidant parameters, anti-stress gene expression, and transcriptome sequencing were used on crayfish following 8 weeks of low-temperature acclimation. The results showed that weight gain, length gain, and molting rates decreased as the temperature decreased. The activity of antioxidant enzymes decreased, while the content of antioxidant substances and the expression of anti-stress genes increased. Transcriptome sequencing identified 589 differentially expressed genes, 279 of which were upregulated and 310 downregulated. The gene functions and pathways for endocrine disorders, glucose metabolism, antioxidant defense, and immune responses were identified. In conclusion, although low-temperature acclimation inhibited the basal metabolism and immune ability of crayfish, it also increased the antioxidant substance content and anti-stress-gene expression to protect the organism from low-temperature damage. This study provided molecular insights into the study of low-temperature responses of low-temperature-tolerant crustacean species.

Hepatopancreas transcriptome analyses provide new insights into the molecular regulatory mechanism of fast ovary maturation in Macrobrachium nipponense

Background

Macrobrachium nipponense is an economically and ecologically important freshwater prawn that is widely farmed in China. In contrast to other species of marine shrimp, M. nipponense has a short sexual maturity period, resulting in not only high stocking densities, but also a reduced survival rate and increased risk of hypoxia. Therefore, there is an urgent need to study the molecular mechanisms underlying fast ovary maturation in this species.

Results

Comparative transcriptome analysis was performed using hepatopancreatic tissue from female M. nipponense across five ovarian maturation stages to explore differentially expressed genes and pathways involved in ovarian maturation. In total, 118.01 Gb of data were generated from 15 transcriptomes. Approximately 90.46% of clean reads were mapped from the M. nipponense reference genome. A comprehensive comparative analysis between successive ovarian maturation stages generated 230–5814 differentially expressed genes. Gene Ontology (GO) enrichment was highly concentrated in the “biological process” category in all four comparison groups, and mainly focused on energy synthesis and accumulation, energy decomposition and transport. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment results showed that, among 20 significantly enriched KEGG pathways, nine were involved in the synthesis, degradation, and metabolism of carbohydrates, lipids, and other nutrient intermediates, suggesting that the hepatopancreas has an important role in energy supply during ovarian maturation. Furthermore, the “Insect hormone biosynthesis” pathway was found to have a dominant role in the development of the ovary from immaturity to maturity, supporting the hypothesis that ecdysteroid- and juvenile hormone-signaling pathways have an important role in hepatopancreas regulation of ovarian maturation.

Conclusion

Taken together, this study sheds light on the role of the hepatopancreas in the molecular regulation of ovary maturation in M. nipponense. The present study provided new insights for understanding the mechanisms of reproductive regulation in crustaceans.

Cathepsin D Plays a Vital Role in Macrobrachium nipponense of Ovary Maturation: Identification, Characterization, and Function Analysis

The oriental river prawn Macrobrachium nipponense is an economically important aquacultural species. However, its aquaculture is negatively impacted by the rapid sexual maturation of female M. nipponense. The fast sexual maturation produces a large number of offspring which leads to a reduction in resilience, a low survival rate, and an increased risk of hypoxia, this in turn, seriously affects the economic benefits of prawn farming. Cathepsin D is a lysosomal protease involved in the ovarian maturation of M. nipponense. In the current study, the cDNA of the gene encoding cathepsin D (Mn-CTSD) was cloned from M. nipponense. The total length was 2391 bp and consisted of an open reading frame (ORF) of 1158 bp encoding 385 amino acids. Sequence analysis confirmed the presence of conserved N-glycosylation sites and characteristic sequences of nondigestive cathepsin D. The qPCR analysis indicated that Mn-CTSD was highly expressed in all tissues tested, most significantly in the ovaries, whereas in situ hybridization showed that expression occurred mainly in oocyte nuclei. Analysis of its expression during development showed that Mn-CTSD peaked during the O-IV stage of ovarian maturation. For the RNAi interference experiment, female M. nipponense specimens in the ovary stage I were selected. Injection of Mn-CTSD double-stranded (ds)RNA into female M. nipponense decreased the expression of Mn-CTSD in the ovaries, such that the Gonad Somatic Index (GSI) of the experimental group was significantly lower than that of the control group (1.79% versus 4.57%; p < 0.05). Ovary development reached the O-III stage in 80% of the control group, compared with 0% in the experimental group. These results suggest that Mn-CTSD dsRNA inhibits ovarian maturation in M. nipponense, highlighting its important role in ovarian maturation in this species and suggesting an approach to controlling ovarian maturation during M. nipponense aquaculture.

Molecular and functional analysis of the microphthalmia-associated transcription factor (mitf) gene duplicates in red tilapia

In vertebrates, the microphthalmia-associated transcription factor (mitf) is at the hub of the melanin synthesis regulation network. However, little information is known about its molecular characterization, expression, location, or function in skin color differentiation and variation of red tilapia. The full-length cDNA sequences (1977 bp and 1999 bp) of mitfa and mitfb, encoding polypeptides of 491 and 514 amino acids, were effectively identified from red tilapia in this study. The Mitfa and Mitfb sequences of red tilapia clustered first with O. aureus, then with other teleost fish, according to phylogenetic analysis. Mitfa and mitfb mRNA were highly expressed in the brain, dorsal skin and eye tissues using quantitative real-time PCR. The mRNA expressions of mitfa and mitfb were the highest in the cleavage stage during the early development of red tilapia. Among three different colors of red tilapia, the expression levels of mitfa and mitfb were highest in the PB (pink with scattered black spots) dorsal skin. After overwintering, the mitfa and mitfb mRNA expressions were high in the dorsal skin of PB (color changed from pink to black). Mitfa and mitfb were mostly found in the epidermal layer of the dorsal skin, according to in situ hybridization (ISH) analysis. After injecting mitf-dsRNA duplicates along the tail vein of red tilapia, the activity of tyrosinase and the level of melanin in the dorsal skin both decreased significantly. The mRNA expressions of mitfa and its downstream genes (tyrb, tyrp1a and dct) decreased, whereas the mRNA expression of mitfb increased after mitfa-dsRNA injection. The mRNA expressions of mitfb, tyrb, tyrp1a and dct decreased, whereas the mRNA expression of mitfa increased after injecting mitfb-dsRNA. These findings suggest that mitf gene duplicates may play an important role in red tilapia skin color differentiation and variation via the melanogenesis pathway.

Identification of ribosomal protein L24 (RPL24) from the oriental river prawn, Macrobrachium nipponense, and its roles in ovarian development

Ribosomal proteins exhibit various extraribosomal functions in addition to their roles in protein synthesis. In this study, complementary DNA (cDNA) of ribosomal protein L24 in Macrobrachium nipponense (MnRPL24) was isolated, and its role in ovarian development was investigated using quantitative real-time PCR (qPCR), immunohistochemistry (IHC), RNA interference (RNAi) and histological observations. The complete cDNA of MnRPL24 is 564 base pairs (bps) and contains a 486 bp open reading frame (ORF) encoding 162 amino acids (aas). The highest expression level of MnRPL24 among eight tissues was found in the ovary, specifically in the stage I ovary. The MnRPL24 protein existed in the cytoplasm and nucleus of developing oocytes, and also existed in the cytoplasm of follicle cells in developing ovaries. After MnRPL24 knockdown by RNAi, the expression levels of vitellogenin (Vg), vitellogenin receptor (Vgr), cyclin-dependent kinase 2 (Cdc2) and M-phase cyclin (Cyclin B) genes and the gonadsomatic index (GSI) did not show the typical trend of gradually elevation with ovarian development and finally decrease in the later stage of ovarian cycle. Moreover, the oviposition rate (OR) was downregulated, and oocyte development was delayed after MnRPL24 knockdown. After eyestalk ablation, the MnRPL24 expression level was considerably elevated in the initial stages and decreased in the late stage of the ovarian development cycle. This investigation illustrates a possible regulatory role of MnRPL24 in the ovarian development of M. nipponense, and MnRPL24 may act as a stimulator of early ovarian development.

Microplastics and bisphenol A hamper gonadal development of whiteleg shrimp (Litopenaeus vannamei) by interfering with metabolism and disrupting hormone regulation

Gonadal development is a prerequisite for the reproductive success of an organism, and might be affected by environmental factors such as emergent pollutants. Although marine crustaceans are threatened by ubiquitous emergent pollutants such as microplastics (MPs) and bisphenol A (BPA) under realistic scenarios, studies about the impacts of these pollutants on the gonadal development of crustacean species are rare. In this study, the effects of MPs and BPA, alone or in combination, on gonadal development were investigated in whiteleg shrimp (Litopenaeus vannamei). The results obtained demonstrated that whiteleg shrimp exposed to MPs and BPA had significantly smaller gonad-somatic index (GSI) and an obvious delay in the gonad developmental stage. In addition, exposure of whiteleg shrimp to pollutants tested resulted in a reduction in the oxygen consumption rate, elevation in the ammonia excretion rate, decline in the O: N ratio, and downregulation in the expression of metabolism-related genes, indicating significant disruptions of shrimp metabolism by the pollutants. Furthermore, in addition to a few exceptions, both the in vivo contents of gonadal development-related hormones (GIH and MIH) and the expression of genes encoding regulatory hormones (GIH, MIH, and CHH) were upregulated by the exposure of whiteleg shrimp to the pollutants investigated, suggesting a significant obstruction of endocrine regulation. Moreover, MP-BPA coexposure was shown to be more toxic to whiteleg shrimp than the corresponding single exposures and significantly greater amount of BPA accumulated in the gonads (both testis and ovaries) of shrimp with the coexistence of MPs, which may be caused by the Trojan horse effect and summation of the toxic impacts on common targets. In general, the data obtained in this study demonstrated that MPs and BPA at environmentally realistic concentrations significantly inhibited the gonadal development of whiteleg shrimp probably by interfering with metabolism and disrupting endocrine regulation.

Genomic structure, expression, and functional characterization of the Fem-1 gene family in the redclaw crayfish, Cherax quadricarinatus

The Fem-1 (Feminization-1) gene, encoding an intracellular protein with conserved ankyrin repeat motifs, has been proven to play a key role in sex differentiation in Caenorhabditis elegans. In the present study, three members of the Fem-1 gene family (designating Fem-1A, Fem-1B, and Fem-1C, respectively) were cloned and characterized in the redclaw crayfish, Cherax quadricarinatus. Sequence analysis showed that all three Fem-1 genes contained the highly conserved ankyrin repeat motifs with variant repeat numbers, which shared similarity with other reported crustaceans. In addition, a phylogenetic tree revealed that the Fem-1 proteins from C. quadricarinatus were clustered with the crustacean Fem-1 homologs, and had the closest evolutionary relationship with Eriocheir sinensis. Quantitative real-time PCR (qRT-PCR) results demonstrated that Fem-1B exhibited a significant higher expression abundance in the ovary than in other tissues. In addition, a regular mRNA expression pattern of the Fem-1B gene appeared in the reproductive cycle of ovarian development. Furthermore, RNA interference experiments were employed to investigate the role of Fem-1B in ovarian development. Moreover, knockdown of Fem-1B by RNAi decreased the expression of VTG in the ovaries and hepatopancreas. In summary, this study pointed out that Fem-1B was involved in the sex differentiation process through regulating VTG expression in C. quadricarinatus, and provided new insights into the role of Fem-1B in ovary development.

Identification of potential functions of polo-like kinase 1 in male reproductive development of the oriental river prawn (Macrobrachium nipponense) by RNA interference analysis

Polo-like kinase 1 (Plk1) has multiple functions in the cell cycle, including in the maturation of centrosomes during the G2/M transition, the separation of centrosomes, and the activation of cyclin-dependent kinase 1 expression and spindle assembly. In this study, we investigated the potential regulatory roles of Plk1 in the reproductive development of the male oriental river prawn (Machrobrachium nipponense). The full cDNA sequence of Mn-Plk1 was 2360 base pairs long, with an open reading frame of 1836 base pairs encoding 611 amino acids. Protein sequence alignment identified a conserved serine/threonine kinase domain and two Polo-boxes. Phylogenetic tree analysis revealed that Mn-Plk1 had the closest evolutionary distance with Plk1s of freshwater prawns and then with those of crustacean species, whereas the evolutionary distance with mollusks was much more distant. Quantitative PCR analysis predicted that Mn-Plk1 plays essential roles in the regulation of gonad development. RNA interference analysis and histological observations showed that expression of insulin-like androgenic gland hormone decreased as the expression of Mn-Plk1 decreased, and fewer than 5% of cells were sperm cells at day 14 in the dsPlk1 injected prawns. This result indicated that Plk1 positively regulated testis development in M. nipponense by affecting the expression of this hormone. Our results highlight the functions of Plk1 in M. nipponense and provide valuable information that can be applied to establish artificial techniques to regulate testis development in this species.

Identification and Characterization of the Pyruvate Dehydrogenase E1 Gene in the Oriental River Prawn, Macrobrachium nipponense

Pyruvate dehydrogenase E1 (PDHE1) is thought to play essential roles in energy metabolism, and a previous study suggested that it also has potential regulatory roles in male sexual development in the oriental river prawn, Macrobrachium nipponense. In this study, we used rapid amplification of cDNA ends, quantitative polymerase chain reaction (qPCR), in situ hybridization, western blotting, RNA interference (RNAi), and histological analyses to assess the potential functions of Mn-PDHE1 in the sexual development of male M. nipponense. The full cDNA sequence of Mn-PDHE1 was 1,614 base pairs long, including a 1,077 base pair open reading frame that encodes 358 amino acids. qPCR analysis revealed the regulatory functions of PDHE1 in male sexual development in M. nipponense and in the metamorphosis process. In situ hybridization and western blot results indicated that PDHE1 was involved in testis development, and RNAi analysis showed that PDHE1 positively regulated the expression of insulin-like androgenic gland factor in M. nipponense. Compared with the cell types in the testes of control prawns, histological analysis showed that the number of sperm was dramatically lower after test subjects were injected with Mn-PDHE1 dsRNA, whereas the numbers of spermatogonia and spermatocytes were higher. Sperm constituted only 1% of cells at 14 days after injection in the RNAi group. This indicated that knockdown of the expression of PDHE1 delayed testis development. Thus, PDHE1 has positive effects on male sexual development in M. nipponense. This study highlights the functions of PDHE1 in M. nipponense and its essential roles in the regulation of testis development.

Identification of candidate genes from androgenic gland in Macrobrachium nipponense regulated by eyestalk ablation

The eyestalk of crustaceans, such as Macrobrachium nipponense, contains many neurosecretory hormones affecting the process of reproduction, molting, metabolism of glucose, and other functions. In this study, important metabolic pathways and candidate genes involved in male sexual development were selected from M. nipponense. The methodology involved performing long-read and next generation transcriptome sequencing of genes from the androgenic gland after eyestalk ablation. qPCR analysis revealed that the mRNA expression of Mn-IAG was significantly increased after ablation of both the single-side (SS) and double-side (DS) eyestalk, compared with the control group (CG). The long-read transcriptome generated 49,840 non-redundant transcripts. A total of 1319, 2092 and 4351 differentially expressed genes (DEGs) were identified between CG versus SS, SS versus DS and CG versus DS, respectively. These data indicated that ablation of the double-sided eyestalk played stronger regulatory roles than the single-side ablation on male sexual development in M. nipponense. This was consistent with the qPCR analysis. Cell Cycle, Cellular Senescence, Oxidative Phosphorylation, Glycolysis/Gluconeogenesis and Steroid Hormone Biosynthesis were the primary enriched metabolic pathways in all three comparisons, and the important genes from these metabolic pathways were also selected. qPCR permitted secondary confirmation of ten DEGs identified through RNA-seq. RNAi-mediated silencing analyses of Hydroxysteroid dehydrogenase like 1 (HSDL1) revealed that HSDL1 has a positive regulatory effect on testes development. This study provides valuable insight into male sexual development in M. nipponense, including metabolic pathways and genes, paving the way for advanced studies on male sexual development in this species and in other crustaceans.

Transcriptome analysis of hepatopancreas from different living states oriental river prawn (Macrobrachium nipponense) in response to hypoxia

As an important economical freshwater prawn, Macrobrachium nipponense has difficulty with adapting to hypoxia. In this study, comparative transcriptome analysis was used for the first time to explore the differences between different living states of Macrobrachium nipponense under hypoxia. A total of 94.22 Gb clean reads were obtained and assembled into 54,688 unigenes. A total of 224, 266, and 750 differently expressed genes were found in the comparison of the control and death groups, the control and moribund groups, and the control and survived groups, respectively. Three signal pathways closely related to hypoxia were found by enriching of the signal pathways in three comparison groups. In addition, much attention was focused on the differential genes in these pathways. Oxidative stress related genes, such as 70 kDa heat shock protein, phosphoenolpyruvate carboxykinase and cyclooxygenase were differentially expressed in different comparisons. After comparing with previous studies, cyclooxygenase was found to be an important hypoxia-related gene that is fully involved in the hypoxic response. Interestingly, two new genes with no Nr annotation were found in this manuscript. This manuscript will enrich our understanding of oxidative stress response to hypoxia and provide a theoretical basis for the subsequent solution of apoptosis caused by hypoxia.

Identification and Characterization of the Succinate Dehydrogenase Complex Iron Sulfur Subunit B Gene in the Oriental River Prawn, Macrobrachium nipponense

Previous studies have revealed that SDHB has potential functions in the male sexual differentiation and development in M. nipponense through providing ATP. In this study, the functions of Mn-SDHB were further analyzed in depth using quantitative polymerase chain reaction (qPCR), in situ hybridization, western-blot, and RNA interference (RNAi), combined with the histological observations. The full-genome sequence of Mn-SDHB was 54,608 bp at Chromosome 34, including 7 introns and 6 exons. The full-length cDNA sequence of Mn-SDHB was 1,268 base pairs (bp) long with an open reading frame of 807 bp, encoding for 268 amino acids. The highest expression level of Mn-SDHB in different tissues was observed in the testis, and male prawns at post-larval developmental stage 25 during different developmental stages, indicating that SDHB was potentially involved in the male sexual development in M. nipponense. In situ hybridization and western-blot analysis indicated that SDHB plays essential roles in the testis development. The in situ hybridization analysis also implies the potential roles of Mn-SDHB in ovarian development. The expressions of Mn-IAG were decreased after Mn-SDHB dsRNA injection, indicating SDHB has the positive regulatory effects on IAG in M. nipponese. Thus, SDHB was involved in the mechanism of the male sexual development. The testis development was inhibited, and sperms were rarely observed after 10 days of Mn-SDHB dsRNA injection, indicating SDHB has positive effects on the male sexual development in M. nipponense. This study highlights the functions of SDHB in M. nipponense, which provide new insights for the future studies of the male sexual development in other crustacean species.

Transcriptome analysis of five ovarian stages reveals gonad maturation in female Macrobrachium nipponense

Background

Macrobrachium nipponense is an economically important species of freshwater shrimp in China. Unlike other marine shrimps, the ovaries in adult female M. nipponense can mature rapidly and periodically during the reproductive period, but the resulting high stocking densities and environmental deterioration can negatively impact the harvest yield and economic benefits. To better understand ovary development in female M. nipponense, we performed systematic transcriptome sequencing of five different stages of ovarian maturation.

Results

We obtained 255,966 Gb of high quality transcriptome data from 15 samples. Of the 105,082 unigenes that were selected, 30,878 were successfully annotated. From these unigenes, we identified 17 differentially expressed genes and identified three distinct gene expression patterns related to different biological processes. We found that cathepins, legumains, and cystatin were enriched in the lysosome pathway, and they are related to vitellogenin hydrolysis. Additionally, we found that myosin heavy chain 67 participated in oocyte excretion.

Conclusions

We provide the first detailed transcriptome data relating to the ovarian maturation cycle in M. nipponense. Our results provide important reference information about the genomics, molecular biology, physiology, and population genetics of M. nipponense and other crustaceans. It is conducive to further solve the problem of M. nipponense rapid ovarian maturation from the aspects of energy supply and cell division.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07737-5.

Transcriptome Profiling Analysis of the Testis After Eyestalk Ablation for Selection of the Candidate Genes Involved in the Male Sexual Development in Macrobrachium nipponense

The eyestalk of crustacean species secretes many hormones, affecting the process of reproduction, molting, metabolism of glucose, and other functions in crustaceans. In this study, important metabolic pathways and candidate genes involved in the male sexual development were identified through performing the transcriptome profiling analysis of the testis after the ablation of eyestalk from Macrobrachium nipponense. The histological observations revealed that the testis development became vigorous after eyestalk ablation, indicating that the hormones secreted by the eyestalk have negative effects on the testis development in M. nipponense. Transcriptome profiling analysis revealed that 1,039, 1,226, and 3,682 differentially expressed genes (DEGs) were identified between normal prawns (CG) vs single-side eyestalk ablation prawns (SS), SS vs double-side eyestalk ablation prawns (DS), and CG vs DS, respectively, indicating that the ablation of double-side eyestalk has more significant regulatory roles on male sexual development than that of single-side ablation, which was consistent with the histological observations. Lysosome, Apoptosis, Glycolysis/Gluconeogenesis, and Insulin signaling pathway were the main enriched metabolic pathways in all of these three comparisons, and the important genes from these metabolic pathways were also selected. The qPCR verifications of 10 DEGs from these metabolic pathways were the same as those of RNA-seq. The qPCR, in situ hybridization, and RNA interference analysis of Mn-NF_kBα revealed that NF_kBα has a positive regulatory effect on testis development. This study provided new insights on male sexual development in M. nipponense, promoting the studies on male sexual development in other crustaceans as well.

Function analysis and molecular characterization of cyclin A in ovary development of oriental river prawn, Macrobrachium nipponense

Macrobrachium nipponense has the characteristics of fast ovarian development cycle, which leads to the coexistence of multiple generations, the reduction of commodity specifications and the low economic benefit. Therefore, the study on the mechanism of ovarian development is of great significance to the development of industry. Cyclin A (CycA)is a key gene regulating ovarian development in vertebrates, but little information was available for its function in crustaceans. In this study, the full-length cDNA of Mn-CycA was obtained from the ovary. The full-length cDNA (2033 bp) with an open reading frame of 1368 bp, encoded a 456-amino acid protein. qRT-PCR revealed tissue-specific expression pattern of Mn-CycA, with abundant expression in the ovary. Results in different developmental stages of ovary indicated that Mn-CycA expression is positively correlated with ovarian maturation. qRT-PCR In different developmental stages, the expression of Mn-CycA mRNA gradually increased during the embryonic stage and decreased significantly on the first day of the hatching stage. At the 25th day of the metamorphosis stage, the expression level of Mn-CycAmRNA in female shrimp was 3.5 times higher than that in male shrimp, which may be related to the proliferation of oogonia and the formation of oocytes. In situ hybridization (ISH) of ovary showed Mn-CycA was examined in all stages and was mainly located in oogonia and oocytes. Compared with the control group, the obvious change of gonad somatic index (GSI) proved that injection of Mn-CycA dsRNA could delay the ovarian development cycle, which provided strong evidence for the involvement of Mn-CycA in ovarian maturation and oogenesis, and expanded a new perspective for studying the fast ovarian development cycle in M. nipponense.

Suppression of a Novel Vitellogenesis-Inhibiting Hormone Significantly Increases Ovarian Vitellogenesis in the Black Tiger Shrimp, Penaeus monodon

In this study, a novel Crustacean Hyperglycemic Hormone-type II gene (CHH-type II) was identified and biologically characterized in a shrimp, Penaeus monodon. Based on its structure and function, this gene was named P. monodon vitellogenesis-inhibiting hormone (PemVIH). The complete cDNA sequence of PemVIH consisted of 1,022 nt with an open reading frame (ORF) of 339 nt encoding a polypeptide of 112 amino acids. It was classified as a member of the CHH-type II family based on conserved cysteine residues, a characteristically positioned glycine residue, and the absence of CHH precursor-related peptide (CPRP) domain. The deduced mature PemVIH shared the highest sequence similarities with giant river prawn sinus gland peptide A. Unlike P. monodon gonad-inhibiting hormone (PemGIH), PemVIH was expressed only in the brain and ventral nerve cord, but not the eyestalks. Whole mount immunofluorescence using a newly generated PemVIH antiserum detected positive signals in neuronal cluster 9/11 and 17 of the brain, commissural ganglion (CoG), and neuronal clusters of ventral nerve cord. The presence of PemVIH-positive neurons in CoG, a part of stomatogastric nervous system, suggested a potential mechanism for crosstalk between nutritional and reproductive signaling. The role of PemVIH in vitellogenesis was evaluated using RNA interference technique. Temporal knockdown of PemVIH in female subadults resulted in a 3-fold increase in ovarian vitellogenin expression, suggesting an inhibitory role of PemVIH in vitellogenesis. This study provided novel insight into the control of vitellogenesis and additional strategies for improving ovarian maturation in P. monodon without the current harmful practice of eyestalk ablation.

Molecular Characterization of a Novel Cathepsin L in Macrobrachium nipponense and Its Function in Ovary Maturation

Cathepsin L genes, which belonged to cysteine proteases, were a series of multifunctional protease and played important roles in a lot of pathological and physiological processes. In this study, we analyzed the characteristics a cathepsin L (named Mn-CL2) in the female oriental river prawn, Macrobrachium nipponense which was involved in ovary maturation. The Mn-CL2 was1,582 bp in length, including a 978 bp open reading frame that encoded 326 amino acids. The Mn-CL2 was classified into the cathepsin L group by phylogenetic analysis. Real-time PCR (qPCR) analysis indicated that Mn-CL2 was highly expressed in the hepatopancreas and ovaries of female prawns. During the different ovarian stages, Mn-CL2 expression in the hepatopancreas and ovaries peaked before ovarian maturation. In situ hybridization studies revealed that Mn-CL2 was localized in the oocyte of the ovary. Injection of Mn-CL2 dsRNA significantly reduced the expression of vitellogenin. Changes in the gonad somatic index also confirmed the inhibitory effects of Mn-CL2 dsRNA on ovary maturation. These results suggest that Mn-CL2 has a key role in promoting ovary maturation.

Three neuroparsin genes from oriental river prawn, Macrobrachium nipponense, involved in ovary maturation

In this study, we identified three neuroparsin (NP) genes in Macrobrachium nipponense: Mn-NP1, Mn-NP2, and Mn-NP3, encoding 99, 100, and 101 amino acid proteins, respectively. Multiple sequence alignments showed that these genes contained 12 cysteine residues, of which 11 were at conserved positions. The total sequence identity between the genes was 47.5%, and they showed a high degree of sequence identity (> 54% similarity) with other crustacean genes. Phylogenetic tree analysis showed that Mn-NPs were clustered at different branches, indicating that Mn-NPs may have different functions. Tissue distribution data revealed that the three genes were present in males and females during the breeding and nonbreeding season, but their expression patterns differed. Mn-NP1 was highly expressed in the breeding season, in the male testis, and highly expressed in the nonbreeding season, in the female ovary. Mn-NP3 exhibited biased female expression in the breeding and nonbreeding season, with dominant expression in the ovary. All Mn-NPs were detected during embryo development, but with different expression patterns. These data indicated that Mn-NP1 may function during embryonic development, and that Mn-NP2 may be expressed during early embryo cell division, and late larval development. Mn-NP3 expression patterns reflected maternal inheritance, and may be associated with ovarian maturation. These expression data suggested that Mn-NP1 and Mn-NP2 are negatively correlated with ovarian development, with inhibition roles during this development. Mn-NP3 may be involved in vitellogenesis.

Molecular cloning, characterization and functional analysis of two neuropeptide F genes from the oriental river prawn (Macrobrachium nipponense)

In invertebrates, neuropeptide F (NPF) has many regulatory functions, similar to NPY, its homologous peptide. In this study, two neuropeptide F genes were identified in Macrobrachium nipponense: Mn-NPF1 and Mn-NPF2. Mn-NPF2 shared the same amino acid sequence with Mn-NPF1, except for a 37 amino acid insert in the middle of the NPF region. The quantitative-PCR (qPCR) results indicated that Mn-NPF1 expression was positively correlated with ovarian maturation, whereas Mn-NPF2 had opposing expression patterns. Both Mn-NPFs were poorly expressed at early embryonic stages, but enhanced expression levels were observed up to day 10 after hatching, when the gonads began to differentiate. Ovary in situ hybridization (ISH) analyses showed that both Mn-NPFs were present at all stages, but were differentially localized to distinct compartments. Temperature gradient studies showed that both Mn-NPFs were implicated in the seasonal regulation of reproduction. A double-stranded (ds) RNA-Mn-NPF2 injection led to a significant 38.5% increase in the vitellogenin (VG) transcript (P < 0.05). These results demonstrated that Mn-NPF2 plays an important role in inhibiting ovarian maturation.

Roles of gonad-inhibiting hormone in the protandric simultaneous hermaphrodite peppermint shrimp†

To date, the molecular mechanisms of the unique gonadal development mode known as protandric simultaneous hermaphroditism (PSH) are unclear in crustaceans. In this study, cDNA of a gonad-inhibiting hormone (Lv-GIH1) was isolated from the PSH peppermint shrimp Lysmata vittata, and its expression was exclusively found in the eyestalk ganglion. Real-time quantitative polymerase chain reaction (qRT-PCR) revealed that the expression of Lv-GIH1 increased during gonadal development of the functional male stages but decreased significantly at subsequent simultaneous hermaphroditism stage. Further in vitro experiment showed that recombinant GIH1 protein (rGIH1) effectively inhibited Vg expression in the cultured hepatopancreas tissues while the short-term injection of GIH1-dsRNA resulted in reduced expression of Lv-GIH1 and upregulated expression of Vg in the hepatopancreas. Moreover, long-term rGIH1 injection led to significantly reduced expression of Lv-Vg, Lv-VgR, and Lv-CFSH1, subdued growth of oocytes, and feathery setae as a secondary sexual characteristic in females. Interestingly, while germ cells in testicular part were suppressed by rGIH1 injection, the expression of Lv-IAGs showed no significant difference; and long-term GIH1-dsRNA injection results were contrary to those of rGIH1 injection. Taken together, the results of this study indicate that Lv-GIH1 is involved in gonadal development and might also participate in controlling secondary sexual characteristic development in L. vittata by inhibiting Lv-CFSH1 expression.

Transcriptomic characterisation of neuropeptides and their putative cognate G protein-coupled receptors during late embryo and stage-1 juvenile development of the Aotearoa-New Zealand crayfish, Paranephrops zealandicus

We de novo assembled a transcriptome for early life-stages of the Aotearoa-New Zealand crayfish, Paranephrops zealandicus, establishing the first genetic resource for this under-developed aquaculture species and for the Paranephrops genus. Mining of this transcriptome for neuropeptides and their putative cognate G protein-coupled receptors (GPCRs) yielded a comprehensive catalogue of neuropeptides, but few putative neuropeptide GPCRs. Of the neuropeptides commonly identified from decapod transcriptomes, only crustacean female sex hormone and insulin-like peptide were absent from our trinity de novo transcriptome assembly, and also RNA-sequence reads. We identified 63 putative neuropeptide precursors from 43 families, predicted to yield 122 active peptides. Transcripts encoding 26 putative neuropeptide GPCRs were identified but were often incomplete. Putative GPCRs for 15 of the neuropeptides identified here were absent from our transcriptome and RNAseq reads. These data highlight the diverse neuropeptide systems already present at the early development life stages sampled here for P. zealandicus.

Molecular cloning, expression, and in situ hybridization analysis of MnGPx-3 and MnGPx-4 from oriental river prawn, Macrobrachium nipponense, in response to hypoxia and reoxygenation

Glutathione peroxidase (GPx) has been the focus of increased research because of its important role as an antioxidant and in reactive oxygen species (ROS) induced damage repair. Studies on GPxs have relevance with Macrobrachium nipponense because it has poor tolerance to hypoxia in Macrobrachium nipponense. The two subunits named as MnGPx-3 and MnGPx-4 according to the glutathione peroxidase nomenclature system. Both full-length cDNAs were cloned from the hepatopancreas. In this study, we analyzed the expression of two GPxs in Macrobrachium nipponense in response to changes in environmental oxygen. Expression levels of MnGPx-3 and MnGPx-4 indicated that both have strong responses to hypoxia. In situ hybridization showed that MnGPx-3 and MnGPx-4 were located in secretory and storage cells in hepatopancreas. These results suggest that GPx gene is expressed and released by secretory cells and released response to hypoxia. In the gill tissue, however, GPxs are located in blood cells, suggesting that they perform different functions in different tissues or organs. The results of in situ hybridization were consistent with those of quantitative Real-time PCR. This study provides a basis for understanding the oxidative stress response in M. nipponense under hypoxia.

Comparative metabolomics analysis of ovarian developmental stages in Macrobrachium nipponense

Rapid sexual maturity of female Macrobrachium nipponense is a severe problem for the aquaculture industry. To date, there have been only transcriptome studies investigating ovarian development, and studies using other tools, such as metabolomics are lacking. Metabolomics reveals changes in the level of metabolites in tissues in relation to current physiological characteristics, and can yield valuable insight into the growth and development of organisms. In this study, we systematically analyzed 15 samples from five different ovarian developmental stages in M. nipponense to learn more about how metabolites change over reproduction. Gas chromatography/time-of-flight mass spectroscopy revealed an array of different compounds and 83-162 pathways depending on the stage. Furthermore, 89 metabolites and 14 pathways were significantly different across stages. It is hypothesized that N-acetyl-N-formyl-5-methoxykynurenamine, ascorbate, fructose-2,6-bisphosphate, cortexolone and other metabolites that significantly differed by stage are regulated by hormones and are closely related to ovarian development. However, for other metabolites that changed with development, such as cytidine and xanthine, an association with ovarian development has yet to be revealed. Quantitative polymerase chain reaction was used to correlate gene changes to metabolites in the pathway for biosynthesis of plant secondary metabolites. We found that the TCA cycle rate may be the cause of female miniaturization during the reproductive period, and that the control of fatty acid content via aquaculture nutrition may be an exogenous tool for regulatory control of maturation. This study provides a systematic and comprehensive metabolomics analysis of ovarian development in M. nipponense and lays a foundation for addressing the problem of rapid sexual maturity.

Insights on Molecular Mechanisms of Ovarian Development in Decapod Crustacea: Focus on Vitellogenesis-Stimulating Factors and Pathways

Vitellogenesis in crustaceans is an energy-consuming process. Though the underlying mechanisms of ovarian maturation in decapod Crustacea are still unclear, evidence indicates the process to be regulated by antagonistically-acting inhibitory and stimulating factors specifically originating from X-organ/sinus gland (XO/SG) complex. Among the reported neuromediators, neuropeptides belonging to the crustacean hyperglycemic hormone (CHH)-family have been studied extensively. The structure and dynamics of inhibitory action of vitellogenesis-inhibiting hormone (VIH) on vitellogenesis have been demonstrated in several species. Similarly, the stimulatory effects of other neuropeptides of the CHH-family on crustacean vitellogenesis have also been validated. Advancement in transcriptomic sequencing and comparative genome analysis has led to the discovery of a large number of neuromediators, peptides, and putative peptide receptors having pleiotropic and novel functions in decapod reproduction. Furthermore, differing research strategies have indicated that neurotransmitters and steroid hormones play an integrative role by stimulating neuropeptide secretion, thus demonstrating the complex intertwining of regulatory factors in reproduction. However, the molecular mechanisms by which the combinatorial effect of eyestalk hormones, neuromediators and other factors coordinate to regulate ovarian maturation remain elusive. These multifunctional substances are speculated to control ovarian maturation possibly via the autocrine/paracrine pathway by acting directly on the gonads or by indirectly exerting their stimulatory effects by triggering the release of a putative gonad stimulating factor from the thoracic ganglion. Acting through receptors, they possibly affect levels of cyclic nucleotides (cAMP and cGMP) and Ca²⁺ in target tissues leading to the regulation of vitellogenesis. The "stimulatory paradox" effect of eyestalk ablation on ovarian maturation continues to be exploited in commercial aquaculture operations, and is outweighed by the detrimental physiological effects of this procedure. In this regard, the development of efficient alternatives to eyestalk ablation based on scientific knowledge is a necessity. In this article, we focus principally on the signaling pathways of positive neuromediators and other factors regulating crustacean reproduction, providing an overview of their proposed receptor-mediated stimulatory mechanisms, intracellular signaling, and probable interaction with other hormonal signals. Finally, we provide insight into future research directions on crustacean reproduction as well as potential applications of such research to aquaculture technology development.

Functional analysis of a SoxE gene in the oriental freshwater prawn, Macrobrachium nipponense by molecular cloning, expression pattern analysis, and in situ hybridization (de Haan, 1849)

In this study, a full-length cDNA sequence of SoxE (subgroup E within the Sox family of transcription factors) was cloned from Macrobrachium nipponense and named MnSoxE1. The full-length cDNA of MnSoxE1 is 1748 bp, consisting of a 110 bp 5' UTR, a 105 bp 3' UTR, and a 1533 bp ORF that encodes 510 amino acids. Conserved domains showed that MnSoxE1 has a high similarity to the SoxE gene of Penaeus vannamei. Phylogenetic tree analysis classified that MnSoxE1 with the SoxE gene of other arthropods into one clade. These results suggested that MnSoxE1 belongs to the SoxE subgroup. During embryonic development, MnSoxE1 was mainly expressed in the gastrula stage, implicating its involvement in tissue cell differentiation and formation. In the post-larval stages, the expression of MnSoxE1 continued to increase on days 1-10. The expression level in males was significantly higher than that in females. Males are clearly distinguishable from females on post-larval day 25, showing that MnSoxE1 may play a role in promoting early development and germ cell and gonadal differentiation, especially for males. qPCR analysis showed that MnSoxE1 may also be involved in oogonium proliferation during ovary development. Further in situ hybridization analysis revealed that MnSoxE1 was mainly located in oocytes and spermatocytes, especially in sertoli cells, and implies that it may be involved in the development of oocytes and spermatocytes, as well as the maintenance of testes in mature prawns. These results indicate that MnSoxE1 is involved in gonadal differentiation and development in M. nipponense, especially testis development.

Identification of Peptides and Their GPCRs in the Peppermint Shrimp Lysmata vittata, a Protandric Simultaneous Hermaphrodite Species

Peptide hormones commonly binding with G-protein coupled receptors (GPCRs) achieve their function in reproduction. The peppermint shrimp Lysmata vittata popular in marine ornamental trade and is known to display protandric simultaneous hermaphrodite (PSH). Knowledge on reproductive biology of this commercial species is critical for resources management and aquaculture. This study employed Illumina sequencing and bioinformatics analysis to identify peptides and their candidate GPCRs from male phase (MP) and euhermaphrodite phase (EP) of L. vittata. A total of 61 peptide and 40 peptide GPCR transcripts derive from 44 peptide families and 13 peptide GPCR families were identified, respectively. Among them, insulin-like androgenic gland hormone and crustacean female sex hormone have two unique mature peptides, respectively, and their transcripts showed higher expression levels in MP than EP, which suggest that these sex differentiation hormones might be involved in sexual characters than spermatogenesis or vitellogenesis. Overall, the first study on identification of peptides and their GPCRs in the genus Lysmata extends our knowledge of peptidergic signaling in PSH species, and provides an important basis for development of aquaculture strategies.

Identification of potentially novel functions of DNA polymerase zeta catalytic subunit in oriental river prawn, Macrobrachium nipoponense: cloning, qPCR, in situ hybridization and RNAi analysis

The goal of this study was to analyze the functions of DNA polymerase zeta catalytic subunit (Rev3) in the oriental river prawn Macrobrachium nipponense (Mn-Rev3) with a focus on its potential roles in sex differentiation and development. The full length of Mn-Rev3 cDNA sequence was 6832 base pairs (bp) with an open reading frame of 6102 bp encoding 2033 amino acids. Mn-Rev3 showed the closest evolutionary relationship with Penaeus vannamei. The highest expression level of Mn-Rev3 occurred in the hepatopancreas and strong signals were observed in hepatopancreas cells, suggesting that Mn-Rev3 played a role in the immune system. Expression levels of Mn-Rev3 also were relatively high in the androgenic gland and testis, suggesting its potential roles in male sexual differentiation and development. During development, expression of Mn-Rev3 was highest on larval day 15 and relatively high from post-larval day 1 (PL1) to PL15, indicating that it played essential roles in promoting metamorphosis and gonad differentiation and development in M. nipponense. Strong Mn-Rev3 signals were detected in spermatids, spermatocytes, and sperm in the testes, and Mn-Rev3 expression was higher in the testes during the reproductive season than in the non-reproductive season. This result indicated that Rev3 promoted whole testis development, and especially sperm development, in M. nipponense. The expression level of Mn-Rev3 was high from ovary V to ovary II stages, indicating that Rev3 may be involved in yolk deposition. The expression level of Mn-insulin-like androgenic gland hormone (Mn-IAG) and the content of testosterone showed the same expression pattern as that of Mn-Rev3 after injection of double-stranded RNA of Mn-Rev3, which indicated that Rev3 had positive effects on male sexual differentiation and development in M. nipponense. The results of this study advance our understanding of male sexual development in M. nipponense and provide the basis for further studies of male sexual differentiation and development in crustaceans.

Integrated analysis of differentially expressed microRNAs and mRNAs to screen miRNAs and genes related to reproduction in Macrobrachium nipponense

Female Macrobrachium nipponense has the characteristic of short sexual maturity during the breeding season, which can increase breeding risk and lead to prevalent female individual miniaturization. In this study, we characterized micro (mi)RNA-seq data of the eyestalk (E) and cerebral ganglia (B) of female M. nipponense during breeding and non-breeding seasons. A total of 393 and 189 differentially expressed miRNAs (DE miRNAs) were identified in BSE vs. NBSE and BSB vs. NBSB, respectively. The most abundant up- and down-regulated DE miRNAs were miR-124, miR-14, and miR-7. Enrichment analysis showed that DE miRNA target genes were mainly involved in 'metabolic process' and 'binding', and were associated with 'neurohormonal regulation' and 'photoreceptor activity' signaling pathways. Integrated analysis of miRNA-mRNA expression showed that the most abundant DE miRNAs were miR-14 and miR-278 in BSE vs. NBSE and BSB vs. NBSB, respectively. Four pairs of DE miRNAs and their corresponding target annotated genes were selected from the DE miRNA-mRNA interaction network (bmo-miR-316-5p/opsin protein, ame-miR-125/skeletal muscle actin 8, dmo-miR-278/sugar transporter, and tca-miR-3885-5p/5-HT1 receptor). Gene expression analysis of these four pairs in different ovary development stages showed their potential regulatory roles in ovary maturation.

Potential Functions of Gem-Associated Protein 2-Like Isoform X1 in the Oriental River Prawn Macrobrachium nipponense: Cloning, qPCR, In Situ Hybridization, and RNAi Analysis

Gem-associated protein 2-like isoform X1 (GEM) was previously predicted to be involved in the sexual development of male Macrobrachium nipponense. In this study, we analyze the GEM functions in depth using quantitative polymerase chain reaction (qPCR), in situ hybridization, and RNA interference (RNAi). The full-length Mn-GEM cDNA sequence was 1018 base pairs (bp) long with an open reading frame of 777 bp encoding 258 amino acids. qPCR analysis of Mn-GEM in different tissues and developmental stages showed that Mn-GEM was highly expressed in the gonad and from post-larval developmental stage day 5 (PL5) to PL15, which indicated that GEM has potential roles in gonad differentiation and development in M. nipponense. In situ hybridization and qPCR analysis of various stages of the reproductive cycle of the testis and ovary indicated that GEM may promote spermatid development and gametogenesis in M. nipponense. After injecting with double-stranded RNA (dsRNA) of Mn-GEM, mRNA expression of Mn-insulin-like androgenic gland hormone (Mn-IAG) and the content of testosterone increased with the decrease of Mn-GEM expression, indicating that GEM has negative effects on the male sexual differentiation and development in M. nipponense. Results of this study highlight the functions of GEM in M. nipponense, which can be applied to future studies of male sexual development in M. nipponense and other crustacean species.

Molecular cloning, expression pattern analysis, and in situ hybridization of a Transformer-2 gene in the oriental freshwater prawn, Macrobrachium nipponense (de Haan, 1849)

In this study, we isolated a full-length cDNA sequence from Macrobrachium nipponense and investigated its gene function. We named the gene Mntra-2a because of high similarities and close evolutionary divergence with arthropod tra-2. The full-length cDNA of Mntra-2a was 1293 bp, consisting of a 212 bp 5' UTR, a 268 bp 3' UTR, and an ORF of 813 bp encoding 270 amino acids. It contained an RNA recognition motif and a linker region. Real-time PCR analysis showed that Mntra-2a was highly expressed in the gonads of both males and females. Further in situ hybridization analysis showed that Mntra-2a was mainly located in oocytes and spermatocytes. During embryogenesis, Mntra-2a expression was higher in the cleavage and nauplius stages. During the ovarian reproductive cycle, Mntra-2a expression reached a peak at OvaryV and decreased to the lowest level at OvaryIV. These results indicated that Mntra-2a probably played important roles in embryonic development and early gonad development in M. nipponense. Our results provide basic information for further functional studies of tra-2 in M. nipponense.

Identification and Characterization of the DMRT11E Gene in the Oriental River Prawn Macrobrachium nipponense

The doublesex and mab-3 related transcription factor (DMRT) gene family involvement in sex development is widely conserved from invertebrates to humans. In this study, we identified a DM (Doublesex/Mab-3)-domain gene in Macrobrachium nipponense, which we named MniDMRT11E because it has many similarities to and phylogenetically close relationships with the arthropod DMRT11E. Amino acid alignments and structural prediction uncovered conservation and putative active sites of the DM domain. Real-time PCR analysis showed that the MniDMRT11E was highly expressed in the ovary and testis in both males and females. Cellular localization analysis showed that DMRT11E was mainly located in the oocytes of the ovary and the spermatocyte of the testis. During embryogenesis, the expression level of MniDMRT11E was higher at the cleavage stage than at other stages. During the different stages of ovarian development, MniDMRT11E expression gradually increased from OI to OIII and decreased to the lowest level at the end of OIV. The results indicated that MniDMRT11E probably played important roles in embryonic development and sex maturity in M. nipponense. MniDMRT11E dsRNA injection also significantly reduced vitellogenin (VG) expression and significantly increased insulin-like androgenic gland factor (IAG) expression, indicating a close relationship in gonad development.

Long noncoding RNA and mRNA expression profiles following igf3 knockdown in common carp, Cyprinus carpio

As a novel IGF system member, igf3 plays an important role in gonadal development of teleost fish. Although studies have reported the unusual expression of igf3 in fish gonad, whether the igf3 affects the expression of long noncoding RNAs (lncRNAs) in gonad remains unknown. In this study, an igf3 knockdown common carp (Cyprinus carpio) model was established by RNA interference. Then RNA sequencing of C. carpio gonad after igf3 knockdown was performed. A total of 327,169,410 and 306,305,018 clean reads were identified from control and igf3-dsRNA interference group, respectively. After a stringent filtering, RNA-seq yielded 14199 lncRNA and 106932 mRNA transcripts with 124 and 353 differentially expressed lncRNAs and mRNAs. Our dataset provides an extensive resource for understanding the potential regulatory molecular mechanism of igf3 in early stage of gonadal development in C. carpio.

Molecular Cloning and Expression of MnGST-1 and MnGST-2 from Oriental River Prawn, Macrobrachium nipponense, in Response to Hypoxia and Reoxygenation

The glutathione-S-transferase (GST) superfamily includes seven classes, and different classes have different functions. GST superfamily members function in various processes including detoxification of xenobiotics, protection against oxidative damage, and intracellular transport of hormones, endogenous metabolites, and exogenous chemicals. Herein, to elucidate the tissue-specific expression pattern of GSTs in response to hypoxia stress, which induces cell death, we investigated the expression of GSTs in response to hypoxia and reoxygenation in oriental river prawn, Macrobrachium nipponense. Full-length cDNAs of two δ class GSTs were cloned from the hepatopancreas, and named MnGST-1 and MnGST-2 based on the established GST nomenclature system. Expression profiles of both GSTs in various tissues were different under acute and chronic experimental hypoxia stress conditions, suggesting that both respond strongly to hypoxia-induced oxidative stress. However, the intensity of responses to hypoxia and reoxygenation were different in different tissues. During acute hypoxia stress, MnGST-1 responds earlier than MnGST-2 in the hepatopancreas and gill, but more slowly in muscle. By contrast, during chronic hypoxia stress, MnGST-2 plays a more important role in the hepatopancreas and gill than MnGST-1.

Myosin Va from Eriocheir sinensis: cDNA cloning, expression and involvement in growth and development

Myosin Va, a member of the myosin superfamily, has been widely identified associated with processes of cellular motility, which include neurotransmitter release and synaptic plasticity during neurodevelopment. However, the function of myosin Va in the growth and development of crustaceans has not yet been reported. In this study, a full-length cDNA of myosin Va (named as EsMyoVa) was cloned from the Chinese mitten crab, Eriocheir sinensis, and the expression patterns were detected in different tissues and larval developmental stages. The full-length cDNA of EsMyoVa was 6037 bp in length. Real time quantitative reverse transcription PCR (qRT-PCR) analysis showed that EsMyoVa transcript has a wide tissue distribution pattern and is expressed in zoeae, megalopa, juvenile crab stages and adults. In order to further study the function of this gene, we used RNAi technology in the muscle, hepatopancreas, gill, and gonad. After double-stranded RNA (dsRNA) injection, the expression level of EsMyoVa was significantly decreased in all tissues in both sexes and the gene knockdown effects of dsRNA persisted for at least 6 days. Subsequently, the role of EsMyoVa was revealed by silencing the transcript through one month injections of Myosin Va dsRNA. Crabs with reduced levels of EsMyoVa transcripts displayed a dramatic slowing in growth rate and considerably higher mortality compared to control groups, which indicated that this gene had important role of regulating growth and development.

Validation and Evaluation of Reference Genes for Quantitative Real-Time PCR in Macrobrachium Nipponense

Quantitative real-time PCR (qPCR) is widely used in molecular biology, although the accuracy of the quantitative results is determined by the stability of the reference genes used. Recent studies have investigated suitable reference genes for some crustaceans under various conditions, but studies in Macrobrachium nipponense are currently lacking. In this study, we selected the following seven genes from among 35 commonly used housekeeping genes as candidate qPCR reference genes for temporal and spatial expression: EIF (eukaryotic translation initiation factor 5A), 18S (18S ribosomal RNA), EF-1α (elongation factor-1α), GAPDH (glyceraldehyde-3-phosphate dehydrogenase), TUB (α-tubulin), β-act (β-actin), and RPL18 (Ribosomal protein L18). The stability of each reference gene was evaluated by GeNorm, NormFinder, BestKeeper, and comparative ∆C _t methods, and was comprehensively ranked using RefFinder. RPL18 was shown to be the most suitable reference gene for adult M. nipponense tissues, while EIF was the most stable in different ovarian and embryo stages and in white spot syndrome virus infection, and β-act was the most stable reference gene under hypoxia stress. The reliability of the rankings was confirmed by RNA interference experiments. To the best of our knowledge, this represents the first systematic analysis of reference genes for qPCR experiments in M. nipponense, and the results will provide invaluable information for future research in closely related crustaceans.

Molecular Cloning and Expression Analysis of Lactate Dehydrogenase from the Oriental River Prawn Macrobrachium nipponense in Response to Hypoxia

Metabolic adaption to hypoxic stress in crustaceans implies a shift from aerobic to anaerobic metabolism. Lactate dehydrogenase (LDH) is a key enzyme in glycolysis in prawns. However, very little is known about the role of LDH in hypoxia inducible factor (HIF) pathways of prawns. In this study, full-length cDNA of LDH (MnLDH) was obtained from the oriental river prawn Macrobrachium nipponense, and was characterized. The full-length cDNA is 2267-bp with an open reading frame of 999 bp coding for a protein of 333 amino acids with conserved domains important for function and regulation. Phylogenetic analysis showed that MnLDH is close to LDHs from other invertebrates. Quantitative real-time PCR revealed that MnLDH is expressed in various tissues with the highest expression level in muscle. MnLDH mRNA transcript and protein abundance in muscle, but not in hepatopancreas, were induced by hypoxia. Silencing of hypoxia-inducible factor 1 (HIF-1) α or HIF-1β subunits blocked the hypoxia-dependent increase of LDH expression and enzyme activity in muscle. A series of MnLDH promoter sequences, especially the full-length promoter, generated an increase in luciferase expression relative to promoterless vector; furthermore, the expression of luciferase was induced by hypoxia. These results demonstrate that MnLDH is probably involved a HIF-1-dependent pathway during hypoxia in the highly active metabolism of muscle.