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

Chromatin loops are an ancestral hallmark of the animal regulatory genome

In bilaterian animals, gene regulation is shaped by a combination of linear and spatial regulatory information. Regulatory elements along the genome are integrated into gene regulatory landscapes through chromatin compartmentalization1,2, insulation of neighbouring genomic regions3,4 and chromatin looping that brings together distal cis-regulatory sequences5. However, the evolution of these regulatory features is unknown because the three-dimensional genome architecture of most animal lineages remains unexplored6,7. To trace the evolutionary origins of animal genome regulation, here we characterized the physical organization of the genome in non-bilaterian animals (sponges, ctenophores, placozoans and cnidarians)8,9 and their closest unicellular relatives (ichthyosporeans, filastereans and choanoflagellates)10 by combining high-resolution chromosome conformation capture11,12 with epigenomic marks and gene expression data. Our comparative analysis showed that chromatin looping is a conserved feature of genome architecture in ctenophores, placozoans and cnidarians. These sequence-determined distal contacts involve both promoter-enhancer and promoter-promoter interactions. By contrast, chromatin loops are absent in the unicellular relatives of animals. Our findings indicate that spatial genome regulation emerged early in animal evolution. This evolutionary innovation introduced regulatory complexity, ultimately facilitating the diversification of animal developmental programmes and cell type repertoires.

The application of irreversible genomic states to define and trace ancient cell type homologies

Homology, or relationship among characters by common descent, has been notoriously difficult to assess for many morphological features, and cell types in particular. The ontogenetic origin of morphological traits means that the only physically inherited information is encoded in the genomes. However, the complexity of the underlying gene regulatory network and often miniscule changes that can impact gene expression, make it practically impossible to postulate a clear demarcation line for what molecular signature should "define" a homologous cell type between two deeply branching animals. In this Hypothesis article, we propose the use of the recently characterized irreversible genomic states, that occur after chromosomal and sub-chromosomal mixing of genes and regulatory elements, to dissect regulatory signatures of each cell type into irreversible and reversible configurations. While many of such states will be non-functional, some may permanently impact gene expression in a given cell type. Our proposal is that such evolutionarily irreversible, and thus synapomorphic, functional genomic states can constitute a criterion for the timing of the origin of deep evolutionary cell type homologies. Our proposal thus aims to close the gap between the clearly defined homology of the individual genomic characters and their genomic states to the homology at the phenotypic level through the identification of the underlying evolutionarily irreversible and regulatory linked states.

Cephalopod sex determination and its ancient evolutionary origin

Octopuses, squids, and cuttlefishes-the coleoid cephalopods-are a remarkable branch in the tree of life whose members exhibit a repertoire of sophisticated behaviors.1 As a clade, coleoids harbor an incredible variety of novel traits, including the most complex nervous system among invertebrates, derived camera-type eyes, and rapid adaptive camouflage abilities.2,3 The burst of evolutionary novelty that distinguishes cephalopods is even more striking in a phylogenetic context; cephalopods are a deeply diverged lineage that last shared a common ancestor with other extant molluscs in the Cambrian period, roughly 550 million years ago.4,5 With recent advances in genome sequencing technologies, we have the capability to explore the genomic foundations of cephalopod novelties. Here, using PacBio long-read sequencing of genomic DNA and Iso-Seq full-length mRNA sequencing, we provide a novel chromosome-scale reference genome and annotation for a female California two-spot octopus (O. bimaculoides). Our assembly reveals that the female octopus has just one sex chromosome, consistent with a ZO karyotype, whereas the male has two (ZZ), providing the first evidence of genetic sex determination in cephalopods. We use our assembly and annotation in combination with existing genomic information from other cephalopods to create the first whole-genome alignments from this group and demonstrate that the sex chromosome is of an ancient origin, before the radiation of extant cephalopods approximately 480 million years ago,4 and has been conserved to the present day in all cephalopod genomes available.

MolluscDB 2.0: a comprehensive functional and evolutionary genomics database for over 1400 molluscan species

Mollusca represents the second-largest animal phylum but remains less explored genomically. The increase in high-quality genomes and diverse functional genomic data holds great promise for advancing our understanding of molluscan biology and evolution. To address the opportunities and challenges facing the molluscan research community in managing vast multi-omics resources, we developed MolluscDB 2.0 (http://mgbase.qnlm.ac), which integrates extensive functional genomic data and offers user-friendly tools for multilevel integrative and comparative analyses. MolluscDB 2.0 covers 1450 species across all eight molluscan classes and compiles ∼4200 datasets, making it the most comprehensive multi-omics resource for molluscs to date. MolluscDB 2.0 expands the layers of multi-omics data, including genomes, bulk transcriptomes, single-cell transcriptomes, proteomes, epigenomes and metagenomes. MolluscDB 2.0 also more than doubles the number of functional modules and analytical tools, updating 14 original modules and introducing 20 new, specialized modules. Overall, MolluscDB 2.0 provides highly valuable, open-access multi-omics platform for the molluscan research community, expediting scientific discoveries and deepening our understanding of molluscan biology and evolution.

Evolution and function of chromatin domains across the tree of life

The genome of all organisms is spatially organized to function efficiently. The advent of genome-wide chromatin conformation capture (Hi-C) methods has revolutionized our ability to probe the three-dimensional (3D) organization of genomes across diverse species. In this Review, we compare 3D chromatin folding from bacteria and archaea to that in mammals and plants, focusing on topology at the level of gene regulatory domains. In doing so, we consider systematic similarities and differences that hint at the origin and evolution of spatial chromatin folding and its relation to gene activity. We discuss the universality of spatial chromatin domains in all kingdoms, each encompassing one to several genes. We also highlight differences between organisms and suggest that similar features in Hi-C matrices do not necessarily reflect the same biological process or function. Furthermore, we discuss the evolution of domain boundaries and boundary-forming proteins, which indicates that structural maintenance of chromosome (SMC) proteins and the transcription machinery are the ancestral sculptors of the genome. Architectural proteins such as CTCF serve as clade-specific determinants of genome organization. Finally, studies in many non-model organisms show that, despite the ancient origin of 3D chromatin folding and its intricate link to gene activity, evolution tolerates substantial changes in genome organization.

Cephalopod Sex Determination and its Ancient Evolutionary Origin

Octopuses, squids, and cuttlefishes - the coleoid cephalopods - are a remarkable branch in the tree of life whose members exhibit a repertoire of sophisticated behaviors (Hanlon and Messenger, 2018). As a clade, coleoids harbor an incredible variety of novel traits including the most complex nervous system among invertebrates, derived camera-type eyes, and rapid adaptive camouflage abilities (Young, 1971; Hanlon, 2007). The burst of evolutionary novelty that distinguishes cephalopods is even more striking in a phylogenetic context; cephalopods are a deeply diverged lineage that last share a common ancestor with other extant molluscs in the Cambrian period, roughly 550 million years ago (Ponder and Lindberg, 2008; Huang et al., 2022). With recent advances in genome sequencing technologies, we have the capability to explore the genomic foundations of cephalopod novelties. Here, using PacBio long-read sequencing of genomic DNA and IsoSeq full-length mRNA sequencing, we provide a novel chromosome-scale reference genome and annotation for a female California two-spot octopus (O. bimaculoides). Our assembly reveals that the female octopus has just one sex chromosome, consistent with a ZO karyotype, while the male has two (ZZ), providing the first evidence of genetic sex determination in cephalopods. We use our assembly and annotation in combination with existing genomic information from other cephalopods to create the first whole genome alignments from this group and demonstrate that the sex chromosome is of an ancient origin, before the radiation of extant cephalopods approximately 480 million years ago (Huang et al., 2022), and has been conserved to the present day in all cephalopod genomes available.

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.

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.

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.

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.

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.