Comparative Transcriptome Analysis Reveals Sex-Biased Gene Expression in Juvenile Chinese Mitten Crab Eriocheir sinensis

Genetic variation mining of the Chinese mitten crab (Eriocheir sinensis) based on transcriptome data from public databases

At present, public databases house an extensive repository of transcriptome data, with the volume continuing to grow at an accelerated pace. Utilizing these data effectively is a shared interest within the scientific community. In this study, we introduced a novel strategy that harnesses SNPs and InDels identified from transcriptome data, combined with sample metadata from databases, to effectively screen for molecular markers correlated with traits. We utilized 228 transcriptome datasets of Eriocheir sinensis from the NCBI database and employed the Genome Analysis Toolkit software to identify 96 388 SNPs and 20 645 InDels. Employing the genome-wide association study analysis, in conjunction with the gender information from databases, we identified 3456 sex-biased SNPs and 639 sex-biased InDels. The KOG and KEGG annotations of the sex-biased SNPs and InDels revealed that these genes were primarily involved in the metabolic processes of E. sinensis. Combined with SnpEff annotation and PCR experimental validation, a highly sex-biased SNP located in the Kelch domain containing 4 (Klhdc4) gene, CHR67-6415071, was found to alter the splicing sites of Klhdc4, generating two splice variants, Klhdc4_a and Klhdc4_b. Additionally, Klhdc4 exhibited robust expression across the ovaries, testes, and accessory glands. The sex-biased SNPs and InDels identified in this study are conducive to the development of unisexual cultivation methods for E. sinensis, and the alternative splicing event caused by the sex-biased SNP in Klhdc4 may serve as a potential mechanism for sex regulation in E. sinensis. The analysis strategy employed in this study represents a new direction for the rational exploitation and utilization of transcriptome data in public databases.

Differential expression of sNPF in male and female eyestalk leading to sex dimorphism of AMP expression in Procambarus clarkii intestine

Antimicrobial peptide (AMP) is an important component of crustaceans' innate immune system. In this study, a short neuropeptide F (sNPF) gene (Pc-sNPF) and a Forkhead box O (FOXO) gene (PcFOXO) from Procambarus clarkii were identified. Analysis findings showed that the expression level of AMP genes differed between male and female P. clarkii. Furthermore, Pc-sNPF and PcFOXO were related to the sex dimorphism of AMP. Knockdown of Pc-sNPF in the eyestalk significantly upregulated the expression of PcFOXO and two anti-lipopolysaccharide factors (PcALF4 and PcALFL) in the intestine of P. clarkii. The expression of PcFOXO in the intestine of female P. clarkii was higher than in that of males. Results from RNA interference revealed that PcFOXO positively regulated the expression of PcALF4 and PcALFL in the intestine of male and female P. clarkii. In summary, our study showed that differences in Pc-sNPF expression in eyestalk of male and female P. clarkii leading to sex dimorphism of AMP expression in the intestine are mediated by the sNPF-FOXO-AMP signal pathway called the eyestalk-intestine axis.

Sex-bias of core intestinal microbiota in different stocks of Chinese mitten crabs (Eriocheir sinensis)

The differences in intestinal microbiota composition are synergistically shaped by internal and external factors of the host. The core microbiota plays a vital role in maintaining intestinal homeostasis. In this study, we conducted 16S rRNA sequencing analysis to investigate the stability of intestinal microbiota and sex-bias of six stocks of Chinese mitten crabs (105 females; and 110 males). The dominant phyla in all six stocks were Proteobacteria, Tenericutes, Bacteroidetes and Firmicutes; however, their relative abundance differed significantly. Twenty-seven core operational taxonomic units (OTUs), corresponding to 18 genera, were screened. Correlation analysis revealed that OTUs of four stocks in the Yangtze River system play important roles in maintaining the stability of intestinal microbiota. Additionally, the core intestinal microbiota was significantly sex-biased, and the top three genera in terms of relative abundance (Acinetobacter, Vibrio, and Candidatus_Hepatoplasma) were significantly dominant in female crabs. Network structure analysis also confirmed gender differences in the association pattern of intestinal microbiota. The intestinal microbiota of male crabs has a higher degree of functional enrichment. This study provided a theoretical basis for further investigating exploring the shaping effect of gender and geographical factors on the intestinal microbiota of Chinese mitten crabs.

Transcriptome analysis of reproductive tract tissues of male river prawn Macrobrachium americanum

BACKGROUND

The river prawn, Macrobrachium americanum (M. americanum), is one of the largest prawns of the genus in Latin America and is an amphidromous species distributed along the Pacific coast of America. This prawn has commercial value due to its size and taste, making it a good option for aquaculture production. Its culture has been attempted in ponds and concrete tanks, but no successful technique can still support commercial production. Understanding the mechanisms that regulate reproduction at the molecular level is very important. This knowledge can provide tools for manipulating transcripts, which could increase the number or size of animals in the culture. Our understanding of the mechanism that regulates the reproduction of M. americanum at the molecular level is limited.

AIM

Perform and analyze the transcriptome assembly of the testes, vas deferens, and terminal ampulla of M. americanum. to provide new molecular information about its reproduction.

METHODS AND RESULTS

The cDNA library was constructed and sequenced for each tissue to identify novel transcripts. A combined transcriptome with the three tissues was assembled using Trinity software. Unigenes were annotated using BLASTx and BLAST2GO. The transcriptome assembly generated 1,059,447 unigenes, of which 7222 genes had significant hits (e-value < 1 × 10-5) when compared against the Swiss-Prot database. Around 75 genes were related to sex determination, testis development, spermatogenesis, spermiogenesis, fertilization, maturation of testicular cells, neuropeptides, hormones, hormone receptors, and/or embryogenesis.

CONCLUSIONS

These results provide new molecular information about M. americanum reproduction, representing a reference point for further genetic studies of this species.

CRISPR/Cas9-mediated gene mutation of EcIAG leads to sex reversal in the male ridgetail white prawn Exopalaemon carinicauda

In the culture of crustaceans, most species show sexual dimorphism. Monosex culture is an effective approach to achieve high yield and economic value, especially for decapods of high value. Previous studies have developed some sex control strategies such as manual segregation, manipulation of male androgenic gland and knockdown of the male sexual differentiation switch gene encoding insulin-like androgenic gland hormone (IAG) in decapods. However, these methods could not generate hereditable changes. Genetic manipulation to achieve sex reversal individuals is absent up to now. In the present study, the gene encoding IAG (EcIAG) was identified in the ridgetail white prawn Exopalaemon carinicauda. Sequence analysis showed that EcIAG encoded conserved amino acid structure like IAGs in other decapod species. CRISPR/Cas9-mediated genome editing technology was used to knock out EcIAG. Two sgRNAs targeting the second exon of EcIAG were designed and microinjected into the prawn zygotes or the embryos at the first cleavage with commercial Cas9 protein. EcIAG in three genetic males was knocked out in both chromosome sets, which successfully generated sex reversal and phenotypic female characters. The results suggest that CRISPR/Cas9-mediated genome editing technology is an effective way to develop sex manipulation technology and contribute to monosex aquaculture in crustaceans.

Sex-specific effects of propiconazole on the molting of the Chinese mitten crab (Eriocheir sinensis)

Given the inevitable exposure of Eriocheir sinensis (E. sinensis) to fungicides in rice-crab co-culture systems, understanding the potential effect of fungisides is important for practical application. Molting is a crucial development process of E. sinensis, which is regulated by endocrine system and genetic factors, and is susceptible to exogenous chemicals. However, the impact of fungicides application on the molting of E. sinensis have been rarely reported. In the present study, propiconazole, a widely used fungicide for rice disease management, was found to exert potential effects on the molting of E. sinensis at residual-related level in the rice-crab co-culture fields. After 14 days of short-term exposure to propiconazole, female crabs exhibited remarkably higher levels of hemolymph ecdysone than males. When the exposure was extended to 28 days, propiconazole markedly accelerated molt-inhibiting hormone expression by 3.3-fold, ecdysone receptor expression by 7.8-fold, and crustacean retinoid X receptor expression by 9.6-fold in male crabs, while it showed the opposite effect in females with suppressed gene expression. Propiconazole also induced the activity of N-acetylglucosaminidase in male crabs rather than females during the experiments. Our study suggests that propiconazole exerts sex-specific effects on the molting of E. sinensis. The impact of propiconazole application in the rice-crab co-culture systems remains more assessment to avoid affecting the growth of cultured E. sinensis.

Interactions Between Temperature Variability and Reproductive Physiology Across Traits in an Intertidal Crab

Thermal extremes alter population processes, which can result in part from temperature-induced movement at different spatial and temporal scales. Thermal thresholds for animal movement likely change based on underlying thermal physiology and life-history stage, a topic that requires greater study. The intertidal porcelain crab Petrolisthes cinctipes currently experiences temperatures that can reach near-lethal levels in the high-intertidal zone at low tide. However, the thermal thresholds that trigger migration to cooler microhabitats, and the extent to which crabs move in response to temperature, remain unknown. Moreover, the influence of reproductive status on these thresholds is rarely investigated. We integrated demographic, molecular, behavioral, and physiological measurements to determine if behavioral thermal limits varied due to reproductive state. Demographic data showed a trend for gravid, egg bearing, crabs to appear more often under rocks in the cooler intertidal zone where crab density is highest. In situ expression of 31 genes related to stress, metabolism, and growth in the field differed significantly based on intertidal elevation, with mid-intertidal crabs expressing the gene for the reproductive yolk protein vitellogenin (vg) earlier in the season. Furthermore, VG protein levels were shown to increase with density for female hemolymph. Testing for temperatures that elicit movement revealed that gravid females engage in heat avoidance behavior at lower temperatures (i.e., have a lower voluntary thermal maximum, VTmax) than non-gravid females. VTmax was positively correlated with the temperature of peak firing rate for distal afferent nerve fibers in the walking leg, a physiological relationship that could correspond to the mechanistic underpinning for temperature dependent movement. The vulnerability of marine organisms to global change is predicated by their ability to utilize and integrate physiological and behavioral strategies in response to temperature to maximize survival and reproduction. Interactions between fine-scale temperature variation and reproductive biology can have important consequences for the ecology of species, and is likely to influence how populations respond to ongoing climate change.

Molecular characterization and expression profiling of transformer 2 and fruitless-like homologs in the black tiger shrimp, Penaeus monodon

Transformer 2 (tra 2) and fruitless (fru) genes have been proven to play a key role in sex determination pathways in many Arthropods, including insects and crustaceans. In this study, a paralog of P. monodon tra 2 (Pmtra 2), P. monodon ovarian associated transformer 2 (PmOvtra 2) and 2 isoforms of P. monodon fruitless-like gene (Pmfru-1 and Pmfru-2) were identified and characterized. The full cDNA sequence of PmOvtra 2 consisted of 1,774 bp with the longest open reading frame (ORF) of 744 bp encoding for 247 amino acids. The PmOvtra 2 exhibited a predicted RNA-recognition motif (RRM) domain and two arginine-serine (RS) regions, suggesting its function in RNA splicing. The full cDNA sequence of Pmfru-1 consisted of 1,306 bp with 1,182 bp ORF encoding for 393 amino acids, whereas the full cDNA sequence of Pmfru-2 consisted of 1,858 bp with 1,437 bp ORF encoding 478 amino acids. The deduced amino acid sequences of Pmfru-1 and Pmfru-2 exhibited highly conserved domains of Fru proteins, including Broad-complex, Tramtrack and Bric-a-brac (BTB), and zinc finger (ZF) domains. In addition, Pmfru-1 and Pmfru-2 were suggestively originated from the same single genomic locus by genomic sequence analysis. Specifically, Pmfru pre-mRNA was alternatively spliced for Pmfru-1 and Pmfru-2 to include mutually exclusive exon 7 and exon 6, respectively. Temporal and spatial expression of PmOvtra 2, Pmfru-1, and Pmfru-2 were also investigated by qPCR. The results showed that all were expressed in early developmental stages with undifferentiated gonads starting from nauplius until postlarvae. The expression of PmOvtra 2 started at nauplius stage and gradually increased from mysis to postlarvae (PL) 1. However, the expression of Pmfru-1 was low at the nauplii stage and slightly increased from protozoea to PL5, whereas the expression of Pmfru-2 maintained a low level from nauplius to mysis and then gradually increased at the PL stages. Expressions of PmOvtra 2, Pmfru-1, and Pmfru-2 were detected in various tissues including nervous tissue, gill, heart, hepatopancreas, gut, and gonads. Interestingly, the sexually dimorphic expression of PmOvtra 2, Pmfru-1, and Pmfru-2 was demonstrated in fully developed gonads in which the ovary showed significantly higher expressions than the testis. The great difference in the expression pattern of PmOvtra 2, Pmfru-1, and Pmfru-2 in the ovary and testis suggested their roles in the female sex determination in P. monodon.

Effect of infection with Aeromonas hydrophila on antioxidant capacity, inflammation response, and apoptosis proteins in Chinese mitten crab (Eriocheir sinensis)

Aeromonas hydrophila (A. hydrophila) as a serious bacterial disease endangering aquaculture and the Chinese mitten crabs (Eriocheir sinensis) industry. The present study was conducted to investigate the effects of A. hydrophila on the antioxidant, inflammation, immunity and apoptosis of the E. sinensis. The E. sinensis (female: 150 crabs and male: 150 crabs; 67.11 ± 0.76 g) were randomly divided into the control group (Foot injection with 200 μl PBS) and infection group (Foot injection with 200 μl A. hydrophila of 10⁶ cfu/mL). The hepatopancreas and serum was collected to detect the related indicators after injection 24 h. The results showed that A. hydrophila significantly reduced the malondialdehyde (MDA) level and gamma-glutamyl-cysteine synthetase (γ-GCS) activity in the hepatopancreas of male and female crabs (P < 0.05). A.hydrophila also significantly decreased the total-superoxide dismutase (T-SOD) activity while the levels of total antioxidant capacity (T-AOC) and total glutathione (T-GSH) were significantly increased in the hepatopancreas and serum of male crabs (P < 0.05). At the transcriptional level, the expression of catalase (CAT) and glutathione peroxidases (GPx), Glutathione S-transferase (GST) in the hepatopancreas of male and female crabs was significantly reduced compared to the control group (P < 0.05). However, A. hydrophila could not significantly change the Kelch-like ECH-associated protein 1 (Keap1) gene expression level in both of male and female carbs. A. hydrophila injection for 24 h, the lysozyme (LZM) and phenoloxidase (PO) activity was significantly increased in the hepatopancreas and serum of the male and female crabs (P < 0.05). Simultaneous increase of immune-related enzyme activity (acid phosphatase and alkaline phosphatase) was found in the serum of male and female crabs (P < 0.05). However, the acid phosphatase (ACP) and alkaline phosphatase (ALP) activity was significantly decreased in the hepatopancreas of male and female crabs (P < 0.05). Meanwhile, the LZM mRNA level was significantly decreased in the hepatopancreas of E. sinensis (P < 0.05). Furthermore, A. hydrophila significantly inhibited the mRNA expression of immune regulated factors (Interleukin enhancer binding factor 2: ILF2, interleukin-16: IL-16, Toll-like receptor: TLR) in the male and female crabs. The levels of inflammatory cytokines (interleukin-1β: IL-1β, interleukin-6: IL-6, interleukin-8: IL-8, interleukin-10: IL-10) were significantly increased in the hepatopancreas of male and female crabs. Moreover, A.hydrophila increased the mRNA expression of apoptosis - related genes in male crabs (p38 mitogen-activated protein kinase: p38, adamalysin 17: ADAM17, Cysteine-aspartic acid protease 3: Caspase 3, and Bcl-2-associated X: BAX), but reduced the expression of p38, ADAM17, Caspase 3 and BAX genes in female crabs. In conclusion, A. hydrophila could induce oxidative stress and the response of inflammation and immunity, and also trigger the mRNA expression changes of apoptosis related-genes in E. sinensis. This study provides a theoretical basis for the study of E. sinensis diseases.

Vitellogenin 2 promotes muscle development and stimulates the browning of white fat

Eggs are rich in nutrients and contain a lot of protein. Although eggs have proved to accelerate the growth of C2C12 cells, the regulatory and mechanism of fertilized egg yolk extract (FEYE) on skeletal muscle development and fat metabolism remains unclearly. The mice were treated with FEYE by gavage for 24 d, we found that FEYE can inhibit the expression of skeletal muscle atrophy genes such as MSTN and Murf-1, and up-regulate the expression levels of MYOD, MYOG and Irisin. In addition, the treatment of FEYE induced UCP1 and PGC1α high expression in WAT, thereby causing WAT browning reaction. In order to confirm the composition of FEYE, we performed protein full spectrum identification (LC MS/MS) analysis and found the most enriched component is vitellogenin 2 (VTG2). Therefore, we added the recombinant protein VTG2 to C2C12 cells and found that VTG2 promoted the proliferation and differentiation of C2C12 cells. After that, we further proved that VTG2 inhibited the expression of MSTN and improved the expression of MYOD and Irisin. Finally, the dual luciferase test proved that VTG2 directly inhibited the transcriptional activity of MSTN. Our results conclude that FEYE inhibits the expression of MSTN in muscle tissues by delivering VTG2, thereby promoting skeletal muscle development, and can also promote the expression level of FNDC5 in serum. Then, FNDC5 acts on the fat through the serum, stimulating the browning reaction of white adipocytes. Therefore, VTG2 can be used to stop muscle consumption, improve skeletal muscle aging, and prevent obesity.

Comparative Analysis of PacBio and Oxford Nanopore Sequencing Technologies for Transcriptomic Landscape Identification of Penaeus monodon

With the advantages that long-read sequencing platforms such as Pacific Biosciences (Menlo Park, CA, USA) (PacBio) and Oxford Nanopore Technologies (Oxford, UK) (ONT) can offer, various research fields such as genomics and transcriptomics can exploit their benefits. Selecting an appropriate sequencing platform is undoubtedly crucial for the success of the research outcome, thus there is a need to compare these long-read sequencing platforms and evaluate them for specific research questions. This study aims to compare the performance of PacBio and ONT platforms for transcriptomic analysis by utilizing transcriptome data from three different tissues (hepatopancreas, intestine, and gonads) of the juvenile black tiger shrimp, Penaeus monodon. We compared three important features: (i) main characteristics of the sequencing libraries and their alignment with the reference genome, (ii) transcript assembly features and isoform identification, and (iii) correlation of the quantification of gene expression levels for both platforms. Our analyses suggest that read-length bias and differences in sequencing throughput are highly influential factors when using long reads in transcriptome studies. These comparisons can provide a guideline when designing a transcriptome study utilizing these two long-read sequencing technologies.

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.

De novo transcriptome assembly of four organs of Collichthys lucidus and identification of genes involved in sex determination and reproduction

The spinyhead croaker (Collichthys lucidus) is a commercially important fish species, which is mainly distributed in the coastal regions of China. However, little is known about the molecular regulatory mechanism underlying reproduction in C. lucidus. A de novo transcriptome assembly in brain, liver, ovary and testis tissues of C. lucidus was performed. Illumina sequencing generated 60,322,004, 57,044,284, 60,867,978 and 57,087,688 clean reads from brain, liver, ovary and testis tissues of C. lucidus, respectively. Totally, 131,168 unigenes with an average length of 644 bp and an N50 value of 1033 bp were assembled. In addition, 1288 genes were differentially expressed between ovary and testis, including 442 up-regulated and 846 down-regulated in ovary. Functional analysis revealed that the differentially expressed genes between ovary and testis were mainly involved in the function of sexual reproduction, sex differentiation, development of primary male sexual characteristics, female gamete generation, and male sex differentiation. A number of genes which might be involved in the regulation of reproduction and sex determination were found, including HYAL and SYCP3 and BMP15. Furthermore, 35,476 simple sequence repeats (SSRs) were identified in this transcriptome dataset, which would contribute to further genetic and mechanism researches. De novo transcriptome sequencing analysis of four organs of C. lucidus provides rich resources for understanding the mechanism of reproductive development of C. lucidus and further investigation of the molecular regulation of sex determination and reproduction of C. lucidus.

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

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

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.

Semi-aquatic spider silks: transcripts, proteins, and silk fibres of the fishing spider, Dolomedes triton (Pisauridae)

To survive in terrestrial and aquatic environments, spiders often rely heavily on their silk. The vast majority of silks that have been studied are from orb-web or cob-web weaving species, leaving the silks of water-associated spiders largely undescribed. We characterize transcripts, proteins, and silk fibres from the semi-aquatic spider Dolomedes triton. From silk gland RNAseq libraries, we report 18 silk transcripts representing four categories of known silk protein types: aciniform, ampullate, pyriform, and tubuliform. Proteomic and structural analyses (scanning electron microscopy, energy dispersive X-ray spectrometry, contact angle) of the D. triton submersible egg sac reveal similarities to silks from aquatic caddisfly larvae. We identified two layers in D. triton egg sacs, notably a highly hydrophobic outer layer with a different elemental composition compared to egg sacs of terrestrial spiders. These features may provide D. triton egg sacs with their water repellent properties.

Transcriptome analysis for identification of candidate genes related to sex determination and growth in Charybdis japonica

Charybdis japonica is an important cultured crab in China and it exhibits sex differences in their growth. Growth is an important economic trait that is controlled by many genes. In order to discover the growth-related regulatory mechanisms, whole-body transcriptomic sequencing and comparative analyses in different genders of C. japonica were conducted based on Illumina RNA-seq technology. In the present study, we obtained 62,879,204 and 60,226,334 clean reads in female and male libraries, respectively. 25,000,000 clean reads of every library were randomly selected and compared with Nt database to examine the possible contamination. Results showed that all clean reads were distributed among C. japonica or other species that were closely relative to this species, indicating no-pollution. De novo assembly was performed and a total of 32,543 and 44,174 unigenes were produced in female and male of C. japonica, respectively. Among all the unigenes, 12,591 and 14,455 unigenes of female and male crabs were annotated based on protein databases. Moreover, a total of 33,926 unigenes were found to contain ORFs and 52,839 SSRs were detected. The contrast between male and female C. japonica identifying 1939 unigenes were significantly differentially expressed. In addition, we specifically discussed some gene functions and pathways potentially associated with sex determination and growth. This is the first systematic report of whole transcriptome in C. japonica. The transcriptome information provides a basic resource for further studies on understanding the molecular basis of biological processes in C. japonica and other crustaceans.

Crustacean Female Sex Hormone From the Mud Crab Scylla paramamosain Is Highly Expressed in Prepubertal Males and Inhibits the Development of Androgenic Gland

Recently, the crustacean female sex hormone (CFSH), which is considered a female-specific hormone, has been shown to play a crucial role in female phenotypes in crustaceans. In this study, two transcripts (Sp-CFSH1 and Sp-CFSH2) encoding the same CFSH precursor were cloned from the mud crab Scylla paramamosain. Homology and phylogenetic analysis showed that CFSHs were homologous to interleukin-17 and highly conserved among brachyuran crabs. PCR analysis revealed that Sp-CFSH was expressed exclusively in the eyestalk ganglion of both prepubertal males and females, and surprisingly, the abundance of Sp-CFSH transcripts detected in the males were not significantly different from that of the females (P > 0.05). In addition, mRNA in situ hybridization showed that Sp-CFSH was localized in the X-organ of the male eyestalk ganglion. During the development of the androgenic gland (AG), the level of Sp-IAG mRNA in AG remained at low levels from stages I to II (early stage) but had a significant increase at stage III (mature stage). In contrast, the level of Sp-CFSH transcripts in the eyestalk ganglion was high in the early stage but extremely low in the mature stage. To investigate the potential function of CFSH in male S. paramamosain, the recombinant protein (∼20 kDa) was expressed in Escherichia coli and was subsequently added to AG explants in vitro. It was demonstrated that recombinant Sp-CFSH protein significantly reduced the expression of Sp-IAG in the AG explants at a concentration of 10-6 M (P < 0.05). In conclusion, our study provides the first piece of evidence that shows CFSH from the eyestalk ganglion acts as a negative regulator inhibiting the development of AG in crustaceans.

Transcriptome-seq provides insights into sex-preference pattern of gene expression between testis and ovary of the crucifix crab (Charybdis feriatus)

The crucifix crab, Charybdis feriatus, which mainly inhabits Indo-Pacific region, is regarded as one of the most high-potential species for domestication and incorporation into the aquaculture sector. However, the regulatory mechanisms of sex determination and differentiation of this species remain unclear. To identify candidate genes involved in sex determination and differentiation, high throughput sequencing of transcriptome from the testis and ovary of C. feriatus was performed by the Illumina platform. After removing adaptor primers, low-quality sequences and very short (<50 nt) reads, we obtained 80.9 million and 66.2 million clean reads from testis and ovary, respectively. A total of 86,433 unigenes were assembled, and ~43% (37,500 unigenes) were successfully annotated to the NR, NT, Swiss-Prot, KEGG, COG, GO databases. By comparing the testis and ovary libraries, we obtained 27,636 differentially expressed genes. Some candidate genes involved in the sex determination and differentiation of C. feriatus were identified, such as vasa, pgds, vgr, hsp90, dsx-f, fem-1, and gpr. In addition, 88,608 simple sequence repeats were obtained, and 61,929 and 77,473 single nucleotide polymorphisms from testis and ovary were detected, respectively. The transcriptome profiling was validated by quantitative real-time PCR in 30 selected genes, which showed a good consistency. The present study is the first high-throughput transcriptome sequencing of C. feriatus. These findings will be useful for future functional analysis of sex-associated genes and molecular marker-assisted selections in C. feriatus.

Proteome and Transcriptome Analysis of Ovary, Intersex Gonads, and Testis Reveals Potential Key Sex Reversal/Differentiation Genes and Mechanism in Scallop Chlamys nobilis

Bivalve mollusks exhibit hermaphroditism and sex reversal/differentiation. Studies generally focus on transcriptional profiling and specific genes related to sex determination and differentiation. Few studies on sex reversal/differentiation have been reported. A combination analysis of gonad proteomics and transcriptomics was conducted on Chlamys nobilis to provide a systematic understanding of sex reversal/differentiation in bivalves. We obtained 4258 unique peptides and 93,731 unigenes with good correlation between messenger RNA and protein levels. Candidate genes in sex reversal/differentiation were found: 15 genes differentially expressed between sexes were identified and 12 had obvious sexual functions. Three novel genes (foxl2, β-catenin, and sry) were expressed highly in intersex individuals and were likely involved in the control of gonadal sex in C. nobilis. High expression of foxl2 or β-catenin may inhibit sry and activate 5-HT receptor and vitellogenin to maintain female development. High expression of sry may inhibit foxl2 and β-catenin and activate dmrt2, fem-1, sfp2, sa6, Amy-1, APCP4, and PLK to maintain male function. High expression of sry, foxl2, and β-catenin in C. nobilis may be involved in promoting and maintaining sex reversal/differentiation. The downstream regulator may not be dimorphic expressed genes, but genes expressed in intersex individuals, males and females. Different expression patterns of sex-related genes and gonadal histological characteristics suggested that C. nobilis may change its sex from male to female. These findings suggest highly conserved sex reversal/differentiation with diverged regulatory pathways during C. nobilis evolution. This study provides valuable genetic resources for understanding sex reversal/differentiation (intersex) mechanisms and pathways underlying bivalve reproductive regulation.

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

The molecular mechanism underlying sex determination and gonadal differentiation of the mud crab (Scylla paramamosain) has received considerable attention, due to the remarkably biological and economic differences between sexes. However, sex-biased genes, especially non-coding genes, which account for these differences, remain elusive in this crustacean species. In this study, the first de novo gonad transcriptome sequencing was performed to identify both differentially expressed genes and long non-coding RNAs (lncRNAs) between male and female S. paramamosain by using Illumina Hiseq2500. A total of 79,282,758 and 79,854,234 reads were generated from ovarian and testicular cDNA libraries, respectively. After filtrating and de novo assembly, 262,688 unigenes were produced from both libraries. Of these unigenes, 41,125 were annotated with known protein sequences in public databases. Homologous genes involved in sex determination and gonadal development pathways (Sxl-Tra/Tra-2-Dsx/Fru, Wnt4, thyroid hormone synthesis pathway, etc.) were identified. Three hundred and sixteen differentially expressed unigenes were further identified between both transcriptomes. Meanwhile, a total of 233,078 putative lncRNAs were predicted. Of these lncRNAs, 147 were differentially expressed between sexes. qRT-PCR results showed that nine lncRNAs negatively regulated the expression of eight genes, suggesting a potential role in sex differentiation. These findings will provide fundamental resources for further investigation on sex differentiation and regulatory mechanism in crustaceans.

De novo transcriptome sequencing and analysis of male and female swimming crab (Portunus trituberculatus) reproductive systems during mating embrace (stage II)

Background

The swimming crab Portunus trituberculatus is one of the most commonly farmed crustaceans in China. As one of the most widely known and high-value edible crabs, it crab supports large crab fishery and aquaculture in China. Only large and sexually mature crabs can provide the greatest economic benefits, suggesting the considerable effect of reproductive system development on fishery. Studies are rarely conducted on the molecular regulatory mechanism underlying the development of the reproductive system during the mating embrace stage in this species. In this study, we used high-throughput sequencing to sequence all transcriptomes of the P. trituberculatus reproductive system.

Results

Transcriptome sequencing of the reproductive system produced 81,688,878 raw reads (38,801,152 and 42,887,726 reads from female and male crabs, respectively). Low-quality (quality <20) reads were trimmed and removed, leaving only high-quality reads (37,020,664 and 41,021,030 from female and male crabs, respectively). A total of 126,188 (female) and 164,616 (male) transcripts were then generated by de novo transcriptome assembly using Trinity. Functional annotation of the obtained unigenes revealed that a large number of key genes and some important pathways may participate in cell proliferation and signal transduction. On the basis of our transcriptome analyses and as confirmed by quantitative real-time PCR, a number of genes potentially involved in the regulation of gonadal development and reproduction of P. trituberculatus were identified: ADRA1B, BAP1, ARL3, and TRPA1.

Conclusion

This study is the first to report on the whole reproductive system transcriptome information in stage II of P. trituberculatus gonadal development and provides rich resources for further studies to elucidate the molecular basis of the development of reproductive systems and reproduction in crabs. The current study can be used to further investigate functional genomics in this species.

Electronic supplementary material

The online version of this article (10.1186/s12863-017-0592-5) contains supplementary material, which is available to authorized users.

Comparative Transcriptome Analysis Reveals Related Regulatory Mechanisms of Androgenic Gland in Eriocheir sinensis

Chinese mitten crab (Eriocheir sinensis) is one of the most commercially important aquaculture species in China. The androgenic gland (AG) of crustaceans plays pivotal roles in the regulation of male differentiation and in maintaining the male sexual characteristics. In order to reveal related mechanisms in AG, we compared transcriptomes of AG between proliferation and secretion phase. A total of 72,000 unigenes and 4,027 differentially expressed genes were obtained. Gene ontology enrichment analysis indicated that biological processes and metabolic pathways related to protein synthesis and secretion such as transcription, translation, and signal transduction were significantly enriched. Critical genes such as IAG, SXL, TRA-2, SRY, FTZ-F1, FOXL2, and FEM-1 were identified and potentially involved in maintaining the testis development and spermatogenesis. Ribosomes pathway revealed the cause of insulin-like androgenic gland hormone secretion increase. Three insulin-like receptors were thought to be associated with growth and spermatogenesis. In the neuroactive ligand-receptor interaction pathway, the expression of octopamine receptor, 5-HT receptor 1, and melatonin receptor was significantly changed, which revealed the key regulation mechanism of aggressive and mating behavior of males. Comparative transcriptome analysis provided new insights into the genome-wide molecular mechanisms of AG development and the regulatory mechanisms of male development.

Sex and tissue specific gene expression patterns identified following de novo transcriptomic analysis of the Norway lobster, Nephrops norvegicus

BACKGROUND

The Norway lobster, Nephrops norvegicus, is economically important in European fisheries and is a key organism in local marine ecosystems. Despite multi-faceted scientific interest in this species, our current knowledge of genetic resources in this species remains very limited. Here, we generated a reference de novo transcriptome for N. norvegicus from multiple tissues in both sexes. Bioinformatic analyses were conducted to detect transcripts that were expressed exclusively in either males or females. Patterns were validated via RT-PCR.

RESULTS

Sixteen N. norvegicus libraries were sequenced from immature and mature ovary, testis and vas deferens (including the masculinizing androgenic gland). In addition, eyestalk, brain, thoracic ganglia and hepatopancreas tissues were screened in males and both immature and mature females. RNA-Sequencing resulted in >600 million reads. De novo assembly that combined the current dataset with two previously published libraries from eyestalk tissue, yielded a reference transcriptome of 333,225 transcripts with an average size of 708 base pairs (bp), with an N50 of 1272 bp. Sex-specific transcripts were detected primarily in gonads followed by hepatopancreas, brain, thoracic ganglia, and eyestalk, respectively. Candidate transcripts that were expressed exclusively either in males or females were highlighted and the 10 most abundant ones were validated via RT-PCR. Among the most highly expressed genes were Serine threonine protein kinase in testis and Vitellogenin in female hepatopancreas. These results align closely with gene annotation results. Moreover, a differential expression heatmap showed that the majority of differentially expressed transcripts were identified in gonad and eyestalk tissues. Results indicate that sex-specific gene expression patterns in Norway lobster are controlled by differences in gene regulation pattern between males and females in somatic tissues.

CONCLUSIONS

The current study presents the first multi-tissue reference transcriptome for the Norway lobster that can be applied to future biological, wild restocking and fisheries studies. Sex-specific markers were mainly expressed in males implying that males may experience stronger selection than females. It is apparent that differential expression is due to sex-specific gene regulatory pathways that are present in somatic tissues and not from effects of genes located on heterogametic sex chromosomes. The N. norvegicus data provide a foundation for future gene-based reproductive studies.

Alternative Splicing Profile and Sex-Preferential Gene Expression in the Female and Male Pacific Abalone Haliotis discus hannai

In order to characterize the female or male transcriptome of the Pacific abalone and further increase genomic resources, we sequenced the mRNA of full-length complementary DNA (cDNA) libraries derived from pooled tissues of female and male Haliotis discus hannai by employing the Iso-Seq protocol of the PacBio RSII platform. We successfully assembled whole full-length cDNA sequences and constructed a transcriptome database that included isoform information. After clustering, a total of 15,110 and 12,145 genes that coded for proteins were identified in female and male abalones, respectively. A total of 13,057 putative orthologs were retained from each transcriptome in abalones. Overall Gene Ontology terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways analyzed in each database showed a similar composition between sexes. In addition, a total of 519 and 391 isoforms were genome-widely identified with at least two isoforms from female and male transcriptome databases. We found that the number of isoforms and their alternatively spliced patterns are variable and sex-dependent. This information represents the first significant contribution to sex-preferential genomic resources of the Pacific abalone. The availability of whole female and male transcriptome database and their isoform information will be useful to improve our understanding of molecular responses and also for the analysis of population dynamics in the Pacific abalone.

Vitellogenin2: spermatozoon specificity and immunoprotection in mud crabs

As the precursor of vitellin (Vn), vitellogenin (Vg) has initially been considered as a female-specific protein involved in vitellogenesis, while it was also present in males induced by hormones or organs manipulation. Distinct from vtg1 we previously found in female mud crab Scylla paramamosain, vtg2 was intriguingly detected in male testis under normal physiological conditions in this study. Sequence analysis showed that vtg2 and vtg1 were actually two isoforms of Vg caused by different types of alternative splicing. PCR and in situ hybridization analysis revealed that vtg2 was localized only in the testicular spermatozoa, while Vn was detected in both the spermatozoa of the testis and seminal vesicle. Therefore, we speculated that Vn was initially translated in testicular spermatozoa, then migrated with spermatozoa, and finally stored in the seminal vesicle, where spermatozoa gradually accomplished maturation. We presumed that vtg2/Vn might act as an immune-relevant molecule in the male reproduction system. In the subsequent experiment, the expression of vtg2/Vn in testis was significantly induced in response to lipopolysaccharide (LPS) and lipoteichoic acid (LTA) injection at both transcriptional and translational levels. In the light of the results presented above, we deemed that vtg2/Vn is a novel candidate of immune-relevant molecules involved in immunoprotection during the spermatozoon maturation, and this research helps to open a new avenue for further exploring the role of Vg.

Applying the Power of Transcriptomics: Understanding Male Sexual Development in Decapod Crustacea

The decapod Crustacea are the most species-rich order of the Crustacea and include some of the most charismatic and highly valued commercial species. Thus the decapods draw a significant research interest in relation to aquaculture, as well as gaining a broader understanding of these species' biology. However, the diverse physiology of the group considered with the lack of a model species have presented an obstacle for comparative analyses. In reflection of this, the recent integration of comparative transcriptomics has rapidly advanced our understanding of key regulatory pathways and developmental phenomena, an example being our understanding of sexual development. We discuss our work in the Eastern spiny lobster, Sagmariasus verreauxi, in the context of what is currently known about male sexual development in the decapods, highlighting the importance of transcriptomic techniques in achieving our recent advancements. We describe the progression made in our understanding of male sexual differentiation and maturation, as mediated by the insulin-like androgenic gland hormone (IAG), integrating the role of regulatory binding proteins (IGFBPs), a tyrosine kinase insulin receptor (TKIR), as well as the upstream effect of neuroendocrine hormones (GIH and MIH). We then consider the less well understood mechanism of male sex determination, with an emphasis on what we believe to be the key regulatory factors, the Dsx- and mab-3-related transcription factors (Dmrts). Finally, we discuss the function of the antennal gland (AnG) in sexual development, relating to the emergence of male-biased upregulation in the AnG in later sexual maturation and the sexually dimorphic expression of two key genes Sv-TKIR and Sv-Dmrt1 We then present the AnG as a case study to illustrate how comparative transcriptomic techniques can be applied to guide preliminary analyses, like the hypothesis that the AnG may function in pheromone biosynthesis. In summary, we describe the power of transcriptomics in facilitating the progress made in our understanding of male sexual development, as illustrated by the commercial decapod species, S. verreauxi Considering future directions, we suggest that the integration of multiple omics-based techniques offers the most powerful tool to ensure we continue to piece together the biology of the important group of decapod Crustacea.

Similarities between decapod and insect neuropeptidomes

Background. Neuropeptides are important regulators of physiological processes and behavior. Although they tend to be generally well conserved, recent results using trancriptome sequencing on decapod crustaceans give the impression of significant differences between species, raising the question whether such differences are real or artefacts. Methods. The BLAST+ program was used to find short reads coding neuropeptides and neurohormons in publicly available short read archives. Such reads were then used to find similar reads in the same archives, and the DNA assembly program Trinity was employed to construct contigs encoding the neuropeptide precursors as completely as possible. Results. The seven decapod species analyzed in this fashion, the crabs Eriocheir sinensis, Carcinus maenas and Scylla paramamosain, the shrimp Litopenaeus vannamei, the lobster Homarus americanus, the fresh water prawn Macrobrachium rosenbergii and the crayfish Procambarus clarkii had remarkably similar neuropeptidomes. Although some neuropeptide precursors could not be assembled, in many cases individual reads pertaining to the missing precursors show unambiguously that these neuropeptides are present in these species. In other cases, the tissues that express those neuropeptides were not used in the construction of the cDNA libraries. One novel neuropeptide was identified: elongated PDH (pigment dispersing hormone), a variation on PDH that has a two-amino-acid insertion in its core sequence. Hyrg is another peptide that is ubiquitously present in decapods and is likely a novel neuropeptide precursor. Discussion. Many insect species have lost one or more neuropeptide genes, but apart from elongated PDH and hyrg all other decapod neuropeptides are present in at least some insect species, and allatotropin is the only insect neuropeptide missing from decapods. This strong similarity between insect and decapod neuropeptidomes makes it possible to predict the receptors for decapod neuropeptides that have been deorphanized in insects. This includes the androgenic insulin-like peptide that seems to be homologous to drosophila insulin-like peptide 8.

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.