Gonadal Transcriptome Sequencing Analysis Reveals the Candidate Sex-Related Genes and Signaling Pathways in the East Asian Common Octopus, Octopus sinensis.

The East Asian common octopus (Octopus sinensis) is an economically important species among cephalopods. This species exhibits a strict dioecious and allogamous reproductive strategy, along with a phenotypic sexual dimorphism, where the third right arm differentiates into hectocotylus in males. However, our understanding of the molecular mechanisms that underlie sex determination and differentiation in this species remains limited. In the present study, we surveyed gene-expression profiles in the immature male and female gonads of O. sinensis based on the RNA-seq, and a total of 47.83 Gb of high-quality data were generated. Compared with the testis, we identified 8302 differentially expressed genes (DEGs) in the ovary, of which 4459 genes were up-regulated and 3843 genes were down-regulated. Based on the GO enrichment, many GO terms related to sex differentiation were identified, such as sex differentiation (GO: 0007548), sexual reproduction (GO: 0019953) and male sex differentiation (GO: 0046661). A KEGG classification analysis identified three conserved signaling pathways that related to sex differentiation, including the Wnt signaling pathway, TGF-ß signaling pathway and Notch signaling pathway. Additionally, 21 sex-related DEGs were selected, of which 13 DEGs were male-biased, including Dmrt1, Foxn5, Foxj1, Sox30, etc., and 8 DEGs were female-biased, including Sox14, Nanos3, ß-tubulin, Suh, etc. Ten DEGs were used to verify the expression patterns in the testis and ovary using the RT-qPCR method, and the results showed that the expression level shown by RT-qPCR was consistent with that from the RNA-seq, which confirmed the reliability of the transcriptome data. The results presented in this study will not only contribute to our understanding of sex-formation mechanisms in O. sinensis but also provide the foundational information for further investigating the molecular mechanisms that underline its gonadal development and facilitate the sustainable development of octopus artificial breeding.

Genome-wide identification of toll-like receptors in Octopus sinensis and expression analysis in response to different PAMPs stimulation.

Toll-like receptors (TLRs) are one of the extensively studied pattern recognition receptors (PRRs) and play crucial roles in the immune responses of vertebrates and invertebrates. In this study, 14 TLR genes were identified from the genome-wide data of Octopus sinensis. Protein structural domain analysis showed that most TLR proteins had three main structural domains: extracellular leucine-rich repeats (LRR), transmembrane structural domains, and intracellular Toll/IL-1 receptor domain (TIR). The results of subcellular localization prediction showed that the TLRs of O. sinensis were mainly located on the plasma membrane. The results of quantitative real-time PCR (qPCR) showed that the detected TLR genes were differentially expressed in the hemolymph, white bodies, hepatopancreas, gills, gill heart, intestine, kidney, and salivary gland of O. sinensis. Furthermore, the present study investigated the expression changes of O. sinensis TLR genes in hemolymph, white bodies, gills, and hepatopancreas in different phases (6 h, 12 h, 24 h, 48 h) after stimulation with PGN, poly(I: C) and Vibrio parahaemolyticus. The expression of most of the TLR genes was upregulated at different time points after infection with pathogens or stimulation with PAMPs, a few genes were unchanged or even down-regulated, and many of the TLR genes were much higher after V. parahaemolyticus infection than after PGN and poly(I:C) stimulation. The results of this study contribute to a better understanding of the molecular immune mechanisms of O. sinensis TLRs genes in resistance to pathogen stimulation.

Genomic evidence for West Antarctic Ice Sheet collapse during the Last Interglacial.

The marine-based West Antarctic Ice Sheet (WAIS) is considered vulnerable to irreversible collapse under future climate trajectories, and its tipping point may lie within the mitigated warming scenarios of 1.5° to 2°C of the United Nations Paris Agreement. Knowledge of ice loss during similarly warm past climates could resolve this uncertainty, including the Last Interglacial when global sea levels were 5 to 10 meters higher than today and global average temperatures were 0.5° to 1.5°C warmer than preindustrial levels. Using a panel of genome-wide, single-nucleotide polymorphisms of a circum-Antarctic octopus, we show persistent, historic signals of gene flow only possible with complete WAIS collapse. Our results provide the first empirical evidence that the tipping point of WAIS loss could be reached even under stringent climate mitigation scenarios.

The complete mitochondrial genome of Octopus vulgaris.

BACKGROUND: The Octopus vulgaris species complex consists of numerous morphologically similar but genetically distinct species. The current publicly available mitogenome of this species has been generated from a specimen collected from Tsukiji Fish Market, Tokyo, Japan. Octopus from the northwestern Pacific Ocean are now considered to be a separate species, Octopus sinensis. For this reason, we hypothesised that the current record of O. vulgaris was sequenced from a specimen of O. sinensis. Here, we sequenced the first complete mitogenome of a specimen of Octopus vulgaris sensu stricto that was collected from the species' confirmed distribution areas in northeastern Atlantic. METHODS AND RESULTS: The complete mitogenome was assembled de novo and annotated using 250 bp paired-end sequences. A single circular contig 15,655 bp in length with a mean read coverage of 1089 reads was reconstructed. The annotation pipeline identified 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNA) and two ribosomal RNAs. A maximum likelihood phylogenetic tree recovered the assembled mitogenome as the sister taxon of a monophyletic group comprising O. sinensis and the previously published mitogenome of "O. vulgaris" from Japan. This confirms that the latter was a Japanese specimen of O. sinensis. CONCLUSION: The mitogenome sequenced here is the first to be published for Octopus vulgaris sensu stricto. It represents an important first step in genetics-informed research on the evolution, conservation, and management of this commercially important species.

Morphology, Histology, and Transcriptome Analysis of Gonadal Development in Octopus minor (Sasaki, 1920).

Octopus minor is an economically important species, but little is known about the histological pattern and regulatory mechanisms during gonadal development. In this study, we investigated the annual changes in total body weight (TW), gonad somatic index (GSI), gonadal histological features, and transcriptome of O. minor. The results indicated that both females and males showed a similar TW trend. The GSI peaked in June in females, while it remained constant at around 3% in males. Nine and four histological stages were observed in ovaries and testes, respectively. Our field sampling results implied that O. minor might have overwintering periods for both eggs and larvae. Transcriptome analysis revealed that a total of 1095 and 2468 genes were significantly expressed during ovarian and testicular development, separately. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis displayed that 126 GO terms and 5 KEGG pathways were significantly enriched in the ovarian group of advanced vitellogenic oocytes vs vitellogenic oocytes (AVO vs VO). The pathways "Ribosomal", "Cell cycle", and "Progesterone-mediated oocyte maturation" were predicted to promote yolk deposition. Additionally, the testicular comparison group of spent vs mature (Spent vs Mature) showed significant enrichment in 674 GO terms and 13 KEGG pathways, suggesting that energy metabolism and cell repair pathways may be involved in the spermatogenesis process. This work revealed the development process of the gonads and shed light on the potential regulatory pathways of O. minor, providing novel insights and laying a molecular basis for artificial breeding.

A chromosome-level reference genome for the common octopus, Octopus vulgaris (Cuvier, 1797).

Cephalopods are emerging animal models and include iconic species for studying the link between genomic innovations and physiological and behavioral complexities. Coleoid cephalopods possess the largest nervous system among invertebrates, both for cell counts and brain-to-body ratio. Octopus vulgaris has been at the center of a long-standing tradition of research into diverse aspects of cephalopod biology, including behavioral and neural plasticity, learning and memory recall, regeneration, and sophisticated cognition. However, no chromosome-scale genome assembly was available for O. vulgaris to aid in functional studies. To fill this gap, we sequenced and assembled a chromosome-scale genome of the common octopus, O. vulgaris. The final assembly spans 2.8 billion basepairs, 99.34% of which are in 30 chromosome-scale scaffolds. Hi-C heatmaps support a karyotype of 1n = 30 chromosomes. Comparisons with other octopus species' genomes show a conserved octopus karyotype and a pattern of local genome rearrangements between species. This new chromosome-scale genome of O. vulgaris will further facilitate research in all aspects of cephalopod biology, including various forms of plasticity and the neural machinery underlying sophisticated cognition, as well as an understanding of cephalopod evolution.

Octopod Hox genes and cephalopod plesiomorphies.

Few other invertebrates captivate our attention as cephalopods do. Octopods, cuttlefish, and squids amaze with their behavior and sophisticated body plans that belong to the most intriguing among mollusks. Little is, however, known about their body plan formation and the role of Hox genes. The latter homeobox genes pattern the anterior-posterior body axis and have only been studied in a single decapod species so far. Here, we study developmental Hox and ParaHox gene expression in Octopus vulgaris. Hox genes are expressed in a near-to-staggered fashion, among others in homologous organs of cephalopods such as the stellate ganglia, the arms, or funnel. As in other mollusks Hox1 is expressed in the nascent octopod shell rudiment. While ParaHox genes are expressed in an evolutionarily conserved fashion, Hox genes are also expressed in some body regions that are considered homologous among mollusks such as the cephalopod arms and funnel with the molluscan foot. We argue that cephalopod Hox genes are recruited to a lesser extent into the formation of non-related organ systems than previously thought and emphasize that despite all morphological innovations molecular data still reveal the ancestral molluscan heritage of cephalopods.

Chilling with cephalopods: Temperature-responsive RNA editing in octopus and squid.

The molecular mechanisms that generate the developmental and physiological complexity found within cephalopods are not well understood. In this issue of Cell, Birk et al. and Rangan and Reck-Peterson show that cephalopods differentially edit their RNA in response to temperature changes and that this editing has consequences on protein function.

Temperature-dependent RNA editing in octopus extensively recodes the neural proteome.

In poikilotherms, temperature changes challenge the integration of physiological function. Within the complex nervous systems of the behaviorally sophisticated coleoid cephalopods, these problems are substantial. RNA editing by adenosine deamination is a well-positioned mechanism for environmental acclimation. We report that the neural proteome of Octopus bimaculoides undergoes massive reconfigurations via RNA editing following a temperature challenge. Over 13,000 codons are affected, and many alter proteins that are vital for neural processes. For two highly temperature-sensitive examples, recoding tunes protein function. For synaptotagmin, a key component of Ca2+-dependent neurotransmitter release, crystal structures and supporting experiments show that editing alters Ca2+ binding. For kinesin-1, a motor protein driving axonal transport, editing regulates transport velocity down microtubules. Seasonal sampling of wild-caught specimens indicates that temperature-dependent editing occurs in the field as well. These data show that A-to-I editing tunes neurophysiological function in response to temperature in octopus and most likely other coleoids.

The Effect of Prohibitins on Mitochondrial Function during Octopus tankahkeei Spermiogenesis.

Mitochondria are essential for spermiogenesis. Prohibitins (PHBs; prohibitin 1, PHB1 or PHB, and prohibitin 2, PHB2) are evolutionarily conserved and ubiquitously expressed mitochondrial proteins that act as scaffolds in the inner mitochondrial membrane. In this study, we analyzed the molecular structure and dynamic expression characteristics of Ot-PHBs, observed the colocalization of Ot-PHB1 with mitochondria and polyubiquitin, and studied the effect of phb1 knockdown on mitochondrial DNA (mtDNA) content, reactive oxygen species (ROS) levels, and apoptosis-related gene expression in spermatids. Our aim was to explore the effect of Ot-PHBs on mitochondrial function during the spermiogenesis of Octopus tankahkeei (O. tankahkeei), an economically important species in China. The predicted Ot-PHB1/PHB2 proteins contained an N-terminal transmembrane, a stomatin/prohibitin/flotillin/HflK/C (SPFH) domain (also known as the prohibitin domain), and a C-terminal coiled-coil domain. Ot-phb1/phb2 mRNA were widely expressed in the different tissues, with elevated expression in the testis. Further, Ot-PHB1 and Ot-PHB2 were highly colocalized, suggesting that they may function primarily as an Ot-PHB compiex in O. tankahkeei. Ot-PHB1 proteins were mainly expressed and localized in mitochondria during spermiogenesis, implying that their function may be localized to the mitochondria. In addition, Ot-PHB1 was colocalized with polyubiquitin during spermiogenesis, suggesting that it may be a polyubiquitin substrate that regulates mitochondrial ubiquitination during spermiogenesis to ensure mitochondrial quality. To further investigate the effect of Ot-PHBs on mitochondrial function, we knocked down Ot-phb1 and observed a decrease in mtDNA content, along with increases in ROS levels and the expressions of mitochondria-induced apoptosis-related genes bax, bcl2, and caspase-3 mRNA. These findings indicate that PHBs might influence mitochondrial function by maintaining mtDNA content and stabilizing ROS levels; in addition, PHBs might affect spermatocyte survival by regulating mitochondria-induced apoptosis during spermiogenesis in O. tankahkeei.

SNP data reveals the complex and diverse evolutionary history of the blue-ringed octopus genus (Octopodidae: Hapalochlaena) in the Asia-Pacific.

The blue-ringed octopus species complex (Hapalochlaena spp.), known to occur from Southern Australia to Japan, currently contains four formally described species (Hapalochlaena maculosa, Hapalochlaena fasciata, Hapalochlaena lunulata and Hapalochlaena nierstraszi). These species are distinguished based on morphological characters (iridescent blue rings and/or lines) along with reproductive strategies. However, the observation of greater morphological diversity than previously captured by the current taxonomic framework indicates that a revision is required. To examine species boundaries within the genus we used mitochondrial (12S rRNA, 16S rRNA, cytochrome c oxidase subunit 1 [COI], cytochrome c oxidase subunit 3 [COIII] and cytochrome b [Cytb]) and genome-wide SNP data (DaRT seq) from specimens collected across its geographic range including variations in depth from 3 m to >100 m. This investigation indicates substantially greater species diversity present within the genus Hapalochlaena than is currently described. We identified 10,346 SNPs across all locations, which when analysed support a minimum of 11 distinct clades. Bayesian phylogenetic analysis of the mitochondrial COI gene on a more limited sample set dates the diversification of the genus to âˆ¼30 mya and corroborates eight of the lineages indicated by the SNP analyses. Furthermore, we demonstrate that the diagnostic lined patterning of H. fasciata found in North Pacific waters and NSW, Australia is polyphyletic and therefore likely the result of convergent evolution. Several "deep water" (>100 m) lineages were also identified in this study with genetic convergence likely to be driven by external selective pressures. Examination of morphological traits, currently being undertaken in a parallel morphological study, is required to describe additional species within the complex.

Brain compartmentalization based on transcriptome analyses and its gene expression in Octopus minor.

Coleoid cephalopods have a high intelligence, complex structures, and large brain. The cephalopod brain is divided into supraesophageal mass, subesophageal mass and optic lobe. Although much is known about the structural organization and connections of various lobes of octopus brain, there are few studies on the brain of cephalopod at the molecular level. In this study, we demonstrated the structure of an adult Octopus minor brain by histomorphological analyses. Through visualization of neuronal and proliferation markers, we found that adult neurogenesis occurred in the vL and posterior svL. We also obtained specific 1015 genes by transcriptome of O. minor brain and selected OLFM3, NPY, GnRH, and GDF8 genes. The expression of genes in the central brain showed the possibility of using NPY and GDF8 as molecular marker of compartmentation in the central brain. This study will provide useful information for establishing a molecular atlas of cephalopod brain.

Gene network analyses of larvae under different egg-protecting behaviors provide novel insights into immune response mechanisms of Amphioctopus fangsiao.

Amphioctopus fangsiao was a representative economic species in cephalopods, which was vulnerable to marine bacteria. Vibrio anguillarum was a highly infectious pathogen that have recently been found to infect A. fangsiao and inhibit its growth and development. There were significant differences in the immune response mechanisms between egg-protected and egg-unprotected larvae. To explore larval immunity under different egg-protecting behaviors, we infected A. fangsiao larvae with V. anguillarum for 24 h and analyzed the transcriptome data about egg-protected and egg-unprotected larvae infected with 0, 4, 12, and 24 h using weighted gene co-expression networks (WGCNA) and protein-protein interaction (PPI) networks. Network analyses revealed a series of immune response processes after infection, and identified six key modules and multiple immune-related hub genes. Meanwhile, we found that ZNF family, such as ZNF32, ZNF160, ZNF271, ZNF479, and ZNF493 might play significant roles in A. fangsiao immune response processes. We first creatively combined WGCNA and PPI network analysis to deeply explore the immune response mechanisms of A. fangsiao larvae with different egg-protecting behaviors. Our results provided further insights into the immunity of V. anguillarum infected invertebrates, and laid the foundation for exploring the immune differences among cephalopods with different egg protecting behaviors.

Genome skimming elucidates the evolutionary history of Octopoda.

Phylogenies for Octopoda have, until now, been based on morphological characters or a few genes. Here we provide the complete mitogenomes and the nuclear 18S and 28S ribosomal genes of twenty Octopoda specimens, comprising 18 species of Cirrata and Incirrata, representing 13 genera and all five putative families of Cirrata (Cirroctopodidae, Cirroteuthidae, Grimpoteuthidae, Opisthoteuthidae and Stauroteuthidae) and six families of Incirrata (Amphitretidae, Argonautidae, Bathypolypodidae, Eledonidae, Enteroctopodidae, and Megaleledonidae) which were assembled using genome skimming. Phylogenetic trees were built using Maximum Likelihood and Bayesian Inference with several alignment matrices. All mitochondrial genomes had the 'typical' genome composition and gene order previously reported for octopodiforms, except Bathypolypus ergasticus, which appears to lack ND5, two tRNA genes that flank ND5 and two other tRNA genes. Argonautoidea was revealed as sister to Octopodidae by the mitochondrial protein-coding gene dataset, however, it was recovered as sister to all other incirrate octopods with strong support in an analysis using nuclear rRNA genes. Within Cirrata, our study supports two existing classifications suggesting neither is likely in conflict with the true evolutionary history of the suborder. Genome skimming is useful in the analysis of phylogenetic relationships within Octopoda; inclusion of both mitochondrial and nuclear data may be key.

Genetic monitoring on the world's first MSC eco-labeled common octopus (O. vulgaris) fishery in western Asturias, Spain.

Octopus vulgaris (Cuvier, 1797) is a cephalopod species with great economic value. In western Asturias (northwest of Spain), O. vulgaris artisanal fisheries are relatively well monitored and conditionally eco-labeled by the Marine Stewardship Council (MSC). Despite this, the Asturian octopus stocks have not been genetically assessed so far. In order to improve the current fishery plan and contrast the octopus eco-label validity in Asturias, 539 individuals from five regions of the O. vulgaris geographic distribution, including temporal samplings in Asturias, were collected and genotyped at thirteen microsatellite loci. All the samples under analysis were in agreement with Hardy-Weinberg expectations. Spatial levels of genetic differentiation were estimated using F-statistics, multidimensional scaling, and Bayesian analyses. Results suggested that the O. vulgaris consists of at least four genetically different stocks coming from two ancestral lineages. In addition, temporal analyses showed stability in terms of genetic variation and high NE (> 50) for several generations in different localities within Asturias, pointing out to indeed sustainable fishery exploitation levels. Even though, the current Asturias fishery plan shows no significant genetic damages to the stocks, the regional-specific management plans need systematic genetic monitoring schemes as part of an efficient and preventive regional fishery regulation strategy.

The gold-ringed octopus (Amphioctopus fangsiao) genome and cerebral single-nucleus transcriptomes provide insights into the evolution of karyotype and neural novelties.

BACKGROUND: Coleoid cephalopods have distinctive neural and morphological characteristics compared to other invertebrates. Early studies reported massive genomic rearrangements occurred before the split of octopus and squid lineages (Proc Natl Acad Sci U S A 116:3030-5, 2019), which might be related to the neural innovations of their brain, yet the details remain elusive. Here we combine genomic and single-nucleus transcriptome analyses to investigate the octopod chromosome evolution and cerebral characteristics. RESULTS: We present a chromosome-level genome assembly of a gold-ringed octopus, Amphioctopus fangsiao, and a single-nucleus transcriptome of its supra-esophageal brain. Chromosome-level synteny analyses estimate that the chromosomes of the ancestral octopods experienced multiple chromosome fission/fusion and loss/gain events by comparing with the nautilus genome as outgroup, and that a conserved genome organization was detected during the evolutionary process from the last common octopod ancestor to their descendants. Besides, protocadherin, GPCR, and C2H2 ZNF genes are thought to be highly related to the neural innovations in cephalopods (Nature 524:220-4, 2015), and the chromosome analyses pinpointed several collinear modes of these genes on the octopod chromosomes, such as the collinearity between PCDH and C2H2 ZNF, as well as between GPCR and C2H2 ZNF. Phylogenetic analyses show that the expansion of the octopod protocadherin genes is driven by a tandem-duplication mechanism on one single chromosome, including two separate expansions at 65 million years ago (Ma) and 8-14 Ma, respectively. Furthermore, we identify eight cell types (i.e., cholinergic and glutamatergic neurons) in the supra-esophageal brain of A. fangsiao, and the single-cell expression analyses reveal the co-expression of protocadherin and GPCR in specific neural cells, which may contribute to the neural development and signal transductions in the octopod brain. CONCLUSIONS: The octopod genome analyses reveal the dynamic evolutionary history of octopod chromosomes and neural-related gene families. The single-nucleus transcriptomes of the supra-esophageal brain indicate their cellular heterogeneities and functional interactions with other tissues (i.e., gill), which provides a foundation for further octopod cerebral studies.

MicroRNAs are deeply linked to the emergence of the complex octopus brain.

Soft-bodied cephalopods such as octopuses are exceptionally intelligent invertebrates with a highly complex nervous system that evolved independently from vertebrates. Because of elevated RNA editing in their nervous tissues, we hypothesized that RNA regulation may play a major role in the cognitive success of this group. We thus profiled messenger RNAs and small RNAs in three cephalopod species including 18 tissues of the Octopus vulgaris. We show that the major RNA innovation of soft-bodied cephalopods is an expansion of the microRNA (miRNA) gene repertoire. These evolutionarily novel miRNAs were primarily expressed in adult neuronal tissues and during the development and had conserved and thus likely functional target sites. The only comparable miRNA expansions happened, notably, in vertebrates. Thus, we propose that miRNAs are intimately linked to the evolution of complex animal brains.

Cell type diversity in a developing octopus brain.

Octopuses are mollusks that have evolved intricate neural systems comparable with vertebrates in terms of cell number, complexity and size. The brain cell types that control their sophisticated behavioral repertoire are still unknown. Here, we profile the cell diversity of the paralarval Octopus vulgaris brain to build a cell type atlas that comprises mostly neural cells, but also multiple glial subtypes, endothelial cells and fibroblasts. We spatially map cell types to the vertical, subesophageal and optic lobes. Investigation of cell type conservation reveals a shared gene signature between glial cells of mouse, fly and octopus. Genes related to learning and memory are enriched in vertical lobe cells, which show molecular similarities with Kenyon cells in Drosophila. We construct a cell type taxonomy revealing transcriptionally related cell types, which tend to appear in the same brain region. Together, our data sheds light on cell type diversity and evolution in the octopus brain.

De novo transcriptome reconstruction in aquacultured early life stages of the cephalopod Octopus vulgaris.

Cephalopods have been considered enigmatic animals that have attracted the attention of scientists from different areas of expertise. However, there are still many questions to elucidate the way of life of these invertebrates. The aim of this study is to construct a reference transcriptome in Octopus vulgaris early life stages to enrich existing databases and provide a new dataset that can be reused by other researchers in the field. For that, samples from different developmental stages were combined including embryos, newly-hatched paralarvae, and paralarvae of 10, 20 and 40 days post-hatching. Additionally, different dietary and rearing conditions and pathogenic infections were tested. At least three biological replicates were analysed per condition and submitted to RNA-seq analysis. All sequencing reads from experimental conditions were combined in a single dataset to generate a reference transcriptome assembly that was functionally annotated. The number of reads aligned to this reference was counted to estimate the transcript abundance in each sample. This dataset compiled a complete reference for future transcriptomic studies in O. vulgaris.