Emergence of novel cephalopod gene regulation and expression through large-scale genome reorganization

Pan-evolutionary and regulatory genome architecture delineated by an integrated macro- and microsynteny approach

The forthcoming massive genome data generated by the Earth BioGenome Project will open up a new era of comparative genomics, for which genome synteny analysis provides an important framework. Profiling genome synteny represents an essential step in elucidating genome architecture, regulatory blocks/elements and their evolutionary history. Here we describe PanSyn, ( https://github.com/yhw320/PanSyn ), the most comprehensive and up-to-date genome synteny pipeline, providing step-by-step instructions and application examples to demonstrate its usage. PanSyn inherits both basic and advanced functions from existing popular tools, offering a user-friendly, highly customized approach for genome macrosynteny analysis and integrated pan-evolutionary and regulatory analysis of genome architecture, which are not yet available in public synteny software or tools. The advantages of PanSyn include: (i) advanced microsynteny analysis by functional profiling of microsynteny genes and associated regulatory elements; (ii) comprehensive macrosynteny analysis, including the inference of karyotype evolution from ancestors to extant species; and (iii) functional integration of microsynteny and macrosynteny for pan-evolutionary profiling of genome architecture and regulatory blocks, as well as integration with external functional genomics datasets from three- or four-dimensional genome and ENCODE projects. PanSyn requires basic knowledge of the Linux environment and Perl programming language and the ability to access a computer cluster, especially for large-scale genomic comparisons. Our protocol can be easily implemented by a competent graduate student or postdoc and takes several days to weeks to execute for dozens to hundreds of genomes. PanSyn provides yet the most comprehensive and powerful tool for integrated evolutionary and functional genomics.

Cephalopod Sex Determination and its Ancient Evolutionary Origin Revealed by Chromosome-level Assembly of the California Two-Spot Octopus

Sex chromosomes are critical elements of sexual reproduction in many animal and plant taxa, however they show incredible diversity and rapid turnover even within clades. Here, using a chromosome-level assembly generated with long read sequencing, we report the first evidence for genetic sex determination in cephalopods. We have uncovered a sex chromosome in California two-spot octopus (Octopus bimaculoides) in which males/females show ZZ/ZO karyotypes respectively. We show that the octopus Z chromosome is an evolutionary outlier with respect to divergence and repetitive element content as compared to other chromosomes and that it is present in all coleoid cephalopods that we have examined. Our results suggest that the cephalopod Z chromosome originated between 455 and 248 million years ago and has been conserved to the present, making it the among the oldest conserved animal sex chromosomes known.

Gene modelling and annotation for the Hawaiian bobtail squid, Euprymna scolopes

Coleoid cephalopods possess numerous complex, species-specific morphological and behavioural adaptations, e.g., a uniquely structured nervous system that is the largest among the invertebrates. The Hawaiian bobtail squid (Euprymna scolopes) is one of the most established cephalopod species. With its recent publication of the chromosomal-scale genome assembly and regulatory genomic data, it also emerges as a key model for cephalopod gene regulation and evolution. However, the latest genome assembly has been lacking a native gene model set. Our manuscript describes the generation of new long-read transcriptomic data and, made using this combined with a plethora of publicly available transcriptomic and protein sequence data, a new reference annotation for E. scolopes.

Cephalopod-omics: Emerging Fields and Technologies in Cephalopod Biology

Few animal groups can claim the level of wonder that cephalopods instill in the minds of researchers and the general public. Much of cephalopod biology, however, remains unexplored: the largest invertebrate brain, difficult husbandry conditions, and complex (meta-)genomes, among many other things, have hindered progress in addressing key questions. However, recent technological advancements in sequencing, imaging, and genetic manipulation have opened new avenues for exploring the biology of these extraordinary animals. The cephalopod molecular biology community is thus experiencing a large influx of researchers, emerging from different fields, accelerating the pace of research in this clade. In the first post-pandemic event at the Cephalopod International Advisory Council (CIAC) conference in April 2022, over 40 participants from all over the world met and discussed key challenges and perspectives for current cephalopod molecular biology and evolution. Our particular focus was on the fields of comparative and regulatory genomics, gene manipulation, single-cell transcriptomics, metagenomics, and microbial interactions. This article is a result of this joint effort, summarizing the latest insights from these emerging fields, their bottlenecks, and potential solutions. The article highlights the interdisciplinary nature of the cephalopod-omics community and provides an emphasis on continuous consolidation of efforts and collaboration in this rapidly evolving field.

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.

Octopus vulgaris Exhibits Interindividual Differences in Behavioural and Problem-Solving Performance

By presenting individual Octopus vulgaris with an extractive foraging problem with a puzzle box, we examined the possible correlation between behavioural performances (e.g., ease of adaptation to captive conditions, prevalence of neophobic and neophilic behaviours, and propensity to learn individually or by observing conspecifics), biotic (body and brain size, age, sex) and abiotic (seasonality and place of origin) factors. We found more neophilic animals showing shorter latencies to approach the puzzle box and higher probability of solving the task; also, shorter times to solve the task were correlated with better performance on the individual learning task. However, the most neophilic octopuses that approached the puzzle box more quickly did not reach the solution earlier than other individuals, suggesting that strong neophilic tendency may lead to suboptimal performance at some stages of the problem-solving process. In addition, seasonal and environmental characteristics of location of origin appear to influence the rate of expression of individual traits central to problem solving. Overall, our analysis provides new insights into the traits associated with problem solving in invertebrates and highlights the presence of adaptive mechanisms that promote population-level changes in octopuses' behavioural traits.

Functional genomics in Spiralia

Our understanding of the mechanisms that modulate gene expression in animals is strongly biased by studying a handful of model species that mainly belong to three groups: Insecta, Nematoda and Vertebrata. However, over half of the animal phyla belong to Spiralia, a morphologically and ecologically diverse animal clade with many species of economic and biomedical importance. Therefore, investigating genome regulation in this group is central to uncovering ancestral and derived features in genome functioning in animals, which can also be of significant societal impact. Here, we focus on five aspects of gene expression regulation to review our current knowledge of functional genomics in Spiralia. Although some fields, such as single-cell transcriptomics, are becoming more common, the study of chromatin accessibility, DNA methylation, histone post-translational modifications and genome architecture are still in their infancy. Recent efforts to generate chromosome-scale reference genome assemblies for greater species diversity and optimise state-of-the-art approaches for emerging spiralian research systems will address the existing knowledge gaps in functional genomics in this animal group.

Emerging questions on the mechanisms and dynamics of 3D genome evolution in spiralians

Information on how 3D genome topology emerged in animal evolution, how stable it is during development, its role in the evolution of phenotypic novelties and how exactly it affects gene expression is highly debated. So far, data to address these questions are lacking with the exception of a few key model species. Several gene regulatory mechanisms have been proposed, including scenarios where genome topology has little to no impact on gene expression, and vice versa. The ancient and diverse clade of spiralians may provide a crucial testing ground for such mechanisms. Sprialians have followed distinct evolutionary trajectories, with some clades experiencing genome expansions and/or large-scale genome rearrangements, and others undergoing genome contraction, substantially impacting their size and organisation. These changes have been associated with many phenotypic innovations in this clade. In this review, we describe how emerging genome topology data, along with functional tools, allow for testing these scenarios and discuss their predicted outcomes.

Comparison of Blastema Formation after Injury in Two Cephalopod Species

Regeneration is the ability to functionally replace significant amounts of lost tissue or whole appendages like arms, limbs or tentacles. The amount of tissue that can be regenerated varies among species, but regeneration is found in both invertebrate and vertebrate animals. Cephalopods have been broadly reported in the literature to regenerate their arms. There are over 800 species of Cephalopod; however, regeneration has only been documented in the literature in a few species (1). Here we compare arm regeneration in two species of cephalopod, the Octopus bimaculoides and the hummingbird bobtail squid Euprymna berryi.

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.

Shifting landscapes: the role of 3D genomic organizations in gene regulatory strategies

3D genome folding enables the physical storage of chromosomes into the compact volume of a cell's nucleus, allows for the accurate segregation of chromatin to daughter cells, and has been shown to be tightly coupled to the way in which genetic information is converted into transcriptional programs [1-3]. Importantly, this link between chromatin architecture and gene regulation is a selectable feature in which modifications to chromatin organization accompany, or perhaps even drive the establishment of new regulatory strategies with enduring impacts on animal body plan complexity. Here, we discuss the nature of different 3D genome folding systems found across the tree of life, with particular emphasis on metazoans, and the relative influence of these systems on gene regulation. We suggest how the properties of these folding systems have influenced regulatory strategies employed by different lineages and may have catalyzed the partitioning and specialization of genetic programs that enabled multicellularity and organ-grade body plan complexity.

Emergence of novel genomic regulatory regions associated with light-organ development in the bobtail squid

Light organs (LO) with symbiotic bioluminescent bacteria are hallmarks of many bobtail squid species. These organs possess structural and functional features to modulate light, analogous to those found in coleoid eyes. Previous studies identified four transcription factors and modulators (SIX, EYA, PAX6, DAC) associated with both eyes and light organ development, suggesting co-option of a highly conserved gene regulatory network. Using available topological, open chromatin, and transcriptomic data, we explore the regulatory landscape around the four transcription factors as well as genes associated with LO and shared LO/eye expression. This analysis revealed several closely associated and putatively co-regulated genes. Comparative genomic analyses identified distinct evolutionary origins of these putative regulatory associations, with the DAC locus showing a unique topological and evolutionarily recent organization. We discuss different scenarios of modifications to genome topology and how these changes may have contributed to the evolutionary emergence of the light organ.

Transcriptome-wide selection and validation of a solid set of reference genes for gene expression studies in the cephalopod mollusk Octopus vulgaris

Octopus vulgaris is a cephalopod mollusk and an active marine predator that has been at the center of a number of studies focused on the understanding of neural and biological plasticity. Studies on the machinery involved in e.g., learning and memory, regeneration, and neuromodulation are required to shed light on the conserved and/or unique mechanisms that these animals have evolved. Analysis of gene expression is one of the most essential means to expand our understanding of biological machinery, and the selection of an appropriate set of reference genes is the prerequisite for the quantitative real-time polymerase chain reaction (qRT-PCR). Here we selected 77 candidate reference genes (RGs) from a pool of stable and relatively high-expressed transcripts identified from the full-length transcriptome of O. vulgaris, and we evaluated their expression stabilities in different tissues through geNorm, NormFinder, Bestkeeper, Delta-CT method, and RefFinder. Although various algorithms provided different assemblages of the most stable reference genes for the different kinds of tissues tested here, a comprehensive ranking revealed RGs specific to the nervous system (Ov-RNF7 and Ov-RIOK2) and Ov-EIF2A and Ov-CUL1 across all considered tissues. Furthermore, we validated RGs by assessing the expression profiles of nine target genes (Ov-Naa15, Ov-Ltv1, Ov-CG9286, Ov-EIF3M, Ov-NOB1, Ov-CSDE1, Ov-Abi2, Ov-Homer2, and Ov-Snx20) in different areas of the octopus nervous system (gastric ganglion, as control). Our study allowed us to identify the most extensive set of stable reference genes currently available for the nervous system and appendages of adult O. vulgaris.

Evolution of 3D chromatin organization at different scales

Most animal genomes fold in 3D chromatin domains called topologically associated domains (TADs) that facilitate interactions between cis-regulatory elements (CREs) and promoters. Owing to their critical role in the control of developmental gene expression, we explore how TADs have shaped animal evolution. In the light of recent studies that profile TADs in disparate animal lineages, we discuss their phylogenetic distribution and the mechanisms that underlie their formation. We present evidence indicating that TADs are plastic entities composed of genomic strata of different ages: ancient cores are combined with newer regions and brought into extant TADs through genomic rearrangements. We highlight that newly incorporated TAD strata enable the establishment of new CRE-promoter interactions and in turn new expression patterns that can drive phenotypical innovation. We further highlight how subtle changes in chromatin folding may fine-tune the expression levels of developmental genes and hold a potential for evolutionary significance.

Deciphering regeneration through non-model animals: A century of experiments on cephalopod mollusks and an outlook at the future

The advent of marine stations in the last quarter of the 19th Century has given biologists the possibility of observing and experimenting upon myriad marine organisms. Among them, cephalopod mollusks have attracted great attention from the onset, thanks to their remarkable adaptability to captivity and a great number of biologically unique features including a sophisticate behavioral repertoire, remarkable body patterning capacities under direct neural control and the complexity of nervous system rivalling vertebrates. Surprisingly, the capacity to regenerate tissues and complex structures, such as appendages, albeit been known for centuries, has been understudied over the decades. Here, we will first review the limited in number, but fundamental studies on the subject published between 1920 and 1970 and discuss what they added to our knowledge of regeneration as a biological phenomenon. We will also speculate on how these relate to their epistemic and disciplinary context, setting the base for the study of regeneration in the taxon. We will then frame the peripherality of cephalopods in regeneration studies in relation with their experimental accessibility, and in comparison, with established models, either simpler (such as planarians), or more promising in terms of translation (urodeles). Last, we will explore the potential and growing relevance of cephalopods as prospective models of regeneration today, in the light of the novel opportunities provided by technological and methodological advances, to reconsider old problems and explore new ones. The recent development of cutting-edge technologies made available for cephalopods, like genome editing, is allowing for a number of important findings and opening the way toward new promising avenues. The contribution offered by cephalopods will increase our knowledge on regenerative mechanisms through cross-species comparison and will lead to a better understanding of the complex cellular and molecular machinery involved, shedding a light on the common pathways but also on the novel strategies different taxa evolved to promote regeneration of tissues and organs. Through the dialogue between biological/experimental and historical/contextual perspectives, this article will stimulate a discussion around the changing relations between availability of animal models and their specificity, technical and methodological developments and scientific trends in contemporary biology and medicine.

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.

Epigenetic machinery is functionally conserved in cephalopods

Background

Epigenetic regulatory mechanisms are divergent across the animal kingdom, yet these mechanisms are not well studied in non-model organisms. Unique features of cephalopods make them attractive for investigating behavioral, sensory, developmental, and regenerative processes, and recent studies have elucidated novel features of genome organization and gene and transposon regulation in these animals. However, it is not known how epigenetics regulates these interesting cephalopod features. We combined bioinformatic and molecular analysis of Octopus bimaculoides to investigate the presence and pattern of DNA methylation and examined the presence of DNA methylation and 3 histone post-translational modifications across tissues of three cephalopod species.

Results

We report a dynamic expression profile of the genes encoding conserved epigenetic regulators, including DNA methylation maintenance factors in octopus tissues. Levels of 5-methyl-cytosine in multiple tissues of octopus, squid, and bobtail squid were lower compared to vertebrates. Whole genome bisulfite sequencing of two regions of the brain and reduced representation bisulfite sequencing from a hatchling of O. bimaculoides revealed that less than 10% of CpGs are methylated in all samples, with a distinct pattern of 5-methyl-cytosine genome distribution characterized by enrichment in the bodies of a subset of 14,000 genes and absence from transposons. Hypermethylated genes have distinct functions and, strikingly, many showed similar expression levels across tissues while hypomethylated genes were silenced or expressed at low levels. Histone marks H3K27me3, H3K9me3, and H3K4me3 were detected at different levels across tissues of all species.

Conclusions

Our results show that the DNA methylation and histone modification epigenetic machinery is conserved in cephalopods, and that, in octopus, 5-methyl-cytosine does not decorate transposable elements, but is enriched on the gene bodies of highly expressed genes and could cooperate with the histone code to regulate tissue-specific gene expression.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01404-1.

Genome and transcriptome mechanisms driving cephalopod evolution

Cephalopods are known for their large nervous systems, complex behaviors and morphological innovations. To investigate the genomic underpinnings of these features, we assembled the chromosomes of the Boston market squid, Doryteuthis (Loligo) pealeii, and the California two-spot octopus, Octopus bimaculoides, and compared them with those of the Hawaiian bobtail squid, Euprymna scolopes. The genomes of the soft-bodied (coleoid) cephalopods are highly rearranged relative to other extant molluscs, indicating an intense, early burst of genome restructuring. The coleoid genomes feature multi-megabase, tandem arrays of genes associated with brain development and cephalopod-specific innovations. We find that a known coleoid hallmark, extensive A-to-I mRNA editing, displays two fundamentally distinct patterns: one exclusive to the nervous system and concentrated in genic sequences, the other widespread and directed toward repetitive elements. We conclude that coleoid novelty is mediated in part by substantial genome reorganization, gene family expansion, and tissue-dependent mRNA editing.

“Cephalopods are known for their large nervous systems, complex behaviors, and morphological innovations. Here, the authors find that soft-bodied cephalopod genomes are more rearranged than other extant molluscs and that mRNA editing patterns are associated with the nervous system and repetitive elements”.

Comparative genomics reveals insight into the evolutionary origin of massively scrambled genomes

Ciliates are microbial eukaryotes that undergo extensive programmed genome rearrangement, a natural genome editing process that converts long germline chromosomes into smaller gene-rich somatic chromosomes. Three well-studied ciliates include Oxytricha trifallax, Tetrahymena thermophila, and Paramecium tetraurelia, but only the Oxytricha lineage has a massively scrambled genome, whose assembly during development requires hundreds of thousands of precisely programmed DNA joining events, representing the most complex genome dynamics of any known organism. Here we study the emergence of such complex genomes by examining the origin and evolution of discontinuous and scrambled genes in the Oxytricha lineage. This study compares six genomes from three species, the germline and somatic genomes for Euplotes woodruffi, Tetmemena sp., and the model ciliate O. trifallax. We sequenced, assembled, and annotated the germline and somatic genomes of E. woodruffi, which provides an outgroup, and the germline genome of Tetmemena sp. We find that the germline genome of Tetmemena is as massively scrambled and interrupted as Oxytricha's: 13.6% of its gene loci require programmed translocations and/or inversions, with some genes requiring hundreds of precise gene editing events during development. This study revealed that the earlier diverged spirotrich, E. woodruffi, also has a scrambled genome, but only roughly half as many loci (7.3%) are scrambled. Furthermore, its scrambled genes are less complex, together supporting the position of Euplotes as a possible evolutionary intermediate in this lineage, in the process of accumulating complex evolutionary genome rearrangements, all of which require extensive repair to assemble functional coding regions. Comparative analysis also reveals that scrambled loci are often associated with local duplications, supporting a gradual model for the origin of complex, scrambled genomes via many small events of DNA duplication and decay.