Molecular characterization and expression profiles of cdc2 and cyclin B during oogenesis and spermatogenesis in green mud crab (Scylla paramamosain)

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.

A novel Dmrt gene of crustacean: functional analysis of idmrt-2 gene in the male reproductive system from Scylla paramamosain

The Dmrt (double-sex and mab-3 related transcription factor) gene family is considered to be a highly conserved gene family related to sex determination and sexual differentiation across species. In order to better understand the role of the idmrt-2 gene in gonad development in Scylla paramamosain, the idmrt-2 gene was cloned and analyzed. The cDNA contains a 1659 bp ORF region encoding 552 amino acids. The qRT-PCR results showed that idmrt-2 was significantly more expressed in the testis than in other tissues (p<0.05). The expression of idmrt-2 was highest in the spermatids stage (T2 stage), followed by the mature sperms stage (T3 stage) and significantly higher than in the spermatocytes stage (T1 stage) (p<0.05) during testicular development and the expression difference was not significant in different stages of ovarian development. RNAi studies revealed that after idmrt-2 was knocked down, the expression of Dmrt-like and foxl-2 genes in the testis decreased, as well as IAG expression in the androgenic gland. The findings suggest that idmrt-2 may be an IAG regulator and involved in testicular development.

Comparative Transcriptome Reveals the Potential Modulation Mechanisms of Spfoxl-2 Affecting Ovarian Development of Scylla paramamosain

Previously, we reported the identification, tissue distribution and confirmed the roles of Spfoxl-2 in regulating vitellogenin (vtg) expression in Scylla paramamosain. Here, we primally analyzed its potential target genes in the ovary with RNAi and RNA-Seq technology. By comparing the transcriptome data of two groups (ovaries that injected with EGFP and Foxl-2 siRNA, respectively), we found 645 DEGs (differentially expressed genes), including several conserved crucial genes involved in ovarian development, such as vtg, vitellogenin receptor (vtgR), adenylate cyclase (AC), cyclinB, and cell division cycle 2 (cdc2). In addition, these DEGs were also enriched in pathways related to ovary development, including relaxin signaling pathway, ovarian steroidogenesis, and progesterone-mediated oocyte maturation. Moreover, several genes were selected for qRT-PCR to validate the accuracy of the bioinformatic result. To the best of our knowledge, the current study was the first report about foxl-2 function through comparative transcriptome analysis in crustacean species, which identified not only relevant genes and pathways involved in ovarian development of S. paramamosain, but also provided new insights into the regulatory mechanisms of foxl-2 at the molecular level in crustacean.

Morphological and molecular features of early regeneration in the marine annelid Ophryotrocha xiamen

Regeneration capability varies in the phylum Annelida making them an excellent group to investigate the differences between closely related organisms. Several studies have described the process of regeneration, while the underlying molecular mechanism remains unclear, especially during the early stage (wound healing and blastema formation). In this study, the newly identified Ophryotrocha xiamen was used to explore the early regeneration. The detailed morphological and molecular analyses positioned O. xiamen within 'labronica' clade. We analyzed the morphological changes during regeneration process (0-3 days post amputation) and molecular changes during the early regeneration stage (1 day post amputation). Wound healing was achieved within one day and a blastema formed one day later. A total of 243 DEGs were mainly involved in metabolism and signal transduction. Currently known regeneration-related genes were identified in O. xiamen which could help with exploring the functions of genes involved in regeneration processes. According to their conserved motif, we identified 8 different Hox gene fragments and Hox5 and Lox2 were found to be absent in early regeneration and during regular growth. Our data can promote further use of O. xiamen which can be used as an experimental model for resolving crucial problems of developmental biology in marine invertebrates.

Identification of Sex-Related Genes from the Three-Spot Swimming Crab Portunus sanguinolentus and Comparative Analysis with the Crucifix Crab Charybdis feriatus

Simple Summary

Crabs within the family Portunidae are important marine species in both aquaculture and fishery sectors. The current aquaculture status of most portunids still relies on wild-caught fisheries due to the lack of essential knowledge regarding their reproductive biology and underlying governing mechanism. In the present study, we compared the differentially expressed genes (DEGs) between the different sexes of Portunus sanguinolentus based on their gonadal transcriptome profiles and subsequently contrasted them with the gonadal DEGs of Charybdis feriatus, the other member of the family Portunidae. In total, 40,964 DEGs between the ovaries and testes of P. sanguinolentus were uncovered, with 27,578 up-regulated and 13,386 down-regulated in females. After comparison, C. feriatus has approximately 63.5% of genes in common with P. sanguinolentus, with 62.6% showing similar expression patterns. Interestingly, the DMRT gene was specifically expressed in male P. sanguinolentus, while its homologous gene—doublesex (DSX)—was specifically expressed in male C. feriatus. The DEGs obtained from the gonadal transcriptome of P. sanguinolentus are a beneficial resource for future genetic and genomic research in P. sanguinolentus and its close species. The transcriptomic comparison analysis might provide references for better understanding the sex determination and differentiation mechanisms among portunids.

Abstract

Crabs within the family Portunidae are important marine species in both aquaculture and fishery sectors. The current aquaculture status of most portunids, however, still relies on wild-caught fisheries due to the lack of essential knowledge regarding their reproductive biology and underlying governing mechanism. With the advancement of sequencing technology, transcriptome sequencing has been progressively used to understand various physiological processes, especially on non-model organisms. In the present study, we compared the differentially expressed genes (DEGs) between sexes of Portunus sanguinolentus based on their gonadal transcriptome profiles and subsequently contrasted them with the gonadal DEGs of Charybdis feriatus, the other member of Family Portunidae. In total, 40,964 DEGs between ovaries and testes were uncovered, with 27,578 up- and 13,386 down-regulated in females. Among those, some sex-related DEGs were identified, including a dmrt-like (DMRT) gene which was specifically expressed in males. C. feriatus has approximately 63.5% of genes common with P. sanguinolentus, with 62.6% showing similar expression patterns. Interestingly, the DMRT gene was specifically expressed in male P. sanguinolentus while its homologous gene—doublesex (DSX)—was specifically expressed in male C. feriatus. The DEGs obtained from the gonadal transcriptome of P. sanguinolentus are a beneficial resource for future genetic and genomic research in P. sanguinolentus and its close species. The transcriptomic comparison analysis might provide references for better understanding the sex determination and differentiation mechanisms among portunids.

Identification of differentially abundant mRNA transcripts and autocrine/paracrine factors in oocytes and follicle cells of mud crabs

The late vitellogenic stage of the mud crab is characterized by large and obvious follicle cells as well as an enlarged oocyte nucleus and a prominent germinal vesicle (GV). The aim of this study was evaluation of functions of oocytes and follicle cells during meiosis as well as at identifying associated ovarian autocrine/paracrine factors using comparative transcriptomics. The results from the KEGG pathway analysis indicated DNA replication, nucleotide excision repair, spliceosome and the ribosome pathways were highly associated with oocyte maturation across both transcriptomes. In addition, there was a larger abundance of mRNA transcripts for cell cycle-related genes in the oocyte, as well as cyclin A, cyclin B and CKS1B in the GV than at the time of germinal vesicle breakdown (GVBD). These findings indicate these cell cycle-related genes might be involved in GVBD induction. Results when there was localization of ligands and the respective receptors of VEGF, TGFβ propeptide and BMP9/10 indicated these proteins might be autocrine/paracrine factors. Results from functional analysis of VEGF, TGFβ propeptide and BMP9/10 in oocyte maturation using RNA interference revealed that these proteins might be involved in oocyte maturation by regulating cyclin abundance. This is the first study on the functions of VEGF in oocyte maturation in invertebrates.

The expression regulation of Cyclins and CDKs in ovary via miR-9c and miR-263a of Scylla paramamosain

Scylla paramamosain is an economically important cultured crab species in China. Cyclins and cyclin-dependent kinases (CDKs) play important roles in regulations of cell cycle and ovarian development. MiRNAs can negatively regulate gene expression at the post-transcriptional level through base-complementary pairing with the 3'-untranslated region (3-UTR) of the target gene. In this study, bioinformatics prediction showed that miR-9c and miR-263a identified from our group's gonad miRNAome of S. paramamosain may bind to the 3' UTR region of cyclin A, cyclin B, cyclin E, cyclin H, CDK1, and CDK2. Furthermore, the results of double luciferase reporter gene assay showed that the luciferase activities of HEK293T cells co-transfected with miR-9c mimics/miR-9c inhibitor and the 3'-UTR plasmid vectors of the five genes (cyclin A, cyclin B, cyclin H, CDK1, and CDK2) were significantly decreased/increased compared with those in the NC (negative control) and BC (blank control) groups. The results in miR-263a were similar to miR-9c, but all of the six genes could be regulated by miR-263a. In in vivo experiments, agomiR-9c (miR-9c enhancer) injection resulted in decreases of cyclin A and CDK1 expression level, and reverse effects were observed by injecting antagomiR-9c. AgomiR-263a decreased the expression of cyclin A, cyclin B, cyclin H, CDK1, and CDK2, but antagomiR-263a increased their expression. Both the in vitro and in vivo experiments confirmed functions of miR-9c and miR-263a in cell cycle progress of ovarian development by expression regulation of cyclin A, cyclin B, cyclin E, cyclin H, CDK1, and CDK2. The findings provide new insights into the reproductive regulation mechanism in mud crab and further enrich the knowledge of cell cycle and ovarian development regulation in invertebrates.

Stimulatory roles of epidermal growth factor receptor (EGFR) in ovarian development of mud crab Scylla paramamosain

The epidermal growth factor receptor (EGFR) is a pleiotropic glycoprotein which plays a role in regulating cell proliferation, migration and differentiation. However, to date little is known about its functions in crustaceans. In this study, we successfully identified SpEGFR from mud crab Scylla paramamosain. RT-PCR result showed that SpEGFR was widely distributed in all tested tissues and highly expressed in ovary. In situ hybridization revealed that SpEGFR mainly localized in oocyte perinuclear region with notably obvious signals. In vitro experiments showed that the expression of SpVgR and SpCyclin B in ovary explants from late vitellogenic stage crabs (summer) were significantly increased when treated with 1 nM human EGF (hEGF) for 1 h, while there was no obvious change towards SpEGFR. Interestingly, as for winter crab at the same vitellogenic stage, the expression of SpVgR and SpCyclin B in ovary explants did not show significant increase until treated with higher concentration of 10 nM hEGF and longer incubation time of 12 h. In addition, the hEGF-induced effect could be suppressed by pre-treated with EGFR inhibitor AG1478 and PD153035, respectively, which further indicated that EGF-EGFR pathway played a vital role in ovarian development in mud crab. In conclusion, SpEGFR might promote ovarian development by stimulating the expression of SpVgR and SpCyclin B under hEGF-induced treatment. The different physiological response to hEGF in the same vitellogenic stage crabs between summer and winter might be attributed to the changes in metabolism and physiological sensitivity.

Identification of MicroRNAs and Their Target Genes Associated with Ovarian Development in Black Tiger Shrimp (Penaeus monodon) Using High-Throughput Sequencing

Plenty of evidence showing that microRNAs (miRNAs) post-transcriptionally regulate gene expression and are involved in a wide range of biological processes. However, the roles of miRNAs in ovarian development process remain largely unknown in shrimp. In the present study, high-throughput sequencing of small RNAs was performed to find specific miRNAs that are involved in ovarian development process in Penaeus monodon. Two small RNA libraries were constructed from undeveloped (UNDEV group) and developed (DEV group) ovarian tissues in P. monodon. In total, 43 differentially expressed miRNAs were identified between the two groups (P ≤ 0.05, |log2 ratio| ≥1), and their expression profiles were validated by qRT-PCR. In order to further clarify the functional roles of these differentially expressed miRNAs during ovarian development process, target gene prediction was performed. In total, 4,102 target genes of 43 miRNAs were predicted, then clustered by the Kyoto Encyclopedia of Genes and Genomes (KEGG) database; only four specific pathways related to ovarian development were obtained (P < 0.05). Dual-luciferase reporter assays and integrated expression analysis were also conducted to further clarify the interaction between the miRNAs and their target mRNAs. This study provides important information about the function of miRNAs involved in ovarian developmental stages in P. monodon.

Molecular Characterization and Expression Profiles of Sp-uchl3 and Sp-uchl5 during Gonad Development of Scylla paramamosain

Ubiquitin C-terminal hydrolases (UCHLs) are a subset of deubiquitinating enzymes, and are involved in numerous physiological processes. However, the role of UCHLs during gonad development has not been studied in crustaceans. In this study, we have first cloned and analyzed expression profiling of Sp-uchl3 and Sp-uchl5 genes from mud crab Scylla paramamosain. The full-length cDNA of Sp-uchl3 is of 1804 bp. Its expression level in the ovary was significantly higher than in other tissues (p < 0.01), and during gonadal development, its expression in both O1 and O5 stages was significantly higher than in the other three stages of ovaries (p < 0.05), while in T3 it was higher than in the former two stages of testes (p < 0.05). Meanwhile, the full-length cDNA of Sp-UCHL5 is 1217 bp. The expression level in the ovary was significantly higher than in other tissues (p < 0.01). Its expression in ovaries was higher than in testes during gonadal development (p < 0.05). The expression level in the O5 stage was the highest, followed by the O3 stage in ovarian development, and with no significant difference in the testis development (p > 0.05). These results provide basic data showing the role of Sp-UCHL3 and Sp-UCHL5 in the gonad development of the crab.

Inhibitory Role of the Mud Crab Short Neuropeptide F in Vitellogenesis and Oocyte Maturation via Autocrine/Paracrine Signaling

Neuropeptides, in addition to their classical role in the nervous system, act on intraovarian factors to regulate reproductive functions in vertebrates. However, this function of neuropeptides has not been characterized in crustaceans. Short neuropeptide F (sNPF), a highly conserved invertebrate neuropeptide, has been reported to be involved in feeding, metabolism, and in differentiation processes including reproduction. Although sNPF and its receptor (sNPFR) have been detected in the ovary in different species, ovarian colocalization of sNPF/sNPFR has not been investigated. In this study, we identified Scylla paramamosain (mud crab) sNPF (Sp-sNPF) as an endogenous ligand for the S. paramamosain orphan G protein-coupled receptor NPY2R in mammalian cell line HEK293T. We designated this receptor as Sp-sNPFR. RNA in situ hybridization in pre-vitellogenic ovary and reverse transcription-PCR on isolated denuded oocytes and follicle layers showed that Sp-sNPF was exclusively localized to the follicle cells, whereas Sp-sNPFR was detected in both follicle cells and oocytes. We also found that Sp-sNPF partly suppressed spontaneous maturation of denuded oocytes and caused intracellular cAMP accumulation and Ca²⁺ mobilization. Moreover, injection of synthetic Sp-sNPF peptides inhibited the expression of vitellogenin and vitellogenin receptor genes in vivo. These combined results suggest for the first time that Sp-sNPF may have inhibitory functions in vitellogenesis and oocyte maturation possibly via the autocrine/paracrine pathway.

Transcriptome analysis of the copepod Eurytemora affinis upon exposure to endocrine disruptor pesticides: Focus on reproduction and development

Copepods-which include freshwater and marine species-represent the most abundant group of aquatic invertebrates. Among them, the calanoid copepod Eurytemora affinis is widely represented in the northern hemisphere estuaries and has become a species of interest in ecotoxicology. Like other non-target organisms, E. affinis may be exposed to a wide range of chemicals such as endocrine disruptors (EDs). This study investigated the gene expression variation in E. affinis after exposure to ED pesticides-chosen as model EDs-in order to (i) improve the knowledge on their effects in crustaceans, and (ii) highlight relevant transcripts for further development of potential biomarkers of ED exposure/effect. The study focused on the reproduction function in response to ED. Copepods were exposed to sublethal concentrations of pyriproxyfen (PXF) and chlordecone (CLD) separately. After 48h, males and females (400 individuals each) were sorted for RNA extraction. Their transcriptome was pyrosequenced using the Illumina(®) technology. Contigs were blasted and functionally annotated using Blast2GO(®). The differential expression analysis between ED- and acetone-exposed organisms was performed according to sexes and contaminants. Half of the 19,721 contigs provided by pyrosequencing were annotated, mostly (80%) from arthropod sequences. Overall, 2,566 different genes were differentially expressed after ED exposures in comparison with controls. As many genes were differentially expressed after PXF exposure as after CLD exposure. In contrast, more genes were differentially expressed in males than in females after both exposures. Ninety-seven genes overlapped in all conditions. Finally, 31 transcripts involved in reproduction, growth and development, and changed in both chemical exposures were selected as potential candidates for future development of biomarkers.

Role of linker histone H1c during the reprogramming of Chinese swamp buffalo (Bubalus Bubalis) embryos produced by somatic cell nuclear transfer

During reprogramming, there is exchange of histone H1c and the oocyte-specific linker histone, and H1c may play a critically important role in the reprogramming process of somatic cell nuclear transfer (SCNT). The aim of the present study was to investigate the role of the H1c gene in SCNT reprogramming in Chinese swamp buffalo (Bubalus bubalis) using RNA interference (RNAi). Chinese swamp buffalo H1c gene sequences were obtained and H1c-RNAi vectors were designed, synthesised and then transfected into a buffalo fetal skin fibroblast cell line. Expression of H1c was determined by real-time polymerase chain reaction to examine the efficiency of vector interference. These cells were then used as a nuclear donor for SCNT so as to observe the further development of SCNT embryos. Inhibition of H1c gene expression in donor cells significantly improved the developmental speed of embryos from the 1-cell to 8-cell stage. Furthermore, compared with the control group, inhibition of H1c gene expression significantly reduced the blastocyst formation rate. It is concluded that linker histone H1c is very important in SCNT reprogramming in Chinese swamp buffalo. Correct expression of the H1c gene plays a significant role in preimplantation embryonic development in B. bubalis.

Transcriptome analysis of the differences in gene expression between testis and ovary in green mud crab (Scylla paramamosain)

Background

The green mud crab (Scylla paramamosain) is the most prevalent crustacean on the southeast coast of China. The molecular regulatory mechanism of sex determination and gonadal differentiation in this species has received considerable attention in recent years because of the huge differences—both biological and economic—between male and female crabs. In this study, next-generation sequencing technology was used to develop deep-coverage transcriptomic sequencing data for the testis and ovary of S. paramamosain.

Results

A total of 365,116 reads (testis 171,962, ovary 193,154) with an average sequence length of 285 bp were produced from testis and ovary cDNA libraries. After filtering out contaminating reads, the clean reads were assembled, producing a total of 21,791 isotigs and leaving 22,814 reads as singlets. Using the BLASTX program, 3,471 unique sequences (2,275 isotigs and 1,196 singletons) were annotated with known protein sequences from the NCBI non-redundant (Nr) protein sequence database. The Gene Ontology and KEGG (Kyoto Encyclopedia of Genes and Genomes) analyses allowed the 224 unique sequences that were annotated with enzyme code (EC) numbers to be mapped into 174 KEGG pathways. After comparing the ovary and testis libraries, 4,021 gonad-differentially, 10,522 ovary-specifically, and 19,013 testis-specifically expressed genes were identified. Moreover, 33 ovary-specific, 14 testis-specific, and 34 gonad-differential transcripts were confirmed by semi-quantitative PCR and quantitative real-time PCR. In addition, 8,610 putative simple sequence repeats (SSRs) and 23,879 potential single nucleotide polymorphisms (SNPs) were identified.

Conclusion

This is the first large-scale RNA sequencing of S. paramamosain to be reported. We have identified many important functional genes and made a preliminary attempt to construct the regulatory network involved in the gonadal development of crustaceans. The annotated transcriptome data will provide fundamental support for future research into the reproduction biology of S. paramamosain. A large number of candidate SSRs and SNPs were detected, which could be used as genetic markers for population genetics and functional genomics in this species.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-585) contains supplementary material, which is available to authorized users.

Characterization and expression of cell division cycle 2 (Cdc2) mRNA and protein during ovarian development of the giant tiger shrimp Penaeus monodon

The meiotic maturation of oocytes is regulated by the maturation-promoting factor (MPF), a complex of Cdc2 (Cdk1) and Cyclin B. Here, the complete open reading frame (ORF) of Cdc2 in Penaeus monodon was characterized. PmCdc2 were 900bp in length corresponding to a polypeptide of 299 amino acids with the conserved Thr14, Tyr15 and Thr161 residues. Quantitative real-time PCR indicated that the expression level of PmCdc2 in wild intact broodstock was significantly increased in stages II (vitellogenic) and III (early cortical rod) ovaries relative to stage I (previtellogenic) ovaries and peaked in stage IV (mature) ovaries (P<0.05). The expression level of PmCdc2 in stages I-IV ovaries of eyestalk-ablated broodstock was greater than that of the same ovarian developmental stages in intact broodstock (P<0.05). Expression levels of PmCdc2 in ovaries of 18-month-old P. monodon upon 5-HT injection (50μg/g body weight) were significantly increased at 1hour post injection (hpi, P<0.05). Recombinant PmCdc2 protein and its polyclonal antibody were successfully produced. Western blot analysis revealed the expected 34kDa band (PmCdc2) along with a smaller band of 23kDa (ribosomal protein S3) in ovaries of juveniles and various ovarian stages of broodstock. Using phospho-Cdc2 (Thr161) polyclonal antibody, the positive signal of 34kDa was observed in all ovarian stages but the most intense signal was found in stage IV ovaries. Results in the present study indicated that PmCdc2 gene/protein plays an important role in the development and maturation of oocytes/ovaries in P. monodon.