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

Digital image processing: A new tool for morphological measurements of freshwater turtles under rehabilitation

Freshwater fauna is facing an uphill task for survival in the Ganga Basin, India, due to a range of factors causing habitat degradation and fragmentation, necessitating conservation interventions. As part of the ongoing efforts to conserve the freshwater fauna of the Basin, we are working on rehabilitating rescued freshwater chelonians. We carry out various interventions to restore rescued individuals to an apparent state of fitness for their release in suitable natural habitats. Morphometric measurements are crucial to managing captive wild animals for assessing their growth and well-being. Measurements are made using manual methods like vernier caliper that are prone to observer error experience and require handling the specimens for extended periods. Digital imaging technology is rapidly progressing at a fast pace and with the advancement of technology. We acquired images of turtles using smartphones along with manual morphometric measurements using vernier calipers of the straight carapace length and straight carapace width. The images were subsequently processed using ImageJ, a freeware and compared with manual morphometric measurements. A significant decrease in the time spent in carrying out morphometric measurements was observed in our study. The difference in error in measurements was, however, not significant. A probable cause for this may have been the extensive experience of the personnel carrying out the measurements using vernier caliper. Digital image processing technology can cause a significant reduction in the stress of the animals exposed to handling during measurements, thereby improving their welfare. Additionally, this can be used in the field to carry out morphometric measurements of free-ranging individuals, where it is often difficult to capture individuals, and challenges are faced in obtaining permission to capture specimens.

Identification and functional analysis of the doublesex gene in the mud crab Scylla paramamosain

Doublesex (Dsx) is a crucial member of the Dmrt gene family and plays a vital role in sex determination and differentiation among the animal kingdom. In the present study, a doublesex (designated as Spdsx) gene was identified and characterized for the first time in the mud crab, Scylla paramamosain. The Spdsx cDNA contains an 801 bp open reading frame (ORF) encoding 266 amino acids with a conserved DM domain. Meanwhile, to elucidate the conservation of Dsx, its orthologus were identified in several crustacean species as well. In addition, the expression pattern of Spdsx in various adult tissues and during embryo development was analyzed with qRT-PCR technology. Finally, the roles of Spdsx might play in the testis, androgenic gland, and ovary were analyzed by RNAi technology. The main results are as follows: (1) the Spdsx gene widely existed in analyzed crustacean species, and the multiple sequences alignment result indicated the conservation of Dsx was low except for the DM domain; (2) only one dsx gene was identified in analyzed crab and lobster, while 2 dsx genes (dsx-1 and dsx-2) were identified in shrimps; (3) the Spdsx gene was widely expressed in analyzed tissues, and the expression level in androgenic gland was obviously higher than that in other tissues. Interestingly, the expression level of Spdsx in the ovary was significantly higher than that in testis (p < 0.05); (4) The expression pattern of Spdsx during embryo development was divided into two groups: remained stable from blastula stage to 5 pairs of appendages stage; after 5 pairs of appendages stage, the expression level increased and remained stable from 7 pairs of appendages stage to hatching stage; (5) After the silencing of Spdsx, the expression level of marker genes in testis, ovary, and androgenic gland significantly changed, among which the expression level of vtg and vtgR in ovary down-regulated, the dmrt-like and dmrt-1a (exclusively expressed in testis) in testis up-regulated and the IAG in androgenic gland down-regulated. All the results above demonstrated that the Spdsx play crucial roles in regulating the reproduction system development of mud crab.

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 Gonadal Transcriptomes Uncovers Reproduction-Related Genes with Sexually Dimorphic Expression Patterns in Diodon hystrix

Diodon hystrix is a new and emerging aquaculture species in south China. However, due to the lack of understanding of reproductive regulation, the management of breeding and reproduction under captivity remains a barrier for the commercial aquaculture of D. hystrix. More genetic information is needed to identify genes critical for gonadal development. Here, the first gonadal transcriptomes of D. hystrix were analyzed and 151.89 million clean reads were generated. All reads were assembled into 57,077 unigenes, and 24,574 could be annotated. By comparing the gonad transcriptomes, 11,487 differentially expressed genes were obtained, of which 4599 were upregulated and 6888 were downregulated in the ovaries. Using enrichment analyses, many functional pathways were found to be associated with reproduction regulation. A set of sex-biased genes putatively involved in gonad development and gametogenesis were identified and their sexually dimorphic expression patterns were characterized. The detailed transcriptomic data provide a useful resource for further research on D. hystrix reproductive manipulation.

Sex Determination and Differentiation in Decapod and Cladoceran Crustaceans: An Overview of Endocrine Regulation

Mechanisms underlying sex determination and differentiation in animals are known to encompass a diverse array of molecular clues. Recent innovations in high-throughput sequencing and mass spectrometry technologies have been widely applied in non-model organisms without reference genomes. Crustaceans are no exception. They are particularly diverse among the Arthropoda and contain a wide variety of commercially important fishery species such as shrimps, lobsters and crabs (Order Decapoda), and keystone species of aquatic ecosystems such as water fleas (Order Branchiopoda). In terms of decapod sex determination and differentiation, previous approaches have attempted to elucidate their molecular components, to establish mono-sex breeding technology. Here, we overview reports describing the physiological functions of sex hormones regulating masculinization and feminization, and gene discovery by transcriptomics in decapod species. Moreover, this review summarizes the recent progresses of studies on the juvenile hormone-driven sex determination system of the branchiopod genus Daphnia, and then compares sex determination and endocrine systems between decapods and branchiopods. This review provides not only substantial insights for aquaculture research, but also the opportunity to re-organize the current and future trends of this field.

Application of quantitative transcriptomics in evaluating the ex vivo effects of per- and polyfluoroalkyl substances on Atlantic cod (Gadus morhua) ovarian physiology

Because of their global consumption and persistence, per- and polyfluoroalkyl substances (PFASs), are ubiquitously distributed in the environment, as well as in wildlife and humans. In the present study, we have employed an ex vivo organ culture technique, based on the floating agarose method, of Atlantic cod ovarian tissue to investigate the effects of three different concentrations of PFOS, PFOA (1, 5 and 25 μM) and PFNA (0.5, 5 and 50 μM), used singly and in also in combination (1×, 20× and 100×). In the 1× exposure mixture, concentrations were decided based on their proportional levels (in molar equivalents) relative to PFOS, which is the most abundant PFAS in cod liver from a 2013 screening project. To investigate the detailed underlying mechanisms and biological processes, transcriptome sequencing was performed on exposed ovarian tissue. The number of differentially expressed genes (DEGs) having at least 0.75 log₂-fold change was elevated in high, compared to low and medium concentration exposures. The highest PFNA, PFOA and PFOS concentrations, and the highest (100×) mixture exposure, showed 40, 68, 1295, and 802 DEGs, respectively. The latter two exposure groups shared a maximum of 438 DEGs. In addition, they both shared the majority of functionally enriched pathways belonging to biological processes such as cellular signaling, cell adhesion, lipid metabolism, immunological responses, cancer, reproduction and metabolism. Shortlisted DEGs that were specifically annotated to reproduction associated gene ontology (GO) terms were observed only in the highest PFOS and mixture exposure groups. These transcripts contributed to ovarian key events such as steroidogenesis (star, cyp19a1a), oocyte growth (amh), maturation (igfbp5b, tgfβ2, tgfβ3), and ovulation (pgr, mmp2). Contrary to other PFAS congeners, the highest PFOS concentration showed almost similar transcript expression patterns compared to the highest mixture exposure group. This indicates that PFOS is the active component of the mixture that significantly altered the normal functioning of female gonads, and possibly leading to serious reproductive consequences in teleosts.

Molecular cloning, characterization, and expression analysis of a sex-biased transcriptional factor sox9 gene of mud crab Scylla paramamosain

Sox9 gene, a crucial member of the Sox gene family, is present in various organisms and involved in many physiological processes, especially in sex determination and gonad development. In this study, we reported a sox9 gene (designated as Spsox9) from Scylla paramamosain through analyzing published gonad transcriptome data. Meanwhile, the accuracy was validated by PCR technology, and the 3' sequences were cloned with 3' RACE technology. The full-length cDNA of Spsox9 is 2843 bp, consisting of a 243 bp 5' UTR, an 1124 bp 3' UTR, and a 1476 bp ORF encoding 491 amino acids. Furthermore, to better understand its conservation among crustacean species, the sox9 gene ortholog was identified in several other crustaceans species with their published transcriptome data, respectively. All of the Sox9 proteins identified in the current study had the common feature of Sox proteins (HMG domain) and were highly conserved among analyzed crustacean species. In all examined tissues, the Spsox9 was mainly expressed in the gonad (testis and ovary), eyestalk, and cerebral ganglion. During embryo development, Spsox9 was highly expressed in 5 pairs of appendages, 7 pairs of appendages, and eye-pigment formation stage. During ovary development, the expression level of Spsox9 remained stable in the first 4 stages (O1-O4) and decreased in the tertiary vitellogenesis (O5) stage. During testis development, the expression level of Spsox9 was highest in the spermatid stage (T2) and was significantly different from that in the spermatocyte stage (T1) and mature sperm stage (T3) (p < 0.05). In addition, Spsox9 exhibited a sex-biased expression pattern in T1 and O1. These present results indicated that the Spsox9 gene might play crucial roles in the gonad and embryo development of mud crab.

Transcriptomic changes across vitellogenesis in the black tiger prawn (Penaeus monodon), neuropeptides and G protein-coupled receptors repertoire curation

The black tiger prawn (Penaeus monodon) is one of the most commercially important prawn species world-wide, yet there are currently key issues that hinder aquaculture of this species, such as low spawning capacity of captive-reared broodstock females and lack of globally available fully domesticated strains. In this study, we analysed the molecular changes that occur from vitellogenesis to spawning of a fully domesticated population of P.monodon (Madagascar) using four tissues [brain and thoracic ganglia (central nervous system - CNS), eyestalks, antennal gland, and ovary] highlighting differentially expressed genes that could be involved in the sexual maturation. In addition, due to their key role in regulating multiple physiological processes including reproduction, transcripts encoding P.monodon neuropeptides and G protein-coupled receptors (GPCRs) were identified and their expression pattern was assessed. A few neuropeptides and their putative GPCRs which were previously implicated in reproduction are discussed. We identified 573 differentially expressed transcripts between previtellogenic and vitellogenic stages, across the four analysed tissues. Multiple transcripts that have been linked to ovarian maturation were highlighted throughout the study, these include vitellogenin, Wnt, heat shock protein 21, heat shock protein 90, teneurin, Fs(1)M3, hemolymph clottable proteins and some other candidates. Seventy neuropeptide transcripts were also characterized from our de novo assembly. In addition, a hybrid approach that involved clustering and phylogenetics analysis was used to annotate all P. monodon GPCRs, revealing 223 Rhodopsin, 100 Secretin and 27 Metabotropic glutamate GPCRs. Given the key commercial significance of P.monodon and the industry requirements for developing better genomic tools to control reproduction in this species, our findings provide a foundation for future gene-based studies, setting the scene for developing innovative tools for reproduction and/or sexual maturation control in P. monodon.

Characterization and antimicrobial evaluation of a new Spgly-AMP, glycine-rich antimicrobial peptide from the mud crab Scylla paramamosain

Antimicrobial peptide (AMP) is a crucial component of the innate immune system in crustaceans. In mud crab, Scylla paramamosain, a commercially important species, a glycine-rich antimicrobial peptide (Spgly-AMP) gene was newly identified and putatively encoded a 26aa signal peptide and 37aa mature peptide. To understand the function of Spgly-AMP, the expression profile of Spgly-amp gene was characterized, which showed Spgly-amp was expressed widely in most tissues of adult crabs with the highest expression level in hemocytes. After Vibrio parahaemolyticus, PGN, or Poly I:C stimulations, the expression level of Spgly-amp was significantly up-regulated in the hemocytes. In antimicrobial assays, chemically synthesized Spgly-AMP peptides exhibited strong antibacterial activities against both Gram-positive and Gram-negative bacteria and high thermal stability after high-temperature heating. These findings in the present study verified the importance of the Spgly-AMP in defense of pathogenic bacteria infection in the mud crab and provided a promising candidate of antimicrobial agents in the crab aquaculture.

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.

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.

Seminal Plasma Exosomes: Promising Biomarkers for Identification of Male and Pseudo-Males in Cynoglossus semilaevis

In mammals, small RNAs enclosed in exosomes have been identified as appropriate signatures for disease diagnosis. However, there is limited information on exosomes derived from seminal plasma, and few studies have reported analyzed the composition of exosomes and enclosed small RNAs in fish. The half-smooth tongue sole (Cynoglossus semilaevis) is an economically important fish for aquaculture, and it exhibits sexual dimorphism: the female gender show higher growth rates and larger body sizes than males. Standard karyotype analysis and artificial gynogenesis tests have revealed that this species uses heterogametic sex determination (ZW/ZZ), and so-called sex-reversed pseudo-males exist. In this study, we successfully identified exosomes in the seminal plasma of C. semilaevis; to the best of our knowledge, this is the first report of exosomes in fish seminal plasma. Analysis of the nucleotide composition showed that miRNAs were dominant in the exosomes, and the miRNAs were sequenced and compared to identify signature miRNAs as sexual biomarkers. Moreover, target genes of the signature miRNAs were predicted by sequence matching and annotation. Finally, four miRNAs (dre-miR-141-3P, dre-miR-10d-5p, ssa-miR-27b-3p, and ssa-miR-23a-3p) with significant differential expression in the males and pseudo-males were selected from the signature candidate miRNAs as markers for sex identification, and their expression profiles were verified using real-time quantitative PCR. Our findings could provide an effective detection method for sex differentiation in fish.

Integrated mRNA and miRNA expression profile analyses reveal the potential roles of sex-biased miRNA-mRNA pairs in gonad tissues of the Chinese concave-eared torrent frog (Odorrana tormota)

The Chinese concave-eared torrent frog (Odorrana tormota) is typically sexually dimorphic. Females are significantly less common than males in the wild. Until now, the molecular mechanisms of reproduction and sex differentiation of frogs remain unclear. Here, we integrated mRNA and microRNA (miRNA) expression profiles to reveal the molecular mechanisms of reproduction and sex differentiation in O. tormota. We identified 234 differentially expressed miRNAs (DEMs) and 18,551 differentially expressed transcripts. Of these, 12,053 mRNAs and 64 miRNAs were upregulated in testes, and 6,498 mRNAs and 170 miRNAs were upregulated in ovaries. Integrated analysis of the miRNA and mRNA expression profiles predicted 75,602 potential miRNA-mRNA interaction sites, with 42,065 negative miRNA-mRNA interactions. We found 36 differentially expressed genes (DEGs) related to reproduction and sex differentiation, of which 15 DEGs formed 92 negative miRNA-mRNA interactions with 34 known DEMs. Thus, miRNAs may play other important roles in O. tormota. Furthermore, Gene Ontology enrichment and Kyoto Encyclopedia of Genes and Genomes pathway analyses showed reproductive-related processes, such as the gonadotropinreleasing hormone signaling pathway and ovarian steroidogenesis. Based on functional annotation and the literature, the retinoic acid signaling pathway, the SOX9-AMH pathway, and the process of spermatogenesis may be involved in the molecular mechanisms of reproduction and sex differentiation in O. tormota, and may be regulated by miRNAs. The miRNA-mRNA pairs described may provide further understanding of the regulatory mechanisms associated with reproduction and sex differentiation, and the molecular mechanism of reproduction in O. tormota.

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.

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

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