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

Transcriptomic Analysis of the Hepatopancreas in the Sex-Related Size Differences of Macrobrachium nipponense

Macrobrachium nipponense, a commercially popular crustacean species within the Chinese context, is recognized for its exceptional nutritional composition and palatability. There are significant differences in growth between male and female M. nipponense. Herein, transcriptomics was used to determine the hepatopancreas transcriptome differences between sex-related size differences in M. nipponense. We identified 974 differentially expressed genes (DEGs) between the SHE (female) and BHE (male) groups, which were validated using RT-qPCR. The genes encoding matrix metalloproteinase-9 (MM9), Ribosome-binding protein 1 (RBP1), Aly/REF export factor 2, and hematological and neurological expressed 1 (HN1) may play a role in modulating the sex-related size differences observed in M. nipponense. Clusters of orthologous groups and gene ontology functional analysis demonstrated that the DEGs for sex-related size in M.nipponense were associated with various biological functions. The Kyoto Encyclopedia of Genes and Genomes pathways analysis demonstrated that upregulated DEGs were mainly enriched in lysine biosynthesis, tryptophan metabolism, and lysine degradation pathways, whereas the downregulated DEGs were mainly enriched in ascorbate and aldarate metabolism, retinol metabolism, and drug metabolism-cytochrome P450 pathways. The results indicated the molecular mechanism underlying the sex-related size differences and identified key genes. This data will be invaluable to support explanations of individual differences between male and female prawns.

A pattern recognition receptor interleukin-1 receptor is involved in reproductive immunity in Macrobrachium nipponense ovary

The family of TIR domain-containing receptors includes numerous proteins involved in innate immunity. In this study, a member of this family was characterized from the ovary of the oriental river prawn Macrobrachium nipponense and identified as interleukin-1 receptor (MnIL-1R). Meanwhile, to elucidate the conservation of IL-1R, its orthologous were identified in several crustacean species as well. In addition, the expression pattern of MnIL-1R in various adult tissues and post different pathogen-associated molecular patterns (PAMPs) challenge in ovary was analyzed with qRT-PCR technology. Finally, the roles of MnIL-1R in the ovary were analyzed by RNAi technology. The main results are as follows: (1) MnIL-1R comprises a 1785 bp ORF encoding 594 amino acids and is structurally composed of five domains: a signal peptide, two immunoglobulin (IG) domains, a transmembrane region, and a TIR-2 domain; (2) the TIR domain showed a high conservation among analyzed crustacean species; (3) MnIL-1R is widely detected in all tested tissues including ovary; (4) MnIL-1R showed a positive response to challenges with LPS, PGN, and polyI:C in the ovary; (5) its IG domain showed strong binding ability to LPS and PGN, confirming its role as a pattern recognition receptor; (6) the expression patterns of several members of the Toll signaling pathway (Myd88, TRAF-6, Dorsal, and Relish) was similar to that of MnIL-1R after challenges with LPS, PGN, and polyI:C in the ovary; (7) the silencing of MnIL-1R resulted in down-regulation of theses gene' (Myd88, TRAF-6, Dorsal, and Relish) expression level in the ovary. These results suggest that MnIL-1R can activate the Toll signaling pathway in the ovary by directly recognizing LPS and PGN through its IG domain, thereby contributing to the immune response in the ovary of M. nipponense.

Microcystin-LR-induced autophagy via miR-282-5p/PIK3R1 pathway in Eriocheir sinensis hepatopancreas

With the intensifying climate warming, blue-green algae blooms have become more frequent and severe, releasing environmental hazards such as microcystin that pose potential threats to human and animal health. Autophagy has been shown to play a crucial role in regulating immune responses induced by environmental hazards, enabling cells to adapt to stress and protect against damage. Although microcystin-LR (MC-LR) has been identified to affect autophagy in mammalian, its impact on aquatic animals has been poorly studied. To investigate the toxicological effects of MC-LR in aquatic ecosystems, we constructed a microRNA profile of acute MC-LR stress in the hepatopancreas of the Chinese mitten crab. Interestingly, we found the MC-LR exposure activated autophagy in the hepatopancreas based on the following evidence. Specifically, mRNA expression level of ATG7, Beclin1 and Gabarap was significantly up-regulated, autophagy regulatory pathways were significantly enriched, and numerous autolysosomes and autophagosomes were observed. Additionally, we found that miR-282-5p and its target gene PIK3R1 played important regulatory roles in autophagy by in vivo and in vitro experiments. Overexpression of miR-282-5p mimicked MC-LR-induced autophagy by inhibiting PIK3R1 expression, while miR-282-5p silencing inhibited autophagy by promoting PIK3R1 expression. Altogether, our findings suggest that MC-LR increases miR-282-5p, which then targets inhibition of PIK3R1 to stimulate autophagy. This study focused on the stress response regulatory mechanisms of juvenile crabs to toxic pollutants in water, offering a potential target for alleviating the toxicity of MC-LR. These findings lay a foundation for reducing the toxicity of MC-LR and environmental hazards in organisms.

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.

PacBio single-molecule long-read sequencing provides new insights into the complexity of full-length transcripts in oriental river prawn, macrobrachium nipponense

BACKGROUND

Oriental river prawn (Macrobrachium nipponense) is one of the most dominant species in shrimp farming in China, which is a rich source of protein and contributes to a significant impact on the quality of human life. Thus, more complete and accurate annotation of gene models are important for the breeding research of oriental river prawn.

RESULTS

A full-length transcriptome of oriental river prawn muscle was obtained using the PacBio Sequel platform. Then, 37.99 Gb of subreads were sequenced, including 584,498 circular consensus sequences, among which 512,216 were full length non-chimeric sequences. After Illumina-based correction of long PacBio reads, 6,599 error-corrected isoforms were identified. Transcriptome structural analysis revealed 2,263 and 2,555 alternative splicing (AS) events and alternative polyadenylation (APA) sites, respectively. In total, 620 novel genes (NGs), 197 putative transcription factors (TFs), and 291 novel long non-coding RNAs (lncRNAs) were identified.

CONCLUSIONS

In summary, this study offers novel insights into the transcriptome complexity and diversity of this prawn species, and provides valuable information for understanding the genomic structure and improving the draft genome annotation of oriental river prawn.

Genome-Wide Association Study of Growth and Sex Traits Provides Insight into Heritable Mechanisms Underlying Growth Development of Macrobrachium nipponense (Oriental River Prawn)

Male hybrid oriental river prawns grow significantly faster than hybrid females. In this study, the growth and sex traits of 181 individuals of Macrobrachium nipponense were recorded, and each individual genotype was evaluated using the 2b-RAD sequencing method. The genetic parameters for growth and sex traits were estimated. A genome-wide association analysis (GWAS) of these traits was performed. In total, 18 growth-related SNPs were detected from 12 chromosomes using a mixed linear model. The most significant loci of weight are located on the position of the SNP (102638935, chromosome 13), which can explain 11.87% of the phenotypic variation. A total of 11 significant SNPs were detected on four chromosomes associated with sex trait (three on chromosome 4, one on chromosome 7 and seven on chromosome 17). The heritability of this trait is 0.8998 and belongs to the range of ultra-high heritability. Genetic correlations were prevalent among the 11 traits examined, the genetic coefficient between sex and body weight reached a significant level of −0.23. This study is the first GWAS for sex of binary and growth traits in oriental river prawn. Our results provide a set of markers for the genetic selection of growth traits and help us to further understand the genetic mechanisms of growth in Macrobrachium 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.

Eyestalk neuropeptide identification in the female red deep-sea crab, Chaceon quinquedens

Eyestalk-derived neuropeptides, primarily the crustacean hyperglycemic hormone (CHH) neuropeptide family, regulate vitellogenesis in decapod crustaceans. The red deep-sea crab, Chaceon quinquedens, a cold-water species inhabiting depths between 200 and 1800 m, has supported a small fishery, mainly harvesting adult males in the eastern US for over 40 years. This study aimed to understand the role of eyestalk-neuropeptides in vitellogenesis in C. quinquedens with an extended intermolt stage. Chromatography shows two CHH and one MIH peak in the sinus gland, with a CHH2 peak area four times larger than CHH1. The cDNA sequence of MIH and CHH of C. quinquedens is isolated from the eyestalk ganglia, and the qPCR assay shows MIH is significantly higher only at ovarian stages 3 than 4 and 5. However, MIH transcript and its neuropeptides do differ between stages 1 and 3. While CHH transcripts remain constant, its neuropeptide levels are higher at stages 3 than 1. Additionally, transcriptomic analysis of the de novo eyestalk ganglia assembly at ovarian stages 1 and 3 found 28 eyestalk neuropeptides. A GIH/VIH or GSH/VSH belonging to the CHH family is absent in the transcriptome. Transcripts per million (TPM) values of ten neuropeptides increase by 1.3 to 2.0-fold at stage 3 compared to stage 1: twofold for Bursicon α, followed by CHH, AKH/corazonin-like, Pyrokinin, CCAP, Glycoprotein B, PDH1, and IDLSRF-like peptide, and 1.3-fold of allatostatin A and short NP-F. WXXXRamide, the only downregulated neuropeptide, decreases TPM by ∼ 2-fold at stage 3, compared to stage 1. Interestingly, neuroparsin with the highest TPM values remains the same in stages 1 and 3. The mandibular organ-inhibiting hormone is not found in de novo assembly. We report that CHH, MIH, and eight other neuropeptides may play a role in vitellogenesis in this species.

Transcriptome analysis of flathead grey mullet (Mugil cephalus) ovarian development induced by recombinant gonadotropin hormones

Background: Treatment with recombinant gonadotropin hormones (rGths), follicle-stimulating hormone (rFsh) and luteinizing hormone (rLh), was shown to induce and complete vitellogenesis to finally obtain viable eggs and larvae in the flathead grey mullet (Mugil cephalus), a teleost arrested at early stages of gametogenesis in intensive captivity conditions. This study aimed to investigate the transcriptomic changes that occur in the ovary of females during the rGths-induced vitellogenesis. Methods: Ovarian samples were collected through biopsies from the same five females at four stages of ovarian development. RNASeq libraries were constructed for all stages studied, sequenced on an Illumina HiSeq4000, and a de novo transcriptome was constructed. Differentially expressed genes (DEGs) were identified between stages and the functional properties of DEGs were characterized by comparison with the gene ontology and Kyoto Encyclopedia. An enrichment analysis of molecular pathways was performed. Results: The de novo transcriptome comprised 287,089 transcripts after filtering. As vitellogenesis progressed, more genes were significantly upregulated than downregulated. The rFsh application induced ovarian development from previtellogenesis to early-to-mid-vitellogenesis with associated pathways enriched from upregulated DEGs related to ovarian steroidogenesis and reproductive development, cholesterol metabolism, ovarian growth and differentiation, lipid accumulation, and cell-to-cell adhesion pathways. The application of rFsh and rLh at early-to-mid-vitellogenesis induced the growth of oocytes to late-vitellogenesis and, with it, the enrichment of pathways from upregulated DEGs related to the production of energy, such as the lysosomes activity. The application of rLh at late-vitellogenesis induced the completion of vitellogenesis with the enrichment of pathways linked with the switch from vitellogenesis to oocyte maturation. Conclusion: The DEGs and enriched molecular pathways described during the induced vitellogenesis of flathead grey mullet with rGths were typical of natural oogenesis reported for other fish species. Present results add new knowledge to the rGths action to further raise the possibility of using rGths in species that present similar reproductive disorders in aquaculture, the aquarium industry as well as the conservation of endangered species.

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.

RNA-seq Provides Novel Insights into Response to Acute Salinity Stress in Oriental River Prawn Macrobrachium nipponense

The oriental river prawn Macrobrachium nipponense is an important aquaculture species in China, Vietnam, and Japan. This species could survive in the salinity ranging from 7 to 20 ppt and accelerate growth in the salinity of 7 ppt. To identify the genes and pathways in response to acute high salinity stress, M. nipponense was exposed to the acute high salinity of 25 ppt. Total RNA from hepatopancreas, gills, and muscle tissues was isolated and then sequenced using high-throughput sequencing method. Differentially expressed genes (DGEs) were identified, and a total of 632, 836, and 1246 DEGs with a cutoff of significant twofold change were differentially expressed in the hepatopancreas, gills, and muscle tissues, respectively. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genome pathway enrichment analyses were conducted. These DEGs were involved in the GO terms of cellular process, metabolic process, membrane, organelle, binding, and catalytic activity. The DEGs of hepatopancreas and gill tissues were mainly enriched in PPAR signaling pathway, longevity regulating pathway, protein digestion and absorption, and the DEGs of muscle tissue in arginine biosynthesis, adrenergic signaling in cardiomyocytes, cardiac muscle contraction, and cGMP-PKG signaling pathway. Real-time PCR conducted with fifteen selected DEGs indicated high reliability of digital analysis using RNA-Seq. The results indicated that the M. nipponense may regulate essential mechanisms such as metabolism, oxidative stress, and ion exchange to adapt the alternation of environment, when exposed to acute high salinity stress. This work reveals the numbers of genes modified by salinity stress and some important pathways, which could provide a comprehensive insight into the molecular responses to high salinity stress in M. nipponense and further boost the understanding of the potential molecular mechanisms of adaptation to salinity stress for euryhaline crustaceans.

Microcystin-LR induces apoptosis in Juvenile Eriocheir sinensis via the mitochondrial pathway

Microcystin-LR (MC-LR), the toxic substance of cyanobacteria secondary metabolism, widely exists in water environments and poses great risks to living organisms. Some toxicological assessments of MC-LR have performed at physiological and biochemical levels. However, plenty of blanks about the potential mechanism in aquatic crustacean remains. In this study, we firstly assessed the exposure toxicity of MC-LR to juvenile E. sinensis and clarified that the 96 h LD50 of MC-LR was 73.23 μg/kg. Then, hepatopancreas transcriptome profiles of MC-LR stressed crabs were constructed at 6 h post-injection and 37 differential expressed genes (DEGs) were identified. These DEGs were enriched in cytoskeleton, peroxisome and apoptosis pathways. To further reveal the toxicity of MC-LR, oxidative stress parameters (SOD, CAT, GSH-px and MDA), apoptosis genes (caspase 3, bcl-2 and bax) and apoptotic cells were detected. Significant accumulated MDA and rise-fall enzyme activities verified the oxidative stress caused by MC-LR. It is noteworthy that quantitative real-time PCR and TUNEL assay indicated that MC-LR stress-induced apoptosis via the mitochondrial pathway. Interestingly, activator protein-1 may play a crucial role in mediating the hepatotoxicity of MC-LR by regulating apoptosis and oxidative stress. Taken together, our study investigated the toxic effects and the potential molecular mechanisms of MC-LR on juvenile E. sinensis. It provided useful data for exploring the toxicity of MC-LR to aquatic crustaceans at molecular levels.

A Novel Legumain-Like Protease in Macrobrachium nipponense: Identification, Characterization, and Function Analysis in Ovary Maturation

Legumain, also called aspartic endopeptidase (AEP), is a member of the cysteine protease family and is involved in various physiological processes. In this study, we analyzed the characteristics of a novel legumain-like (named Mn-Lel) in the female oriental river prawn, Macrobrachium nipponense, which is involved in ovary maturation. The Mn-Lel is 1,454 bp in length, including a 1,290-bp open reading frame that encodes 430 amino acids. qPCR analysis indicated that Mn-Lel is specifically highly expressed in the hepatopancreas and ovaries of female prawns. It is rarely expressed in embryogenesis, weakly expressed in early larval development stages, and then significantly increased after metamorphosis, which indicated that Mn-Lel is not a maternal gene and mainly plays a role in adults. During the different ovarian stages, Mn-Lel expression in the hepatopancreas had no obvious rules, while its expression in the ovaries had a significant peak in stage III. In situ hybridization studies revealed that Mn-Lel is localized in the oocyte of the ovary. Changes in the gonadosomatic index confirmed the inhibitory effects of Mn-Lel dsRNA on ovary maturation. These results suggest that Mn-Lel has a key role in promoting ovary maturation.