Transcriptome Sequencing and De Novo Assembly of Golden Cuttlefish Sepia esculenta Hoyle

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.

Preliminary study on toxicological mechanism of golden cuttlefish (Sepia esculenta) larvae exposed to cd

BACKGROUND

Cadmium (Cd) flows into the ocean with industrial and agricultural pollution and significantly affects the growth and development of economic cephalopods such as Sepia esculenta, Amphioctopus fangsiao, and Loligo japonica. As of now, the reasons why Cd affects the growth and development of S. esculenta are not yet clear.

RESULTS

In this study, transcriptome and four oxidation and toxicity indicators are used to analyze the toxicological mechanism of Cd-exposed S. esculenta larvae. Indicator results indicate that Cd induces oxidative stress and metal toxicity. Functional enrichment analysis results suggest that larval ion transport, cell adhesion, and some digestion and absorption processes are inhibited, and the cell function is damaged. Comprehensive analysis of protein-protein interaction network and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis was used to explore S. esculenta larval toxicological mechanisms, and we find that among the 20 identified key genes, 14 genes are associated with neurotoxicity. Most of them are down-regulated and enriched to the neuroactive ligand-receptor interaction signaling pathway, suggesting that larval nervous system might be destroyed, and the growth, development, and movement process are significantly affected after Cd exposure.

CONCLUSIONS

S. esculenta larvae suffered severe oxidative damage after Cd exposure, which may inhibit digestion and absorption functions, and disrupt the stability of the nervous system. Our results lay a function for understanding larval toxicological mechanisms exposed to heavy metals, promoting the development of invertebrate environmental toxicology, and providing theoretical support for S. esculenta artificial culture.

Preliminary Exploration of Metabolic Mechanisms in Copper-Exposed Sepia esculenta Based on Transcriptome Analysis

As a common and high-concentration heavy metal in the ocean, Cu can induce metal toxicity and significantly affect the metabolic function of marine organisms. Sepia esculenta is an important economic cephalopod found along the east coast of China, the growth, movement, and reproduction of which are all affected by heavy metals. Hitherto, the specific metabolic mechanism of heavy-metal exposure in S. esculenta is still unclear. In this study, we identified 1131 DEGs through transcriptome analysis of larval S. esculenta within 24 h of Cu exposure. GO and KEGG functional enrichment analysis results indicated that Cu exposure may affect purine metabolism, protein digestion and absorption, cholesterol metabolism, and other metabolic processes in S. esculenta larvae. It is worth noting that in this study we explore metabolic mechanism of Cu-exposed S. esculenta larvae through the comprehensive analysis of protein–protein interaction network and KEGG enrichment analysis for the first time and find 20 identified key and hub genes such as CYP7A1, CYP3A11, and ABCA1. Based on their expression, we preliminarily speculate that Cu exposure may inhibit multiple metabolic processes and induce metabolic disorders. Our results lay a foundation for further understanding the metabolic mechanism of S. esculenta against heavy metals and provide theoretical help for S. esculenta artificial breeding.

Effect of acute Cu exposure on immune response mechanisms of golden cuttlefish (Sepia esculenta)

Sepia esculenta is a common economic cephalopod that has received extensive attention due to the tender meat, rich protein content and certain medicinal value thereof. Over the past decade, multiple industries have discharged waste into the ocean in large quantities, thereby significantly increasing the concentration of heavy metals in the ocean. Copper (Cu) is a common heavy metal in the ocean. The increase of Cu content will affect numerous biological processes such as immunity and metabolism of marine organisms. High concentrations of Cu may inhibit S. esculenta growth, development, swimming, and other processes, which would significantly affect its culture. In this research, transcriptome analysis is used to initially explore Cu-exposed S. esculenta larval immune response mechanisms. And compared to control group with normally growing larvae, 2056 differentially expressed genes (DEGs) are identified in experimental group with Cu-exposed larvae. The results of DEGs functional enrichment analyses including GO and KEGG indicate that Cu exposure might promote inflammatory and innate immune responses in cuttlefish larvae. Then, 10 key genes that might regulate larval immunity are identified using a comprehensive analysis that combines protein-protein interaction (PPI) network and KEGG functional enrichment analyses, of which three genes with the highest number of protein interactions or involve in more KEGG signaling pathways are identified as hub genes that might significantly affect larval immune response processes. Comprehensive analysis of PPI network and KEGG signaling pathway are used for the first time to explore Cu-exposed S. esculenta larval immune response mechanisms. Our results preliminarily reveal immune response mechanisms of cephalopods exposed to heavy metals and provide valuable resources for further understanding mollusk immunity.

Exploration of immune response mechanisms in cadmium and copper co-exposed juvenile golden cuttlefish (Sepia esculenta) based on transcriptome profiling

Sepia esculenta is a popular economic cephalopod with high yield, delicious meat, and rich nutrition. With the rapid development of heavy industry and medical industry, a large amount of waste has been released into the ocean recklessly in recent years, inducing a significant increase in the content of heavy metals, especially cadmium (Cd) and copper (Cu), in the ocean. This phenomenon significantly affects the growth and development of S. esculenta, causing a serious blow to its artificial breeding. In this study, transcriptome analysis is used to initially explore immune response mechanisms of Cd and Cu co-exposed juvenile S. esculenta. The results show that 1,088 differentially expressed genes (DEGs) are identified. And DEGs functional enrichment analysis results suggests that co-exposure may promote inflammatory and innate immune responses in juvenile S. esculenta. Fifteen key genes that might regulate the immunity of S. esculenta are identified using protein-protein interaction (PPI) network and KEGG enrichment analyses, of which the three genes with the highest number of interactions or involve in more KEGG pathways are identified as hub genes that might significantly affect the immune response processes. Comprehensive analysis of PPI network and KEGG signaling pathway is used for the first time to explore co-exposed S. esculenta juvenile immune response processes. Our results preliminarily reveal immune response mechanisms of cephalopods exposed to heavy metals and provide a valuable resource for further understanding of mollusk immunity.

De novo transcriptome sequencing and analysis of the cuttlefish (Sepiella japonica) with different embryonic developmental stages

In this study, embryos of Sepiella japonica from eye primordium formation to the larval growing stage were collected and used for RNA-Seq analysis. A total of 183,542,186 clean reads were assembled de novo into 58,054 unigenes consisting of 54,118,228 bp, with the average length at 932 bp and the N50 at 1667 bp. 21,469 (36.98%) unigenes were annotated at least in one of four databases including non-redundant protein (NR), Swiss-Prot, clusters of orthologous groups of proteins (KOG) and Kyoto Encyclopedia of Genes and Genomes (KEGG). 4460 (7.68%) unigenes were annotated in all databases. Analysis of differentially expressed genes (DEGs) was carried out on embryos at Eye primordium formation stage (SJ1), organ differentiation stage (SJ2), and hatching stage (SJ3). Overall, the current study provided the de novo assembly of S. japonica transcriptome and identified the DEGs and pathways during embryonic development, which will provide a fundamental genetic resource for further functional research.

Global impact of diet and temperature over aquaculture of Octopus vulgaris paralarvae from a transcriptomic approach

Common octopus, Octopus vulgaris, is an economically important cephalopod species. However, its rearing under captivity is currently challenged by massive mortalities previous to their juvenile stage due to nutritional and environmental factors. Dissecting the genetic basis and regulatory mechanism behind this mortality requires genomic background knowledge. A transcriptomic sequencing of 10 dph octopus paralarvae from different experimental conditions was constructed via RNA-seq. A total of 613,767,530 raw reads were filtered and de novo assembled into 363,527 contigs of which 82,513 were annotated in UniProt carrying also their GO and KEGG information. Differential gene expression analysis was carried out on paralarvae reared under different diet regimes and temperatures, also including wild paralarvae. Genes related to lipid metabolism exhibited higher transcriptional levels in individuals whose diet includes crustacean zoeas, which had an impact over their development and immune response capability. High temperature induces acclimation processes at the time that increase metabolic demands and oxidative stress. Wild individuals show an expression profile unexpectedly similar to Artemia fed individuals. Proteomic results support the hypothesis revealed by transcriptional analysis. The comparative study of the O. vulgaris transcriptomic profiles allowed the identification of genes that deserve to be further studied as candidates for biomarkers of development and health. The results obtained here on the transcriptional variations of genes caused by diet and temperature will provide new perspectives in understanding the molecular mechanisms behind nutritional and temperature requirements of common octopus that will open new opportunities to deepen in paralarvae rearing requirements.

Transcriptome analysis for identifying possible causes of post-reproductive death of Sepia esculenta based on brain tissue

BACKGROUND

The subpeduncle lobe/olfactory lobe-optic gland axis is called the endocrine regulation center of cephalopods. However, little is known about the mechanism of the subpeduncle lobe/olfactory lobe-optic gland axis regulate the sexual maturation and post-reproductive death of Sepia esculenta Hoyle.

OBJECTIVES

The primary objective of this study was to provide basic information for revealing the mechanism of the subpeduncle lobe/olfactory lobe-optic axis regulating the rapid post-reproductive death of S. esculenta.

METHODS

In this paper, Illumina sequencing based transcriptome analysis was performed on the brain tissue of female S. esculenta in the three key developmental stages: growth stage (BG), spawning stage (BS), and post-reproductive death stage (BA).

RESULTS

A total of 66.19 Gb Illumina sequencing data were obtained. A comparative analysis of the three stages showed 2609, 3333, and 170 differentially expressed genes (DEGs) in BG-vs-BA, BG-vs-BA, and BS-vs-BA, respectively. The Gene Ontology (GO) enrichment analysis of DEGs revealed that the regulation of cyclin-dependent protein serine/threonine kinase activity, oxidative phosphorylation, and respiratory chain were significantly enriched. The significant enrichment analysis of the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway identified pathways associated with the regulation of death, such as the mammalian target of rapamycin (mTOR) signaling pathway, AMPK signaling pathway, oxidative phosphorylation, and cell cycle.

CONCLUSION

The post-reproductive death of S. esculenta was found to be a complex energy steady-state regulation network system. The mTOR acted as an energy receptor and had a key role in regulating energy homeostasis.

De Novo Transcriptome Assembly and Population Genetic Analyses for an Endangered Chinese Endemic Acer miaotaiense (Aceraceae)

Acer miaotaiense (P. C. Tsoong) is a rare and highly endangered plant in China. Because of the lack of genomic information and the limited number of available molecular markers, there are insufficient tools to determine the genetic diversity of this species. Here, 93,305 unigenes were obtained by multiple assembled contigs with a transcriptome sequencing program. Furthermore, 12,819 expressed sequence tag-derived simple sequence repeat (EST-SSR) markers were generated, 300 were randomly selected and synthesized, 19 primer pairs were identified as highly polymorphic (average number of alleles (N_a) = 8, expected heterozygosity (H_e) = 0.635, polymorphism information content (PIC) = 0.604) and were further used for population genetic analysis. All 261 samples were grouped into two genetic clusters by UPGMA, a principal component analyses and a STRUCTURE analyses. A moderate level of genetic differentiation (genetic differentiation index (F_st) = 0.059–0.116, gene flow = 1.904–3.993) among the populations and the major genetic variance (81.01%) within populations were revealed by the AMOVA. Based on the results, scientific conservation strategies should be established using in situ and ex situ conservation strategies. The study provides useful genetic information for the protection of precious wild resources and for further research on the origin and evolution of this endangered plant and its related species.

Transcriptome analysis of gene expression patterns during embryonic development in golden cuttlefish (Sepia esculenta)

Golden cuttlefish (Sepia esculenta) is an important economic species in China. Because of the rapid decline of its natural resource, researchers are exploring breeding technique for this species. The major obstacle that hinders artificial breeding of S. esculenta is the low larvae survival rate. Mortality is especially high during the mouth-opening stage. Investigating the embryogenesis before the first feed could provide theoretical guidance for reproduction control and breeding of S. esculenta and other Sepia species. In this study, we analyzed the dynamics of the S. esculenta transcriptome along different stages of embryonic development by mRNA-sEq. Our bioinformatics protocol identified 1492 differentially expressed genes (DEGs) across the early developmental stages. Gene ontology enrichment analysis showed that the DEGs were significantly involved in developmental processes and molecular functions, including chitin metabolic process, peptidase activity, catalytic activity, and calcium ion binding. Our results indicated that genes related to cuttlebone development and gene regulation functions were active during the early life phase of S. esculenta. Hierarchical clustering of the DEGs reflected the successiveness of the developmental stages, revealing that gene expression patterns of neighboring stages were similar. The DEG analysis allowed us to identify specific genes and relevant biological pathways to better understand the molecular mechanisms during each developmental stage. This study provides novel insights into the processes underlying the early developmental stages of S. esculenta. The transcriptomic data and identified genes will serve as valuable references for the developmental biology of this species and will help promote its aquaculture research.

The Gastric Ganglion of Octopus vulgaris: Preliminary Characterization of Gene- and Putative Neurochemical-Complexity, and the Effect of Aggregata octopiana Digestive Tract Infection on Gene Expression

The gastric ganglion is the largest visceral ganglion in cephalopods. It is connected to the brain and is implicated in regulation of digestive tract functions. Here we have investigated the neurochemical complexity (through in silico gene expression analysis and immunohistochemistry) of the gastric ganglion in Octopus vulgaris and tested whether the expression of a selected number of genes was influenced by the magnitude of digestive tract parasitic infection by Aggregata octopiana. Novel evidence was obtained for putative peptide and non-peptide neurotransmitters in the gastric ganglion: cephalotocin, corticotrophin releasing factor, FMRFamide, gamma amino butyric acid, 5-hydroxytryptamine, molluscan insulin-related peptide 3, peptide PRQFV-amide, and tachykinin-related peptide. Receptors for cholecystokininA and cholecystokininB, and orexin2 were also identified in this context for the first time. We report evidence for acetylcholine, dopamine, noradrenaline, octopamine, small cardioactive peptide related peptide, and receptors for cephalotocin and octopressin, confirming previous publications. The effects of Aggregata observed here extend those previously described by showing effects on the gastric ganglion; in animals with a higher level of infection, genes implicated in inflammation (NFκB, fascin, serpinB10 and the toll-like 3 receptor) increased their relative expression, but TNF-α gene expression was lower as was expression of other genes implicated in oxidative stress (i.e., superoxide dismutase, peroxiredoxin 6, and glutathione peroxidase). Elevated Aggregata levels in the octopuses corresponded to an increase in the expression of the cholecystokininA receptor and the small cardioactive peptide-related peptide. In contrast, we observed decreased relative expression of cephalotocin, dopamine β-hydroxylase, peptide PRQFV-amide, and tachykinin-related peptide genes. A discussion is provided on (i) potential roles of the various molecules in food intake regulation and digestive tract motility control and (ii) the difference in relative gene expression in the gastric ganglion in octopus with relatively high and low parasitic loads and the similarities to changes in the enteric innervation of mammals with digestive tract parasites. Our results provide additional data to the described neurochemical complexity of O. vulgaris gastric ganglion.