Xiwei Jia

Transcriptome analysis of the differences in gene expression between testis and ovary in green mud crab (Scylla paramamosain)

BackgroundThe green mud crab (Scylla paramamosain) is the most prevalent crustacean on the southeast coast of China. The molecular regulatory mechanism of sex determination and gonadal differentiation in this species has received considerable attention in recent years because of the huge differences—both biological and economic—between male and female crabs. In this study, next-generation sequencing technology was used to develop deep-coverage transcriptomic sequencing data for the testis and ovary of S. paramamosain.ResultsA total of 365,116 reads (testis 171,962, ovary 193,154) with an average sequence length of 285 bp were produced from testis and ovary cDNA libraries. After filtering out contaminating reads, the clean reads were assembled, producing a total of 21,791 isotigs and leaving 22,814 reads as singlets. Using the BLASTX program, 3,471 unique sequences (2,275 isotigs and 1,196 singletons) were annotated with known protein sequences from the NCBI non-redundant (Nr) protein sequence database. The Gene Ontology and KEGG (Kyoto Encyclopedia of Genes and Genomes) analyses allowed the 224 unique sequences that were annotated with enzyme code (EC) numbers to be mapped into 174 KEGG pathways. After comparing the ovary and testis libraries, 4,021 gonad-differentially, 10,522 ovary-specifically, and 19,013 testis-specifically expressed genes were identified. Moreover, 33 ovary-specific, 14 testis-specific, and 34 gonad-differential transcripts were confirmed by semi-quantitative PCR and quantitative real-time PCR. In addition, 8,610 putative simple sequence repeats (SSRs) and 23,879 potential single nucleotide polymorphisms (SNPs) were identified.ConclusionThis is the first large-scale RNA sequencing of S. paramamosain to be reported. We have identified many important functional genes and made a preliminary attempt to construct the regulatory network involved in the gonadal development of crustaceans. The annotated transcriptome data will provide fundamental support for future research into the reproduction biology of S. paramamosain. A large number of candidate SSRs and SNPs were detected, which could be used as genetic markers for population genetics and functional genomics in this species.

Characterization and expression profile of Vitellogenin gene from Scylla paramamosain

The full-length (7816 bp) cDNA of Vitellogenin (Vtg) encoding 2560 aa with an estimated molecular mass of 287.743 kDa was cloned from the green mud crab Scylla paramamosain. Semi-quantitative PCR (sq-PCR) revealed a specific expression pattern of Sp-vtg gene in ovaries and hepatopancreas. With the development of ovaries, the expression level of Sp-vtg gene showed an increasing trend both in ovaries and hepatopancreas, and the expression level of Sp-vtg gene in hepatopancreas and ovary was stable after stage IV. By in situ hybridization, the positive signals of Sp-vtg gene were detected in the cytoplasm of oocytes in stage I, in the follicle cell and the surrounding of the nucleus in stage III, and in the nucleus in stage V. Furthermore, the signals become stronger with the later development stages of ovary. Moreover, in situ hybridization analysis revealed that positive signals of Sp-vtg gene are present in the hepatopancreatic tubule, and the signals increase during the development, becoming the strongest in stage V. Our results indicate that both ovaries and hepatopancreas are sites of Vitellogenin gene synthesis in S. paramamosain.

Identification and comparative analysis of the ovary and testis microRNAome of mud crab Scylla paramamosain : JIA et al.

The role of microRNA (miRNA) in reproductive regulation is attracting increasingly more attention. In this study, we obtained 9,643,114 and 15,498,999 raw reads from the ovary and testis library of important farmed mud crab Scylla paramamosain, respectively. After data mining, a total of 4,096,464 and 11,737,973 mappable small RNA sequences remained for analysis. By mapping to the reference genome and expressed sequence tag (EST) of Daphnia pulex and other crabs, a total of 1,417 miRNAs were identified. On the basis of 1,417 miRNAs, 514 (36.3%) unique miRNAs coexpressed in the gonad of female and male libraries, and 336 (23.7%) and 567 (40%) expressed preferentially in female and male libraries, respectively. Analysis of library sequencing data resulted in the identification of 108 miRNAs (out of 1,417; 7.6%) that showed significant differential expression between the two samples. Of these, 13 miRNAs were expressed only in the testis, two miRNAs were expressed only in the ovary, and 93 miRNAs were coexpressed: 57 (61.3%) were upregulated (ovary/testis) and 36 (38.7%) were downregulated (ovary/testis). To confirm the expression patterns of the predicted miRNAs, we randomly selected 14 candidate miRNAs from 108 differentially expressed miRNAs and performed stem–loop real time quantitative PCR (RT‐qPCR) assays in five ovary developing stages. Five miRNAs showed similar expression patterns in almost every stage as those revealed by identification of differentially expressed genes (IDEG6) analysis. The above five miRNAs were predicted to match the 3′‐untranslated region of the published S. paramamosain gene. Four out of five miRNA had a regulation effect on many genes, especially the genes related to gonadal development.

miR-34 regulates reproduction by inhibiting the expression of MIH, CHH, EcR and FAMeT genes in mud crab Scylla paramamosain : ZHOU et al.

Mud crab Scylla paramamosain is a commercially important species widely cultured in China. It is well known that the eyestalk regulates reproductive activities in crustaceans. In our previous research, we found that the miR‐34 expression level in male eyestalk was significantly higher than that in females. Thus, we assumed that it may play an important role in regulating reproduction. In this study, we used bioinformatic tools to identify the target genes of miR‐34 in eyestalk. Six reproduction‐related genes with an intact 3′‐untranslated region (UTR), including molt‐inhibiting hormone (MIH), crustacean hyperglycemic hormone (CHH), vitellogenesis‐inhibiting hormone, red pigment concentrating hormone, ecdysone receptor (EcR), and farnesoic acid methyltransferase (FAMeT) were identified. When the 3′‐UTR plasmid vectors of the six genes were cotransfected with miR‐34 mimics into 293FT cells, respectively, the luciferase activities of four genes (MIH, CHH, EcR, and FAMeT) were significantly decreased compared with that in the control group; on the contrary, when the six plasmid vectors were cotransfected with the miR‐34 inhibitor respectively, the luciferase activities of four genes (MIH, CHH, EcR, and FAMeT) were significantly higher than that in the control group. When agomiR‐34 and antagomiR‐34 were injected into the eyestalk respectively in vivo, the expression levels of the MIH, CHH, EcR, and FAMeT genes were detected by a quantitative real‐time polymerase chain reaction. The results showed that agomiR‐34 suppressed the expression of the four genes, whereas antagomiR‐34 enhanced their expression. These experimental results confirmed our hypothesis that miR‐34 may indirectly regulate reproduction via binding to the 3′‐UTRs of MIH, CHH, EcR, and FAMeT genes and suppressing their expression.