Constructing and random sequencing analysis of normalized cDNA library of testis tissue from oriental river prawn (Macrobrachium nipponense)

Comparative transcriptome analysis of brain and gonad reveals reproduction-related miRNAs in the giant prawn, Macrobrachium rosenbergii

Macrobrachium rosenbergii (M. rosenbergii), as a species of common prawn, is a delicacy that is consumed all over the world. By interacting with the target gene 3′-untranslated region (3'-UTR), microRNAs (miRNAs) regulate its expression and ultimately participate in the regulation of reproductive development. However, research focusing on miRNA regulation during gonadal development in M. rosenbergii received very little attention. To explore the association between miRNA and reproduction, we performed RNA sequencing (RNA-seq) on brain and gonad organs in male and female M. rosenbergii. A total of 494 miRNAs were obtained in RNA-seq, including 31 and 59 differentially expressed (DE) miRNAs in the brain and gonads, respectively. Furthermore, 9 DE miRNAs were randomly selected from the brain and gonads, and qRT-PCR was conducted to validate the results of RNA-seq. Interestingly, dpu-miR-133 was found to be substantially expressed in the male brain and testis but poorly expressed in the female brain, ovary, and other organs. Analysis of dpu-miR-133 by Targetscan and MiRanda predicted to target 5-HT1. Furthermore, the dual-luciferase reporter assay manifested that dpu-miR-133 can combine with 5-HT1. Overall, our research work provides basic data for further study on the miRNA-mediated regulation of brain, gonad, and reproductive development of study M. rosenbergii.

Cathepsin D Plays a Vital Role in Macrobrachium nipponense of Ovary Maturation: Identification, Characterization, and Function Analysis

The oriental river prawn Macrobrachium nipponense is an economically important aquacultural species. However, its aquaculture is negatively impacted by the rapid sexual maturation of female M. nipponense. The fast sexual maturation produces a large number of offspring which leads to a reduction in resilience, a low survival rate, and an increased risk of hypoxia, this in turn, seriously affects the economic benefits of prawn farming. Cathepsin D is a lysosomal protease involved in the ovarian maturation of M. nipponense. In the current study, the cDNA of the gene encoding cathepsin D (Mn-CTSD) was cloned from M. nipponense. The total length was 2391 bp and consisted of an open reading frame (ORF) of 1158 bp encoding 385 amino acids. Sequence analysis confirmed the presence of conserved N-glycosylation sites and characteristic sequences of nondigestive cathepsin D. The qPCR analysis indicated that Mn-CTSD was highly expressed in all tissues tested, most significantly in the ovaries, whereas in situ hybridization showed that expression occurred mainly in oocyte nuclei. Analysis of its expression during development showed that Mn-CTSD peaked during the O-IV stage of ovarian maturation. For the RNAi interference experiment, female M. nipponense specimens in the ovary stage I were selected. Injection of Mn-CTSD double-stranded (ds)RNA into female M. nipponense decreased the expression of Mn-CTSD in the ovaries, such that the Gonad Somatic Index (GSI) of the experimental group was significantly lower than that of the control group (1.79% versus 4.57%; p < 0.05). Ovary development reached the O-III stage in 80% of the control group, compared with 0% in the experimental group. These results suggest that Mn-CTSD dsRNA inhibits ovarian maturation in M. nipponense, highlighting its important role in ovarian maturation in this species and suggesting an approach to controlling ovarian maturation during M. nipponense aquaculture.

Transcriptome Profiling Analysis of the Testis After Eyestalk Ablation for Selection of the Candidate Genes Involved in the Male Sexual Development in Macrobrachium nipponense

The eyestalk of crustacean species secretes many hormones, affecting the process of reproduction, molting, metabolism of glucose, and other functions in crustaceans. In this study, important metabolic pathways and candidate genes involved in the male sexual development were identified through performing the transcriptome profiling analysis of the testis after the ablation of eyestalk from Macrobrachium nipponense. The histological observations revealed that the testis development became vigorous after eyestalk ablation, indicating that the hormones secreted by the eyestalk have negative effects on the testis development in M. nipponense. Transcriptome profiling analysis revealed that 1,039, 1,226, and 3,682 differentially expressed genes (DEGs) were identified between normal prawns (CG) vs single-side eyestalk ablation prawns (SS), SS vs double-side eyestalk ablation prawns (DS), and CG vs DS, respectively, indicating that the ablation of double-side eyestalk has more significant regulatory roles on male sexual development than that of single-side ablation, which was consistent with the histological observations. Lysosome, Apoptosis, Glycolysis/Gluconeogenesis, and Insulin signaling pathway were the main enriched metabolic pathways in all of these three comparisons, and the important genes from these metabolic pathways were also selected. The qPCR verifications of 10 DEGs from these metabolic pathways were the same as those of RNA-seq. The qPCR, in situ hybridization, and RNA interference analysis of Mn-NF_kBα revealed that NF_kBα has a positive regulatory effect on testis development. This study provided new insights on male sexual development in M. nipponense, promoting the studies on male sexual development in other crustaceans as well.

De novo gonad transcriptome analysis of the common littoral shrimp Palaemon serratus: novel insights into sex-related genes

BACKGROUND

The common littoral shrimp Palaemon serratus is an economically important decapod resource in some European communities. Aquaculture practices prevent the genetic deterioration of wild stocks caused by overfishing and at the same time enhance the production. The biotechnological manipulation of sex-related genes has the proved potential to improve the aquaculture production but the scarcity of genomic data about P. serratus hinders these applications. RNA-Seq analysis has been performed on ovary and testis samples to generate a reference gonadal transcriptome. Differential expression analyses were conducted between three ovary and three testis samples sequenced by Illumina HiSeq 4000 PE100 to reveal sex-related genes with sex-biased or sex-specific expression patterns.

RESULTS

A total of 224.5 and 281.1 million paired-end reads were produced from ovary and testis samples, respectively. De novo assembly of ovary and testis trimmed reads yielded a transcriptome with 39,186 transcripts. The 29.57% of the transcriptome retrieved at least one annotation and 11,087 differentially expressed genes (DEGs) were detected between ovary and testis replicates. Six thousand two hundred seven genes were up-regulated in ovaries meanwhile 4880 genes were up-regulated in testes. Candidate genes to be involved in sexual development and gonadal development processes were retrieved from the transcriptome. These sex-related genes were discussed taking into account whether they were up-regulated in ovary, up-regulated in testis or not differentially expressed between gonads and in the framework of previous findings in other crustacean species.

CONCLUSIONS

This is the first transcriptome analysis of P. serratus gonads using RNA-Seq technology. Interesting findings about sex-related genes from an evolutionary perspective (such as Dmrt1) and for putative future aquaculture applications (Iag or vitellogenesis genes) are reported here. We provide a valuable dataset that will facilitate further research into the reproductive biology of this shrimp.

Integrated analysis of differentially expressed microRNAs and mRNAs to screen miRNAs and genes related to reproduction in Macrobrachium nipponense

Female Macrobrachium nipponense has the characteristic of short sexual maturity during the breeding season, which can increase breeding risk and lead to prevalent female individual miniaturization. In this study, we characterized micro (mi)RNA-seq data of the eyestalk (E) and cerebral ganglia (B) of female M. nipponense during breeding and non-breeding seasons. A total of 393 and 189 differentially expressed miRNAs (DE miRNAs) were identified in BSE vs. NBSE and BSB vs. NBSB, respectively. The most abundant up- and down-regulated DE miRNAs were miR-124, miR-14, and miR-7. Enrichment analysis showed that DE miRNA target genes were mainly involved in 'metabolic process' and 'binding', and were associated with 'neurohormonal regulation' and 'photoreceptor activity' signaling pathways. Integrated analysis of miRNA-mRNA expression showed that the most abundant DE miRNAs were miR-14 and miR-278 in BSE vs. NBSE and BSB vs. NBSB, respectively. Four pairs of DE miRNAs and their corresponding target annotated genes were selected from the DE miRNA-mRNA interaction network (bmo-miR-316-5p/opsin protein, ame-miR-125/skeletal muscle actin 8, dmo-miR-278/sugar transporter, and tca-miR-3885-5p/5-HT1 receptor). Gene expression analysis of these four pairs in different ovary development stages showed their potential regulatory roles in ovary maturation.

Validation and Evaluation of Reference Genes for Quantitative Real-Time PCR in Macrobrachium Nipponense

Quantitative real-time PCR (qPCR) is widely used in molecular biology, although the accuracy of the quantitative results is determined by the stability of the reference genes used. Recent studies have investigated suitable reference genes for some crustaceans under various conditions, but studies in Macrobrachium nipponense are currently lacking. In this study, we selected the following seven genes from among 35 commonly used housekeeping genes as candidate qPCR reference genes for temporal and spatial expression: EIF (eukaryotic translation initiation factor 5A), 18S (18S ribosomal RNA), EF-1α (elongation factor-1α), GAPDH (glyceraldehyde-3-phosphate dehydrogenase), TUB (α-tubulin), β-act (β-actin), and RPL18 (Ribosomal protein L18). The stability of each reference gene was evaluated by GeNorm, NormFinder, BestKeeper, and comparative ∆C _t methods, and was comprehensively ranked using RefFinder. RPL18 was shown to be the most suitable reference gene for adult M. nipponense tissues, while EIF was the most stable in different ovarian and embryo stages and in white spot syndrome virus infection, and β-act was the most stable reference gene under hypoxia stress. The reliability of the rankings was confirmed by RNA interference experiments. To the best of our knowledge, this represents the first systematic analysis of reference genes for qPCR experiments in M. nipponense, and the results will provide invaluable information for future research in closely related crustaceans.

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.

Molecular insights into reproduction regulation of female Oriental River prawns Macrobrachium nipponense through comparative transcriptomic analysis

The oriental river prawn, Macrobrachium nipponense, is an important commercial aquaculture resource in China. During breeding season, short ovary maturation cycles of female prawns cause multi-generation reunions in ponds and affect the growth of females representing individual miniaturization (known as autumn -propagation). These reproductive characteristics pose problems for in large - scale farming. To date, the molecular mechanisms of reproduction regulation of M. nipponense remain unclear. To address this issue, we performed transcriptome sequencing and gene expression analyses of eyestalk and cerebral ganglia of female M. nipponense during breeding and non-breeding seasons. Differentially expressed gene enrichment analysis results revealed several important reproduction related terms and signaling pathways, such as "photoreceptor activity", "structural constituent of cuticle" and "G-protein coupled receptor activity". The following six key genes from the transcriptome were predicted to mediate environmental factors regulating reproduction of M. nipponense: neuroparsin, neuropeptide F II, orcokinin II, crustacean cardioactive peptide, pigment-dispersing hormone 3 and tachykinin. These results will contribute to a better understanding of the molecular mechanisms of reproduction of oriental river prawns. Further detailed functional analyses of the candidate reproduction regulation related neuropeptides are needed to shed light on the mechanisms of reproduction of crustacean.

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.

Insights into Sexual Precocity of Female Oriental River Prawn Macrobrachium nipponense through Transcriptome Analysis

Background

The oriental river prawn (Macrobrachium nipponense) is the most prevalent aquaculture species in China. The sexual precocity in this species has received considerable attention in recent years because more and more individuals matured at a small size, which devalues the commercial production. In this study, we developed deep-coverage transcriptomic sequencing data for the ovaries of sexually precocious and normal sexually mature M. nipponense using next-generation RNA sequencing technology and attempted to provide the first insight into the molecular regulatory mechanism of sexual precocity in this species.

Results

A total of 63,336 unigenes were produced from the ovarian cDNA libraries of sexually precocious and normal sexually mature M. nipponense using Illumina HiSeq 2500 platform. Through BLASTX searches against the NR, STRING, Pfam, Swissprot and KEGG databases, 15,134 unigenes were annotated, accounting for 23.89% of the total unigenes. 5,195 and 3,227 matched unigenes were categorized by GO and COG analysis respectively. 15,908 unigenes were consequently mapped into 332 KEGG pathways, and many reproduction-related pathways and genes were identified. Moreover, 26,008 SSRs were identified from 18,133 unigenes. 80,529 and 80,516 SNPs were yielded from ovarian libraries of sexually precocious and normal sexually mature prawn, respectively, and 29,851 potential SNPs between these two groups were also predicted. After comparing the ovarian libraries of sexually precocious and normal sexually mature prawn, 549 differentially expressed genes (DEGs) and 9 key DEGs that may be related to sexual precocity of M. nipponense were identified. 20 DEGs were selected for validation by quantitative real-time PCR (QPCR) and 19 DEGs show consistent expression between QPCR and RNAseq-based differential expression analysis datasets.

Conclusion

This is the first report on the large-scale RNA sequencing of ovaries of sexually precocious and normal sexually mature M. nipponense. The annotated transcriptome data will provide fundamental support for future research into the reproduction biology of M. nipponense. The large number of candidate SNPs and SSRs detected in this study could be used as genetic markers for population genetics and functional genomics in this species. More importantly, many DEGs, especially nine key DEGs between sexually precocious and normal sexually mature prawns were identified, which will dramatically improve understanding of molecular regulatory mechanism of sexual precocity of this species.

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.

Characterization of MMP-9 gene from a normalized cDNA library of kidney tissue of yellow catfish (Pelteobagrus fulvidraco)

Matrix metalloproteinase-9 (MMP-9), one of members of the MMP family, is important for the cleaving of structural extracellular matrix (ECM) molecules and involved in inflammatory processes. In this study, MMP-9 cDNA was isolated and characterized from a normalized cDNA library of kidney tissue of yellow catfish (designated as YcMMP-9). The complete sequence of YcMMP-9 cDNA consisted of 2561 nucleotides. The open reading frame potentially encoded a protein of 685 amino acids with a calculated molecular mass of approximately 77.182 kDa. Amino acid sequence of YcMMP-9 have typical characteristics of MMP-9 family and showed highest identity (85.3%) to channel catfish MMP-9. The YcMMP-9 genomic DNA contains 13 exons and 12 introns. Quantitative RT-PCR (qRT-PCR) analysis showed that YcMMP-9 mRNA was constitutively expressed in all examined tissues in normal fish with high expression in head kidney, trunk kidney, blood, and spleen. However, expression of YcMMP-9 mRNA was induced by Aeromonas hydrophila stimulation, especially in these four tissues mentioned above. It indicated that YcMMP-9 was involved in innate immune responses against bacterial infection.

Transciptomic and histological analysis of hepatopancreas, muscle and gill tissues of oriental river prawn (Macrobrachium nipponense) in response to chronic hypoxia

Background

Oriental river prawn, Macrobrachium nipponense, is a commercially important species found in brackish and fresh waters throughout China. Chronic hypoxia is a major physiological challenge for prawns in culture, and the hepatopancreas, muscle and gill tissues play important roles in adaptive processes. However, the effects of dissolved oxygen availability on gene expression and physiological functions of those tissues of prawns are unknown. Adaptation to hypoxia is a complex process, to help us understand stress-sensing mechanism and ultimately permit selection for hypoxia- tolerant prawns, we performed transcriptomic analysis of juvenile M. nipponense hepatopancreas, gill and muscle tissues by RNA-Seq.

Results

Approximately 46,472,741; 52,773,612 and 58,195,908 raw sequence reads were generated from hepatopancreas, muscle and gill tissues, respectively. A total of 62,722 unigenes were generated, of the assembled unigenes, we identified 8,892 genes that were significantly up-regulated, while 5,760 genes were significantly down-regulated in response to chronic hypoxia. Genes from well known functional categories and signaling pathways associated with stress responses and adaptation to extreme environments were significantly enriched, including genes in the functional categories “response to stimulus”, “transferase activity” and “oxidoreductase activity”, and the signaling pathways “oxidative phosphorylation”, “glycolysis/gluconeogenesis” and “MAPK signaling”. The expression patterns of 18 DEGs involved in hypoxic regulation of M. nipponense were validated by quantitative real-time reverse-transcriptase polymerase chain reactions (qRT-PCR; average correlation coefficient = 0.94). In addition, the hepatopancreas and gills exhibited histological differences between hypoxia and normoxia groups. These structural alterations could affect the vital physiological functions of prawns in response to chronic hypoxia, which could adversely affect growth and survival of M. nipponense.

Conclusions

Gene expression changes in tissues from the oriental river prawn provide a preliminary basis to better understand the molecular responses of M. nipponense to chronic hypoxia. The differentially expressed genes (DEGs) identified in M. nipponense under hypoxia stress may be important for future genetic improvement of cultivated prawns or other crustaceans through transgenic approaches aimed at increasing hypoxia tolerance.

Electronic supplementary material

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

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.

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.