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

Transcriptome profiling of banana shrimp (Fenneropenaeus merguiensis) ovaries and testes: Insights into FoxL2

The banana shrimp is found in the Pacific and Indian Oceans. Female shrimp are preferred for consumption because they are larger than males. Understanding the mechanism of sex differentiation is important for developing techniques to increase the number of female shrimp for economic benefits. This study investigates the reproductive development of F. merguiensis using transcriptome analysis. Sxl2, dsx, AGH, FEM-1, and Nrg-X2 were classified as essential genes for testes development during the juvenile stage. Several genes were required for both juvenile and adult male development. Additionally, the expression of several genes was shown to be required for juvenile and adult ovarian development, including SOP1, SOP2, Ptgr1, EST, Vgr, Vmol1, and TR-beta A. Interestingly, high levels of FoxL2 expression were observed in the testes, in contrast to previous studies in humans and other mammals. The binding of FoxL2 to the Vtg promoter was demonstrated in silico with the highest relative binding score (RS = 0.89) using the JASPAR program. Knock-down of the FoxL2 gene with dsRNA significantly suppressed FoxL2 at 2, 4, and 6 d. As a result, Vtg expression increased when compared with the control at 2, 4, and 6 d, indicating that FoxL2 plays an important role in Vtg expression in the ovary. Our findings highlight the role of FoxL2 in banana shrimp reproduction and provide valuable information on the genes associated with the F. merguiensis reproductive system.

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.

Comparative Transcriptome Analysis of Differentially Expressed Genes in the Testis and Ovary of Sea Urchin (Strongylocentrotus intermedius)

The quality of sea urchin gonad is important to consumers with high standards for nutrition and taste. However, few studies have been conductedon the molecular mechanisms that determine the quality of male and female sea urchins. In this study, our goal was to understand the differences and characteristics of gonad quality between sea urchin (Strongylocentrotus intermedius) males and females. The transcriptomes of males and females were obtained, with totals of 43,797,146 and 56,222,782 raw reads, respectively, comprising 128,979 transcripts and 85,745 unigenes. After comparative transcriptome analysis, a total of 6736 differentially expressed genes (DEGs) between the males and females were identified, of which 2950 genes were up-regulated and 3786 genes were down-regulated in males. We compared the expression of twelve DEGs with significant differences their expression levels and functional annotations to confirm the reliability of the RNA-Seq data. Five DEGs related to gonadal quality were found through enrichment analysis of KEGG pathways: 17β-HSD8, PGDH, FAXDC2, C4MO, and PNPLA7. Our study analyzes genes related to the taste and flavor of sea urchin gonads among the sexes and provides reference sequences and fundamental information concerning the nutrition and taste of S. intermedius gonads.

Metabolomic Alterations in the Digestive System of the Mantis Shrimp Oratosquilla oratoria Following Short-Term Exposure to Cadmium

Mantis shrimp Oratosquilla oratoria is an economically critical aquatic species along the coast of China but strongly accumulates marine pollutant cadmium (Cd) in its digestive system. It is necessary to characterize the toxicity of Cd in the digestive system of mantis shrimp. The metabolic process is an essential target of Cd toxicity response. In this work, we used ultra-performance liquid chromatography coupled with time-of-flight mass spectrometry (UPLC-TOF-MS) for untargeted metabolomics to characterize the metabolic changes in the digestive system of O. oratoria, exposed to 0.05 mg/L for 96 h. The aim of this study was to further investigate the effect of O. oratoria on Cd response to toxicity and develop biomarkers. Metabolomics analysis showed the alteration of metabolism in the digestive system of mantis shrimp under Cd stress. A total of 91 metabolites were differentially expressed and their main functions were classified into amino acids, phospholipids, and fatty acid esters. The enrichment results of differential metabolite functional pathways showed that biological processes such as amino acid metabolism, transmembrane transport, energy metabolism, and signal transduction are significantly affected. Based on the above results, the Cd-induced oxidative stress and energy metabolism disorders were characterized by the differential expression of amino acids and ADP in mantis shrimp, while the interference of transmembrane transport and signal transduction was due to the differential expression of phospholipids. Overall, this work initially discussed the toxicological response of Cd stress to O. oratoria from the metabolic level and provided new insights into the mechanism.

Analyzing the gonadal transcriptome of the frog Hoplobatrachus rugulosus to identify genes involved in sex development

BACKGROUND

The tiger frog (Hoplobatrachus rugulosus) is listed as a national Class II protected species in China. In the context of global warming, the sex ratio of amphibians will be affected, and the development of the population will be limited. Therefore, considering the potential for a decrease in the number of amphibians, studying sex evolution and molecular regulation of gonadal development in H. rugulosus, phenomenon that are currently unclear, is of great significance.

RESULTS

Here, H. rugulosus was used to explore the mechanisms regulating gonadal development in amphibians. Illumina HiSeq 3000 was used to sequence the gonadal transcriptome of male and female H. rugulosus at two growth stages to identify genes related to gonadal development and analyze expression differences in the gonads. This analysis indicated that cyp17α, hsd3β, hsd11β1, cyp19α, and hsd17β12 perform vital functions in sex development in amphibians. Specifically, the expression of cyp3α, cyp17α, hsd3β, hsd11β1, sox2, sox9, sox30, soat, cyp19α, hsd17β12, and hspα1s was correlated with gonadal development and differentiation in H. rugulosus, as determined using the quantitative reverse transcriptase-polymerase chain reaction.

CONCLUSION

Significant differences were found in the gonadal gene expression levels in H. rugulosus of both sexes, and we identified a steroid hormone synthesis pathway in this species and analyzed related gene expression, but the changes during sex differentiation were still unclear. To our knowledge, this report presents the first analysis of the H. rugulosus gonadal transcriptome and lays the foundation for future research.

Characterization of the foxl2 gene involved in the vtg expression in mud crab (Scylla paramamosain)

Forkhead box protein L2 (Foxl2) is involved in multiple physiological processes, such as ovarian development, granulosa cell differentiation, ovarian follicle development, and oocyte growth. In this study, a Spfoxl2 gene encoded 530 amino acid protein with characteristic forkhead (FH) domain was identified from transcriptome data of mud crab Scylla paramamosain and validated the accuracy by PCR technology. Meanwhile, the orthologues of the Spfoxl2 gene in other 14 crustacean species were identified with the same method. Further multiple sequence alignment analysis revealed the Foxl2 was highly conserved, especially in the FH domain, even completely identical in several species. Besides, the semi-quantitative PCR (Sq-PCR) result showed Spfoxl2 gene was mainly expressed in the gonad (testis and ovary). Further quantitative real-time PCR (qRT-PCR) result demonstrated its expression level in the testis was significantly higher than that in the ovary (p < 0.01). In addition, the qRT-PCR result showed that in zoea V, megalopa, and larval I, the expression level of Spfoxl2 in megalopa is the highest. In addition, a putative Foxl2 binding site was identified on the promoter region of Spvtg, and knockdown of Spfoxl2 mediated by RNAi technology increased the expression of Spvtg in the ovary, suggesting Spfoxl2 might be the upstream negative regulator of Spvtg. Overall, this study provided new insights into the role of Spfoxl2 in ovary development through regulating Spvtg expression in S. paramamosain.

Genome-Wide DNA Methylation and Transcriptome Analyses Reveal Epigenetic and Genetic Mechanisms Underlying Sex Maintenance of Adult Chinese Alligator

The sexes of Chinese alligators are determined during embryonic development and remain fixed thereafter. In this study, we investigated the genetic and epigenetic mechanisms underlying sex maintenance in Chinese alligators through RNA sequencing and bisulfite sequencing data analyses of the adult gonads. We identified the genes and pathways (e. g., DMRT1-SOX9-AMH pathway for males and oocyte meiotic maturation pathway for females) involved in male and female sex maintenance and gonadal development of adult Chinese alligators. In contrast to their expression patterns in the embryo, both DMRT1 and the steroid hormone biosynthesis related genes showed a male-biased expression in adult gonads. The overall DNA methylation density and level were higher in testes than in ovaries. Hypermethylation in the gene bodies enhanced the expression of male-biased genes (such as DMRT1-SOX9-AMH and steroid hormone biosynthesis related genes) in the testis, as opposed to the normalization of gene expression. Our results provide insights into the genetic and epigenetic mechanisms underlying sex maintenance in adult Chinese alligators, and are expected to contribute to the development of scientific programs for the successful conservation of this endangered species.

Full Functional Sex Reversal Achieved Through Silencing of MroDmrt11E Gene in Macrobrachium rosenbergii: Production of All-Male Monosex Freshwater Prawn

The freshwater prawn Macrobrachium rosenbergii is one kind of important economic aquaculture species and displays remarkable sexual dimorphism. The molecular mechanism of sexual differentiation in M. rosenbergii has been primarily unraveled through the research efforts of the androgenic gland and its related genes. However, the understanding of conserved genes involved in the molecular mechanism underpinning sex determination and sexual differentiation of M. rosenbergii is still fragmentary. MroDmrt11E is a member of the doublesex and mab-3-related transcription factor (Dmrt) gene family and is prominently expressed in the testis. In the present study, in vivo knockdown of MroDmrt11E at the postlarva stage in male prawn induced a complete and functional sex reversal and achieved the production of an all-male monosex population. Furthermore, a great deal of new information of upregulated and downregulated transcriptions involved in sexual differentiation of MroDmrt11E knockdown was enriched by comparative transcriptomic analysis. The effects of RNAi-mediated gene knockdown of MroDmrt11E on the differentially expressed and sex-related candidate genes, such as transformer, fruitless, feminization, insulin-like androgenic gland gene, Dmrt gene family, were primarily focused on, and their possible molecular regulatory relationships in sexual differentiation were analyzed. Meanwhile, the response of primary Kyoto Encyclopedia of Genes and Genomes (KEGG) biological pathways was investigated to expound the potential roles of MroDmrt11E in male sexual differentiation, which provided a deeper understanding of the molecular regulatory network underlying sexual differentiation of M. rosenbergii. The finding provided a novel sexual manipulation technique through silencing of Dmrt gene family for achieving a complete and functional sex reversal and offered a new insight regarding the mechanism of the Dmrt gene family in the sexual differentiation of crustaceans.

De novo gonad transcriptome analysis of the common littoral shrimp Palaemon serratus: novel insights into sex-related genes

BACKGROUND

The common littoral shrimp Palaemon serratus is an economically important decapod resource in some European communities. Aquaculture practices prevent the genetic deterioration of wild stocks caused by overfishing and at the same time enhance the production. The biotechnological manipulation of sex-related genes has the proved potential to improve the aquaculture production but the scarcity of genomic data about P. serratus hinders these applications. RNA-Seq analysis has been performed on ovary and testis samples to generate a reference gonadal transcriptome. Differential expression analyses were conducted between three ovary and three testis samples sequenced by Illumina HiSeq 4000 PE100 to reveal sex-related genes with sex-biased or sex-specific expression patterns.

RESULTS

A total of 224.5 and 281.1 million paired-end reads were produced from ovary and testis samples, respectively. De novo assembly of ovary and testis trimmed reads yielded a transcriptome with 39,186 transcripts. The 29.57% of the transcriptome retrieved at least one annotation and 11,087 differentially expressed genes (DEGs) were detected between ovary and testis replicates. Six thousand two hundred seven genes were up-regulated in ovaries meanwhile 4880 genes were up-regulated in testes. Candidate genes to be involved in sexual development and gonadal development processes were retrieved from the transcriptome. These sex-related genes were discussed taking into account whether they were up-regulated in ovary, up-regulated in testis or not differentially expressed between gonads and in the framework of previous findings in other crustacean species.

CONCLUSIONS

This is the first transcriptome analysis of P. serratus gonads using RNA-Seq technology. Interesting findings about sex-related genes from an evolutionary perspective (such as Dmrt1) and for putative future aquaculture applications (Iag or vitellogenesis genes) are reported here. We provide a valuable dataset that will facilitate further research into the reproductive biology of this shrimp.

The Single-molecule long-read sequencing of Scylla paramamosain

Scylla paramamosain is an important aquaculture crab, which has great economical and nutritional value. To the best of our knowledge, few full-length crab transcriptomes are available. In this study, a library composed of 12 different tissues including gill, hepatopancreas, muscle, cerebral ganglion, eyestalk, thoracic ganglia, intestine, heart, testis, ovary, sperm reservoir, and hemocyte was constructed and sequenced using Pacific Biosciences single-molecule real-time (SMRT) long-read sequencing technology. A total of 284803 full-length non-chimeric reads were obtained, from which 79005 high-quality unique transcripts were obtained after error correction and sequence clustering and redundant. Additionally, a total of 52544 transcripts were annotated against protein database (NCBI nonredundant, Swiss-Prot, KOG, and KEGG database). A total of 23644 long non-coding RNAs (lncRNAs) and 131561 simple sequence repeats (SSRs) were identified. Meanwhile, the isoforms of many genes were also identified in this study. Our study provides a rich set of full-length cDNA sequences for S. paramamosain, which will greatly facilitate S. paramamosain research.

Ultrastructure of spermiogenesis and the distribution of spermatozoal nuclear histones in the Japanese mantis shrimp, Oratosquilla oratoria (Crustacea: Stomatopoda)

The Japanese mantis shrimp Oratosquilla oratoria (Stomatopoda; Crustacea) is one of the most economically important aquatic species of Pacific shrimp and it is distributed from Japan to the coast of China, the Philippines, the Malay Peninsula, and the Hawaiian Islands. Early studies described certain characteristics of spermatogenesis and the sperm ultrastructure in Stomatopoda, but the composition of sperm basic nuclear proteins (SBNPs) remains completely unknown. We studied the sperm ultrastructure of O. oratoria using transmission electron microscopy and the histone composition using immunofluorescence and immunoelectron microscopy. We found that the spherical nucleus is adjacent to the electron translucent external coat, which occurs in early spermatids. The acrosomal structure begins to form at the junction of the nucleus and the external coat. At the mid-spermatid stage, part of the chromatin appears to be more electron-dense than the external coat side. The aflagellate sperm of O. oratoria, are rounded or slightly ovoid in shape and have a consistent granular nucleus, an acrosome structure of pushpin shape and a spherical vesicular body in which faintly granular material is scattered. The acrosome consists of an acrosomal vesicle, perforatorium, and subacrosomal material. The sperm contains histones H2A, H2B, H3, H4, H3.3, H2AX, and H2AZ as well as some histone modifications, that is, H3K9me3, H3K4me2, H3S10ph, H4Kac, and H2A + H4S1ph. Histones are localized not only in the nucleus of the sperm but also in other structures outside the nucleus. The results may provide new perspectives for systematic studies of crustaceans and their sperm chromatin components. These findings extend the study of the sperm structure of Stomatopoda and provide basic data to elucidate the epigenetic mechanism of fertilization.

Transcriptome Analysis of Male and Female Mature Gonads of Silver Sillago (Sillago sihama)

Silver sillago (Sillago sihama) is an emerging commercial marine aquaculture species in China. To date, fundamental information on S. sihama, such as genomic information, is lacking, and no data are available on the gonad transcriptome of S. sihama. Here, the first gonadal transcriptomes of S. sihama have been constructed and genes potentially involved in gonadal development and reproduction identified. Illumina sequencing generated 60.18 million clean reads for the testis and 59.10 million for the ovary. All reads were assembled into 74,038 unigenes with a mean length of 1,004 bp and N50 value of 2,190 bp. Among all the predictable unigenes, a total of 34,104 unigenes (46%) were searched against multiple databases, including 33,244 unigenes annotated in the RefSeq Non- Redundant database at NCBI, and 28,924 in Swiss-Prot. By comparing the ovary and testis, 35,367 unigenes were identified as being differentially expressed between males and females, of which 29,127 were upregulated in the testis and 6,240 were upregulated in the ovary. Numerous differentially expressed genes (DEGs) known to be involved in gonadal development and gametogenesis were identified, including amh, dmrt1, gsdf, cyp19a1a, gnrhr, and zps. Using gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses, the top 20 KEGG pathways with highest number of DEGs were found to be involved in regulating gonadal development and gametogenesis in S. sihama. Moreover, 22,666 simple sequence repeats (SSRs) were identified in 14,577 SSR-containing sequences. The findings provide a valuable dataset for future functional analyses of sex-associated genes and molecular marker assisted selection in S. sihama.