Comparative transcriptomic analysis of the different developmental stages of ovary in the cuttlefish Sepia pharaonis

The cuttlefish, Sepia pharaonis, is characterized by rapid growth and strong disease resistance, making it an important commercially farmed cephalopod species in the southeastern coastal regions of China. However, in the reproductive process of S. pharaonis, there are challenges such as a low output of eggs, poor quality, and low survival rates of newly hatched juveniles. Therefore, there is an urgent need to study the molecular mechanisms underlying ovarian development in this species. In this study, we conducted the first transcriptomic analysis of the ovary at four developmental stages: the undeveloped stage, developing stage, nearly-ripe stage, and ripe stage, and compared the transcriptomics among these four stages using Illumina sequencing technology. The total numbers of clean reads of the four stages ranged from 40,890,772 to 52,055,714 reads. A total of 136,829 DEGs were obtained, GC base ratios of raw data were between 38.44 and 44.59%, and the number of uniquely mapped reads spanned from 88.08 to 95.90%. The Pearson correlation coefficient demonstrated a strong correlation among different samples within the same group, PCA and Anosim analysis also revealed that the grouping of these four stages was feasible, and each stage could be distinguished from the others. GO enrichment analysis demonstrated that ovarian follicle growth, sex differentiation, and transforming growth factor beta receptor, played a foreshadowing role at the early ovarian development stage, and the terms of small molecule metabolic process, peptide metabolic process, and catalytic activity were prominent at the mature stage. Meanwhile, KEGG analysis showed that the early ovarian development of S. pharaonis was mainly associated with the cell cycle, DNA replication, and carbon metabolism, while the mid-late ovarian development was involved with the signal transduction, endocrine system, and reproduction pathway. RT-qPCR further confirmed the consistent expression patterns of genes such as 17β-HSD, GH, VGS, NFR, and NYR in the ovaries of S. pharaonis, exhibiting elevated levels of expression during the maturation stage. Conversely, ER and OM exhibited high expression levels during the early stages of ovarian development. These transcriptomic data provide insights into the molecular mechanisms of S. pharaonis ovarian development. The findings of this study will contribute to improving the reproduction and development of cuttlefish and enriching the bioinformatics knowledge of cephalopods.