Comparative Transcriptome Analysis Provides Insights into Differentially Expressed Genes and Long Non-Coding RNAs between Ovary and Testis of the Mud Crab (Scylla paramamosain)

Whole transcriptome RNA sequencing provides novel insights into the molecular dynamics of ovarian development in mud crab, Scylla paramamosain after mating

Ovarian development in animals is a complicated biological process, requiring the simultaneous coordination among various genes and pathways. To understand the dynamic changes and molecular regulatory mechanisms of ovarian development in mud crab (Scylla paramamosain), both histological observation and whole transcriptome sequencing of ovarian tissues at different mating stages were implemented in this study. The histological results revealed that ovarian development was delayed in unmated females (60 days after courtship behavior but not mating), who exhibited an oocyte diameter of 56.38 ± 15.17 μm. Conversely, mated females exhibited accelerated the ovarian maturation process, with females reaching ovarian stage III (proliferative stage) 23 days after mating and attained an average oocyte diameter of 132.19 ± 15.07 μm. Thus, mating process is essential in promoting the rapid ovarian development in mud crab. Based on the whole transcriptome sequencing analysis, a total of 518 mRNAs, 1502 lncRNAs, 18 circRNAs and 151 miRNAs were identified to be differentially expressed between ovarian tissues at different mating stages. Notably, six differentially expressed genes (DEGs) associated with ovarian development were identified, including ovary development-related protein, red pigment concentrating hormone receptor, G2/mitotic-specific cyclin-B3-like, lutropin-chorio gonadotropic hormone receptor, renin receptor, and SoxB2. More importantly, both DEGs and targets of differentially expressed non-coding RNAs (DEncRNAs) were enriched in renin-angiotensin system, TGF-β signaling, cell adhesion molecules, MAPK signaling pathway, and ECM-receptor interaction, suggesting that these pathways may play significant roles in the ovarian development of mud crabs. Moreover, competition endogenous RNA (ceRNA) networks were constructed while mRNAs were differentially expressed between mating stages were involved in Gene Ontology (GO) biological processes such as developmental process, reproduction, and growth. These findings could provide solid foundations for the future development of female mud crab maturation enhancement strategy, and improve the understanding of the ovarian maturation process in crustaceans.

Global expression pattern of genes containing positively selected sites in European anchovy (Engraulis encrasicolus L.) may shed light on teleost reproduction

European anchovy is a multiple-spawning and highly fecundate pelagic fish with high economic and ecological significance. Although fecundity is influenced by nutrition, temperature and weight of spawners, high reproductive capacity is related to molecular processes in the ovary. The ovary is an essential and complex reproductive organ composed of various somatic and germ cells, which interact to facilitate the development of the ovary and functional oocytes. Revealing the ovarian transcriptome profile of highly fecundate fishes provides insights into oocyte production in teleosts. Here we use a comprehensive tissue-specific RNA sequencing which yielded 102.3 billion clean bases to analyze the transcriptional profiles of the ovary compared with other organs (liver, kidney, ovary, testis, fin, cauda and gill) and juvenile tissues of European anchovy. We conducted a comparative transcriptome and positive selection analysis of seven teleost species with varying fecundity rates to identify genes potentially involved in oogenesis and oocyte development. Of the 2,272 single copies of orthologous genes found, up to 535 genes were under positive selection in European anchovy and these genes are associated with a wide spectrum of cellular and molecular functions, with enrichments such as RNA methylation and modification, ribosome biogenesis, DNA repair, cell cycle processing and peptide/amide biosynthesis. Of the 535 positively selected genes, 55 were upregulated, and 45 were downregulated in the ovary, most of which were related to RNA and DNA transferase, developmental transcription factors, protein kinases and replication factors. Overall, our analysis of the transcriptome level in the ovarian tissue of a teleost will provide further insights into molecular processes and deepen our genetic understanding of egg production in highly fecund fish.

Effects of berberine hydrochloride on immune response in the crab Charybdis japonica

Berberine hydrochloride is the main effective component of Coptis spp. used in Chinese herbal medicine and its underlying molecular mechanisms, responsible for inducing effects in crustacean species, are not fully understood. In this study, the molecular response of the crab Charybdis japonica to berberine hydrochloride exposure was studied using transcriptome sequencing. The survival rate, gene expression and activities of several immune enzymes were measured after berberine hydrochloride treatments, with or without injection of the pathogenic bacterium Aeromonas hydrophila. A total of 962 differentially expressed genes (464 up-regulated and 498 down-regulated) were observed during exposure to 100 mg/L of berberine hydrochloride and in the control group after 48 h. Enrichment analysis revealed that these genes are involved in metabolism, cellular processes, signal transduction and immune functions, indicating that exposure to berberine hydrochloride activated the immune complement system. This bioactive compound simultaneously activated fibrinogen beta (FGB), fibrinogen alpha (FGA), alpha-2-macroglobulin (A2M), kininogen (KNG), fibrinogen gamma chain (FGB), alpha-2-HS-glycoprotein (AHSG), caspase-8 (CASP8), cathepsin L (CTSL), adenylate cyclase 3 (Adcy3) and MMP1. Its action could significantly increase the survival rate of the crabs injected with A. hydrophila and promote the activity of LZM, Caspas8, FGA, ACP and AKP in the hepatopancreas. When A. hydrophila was added, the neutralization of 300 mg/L berberine hydrochloride maximized the activities of Caspas8, LZM, ACP and AKP. Our results provide a new understanding of the potential effects of berberine hydrochloride on the immune system mechanisms in crustaceans.

Full-Length Transcriptome Reconstruction Reveals the Genetic Mechanisms of Eyestalk Displacement and Its Potential Implications on the Interspecific Hybrid Crab (Scylla serrata ♀ × S. paramamosain ♂)

Simple Summary

The eyestalk is a key organ in crustaceans that produces neurohormones and regulates a range of physiological functions. Eyestalk displacement was discovered in some first-generation (F1) offspring of the novel interspecific hybrid crab (Scylla serrata ♀ × S. paramamosain ♂). To uncover the genetic mechanism underlying eyestalk displacement and its potential implications, high-quality transcriptome was reconstructed using single-molecule real-time (SMRT) sequencing. A total of 37 significantly differential alternative splicing (DAS) events (17 up-regulated and 20 down-regulated) and 1475 significantly differential expressed transcripts (DETs) (492 up-regulated and 983 down-regulated) were detected in hybrid crabs with displaced eyestalks (DH). The most significant DAS events and DETs were annotated as being endoplasmic reticulum chaperone BiP and leucine-rich repeat protein lrrA-like isoform X2. In addition, the top ten significant gene ontology (GO) terms were related to the cuticle or chitin. Overall, this study highlights the underlying genetic mechanisms of eyestalk displacement and provide useful knowledge for mud crab (Scylla spp.) crossbreeding.

Abstract

The lack of high-quality juvenile crabs is the greatest impediment to the growth of the mud crab (Scylla paramamosain) industry. To obtain high-quality hybrid offspring, a novel hybrid mud crab (S. serrata ♀ × S. paramamosain ♂) was successfully produced in our previous study. Meanwhile, an interesting phenomenon was discovered, that some first-generation (F1) hybrid offspring’s eyestalks were displaced during the crablet stage I. To uncover the genetic mechanism underlying eyestalk displacement and its potential implications, both single-molecule real-time (SMRT) and Illumina RNA sequencing were implemented. Using a two-step collapsing strategy, three high-quality reconstructed transcriptomes were obtained from purebred mud crabs (S. paramamosain) with normal eyestalks (SPA), hybrid crabs with normal eyestalks (NH), and hybrid crabs with displaced eyestalks (DH). In total, 37 significantly differential alternative splicing (DAS) events (17 up-regulated and 20 down-regulated) and 1475 significantly differential expressed transcripts (DETs) (492 up-regulated and 983 down-regulated) were detected in DH. The most significant DAS events and DETs were annotated as being endoplasmic reticulum chaperone BiP and leucine-rich repeat protein lrrA-like isoform X2. In addition, the top ten significant GO terms were related to the cuticle or chitin. Overall, high-quality reconstructed transcriptomes were obtained for the novel interspecific hybrid crab and provided valuable insights into the genetic mechanisms of eyestalk displacement in mud crab (Scylla spp.) crossbreeding.

Identification, characterization and functional analysis of gonadal long noncoding RNAs in a protogynous hermaphroditic teleost fish, the ricefield eel (Monopterus albus)

Background

An increasing number of long noncoding RNAs (lncRNAs) have been found to play important roles in sex differentiation and gonad development by regulating gene expression at the epigenetic, transcriptional and posttranscriptional levels. The ricefield eel, Monopterus albus, is a protogynous hermaphroditic fish that undergoes a sequential sex change from female to male. However, the roles of lncRNA in the sex change is unclear.

Results

Herein, we performed RNA sequencing to analyse lncRNA expression patterns in five different stages of M. albus development to investigate the roles of lncRNAs in the sex change process. A total of 12,746 lncRNAs (1503 known lncRNAs and 11,243 new lncRNAs) and 2901 differentially expressed lncRNAs (DE-lncRNAs) were identified in the gonads. The target genes of the DE-lncRNAs included foxo1, foxm1, smad3, foxr1, camk4, ar and tgfb3, which were mainly enriched in signalling pathways related to gonadal development, such as the insulin signalling pathway, MAPK signalling pathway, and calcium signalling pathway. We selected 5 highly expressed DE-lncRNAs (LOC109952131, LOC109953466, LOC109954337, LOC109954360 and LOC109958454) for full length amplification and expression pattern verification. They were all expressed at higher levels in ovaries and intersex gonads than in testes, and exhibited specific time-dependent expression in ovarian tissue incubated with follicle-stimulating hormone (FSH) and human chorionic gonadotropin (hCG). The results of quantitative real-time PCR (qRT-PCR) analysis and a dual-luciferase assay showed that znf207, as the gene targeted by LOC109958454, was expressed in multiple tissues and gonadal developmental stages of M. albus, and its expression was also inhibited by the hormones FSH and hCG.

Conclusions

These results provide new insights into the role of lncRNAs in gonad development, especially regarding natural sex changes in fish, which will be useful for enhancing our understanding of sequential hermaphroditism and sex changes in the ricefield eel (M. albus) and other teleosts.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08679-2.

Whole-transcriptome RNA sequencing revealed the roles of chitin-related genes in the eyestalk abnormality of a novel mud crab hybrid (Scylla serrata ♀ × S. paramamosain ♂)

Chitin is a kind of insoluble structural polysaccharide and plays different roles in different species. In crustaceans, it forms the structural components in the exoskeleton. In our previous studies, novel mud crab hybrids have been produced from the interspecific hybridization of Scylla serrata ♀ × S. paramamosain ♂. Some of the hybrid crabs have been found to be morphologically (eyestalk) abnormal, but the genetic mechanism remains unknown. To address this question, we performed whole-transcriptome RNA sequencing on the control group (normal hybrids), abnormal hybrids, and S. paramamosain to uncover the genetic basis underlying this morphological abnormality. A total of 695 mRNAs, 10 miRNAs, 44 circRNAs, and 1957 lncRNAs were differentially expressed between normal and abnormal hybrids. Several differentially expressed genes (DEGs) associated with chitin and cuticle metabolism were identified, including chitin synthase, chitinase, chitin deacetylase, β-N-acetylglucosaminidase, β-1,4-endoglucanase, N-alpha-acetyltransferase, cuticle proprotein, early cuticle protein, and arthrodial cuticle protein. Functional analysis showed that DE miRNAs, DE circRNAs, DE lncRNAs, and lncRNA/circRNA-miRNA-mRNA network were enriched in pathways related to the amino acid, carbohydrate, and glycogen metabolism. Considering the importance of the chitin and cuticle in exoskeleton formation, it can be concluded that the changes in the chitin and cuticle biosynthesis might have caused the eyestalk abnormality in hybrid crabs. These findings can lay the solid foundation for a better understanding of the important roles of chitin and cuticle related genes and the development of hybridization techniques in crustaceans.

Cloning, expression, and function of the Spdmrt-like gene in Scylla paramamosain

BACKGROUND

The mud crab Scylla paramamosain is an economically important species for aquaculture in China and has sexually dimorphic between females and males. Understanding sex differentiation in this species is essential for the development of monosex aquaculture. The Dmrt genes play a vital role in sex differentiation in animals.

METHODS AND RESULTS

In this study, two dmrt-like transcript variants, Spdmrt-like-tv1 and Spdmrt-like-v2, were cloned. SpDmrt-like-tv1 contained a DM domain, while SpDmrt-like-tv2 contained a DM and a DMA domain. Spdmrt-like-tv1 and Spdmrt-like-tv2 were both specifically expressed in testis. During testicular development, the expression level of Spdmrt-like-tv1 increased from stage I to stage II (P > 0.05) and then decreased from stage II to stage III (P < 0.05). The expression level of Spdmrt-like-tv2 in stages I and II was significantly higher than that in stage III (P < 0.05). During embryonic development, the expression level of Spdmrt-like-tv1 was higher in the mid-embryonic stage compared with the early and late stages, but the differences were not significant. Moreover, the expression level of Spdmrt-like-tv2 was stable and remained high throughout embryonic development. Furthermore, the expression level of Spdmrt-like-tv2 was significantly higher than that of Spdmrt-like-tv1. Knockdown of Spdmrt-like variants indicated that the regulative target gene of Spdmrt-like-tv1 was Spsox21, and the regulative target genes of Spdmrt-like-tv2 were Spfoxl2 and Spsox21. Combined with the results in our previously published peer-reviewed articles that the expression of Spfoxl2 in the testis was significantly higher than that in the ovary, and Spfoxl2 negatively regulated Spvtg expression. Spsox21 played a role in the development and maintenance of testis as well as in the process of neural development and regulation of body segmentation.

CONCLUSION

Therefore, we suggest that Spdmrt-like-tv1 and Spdmrt-like-tv2 might be involved in testicular development and embryonic development, and Spdmrt-like-tv2 might play more important roles in these two developmental processes by regulating the expression of Spfoxl2 and Spsox21 due to its high expression.

Long non-coding RNA profile in banana shrimp, Fenneropenaeus merguiensis and the potential role of lncPV13 in vitellogenesis

The long non-coding RNAs (lncRNAs) have been known to play important roles in several biological processes as well as in reproduction. This study aimed to identify lncRNA in ovary female banana shrimp, Fenneropenaeus merguiensis, and investigate the potential role of lncPV13 in the vitellogenesis. After the in silico identification of the ovarian transcriptome, a total of 24,733 putative lncRNAs were obtained, and only 147 putative lncRNAs were significantly differentially expressed among the ovarian development stages. To validate the in silico identification of lncRNAs, the 16 lncRNAs with the highest differential expression in the transcriptome analysis were evaluated by RT-qPCR. The 6 lncRNAs showed higher expression levels in the mature stage than in the previtellogenic stage and were found in several tissues such as in eyestalks, brains, thoracic ganglia, gills, and muscle. Furthermore, most candidate lncRNAs were amplifiable in Litopenaeus vannamei's and Penaeus monodon's DNA but not in Macrobrachium rosenbergii's DNA, suggesting some lncRNAs are expressed in a species-specific manner among penaeid shrimp. In this study, the lncPV13 was investigated for its vitellogenin regulating function by RNA interference. The result indicates that the lncPV13 expression was suppressed in the ovary on day 7 after the injection of double-stranded RNA specific to lncPV13 (dslncPV13), while vitellogenin (Vg) expression was significantly decreased. In contrast, the gonad inhibiting hormone (GIH) expression was significantly increased in the lncPV13 knockdown shrimp. However, the oocyte proliferation was not significantly different between control and lncPV13 knockdown shrimp. This suggests that lncPV13 regulate Vg synthesis through GIH inhibition. Finally, our findings provide lncRNA information and potential lncRNAs involved in the vitellogenesis of female banana shrimp.

Evidence of Sex Differentiation Based on Morphological Traits During the Early Development Stage of Mud Crab Scylla paramamosain

In order to uncover the sexual difference in morphology and how early they appear during the development stage of mud crab Scylla paramamosain, we measured, observed, and biostatistically analyzed morphological traits related to sex. For unveiling the morphological differences between sexes, morphological traits involving abdomen width (AW), carapace length (CL), and carapace width (CW) were first measured during the crablet development stage of S. paramamosain in the present study. The correlation analyses and path analyses exhibited that sexual dimorphism in the third abdomen width (AW3) and fourth abdomen width (AW4) could be used for sex identification from stage C VI (stage VI of crablet). Based on the stepwise discriminant analysis and standardized traits, a sex discriminant equation was constructed, which is capable for sex identification in crablets from stage C VI. Observations for secondary sexual traits and abdomen morphology (shape and pleopods) using a dissecting microscope or scanning electron microscope indicated that sexes are easily identified at stage C VIII according to the abdomen shape; meanwhile, at stage C II based on pleopod difference, and at stage C I by the presence or absence of gonopores. The findings in this study contribute greatly to the accuracy of sex identification of S. paramamosain during the early development stage, which promotes the understanding of the morphological differentiation mechanism of sex.

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.

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.

Comparison of morphology and genetic diversity between broodstock and hybrid offspring of oriental river prawn, Macrobrachium nipponense based on morphological analysis and SNP markers

Hybridization is an effective method for the genetic improvement of farmed species. In this study, three broodstock populations (Changjiang, CJ, Dongting, DT, and Dianshan, DS) of oriental river prawn, Macrobrachium nipponense were used, and DS was used as the female broodstock. Through three-line hybridization, two hybrid populations were finally obtained. The F3 generation of the broodstock population and the F1 generation of the hybrid population were cultured indoors for 3 months. Through morphological analysis (cluster analysis, discriminant analysis and path analysis) it was found that the hybrid population and the broodstock had some differences, but not enough to reach the subspecies level, and the dominant traits exhibited differentiation and reorganization. This study identified SNP genetic markers, carried out systematic evolution analysis and genetic diversity analysis and found that the nucleotide diversity π and heterozygosity Het of the hybrid population were higher than those of the broodstock. Among broodstocks, the differentiation index (F_st ) of SCD and SDC was smallest (0.055). This research provides some valuable reference for genetic breeding.

The significant sex-biased expression pattern of Sp-Wnt4 provides novel insights into the ovarian development of mud crab (Scylla Paramamosain)

The wingless-type MMTV integration site family member-4 (Wnt4), a member of the wingless-related integration site (Wnt) family, is widely accepted as a key regulator of ovarian development in mammals. In this study, a full-length cDNA of Wnt4 (designated as Sp-Wnt4) was cloned, characterized, and functionally studied in mud crab (Scylla paramamosain). The full-length cDNA of Sp-Wnt4 consists of 2659 bp with an open reading frame (ORF) encoding 359 amino acids, a 907 bp 5'-UTR and a 672 bp 3'-UTR. Sp-Wnt4 contains 25 cysteine (Cys) residues and three potential N-glycosylation sites. Sp-Wnt4 protein shared the highest identity (98.9%) to the Wnt4 protein of Portunus trituberculatus. The phylogenetic tree showed that Sp-Wnt4 and Wnt4 protein of Malacostracan crustaceans clustered together, indicating that they had a close genetic distance. Sp-Wnt4 was expressed at a higher level in the ovary compared to other tissues, with the highest expression level at the third stage (O-III) of the ovarian development (P < 0.05). A downward trend was observed in the expression level of Sp-Wnt4 from the embryo stage to crablet stages (P < 0.05). After unilateral eyestalk ablation, the expression level of Sp-Wnt4 significantly increased in testis (14-fold) and downregulated (3.1-fold) in the gill (P < 0.05) of females. In situ hybridization (ISH) assay revealed that Sp-Wnt4 transcripts were mainly localized in the cytoplasm of oocyte cells. These findings showed that Sp-Wnt4 play crucial roles in the ovarian development of S. paramamosain. In conclusion, our study provides novel insights into the evolution and roles of the Wnt4 gene.

Discovery of the Dmrt gene family members based on transcriptome analysis in mud crab Scylla paramamosain

Doublesex and mab-3 related transcription factors (Dmrts) play crucial roles in sex determination/differentiation and gonad development. The information on Dmrts and their functions are still scarce in mud crab Scylla paramamosain. In this study, 12 published transcriptome data of S. paramamosain were retrieved, pooled, and assembled. From the assembly, 7 Dmrt gene family members were identified and consisted of Spdmrt-like, Spdmrt-1a, Spdmrt-3, Spdmrt-11E, Spidmrt-1, Spdoublesex (Spdsx), and Spidmrt-2. These dmrt genes were predicted to encode 224 aa, 465 aa, 435 aa, 276 aa, 520 aa, 552 aa, and 266 aa protein precursors, respectively. The expression patterns of the dmrt genes were characterized by semi-quantitative PCR. The Spdmrt-like and Spdmrt-1a were exclusively detected in gonads, of which both expression levels in the testis were higher than that in the ovary. The Spdmrt-3, Spdmrt-11E, Spidmrt-1, Spdsx, and Spidmrt-2 were observed in various tissues; all these genes were sexually dimorphic except for dmrt-11E. Specifically, the expression level of Spdmrt-3 and Spidmrt-2 were higher in the testis than that in the ovary. On the contrary, the Spdsx and Spidmrt-1 expression level were higher in ovary than that in testis. The present study's findings provided a fundamental understanding of Dmrt gene family members involving sex determination/differentiation and gonad development in the S. paramamosain.

De novo transcriptome assembly of the Southern Ocean copepod Rhincalanus gigas sheds light on developmental changes in gene expression

Copepods are small crustaceans that dominate most zooplankton communities in terms of both abundance and biomass. In the polar oceans, a subset of large lipid-storing copepods occupy central positions in the food web because of their important role in linking phytoplankton and microzooplankton with higher trophic levels. In this paper, we generated a high-quality de novo transcriptome for Rhincalanus gigas, the largest-and among the most abundant-of the Southern Ocean copepods. We then conducted transcriptional profiling to characterize the developmental transition between late-stage juveniles and adult females. We found that juvenile R. gigas substantially upregulate lipid synthesis and glycolysis pathways relative to females, as part of a developmental gene expression program that also implicates processes such as muscle growth, chitin formation, and ion transport. This study provides the first transcriptional profile of a developmental transition within Rhincalanus gigas or any endemic Southern Ocean copepod, thereby extending our understanding of copepod molecular physiology.

Comparative transcriptomic analysis for identification of candidate sex-related genes and pathways in Crimson seabream (Parargyrops edita)

Teleost fishes display the largest array of sex-determining systems among animals, resulting in various reproductive strategies. Research on sex-related genes in teleosts will broaden our understanding of the process, and provide important insight into the plasticity of the sex determination process in vertebrates in general. Crimson seabream (Parargyrops edita Tanaka, 1916) is one of the most valuable and abundant fish resources throughout Asia. However, little genomic information on P. edita is available. In the present study, the transcriptomes of male and female P. edita were sequenced with RNA-seq technology. A total of 388,683,472 reads were generated from the libraries. After filtering and assembling, a total of 79,775 non redundant unigenes were obtained with an N50 of 2,921 bp. The unigenes were annotated with multiple public databases, including NT (53,556, 67.13%), NR (54,092, 67.81%), Swiss-Prot (45,265, 56.74%), KOG (41,274, 51.74%), KEGG (46,302, 58.04%), and GO (11,056, 13.86%) databases. Comparison of the unigenes of different sexes of P. edita revealed that 11,676 unigenes (9,335 in females, 2,341 in males) were differentially expressed between males and females. Of these, 5,463 were specifically expressed in females, and 1,134 were specifically expressed in males. In addition, the expression levels of ten unigenes were confirmed to validate the transcriptomic data by qRT-PCR. Moreover, 34,473 simple sequence repeats (SSRs) were identified in SSR-containing sequences, and 50 loci were randomly selected for primer development. Of these, 36 loci were successfully amplified, and 19 loci were polymorphic. Finally, our comparative analysis identified many sex-related genes (zps, amh, gsdf, sox4, cyp19a, etc.) and pathways (MAPK signaling pathway, p53 signaling pathway, etc.) of P. edita. This informative transcriptomic analysis provides valuable data to increase genomic resources of P. edita. The results will be useful for clarifying the molecular mechanism of sex determination and for future functional analyses of sex-associated genes.

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.

The long non-coding RNA TSLC8 inhibits colorectal cancer by stabilizing puma

The colorectal cancer (CRC) dictates a common malignancy with high recurrence rate. Long non-coding RNAs (lncRNAs) belong to a class of regulatory factors involved in multiple cancers. In current work, we have uncovered a novel lncRNA named TSLC8. TSLC8 was dramatically downregulated in CRC samples and cell lines. Reintroduction of TSLC8 inhibited tumor sphere formation and viability in CRC cells. In vivo experiments further confirmed the tumor suppressive function of TSLC8. Ectopic TSLC8 expression elevates puma abundance whereas this effect is mediated by TSLC8-puma binding and stabilization. FOXO1 can transcriptionally induce TSLC8 expression. Epigenetic investigation suggested that TSLC8 locus was hypermethylated in CRC leading to diminished TSLC8 expression. Our current work has identified a novel tumor suppressive function of TSLC8, whose reduced expression may facilitate malignant phenotypes during CRC progression.

Genome-Wide Characterization of Alternative Splicing Events and Their Responses to Cold Stress in Tilapia

Alternative splicing (AS) is an important post-transcriptional regulatory mechanism for cells to generate transcript variability and proteome diversity. No systematic investigation of AS events among different tissues in response to stressors is available for tilapia currently. In this study, AS among different tissues was identified and the cold stress-related AS events were explored in a Nile tilapia (Oreochromis niloticus) line based on 42 RNA-seq datasets using a bioinformatics pipeline. 14,796 (82.76%; SD = 2,840) of the expression genes showed AS events. The two most abundant AS types were alternative transcription start site (TSS) and terminal site (TTS) in tilapia. Testis, brain and kidney possess the most abundant AS gene number, while the blood, muscle and liver possess the least number in each tissue. Furthermore, 208 differentially alternative splicing (DAS) genes in heart and 483 DAS in brain in response to cold stress. The number of AS types for alternative exon end, exon skipping and retention of single intron increased significantly under cold stress. GO enrichment and pathway overrepresentation analysis indicated that many DAS genes, e.g., genes in circadian clock pathway, may influence expression of downstream genes under cold stress. Our study revealed that AS exists extensively in tilapia and plays an important role in cold adaption.

Transcriptome analysis and histopathology of the mud crab (Scylla paramamosain) after air exposure

The mud crab, Scylla paramamosain, is an economically-important crab in China. Air exposure is an important environmental stressor during mud crab culture and transportation. Adaptive mechanisms responding to air exposure in mud crabs are still poorly understood. In this study, mud crabs were exposed to air for 120 h. Air exposure decreased total hemocyte counts, led to cytological damage, and caused high mortality. Transcriptomic analysis was conducted at 0, 6 and 96 h after air exposure. A total of 3530 differentially expressed genes (DEGs) were identified. DEGs were mainly involved in the oxidative stress response, metabolism, cellular processes, signal transduction, and immune functions. Transcriptomic analysis also revealed that genes of glycolysis and of the tricarboxylic acid cycle were key factors in regulating the mud crab adaptation to air exposure.

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.

Gonadal microRNA Expression Profiles and Their Potential Role in Sex Differentiation and Gonadal Maturation of Mud Crab Scylla paramamosain

Although the sexual dimorphism in terms of gonadal development and gametogenesis of mud crab has been described, the internal regulating mechanism and sex differentiation process remain unclear. A comparative gonadal miRNA transcriptomic study was conducted to identify miRNAs that are differentially expressed between testes and ovaries, and potentially uncover miRNAs that might be involved in sex differentiation and gonadal maturation mechanisms of mud crabs (Scylla paramamosain). A total of 10 known miRNAs and 130 novel miRNAs were identified, among which 54 were differentially expressed. Target gene prediction revealed a significant enrichment in 30 KEGG pathways, including some reproduction-related pathways, e.g. phosphatidylinositol signalling system and inositol phosphate metabolism pathways. Further analysis on six differentially expressed known miRNAs, six differentially expressed novel miRNAs and their reproduction-related putative target genes shows that both miRNAs and putative target genes showed stage-specific expression during gonadal maturation, suggesting their potential regulatory roles in sex differentiation and reproductive development. This study reveals the sex-biased miRNA profile and establishes a solid foundation for understanding the sex differentiation and gonadal maturation mechanisms of S. paramamosain.

High-Density Genetic Linkage Maps Provide Novel Insights Into ZW/ZZ Sex Determination System and Growth Performance in Mud Crab (Scylla paramamosain)

Mud crab, Scylla paramamosain is one of the most important crustacean species in global aquaculture. To determine the genetic basis of sex and growth-related traits in S. paramamosain, a high-density genetic linkage map with 16,701 single nucleotide polymorphisms (SNPs) was constructed using SLAF-seq and a full-sib family. The consensus map has 49 linkage groups, spanning 5,996.66 cM with an average marker-interval of 0.81 cM. A total of 516 SNP markers, including 8 female-specific SNPs segregated in two quantitative trait loci (QTLs) for phenotypic sex were located on LG32. The presence of female-specific SNP markers only on female linkage map, their segregation patterns and lower female: male recombination rate strongly suggest the conformation of a ZW/ZZ sex determination system in S. paramamosain. The QTLs of most (90%) growth-related traits were found within a small interval (25.18–33.74 cM) on LG46, highlighting the potential involvement of LG46 in growth. Four markers on LG46 were significantly associated with 10–16 growth-related traits. BW was only associated with marker 3846. Based on the annotation of transcriptome data, 11 and 2 candidate genes were identified within the QTL regions of sex and growth-related traits, respectively. The newly constructed high-density genetic linkage map with sex-specific SNPs, and the identified QTLs of sex- and growth-related traits serve as a valuable genetic resource and solid foundation for marker-assisted selection and genetic improvement of crustaceans.

Effect of fasting and subsequent refeeding on the transcriptional profiles of brain in juvenile Spinibarbus hollandi

Starvation is a common stress in fish. The underlying molecular mechanisms associated with growth depression caused by feeding restriction and compensatory growth are not well understood. We investigated the effect of fasting and refeeding on the transcriptome profiles of brain in juvenile S. hollandi using RNA-seq. A total of 4.73 × 10⁸ raw reads were obtained from nine brain samples. De novo transcriptome assembly identified 387,085 unigenes with 2.1×10⁹ nucleotides. A total of 936 annotated unigenes showed significantly differential expression among the control, fasting, and fasting-refeeding groups. The down-regulated differentially expressed genes (DEGs) during fasting were mainly associated with cell cycle, DNA replication, and mitosis. The up-regulated DEGs were mainly related to glucose and lipid metabolism, material transportation, and transcription factors. Most decreased DEGs during fasting were restored to normal levels after refeeding. Comparing with the control group, genes associated with protein synthesis, stimulus response, and carbohydrate metabolism were significantly over-expressed and pro-opio melanocortin (POMC) was down-regulated during the refeeding period. In conclusion, fish mobilized stored energetic materials and reduced energy consumption to prolong survival during fasting. After refeeding, the down-regulation of DEGs, e.g., POMC may be associated with compensatory growth. Up-regulation of DEGs related to ribosomal protein, stimulus response, and carbohydrate metabolism may contribute to eliminate negative effect of starvation on brain. This study provided the first transcriptome data related with impact of short-time starvation and refeeding in S. hollandi brains.

Identification of male-specific SNP markers and development of PCR-based genetic sex identification technique in crucifix crab (Charybdis feriatus) with implication of an XX/XY sex determination system

In this study, we first identified male-specific SNP markers using restriction site-associated DNA sequencing, and further developed a PCR-based sex identification technique for Charybdis feriatus. A total of 296.96 million clean reads were obtained, with 114.95 and 182.01 million from females and males. After assembly and alignment, 10 SNP markers were identified being heterozygous in males but homozygous in females. Five markers were further confirmed to be male-specific in a large number of individuals. Moreover, two male-specific sense primers and a common antisense primer were designed, using which, a PCR-based genetic sex identification method was successfully developed and used to identify the sex of 103 individuals, with a result of 49 females and 54 males. The presence of male-specific SNP markers suggests an XX/XY sex determination system for C. feriatus. These findings should be helpful for better understanding sex determination mechanism, and drafting artificial breeding program in crustaceans.

Gonadal Transcriptome Analysis of Pacific Abalone Haliotis discus discus: Identification of Genes Involved in Germ Cell Development

Little is known about the molecular mechanisms governing gonadal developmental processes in abalones. Here, we conducted transcriptome analysis of Pacific abalone Haliotis discus discus for gene discovery in the brain, ovary, testis, and unfertilized eggs. Among the annotated unigenes, 48.6% of unigenes were identified by Venn diagram analysis as having universal or tissue-specific expression. Twenty-three genes with gonad-biased gene ontology (GO) terms were first obtained. Secondly, 36 genes were found by screening known gene names related to germ cell development. Finally, 17 genes were obtained by querying the annotated unigene database for zygotically expressed gonadal genes (ovary and testis) and maternally expressed gonadal genes (ovary, testis, and unfertilized eggs) using keywords related to reproduction. To further verify tissue distribution pattern and subcellular localization of these genes, RT-PCR and in situ hybridization were performed using a unigene encoding a germ cell marker, vasa, as control. The results showed that vasa was expressed mainly in the early developmental stages of germ cells in both sexes. One of the candidate genes, vitelline envelope zona pellucida domain protein 12 (ZP12), was expressed in the primordial germ cells of immature gonad and early developmental stages of germ cells of the adult female. The results obtained from the present study suggest that vasa and ZP12 are involved in germ cell development of Pacific abalone and that ZP12 is an especially useful germ cell-specific marker in immature adults. The current gonadal transcriptome profile is an extensive resource for future reproductive molecular biology studies of this species.

Transcriptome sequencing and molecular markers discovery in the gonads of Portunus sanguinolentus

Crab culture has gained prominence in the last decade due to the large global market demand for live crabs and crab products. Portunus sanguinolentus is one of the economically important crab species in the Indo-Pacific region, with distinct differences in growth and size between male and female crabs, thus, leading to huge difference in their market values. The culture of P. sanguinolentus is still in its infancy, with crab supplies heavily dependent on wild catch. In order to unravel the molecular differences between male and female crabs, we generated a comprehensive transcriptomic dataset for P. sanguinolentus by sequencing the gonads of both sexes using the Illumina Hiseq 2500 system. Transcriptomes were assembled using Trinity de novo assembly followed by annotation. This transcriptomic data set for P. sanguinolentus would serve as an important reference data for genomic and genetic studies in this crab and related species.