Characterization of the vasa gene in the Chinese mitten crab Eriocheir sinensis: a germ line molecular marker

Characteristics of the Vasa Gene in Silurus asotus and Its Expression Response to Letrozole Treatment

The identification and expression of germ cells are important for studying sex-related mechanisms in fish. The vasa gene, encoding an ATP-dependent RNA helicase, is recognized as a molecular marker of germ cells and plays a crucial role in germ cell development. Silurus asotus, an important freshwater economic fish species in China, shows significant sex dimorphism with the female growing faster than the male. However, the molecular mechanisms underlying these sex differences especially involving in the vasa gene in this fish remain poorly understood. In this work, the vasa gene sequence of S. asotus (named as Savasa) was obtained through RT-PCR and rapid amplification of cDNA end (RACE), and its expression in embryos and tissues was analyzed using qRT-PCR and an in situ hybridization method. Letrozole (LT) treatment on the larvae fish was also conducted to investigate its influence on the gene. The results revealed that the open reading frame (ORF) of Savasa was 1989 bp, encoding 662 amino acids. The SaVasa protein contains 10 conserved domains unique to the DEAD-box protein family, showing the highest sequence identity of 95.92% with that of Silurus meridionalis. In embryos, Savasa is highly expressed from the two-cell stage to the blastula stage in early embryos, with a gradually decreasing trend from the gastrula stage to the heart-beating stage. Furthermore, Savasa was initially detected at the end of the cleavage furrow during the two-cell stage, later condensing into four symmetrical cell clusters with embryonic development. At the gastrula stage, Savasa-positive cells increased and began to migrate towards the dorsal side of the embryo. In tissues, Savasa is predominantly expressed in the ovaries, with almost no or lower expression in other detected tissues. Moreover, Savasa was expressed in phase I-V oocytes in the ovaries, as well as in spermatogonia and spermatocytes in the testis, implying a specific expression pattern of germ cells. In addition, LT significantly upregulated the expression of Savasa in a concentration-dependent manner during the key gonadal differentiation period of the fish. Notably, at 120 dph after LT treatment, Savasa expression was the lowest in the testis and ovary of the high concentration group. Collectively, findings from gene structure, protein sequence, phylogenetic analysis, RNA expression patterns, and response to LT suggest that Savasa is maternally inherited with conserved features, serving as a potential marker gene for germ cells in S.asotus, and might participate in LT-induced early embryonic development and gonadal development processes of the fish. This would provide a basis for further research on the application of germ cell markers and the molecular mechanisms of sex differences in S. asotus.

Identification and characterization of a new germline-specific marker vasa gene and its promoter in the giant freshwater prawn Macrobrachium rosenbergii

Vasa gene encodes a protein member of DEAD-box superfamily of ATP-dependent RNA helicases, which plays a key role in germline development in metazoans. In present study, we identified a new germline-specific marker Mrvasa in the prawn Macrobrachium rosenbergii, whose genomic DNA sequence consists of 14 exons and 13 introns. A 2516 bp of full-length Mrvasa cDNA encodes a protein of 603 amino acids. It contains nine conserved motifs, a zinc-finger motif, and RGG repeats. RT-PCR indicated that Mrvasa mRNA was specifically expressed in gonads. QPCR analysis further revealed that the expression of Mrvasa mRNA is much higher in testis than in ovary. In testis, the relative expression level of Mrvasa mRNA in late developing stage is significantly higher than that in early-middle developing stage. During ovarian development, no significant difference in expression was found. In situ hybridization demonstrated that Mrvasa mRNA was localized in germline cells including spermatogonia, spermatocytes, and spermatozoa in testes, and previtellogenic and vitellogenic oocytes in ovary. We then isolated the Mrvasa promoter and determined the transcription core region of this promoter. This is the first report on identification of vasa core promoter in crustaceans. Our results will provide a useful germline-specific marker Mrvasa for tracing germline cell formation and development in M. rosenbergii.

Cloning and expression analysis of Shvasa and the molecular regulatory pathways implicated in Cd-induced reproductive toxicity in the freshwater crab Sinopotamon henanense

Cadmium (Cd), a widespread, severely toxic heavy metal, can cause serious reproductive toxicity in animals. However, the molecular pathways associated with Cd-induced effects remain unknown. In this study, we first cloned the vasa gene (Shvasa) and characterized the VASA protein (ShVASA) in Sinopotamon henanense. We then investigated the molecular mechanisms of Cd-induced reproductive toxicity. Shvasa was specifically expressed in the ovary and testis. ShVASA was abundant in early ovarian development and significantly less abundant in mature ovaries. During oogenesis, ShVASA was abundant and evenly distributed in the cytoplasm of the oogonium and previtellogenic oocytes, but gradually accumulated in the nuclear periphery of vitellogenic and mature oocytes. As Cd concentration increased, ShVASA abundance decreased gradually in proliferation-stage ovaries, and increased gradually in mature ovaries. Notably, at the small and large growth stages, ShVASA was upregulated following exposure to 14.5 mg/L Cd and downregulated following exposure to 29 mg/L Cd. In contrast to the unexposed control, ShVASA accumulated around the nuclear periphery in Cd-exposed previtellogenic oocytes and scattered gradually into the cytoplasm in Cd-exposed vitellogenic and mature oocytes. Shvasa RNA interference (RNAi) downregulated Shnanos and Shpiwi, but simultaneous Cd exposure and Shvasa RNAi significantly upregulated Shnanos and downregulated Shpiwi. These data suggested that Cd disrupted Shvasa expression and function, as well as the functions of Shnanos and Shpiwi, leading to severe reproductive toxicity in S. henanense.

Comparative transcriptome reveals the potential modulation mechanisms of estradiol affecting ovarian development of female Portunus trituberculatus

Estradiol is an important sex steroid hormone that is involved in the regulation of crustacean ovarian development. However, the molecular regulatory mechanisms of estradiol on ovarian development are largely unknown. This study performed transcriptome sequencing of ovary, hepatopancreas, brain ganglion, eyestalk, and mandibular organ of crabs after estradiol treatment (0.1μg g^-1 crab weight). A total of 23, 806 genes were annotated, and 316, 1300, 669, 142, 383 genes were expressed differently in ovary, hepatopancreas, brain ganglion, eyestalk, and mandibular organ respectively. Differentially expressed gene enrichment analysis revealed several crucial pathways including protein digestion and absorption, pancreatic secretion, insect hormone biosynthesis, drug metabolism-cytochrome P450 and signal transduction pathway. Through this study, some key genes in correlation with the ovarian development and nutrition metabolism were significantly affected by estradiol, such as vitelline membrane outer layer 1-like protein, heat shock protein 70, Wnt5, JHE-like carboxylesterase 1, cytochrome P302a1, crustacean hyperglycemic hormone, neuropeptide F2, trypsin, carboxypeptidase B, pancreatic triacylglycerol lipase-like, and lipid storage droplet protein. Moreover, RT-qPCR validation demonstrated that expression of transcripts related to ovarian development (vitelline membrane outer layer 1-like protein and cytochrome P302a1) and nutrition metabolism (trypsin, glucose dehydrogenase and lipid storage droplet protein) were significantly affected by estradiol treatment. This study not only has identified relevant genes and several pathways that are involved in estradiol regulation on ovarian development of P. trituberculatus, but also provided new insight into the understanding of the molecular function mechanisms of estradiol in crustacean.

Alternative splice variants and differential relative abundance patterns of vasa mRNAs during gonadal development in the Chinese mitten crab Eriocheir sinensis

Gonadal development usually involves alternative splicing of sex-related genes. Vasa, a highly conserved ATP-dependent RNA helicase present mainly in germ cells, has an important function in gonadal development. As an important sex-related gene, recent evidence indicates that different splice variants of vasa exist in many species. In this study, there was identification of two types of vasa splice variants in the Chinese mitten crab Eriocheir sinensis, termed Esvasa-l and Esvasa-s, respectively. Furthermore, splice variants of Esvasa-s were sub-divided into Esvasa-s1, Esvasa-s2, Esvasa-s3, Esvasa-s4, and Esvasa-s5, based on differing numbers of TGG repeats. Results from genomic structure analyses indicated that these forms are alternatively spliced transcripts from a single vasa gene. Results from tissue distribution assessments indicate the vasa splice variants were exclusively expressed in the gonads of male and female adult crabs. In situ hybridization results indicate Esvasa mRNA was mainly present in the cytoplasm of previtellogenic oocytes. As oocyte size increased, relative abundance of Esvasa mRNA decreased and became distributed near the cellular membrane. The Esvasa mRNA was not detectable in mature oocytes. In testis, Esvasa mRNA was detected in spermatids and spermatozoa, but not in spermatogonia and spermatocytes. Notably, results from qPCR analysis of Esvasa-l and Esvasa-s indicate there are different relative proportions during gametogenesis, implying that splice variants of the Esvasa gene may have different biological functions during crab gonadal development.

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.

PL10 DEAD-Box Protein is Expressed during Germ Cell Differentiation in the Reptile Podarcis sicula (Family Lacertidae)

Among genes involved in the regulation of germ cell differentiation, those of DDX4/Vasa and the Ded1/DDX3 subfamilies encode for DEAD-box ATP-dependent RNA helicases, proteins involved in many mechanisms related to RNA processing. For the first time in reptiles, using specific antibodies at confocal microscopy, we analysed the localization pattern of a Ded1/DDX3 subfamily member in testis and ovary of Podarcis sicula (Ps-PL10) during the reproductive cycle. In testis, Ps-PL10 is expressed in the cytoplasm of spermatocytes and it is not detected in spermatogonia. Differently from Ps-VASA, in round spermatids, Ps-PL10 is not segregated in the chromatoid body but it accumulates in the cytoplasm of residual bodies, and mature spermatozoa are unstained. These observations suggest that in males, Ps-PL10 (1) is involved in spermatogenesis and (2) is then eliminated with residual bodies. In the ovary, Ps-PL10 is present with granules in the cytoplasm of early meiotic cells of the germinal bed (GB), while it is not present in oogonia and somatic cells of the GB stroma. In follicular cells of ovarian follicles, Ps-PL10 expression starts after their fusion with the oocyte. Numerous Ps-PL10 spots are visible in pyriform (nurse-like) cells concomitantly with the protein accumulation in the cytoplasm of differentiating oocyte. In pyriform cells, Ps-PL10 spots are present in the cytoplasm and nuclei, as observed for Ps-VASA, and in the nucleoli, suggesting for Ps-PL10 a role in rRNA processing and in the transport of molecules from the nucleus to cytoplasm and from nurse cells to the oocyte.

Transcriptome Analysis and Differential Gene Expression on the Testis of Orange Mud Crab, Scylla olivacea, during Sexual Maturation

Adequate genetic information is essential for sustainable crustacean fisheries and aquaculture management. The commercially important orange mud crab, Scylla olivacea, is prevalent in Southeast Asia region and is highly sought after. Although it is a suitable aquaculture candidate, full domestication of this species is hampered by the lack of knowledge about the sexual maturation process and the molecular mechanisms behind it, especially in males. To date, data on its whole genome is yet to be reported for S. olivacea. The available transcriptome data published previously on this species focus primarily on females and the role of central nervous system in reproductive development. De novo transcriptome sequencing for the testes of S. olivacea from immature, maturing and mature stages were performed. A total of approximately 144 million high-quality reads were generated and de novo assembled into 160,569 transcripts with a total length of 142.2 Mb. Approximately 15–23% of the total assembled transcripts were annotated when compared to public protein sequence databases (i.e. UniProt database, Interpro database, Pfam database and Drosophila melanogaster protein database), and GO-categorised with GO Ontology terms. A total of 156,181 high-quality Single-Nucleotide Polymorphisms (SNPs) were mined from the transcriptome data of present study. Transcriptome comparison among the testes of different maturation stages revealed one gene (beta crystallin like gene) with the most significant differential expression—up-regulated in immature stage and down-regulated in maturing and mature stages. This was further validated by qRT-PCR. In conclusion, a comprehensive transcriptome of the testis of orange mud crabs from different maturation stages were obtained. This report provides an invaluable resource for enhancing our understanding of this species’ genome structure and biology, as expressed and controlled by their gonads.

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

Sex-biased genes are considered to account for most of phenotypic differences between males and females. In order to explore the sex-biased gene expression in crab, we performed the whole-body transcriptome analysis in male and female juveniles of the Chinese mitten crab Eriocheir sinensis using next-generation sequencing technology. Of the 23,349 annotated unigenes, 148 were identified as sex-related genes. A total of 29 candidate genes involved in primary sex determination pathways were detected, indicating the sex determination cascade of the mitten crab might be more complex than previously supposed. Differential expression analysis showed 448 differentially expressed genes (DEGs) between the two transcriptomes. Most of DEGs were involved in processes such as metabolism and immunity, and not associated with obvious sexual function. The pathway predominantly enriched for DEGs were related to lysosome, which might reflect the differences in metabolism between males and females. Of the immune DGEs, 18 up-regulated genes in females were humoral immune factors, and eight up-regulated genes in males were pattern recognition receptors, suggesting sex differences of immune defense might exist in the mitten crab. In addition, two reproduction-related genes, vitellogenin and insulin-like androgenic gland factor, were identified to express in both sexes but with significantly higher level in males. Our research provides the first whole-body RNA sequencing of sex-specific transcriptomes for juvenile E. sinensis and will facilitate further studies on molecular mechanisms of crab sexual dimorphism.

High-density linkage mapping aided by transcriptomics documents ZW sex determination system in the Chinese mitten crab Eriocheir sinensis

The sex determination system in crabs is believed to be XY-XX from karyotypy, but centromeres could not be identified in some chromosomes and their morphology is not completely clear. Using quantitative trait locus mapping of the gender phenotype, we revealed a ZW-ZZ sex determination system in Eriocheir sinensis and presented a high-density linkage map covering ~98.5% of the genome, with 73 linkage groups corresponding to the haploid chromosome number. All sex-linked markers in the family we used were located on a single linkage group, LG60, and sex linkage was confirmed by genome-wide association studies (GWAS). Forty-six markers detected by GWAS were heterozygous and segregated only in the female parent. The female LG60 was thus the putative W chromosome, with the homologous male LG60 as the Z chromosome. The putative Z and W sex chromosomes were identical in size and carried many homologous loci. Sex ratio (5:1) skewing towards females in induced triploids using unrelated animals also supported a ZW-ZZ system. Transcriptome data were used to search for candidate sex-determining loci, but only one LG60 gene was identified as an ankyrin-2 gene. Double sex- and mab3-related transcription factor 1 (Dmrt1), a Z-linked gene in birds, was located on a putative autosome. With complete genome sequencing and transcriptomic data, more genes on putative sex chromosomes will be characterised, thus leading towards a comprehensive understanding of the sex determination and differentiation mechanisms of E. sinensis, and decapod crustaceans in general.

Characteristic of PGDS potential regulation role on spermatogenesis in the Chinese mitten crab Eriocheir sinensis

Prostaglandin D synthase (PGDS) catalyzes the isomerization of PGH2 to produce PGD2 in the presence of sulfhydryl compounds. In this study, a full length PGDS gene comprising 1250 nucleotides from the Chinese mitten crab Eriocheir sinensis (Es-PGDS) was characterized, with a 615 bp open reading frame encoding 204 amino acid residues. Its deduced peptide has high homology with other species' PGDS protein. The Es-PGDS mRNA expression was tissue-related, with the highest expression observed in the hepatopancreas, accessory sex gland, testis and ovaries. We also detected the different stages of tissue expression and the enzyme activity for Es-PGDS in the testis and male crab hepatopancreas. The different expression patterns and its corresponding enzyme activity level indicated that PGDS is involving in the regulation of reproductive action during the period of rapid development in E. sinensis. Furthermore our research could arouse a heat debate on the PGDS reproductive function in invertebrate and further study will be needed to determine the molecular mechanism(s) linking PGDS functions to spermatogenesis and ontogenesis if this gene is to be exploited as a molecular biomarker in further studies of development.

Localization of germline marker vasa homolog RNA to a single blastomere at early cleavage stages in the oriental river prawn Macrobrachium nipponense: evidence for germ cell specification by preformation

Germ cells are specified by the inheritance of maternal germline determinants (preformation mode) or inductive signals from somatic cells (epigenesis mode) during embryogenesis. However, the germline specification in decapod crustaceans is unclear so far. Using vasa homolog (MnVasa) as a germ cell marker, here we probed the early events of germline specification in the oriental river prawn Macrobrachium nipponense. Quantitative RT-PCR analysis of unfertilized eggs and embryos demonstrated that the prawn MnVasa mRNA is a maternal factor. Whole-mount in situ hybridization further indicated that MnVasa transcripts are maternally supplied to only one blastomere at the very early cleavage stages. As cleavage proceeds, the MnVasa-positive blastomere undergoes proliferation and increases in number. During gastrulation, the MnVasa-positive cells are found to be around a blastopore and could migrate into an embryo through the blastopore. At the zoea stage, clusters of the MnVasa-positive cells distribute not only in the gonad rudiment in the cephalothorax but also at an extragonadic site, dorsal to the posterior hindgut in the abdomen, suggesting that MnVasa-positive cells could migrate anteriorly to the genital rudiment through the hindgut. Based on the dynamic localization and number of MnVasa-positive cells during embryogenesis, we concluded that the MnVasa-positive cells are primordial germ cells (PGC) or founder cells of PGC that are separated from soma at the early cleavage stage. MnVasa mRNA might have a key function in the specification of the prawn germline cells as a maternal determinant. These results provide the first evidence that the germline specification in decapod crustaceans follows a preformation mode.