Differential expression of ubiquitin-conjugating enzyme E2r in the developing ovary and testis of penaeid shrimp Marsupenaeus japonicus

Gonadal transcriptome analysis of genes related to sex differentiation and sex development in the Pomacea canaliculata

As an invasive alien animal, Pomacea canaliculata poses a great danger to the ecology and human beings. Recently, there has been a gradual shift towards bio-friendly control. Based on the development of RNA interference and CRISPR technology as molecular regulatory techniques for pest control, it was determined if the knockout of genes related to sex differentiation in P. canaliculata could induce sterility, thereby helping in population control. However, the knowledge of sex differentiation- and development-related genes in P. canaliculata is currently lacking. Here, transcriptomic approaches were used to study the genes expressed in the two genders of P. canaliculata at various developmental stages. Gonad transcriptomes of immature or mature males and females were compared, revealing 12,063 genes with sex-specific expression, of which 6066 were male- and 5997 were female-specific. Among the latter, 581 and 235 genes were up-regulated in immature and mature females, respectively. The sex-specific expressed genes identified included GnRHR2 and TSSK3 in males and ZAR1 and WNT4 in females. Of the genes, six were involved in reproduction: CCNBLIP1, MND1, DMC1, DLC1, MRE11, and E(sev)2B. Compared to immature snail gonads, the expression of HSP90 and CDK1 was markedly reduced in gonadal. It was hypothesized that the two were associated with the development of females. These findings provided new insights into crucial genetic information on sex differentiation and development in P. canaliculata. Additionally, some candidate genes were explored, which can contribute to future studies on controlling P. canaliculata using molecular regulatory techniques.

Genome-wide identification and transcriptome-based expression profiling of E2 gene family: Implication for potential roles in gonad development of Crassostrea gigas

In this study, we investigated the role of E2 ubiquitin conjugating enzymes (E2) in the Pacific oyster Crassostrea gigas, with a focus on their involvement in gonad development. We identified 34 E2 genes clustered into nine subgroups and 24 subfamilies. The gene structure and intron-exon location were conserved within the same subfamily, but motif variation suggested functional diversity. Tissue transcriptome analyses revealed that most E2 genes were broadly expressed, with UBE2CL showing specific expression in the female gonad. Expression profiling of E2 genes during early embryo-larvae development stages suggested that five E2 genes were highly expressed in early embryo development, indicating their involvement in cell division processes. Furthermore, by profiling the expression of E2 genes in different gonadal developmental stages, we observed a gradual increase in expression for four genes during gametogenesis, with significantly higher expression in the female gonad at the maturation stage. Similarly, five E2 genes displayed elevated expression levels in the male gonad at the maturation stage, indicating their crucial roles in gonadal development and gametogenesis. Our study provides valuable insights into the potential functions of the E2 gene family in C. gigas, shedding light on the molecular mechanisms underlying gonad development in oysters.

miRNA-mRNA analysis of sheep adrenal glands reveals the network regulating reproduction

Background

The adrenal gland participates in the process of sheep reproduction. MicroRNAs (miRNAs), endogenous small noncoding RNAs, regulate gene expression at the posttranscriptional level. However, the miRNA-mRNA network profile of adrenal glands relating to reproduction in sheep is still not well-studied. As sheep with FecB^BB genotype show higher lambing number compare with the sheep with FecB⁺⁺ genotype. This research aims to compare gene expression by small RNA-seq in adrenal tissues at follicular (F) and luteal (L) phases in FecB^BB (MM) and FecB⁺⁺ (ww) sheep. After analysis of gene expression, significant differentially expressed microRNAs (DEMs) and corresponding target genes were identified.

Results

A total of 180 miRNAs were found in this study, of which 19 DEMs were expressed in the four comparison groups (MM_F_A vs. MM_L_A, MM_F_A vs. ww_F_A, MM_L_A vs. ww_L_A, ww_F_A vs. ww_L_A). Subsequently, 354 target genes of 19 DEMs were predicted by integrated analysis. Cluster analysis was performed by K_means_cluster, and the expression patterns of these DEMs were separated into four subclusters. Functional analysis of target genes was performed by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG). The results indicated that the target genes were involved mainly in the Notch signaling pathway, signal transduction, cell communication, innate immune response and amino acid metabolism. Specifically, the Notch signaling pathway, biosynthetic process and metabolic process of pyrimidine nucleotide and amino acid metabolism appear to play key regulatory roles in the sheep fertility trait. Furthermore, miRNA-mRNA interaction networks were constructed by differentially expressed genes combined with our previous study of transcriptome data. The results showed that several key genes, including TDRD3, ANAPC7, CCNL2, BRD2 and MUT, were related to the transformation from the follicular phase to the luteal phase. PLAC8L1, NFAT5, DDX24 and MBD1 were related to the high fecundity of small tail Han sheep.

Conclusions

In this study, the miRNA transcriptome profile was identified, and miRNA-mRNA interaction networks were constructed in adrenal gland tissue of small tail Han sheep, the interaction between miR-370-3p and its targets were considered to play a major role in the reproduction regulation process. The results enriched the number of known miRNAs in adrenal glands and provided novel ideas and further information to demonstrate posttranscriptional regulation mechanisms at follicular and luteal phases in different genotypes of small tail Han sheep.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12863-022-01060-y.

Knockdown of E1- and E2-ubiquitin enzymes triggers defective chorion biogenesis and modulation of autophagy-related genes in the follicle cells of the vector Rhodnius prolixus

In insects, the last stage of oogenesis is the process where the chorion layers (eggshell) are synthesized and deposited on the surface of the oocytes by the follicle cells. Protein homeostasis is determined by the fine-tuning of translation and degradation pathways, and the ubiquitin-proteasome system is one of the major degradative routes in eukaryotic cells. The conjugation of ubiquitin to targeted substrates is mediated by the ordered action of E1-activating, E2-conjugating, and E3-ligase enzymes, which covalently link ubiquitin to degradation-targeted proteins delivering them to the proteolytic complex proteasome. Here, we found that the mRNAs encoding polyubiquitin (pUbq), E1, and E2 enzymes are highly expressed in the ovaries of the insect vector of Chagas Disease Rhodnius prolixus. RNAi silencing of pUbq was lethal whereas the silencing of E1 and E2 enzymes resulted in drastic decreases in oviposition and embryo viability. Eggs produced by the E1- and E2-silenced insects presented particular phenotypes of altered chorion ultrastructure observed by high-resolution scanning electron microscopy as well as readings for dityrosine cross-linking and X-ray elemental microanalysis, suggesting a disruption in the secretory routes responsible for the chorion biogenesis. In addition, the ovaries from silenced insects presented altered levels of autophagy-related genes as well as a tendency of upregulation in ER chaperones, indicating a disturbance in the general biosynthetic-secretory pathway. Altogether, we found that E1 and E2 enzymes are essential for chorion biogenesis and that their silencing triggers the modulation of autophagy genes suggesting a coordinated function of both pathways for the progression of choriogenesis.

Comparative Transcriptomics of Gonads Reveals the Molecular Mechanisms Underlying Gonadal Development in Giant Freshwater Prawns (Macrobrachium rosenbergii)

The giant freshwater prawn, Macrobrachium rosenbergii, is a prawn that has economic significance throughout the world. It exhibits sex-related growth dimorphism, whereby the males grow significantly more rapidly than the females. Therefore, a study on the molecular regulatory mechanism, which underlies the sexual differentiation of M. rosenbergii, is of both scientific and commercial importance. However, a scarcity of genomic and transcriptomic resources severely limits our knowledge of the sexual differentiation mechanisms in M. rosenbergii. Here, transcriptome sequencing of several gonadic samples of males and females in M. rosenbergii was performed to investigate the molecular basis underlying gonadal development. Our results showed that 2149 unigenes presented as differentially expressed genes (DEGs) in the ovaries of females compared to the testes of males, which contained 484 down-regulated and 1665 up-regulated genes. Enrichment analysis of DEGs revealed many of these genes to be related to sexual differentiation and gonadal development. From our transcriptome analyses, and as confirmed by quantitative real-time PCR, male-related genes (Mrr, MRPINK, IR, IAGBP, TESK1, and dsx) in the testes were significantly up-regulated, and female-related genes (ERR, Sxl3, cyclinB, Dmrt99B, PPP2A, and ADCY9) in the ovaries were also significantly up-regulated. This indicates the potential role these genes play in the gonadal development of M. rosenbergii. Furthermore, multiple signal transduction pathways relating to gonadal maturation and spermatogenesis, including MAPK, were identified herein. Our data also supports previous ideas that IAG and IAGBP-IR signaling schemes could help in the regulation of testis’ development in M. rosenbergii and the ERR gene could regulate ovarian development by affecting the expression of cyclinB, PPP2A, and ADCY9. The data from this study provides incredibly usefully genomic resources for future research on the sexual differentiation and practical aquaculture of M. rosenbergii.

Transcriptome analysis of genes related to gonad differentiation and development in Muscovy ducks

BACKGROUND

Sex-related genes play a crucial role in gonadal differentiation into testes or ovaries. However, the genetic control of gonadal differentiation in Muscovy ducks remains unknown. Therefore, the objective of our study was to screen new candidate genes associated with ovarian and testicular development.

RESULTS

In this study, 24 males before gonadal differentiation (MB), 24 females before gonadal differentiation (FB), 24 males after gonadal differentiation (MA) and 24 females after gonadal differentiation (FA) were selected from Putian Muscovy ducks, forming 4 groups. RNA-Seq revealed 101.76 Gb of clean reads and 2800 differentially expressed genes (DEGs), including 46 in MB vs FB, 609 in MA vs FA, 1027 in FA vs FB, and 1118 in MA vs MB. A total of 146 signalling pathways were enriched by KEGG analysis, among which 20, 108, 108 and 116 signalling pathways were obtained in MB vs FB, MA vs MB, MA vs FA and FA vs FB, respectively. In further GO and KEGG analyses, a total of 21 candidate genes related to gonad differentiation and development in Muscovy ducks were screened. Among these, 9 genes were involved in the differentiation and development of the testes, and 12 genes were involved in the differentiation and development of the ovaries. In addition, RNA-Seq data revealed 2744 novel genes.

CONCLUSIONS

RNA-Seq data revealed 21 genes related to gonadal differentiation and development in Muscovy ducks. We further identified 12 genes, namely, WNT5B, HTRA3, RSPO3, BMP3, HNRNPK, NIPBL, CREB3L4, DKK3, UBE2R2, UBPL3KCMF1, ANXA2, and OSR1, involved in the differentiation and development of ovaries. Moreover, 9 genes, namely, TTN, ATP5A1, DMRT1, DMRT3, AMH, MAP3K1, PIK3R1, AGT and ADAMTSL1, were related to the differentiation and development of testes. Moreover, after gonadal differentiation, DMRT3, AMH, PIK3R1, ADAMTSL1, AGT and TTN were specifically highly expressed in males. WNT5B, ANXA2 and OSR1 were specifically highly expressed in females. These results provide valuable information for studies on the sex control of Muscovy ducks and reveal novel candidate genes for the differentiation and development of testes and ovaries.

Identification of Candidate Genes for the Plateau Adaptation of a Tibetan Amphipod, Gammarus lacustris, Through Integration of Genome and Transcriptome Sequencing

The amphipod Gammarus lacustris has been distributing in the Tibetan region with well-known uplifts of the Tibetan plateau. It is hence considered as a good model for investigating stress adaptations of the plateau. Here, we sequenced the whole-genome and full-length transcriptome of G. lacustris, and compared the transcriptome results with its counterpart Gammarus pisinnus from a nearby plain. Our main goal was to provide a genomic resource for investigation of genetic mechanisms, by which G. lacustris adapted to living on the plateau. The final draft genome assembly of G. lacustris was 5.07 gigabases (Gb), and it contained 443,304 scaffolds (>2 kb) with an N50 of 2,578 bp. A total of 8,858 unigenes were predicted in the full-length transcriptome of G. lacustris, with an average gene length of 1,811 bp. Compared with the G. pisinnus transcriptome, 2,672 differentially expressed genes (DEGs) were up-regulated and 2,881 DEGs were down-regulated in the G. lacustris transcriptome. Along with these critical DEGs, several enriched metabolic pathways, such as oxidative phosphorylation, ribosome, cell energy homeostasis, glycolysis and gluconeogenesis, were predicted to play essential roles in the plateau adaptation. In summary, the present study provides a genomic basis for understanding the plateau adaption of G. lacustris, which lays a fundamental basis for further biological and ecological studies on other resident aquatic species in the Tibetan plateau.

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.

Integrated analysis of microRNA and mRNA expression profiles during the sex-differentiation sensitive period in oriental river prawn, Macrobrachium nipponense

Male oriental river prawns (Macrobrachium nipponense) grow faster than females, and therefore, reach larger sizes by harvest time. Histological observations have indicated that the sex-differentiation sensitive period (which includes the formation of the androgenic gland, the testis, and the ovary) is from post-larvae (PL) developmental stage for M. nipponense. In this study, we prepared four microRNA (miRNA) and mRNA libraries using samples collected from sex-differentiation sensitive period (PL7 to PL16) to perform RNA-sequencing for identifying sex-related candidate miRNAs, genes, and metabolic pathways. A total of nine intersection miRNAs were identified, of which three were highly expressed in the androgenic gland, and their expression was verified by quantitative Real-Time PCR (qPCR). These three miRNAs and their 11 predicted target genes may be strong candidates for sex-related miRNAs and sex-related genes in M. nipponense. Five vital sex-related metabolic pathways were also identified that may regulate other sex-differentiation and sex-determination mechanisms. Finding of the study provide important insights to enhance our understanding on sex-differentiation and sex-determination mechanisms for M. nipponense.

Ovarian Transcriptome Analysis of Vitellogenic and Non-Vitellogenic Female Banana Shrimp (Fenneropenaeus merguiensis)

The banana shrimp (Fenneropenaeus merguiensis) is one of the most commercially important penaeid species in the world. Its numbers are declining in the wild, leading to a loss of broodstock for farmers of the shrimp and a need for more successful breeding programs. However, the molecular mechanism of the genes involved in this shrimp’s ovarian maturation is still unclear. Consequently, we compared transcriptomic profiles of ovarian tissue from females in both the vitellogenic stage and the non-vitellogenic stage. Using RNA-Seq technology to prepare the transcriptome libraries, a total of 12,187,412 and 11,694,326 sequencing reads were acquired from the non-vitellogenic and vitellogenic stages respectively. The analysis of the differentially expressed genes identified 1,025 which were significantly differentially expressed between the two stages, of which 694 were up-regulated and 331 down-regulated. Four genes putatively involved in the ovarian maturation pathway were chosen for validation by quantitative real-time PCR (RT-qPCR). The data from this study provided information about gene expression in ovarian tissue of the banana shrimp which could be useful for a better understanding of the regulation of this species’ reproductive cycle.

Characterisation and expression analysis of UBC9 and UBS27 genes in developing gonads of cicindelids (Coleoptera: Cicindelidae)

Ubiquitin and small ubiquitin-like modifiers (SUMO) are post-translational modifiers essential in a variety of cellular processes, including gametogenesis. SUMO-conjugating enzyme (UBC9) and the ubiquitin ribosomal fusion protein UBS27 have been characterised in several model species. However, their expression in coleopteran remains unstudied. In this study, UBC9 and UBS27 genes have been characterised in the tiger beetle Cicindela campestris for the first time. Bioinformatic analysis showed that the Cc-UBC9 gene encoded a 159 amino acid protein with a predicted molecular weight of 18.18kDa, and the Cc-UBS27 gene encoded a 156 amino acid protein with a predicted molecular weight of 17.71kDa. Selection analyses carried out in several cicindelid species revealed that both genes were affected by purifying selection. Real time quantitative PCR analysis demonstrated that Cc-UBC9 and Cc-UBS27 were expressed in different tissues. The highest expression on both genes was found in the ovary and testis, and there were differential expression levels between immature and mature stages of testis development. The expression patterns of Cc-UBC9 and Cc-UBS27 suggest that these genes play important roles in gametogenesis in C. campestris. This information is relevant to better understand the reproductive process in cicindelids and the function of ubiquitin and small ubiquitin-related modifier genes in the Coleoptera.

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.

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

Background

The 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.

Results

A 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.

Conclusion

This 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.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-585) contains supplementary material, which is available to authorized users.

Identification and expression analysis of differentially expressed genes from shrimp (Penaeus monodon) in response to low salinity stress

Four suppression subtractive hybridization (SSH) cDNA libraries were constructed to identify differentially expressed salinity stress responsive genes of black tiger shrimp, Penaeus monodon exposed to low (3 ppt) salinity conditions. Forward and reverse SSH cDNA libraries were developed from the gill and gut tissues of shrimp and clones having inserts larger than 300 bp were unidirectionally sequenced. Based on the sequence homology search, the identified genes were categorized for their putative functions related to a wide range of biological roles, such as nucleic acid regulation and replication, immune response, energy and metabolism, cell signaling, cellular process, cytoskeleton and membrane structure, stress and osmoregulation. Gene expression levels in response to low salinity conditions at 2 weeks post salinity stress of thirteen selected differentially expressed genes identified from SSH cDNA libraries (14-3-3 like protein, crust in, lysozyme, arginine kinase, Naþ/Kþ-ATPase a-subunit, intracellular fatty acid binding protein, cathepsin B, anti-lipopolysaccharide factor, ferritin, ubiquitin conjugating enzyme E2, calreticulin, innexin 2 and heat shock protein 21) were analyzed by RT-PCR. The highest gene expression levels were observed for Naþ/Kþ-ATPase a-subunit (34.28-folds) in gill tissues, intracellular fatty acid binding protein (13.30-folds) in gut tissues and innexin 2 (14.43-folds) in muscle tissues respectively. The differential and significant levels of gene expression indicate the functional role of these genes in shrimp salinity stress adaptive mechanisms.

Transcriptome analysis of androgenic gland for discovery of novel genes from the oriental river prawn, Macrobrachium nipponense, using Illumina Hiseq 2000

Background

The oriental river prawn, Macrobrachium nipponense, is an important aquaculture species in China, even in whole of Asia. The androgenic gland produces hormones that play crucial roles in sexual differentiation to maleness. This study is the first de novo M. nipponense transcriptome analysis using cDNA prepared from mRNA isolated from the androgenic gland. Illumina/Solexa was used for sequencing.

Methodology and Principal Finding

The total volume of RNA sample was more than 5 ug. We generated 70,853,361 high quality reads after eliminating adapter sequences and filtering out low-quality reads. A total of 78,408 isosequences were obtained by clustering and assembly of the clean reads, producing 57,619 non-redundant transcripts with an average length of 1244.19 bp. In total 70,702 isosequences were matched to the Nr database, additional analyses were performed by GO (33,203), KEGG (17,868), and COG analyses (13,817), identifying the potential genes and their functions. A total of 47 sex-determination related gene families were identified from the M. nipponense androgenic gland transcriptome based on the functional annotation of non-redundant transcripts and comparisons with the published literature. Furthermore, a total of 40 candidate novel genes were found, that may contribute to sex-determination based on their extremely high expression levels in the androgenic compared to other sex glands,. Further, 437 SSRs and 65,535 high-confidence SNPs were identified in this EST dataset from which 14 EST-SSR markers have been isolated.

Conclusion

Our study provides new sequence information for M. nipponense, which will be the basis for further genetic studies on decapods crustaceans. More importantly, this study dramatically improves understanding of sex-determination mechanisms, and advances sex-determination research in all crustacean species. The huge number of potential SSR and SNP markers isolated from the transcriptome may shed the lights on research in many fields, including the evolution and molecular ecology of Macrobrachium species.

Molecular cloning, characterization, and gene expression of the androgen receptor in the large yellow croaker, Larimichthys crocea

Androgens mediate a wide range of physiological responses and developmental processes in vertebrates, involving both reproductive and nonreproductive systems. The activity of androgens is mediated by the androgen receptor (AR), a member of the nuclear receptor superfamily. In this study, an AR gene was cloned from the large yellow croaker (Larimichthys crocea) for the first time. qRT-PCR revealed ubiquitous expression of AR in all adult tissues examined, with higher expression in the gonad and liver of both sexes and highest expression in the blastula stage of embryonic development. Using in situ hybridization, we detected positive signals of AR in the spermatogonium, spermatocyte, spermatid, and spermatozoon during spermatogenesis, in the cytoplasm of all oocytes during oogenesis and in the follicle cells of stage IV oocytes. Our findings support the important role that AR plays in gametogenesis, gonadal development, and the early stages of embryonic development.

Identification of reproduction-related proteins and characterization of proteasome alpha 3 and proteasome beta 6 cDNAs in testes of the giant tiger shrimp Penaeus monodon

Cellular proteomic analysis was carried out to identify reproduction-related proteins in testes of wild and domesticated broodstock of Penaeus monodon. In total, 642 protein spots were characterized and 287 spots (44.70%) significantly matched protein sequences in the databases (P<0.05). To examine a role of the proteasome system in testicular development of P. monodon, the expression profiles of proteasome alpha 3 subunit (PmPsma3) and proteasome beta 6 (PmPsmb6) mRNA in different groups of domesticated shrimp and in wild broodstock were examined. The expression levels of these transcripts in testes of 18-month-old domesticated shrimp were significantly lower than those of wild broodstock (P<0.05). Interestingly, the expression levels of testicular PmPsma3 and PmPsmb6 in 18-month-old shrimp were significantly increased at 24 h following serotonin injection (50 μg/g body weight). Results suggested that reduced degrees of maturation in captive P. monodon males may be partially resolved by exogenous 5-HT administration.

Transcriptome profiling of testis during sexual maturation stages in Eriocheir sinensis using Illumina sequencing

The testis is a highly specialized tissue that plays dual roles in ensuring fertility by producing spermatozoa and hormones. Spermatogenesis is a complex process, resulting in the production of mature sperm from primordial germ cells. Significant structural and biochemical changes take place in the seminiferous epithelium of the adult testis during spermatogenesis. The gene expression pattern of testis in Chinese mitten crab (Eriocheir sinensis) has not been extensively studied, and limited genetic research has been performed on this species. The advent of high-throughput sequencing technologies enables the generation of genomic resources within a short period of time and at minimal cost. In the present study, we performed de novo transcriptome sequencing to produce a comprehensive transcript dataset for testis of E. sinensis. In two runs, we produced 25,698,778 sequencing reads corresponding with 2.31 Gb total nucleotides. These reads were assembled into 342,753 contigs or 141,861 scaffold sequences, which identified 96,311 unigenes. Based on similarity searches with known proteins, 39,995 unigenes were annotated based on having a Blast hit in the non-redundant database or ESTscan results with a cut-off E-value above 10⁻⁵. This is the first report of a mitten crab transcriptome using high-throughput sequencing technology, and all these testes transcripts can help us understand the molecular mechanisms involved in spermatogenesis and testis maturation.

Identification of robust hypoxia biomarker candidates from fin of medaka (Oryzias latipes)

Aquatic hypoxia caused by organic pollution and eutrophication is a pressing worldwide water pollution problem. Better methods for monitoring oxygen levels are needed to assist efforts to maintain and protect the health of natural aquatic environments. In this project, we used a Japanese ricefish (medaka, Oryzias latipes) 8K oligonucleotide array as a platform to identify potential hypoxic biomarkers in different organs (fin, gill, liver and brain) upon exposure to hypoxia. The microarray results were validated by qRT-PCR employing a subset of candidate biomarkers. Interestingly, the largest number and most significant of hypoxia responding array features were detected in hypoxia exposed fin tissues. We identified 173 array features that exhibited a significant response (over 2 fold change in expression) upon exposure to hypoxic conditions and validated a subset of these by quantitative RT-PCR. These gene targets were subjected to annotation and gene ontology mining. Positively identifiable gene targets that may be useful for development of a rapid and accurate biomarker test using fin clips are discussed in relation to previous reports on hypoxia responsive genes.

Enzyme E2 from Chinese white shrimp inhibits replication of white spot syndrome virus and ubiquitinates its RING domain proteins

Recent studies have shown that the ubiquitin (Ub) proteasome pathway (UPP) is closely related to immune defense. We have identified a ubiquitin-conjugating enzyme, E2, from the Chinese white shrimp, Fenneropenaeus chinensis (FcUbc). Injection of recombinant FcUbc protein (rFcUbc) reduced the mortality of shrimp infected with white spot syndrome virus (WSSV) and inhibited replication of WSSV. rFcUbc, but not a mutant FcUbc (mFcUbc), bound to WSSV RING domains (WRDs) from four potential E3 ligase proteins of WSSV in vitro. Importantly, rFcUbc could ubiquitinate the RING domains (named WRD2 and WRD3) of WSSV277 and WSSV304 proteins in vitro and the two proteins in WSSV-infected Drosophila melanogaster Schneider 2 (S2) cells. Furthermore, overexpression of FcUbc increased ubiquitination of WSSV277 and WSSV304 during WSSV infection. In summary, our study demonstrates that FcUbc from Chinese white shrimp inhibited WSSV replication and could ubiquitinate WSSV RING domain-containing proteins. This is the first report about antiviral function of Ubc E2 in shrimp.

EST analysis on the gonad development related organs and microarray screen for differentially expressed genes in mature ovary and testis of Scylla paramamosain

A total of 5160 high quality ESTs (expressed sequence tags) averaging 357 bp were collected from normalized cDNA libraries created from testis, ovary and mixed organs of mud crab Scylla paramamosain. Clustering and assembly of these ESTs resulted in a total of 3837 unique sequences with 576 overlapping contigs and 3261 singletons. Comparisons with the GenBank non-redundant (Nr) protein database (BLASTx, e-values <10(-5)) revealed putative functions or matched homologs from other organisms for 847 (22%) of the ESTs. Several gonad development related genes such as cathepsin C, thioredoxin peroxidase, vitellogenin receptor precursor, 50S ribosomal protein L24 and ubiquitin-conjugating enzyme E2 isoform 2 were identified from this EST project and demonstrated as gonad differential expression genes by rqRT-PCR. Sixty five different types of SSRs (simple sequence repeats) were identified from the total 411 EST-SSR motifs. A home-made cDNA microarray containing 5664 spots was developed and the hybridization results indicated that 39 unique transcripts were differentially expressed in testis and ovaries (P<0.05). The expression levels of eleven unique transcripts examined by rqRT-PCR were matched with microarray fairly. These results will provide a useful resource for functional genomic studies on the biology of reproduction of mud crab.

Identification of testis-relevant genes using in silico analysis from testis ESTs and cDNA microarray in the black tiger shrimp (Penaeus monodon)

Background

Poor reproductive maturation of the black tiger shrimp (Penaeus monodon) in captivity is one of the serious threats to sustainability of the shrimp farming industry. Understanding molecular mechanisms governing reproductive maturation processes requires the fundamental knowledge of integrated expression profiles in gonads of this economically important species. In P. monodon, a non-model species for which the genome sequence is not available, expressed sequence tag (EST) and cDNA microarray analyses can help reveal important transcripts relevant to reproduction and facilitate functional characterization of transcripts with important roles in male reproductive development and maturation.

Results

In this study, a conventional testis EST library was exploited to reveal novel transcripts. A total of 4,803 ESTs were unidirectionally sequenced and analyzed in silico using a customizable data analysis package, ESTplus. After sequence assembly, 2,702 unique sequences comprised of 424 contigs and 2,278 singletons were identified; of these, 1,133 sequences are homologous to genes with known functions. The sequences were further characterized according to gene ontology categories (41% biological process, 24% molecular function, 35% cellular component). Through comparison with EST libraries of other tissues of P. monodon, 1,579 transcripts found only in the testis cDNA library were identified. A total of 621 ESTs have not been identified in penaeid shrimp. Furthermore, cDNA microarray analysis revealed several ESTs homologous to testis-relevant genes were more preferentially expressed in testis than in ovary. Representatives of these transcripts, homologs of saposin (PmSap) and Dmc1 (PmDmc1), were further characterized by RACE-PCR. The more abundant expression levels in testis than ovary of PmSap and PmDmc1 were verified by quantitative real-time PCR in juveniles and wild broodstock of P. monodon.

Conclusions

Without a genome sequence, a combination of EST analysis and high-throughput cDNA microarray technology can be a useful integrated tool as an initial step towards the identification of transcripts with important biological functions. Identification and expression analysis of saposin and Dmc1 homologs demonstrate the power of these methods for characterizing functionally important genes in P. monodon.

Molecular cloning of proliferating cell nuclear antigen and its differential expression analysis in the developing ovary and testis of penaeid shrimp Marsupenaeus japonicus

To understand the molecular mechanisms of gonadal development and maturation in penaeid shrimp Marsupenaeus japonicus, eight differentially expressed genes were obtained using a modified annealing control primer system. One of these genes is a proliferating cell nuclear antigen (PCNA). Bioinformatics analyses showed that full-length cDNA of M. japonicus PCNA (mjPCNA) consists of 75 bp of 5' untranslated region, 783 bp of coding region, and 65 bp of 3' untranslated region (excluding the polyA tail), encoding a protein of 260 amino acids with a predicted molecular mass of 28.85 kDa and an isoelectric point of 4.59. Real-time polymerase chain reaction analyses demonstrated that the gene expression level changed significantly in the developing testis and ovary. In stage 1 of ovary and testis, mjPCNA showed its lowest level during development and reached its highest expression level in stage 2 of ovary and testis. In stages 4 and 5 of ovary and the stage 3 of testis, mjPCNA held a steady expression level. Data suggest that PCNA plays an important role in the testis and ovary development, especially in the process of mitosis and meiosis.

EST analysis and identification of gonad-related genes from the normalized cDNA library of large yellow croaker, Larimichthys crocea

On grounds of the especially limited numbers of identified gonad-specific or gonad-related genes of large yellow croaker Larimichthys crocea which may represent a major obstacle for the study of gonad development and sex differentiation, we initiated a sequencing program of Expressed Sequence Tags (ESTs) in large yellow croaker. In this study, we firstly constructed a normalized gonad cDNA library using the combination of SMART technique and DSN treatment. The titer of amplified cDNA library was 4.8x10(11) and the percentage of unique cDNA sequences of the library was 82.49%. 2916 unique cDNAs were clustered from the 3535 high quality ESTs. Among the 1785 ESTs which had significant homology with known genes in the NCBI database, about 64 significant gonad-related genes were found, accounting for 3.59% of the total unique cDNAs. Specifically, the testis-specific LRR gene and testis-specific chromodomain Y-like protein gene were identified from fish for the first time. Six gonad-related microsatellite-containing ESTs were identified from the 129 ESTs containing 149 microsatellites. Expression patterns of 10 of these gonad-related gene homologues in ovaries and testes were examined by qRT-PCR. The results will be powerful resources for our further investigation to establish the molecular mechanisms of gonad development and sex differentiation in large yellow croaker.

cDNA cloning and expression of Ubc9 in the developing embryo and ovary of Oriental river prawn, Macrobrachium nipponense

The small ubiquitin-like modifier (SUMO) pathway in eukaryotes is an essential biological process involving cellular processes, development and organelle biogenesis. In a sequential enzymatic action, Ubc9 is an important conjunction enzyme in the SUMO pathway. Although the Ubc9 has been found in vertebrates, its expression in crustaceans is little known. In this study, the Ubc9 was identified in the embryo and ovary of a freshwater prawn Macrobrachium nipponense for the first time and it was denoted as MnUbc9. Bioinformatics analyses showed that this gene encodes a protein of 161 amino acids with predicted molecular mass of 18.32kDa. Real-time quantitative PCR analyses demonstrated that the expression levels varied significantly in the developing embryo and ovary. In the embryo, the expression level of MnUbc9 was higher at the cleavage stage (CS) than at the blastula stage (BS), and reached even higher levels at the protozoea stage (PS) and the zoea stage (ZS). In the ovary, the MuUbc9 expression was low at the early stage, but reached the highest at the yolk granule stage (YG), and then abruptly declined at the maturation stage (MA). The differential expressions of MnUbc9 in the embryo and ovary suggest that MnUbc9 may play an important role in embryogenesis and oogenesis of M. nipponense.