Feeding rhythm of the zoea larvae of Scylla paramamosain: The dynamic feeding rhythm is not completely synchronized with photoperiod.

The feeding rhythm is one of the key factors determining the success of artificial breeding of S. paramamosain. To understand the feeding rhythm of the different zoea larva developmental stages of S. paramamosain, the feeding rate, digestive enzyme activity, and expression of metabolism-related genes were investigated in the present study. The results showed that the S. paramamosain feeding rate has strong diurnal feeding rhythm, being significantly higher at 10:00-14:00 from stages ZI to ZIV. While the feeding rate peaked at 14:00 on Days 10 and 11, the peak shifted to 18:00 on Day 12. The activity of digestive enzymes amylase, pepsin and lipase decreased at night but increased in the daytime, showing a single-phase rhythm similar to that of the feeding rate, suggesting that the digestive enzyme activity was closely associated with the feeding rate during the larval development. Compared to pepsin and lipase, the activity of amylase was the most consistent with feeding rate. In particular, amylase activity peaked at 18:00 on Day 12. Due to its synchronicity with feeding activity, the activity of amylase could provide a potential reference for determining the best feeding time during zoea stages in S. paramamosain breeding. Moreover, the relative mRNA expression of metabolism-related genes SpCHH and SpFAS at most tested points was lower from 10:00 to 14:00, but higher at 18:00 to 6:00 of the next day. On the other hand, the expression patterns of SpHSL and SpTryp were converse to those of SpCHH and SpFAS. Our findings revealed that the S. paramamosain zoea has an obvious feeding rhythm, and the most suitable feeding time was 10:00-18:00 depending on different stages. The feeding rhythm is a critical aspect in aquaculture, influencing a series of physiological functions in aquatic animals. This study provides insights into the feeding rhythm during the zoea development of S. paramamosain, making a significant contribution to optimizing feeding strategy, improving aquafeed utilization, and reducing the impact of residual feed on water environment.

Temporal and Spatial Signatures of Scylla paramamosain Transcriptome Reveal Mechanistic Insights into Endogenous Ovarian Maturation under Risk of Starvation.

Variability in food availability leads to condition-dependent investments in reproduction. This study is aimed at understanding the metabolic response and regulatory mechanism of female Scylla paramamosain in response to starvation in a temporal- and tissue-specific manner. The mud crabs were starved for 7 (control), 14, 28, and 40 days for histological and biochemical analysis in the hepatopancreas, ovary, and serum, as well as for RNA sequencing on the hepatopancreas and ovary. We further highlighted candidate gene modules highly linked to physiological traits. Collectively, our observations suggested that starvation triggered endogenous ovarian maturation at the expense of hepatopancreas mass, with both metabolic adjustments to optimize energy and fatty acid supply from hepatopancreas to ovary in the early phase, followed by the activation of autophagy-related pathways in both organs over prolonged starvation. These specific adaptive responses might be considered efficient strategies to stimulate ovarian maturation of Scylla paramamosain under fasting stress, which improves the nutritional value of female mud crabs and other economically important crustaceans.

Enzymatic Stress Responses of Coreius guichenoti to Microplastics with Different Particle Sizes.

The wild population resources of Coreius guichenoti have sharply declined in recent decades, and any negative factors may have a significant impact on their survival. In this study, the enzymatic stress responses of C. guichenoti to 25 and 48 µm polyethylene fragments were explored for the first time. This was achieved by evaluating the changes in physiological and biochemical indicators of the species in response to the environmental stimuli of microplastics. In this study, we observed an early stress response in the external tissues of C. guichenoti following exposure to microplastics. The TP content in skin and muscle and the MDA content in skin, gill and muscle initially showed a significant increase. The skin, gill, and muscle exhibited greater stress responses to M5 particles, whereas M3 particles caused a greater response in the intestine and especially the liver. After the removal of microplastic exposure, the stress state of the C. guichenoti would be alleviated in a short period, but it could not fully recover to the pre-exposure level. In summary, microplastics pose a significant threat to C. guichenoti. While their negative effects can be alleviated by the removal of microplastics exposure, full recovery does not occur in a short period. Continuous monitoring of microplastics in natural waters and targeted aquatic ecological restoration are essential to ensure the normal growth and reproduction of the wild population of C. guichenoti.

Molecular characterization and transcriptional response to TiO2-GO nanomaterial exposure of two molt-related genes in the juvenile prawn, Macrobrachium rosenbergii.

In recent years, with the widespread use of TiO2-GO nanocomposite in industry, especially in the remediation of water environments, its toxic effects on aquatic organisms have received increasing attention. As molting is extremely important for crustaceans in their growth, in this study, we cloned the full-length cDNA sequences of two key genes related to molting, nuclear hormone receptor E75 (E75) and nuclear hormone receptor HR3 (HR3), in Macrobrachium rosenbergii, examined the gene expression profile, and investigated their toxicological effects on crustacean molting through nanomaterial exposure. The amino acid sequences for E75 and HR3 were respectively determined to encode 1138 and 363 acid residues. Sequence analysis showed that both E75 and HR3 contain a HOLI domain, with the E75 of M. rosenbergii being more closely related to the E75 of Palaemon carinicauda. These two genes were expressed at the highest levels in muscle, followed by hepatopancreas. The results showed that the expressions of E75 and HR3 in hepatopancreas and muscle tissues were significantly decreased after exposure to 0.1 mg/L of TiO2-GO composite nanoparticles (P < 0.05). This study will serve as a foundation for subsequent research into the evaluation of nanomaterial toxicity on crustacean species.

Impact of Titanium Dioxide-Graphene Oxide (TiO2-GO) Composite Nanoparticle on the Juveniles of the Giant River Prawn, Macrobrachium rosenbergii: Physio-Biochemistry and Transcriptional Response.

Nanomaterials are used in many fields, resulting in inevitably releasing into the aquatic environment. The presence of nanomaterials, including TiO2-GO in the aquatic environment, can be toxic to aquatic organisms. However, few studies have focused on the effects of TiO2-GO composite nanoparticle on crustaceans. In the present study, the giant river prawn Macrobrachium rosenbergii juveniles were exposed to two concentrations of TiO2-GO composite nanoparticle (0.1 and 0.5 mg/L). The effects of TiO2-GO composite exposure on activities of digestive and antioxidant-related enzymes and expressions of growth and immune-related genes at the transcriptome were studied. The results showed that the survival rate and growth performance were not negatively affected by TiO2-GO composite at the two exposure levels. Nevertheless, exposure to TiO2-GO composite causes an effect on the activities of digestive and antioxidant enzymes in the juvenile prawns. The enzyme activities of CAT, SOD, GSH-Px, AMS, TPS, and LPS in the 0.1 mg/L TiO2-GO composite experimental group were markedly reduced than those in the control group. Additionally, the expression level of genes involved in growth and immunity was significantly affected by TiO2-GO composite. After exposure to the 0.1 mg/L TiO2-GO composite, the mRNA expression level of MSTN was significantly increased, but the level of EcR, Raptor, and CaBP was significantly decreased. However, the mRNA levels of the CTL, TLR, JAK, and STAT were significantly increased after exposure to the 0.5 mg/L concentration of TiO2-GO composite. Furthermore, to understand the molecular mechanism of M. rosenbergii under TiO2-GO composite exposure, RNA-Seq was employed to analyze the changes of the muscle and hepatopancreas transcriptome. Compared with the control group, we identified 5166 and 4784 differentially expressed genes (DEGs) in the muscle and hepatopancreas, respectively (p < 0.05). Based on gene ontology and KEGG analysis, significant differences were observed in the DEGs involved in activity and binding, metabolism, immune response, and environmental information processing. These results showed that exposure to TiO2-GO composite nanoparticle led to the changes of enzyme activity and gene expression, suggesting that TiO2-GO composite existing in aquatic environments would disrupt the physiology of M. rosenbergii. This study will serve as a foundation for subsequent research into the evaluation of nanomaterial toxicity on crustacean species.

Effects of titanium dioxide (TiO2)/activated carbon (AC) nanoparticle on the growth and immunity of the giant freshwater prawn, Macrobrachium rosenbergii: potential toxicological risks to the aquatic crustaceans.

Due to their unique physicochemical characteristics, nanomaterials exhibit many excellent properties and functions, leading to their applications in numerous fields. The large-scale production and widespread application of nanomaterials have inevitably resulted in their release into the environment, especially the water environment. Several studies have confirmed that exposure to nanomaterials can be toxic to aquatic organisms. However, few studies have focused on the effects of nanomaterial exposure on growth and immunity in crustaceans. In the present study, juvenile Macrobrachium rosenbergii were exposed to different concentrations of titanium dioxide (TiO2)/activated carbon (AC) composite nanomaterial (0.1 and 0.5 mg/L) for 45 days. The effects of nanoparticle exposure on digestion and antioxidant-related enzyme activities, as well as the expression of growth and immunity-related genes and signaling pathway, were evaluated. Our results show that in response to low concentration of TiO2/AC nanoparticle (0.1 mg/L), most of the enzyme activities related to digestion and antioxidation (TPS, LPS, AMS, SOD, and CAT) were diminished. On the contrary, the GSH-Px activity increased under the 0.1 mg/L group of TiO2/AC nanoparticle concentration. Additionally, the level of digestive and antioxidant enzyme activities we detected was increased when exposed to 0.5 mg/L TiO2/AC nanoparticle. By comparison to the expression level of growth-related genes in the control group, MSTN, CaBP, E75, Raptor, EcR, and EGF were significantly inhibited at 0.1 and 0.5 mg/L concentrations of TiO2/AC nanoparticle, whereas the expression level of genes (TLR, JAK, STAT, PPAF, ACP, and AKP) related to immunity was increased when exposed to different concentrations of TiO2/AC nanoparticle. Compared with the control group (0 mg/L concentration), 5166 DEGs were identified in the TiO2/AC nanoparticle group, and a large number of DEGs were involved in molting, energy metabolism, stress tolerance, and germ cell development. Moreover, KEGG analysis revealed that many DEGs were assigned into signaling pathways related to metabolic growth and immune stress. These results showed that exposure to TiO2/AC nanoparticle will result in the changes of enzyme activity and routine mRNA expression, suggesting that TiO2/AC nanoparticle which existed in aquatic environment might affect the physiology of M. rosenbergii. This study will provide significant information for the evaluation of nanomaterial toxicity on aquatic crustaceans.

Effects of green tide on microbial communities in waters of the Jiangsu coastal area, China.

Recently, green tide outbreaks have resulted in severe coastal ecology and economic effects in China. Jiangsu coastal areas are usually the site of early green tide outbreaks. To clarify the effects of green tide outbreaks in Jiangsu coastal areas, this study analyzed microbial communities during green tide-free and green tide outbreak periods (May and July, respectively) through 16S rDNA sequencing. Sequences were clustered into 4117 operational taxonomic units (OTUs), 1044 and 3834 of which were obtained from the May and July groups, respectively. Redundancy analysis indicated that green tide occurrence was closely associated with the temperature, pH, and concentrations of various nutrients. Diversity analysis revealed that the July group had a richer microbial community than the May group, in agreement with the results of propagule culture. Moreover, comparative analysis revealed that samples in the May and July groups clustered together. According to Megan analysis, the May group had much more Psychrobacter, Sulfitobacter, and Marinomonas than the July group, whereas the other genera were predominantly found in July, such as Ascidiacerhabitans, Synechococcus Hydrotalea, and Burkholderia-Paraburkholderia. These findings suggest that green tide outbreaks affect marine microbial communities, and detecting the changes in the identified genera during green tide outbreaks may contribute to green tide forecasting. PRACTITIONER POINTS: Jiangsu coastal areas are usually the site of early green tide outbreaks. Green tide occurrence was related to the concentrations of various nutrients. Microbial species and community structure significantly changed after green tide outbreak.

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.

Transcriptome analysis of the post-larvae of giant freshwater prawn (Macrobrachium rosenbergii) after IAG gene knockdown with microRNA interference.

The insulin-like androgenic gland hormone gene (IAG) of crustaceans plays pivotal roles in the regulation of sex differentiation. MicroRNAs (miRNAs) are a class of short, non-coding RNAs that function as post-transcriptional gene regulators. However, little information about the regulatory relationship between miRNA and Macrobrachium rosenbergii IAG (MrIAG) were exposed. In this study, we used the 3' untranslated region (UTR) of MrIAG to predict potential target sites of miRNAs. The results showed that miR-184 has one target site in the 3'UTR of MrIAG. Dual-luciferase report assay in vitro confirmed that miR-184 can significantly down-regulate MrIAG expression. Besides, we constructed mutant plasmids of 3'UTR of MrIAG. The result displayed that after co-transfection of mutant plasmids and miR-184 agomir, the activity of luciferase was not affected compared to the control. These results indicated that miR-184 could directly regulate MrIAG. In addition, we found that overexpression of miR-184 in M. rosenbergii can lead to significant changes in the transcription level of genes. Compared with control group, we identified 1510 differentially expressed genes (DEGs) in the miR-184 injection group. Some DEGs were involved in sex differentiation, gonad development, growth and molting were found. qRT-PCR verification was performed on eight DEGs randomly, and the results showed that the expression level of sex-, growth-, and metabolism-related genes changed significantly after MrIAG gene knockdown. Collectively, findings from this study suggest that miR-184, by mediating IAG expression, may be involved in many physiological processes in M. rosenbergii. The current study lays a basic understanding for short-term silencing of MrIAG with miR-184, and facilitates miRNA function analysis in M. rosenbergii in future.

Comparative transcriptome analysis of the gills and hepatopancreas from Macrobrachium rosenbergii exposed to the heavy metal Cadmium (Cd2+).

Heavy metal Cadmium (Cd2+) pollution has become a severe environmental problem for aquatic organisms. In crustaceans, gills (Gi) and hepatopancreas (Hp) play a vital role in the toxicology. However, in Macrobrachium rosenbergill, there are few researches about gill and hepatopancreases responding to Cd2+ stress at a molecular level. In this study, transcriptomic analysis was applied to characterize gene expression profiles of gills and hepatopancreas of M. rosenbergill after Cd2+ exposure for 0 h, 3 h and 3 d. Six cDNA libraries (Gi 0 h, Gi 3 h, Gi 3 d, Hp 0 h, Hp 3 h, and Hp 3 d) were constructed and a total of 66,676 transcripts and 48,991 unigenes were annotated. Furthermore, differentially expressed genes (DEGs) were isolated by comparing the Cd2+ treated time-point libraries (3 h and 3 d group) with the control library (0 h group). The results showed that most of the DEGs were down-regulated after Cd2+ exposure and the number of DEGs among gill groups were significantly higher than those among hepatopancreas groups. GO functional and KEGG pathway analysis suggested many key DEGs in response to the Cd2+ stress, such as metallothionein and Hemocyanin. Additionally, a total of six DEGs were randomly selected to further identify their expressional profile by qPCR. The results indicated that these DEGs were involved in the response to Cd2+. This comparative transcriptome provides valuable molecular information on the mechanisms of responding to Cd2+ stress in M. rosenbergii, which lays the foundation for further understanding of heavy metal stress.

Transcriptional changes revealed water acidification leads to the immune response and ovary maturation delay in the Chinese mitten crab Eriocheir sinensis.

Nowadays, due to increasing carbon dioxide released, water acidification poses a series of serious impacts on aquatic organisms. To evaluate the effects of water acidification on crustaceans, we focused on the Chinese mitten crab Eriocheir sinensis, which is a spawning migration and farmed species in China. Based on histological and oocyte transparent liquid observation, we found that the acidified environment significantly delayed the ovarian maturation of E. sinensis. Moreover, RNA-seq was applied to obtain gene expression profile from the crab's gills and ovaries in response to acidified environment. Compared with control groups, a total of 5471 differentially expressed genes (DEGs) were identified in acidified gills and 485 DEGs were identified in acidified ovaries. Enrichment analysis indicated that some pathways also responded to the acidified environment, such as PI3K-Akt signaling pathway, Chemokine signaling pathway, apoptosis, and toll-like receptor signaling pathway. Subsequently, some DEGs involved in immune response (ALF, Cathepsin A, HSP70, HSP90, and catalase) and ovarian maturation (Cyclin B, Fem-1a, Fem-1b, and Fem-1c) were selected to further validate the influence of water acidification on gene expression by qRT-PCR. The results showed that the expression level of immune-related genes was significantly increased to response to the water acidification, while the ovarian maturation-related genes were significantly decreased. Overall, our data suggested that E. sinensis was sensitive to the reduced pH. This comparative transcriptome also provides valuable molecular information on the mechanisms of the crustaceans responding to acidified environment.

Identification, characterization and mRNA transcript abundance profiles of the carboxylesterase (CXE5) gene in Eriocheir sinensis suggest that it may play a role in methyl farnesoate degradation.

The sesquiterpenoid methyl farnesoate (MF) is a de-epoxidized form of insect juvenile hormone (JH) III in crustaceans, and its precise titer plays important roles in regulating many critical physiological processes, including reproduction and ovarian maturation. Understanding the synthetic and degradation pathways of MF is equally important for determining how to maintain MF titers at appropriate levels and thus for potential applications in crab aquaculture. Although the synthetic pathway of MF has been well established, little is known about MF degradation. Previous research proposed that specific carboxylesterases (CXEs) that degrade MF in crustaceans are conserved from those of JH III. In this study, we identified a novel Es-CXE5 gene from Eriocheir sinensis. The Es-CXE5 protein contains some conserved motifs, including catalytic triad and oxyanion hole, which are characteristics of the biologically active CXE family. The phylogenetic analysis showed that Es-CXE5 belongs to the hormone/semiochemical processing group of the CXE family. Moreover, Tissue and stage-specific expression results suggested that Es-CXE5 expression in hepatopancreas was highest and associated with the hemolymph MF titer. Furthermore, Es-CXE5 mRNA transcripts were detected in both in vitro and in vivo experiments and ESA experiment in the hepatopancreas and ovary. The results of this study showed that Es-CXE5 mRNA abundance in the hepatopancreas was notably induced by MF addition but had no effect on the ovary. Taken together, our results suggest that Es-CXE5 may degrade MF in the hepatopancreas and may thus be involved in ovarian development in E. sinensis.

Identification of a novel germ cell marker MnTdrd from the oriental river prawn Macrobrachium nipponense.

Germ cell-specific genes play an important role in establishing the reproductive system in sexual organisms and have been used as valuable markers for studying gametogenesis and sex differentiation. Previously, we isolated a vasa transcript as a germ cell marker to trace the origin and migration of germ cells in the oriental river prawn Macrobrachium nipponense. Here, we identified a new germ cell-specific marker MnTdrd RNA and assessed its temporal and spatial expression during oogenesis and embryogenesis. MnTdrd transcripts were expressed in high abundance in unfertilized eggs and embryos at cleavage stage and then dropped significantly during late embryogenesis, suggesting that MnTdrd mRNA is maternally inherited. In situ hybridization of ovarian tissue showed that MnTdrd mRNA was initially present in the cytoplasm of previtellogenic oocyte and localized to the perinuclear region as the accumulation of yolk in vitellogenic oocyte. Whole-mount in situ hybridization of embryos showed that MnTdrd-positive signals were only localized in one blastomere until 16-cell stage. In the blastula, there were approximately 16 MnTdrd-positive blastomeres. During embryonized-zoea stage, the MnTdrd-positive cells aggregated as a cluster and migrated to the genital rudiment which would develop into primordial germ cells (PGCs). The localized expression pattern of MnTdrd transcripts resembled that of the previously identified germ cell marker vasa, supporting the preformation mode of germ cell specification. Therefore, we concluded that MnTdrd, together with vasa, is a component of the germ plasm and might have critical roles in germ cell formation and differentiation in the prawn. Thus, MnTdrd can be used as a novel germ cell marker to trace the origin and migration of germ cells.

Molecular characterization of ovary-specific gene Mrfem-1 and siRNA-mediated regulation on targeting Mrfem-1 in the giant freshwater prawn, Macrobrachium rosenbergii.

Characterized by ankyrin repeat motifs, the feminization-1 (fem-1) gene plays an essential role in sex determination/differentiation in Caenorhabditis elegans. However, there are only a few reports on fem-1 in crustaceans. In this study, a fem-1 gene (Mrfem-1) was first isolated from the giant freshwater prawn Macrobrachium rosenbergii. The full-length cDNA of Mrfem-1 was 2607 bp long, containing an open reading frame encoding 615 amino acids, and presenting eight ankyrin repeats. The full-length cDNA has been submitted to GenBank with the accession no. MT160093. According to the RT-PCR results, Mrfem-1 was exclusively expressed in the ovary. The expression level of Mrfem-1 had increased with ovarian maturation and reached the highest peak at vitellogenic stage. In situ hybridization results showed that positive signals were concentrated in the cytoplasm of previtellogenic stage, and scattered in the cytoplasm and follicular cells at vitellogenic stage, suggesting that Mrfem-1 might be associated with ovarian maturation. Moreover, two effective siRNAs targeting Mrfem-1 were found and their effectiveness verified in vitro. These results on Mrfem-1 will help us to better understand the fem family and provide a new resource for subsequent investigations of siRNA-mediated regulation on ovarian development in M. rosenbergii.

Alternative splice variants and differential relative abundance patterns of vasa mRNAs during gonadal development in the Chinese mitten crab Eriocheir sinensis.

Gonadal development usually involves alternative splicing of sex-related genes. Vasa, a highly conserved ATP-dependent RNA helicase present mainly in germ cells, has an important function in gonadal development. As an important sex-related gene, recent evidence indicates that different splice variants of vasa exist in many species. In this study, there was identification of two types of vasa splice variants in the Chinese mitten crab Eriocheir sinensis, termed Esvasa-l and Esvasa-s, respectively. Furthermore, splice variants of Esvasa-s were sub-divided into Esvasa-s1, Esvasa-s2, Esvasa-s3, Esvasa-s4, and Esvasa-s5, based on differing numbers of TGG repeats. Results from genomic structure analyses indicated that these forms are alternatively spliced transcripts from a single vasa gene. Results from tissue distribution assessments indicate the vasa splice variants were exclusively expressed in the gonads of male and female adult crabs. In situ hybridization results indicate Esvasa mRNA was mainly present in the cytoplasm of previtellogenic oocytes. As oocyte size increased, relative abundance of Esvasa mRNA decreased and became distributed near the cellular membrane. The Esvasa mRNA was not detectable in mature oocytes. In testis, Esvasa mRNA was detected in spermatids and spermatozoa, but not in spermatogonia and spermatocytes. Notably, results from qPCR analysis of Esvasa-l and Esvasa-s indicate there are different relative proportions during gametogenesis, implying that splice variants of the Esvasa gene may have different biological functions during crab gonadal development.

Identification and profiling of microRNAs during gonadal development in the giant freshwater prawn Macrobrachium rosenbergii.

As post-transcriptional regulators, microRNAs (miRNAs) play an important role in growth and reproductive processes. So far, there is limited information regarding crustacean miRNAs. To explore the potential role of miRNAs in the gonadal development of the prawn Macrobrachium rosenbergii, we constructed seven small RNA libraries from ovarian and testicular tissues at various stages using somatic tissue as the control. A total of 1,954 known and 129 novel miRNAs were retrieved. By comparing differentially expressed miRNAs (DEMs) between testes and ovaries, forty-one miRNAs were identified with sex-biased expression patterns, including 17 ovary-biased and 24 testis-biased patterns. Furthermore, the putative target genes of the sex-biased miRNAs, such as cyclin L1, mitogen-activated protein kinase 7 (MAPK 7), heat shock protein (HSP), and zinc finger protein, were significantly enriched in many reproduction-related pathways including the Gonadotropin-releasing hormone (GnRH) pathway, glycolysis, gluconeogenesis pathway, ovarian steroidogenesis, estrogen signaling pathway, MAPK pathway, Wnt pathway, and insulin signaling pathway, implicating potential regulatory roles of miRNAs in reproduction. These data aid in the further investigation of the mechanism of gonadal development and reproductive regulation mediated by miRNA in M. rosenbergii.

Construction of a Genomic Bacterial Artificial Chromosome (BAC) Library for the Prawn Macrobrachium rosenbergii and Initial Analysis of ZW Chromosome-Derived BAC Inserts.

Knowledge on sex determination has proven valuable for commercial production of the prawn Macrobrachium rosenbergii due to sex dimorphism of the male and female individuals. Previous studies indicated that prawn sex is determined by a ZW-ZZ chromosomal system, but no genomic information is available for the sex chromosome. Herein, we constructed a genomic bacterial artificial chromosome (BAC) library and identified the ZW-derived BAC clones for initial analysis of the sex chromosomal DNA sequence. The arrayed BAC library contains 200,448 clones with average insert size of 115.4 kb, corresponding to ∼ 4× coverage of the estimated 5.38 Gb genome. Based on a short female-specific marker, a Z- and a W-fragment were retrieved with the genomic walking method. Screening the BAC library using a ZW-specific marker as probe resulted in 12 positive clones. From these, a Z-derived (P331M17) and a W-derived (P122G2) BAC clones were randomly selected and sequenced by PacBio method. We report the construction of a large insert, deep-coverage, and high-quality BAC library for M. rosenbergii that provides a useful resource for positional cloning of target genes, genomic organization, and comparative genomics analysis. Our study not only confirmed the ZW/ZZ system but also discovered sex-linked genes on ZW chromosomes for the first time, contributing to a comprehensive understanding of the genomic structure of sex chromosomes in M. rosenbergii.

Transcriptome analysis of the brain of the Chinese mitten crab, Eriocheir sinensis, for neuropeptide abundance profiles during ovarian development.

Neuropeptides, important messenger molecules, regulate various physiological processes, such as growth, development, and reproduction. In the present study, cDNA libraries from brains of E. sinensis were constructed and sequenced using the Illumina technique for transcript analysis and neuropeptides discovery. There were 233,887 transcripts assembled for 194,286 unigenes. According to the annotations of NCBI non-redundant protein (NR) database, 2487 (11.31%) unigenes were annotated successfully. In total, 1273 transcripts were assigned to the "signal transduction mechanisms" category using KOG analysis. The results of KEGG indicate signal transduction and translation pathways were the dominant and enriched signal pathways. Additionally, results indicated C2H2 was the main transcription factor (TF) family. Analysis of the assembled transcripts indicated there were 22 neuropeptide transcripts, such as allatostatin, crustacean female sex hormone, crustacean hyperglycemic hormone, diuretic hormone 31, and eclosion hormone. The detection of these neuropeptides provide for a basic understanding for future study of functions in development, reproduction, and sexual maturation in crustaceans.

Molecular cloning and characterization of a gonadotropin-releasing hormone receptor homolog in the Chinese mitten crab, Eriocheir sinensis.

As an essential mediator in the Gonadotropin-releasing hormone (GnRH) signaling pathway, GnRH receptor (GnRHR) coupled to GnRH, plays an important role in activating the downstream pathway after stimulating a series of cascades to regulate reproduction. To detect the existence of GnRHR and potential GnRH signaling pathway, we cloned and characterized GnRHR in the Chinese mitten crab, Eriocheir sinensis (named EsGnRHR). The full-length EsGnRHR cDNA is 2038 bp in length, including an open reading frame (ORF) of 1566 bp, a 57 bp 5'-untranslated region (5'-UTR) and a 415 bp 3'-UTR. Prediction of transmembrane domains in protein sequence revealed that the EsGnRHR protein contained seven hydrophobic transmembrane regions (TMs). Reverse transcription PCR revealed that EsGnRHR was mainly expressed in the thoracic nerve group and ovary, and weakly distributed in the testis and brain. In situ hybridization further demonstrated that EsGnRHR mRNA was localized at the protocerebrum and deutocerebrum. In the ovary and testis, the hybridization signal was dominantly at the earlier developmental stages. The signal was mainly localized in the cytoplasm cell in the ovary, and in the epithelium cell in the testis. During the different stages of gonadal development, EsGnRHR displayed increasing trends in both female and male when analyzed by quantitative real-time PCR, suggesting that EsGnRHR was involved in controlling gonadal development. Our study provides important information for further research on the molecular mechanisms underlying crab development.

Characterization of the complete mitochondrial genome sequence of Charybdis bimaculata and phylogenetic analysis.

In this study, the complete mitogenome of Charybdis bimaculata was sequenced and annotated. The mitogenome is 15,441 bp in length, containing 37 classical eukaryotic mitochondrial regions (13 typical protein-coding genes (PCGs), 22 tRNA genes, two rRNA genes) and a non-coding control region. Most of the genes are initiated with ATA, ATG, and ATT, though GTG is also used as an initiation codon. Twelve PCGs stop with complete termination codon TAA and TAG, while Cob uses incomplete codon (T-). The phylogenetic relationships based on 13 PCGs show that C. bimaculata clusters closest to C. fariata and C. natator.