Identification, characterization and expression of sex-related genes in testes of the giant tiger shrimp Penaeus monodon

Let's talk about sex: Mechanisms of neural sexual differentiation in Bilateria

In multicellular organisms, sexed gonads have evolved that facilitate release of sperm versus eggs, and bilaterian animals purposefully combine their gametes via mating behaviors. Distinct neural circuits have evolved that control these physically different mating events for animals producing eggs from ovaries versus sperm from testis. In this review, we will describe the developmental mechanisms that sexually differentiate neural circuits across three major clades of bilaterian animals-Ecdysozoa, Deuterosomia, and Lophotrochozoa. While many of the mechanisms inducing somatic and neuronal sex differentiation across these diverse organisms are clade-specific rather than evolutionarily conserved, we develop a common framework for considering the developmental logic of these events and the types of neuronal differences that produce sex-differentiated behaviors. This article is categorized under: Congenital Diseases > Stem Cells and Development Neurological Diseases > Stem Cells and Development.

Differentially expressed genes in the testes from early to mature development of banana shrimp (Fenneropenaeus merguiensis)

Banana shrimp (Fenneropenaeus merguiensis) is an economically important species in Thailand owing to the high value of globally exported frozen brine shrimps. However, the regulatory mechanisms governing spermatogenesis and testicular development in this species are poorly understood. High-throughput RNA sequencing was used to investigate the mechanisms and regulated genes involved in testis development using transcriptome profiling of juvenile and adult banana shrimp testes. Differentially expressed genes (DEGs) in these two libraries were identified and quantified to confirm gene expression. DEGs were found in 7,347 genes, with 4,465 upregulated and 2,882 downregulated. Some of these genes were designated as candidate genes, and six specific DEGs, including PRM1, SPATA20, Sry, SSRF, Sxl, and Tra-2c, were selected to confirm the reliability of the RNA-seq data using qPCR. Moreover, six non-DEGs were chosen based on testis-specific and regulatory genes that support a specific function in spermatogenesis and testis development in this species, including Dsx, Gfra2, IAG, Sox9, Sox13, and Sox14A. Furthermore, Sry, Sox14A, Sox14B and SPATA20 were identified in early stages (nauplius-postlarvae) of shrimp development to provide more information involving testes formation and development. The transcript data from this study could differentiate a group of genes required at the early and late stages of testis development and both sets of testis development. Therefore, this information would help in manipulating each stage of testicular development.

IAG Regulates the Expression of Cytoskeletal Protein-Encoding Genes in Shrimp Testis

Insulin-like androgenic gland hormone (IAG) is the master regulator of sexual differentiation and testis development in male crustaceans. However, the molecular mechanism on how IAG functions during testis development is still largely unknown. Here, the transcriptional changes were analyzed in the testes of shrimp after LvIAG knockdown in Litopenaeus vannamei. Differential expression analysis identified 111 differentially expressed genes (DEGs), including 48 upregulated DEGs and 63 downregulated DEGs, in testes of shrimp after LvIAG knockdown. Gene ontology (GO) analysis showed that these DEGs were apparently enriched in cytoskeleton-related GO items. Gene function analysis showed that genes enriched in these GO items mainly encoded actin, myosin, and heat shock protein. Interestingly, these genes were all downregulated in testis after LvIAG knockdown, which was confirmed by qRT-PCR detection. Furthermore, injection of LvIAG protein that was recombinantly expressed in insect cells upregulated the expression levels of these genes. The present study revealed that shrimp IAG might function in testis development through regulating the expression of cytoskeletal protein-encoding genes, which would provide new insights into understanding the functional mechanisms of IAG on male sexual development of crustaceans.

Application of proteomics in shrimp and shrimp aquaculture

Since proteins play an important role in the life of an organism, many researchers are now looking at how genes and proteins interact to form different proteins. It is anticipated that the creation of adequate tools for rapid analysis of proteins will accelerate the determination of functional aspects of these biomolecules and develop new biomarkers and therapeutic targets for the diagnosis and treatment of various diseases. Though shrimp contains high-quality marine proteins, there are reports about the heavy losses to the shrimp industry due to the poor quality of shrimp production and many times due to mass mortality also. Frequent outbreaks of diseases, water pollution, and quality of feed are some of the most recognized reasons for such losses. In the seafood export market, shrimp occupies the top position in currency earnings and strengthens the economy of many developing nations. Therefore, it is vital for shrimp-producing companies they produce healthy shrimp with high-quality protein. Though aquaculture is a very competitive market, global awareness regarding the use of scientific knowledge and emerging technologies to obtain better-farmed organisms through sustainable production has enhanced the importance of proteomics in seafood biology research. Proteomics, as a powerful tool, has therefore been increasingly used to address several issues in shrimp aquaculture. In the present paper, efforts have been made to address some of them, particularly the role of proteomics in reproduction, breeding and spawning, immunological responses and disease resistance capacity, nutrition and health, microbiome and probiotics, quality and safety of shrimp production, bioinformatics applications in proteomics, the discovery of protein biomarkers, and mitigating biotic and abiotic stresses. Future challenges and research directions on proteomics in shrimp aquaculture have also been discussed.

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

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

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

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

Association between ALFPm3 single nucleotide polymorphism and white spot syndrome virus resistance in black tiger shrimp Penaeus monodon

Here single nucleotide polymorphisms (SNPs) were associated with white spot syndrome virus (WSSV) resistance in black tiger shrimp Penaeus monodon. SNPs were identified by single-strand conformation polymorphism (SSCP) screening and DNA sequencing of shrimp sampled from 3 families (100 shrimp per family) challenged with WSSV. Shrimp that died over the 14 d challenge trial were designated susceptible, with those remaining alive on Day 14 designated resistant. To compare SNPs, 10 samples from the susceptible and resistant groups, each comprising DNA pooled from 3 shrimp, were amplified by polymerase chain reaction (PCR) using primers to 12 selected genes and screened by SSCP. SNPs were only identified in the anti-lipopolysaccharide factor 3 (ALFPm3) gene product. Analysis of complete ALFPm3 gene sequences confirmed the existence of 3 SNPs (g.934C>G, g.1186A>G, and g.1898C>G) that were polymorphic between the susceptible and resistant groups. Further analyses using specific tetra-primer amplification refractory mutation system PCR primer sets associated these 3 SNPS, and particularly the g.1186A>G SNP, with WSSV resistance. This SNP thus has potential for use as a DNA marker to select for WSSV resistance in P. monodon breeding programs.

Penaeus stylirostris densovirus proteins CP and NS1 interact with peritrophin of Litopenaeus vannamei

The Penaeus stylirostris densovirus (PstDNV) is a major virus of shrimps that severely harms the shrimp farming industry. Peritrophin is a peritrophic membrane protein with chitin binding activity. To examine the roles of peritrophin in viral infection, we used yeast two-hybrid to analyze the interaction between the Pacific white shrimp (Litopenaeus vannamei) peritrophin and PstDNV proteins (CP, NS1 and NS2). The yeast two-hybrid results showed that NS1 and peritrophin had an interaction, CP and peritrophin had an interaction as well, and NS2 had no interaction with peritrophin. We validated the interactions with GST pull-down assays. We then conducted RNA interference and qRT-PCR. The results showed that when pre-injection of dsRNA-peritrophin, the quantity of PstDNV in the shrimps injected with viruses was significantly lower than in the control group (P < 0.01), indicating the viral infection was decreased when the peritrophin gene expression was inhibited. The results indicated that peritrophin of L. vannamei participated in the PstDNV infection.

Molecular Characterization and Expression Profiling of Three Transformer-2 Splice Isoforms in the Redclaw Crayfish, Cherax quadricarinatus

Sex determination/sex differentiation is determined by genetics, environmental factors, or the interactions of the two. The Transformer-2 (Tra-2) gene plays an important role in the sex determination cascade signal pathway in insects. In this study, the Tra-2 gene was isolated and characterized from the cDNA library of gonad tissues in the redclaw crayfish, Cherax quadricarinatus. Three splice variants were identified, designated as CqTra-2A, CqTra-2B, and CqTra-2C, and sequence analysis showed that they had a highly conserved RRM domain. Phylogenetic analysis was performed by the NJ method, and the results revealed that the Tra-2 protein of the redclaw crayfish was very closely related to those of Macrobrachium rosenbergii, Fenneropenaeus chinensis, and Macrobrachium nipponense. Real-time PCR analysis showed that the three isoforms were predominantly expressed in the ovary and gradually increased with embryonic development. Additionally, the expression pattern of CqTra-2 at different developmental stages was analyzed by qPCR and revealed that the phase of having a body length of 3 cm may be the key period for the sex differentiation of C. quadricarinatus. RNAi-targeting gene silencing further confirmed the function of CqTra-2 in sexual differentiation in redclaw crayfish. Our experimental data will contribute to understanding the mechanism of sex determination in crustaceans.

Molecular cloning, expression pattern analysis, and in situ hybridization of a Transformer-2 gene in the oriental freshwater prawn, Macrobrachium nipponense (de Haan, 1849)

In this study, we isolated a full-length cDNA sequence from Macrobrachium nipponense and investigated its gene function. We named the gene Mntra-2a because of high similarities and close evolutionary divergence with arthropod tra-2. The full-length cDNA of Mntra-2a was 1293 bp, consisting of a 212 bp 5' UTR, a 268 bp 3' UTR, and an ORF of 813 bp encoding 270 amino acids. It contained an RNA recognition motif and a linker region. Real-time PCR analysis showed that Mntra-2a was highly expressed in the gonads of both males and females. Further in situ hybridization analysis showed that Mntra-2a was mainly located in oocytes and spermatocytes. During embryogenesis, Mntra-2a expression was higher in the cleavage and nauplius stages. During the ovarian reproductive cycle, Mntra-2a expression reached a peak at OvaryV and decreased to the lowest level at OvaryIV. These results indicated that Mntra-2a probably played important roles in embryonic development and early gonad development in M. nipponense. Our results provide basic information for further functional studies of tra-2 in M. nipponense.

A natural antisense transcript of the fem-1 gene was found expressed in female gonads during the characterization, expression profile, and cellular localization of the fem-1 gene in Pacific white shrimp Penaeus vannamei

The fem-1 gene in Caenorhabditis elegans is involved in sex differentiation; it is specifically required for all aspects of male development. In this study, the full-length cDNA of the fem-1 (Pvfem-1) gene was isolated from the Pacific whiteleg shrimp Penaeus vannamei. The Pvfem-1 transcript is 3778 nt long and encodes a putative protein (PvFEM-1) of 638 amino acids that presented eight ankyrin repeats. The translated protein showed a significant (P < 0.05) structural similitude by superposition with C. elegans FEM-1 protein. Pvfem-1 expression was evaluated by qPCR and in situ hybridization (ISH) during embryogenesis, larval development, and gonads of both genders in subadult and adult life stages. Pvfem-1 was found expressed in brain, intestine, hepatopancreas, and in the gonads of both genders in subadults and adults when quantified by RT-qPCR. A significant finding was the discovery of a natural antisense transcript (NAT) of Pvfem-1 by ISH. It was present in the oocyte nucleus of subadult female shrimp gonads but was not seen within oocytes from adult females, although it was detected in follicular cells, suggesting a possible post-transcriptional regulation of Pvfem-1 in female gonad. Conversely, in males, no NAT was observed, and Pvfem-1 was found expressed in spermatogonia of both, subadult and adult shrimps indicating a function in male sexual differentiation and gametes generation. This study represents the first step for future functional analysis that is expected to contribute to clarifying the role of Pvfem-1 in sex differentiation and determination.

The State of "Omics" Research for Farmed Penaeids: Advances in Research and Impediments to Industry Utilization

Elucidating the underlying genetic drivers of production traits in agricultural and aquaculture species is critical to efforts to maximize farming efficiency. "Omics" based methods (i.e., transcriptomics, genomics, proteomics, and metabolomics) are increasingly being applied to gain unprecedented insight into the biology of many aquaculture species. While the culture of penaeid shrimp has increased markedly, the industry continues to be impeded in many regards by disease, reproductive dysfunction, and a poor understanding of production traits. Extensive effort has been, and continues to be, applied to develop critical genomic resources for many commercially important penaeids. However, the industry application of these genomic resources, and the translation of the knowledge derived from "omics" studies has not yet been completely realized. Integration between the multiple "omics" resources now available (i.e., genome assemblies, transcriptomes, linkage maps, optical maps, and proteomes) will prove critical to unlocking the full utility of these otherwise independently developed and isolated resources. Furthermore, emerging "omics" based techniques are now available to address longstanding issues with completing keystone genome assemblies (e.g., through long-read sequencing), and can provide cost-effective industrial scale genotyping tools (e.g., through low density SNP chips and genotype-by-sequencing) to undertake advanced selective breeding programs (i.e., genomic selection) and powerful genome-wide association studies. In particular, this review highlights the status, utility and suggested path forward for continued development, and improved use of "omics" resources in penaeid aquaculture.

Identification of MicroRNAs and Their Target Genes Associated with Ovarian Development in Black Tiger Shrimp (Penaeus monodon) Using High-Throughput Sequencing

Plenty of evidence showing that microRNAs (miRNAs) post-transcriptionally regulate gene expression and are involved in a wide range of biological processes. However, the roles of miRNAs in ovarian development process remain largely unknown in shrimp. In the present study, high-throughput sequencing of small RNAs was performed to find specific miRNAs that are involved in ovarian development process in Penaeus monodon. Two small RNA libraries were constructed from undeveloped (UNDEV group) and developed (DEV group) ovarian tissues in P. monodon. In total, 43 differentially expressed miRNAs were identified between the two groups (P ≤ 0.05, |log2 ratio| ≥1), and their expression profiles were validated by qRT-PCR. In order to further clarify the functional roles of these differentially expressed miRNAs during ovarian development process, target gene prediction was performed. In total, 4,102 target genes of 43 miRNAs were predicted, then clustered by the Kyoto Encyclopedia of Genes and Genomes (KEGG) database; only four specific pathways related to ovarian development were obtained (P < 0.05). Dual-luciferase reporter assays and integrated expression analysis were also conducted to further clarify the interaction between the miRNAs and their target mRNAs. This study provides important information about the function of miRNAs involved in ovarian developmental stages in P. monodon.

Characterization of transforming growth factor beta regulator 1-like and association between its expression levels and growth of the giant tiger shrimp Penaeus monodon

Molecular markers that allow selection of juveniles and broodstock with improved growth performances are useful for the shrimp industry. Here, the full-length cDNA of transforming growth factor beta regulator 1-like (PmTbrg1-l) in the giant tiger shrimp Penaeus monodon was determined. It was 1184 bp in length and contained an open reading frame (ORF) of 975 bp corresponding to a deduced polypeptide of 324 amino acids. Successful RNA interference (RNAi) carried out using juveniles injected with PmTbrg1-l dsRNA revealed reduced levels of PmTbrg1-l and myostatin (PmMstm) in hemocytes when compared to shrimp injected with saline solution and GFP dsRNA (P < .05). Associations between single-strand conformational polymorphism (SSCP) patterns or single nucleotide polymorphism (SNP) patterns and growth-related parameters (average body weight and total length) were examined. Juveniles with pattern III (corresponding to A/A₉₁₈; N = 37) showed a trend for greater average body weight and total length than those with patterns II (G/G₉₁₈; N = 42) and IV (A/G₉₁₈; N = 75). The expression level of PmTbrg1-l in the hepatopancreas of females was significantly higher than that in males (P < .05) in two sample sets of three-month-old domesticated juveniles (N = 59 and 50; P < .05). Moreover, its expression level in large-size juveniles was significantly higher than that in medium-size and small-size juveniles in both groups of samples (P < .05). Results indicated that PmTbrg1-l is functionally related with growth of P. monodon.

Isolation of the sex-determining gene Sex-lethal (Sxl) in Penaeus (Litopenaeus) vannamei (Boone, 1931) and characterization of its embryogenic, gametogenic, and tissue-specific expression

The Pacific white shrimp Penaeus vannamei is the most cultured shrimp species around the world. Because females grow larger than males, the culture of 'only females' is of great interest, but knowledge on sex determination and differentiation is required for producing only females. In an effort to obtain information associated with reproduction in P. vannamei, transcriptomic data from female gonads was generated, and partial sequences of a transcript were identified as Sex-lethal (Sxl). Its characterization indicated that, differently from other penaeids in which this gene has been isolated, there are six isoforms of the Sxl transcript in P. vannamei (PvanSxl 1-6). These isoforms result from alternative splicing at three splice sites (SS1, SS2, SS3). The first splice-site is unique to P. vannamei, as it has not been reported for other Arthropod species; the second splice-site (SS2) is common among crustaceans, and the third splice-site (SS3) is also unique to P. vannamei and when spliced-out, it is always together with SS2. All isoforms are expressed during embryogenesis as well as gametogenesis of both genders. The two shorter isoforms, PvanSxl-5 and PvanSxl-6, which result from the splicing of SS2 and SS3, were found mostly expressed in adult testis, but PvanSxl-6 was also expressed in oocytes during gametogenesis. During oogenesis, the second largest isoform, PvanSxl-2, which splices-out only SS1, and PvanSxl-4 that splices-out SS1 and SS2 were highly expressed. These two isoforms were also highly expressed during embryonic development. In situ hybridization allowed pinpointing more specifically the cells where the PvanSxl transcripts were expressed. During embryogenesis, hybridization was observed from the one-cell stage embryo to late gastrula. In the female gonad in previtellogenesis, hybridization occurred in the nucleus of oocytes, whereas in secondary vitellogenesis the transcript also hybridized cytoplasmic granules and cortical crypts. Finally, in situ hybridization corroborated the expression of PvanSxl also in the male gonad during spermatogenesis, mostly occurring in the cytoplasm from spermatogonia and spermatocytes.

De novo transcriptome analysis and differentially expressed genes in the ovary and testis of the Japanese mantis shrimp Oratosquilla oratoria by RNA-Seq

The mantis shrimp Oratosquilla oratoria is a widely distributed, commercially important crustacean species. Although its conservation and the development of successful artificial breeding technologies have recently received considerable attention, there are currently no available data regarding the molecular mechanisms in controlling reproduction. In this study, we performed transcriptome sequencing of the testis, ovary, female and male eyestalks and the androgenic gland of O. oratoria, and compared the expression pattern of transcripts from the testis and ovary libraries to identify genes involved in gonadal development. A total of 147,130,937 clean reads were retrieved after removing the adapters in reads and filtering out low-quality data. All the reads were assembled into 94,990 unigenes (23,133 in testis and ovary) with an average length of 783 base pairs (bp) and N50 of 1502 bp. A search of all-unigenes against COG, GO, KEGG, KOG, Pfam, Swiss-Prot and Nr databases resulted in a total of 19,404 annotated unigenes. Comparison of the sequences in the ovary and testis libraries revealed that 1188 unigenes were up-regulated in the ovary and 2732 were up-regulated in the testis. Twenty ovary-up-regulated and 21 testis-up-regulated unigenes were confirmed by quantitative real-time PCR. Additionally, 13,437 simple sequence repeats (SSRs) and 275,799 putative single nucleotide polymorphisms (SNPs) were identified. The important functional genes and pathways identified here provide a valuable dataset for understanding the molecular mechanisms controlling gonad development in O. oratoria, and the numerous (13,437 SSRs and 275,799 SNPs) molecular markers obtained here will provide fundamental basis for functional genomic and population genetic studies of O. oratoria.

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

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

Molecular characterization and function of the Prohibitin2 gene in Litopenaeus vannamei responses to Vibrio alginolyticus

Prohibitin2 (PHB2), a potential tumor suppressor protein, plays important roles in inhibition of cell cycle progression, transcriptional regulation, apoptosis and the mitochondrial respiratory chain. To explore its potential roles in crustaceans' immune responses we have identified and characterized LvPHB2, a 891 bp gene encoding a 297 amino acids protein in the shrimp Litopenaeus vannamei. Expression analyses showed that LvPHB2 is expressed in all examined tissues, and largely present in cytoplasm, correlating with its known anti-oxidation function in mitochondria. Luciferase reporter assays showed that over-expression of LvPHB2 could activate the p53 pathway, indicating that it might participate in apoptosis regulation. Quantitative real-time PCR revealed that infection with Vibrio alginolyticus induces its up-regulation in hepatopancreas. Moreover, RNAi knock-down of LvPHB2 in vivo raises mortality rates of L. vannamei infected by V. alginolyticus, and affects expression of STAT3, Caspase3 and p53 genes. We found significantly higher reactive oxygen species production, DNA damage and apoptosis rates in LvPHB2-silenced shrimp challenged with V. alginolyticus than in controls injected with a Green Fluorescent Protein-silencing construct. Our results suggest that LvPHB2 plays a vital role in shrimp responses to V. alginolyticus infection through its participation in regulation of oxidants and apoptosis.

Identification of differentially expressed reproductive and metabolic proteins in the female abalone (Haliotis laevigata) gonad following artificial induction of spawning

Inefficient control of temperate abalone spawning prevents pair-wise breeding and production of abalone with highly marketable traits. Traditionally, abalone farmers have used a combination of UV irradiation and application of temperature gradients to the tank water to artificially induce spawning. Proteins are known to regulate crucial processes such as respiration, muscle contraction, feeding, growth and reproduction. Spawning as a pre-requisite of abalone reproduction is likely to be regulated, in part, by endogenous proteins. A first step in elucidating the mechanisms that regulate spawning is to identify which proteins are directly involved during spawning. The present study examined protein expression following traditional spawning induction in the Haliotis laevigata female. Gonads were collected from abalone in the following physiological states: (1) spawning; (2) post-spawning; and (3) failed-to-spawn. Differential protein abundance was initially assessed using two-dimensional difference in-gel electrophoresis coupled with mass spectrometry for protein identification. A number of reproductive proteins such as vitellogenin, vitelline envelope zona pellucida domain 29 and prohibitin, and metabolic proteins such as thioredoxin peroxidase, superoxide dismutase and heat shock proteins were identified. Differences in protein abundance levels between physiological states were further assessed using scheduled multiple reaction monitoring mass spectrometry. Positive associations were observed between the abundance of specific proteins, such as heat shock cognate 70 and peroxiredoxin 6, and the propensity or failure to spawn in abalone. These findings have contributed to better understand both the effects of oxidative and heat stress over abalone physiology and their influence on abalone spawning.

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

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

Male Sexual Development and the Androgenic Gland: Novel Insights through the de novo Assembled Transcriptome of the Eastern Spiny Lobster, Sagmariasus verreauxi

The Eastern spiny lobster, Sagmariasus verreauxi, is commercially important in fisheries, with growing aquaculture potential, driving an interest to better understand male sexual differentiation. Amongst the Decapoda, the androgenic gland (AG) and the insulin-like androgenic gland hormone (IAG) have a well-defined function in male sexual differentiation. However, IAG is not a sex determinant and therefore must be considered as part of a broader, integrated pathway. This work uses a transcriptomic, multi-tissue approach to provide an integrated description of male-biased expression as mediated through the AG. Transcriptomic analyses demonstrate that IAG expression is stage- and eyestalk-regulated (low in immature, high in mature and 6-times higher in hypertrophied glands), with IAG being the predominant AG-specific factor. The low expression of this key factor in immature males suggests the involvement of other tissues in male sexual differentiation. Across tissues, the gonad (87.8%) and antennal gland (73.5%) show the highest male-biased differential expression of transcripts and also express 4 sex-determination regulators, known as Dmrts, with broader expression of Sv-Sxl and Sv-TRA-2. In order to better understand male sexual differentiation, tissues other than the AG must also be considered. This research highlights the gonad and antennal gland as showing significant male-biased expression patterns, including the Sv-Dmrts.

Cellular and molecular markers in monitoring the fate of lymphoid cell culture from Penaeus monodon Fabricius (1798)

Lymphoid cell culture from penaeid shrimps has gained much acceptance as an in vitro platform to facilitate research on the development of prophylaxis, and therapeutic strategies against viruses and for cell line development. However, lymphoid cells can be used as platform for in vitro research, only if they are in metabolically and mitotically active state in vitro with unaltered cell surface receptors. Through this study, we addressed the response of lymphoid cells to a new microenvironment at cellular and molecular levels; including the study of mitotic events, DNA synthesis, expression profile of cell cycle genes, cytoskeleton organization, metabolic activity and viral susceptibility. The S-phase entry and synthesis of new DNA was recorded by immunoflourescent technique. Cdc2, CycA, CycB, EF-1α and BUB3 genes involved in cell cycle were studied in both the cells and tissue, of which EF-1α showed an elevated expression in cells in vitro (∼ 19.7%). Cytoskeleton network of the cell was examined by studying the organization of actin filaments. As the markers for metabolic status, mitochondrial dehydrogenase, protein synthesis and glucose assimilation by the cells were also assessed. Viral susceptibility of the cell was determined using WSSV to confirm the preservation of cellular receptors. This study envisages to strengthen the shrimp cell line research and to bring forth lymphoid cell culture system as a 'model' in vitro system for shrimp and crustaceans altogether.

Gonadal transcriptomic analysis and differentially expressed genes in the testis and ovary of the Pacific white shrimp (Litopenaeus vannamei)

Background

The Pacific white shrimp (Litopenaeus vannamei) is the world’s most prevalent cultured crustacean species. However, the supply of high-quality broodstocks is limited and baseline information related to its reproductive activity and molecular issues related to gonad development are scarce. In this study, we performed transcriptome sequencing on the gonads of adult male and female L. vannamei to identify sex-related genes.

Results

A total of 25.16 gigabases (Gb) of sequences were generated from four L. vannamei gonadal tissue libraries. After quality control, 24.11 Gb of clean reads were selected from the gonadal libraries. De-novo assembly of all the clean reads generated a total of 65,218 unigenes with a mean size of 1021 bp and a N50 of 2000 bp. A search of all-unigene against Nr, SwissProt, KEGG, COG and NT databases resulted in 26,482, 23,062, 20,659, 11,935 and 14,626 annotations, respectively, providing a total of 30,304 annotated unigenes. Among annotated unigenes, 12,320 unigenes were assigned to gene ontology categories and 20,659 unigenes were mapped to 258 KEGG pathways. By comparing the ovary and testis libraries, 19,279 testicular up-regulated and 3,529 ovarian up-regulated unigenes were identified. Enrichment analysis of differentially expressed unigenes resulted in 1060 significantly enriched GO terms and 34 significantly enriched KEGG pathways. Nine ovary-specific, 6 testis-specific, 45 testicular up-regulated and 39 ovarian up-regulated unigenes were then confirmed by semi-quantitative PCR and quantitative real-time PCR. In addition, using all-unigenes as a reference, a total of 13,233 simple sequence repeats (SSRs) were identified in 10,411 unigene sequences.

Conclusions

The present study depicts the first large-scale RNA sequencing of shrimp gonads. We have identified many important sex-related functional genes, GO terms and pathways, all of which will facilitate future research into the reproductive biology of shrimp. We expect that the SSRs detected in this study can then be used as genetic markers for germplasm evaluation of breeding and imported populations.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2219-4) contains supplementary material, which is available to authorized users.

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.

Transcriptome profiles of Penaeus (Marsupenaeus) japonicus animal and vegetal half-embryos: identification of sex determination, germ line, mesoderm, and other developmental genes

There is virtually no knowledge of the molecular events controlling early embryogenesis in Penaeid shrimp. A combination of controlled spawning environment, shrimp embryo micro-dissection techniques, and next-generation sequencing was used to produce transcriptome EST datasets of Penaeus japonicus animal and vegetal half-embryos. Embryos were collected immediately after spawning, and then blastomeres were separated at the two-cell stage and allowed to develop to late gastrulation, then pooled for RNA isolation and cDNA synthesis. Ion Torrent sequencing of cDNA from approximately 500 pooled animal and vegetal half-embryos from multiple spawnings resulted in 560,516 and 493,703 reads, respectively. Reads from each library were assembled and Gene Ontogeny analysis produced 3479 annotated animal contigs and 4173 annotated vegetal contigs, with 159/139 hits for developmental processes in the animal/vegetal contigs, respectively. Contigs were subject to BLAST for selected developmental toolbox genes. Some of the genes found included the sex determination genes sex-lethal and transformer; the germ line genes argonaute 1, boule, germ cell-less, gustavus, maelstrom, mex-3, par-1, pumilio, SmB, staufen, and tudor; the mesoderm genes brachyury, mef2, snail, and twist; the axis determination/segmentation genes β-catenin, deformed, distal-less, engrailed, giant, hairy, hunchback, kruppel, orthodenticle, patched, tailless, and wingless/wnt-8c; and a number of cell-cycle regulators. Animal and vegetal contigs were computationally subtracted from each other to produce sets unique to either half-embryo library. Genes expressed only in the animal half included bmp1, kruppel, maelstrom, and orthodenticle. Genes expressed only in the vegetal half included boule, brachyury, deformed, dorsal, engrailed, hunchback, spalt, twist, and wingless/wnt-8c.

Identification of genes involved in reproduction and lipid pathway metabolism in wild and domesticated shrimps

The aims of this study were to identify genes involved in reproduction and lipid pathway metabolism in Penaeus monodon and correlate their expression with reproductive performance. Samples of the hepatopancreas and ovaries were obtained from a previous study of the reproductive performance of wild and domesticated P. monodon broodstock. Total mRNA from the domesticated broodstock was used to create two next generation sequencing cDNA libraries enabling the identification of 11 orthologs of key genes in reproductive and nutritional metabolic pathways in P. monodon. These were identified from the library of de novo assembled contigs, including the description of 6 newly identified genes. Quantitative RT-PCR of these genes in the hepatopancreas prior to spawning showed that the domesticated mature females significantly showed higher expression of the Pm Elovl4, Pm COX and Pm SUMO genes. The ovaries of domesticated females had a significantly decreased expression of the Pm Elovl4 genes. In the ovaries of newly spawned females, a significant correlation was observed between hepatosomatic index and the expression of Pm FABP and also between total lipid content and the expression of Pm CYP4. Although not significant, the highest levels of correlation were found between relative fecundity and Pm CRP and Pm CYP4 expression, and between hatching rate and Pm Nvd and Pm RXR expression. This study reports the discovery of genes involved in lipid synthesis, steroid biosynthesis and reproduction in P. monodon. These results indicate that genes encoding enzymes involved in lipid metabolism pathways might be potential biomarkers to assess reproductive performance.

RNAi-mediated knockdown of transformer-2 in the predatory mite Metaseiulus occidentalis via oral delivery of double-stranded RNA

Little is known about the process of sex determination at the molecular level in Metaseiulus occidentalis, a parahaploid species and natural enemy of phytophagous pest mites. Detailed knowledge of the sex-determination pathway could allow genetic manipulation of M. occidentalis to produce more female offspring, which could improve its effectiveness as a biological control agent. RNA interference is useful for assessing the function of putative sex-determination genes by reducing or eliminating gene expression. In many insect species the transformer-2 (tra-2) gene is an upstream regulatory element in the sex-determination cascade, and knockdown of tra-2 expression can alter the sex ratio. We assessed whether oral delivery of tra-2 double-stranded RNA to M. occidentalis virgin females would affect the sex of her progeny. Females that ingested tra-2 dsRNA produced significantly fewer eggs compared to control females suggesting that tra-2 is somehow involved in reproduction by females. However, the sex ratio of the few progeny that were laid was not altered, so it is unclear whether tra-2 is involved in sex determination. This is an initial step towards elucidating the molecular components of sex determination in M. occidentalis.

Identification of genes associated with reproduction in the Mud Crab (Scylla olivacea) and their differential expression following serotonin stimulation

The central nervous system (CNS) is often intimately involved in reproduction control and is therefore a target organ for transcriptomic investigations to identify reproduction-associated genes. In this study, 454 transcriptome sequencing was performed on pooled brain and ventral nerve cord of the female mud crab (Scylla olivacea) following serotonin injection (5 µg/g BW). A total of 197,468 sequence reads was obtained with an average length of 828 bp. Approximately 38.7% of 2,183 isotigs matched with significant similarity (E value < 1e⁻⁴) to sequences within the Genbank non-redundant (nr) database, with most significant matches being to crustacean and insect sequences. Approximately 32 putative neuropeptide genes were identified from nonmatching blast sequences. In addition, we identified full-length transcripts for crustacean reproductive-related genes, namely farnesoic acid o-methyltransferase (FAMeT), estrogen sulfotransferase (ESULT) and prostaglandin F synthase (PGFS). Following serotonin injection, which would normally initiate reproductive processes, we found up-regulation of FAMeT, ESULT and PGFS expression in the female CNS and ovary. Our data here provides an invaluable new resource for understanding the molecular role of the CNS on reproduction in S. olivacea.

Transcriptome analysis of red swamp crawfish Procambarus clarkii reveals genes involved in gonadal development

Background

The red swamp crawfish, Procambarus clarkii, has become one of the most economically important cultured species in China. Currently, little is known about the gonadal development of this species. Isolation and characterization of genes are an initial step towards understanding gonadal development of P. clarkii.

Results

Using the 454 pyrosequencing technology, we obtained a total of 1,134,993 high quality sequence reads from the crawfish testis and ovary libraries. We aimed to identify different genes with a potential role in gonad development. The assembly formed into 22,652 isotigs, distributed by GO analysis across 55 categories in the three ontologies, ‘molecular function’, ‘cellular component’, and ‘biological processes’. Comparative transcript analysis showed that 1,720 isotigs in the ovary were up-regulated and 2138 isotigs were down-regulated. Several gonad development related genes, such as vitellogenin, cyclin B, cyclin-dependent kinases 2, Dmc1 and ubiquitin were identified. Quantitative real-time PCR verified the expression profiles of 14 differentially expressed genes, and confirmed the reliability of the 454 pyrosequencing.

Conclusions

Our findings provide an archive for future research on gonadal development at a molecular level in P. clarkii and other crustacean. This data will be helpful to develop new ideas for artificial regulation of the reproductive process in crawfish aquaculture.

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 of sex-linked DNA markers by AFLP in Portunus trituberculatus

To investigate the mechanism of sex determination of crabs, amplified fragment-length polymorphism DNA analysis was applied to identify sex-specific markers for Portunus trituberculatus. Eleven pairs of primers were used to amplify DNA isolated from male and female crabs. A total of 481 and 499 fragments were amplified for females and males, respectively, and the ratios of polymorphic fragments were 20.59-65.52% and 27.27-62.22%, respectively. Distribution difference values between female and male groups for 10 polymorphisms were greater than 40%; Cluster analysis revealed that Nei's genetic distances for 57 polymorphisms clustered according to sex. The results indicate that a difference at the DNA level exists between female and male crabs and provide data for further study on sex determination.

Proteomic analysis of ovarian proteins and characterization of thymosin-β and RAC-GTPase activating protein 1 of the giant tiger shrimp Penaeus monodon

Cellular proteomics of total proteins in ovaries of domesticated and wild giant tiger shrimp (Penaeus monodon) were examined using GeLC-MS/MS. In total, 1638 proteins matched those previously deposited in databases and 1253 (76.50%) of these significantly matched known proteins. Several reproduction-related proteins (e.g. Cdc2, Cyclin B, Cdc25, 14-3-3, thymosin-β and Rac-GTPase activating protein 1) were identified. In addition, the full-length cDNA of P. monodon thymosin-β (PmTmsb; 1084 bp with an ORF of 387 bp and 128 deduced aa) and Rac-GTPase activating protein 1 (PmRacgap1; an ORF of 1881 bp and 626 deduced aa) were further characterized. PmTmsb was constitutively expressed in all tissues. In contrast, PmRacgap1 was more abundantly expressed in gonads than in several non-reproductive tissues (e.g. subcuticular epithelium, hepatopancreas, intestine, pleopods, stomach and thoracic ganglion). The expression levels of PmTmsb and PmRacgap1 in ovaries of wild adult broodstock were significantly greater than those in ovaries of juveniles (P<0.05). However, their expression levels did not vary significantly during ovarian development stages in intact broodstock. However, eyestalk ablation resulted in a significant reduction in PmTmsb expression at stages I and III ovaries (P<0.05), although it did not affect PmRacgap1 transcription significantly at these stages. On the other hand, use of polyclonal antibodies derived from recombinant PmTmsb and PmRacgap1 revealed that levels of both proteins decreased at the late stage (IV) of ovarian development. Our results suggested that PmTmsb and PmRacgap1 may act as negative effectors during ovarian development in P. monodon.

Characterization and expression of cell division cycle 2 (Cdc2) mRNA and protein during ovarian development of the giant tiger shrimp Penaeus monodon

The meiotic maturation of oocytes is regulated by the maturation-promoting factor (MPF), a complex of Cdc2 (Cdk1) and Cyclin B. Here, the complete open reading frame (ORF) of Cdc2 in Penaeus monodon was characterized. PmCdc2 were 900bp in length corresponding to a polypeptide of 299 amino acids with the conserved Thr14, Tyr15 and Thr161 residues. Quantitative real-time PCR indicated that the expression level of PmCdc2 in wild intact broodstock was significantly increased in stages II (vitellogenic) and III (early cortical rod) ovaries relative to stage I (previtellogenic) ovaries and peaked in stage IV (mature) ovaries (P<0.05). The expression level of PmCdc2 in stages I-IV ovaries of eyestalk-ablated broodstock was greater than that of the same ovarian developmental stages in intact broodstock (P<0.05). Expression levels of PmCdc2 in ovaries of 18-month-old P. monodon upon 5-HT injection (50μg/g body weight) were significantly increased at 1hour post injection (hpi, P<0.05). Recombinant PmCdc2 protein and its polyclonal antibody were successfully produced. Western blot analysis revealed the expected 34kDa band (PmCdc2) along with a smaller band of 23kDa (ribosomal protein S3) in ovaries of juveniles and various ovarian stages of broodstock. Using phospho-Cdc2 (Thr161) polyclonal antibody, the positive signal of 34kDa was observed in all ovarian stages but the most intense signal was found in stage IV ovaries. Results in the present study indicated that PmCdc2 gene/protein plays an important role in the development and maturation of oocytes/ovaries in P. monodon.

Identification of a sex-linked SNP marker in the salmon louse (Lepeophtheirus salmonis) using RAD sequencing

The salmon louse (Lepeophtheirus salmonis (Krøyer, 1837)) is a parasitic copepod that can, if untreated, cause considerable damage to Atlantic salmon (Salmo salar Linnaeus, 1758) and incurs significant costs to the Atlantic salmon mariculture industry. Salmon lice are gonochoristic and normally show sex ratios close to 1:1. While this observation suggests that sex determination in salmon lice is genetic, with only minor environmental influences, the mechanism of sex determination in the salmon louse is unknown. This paper describes the identification of a sex-linked Single Nucleotide Polymorphism (SNP) marker, providing the first evidence for a genetic mechanism of sex determination in the salmon louse. Restriction site-associated DNA sequencing (RAD-seq) was used to isolate SNP markers in a laboratory-maintained salmon louse strain. A total of 85 million raw Illumina 100 base paired-end reads produced 281,838 unique RAD-tags across 24 unrelated individuals. RAD marker Lsa101901 showed complete association with phenotypic sex for all individuals analysed, being heterozygous in females and homozygous in males. Using an allele-specific PCR assay for genotyping, this SNP association pattern was further confirmed for three unrelated salmon louse strains, displaying complete association with phenotypic sex in a total of 96 genotyped individuals. The marker Lsa101901 was located in the coding region of the prohibitin-2 gene, which showed a sex-dependent differential expression, with mRNA levels determined by RT-qPCR about 1.8-fold higher in adult female than adult male salmon lice. This study’s observations of a novel sex-linked SNP marker are consistent with sex determination in the salmon louse being genetic and following a female heterozygous system. Marker Lsa101901 provides a tool to determine the genetic sex of salmon lice, and could be useful in the development of control strategies.

Identification of differentially expressed genes in American cockroach ovaries and testes by suppression subtractive hybridization and the prediction of its miRNAs

Studies on the cockroach have contributed to our understanding of several important developmental processes, especially those that can be easily studied in the embryo. However, our knowledge on late events such as gonad differentiation in the cockroach is still limited. The major aim of the present study was to identify sex-specific genes between adult female and male Periplaneta americana. Two cDNA libraries were constructed using the suppression subtractive hybridization method; a total of 433 and 599 unique sequences were obtained from the forward library and the reverse library, respectively, by cluster assembly, and sequence alignment of 1,032 expressed sequence tags. The analysis of the differentially expressed gene functions allowed these genes to be categorized into three groups: biological process, molecular function, and cellular component. The differentially expressed genes were suggested to be related to the development of the gonads of P. americana. Twelve differentially expressed genes were randomly selected and verified using relative quantitative real-time polymerase chain reaction (qRT-PCR). Meanwhile, by adopting a range of filtering criteria, we predicted two potential microRNA sequences for P. americana, pam-miR100-3p and pam-miR7. To confirm the expression of potential microRNAs (miRNAs) in American cockroach, a qRT-PCR approach was also employed. The data presented here offer the insights into the molecular foundation of sex differences in American cockroach, and the first report for the miRNAs in this species. In addition, the results can be used as a reference for unraveling candidate genes associated with the sex and reproduction of cockroaches.

Validation of reference genes for real-time PCR of reproductive system in the black tiger shrimp

Gene expression of reproductive system of the black tiger shrimp (Peneaus monodon) has been widely studied to address poor maturation problem in captivity. However, a systematic evaluation of reference genes in quantitative real-time PCR (qPCR) for P. monodon reproductive organs is lacking. In this study, the stability of four potential reference genes (18s rRNA, GAPDH, β-actin, and EF1-α) was examined in the reproductive tissues in various conditions using bioinformatic tools: NormFinder and geNorm. For NormFinder, EF1-α and GAPDH ranked first and second as the most stable genes in testis groups whereas GAPDH and EF1-α were for ovaries from wild-caught broodstock and domesticated groups. EF1-α and β-actin ranked first and second for the eyestalk ablated ovaries. For geNorm, EF1-α and GAPDH had the best stability in all testis and ovaries from domesticated groups whereas EF1-α and β-actin were the best for ovaries from wild-caught and eyestalk ablated groups. Moreover, the expression levels of two well-known reproductive genes, Dmc1 and Vitellogenin, were used to validate these reference genes. When normalized to EF1-α, the expected expression patterns were obtained in all cases. Therefore, this work suggests that EF1-α is more versatile as reference genes in qPCR analysis for reproductive system in P. monodon.

Identification and characterization of the sex-determiner transformer-2 homologue in Chinese shrimp, Fenneropenaeus chinensis

The transformer-2 gene, encoding a protein (Tra-2) which directs sex-specifically alternative splicing of doublesex (dsx) pre-mRNA in combination with the transformer (Tra) protein, has been proved to play important roles on sex differentiation and sex development in Drosophila melanogaster. In the present study, a tra-2 homologue (FcTra-2) was cloned and characterized in the Chinese shrimp, Fenneropenaeus chinensis. A FcTra-2 genomic DNA sequence with a length of 8,871 bp was obtained and verified to consist of 7 exons and 6 introns. Three alternatively spliced mRNA transcripts, designated as FcTra-2a, FcTra-2b and FcTra-2c, were isolated and characterized. Sequence analysis showed that FcTra-2 included a RNA recognition motif and a linker region, which shared high sequence identities with Tra-2 from other species and 2 arginine/serine rich regions. Further studies were performed on the isoform FcTra-2c, since it exhibited a significantly higher expression level in ovary than in other tissues. In early developmental stages of the shrimp, FcTra-2c was detected to suddenly increase its expression level at the mysis stage. In juvenile stage, FcTra-2c displayed a significantly higher expression level in female Chinese shrimp than in males. These data indicated that FcTra-2 might be involved in female sex determination in Chinese shrimp.

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.

SUMO-1 of mud crab (Scylla paramamosain) in gametogenesis

The small ubiquitin-related modifier-1 (SUMO-1) is a member of a family of ubiquitin-related proteins. SUMO pathway, which is involved in gene expression in eukaryotic posttranslational processing, plays important roles in gene expression, genomic stability and the occurrence of cells, development and other biological processes. Scylla paramamosain is one of the important economic breeding crabs in the southeast coast of China. To date, little is known about the distinct roles of SUMO in crustacean, especially in crabs. In the present study, we report the identification and characterization of mud crab, S. paramamosain SUMO-1 (SpSUMO-1) gene using an approach which combines expressed sequence tag (EST) and rapid amplification cDNA end (RACE). The full length cDNA of SpSUMO-1 gene (GenBank: HM581660) is of 732 bp, including a 282 bp open reading frame which encodes a protein of 93 amino acids. Tissue distribution analysis showed that SpSUMO-1 was expressed more abundantly in the ovary than in other tissues (P<0.01). And the expression profiles of SpSUMO-1 in the different gonad developing stages revealed that the highest expression of SpSUMO-1 occurred at proliferation stage, and then decreased gradually as the ovarian development progressed, while in the testis, the expression level of SpSUMO-1 was relatively stable at different stages of testis development. The distribution of SpSUMO-1 mRNA and its protein was observed in the crab gametogenesis by in situ hybridization and immunocytochemical method respectively. In oogenesis, SpSUMO-1 transcripts presented at the cytoplasm and nucleus of oocytes from proliferation stage to primary vitellogenesis stage, but only appeared in the nucleus of oocytes in secondary and tertiary vitellogenesis stages. Meanwhile, SpSUMO-1 protein was localized in the cytoplasm of oogonia and different developing oocytes. On the other hand, the SpSUMO-1 transcript was detected throughout the spermatogenesis, with the strong positive signals of SpSUMO-1 presented at the nuclei of primary and secondary spermatocytes, spermatids and spermatozoa. Interestingly, the positive signals of acrosomal tubules of spermatozoa were also detected. SpSUMO-1 protein was localized in spermatogonium, primary spermatocyte, secondary spermatocyte and spermatid, but the positive signal was only detected in the nucleus of spermatozoa. All these results suggested that SUMO-1 may play essential roles in the gametogenesis of the crustacea.

Transcriptomics of a giant freshwater prawn (Macrobrachium rosenbergii): de novo assembly, annotation and marker discovery

BACKGROUND

Giant freshwater prawn (Macrobrachium rosenbergii or GFP), is the most economically important freshwater crustacean species. However, as little is known about its genome, 454 pyrosequencing of cDNA was undertaken to characterise its transcriptome and identify genes important for growth.

METHODOLOGY AND PRINCIPAL FINDINGS

A collection of 787,731 sequence reads (244.37 Mb) obtained from 454 pyrosequencing analysis of cDNA prepared from muscle, ovary and testis tissues taken from 18 adult prawns was assembled into 123,534 expressed sequence tags (ESTs). Of these, 46% of the 8,411 contigs and 19% of 115,123 singletons possessed high similarity to sequences in the GenBank non-redundant database, with most significant (E value < 1e(-5)) contig (80%) and singleton (84%) matches occurring with crustacean and insect sequences. KEGG analysis of the contig open reading frames identified putative members of several biological pathways potentially important for growth. The top InterProScan domains detected included RNA recognition motifs, serine/threonine-protein kinase-like domains, actin-like families, and zinc finger domains. Transcripts derived from genes such as actin, myosin heavy and light chain, tropomyosin and troponin with fundamental roles in muscle development and construction were abundant. Amongst the contigs, 834 single nucleotide polymorphisms, 1198 indels and 658 simple sequence repeats motifs were also identified.

CONCLUSIONS

The M. rosenbergii transcriptome data reported here should provide an invaluable resource for improving our understanding of this species' genome structure and biology. The data will also instruct future functional studies to manipulate or select for genes influencing growth that should find practical applications in aquaculture breeding programs.

A review of the major penaeid shrimp EST studies and the construction of a shrimp transcriptome database based on the ESTs from four penaeid shrimp

By economic value, shrimp is currently the most important seafood commodity worldwide, and these animals are often the subject of scientific research in shrimp farming countries. High throughput methods, such as expressed sequence tags (ESTs), were originally developed to study human genomics, but they are now available for studying other important organisms, including shrimp. ESTs are short sequences generated by sequencing randomly selected cDNA clones from a cDNA library. This is currently the most efficient and powerful method for providing transcriptomic data for organisms with an uncharacterized genome. This review will summarize the sixteen major shrimp EST studies that have been conducted to date. In addition, we analyzed the EST data downloaded from NCBI dbEST for the four major penaeid shrimp species and constructed a database to host all of these EST data as well as our own analysis results. This database provides the shrimp aquaculture research community with an outline of the shrimp transcriptome as well as a tool for shrimp gene identification.

A cDNA microarray, UniShrimpChip, for identification of genes relevant to testicular development in the black tiger shrimp (Penaeus monodon)

Background

Poor reproductive maturation in captive male broodstock of the black tiger shrimp (Penaeus monodon) is one of the serious problems to the farming industries. Without genome sequence, EST libraries of P. monodon were previously constructed to identify transcripts with important biological functions. In this study, a new version of cDNA microarray, UniShrimpChip, was constructed from the Peneaus monodon EST libraries of 12 tissues, containing 5,568 non-redundant cDNA clones from 10,536 unique cDNA in the P. monodon EST database. UniShrimpChip was used to study testicular development by comparing gene expression levels of wild brooders from the West and East coasts of Thailand and domesticated brooders with different ages (10-, 14-, 18-month-old).

Results

The overall gene expression patterns from the microarray experiments revealed distinct transcriptomic patterns between the wild and domesticated groups. Moreover, differentially expressed genes from the microarray comparisons were identified, and the expression patterns of eight selected transcripts were subsequently confirmed by reverse-transcriptase quantitative PCR (RT-qPCR). Among these, expression levels of six subunits (CSN2, 4, 5, 6, 7a, and 8) of the COP9 signalosome (CSN) gene family in wild and different ages of domesticated brooders were examined by RT-qPCR. Among the six subunits, CSN5 and CSN6 were most highly expressed in wild brooders and least expressed in the 18-month-old domesticated group; therefore, their full-length cDNA sequences were characterized.

Conclusions

This study is the first report to employ cDNA microarray to study testicular development in the black tiger shrimp. We show that there are obvious differences between the wild and domesticated shrimp at the transcriptomic level. Furthermore, our study is the first to investigate the feasibility that the CSN gene family might have involved in reproduction and development of this economically important species.

Isolation and characterization of genes functionally involved in ovarian development of the giant tiger shrimp Penaeus monodon by suppression subtractive hybridization (SSH)

Suppression subtractive hybridization (SSH) libraries between cDNA in stages I (previtellogenic) and III (cortical rod) ovaries of the giant tiger shrimp (Penaeus monodon) were established. In all, 452 ESTs were unidirectionally sequenced. Sequence assembly generated 28 contigs and 201 singletons, 109 of which (48.0%) corresponding to known sequences previously deposited in GenBank. Several reproduction-related transcripts were identified. The full-length cDNA of anaphase promoting complex subunit 11 (PmAPC11; 600 bp with an ORF of 255 bp corresponding to a polypeptide of 84 amino acids) and selenoprotein Mprecursor (PmSePM; 904 bp with an ORF of 396 bp corresponding to a polypeptide of 131 amino acids) were characterized and reported for the first time in penaeid shrimp. Semiquantitative RT-PCR revealed that the expression levels of PmSePM and keratinocyte-associated protein 2 significantly diminished throughout ovarian development, whereas Ser/Thrcheckpoint kinase 1 (Chk1), DNA replication licensing factor mcm2 and egalitarian were down-regulated in mature ovaries of wild P. monodon (p < 0.05). Accordingly, the expression profiles of PmSePM and keratinocyte-associated protein 2 could be used as biomarkers for evaluating the degree of reproductive maturation in domesticated P. monodon.

Expression of sex and reproduction-related genes in Marsupenaeus japonicus

Expressed sequence tags (ESTs) were identified from reciprocal suppression subtractive hybridization cDNA libraries from Marsupenaeus japonicus (Kuruma shrimp) female and male gonads. The expression profiles of 24 of these ESTs were determined in female and male gonads and developing postlarvae by real-time quantitative reverse transcription-polymerase chain reaction. When expression was determined in gonads, six of the ESTs were expressed in ovaries only, and five of the ESTs were expressed in testes only. When expression was determined in whole individuals during postlarval development, expression of the ESTs was low and inconsistent until stage PL110 (110 days since metamorphosis from mysis stage to the first postlarval stage). At PL110, seven of the ESTs were detected in females only, and seven ESTs were detected in males only. Sex-specific expression at this developmental stage indicates that these ESTs act as important gonadal development markers and may have a role in gametogenesis.