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 Sex-Specific Splicing of Doublesex in Brine Shrimp Artemia franciscana

The understanding of sex determination and differentiation in animals has recently made remarkable strides through the use of advanced research tools. At the gene level, the Mab-3-related transcription factor (Dmrt) gene family, which encodes for the typical DNA-binding doublesex/Mab-3 (DM) domain in their protein, is known for its contribution to sex determination and differentiation in insects. In this study, DNA-binding DM domain screening has identified eight transcripts from Artemia franciscana transcriptomic that encode proteins containing one conserved DNA-binding DM domain. The genome mapping confirmed that these eight transcripts are transcribed from six different loci on the A. franciscana genome assembly. One of those loci, the Af.dsx-4 locus, is closely related to Doublesex, a gene belonging to the Dmrt gene family. This locus could be transcribed into three alternative transcripts, namely Af.dsx⁴, Af.dsx^F and Af.dsx^M. While Af.dsx⁴ and Af.dsx^F could putatively be translated to form an identical Af.dsx^F protein of 186 aa long, Af.dsx^M translates for an Af.dsx^M protein of 289 aa long but shares a DNA-binding DM domain. Interestingly, Af.dsx^F and Af.dsx^M are confirmed as sex-specific transcripts, Af.dsx^F is only present in females, and Af.dsx^M is only present in male individuals. The results suggest that the sex-specific splicing mechanism of the doublesex described in insects is also present in A. franciscana. Af.dxs-4 locus can be used in further studies to clarify the sex determination pathways in A. fracnciscana.

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.

Comparative Transcriptome Analysis of Gonads for the Identification of Sex-Related Genes in Giant Freshwater Prawns (Macrobrachium Rosenbergii) Using RNA Sequencing

The giant freshwater prawn (Macrobrachium rosenbergii) exhibits sex dimorphism between the male and female individuals. To date, the molecular mechanism governing gonadal development was unclear, and limited data were available on the gonad transcriptome of M. rosenbergii. Here, we conducted comprehensive gonadal transcriptomic analysis of female (ZW), super female (WW), and male (ZZ) M. rosenbergii for gene discovery. A total of 70.33 gigabases (Gb) of sequences were generated. There were 115,338 unigenes assembled with a mean size of 1196 base pair (bp) and N50 of 2195 bp. Alignment against the National Center for Biotechnology Information (NCBI) non-redundant nucleotide/protein sequence database (NR and NT), the Kyoto Encyclopedia of Genes and Genomes (KEGG) database, SwissProt database, Protein family (Pfam), Gene ontology (GO), and the eukaryotic orthologous group (KOG) database, 36,282 unigenes were annotated at least in one database. Comparative transcriptome analysis observed that 10,641, 16,903, and 3393 genes were significantly differentially expressed in ZW vs. ZZ, WW vs. ZZ, and WW vs. ZW samples, respectively. Enrichment analysis of differentially expressed genes (DEGs) resulted in 268, 153, and 42 significantly enriched GO terms, respectively, and a total of 56 significantly enriched KEGG pathways. Additionally, 23 putative sex-related genes, including Gtsf1, IR, HSP21, MRPINK, Mrr, and other potentially promising candidate genes were identified. Moreover, 56,241 simple sequence repeats (SSRs) were identified. Our findings provide a valuable archive for further functional analyses of sex-related genes and future discoveries of underlying molecular mechanisms of gonadal development and sex determination.

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.

Gill transcriptome alterations in Macrobrachium rosenbergii under copper exposure

Copper is one of common contaminants in estuaries and coastal zones, which may cause physiological dysfunction in aquatic organisms. However, molecular response triggered by Cu have remained largely unknown in freshwater prawn Macrobrachium rosenbergii. In the present study, we performed transcriptomic analysis to characterize molecular mechanisms of copper immunotoxicity in gills from M. rosenbergii. A large number of potential simple sequence repeats (SSRs) and single nucleotide polymorphisms (SNPs) loci in the transcriptome were identified. 19,417 and 8989 differentially expressed genes (DEGs) were obtained at 3 h and 48 h after exposure, respectively. Most of these DEGs were down-regulated implying that gene expressions were largely inhibited by Cu, which might lead to impairments of biological functions. Functional enrichment analysis of these DEGs revealed that immune, detoxification and apoptosis were the differentially regulated processes by Cu stress. 12 DGEs involved in immune response and heavy metal detoxification were discovered and validated by qRT-PCR. The results indicated that the M. rosenbergii might counteract the toxicity of Cu at the transcriptomic level by increasing expressions of immune- and heavy metal detoxification-related genes, and these selected genes could be used as molecular indicators for Cu stress. Our study firstly reported the stress response at transcriptional level in M. rosenbergii during Cu exposure. The genes and pathways identified here not only give us new insight into molecular mechanisms underlying Cu toxicity effects in prawn, but facilitate biomarker identification and stress-resistant breeding studies.

Screening for Dmrt genes from embryo to mature Macrobrachium rosenbergii prawns

The doublesex and mab-3 related transcription factor (Dmrt) gene family is known to be related to the sexual regulators doublesex of arthropods and mab-3 of annelids and to hold highly conserved functions in sexual determination and differentiation across phyla. Here, we report a study of the Dmrt gene family in the freshwater prawn Macrobrachium rosenbergii, a crustacean whose sexual differentiation has been widely researched. A wide transcriptomic screen, from the embryo to the adult M. rosenbergii, identified five novel Dmrt genes (MroDmrts) and confirmed two known MroDmrts. The seven MroDmrts encode proteins of 275-855 amino acids; each protein contained at least one conserved DNA-binding DM domain, which is typical of Dmrt proteins, and five proteins contained 1-4 transactivation domains (TADs). Importantly, in the embryonic, larval and post-larval stages, MroDmrt genes exhibited time-dependent expression patterns rather than sex-specific expression. In-silico screening of the expression of the MroDmrt genes in adult males revealed the enrichment of MroiDmrt1b and MroiDmrt1c in the androgenic gland (AG) as compared to the eyestalks. In vivo silencing of the androgenic gland insulin-like (IAG) encoding gene significantly decreased the expression of the above two Dmrt genes, while not affecting the expression of control genes, thereby suggesting the possible role of these two genes in the IAG-switch and in sex-differentiation processes.

A high-throughput SNP discovery strategy for RNA-seq data

Background

Single nucleotide polymorphisms (SNP) have been applied as important molecular markers in genetics and breeding studies. The rapid advance of next generation sequencing (NGS) provides a high-throughput means of SNP discovery. However, SNP development is limited by the availability of reliable SNP discovery methods. Especially, the optimum assembler and SNP caller for accurate SNP prediction from next generation sequencing data are not known.

Results

Herein we performed SNP prediction based on RNA-seq data of peach and mandarin peel tissue under a comprehensive comparison of two paired-end read lengths (125 bp and 150 bp), five assemblers (Trinity, IDBA, oases, SOAPdenovo, Trans-abyss) and two SNP callers (GATK and GBS). The predicted SNPs were compared with the authentic SNPs identified via PCR amplification followed by gene cloning and sequencing procedures. A total of 40 and 240 authentic SNPs were presented in five anthocyanin biosynthesis related genes in peach and in nine carotenogenic genes in mandarin. Putative SNPs predicted from the same RNA-seq data with different strategies led to quite divergent results. The rate of false positive SNPs was significantly lower when the paired-end read length was 150 bp compared with 125 bp. Trinity was superior to the other four assemblers and GATK was substantially superior to GBS due to a low rate of missing authentic SNPs. The combination of assembler Trinity, SNP caller GATK, and the paired-end read length 150 bp had the best performance in SNP discovery with 100% accuracy both in peach and in mandarin cases. This strategy was applied to the characterization of SNPs in peach and mandarin transcriptomes.

Conclusions

Through comparison of authentic SNPs obtained by PCR cloning strategy and putative SNPs predicted from different combinations of five assemblers, two SNP callers, and two paired-end read lengths, we provided a reliable and efficient strategy, Trinity-GATK with 150 bp paired-end read length, for SNP discovery from RNA-seq data. This strategy discovered SNP at 100% accuracy in peach and mandarin cases and might be applicable to a wide range of plants and other organisms.

Electronic supplementary material

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

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.

Enrichment of genomic resources and identification of simple sequence repeats from medicinally important Clausena excavata

To broaden and delve into the genomic information of Clausena excavata, an important medicinal plant in many Asian countries, RNA sequencing (RNA-seq) analysis was performed and a total of 16,638 non-redundant unigenes (≥ 300 bp) with an average length of 755 bp were generated by de novo assembly from 17,580,456 trimmed clear reads. The functional categorization of the identified unigenes by a gene ontology (GO) term resulted in 2305 genes in the cellular component, 5577 in the biological processes, and 8056 in the molecular functions, respectively. The top sub-category in biological processes was the metabolic process with 4374 genes. Among annotated genes, 3006 were mapped to 123 metabolic pathways by the Kyoto Encyclopedia of Genes and Genomes (KEGG) metabolic pathway analysis tool. The search for simple sequence repeats (SSRs) resulted in 845 SSRs from 749 SSR-containing unigenes and the most abundant SSR motifs was AAG/CTT with 179 occurrences. Twelve SSR markers were tested for cross transferability among five Clausena species; eight of them exhibited polymorphism. Taken together, these data provide valuable resources for genomic or genetic studies of Clausena species and other relative studies. The transcriptome shotgun assembly data have been deposited at DDBJ/EMBL/GenBank under the accession GGEM00000000.

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.

Transcriptome profiling of the Macrobrachium rosenbergii lymphoid organ under the white spot syndrome virus challenge

Macrobrachium rosenbergii is a crustacean with economic importance, and adult prawns are generally thought to be tolerant to white spot syndrome virus (WSSV) infection. Although certain genes are known to respond to WSSV infection and lymphoid tissue is an important immune organ, the response of lymphoid organ to WSSV infection is unclear. Next-generation sequencing was employed in this study to determine the transcriptome differences between WSSV infection and mock lymphoid organs. A total of 44,606,694 and 40,384,856 clean reads were generated and assembled into 73,658 and 72,374 unigenes from the control sample and the WSSV infection sample, respectively. Based on homology searches, KEGG, GO, and COG analysis, 21,323 unigenes were annotated. Among them, 4951 differential expression genes were identified and categorized into 244 metabolic pathways. Coagulation cascades, and pattern recognition receptor signaling pathways were used as examples to discuss the response of host to WSSV infection. We also identified 12,308 simple sequence repeats, which can be further used as functional markers. Results contribute to a better understanding of the immune response of prawn lymphoid organ to WSSV and provide information for identifying novel genes in the absence of the prawn genome.