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.