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

Study on reproductive endocrine disturbance and DNA damage mechanism of female Ruditapes philippinarum under Benzo[a]pyrene stress

The reproductive toxicity of polycyclic aromatic hydrocarbons (PAHs) in aquatic organisms has attracted increasing attention from scholars. Currently, research in this field primarily focuses on vertebrates such as zebrafish and other model species. However, there is still a significant knowledge gap in the toxicity of PAHs to invertebrates and its potential mechanisms. Benzo[a]pyrene (B[a]P) is one of the most representative PAHs. In this study, female Ruditapes philippinarum (R. philippinarum) was treated with B[a]P concentrations of 0, 0.8, 4, and 20 μg/L to investigate reproductive indicators in the proliferative, growth, mature, and spawn stages. Transcriptomics was used to investigate the expression of genes associated with the reproductive endocrine system, DNA repair, autophagy, apoptosis, and ovarian development at different reproductive stages. Our results suggested that B[a]P disrupted the endocrine system by interfering with the production of steroid hormones and the transmission of estrogen signals in female R. philippinarum. The structure of the ovarian DNA duplex is severely damaged under the stress of B[a]P, and a series of cellular responses caused by DNA damage are also interfered. Additionally, we observed a reduction in the gonadosomatic index (GSI) and mature oocytes numbers after B[a]P exposed. Tissue section indicated that severe damage to the ovarian structure at mature and spawn stages. In conclusion, this study combined transcriptomic and toxicological to explore the negative effects on ovarian development induced by B[a]P, focusing on reproductive endocrine disturbance and DNA damage.

The miRNAs let-7b and miR-141 Coordinately Regulate Vitellogenesis by Modulating Methyl Farnesoate Degradation in the Swimming Crab Portunus trituberculatus

Methyl farnesoate (MF), a crucial sesquiterpenoid hormone, plays a pivotal role in the reproduction of female crustaceans, particularly in the vitellogenesis process. Despite extensive research on its functions, the molecular mechanisms that regulate MF levels during the vitellogenic phase remain largely elusive. This study investigates the roles of microRNAs (miRNAs), significant post-transcriptional regulators of gene expression, in controlling MF levels in the swimming crab Portunus trituberculatus. Through bioinformatic analysis, four miRNAs were identified as potential regulators targeting two genes encoding Carboxylesterases (CXEs), which are key enzymes in MF degradation. Dual luciferase reporter assays revealed that let-7b and miR-141 suppress CXE1 and CXE2 expression by directly binding to their 3' UTRs. In vivo overexpression of let-7b and miR-141 significantly diminished CXE1 and CXE2 levels, consequently elevating hemolymph MF and enhancing vitellogenin expression. Spatiotemporal expression profile analysis showed that these two miRNAs and their targets exhibited generally opposite patterns during ovarian development. These findings demonstrate that let-7b and miR-141 collaboratively modulate MF levels by targeting CXEs, thus influencing vitellogenesis in P. trituberculatus. Additionally, we found that the expression of let-7b and miR-141 were suppressed by MF, constituting a regulatory loop for the regulation of MF levels. The findings contribute novel insights into miRNA-mediated ovarian development regulation in crustaceans and offer valuable information for developing innovative reproduction manipulation techniques for P. trituberculatus.

17β-Estradiol Induced Sex Reversal and Gonadal Transcriptome Analysis in the Oriental River Prawn (Macrobrachium nipponense): Mechanisms, Pathways, and Potential Harm

Sex reversal induced by 17β-estradiol (E₂) has shown the potential possibility for monoculture technology development. The present study aimed to determine whether dietary supplementation with different concentrations of E₂ could induce sex reversal in M. nipponense, and select the sex-related genes by performing the gonadal transcriptome analysis of normal male (M), normal female (FM), sex-reversed male prawns (RM), and unreversed male prawns (NRM). Histology, transcriptome analysis, and qPCR were performed to compare differences in gonad development, key metabolic pathways, and genes. Compared with the control, after 40 days, feeding E₂ with 200 mg/kg at PL25 (PL: post-larvae developmental stage) resulted in the highest sex ratio (female: male) of 2.22:1. Histological observations demonstrated the co-existence of testis and ovaries in the same prawn. Male prawns from the NRM group exhibited slower testis development without mature sperm. RNA sequencing revealed 3702 differentially expressed genes (DEGs) between M vs. FM, 3111 between M vs. RM, and 4978 between FM vs. NRM. Retinol metabolism and nucleotide excision repair pathways were identified as the key pathways for sex reversal and sperm maturation, respectively. Sperm gelatinase (SG) was not screened in M vs. NRM, corroborating the results of the slice D. In M vs. RM, reproduction-related genes such as cathepsin C (CatC), heat shock protein cognate (HSP), double-sex (Dsx), and gonadotropin-releasing hormone receptor (GnRH) were expressed differently from the other two groups, indicating that these are involved in the process of sex reversal. Exogenous E₂ can induce sex reversal, providing valuable evidence for the establishment of monoculture in this species.

Small RNA-seq and hormones in the testes of dwarf hamsters (Cricetulus barabensis) reveal the potential pathways in photoperiod regulated reproduction

Photoperiod regulates the functions and development of gonadal organs of seasonally breeding animals, resulting in breeding peaks in specific seasons. miRNA plays an important role in the regulation of testicular physiological functions. However, the relationship between photoperiods and miRNA levels in testes has yet to be conclusively determined. We investigated testicular miRNA of striped dwarf hamster (Cricetulus barabensis) responses to different photoperiods (long daylength [LD], moderate daylength [MD], and short daylength [SD]) and the potential pathways involved in photoperiod regulated reproduction. Testicular weights and reproductive hormone levels were measured in each of photoperiod treatments after 30 days. The concentrations of testosterone (T) and dihydrogen testosterone (DHT) in testes and Gonadotropin-releasing hormone (GnRH), follicle-stimulating hormone (FSH), and luteinizing hormone (LH) in serum were higher in MD than in the other two groups. Testicular weights were heaviest in MD. Small RNA-seq was performed for the testes of hamsters in three groups. A total of 769 miRNAs were identified, of which 83 were differentially expressed between LD, MD, and SD. GO and KEGG analysis of target genes revealed that some miRNAs influence testicular activities by regulating the pathways related to cell apoptosis and metabolism. Gene expression pattern analysis showed that the MAPK signaling pathway may be the core pathway for photoperiodic regulation of reproduction. These results suggest that moderate daylength is more suitable for hamster reproduction while long daylength and short daylength may regulate reproduction through different molecular pathways.

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.

Integrated transcriptome analysis of immune-related mRNAs and microRNAs in Macrobrachium rosenbergii infected with Spiroplasma eriocheiris

Macrobrachium rosenbergii (M. rosenbergii), is a major aquaculture species in China and Southeast Asia. However, infection with Spiroplasma eriocheiris (S. eriocheiris) has caused huge economic losses to the cultivation of M. rosenbergii. Currently, there are few reports on the immune response mechanism of M. rosenbergii that are infected with S. eriocheiris. To clarify the immune response mechanism of M. rosenbergii infected with S. eriocheiris, the key immune genes which respond to the infection with the pathogen and the regulation of related microRNAs (miRNAs) on them were identified. In this study, the mRNA and miRNA transcriptome of hepatopancreas of M. rosenbergii at different infection stages were analyzed using high-throughput sequencing and qRT-PCR. In the mRNA transcriptome, 27,703 and 33,402 genes were expressed in healthy and susceptible M. rosenbergii, respectively. By digital gene-expression profiling analysis, 23,929 and 24,325 genes were expressed, and 223 and 373 genes were significantly up-regulated and down-regulated, respectively. A total of 145 key genes related to Toll, IMD, JAK/STAT and MAPK were excavated from the transcriptome. In the miRNA transcriptome, 549 miRNAs (Conserved: 41, PN-type: 83, PC-type: 425) were sequenced, of which 87 were significantly up-regulated and 23 were significantly down-regulated. Among the related immune pathways, there are 259 miRNAs involved in the regulation of target genes in the Toll and IMD pathways, 231 JAK/STAT pathways and 122 MAPK pathways. qRT-PCR differential detection of immune-related miRNAs and mRNAs showed that 22 miRNAs with significant differences (P < 0.05) such as mro-miR-100, PC-mro-3p-27 and PN-mro-miR-316 had corresponding regulatory relationships with 22 important immune genes such as TLR2, TLR3, TLR4, TLR5, MyD88, Pelle and Relish in different stages after infection. In this study, the immune genes and related regulatory miRNAs of M. rosenbergii in response to S. eriocheiris infection were obtained. The results can provide basic data to further reveal the immune defense mechanism of M. rosenbergii against S. eriocheiris infection.

Co-expression network analysis identifies novel molecular pathways associated with cadmium and pyriproxyfen testicular toxicity in Gammarus fossarum

Omics approaches are continuously providing new clues on the mechanisms of action of contaminants in species of environmental relevance, contributing to the emergence of molecular ecotoxicology. Co-expression network approaches represent a suitable methodological framework for studying the rich content of omics datasets. This study aimed to find evidence of key pathways and proteins related to the testicular toxicity in the sentinel crustacean species Gammarus fossarum exposed to endocrine disruptors using a weighted protein co-expression network analysis. From a shotgun proteomics dataset of male gonads of G. fossarum organisms exposed to cadmium (Cd), pyriproxyfen (Pyr) and methoxyfenozide (Met) in laboratory conditions, four distinct modules were identified as significantly correlated to contaminants' exposure. Protein set enrichment analysis identified modules involved in cytoskeleton organization and oxidative stress response associated with the Cd exposure. The module associated with Pyr exposure was associated with endoplasmic reticulum stress (ER) response, and the module correlated with Met exposure was characterized by a significant proportion of amphipod-restricted proteins whose functions are still not characterized. Our results show that co-expression networks are efficient and adapted tools to identify new potential mode of actions from environmental sentinel species, such as G. fossarum, using a proteogenomic approach, even without an annotated genome.

Identification and characterization of sex-biased and differentially expressed miRNAs in gonadal developments of the Chinese mitten crab, Eriocheir sinensis

MicroRNA (miRNA) is a posttranscriptional downregulator that plays a vital role in a wide variety of biological processes. In this study, we constructed five ovarian and testicular small RNA libraries using two somatic libraries as reference controls for the identification of sex-biased miRNAs and gonadal differentially expressed miRNAs (DEMs) of the Chinese mitten crab, Eriocheir sinensis. A total of 535 known and 243 novel miRNAs were identified, including 312 sex-biased miRNAs and 402 gonadal DEMs. KEGG pathway analysis showed that DEM target genes were statistically enriched in MAPK, Wnt, and GnRH signaling pathway, and so on. A number of the sex-biased miRNAs target genes associated with sex determination/differentiation, such as IAG, Dsx, Dmrt1, and Fem1, while others target the genes related to gonadal development, such as P450s, Wnt, Ef1, and Tra-2c. Dual-luciferase reporter assay in vitro further confirmed that miR-34 and let-7b can downregulate IAG expression, miR-9-5p, let-7d, let-7b, and miR-8915 can downregulate Dsx. Taken together, these data strongly suggest a potential role for the sex-biased miRNAs in sex determination/differentiation and gonadal development in the crab.

Integrated mRNA and miRNA expression profile analysis of female and male gonads in Hyriopsis cumingii

Hyriopsis cumingii is an important species for freshwater pearl cultivation in China. In terms of pearl production, males have larger pearls and better glossiness than females, but there are few reports focusing on the sex of H. cumingii. In this study, six mRNA and six microRNA (miRNA) libraries were prepared from ovaries and testes. Additionally, 28,502 differentially expressed genes (DEGs) and 32 differentially expressed miRNAs (DEMs) were identified. Compared with testis, 14,360 mRNAs and 20 miRNAs were up-regulated in ovary, 14,142 mRNAs and 12 miRNAs were down-regulated. In DEGs, the known genes related to sex determinism and/or differentiation were also identified, such as DMRT1, SOX9, SF1 for males, FOXL2 for females, and other potentially significant candidate genes. Three sex-related pathways have also been identified, which are Wnt, Notch, and TGF-beta. In 32 DEMs, the three miRNAs (miR-9-5p, miR-92, miR-184) were paid more attention, they predicted 28 target genes, which may also be candidates for sex-related miRNAs and genes. Differential miRNAs target genes analysis reveals the pathway associated with oocyte meiosis and spermatogenesis. Overall, the findings of the study provide significant insights to enhance our understanding of sex differentiation and/or sex determination mechanisms for H. cumingii.

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.

Identification and Integrated Analysis of MicroRNA and mRNA Expression Profiles During Agonistic Behavior in Chinese Mitten Crab (Eriocheir sinensis) Using a Deep Sequencing Approach

As a commercially important species, the Chinese mitten crab (Eriocheir sinensis) has been cultured for a long time in China. Agonistic behavior often causes limb disability and requires much energy, which is harmful to the growth and survival of crabs. In this paper, we divided crabs into a control group (control, no treatment) and an experimental group (fight, agonistic behavior after 1 h) and then collected the thoracic ganglia (TG) to extract RNA. Subsequently, we first used a deep sequencing approach to examine the transcripts of microRNAs (miRNAs) and messenger RNAs (mRNAs) in E. sinensis displaying agonistic behavior. According to the results, we found 29 significant differentially expressed miRNAs (DEMs) and 116 significant differentially expressed unigenes (DEGs). The DEMs esi-miR-199a-5p, esi-let-7d, esi-miR-200a, and esi-miR-200b might participate in the regulation of agonistic behavior by mediating neuroregulation and energy metabolism. Focusing on the transcripts of the mRNAs, the renin–angiotensin system (RAS) Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway might be involved in the regulation of agonistic behavior through glucose metabolism as this pathway was significantly enriched with DEGs. Besides, an integrated analysis of the miRNA and mRNA profiles revealed that the retinoid X receptor (RXR) was also involved in visual signal transduction, which was important for agonistic behavior. In addition, four vital agonistic behavior-related metabolic pathways, including the cAMP signaling, MAPK, protein digestion and absorption, and fatty acid metabolism pathways, were significantly enriched with the predicted target unigenes. In conclusion, the findings of this study might provide important insight enhancing our understanding of the underlying molecular mechanisms of agonistic behavior in E. sinensis.

Integrative Proteomic and MicroRNA Analysis: Insights Into Mechanisms of Eyestalk Ablation-Induced Ovarian Maturation in the Swimming Crab Portunus trituberculatus

Eyestalk ablation is the most common method to induce ovarian maturation in decapod crustacean aquaculture, but it jeopardizes broodstock survival and larvae production. It is important to understand the molecular basis underlying the maturation triggered by ablation and thereby develop an alternative measure for maturation manipulation. In this study, we investigate alterations of ovarian proteome and miRNA profile after ablation in a commercially important marine crab Portunus trituberculatus. Quantitative proteomic analysis using iTRAQ reveals that 163 proteins are differentially expressed following ablation, and modulation of methyl farnesoate metabolism and activation of calcium signaling may play important roles in the ovarian maturation induced by ablation. miRNA expression profiling identifies 31 miRNAs that show statistically significant changes. Integration of miRNA and proteome expression data with miRNA target prediction algorithms generates a potential regulatory network consisting of 26 miRNAs and 30 proteins linked by 71 possible functional associations. The miRNA-protein network analysis suggests that miRNAs are involved in promoting ovarian maturation by controlling expression of proteins related to methyl farnesoate synthesis, calcium signals, and energy metabolism. Experimental validation and temporal expression analysis indicate multiple miRNAs can act synergistically to regulate expression of Farnesoic acid O-methyltransferase and Calmodulin. Our findings provide new insights for elucidating the mechanisms underlying eyestalk ablation-induced ovarian maturation and could be useful for devising an alternative technique for manipulating reproduction in P. trituberculatus and other decapods.

Comparative Transcriptome Analysis of Gonads for the Identification of Sex-Related Genes in Giant Freshwater Prawns (Macrobrachium Rosenbergii) Using RNA Sequencing

The giant freshwater prawn (Macrobrachium rosenbergii) exhibits sex dimorphism between the male and female individuals. To date, the molecular mechanism governing gonadal development was unclear, and limited data were available on the gonad transcriptome of M. rosenbergii. Here, we conducted comprehensive gonadal transcriptomic analysis of female (ZW), super female (WW), and male (ZZ) M. rosenbergii for gene discovery. A total of 70.33 gigabases (Gb) of sequences were generated. There were 115,338 unigenes assembled with a mean size of 1196 base pair (bp) and N50 of 2195 bp. Alignment against the National Center for Biotechnology Information (NCBI) non-redundant nucleotide/protein sequence database (NR and NT), the Kyoto Encyclopedia of Genes and Genomes (KEGG) database, SwissProt database, Protein family (Pfam), Gene ontology (GO), and the eukaryotic orthologous group (KOG) database, 36,282 unigenes were annotated at least in one database. Comparative transcriptome analysis observed that 10,641, 16,903, and 3393 genes were significantly differentially expressed in ZW vs. ZZ, WW vs. ZZ, and WW vs. ZW samples, respectively. Enrichment analysis of differentially expressed genes (DEGs) resulted in 268, 153, and 42 significantly enriched GO terms, respectively, and a total of 56 significantly enriched KEGG pathways. Additionally, 23 putative sex-related genes, including Gtsf1, IR, HSP21, MRPINK, Mrr, and other potentially promising candidate genes were identified. Moreover, 56,241 simple sequence repeats (SSRs) were identified. Our findings provide a valuable archive for further functional analyses of sex-related genes and future discoveries of underlying molecular mechanisms of gonadal development and sex determination.