Molecular characterization and expression profiles of cyclin A and cyclin B during ovarian development of the giant tiger shrimp Penaeus monodon

Proteome analysis provides insights into sex differences in Holothuria Scabra

Sex-determining mechanism is still ambiguous for sea cucumber Holothuria scabra which only manifests gonochorism in gonad. In this study, proteomic analysis was employed to delineate sex-related proteins and genes in gonads of H. scabra, subsequently validated through Quantitative real-time polymerase chain reaction (qRT-PCR). A total of 5,313 proteins were identified via proteome sequencing. Among these, 817 proteins exhibited expression in both the ovary and testis, with 445 proteins displaying up-regulation and 372 proteins showing down-regulation (ovary vs testis). Furthermore, 136 and 69 proteins were identified as ovary-specific and testis-specific Differentially Abundant Proteins (DAPs), respectively. And 9 DAP coding genes which play crucial role in ovary and testis were verified by qRT-PCR. Notably, 24 ovary-bias proteins enriched in ribosome pathway strongly indicated the crucial role of ribosome in ovary. This study serves to furnish novel evidence pertaining to sex differences in H. scabra.

Differentially expressed genes in the testes from early to mature development of banana shrimp (Fenneropenaeus merguiensis)

Banana shrimp (Fenneropenaeus merguiensis) is an economically important species in Thailand owing to the high value of globally exported frozen brine shrimps. However, the regulatory mechanisms governing spermatogenesis and testicular development in this species are poorly understood. High-throughput RNA sequencing was used to investigate the mechanisms and regulated genes involved in testis development using transcriptome profiling of juvenile and adult banana shrimp testes. Differentially expressed genes (DEGs) in these two libraries were identified and quantified to confirm gene expression. DEGs were found in 7,347 genes, with 4,465 upregulated and 2,882 downregulated. Some of these genes were designated as candidate genes, and six specific DEGs, including PRM1, SPATA20, Sry, SSRF, Sxl, and Tra-2c, were selected to confirm the reliability of the RNA-seq data using qPCR. Moreover, six non-DEGs were chosen based on testis-specific and regulatory genes that support a specific function in spermatogenesis and testis development in this species, including Dsx, Gfra2, IAG, Sox9, Sox13, and Sox14A. Furthermore, Sry, Sox14A, Sox14B and SPATA20 were identified in early stages (nauplius-postlarvae) of shrimp development to provide more information involving testes formation and development. The transcript data from this study could differentiate a group of genes required at the early and late stages of testis development and both sets of testis development. Therefore, this information would help in manipulating each stage of testicular development.

Molecular Markers of Ovarian Germ Cells of Banana Prawn (Fenneropenaeus merguiensis)

The banana prawn (Fenneropenaeus merguiensis) is a valuable prawn in the worldwide market. However, cultivation of this species is limited owing to the difficulty in culture management and limited knowledge of reproduction. Therefore, we studied the gene expression and molecular mechanisms involved in oogenesis for elucidating ovarian germ cell development in banana prawns. The tissue-specific distribution of certain genes identified from previous transcriptome data showed that FmCyclinB, FmNanos, and nuclear autoantigenic sperm protein (FmNASP) were only expressed in gonads. The in situ hybridization (ISH) of these three genes showed different expression patterns throughout oogenesis. FmCyclinB was highly expressed in pre-vitellogenic oocytes. FmNanos was expressed at almost the same level during oogenesis but showed the most expression in late pre-vitellogenic stages. Based on the highest expression of FmCyclinB and FmNanos in mid pre-vitellogenic and late pre-vitellogenic oocytes, respectively, we suggested that FmNanos may suppress FmCyclinB expression before initiation of vitellogenesis. Meanwhile, FmNASP expression was detected only in pre-vitellogenesis. Moreover, quantitative real-time polymerase chain reaction (qRT-PCR) analysis of FmNASP expression was supported by FmNASP ISH analysis based on high expression of FmNASP in sub-adult ovaries, which contain most of pre-vitellogenic oocytes. In this study, we found three reliable ovarian markers for banana prawns and also found a dynamic change of molecular mechanism during the sub-stage of pre-vitellogenesis. We determined the expression levels of these genes involved in oogenesis. Our findings provide information for further studies on banana prawn reproduction which may assist in their cultivation.

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.

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.

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

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

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.

Proteomic and deep sequencing analysis of extracellular vesicles isolated from adult male and female Schistosoma japonicum

Schistosomes are the causative agent of schistosomiasis, which affects more than 200 million people worldwide. Unlike other trematode parasites, schistosomes (along with the Didymozoidae) have evolved separate sexes. Pairing of males and females is a prerequisite for female sexual development and subsequent egg production. However, the mechanisms underlying these processes remain poorly understood. Extracellular vesicles (EVs) have been shown to play important roles in many biological processes. In the present study, we characterized EVs isolated from adult male and female Schistosoma japonicum. Proteomic analyses of the isolated EVs revealed that some proteins are significantly enriched in male or female EVs. RNA-sequencing analysis of a small RNA population associated with EVs identified 18 miRNAs enriched in male and female S. japonicum EVs. Among these, miR-750 was specifically enriched in female EVs. Additionally, the inhibition of miR-750 by a miRNA inhibitor led to decreased egg production in female schistosomes cultured in vitro. Collectively, our results suggest that miR-750 within female EV cargo may be involved in regulating ovary development and egg production in S. japonicum females.

Schistosomiasis is a neglected tropical disease caused by the genus Schistosoma and affects more than 200 million people worldwide. Previously, we and other groups found that Schistosoma japonicum can secrete extracellular vesicles (EVs) that are taken up by mammalian cells. Here, we characterized EVs isolated from adult male and female S. japonicum and found that some proteins and microRNAs (miRNAs) were significantly enriched in male or female EVs. More importantly, the inhibition of miR-750, which is specifically enriched in female EVs, resulted in decreased egg production. Overall, our study suggests that female EV miRNA cargo may play important roles in regulating female ovary development and egg production during male-female pairing in S. japonicum.

Identification of a Growth-Associated Single Nucleotide Polymorphism (SNP) in Cyclin C of the Giant Tiger Shrimp Penaeus monodon

The full-length cDNA of cyclin C of the giant tiger shrimp Penaeus monodon (PmCyC) was isolated by RACE-PCR. It was 1443 bp in length containing an open reading frame (ORF) of 804 bp and 267 deduced amino acids. Tissue distribution analysis indicated that PmCyC was more abundantly expressed in ovaries and testes than other tissues of female and male juveniles (P < 0.05). A pair of primers was designed, and an amplification product of 403 bp containing an intron of 123 bp was obtained. Polymorphism of amplified PmCyC gene segments of the 5th (3-month-old G5, N = 30) and 7th (5-month-old G7, N = 18) generations of domesticated juveniles was analyzed. Four conserved SNPs (T>C₁₃₄, T>C₁₈₈, G>A₃₇₉, and T>C₃₈₂) were found within the examined sequences. A TaqMan genotyping assay was developed for detection of a T>C₁₃₄ SNP. Association analysis indicated that this SNP displayed significant association with body weight (P < 4.2e-10) and total length (P < 2e-09) of the examined G7 P. monodon (N = 419) with an allele substitution effect of 5.02 ± 0.78 g and 1.41 ± 0.19 cm, respectively. Juveniles with C/C₁₃₄ (22.80 ± 2.51 g and 12.97 ± 0.53 cm, N = 19) and T/C₁₃₄ (20.41 ± 0.93 g and 12.77 ± 0.21 cm, N = 129) genotypes exhibited a significantly greater average body weight and total length than those with a T/T₁₃₄ genotype (14.72 ± 0.53 g and 11.37 ± 0.13 cm, N = 271) (P < 0.05).

Exploring potential proteins associated with cyclin B 3'UTR in Procambarus clarkii oocytes

Cyclin B is a ubiquitous regulatory molecule and has been implicated in mitosis and meiosis in oocytes. Phenomenon that differ in the length of cyclin B 3'UTR in crustacean has attracted much attention, although molecular details are poorly understood. The study of 3'UTR-interacting proteins could yield much information in translational regulation and the mRNA localization process. Previous studies on crayfish suggested that the 3'UTR (1300 bp) probably contains the potential regulatory sequence/motifs such as CPEs and K-box et al. In present study, using pull-down assay coupled with mass spectrometry approach allowing us to explore the potential proteins associated with the 3'UTR. We finally identified four candidate proteins including Hspg 2, Vtg, eef1a and Tuba1a, which annotated as significant roles involved in cell differentiation, lipid transporter activity, and meiotic cell cycle process. The preliminary results will contribute to the advance in understanding the translational activation of cyclin B in oocyte maturation regulation in crustacean.

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.

Genomic structure, expression, and functional characterization of checkpoint kinase 1 from Penaeus monodon

Chk1 is a cell-cycle regulator. Chk1 has been identified in organisms ranging from yeast to humans, but few researchers have studied Chk1 in shrimps. We cloned Chk1 from the black tiger shrimp (Penaeus monodon). The full-length cDNA sequence of PmChk1 had 3,334 base pairs (bp), with an open reading frame of 1,455 bp. The complete genomic sequence of PmChk1 (11,081 bp) contained 10 exons separated by nine introns. qRT-PCR showed that PmChk1 was highly expressed in the ovaries and gills of P. monodon. The lowest PmChk1 expression was noted in stage III of ovarian development in P. monodon. PmChk1 expression decreased significantly after injection of 5-hydroxytryptamine and eyestalk ablation in P. monodon ovaries. RNA interference experiments were undertaken to examine the expression of PmChk1, PmCDC2, and PmCyclin B. PmChk1 knockdown in the ovaries and hepatopancreas by dsRNA-Chk1 was successful. The localization and level of PmChk1 expression in the hepatopancreas was studied using in situ hybridization, which showed that data were in accordance with those of qRT-PCR. The Gonadosomatic Index of P. monodon after dsRNA-Chk1 injection was significantly higher than that after injection of dsRNA-GFP or phosphate-buffered saline. These data suggest that PmChk1 may have important roles in the ovarian maturation of P. monodon.

De novo transcriptome analysis and differentially expressed genes in the ovary and testis of the Japanese mantis shrimp Oratosquilla oratoria by RNA-Seq

The mantis shrimp Oratosquilla oratoria is a widely distributed, commercially important crustacean species. Although its conservation and the development of successful artificial breeding technologies have recently received considerable attention, there are currently no available data regarding the molecular mechanisms in controlling reproduction. In this study, we performed transcriptome sequencing of the testis, ovary, female and male eyestalks and the androgenic gland of O. oratoria, and compared the expression pattern of transcripts from the testis and ovary libraries to identify genes involved in gonadal development. A total of 147,130,937 clean reads were retrieved after removing the adapters in reads and filtering out low-quality data. All the reads were assembled into 94,990 unigenes (23,133 in testis and ovary) with an average length of 783 base pairs (bp) and N50 of 1502 bp. A search of all-unigenes against COG, GO, KEGG, KOG, Pfam, Swiss-Prot and Nr databases resulted in a total of 19,404 annotated unigenes. Comparison of the sequences in the ovary and testis libraries revealed that 1188 unigenes were up-regulated in the ovary and 2732 were up-regulated in the testis. Twenty ovary-up-regulated and 21 testis-up-regulated unigenes were confirmed by quantitative real-time PCR. Additionally, 13,437 simple sequence repeats (SSRs) and 275,799 putative single nucleotide polymorphisms (SNPs) were identified. The important functional genes and pathways identified here provide a valuable dataset for understanding the molecular mechanisms controlling gonad development in O. oratoria, and the numerous (13,437 SSRs and 275,799 SNPs) molecular markers obtained here will provide fundamental basis for functional genomic and population genetic studies of O. oratoria.

Molecular cloning, characterization and expression analysis of a novel PDRG1 gene from black tiger shrimp (Penaeus monodon)

P53 And DNA Damage-Regulated Gene 1 (PDRG1) is a novel gene which plays an important role in chaperone-mediated protein folding. In the present study, the full-length complementary DNA (cDNA) sequence of the PDRG1 gene from Penaeus monodon (PmPDRG1) was cloned by the rapid amplification of cDNA ends (RACE) method. The cDNA of PmPDRG1 spans 1,613 bp, interrupted by only one short intron, and encodes a protein of 136 amino acids with calculated molecular weight of 15.49 kDa. The temporal expression profile of PmPDRG1 in different tissues and in different developmental stages of the ovary was investigated by real-time quantitative PCR (RT-qPCR). An RNA interference (RNAi) experiment was performed to study the relationship between P. monodon p53 (Pmp53) and PmPDRG1, and the results showed that the relative expression level of PmPDRG1 mRNA was notably up-regulated from 12 h to 96 h after Pmp53 was silenced both in ovary and hepatopancreas. To further explore the role of PmPDRG1 in ovarian development, dopamine (DA) and 5-hydroxytryptamine (5-HT)-injected shrimps were analyzed by RT-qPCR, indicating that PmPDRG1 may be involved in the regulation of ovarian development of P. monodon.

Optimizing Hybrid de Novo Transcriptome Assembly and Extending Genomic Resources for Giant Freshwater Prawns (Macrobrachium rosenbergii): The Identification of Genes and Markers Associated with Reproduction

The giant freshwater prawn, Macrobrachium rosenbergii, a sexually dimorphic decapod crustacean is currently the world’s most economically important cultured freshwater crustacean species. Despite its economic importance, there is currently a lack of genomic resources available for this species, and this has limited exploration of the molecular mechanisms that control the M. rosenbergii sex-differentiation system more widely in freshwater prawns. Here, we present the first hybrid transcriptome from M. rosenbergii applying RNA-Seq technologies directed at identifying genes that have potential functional roles in reproductive-related traits. A total of 13,733,210 combined raw reads (1720 Mbp) were obtained from Ion-Torrent PGM and 454 FLX. Bioinformatic analyses based on three state-of-the-art assemblers, the CLC Genomic Workbench, Trans-ABySS, and Trinity, that use single and multiple k-mer methods respectively, were used to analyse the data. The influence of multiple k-mers on assembly performance was assessed to gain insight into transcriptome assembly from short reads. After optimisation, de novo assembly resulted in 44,407 contigs with a mean length of 437 bp, and the assembled transcripts were further functionally annotated to detect single nucleotide polymorphisms and simple sequence repeat motifs. Gene expression analysis was also used to compare expression patterns from ovary and testis tissue libraries to identify genes with potential roles in reproduction and sex differentiation. The large transcript set assembled here represents the most comprehensive set of transcriptomic resources ever developed for reproduction traits in M. rosenbergii, and the large number of genetic markers predicted should constitute an invaluable resource for future genetic research studies on M. rosenbergii and can be applied more widely on other freshwater prawn species in the genus Macrobrachium.

Gonadal transcriptomic analysis and differentially expressed genes in the testis and ovary of the Pacific white shrimp (Litopenaeus vannamei)

Background

The Pacific white shrimp (Litopenaeus vannamei) is the world’s most prevalent cultured crustacean species. However, the supply of high-quality broodstocks is limited and baseline information related to its reproductive activity and molecular issues related to gonad development are scarce. In this study, we performed transcriptome sequencing on the gonads of adult male and female L. vannamei to identify sex-related genes.

Results

A total of 25.16 gigabases (Gb) of sequences were generated from four L. vannamei gonadal tissue libraries. After quality control, 24.11 Gb of clean reads were selected from the gonadal libraries. De-novo assembly of all the clean reads generated a total of 65,218 unigenes with a mean size of 1021 bp and a N50 of 2000 bp. A search of all-unigene against Nr, SwissProt, KEGG, COG and NT databases resulted in 26,482, 23,062, 20,659, 11,935 and 14,626 annotations, respectively, providing a total of 30,304 annotated unigenes. Among annotated unigenes, 12,320 unigenes were assigned to gene ontology categories and 20,659 unigenes were mapped to 258 KEGG pathways. By comparing the ovary and testis libraries, 19,279 testicular up-regulated and 3,529 ovarian up-regulated unigenes were identified. Enrichment analysis of differentially expressed unigenes resulted in 1060 significantly enriched GO terms and 34 significantly enriched KEGG pathways. Nine ovary-specific, 6 testis-specific, 45 testicular up-regulated and 39 ovarian up-regulated unigenes were then confirmed by semi-quantitative PCR and quantitative real-time PCR. In addition, using all-unigenes as a reference, a total of 13,233 simple sequence repeats (SSRs) were identified in 10,411 unigene sequences.

Conclusions

The present study depicts the first large-scale RNA sequencing of shrimp gonads. We have identified many important sex-related functional genes, GO terms and pathways, all of which will facilitate future research into the reproductive biology of shrimp. We expect that the SSRs detected in this study can then be used as genetic markers for germplasm evaluation of breeding and imported populations.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2219-4) contains supplementary material, which is available to authorized users.

Transcriptional changes during Daphnia pulex development indicate that the maturation decision resembles a rate more than a threshold

Maturation is a critical developmental process, and the age and size at which it occurs have important fitness consequences. Although maturation is remarkably variable, certain mechanisms, including a minimum size or state threshold, are proposed to underlie the process across a broad diversity of taxa. Recent evidence suggests that thresholds may themselves be developmentally plastic, and in the crustacean Daphnia pulex it is unclear whether maturation follows a threshold or is a gradual process more akin to a rate. Changes in gene expression across four instars before and during maturation were compared in a cDNA microarray experiment. Developmental stage was treated statistically both as a discontinuous and as a continuous variable, to determine whether genes showed gradual or discrete changes in expression. The continuous analysis identified a greater number of genes with significant differential expression (45) than the discontinuous analysis (11). The majority of genes, including those coding for histones, factors relating to transcription and cell cycle processes, and a putative developmental hormone showed continuous increases or decreases in expression from the first to the fourth instars that were studied, suggestive of a prolonged and gradual maturation process. Three genes coding for a fused vitellogenin/superoxide dismutase showed increases in expression following the second instar and coincided with the posited maturation threshold, but even their expression increased in a continuous fashion.

Transcriptome analysis of red swamp crawfish Procambarus clarkii reveals genes involved in gonadal development

Background

The red swamp crawfish, Procambarus clarkii, has become one of the most economically important cultured species in China. Currently, little is known about the gonadal development of this species. Isolation and characterization of genes are an initial step towards understanding gonadal development of P. clarkii.

Results

Using the 454 pyrosequencing technology, we obtained a total of 1,134,993 high quality sequence reads from the crawfish testis and ovary libraries. We aimed to identify different genes with a potential role in gonad development. The assembly formed into 22,652 isotigs, distributed by GO analysis across 55 categories in the three ontologies, ‘molecular function’, ‘cellular component’, and ‘biological processes’. Comparative transcript analysis showed that 1,720 isotigs in the ovary were up-regulated and 2138 isotigs were down-regulated. Several gonad development related genes, such as vitellogenin, cyclin B, cyclin-dependent kinases 2, Dmc1 and ubiquitin were identified. Quantitative real-time PCR verified the expression profiles of 14 differentially expressed genes, and confirmed the reliability of the 454 pyrosequencing.

Conclusions

Our findings provide an archive for future research on gonadal development at a molecular level in P. clarkii and other crustacean. This data will be helpful to develop new ideas for artificial regulation of the reproductive process in crawfish aquaculture.

Proteomic analysis of ovarian proteins and characterization of thymosin-β and RAC-GTPase activating protein 1 of the giant tiger shrimp Penaeus monodon

Cellular proteomics of total proteins in ovaries of domesticated and wild giant tiger shrimp (Penaeus monodon) were examined using GeLC-MS/MS. In total, 1638 proteins matched those previously deposited in databases and 1253 (76.50%) of these significantly matched known proteins. Several reproduction-related proteins (e.g. Cdc2, Cyclin B, Cdc25, 14-3-3, thymosin-β and Rac-GTPase activating protein 1) were identified. In addition, the full-length cDNA of P. monodon thymosin-β (PmTmsb; 1084 bp with an ORF of 387 bp and 128 deduced aa) and Rac-GTPase activating protein 1 (PmRacgap1; an ORF of 1881 bp and 626 deduced aa) were further characterized. PmTmsb was constitutively expressed in all tissues. In contrast, PmRacgap1 was more abundantly expressed in gonads than in several non-reproductive tissues (e.g. subcuticular epithelium, hepatopancreas, intestine, pleopods, stomach and thoracic ganglion). The expression levels of PmTmsb and PmRacgap1 in ovaries of wild adult broodstock were significantly greater than those in ovaries of juveniles (P<0.05). However, their expression levels did not vary significantly during ovarian development stages in intact broodstock. However, eyestalk ablation resulted in a significant reduction in PmTmsb expression at stages I and III ovaries (P<0.05), although it did not affect PmRacgap1 transcription significantly at these stages. On the other hand, use of polyclonal antibodies derived from recombinant PmTmsb and PmRacgap1 revealed that levels of both proteins decreased at the late stage (IV) of ovarian development. Our results suggested that PmTmsb and PmRacgap1 may act as negative effectors during ovarian development in P. monodon.

Comparative proteomic profiling during ovarian development of the shrimp Metapenaeus ensis

Two-dimensional electrophoresis and mass spectrometry were used to identify proteins that are differentially expressed during ovarian maturation in Metapenaeus ensis. 87 spots with consistently significant quantitative differences (≥ 1.5-fold for vol%) among stage I, III and V ovaries were chosen for MS/MS analysis. 45 spots were significantly matched to known proteins in the database (Mascot score >40). Half of them were down-regulated, in contrast to 9 out of 45 proteins that were up-regulated as ovarian maturation proceeded. Functionally, these identified proteins could be classified into five major groups, including cytoskeleton (11 %), metabolism (18 %), signal transduction (32 %), gene expression (14 %) and immune response (7 %). Among the differentially expressed reproduction-related proteins, the mRNA expression level of cellular retinoic acid/retinol binding protein in M. ensis (MeCRABP) during ovarian maturation was further characterized by quantitative real-time PCR. It was down-regulated during ovarian maturation. In situ hybridization further revealed that MeCRABP transcript was localized in ooplasm of previtellogenic oocytes but not in vitellogenic oocytes. These results demonstrate the application of proteomic analysis for identification of proteins involved in shrimp ovarian maturation and they provide new insights into ovarian development.

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.

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

The maturation promoting factor (MPF) is a key regulator of controlling G2/M phase transition in the meiotic maturation of oocyte and spermatocyte in animals, which is a complex of CDC2 (CDK1) and cyclin B. To better understand the molecular mechanism of oocyte and spermatocyte maturation in mud crab (Scylla paramamosain), the full length cDNA of cdc2 (Sp-cdc2) and cyclin B (Sp-cyclin B) were cloned and characterized. The full length cDNA of Sp-cdc2 gene is of 1593 bp encoding a protein of 299 amino acids. Real-time quantitative PCR analysis revealed that the expression level of Sp-cdc2 in the ovary was higher than in other tissues (P<0.01); and its expression level was not significantly different in different stages of ovary development (P>0.05), meanwhile there was higher expression in T3 stage than in T1 and T2 stages (P<0.05). The full length cDNA of Sp-cyclin B is 1492 bp encoding a protein of 391 amino acids. The real-time PCR results showed that its expression level in the ovary was the highest in all examined tissues (P<0.01), and the gonad expression level in O5 stage was significantly higher than in previous 4 stages and the testis (P<0.05), and was also significantly higher in T2 stage than in T1 stage (P<0.05). In situ hybridization analysis showed that the expressions of Sp-cdc2 and Sp-cyclin B transcripts were presented in similar distribution patterns in different developing stages of ovary and testis. The positive signals of Sp-cdc2 and Sp-cyclin B mRNA were detected in the oocytoplasm of oogonia and pre-vitellogenic and primary vitellogenic oocytes, while these two genes had higher expression level in the spermatid and secondary spermatocyte following primary spermatocyte. These results suggested that Sp-cdc2 and Sp-cyclin B may play essential roles in the oogenesis and spermatogenesis of the crab.

Molecular cloning and expression analysis of the Mitogen-activating protein kinase 1 (MAPK1) gene and protein during ovarian development of the giant tiger shrimp Penaeus monodon

Isolation and characterization of genes and/or proteins differentially expressed in ovaries are necessary for understanding ovarian development in the giant tiger shrimp (Penaeus monodon). In this study, the full-length cDNA of P. monodon mitogen-activating protein kinase 1 (PmMAPK1) was characterized. PmMAPK1 was 1,398 bp in length containing an open reading frame of 1,098 bp that corresponded to a polypeptide of 365 amino acids. PmMAPK1 was more abundantly expressed in ovaries than in testes of P. monodon. Quantitative real-time PCR revealed differential expression levels of PmMAPK1 mRNA during ovarian development of intact broodstock, where it peaked in early cortical rod (stage III) ovaries (P < 0.05) and slightly decreased afterwards (P > 0.05). Likewise, the expression level of PmMAPK1 in early cortical rod and mature (IV) ovaries was significantly greater than that in previtellogenic (I) and vitellogenic (II) ovaries of eyestalk-ablated broodstock (P < 0.05). The PmMAPK1 transcript was localized in ooplasm of previtellogenic oocytes. In intact broodstock, the expression of the PmMAPK1 protein was clearly increased from previtellogenic ovaries in subsequent stages of ovarian development (P < 0.05). In contrast, the level of ovarian PmMAPK1 protein was comparable during oogenesis in eyestalk-ablated broodstock (P > 0.05). The PmMAPK1 protein was localized in ooplasm of previtellogenic and vitellogenic oocytes. It was also detected around the nuclear membrane of early cortical rod oocytes in both intact and eyestalk-ablated broodstock. Results indicated that PmMAPK1 gene products seem to play functional roles in the development and maturation of oocytes/ovaries in P. monodon.

A cDNA microarray, UniShrimpChip, for identification of genes relevant to testicular development in the black tiger shrimp (Penaeus monodon)

Background

Poor reproductive maturation in captive male broodstock of the black tiger shrimp (Penaeus monodon) is one of the serious problems to the farming industries. Without genome sequence, EST libraries of P. monodon were previously constructed to identify transcripts with important biological functions. In this study, a new version of cDNA microarray, UniShrimpChip, was constructed from the Peneaus monodon EST libraries of 12 tissues, containing 5,568 non-redundant cDNA clones from 10,536 unique cDNA in the P. monodon EST database. UniShrimpChip was used to study testicular development by comparing gene expression levels of wild brooders from the West and East coasts of Thailand and domesticated brooders with different ages (10-, 14-, 18-month-old).

Results

The overall gene expression patterns from the microarray experiments revealed distinct transcriptomic patterns between the wild and domesticated groups. Moreover, differentially expressed genes from the microarray comparisons were identified, and the expression patterns of eight selected transcripts were subsequently confirmed by reverse-transcriptase quantitative PCR (RT-qPCR). Among these, expression levels of six subunits (CSN2, 4, 5, 6, 7a, and 8) of the COP9 signalosome (CSN) gene family in wild and different ages of domesticated brooders were examined by RT-qPCR. Among the six subunits, CSN5 and CSN6 were most highly expressed in wild brooders and least expressed in the 18-month-old domesticated group; therefore, their full-length cDNA sequences were characterized.

Conclusions

This study is the first report to employ cDNA microarray to study testicular development in the black tiger shrimp. We show that there are obvious differences between the wild and domesticated shrimp at the transcriptomic level. Furthermore, our study is the first to investigate the feasibility that the CSN gene family might have involved in reproduction and development of this economically important species.

Isolation and characterization of genes functionally involved in ovarian development of the giant tiger shrimp Penaeus monodon by suppression subtractive hybridization (SSH)

Suppression subtractive hybridization (SSH) libraries between cDNA in stages I (previtellogenic) and III (cortical rod) ovaries of the giant tiger shrimp (Penaeus monodon) were established. In all, 452 ESTs were unidirectionally sequenced. Sequence assembly generated 28 contigs and 201 singletons, 109 of which (48.0%) corresponding to known sequences previously deposited in GenBank. Several reproduction-related transcripts were identified. The full-length cDNA of anaphase promoting complex subunit 11 (PmAPC11; 600 bp with an ORF of 255 bp corresponding to a polypeptide of 84 amino acids) and selenoprotein Mprecursor (PmSePM; 904 bp with an ORF of 396 bp corresponding to a polypeptide of 131 amino acids) were characterized and reported for the first time in penaeid shrimp. Semiquantitative RT-PCR revealed that the expression levels of PmSePM and keratinocyte-associated protein 2 significantly diminished throughout ovarian development, whereas Ser/Thrcheckpoint kinase 1 (Chk1), DNA replication licensing factor mcm2 and egalitarian were down-regulated in mature ovaries of wild P. monodon (p < 0.05). Accordingly, the expression profiles of PmSePM and keratinocyte-associated protein 2 could be used as biomarkers for evaluating the degree of reproductive maturation in domesticated P. monodon.

Molecular cloning and expression of progestin membrane receptor component 1 (Pgmrc1) of the giant tiger shrimp Penaeus monodon

Knowledge on molecular mechanisms of steroid hormonal induction on oocyte development may lead to the possible ways to effectively induce ovarian maturation in shrimp. In this study, progestin membrane receptor component 1 (Pgmrc1) of the giant tiger shrimp (Penaeus monodon) initially identified by EST analysis was further characterized. The full-length cDNA of Pgmrc1 was 2015bp in length containing an ORF of 573bp corresponding to a polypeptide of 190 amino acids. Northern blot analysis revealed a single form of Pgmrc1 in ovaries of P. monodon. Quantitative real-time PCR indicated that the expression level of Pgmrc1 mRNA in ovaries of both intact and eyestalk-ablated broodstock was greater than that of juveniles (P<0.05). Pgmrc1 was up-regulated in mature (stage IV) ovaries of intact broodstock (P<0.05). Unilateral eyestalk ablation resulted in an earlier up-regulation of Pgmrc1 since the vitellogenic (II) ovarian stage. Moreover, the expression level of Pgmrc1 in vitellogenic, early cortical rod and mature (II-IV) ovaries of eyestalk-ablated broodstock was greater than that of the same ovarian stages in intact broodstock (P<0.05). Pgmrc1 mRNA was clearly localized in the cytoplasm of follicular cells, previtellogenic and early vitellogenic oocytes. Immunohistochemistry revealed the positive signals of the Pgmrc1 protein in the follicular layers and cell membrane of follicular cells and various stages of oocytes. Taken the information together, Pgmrc1 gene products seem to play the important role on ovarian development and may be used as the bioindicator for monitoring progression of oocyte maturation of P. monodon.

EST analysis and identification of gonad-related genes from the normalized cDNA library of large yellow croaker, Larimichthys crocea

On grounds of the especially limited numbers of identified gonad-specific or gonad-related genes of large yellow croaker Larimichthys crocea which may represent a major obstacle for the study of gonad development and sex differentiation, we initiated a sequencing program of Expressed Sequence Tags (ESTs) in large yellow croaker. In this study, we firstly constructed a normalized gonad cDNA library using the combination of SMART technique and DSN treatment. The titer of amplified cDNA library was 4.8x10(11) and the percentage of unique cDNA sequences of the library was 82.49%. 2916 unique cDNAs were clustered from the 3535 high quality ESTs. Among the 1785 ESTs which had significant homology with known genes in the NCBI database, about 64 significant gonad-related genes were found, accounting for 3.59% of the total unique cDNAs. Specifically, the testis-specific LRR gene and testis-specific chromodomain Y-like protein gene were identified from fish for the first time. Six gonad-related microsatellite-containing ESTs were identified from the 129 ESTs containing 149 microsatellites. Expression patterns of 10 of these gonad-related gene homologues in ovaries and testes were examined by qRT-PCR. The results will be powerful resources for our further investigation to establish the molecular mechanisms of gonad development and sex differentiation in large yellow croaker.