Genome and transcriptome mechanisms driving cephalopod evolution

Neuronal segmentation in cephalopod arms

The prehensile arms of the cephalopod are among these animals most remarkable features, but the neural circuitry governing arm and sucker movements remains largely unknown. We studied the neuronal organization of the adult axial nerve cord (ANC) of Octopus bimaculoides with molecular and cellular methods. The ANCs, which lie in the center of every arm, are the largest neuronal structures in the octopus, containing four times as many neurons as found in the central brain. In transverse cross section, the cell body layer (CBL) of the ANC wraps around its neuropil (NP) with little apparent segregation of sensory and motor neurons or nerve exits. Strikingly, when studied in longitudinal sections, the ANC is segmented. ANC neuronal cell bodies form columns separated by septa, with 15 segments overlying each pair of suckers. The segments underlie a modular organization to the ANC neuropil: neuronal cell bodies within each segment send the bulk of their processes directly into the adjoining neuropil, with some reaching the contralateral side. In addition, some nerve processes branch upon entering the NP, forming short-range projections to neighboring segments and mid-range projections to the ANC segments of adjoining suckers. The septa between the segments are employed as ANC nerve exits and as channels for ANC vasculature. Cellular analysis establishes that adjoining septa issue nerves with distinct fiber trajectories, which across two segments (or three septa) fully innervate the arm musculature. Sucker nerves also use the septa, setting up a nerve fiber "suckerotopy" in the sucker-side of the ANC. Comparative anatomy suggests a strong link between segmentation and flexible sucker-laden arms. In the squid Doryteuthis pealeii, the arms and the sucker-rich club of the tentacles have segments, but the sucker-poor stalk of the tentacles does not. The neural modules described here provide a new template for understanding the motor control of octopus soft tissues. In addition, this finding represents the first demonstration of nervous system segmentation in a mollusc.

Dynamic skin behaviors in cephalopods

The coleoid cephalopods (cuttlefish, octopus, and squid) are a group of soft-bodied mollusks that exhibit a wealth of complex behaviors, including dynamic camouflage, object mimicry, skin-based visual communication, and dynamic body patterns during sleep. Many of these behaviors are visually driven and engage the animals' color changing skin, a pixelated display that is directly controlled by neurons projecting from the brain. Thus, cephalopod skin provides a direct readout of neural activity in the brain. During camouflage, cephalopods recreate on their skin an approximation of what they see, providing a window into perceptual processes in the brain. Additionally, cephalopods communicate their internal state during social encounters using innate skin patterns, and create waves of pigmentation on their skin during periods of arousal. Thus, by leveraging the visual displays of cephalopods, we can gain insight into how the external world is represented in the brain and how this representation is transformed into a recapitulation of the world on the skin. Here, we describe the rich skin behaviors of the coleoid cephalopods, what is known about cephalopod neuroanatomy, and how advancements in gene editing, machine learning, optical imaging, and electrophysiological tools may provide an opportunity to explore the neural bases of these fascinating behaviors.

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.

Acceleration of genome rearrangement in clitellate annelids

Comparisons of multiple metazoan genomes have revealed the existence of ancestral linkage groups (ALGs), genomic scaffolds sharing sets of orthologous genes that have been inherited from ancestral animals for hundreds of millions of years (Simakov et al. 2022; Schultz et al. 2023) These ALGs have persisted across major animal taxa including Cnidaria, Deuterostomia, Ecdysozoa and Spiralia. Notwithstanding this general trend of chromosome-scale conservation, ALGs have been obliterated by extensive genome rearrangements in certain groups, most notably including Clitellata (oligochaetes and leeches), a group of easily overlooked invertebrates that is of tremendous ecological, agricultural and economic importance (Charles 2019; Barrett 2016). To further investigate these rearrangements, we have undertaken a comparison of 12 clitellate genomes (including four newly sequenced species) and 11 outgroup representatives. We show that these rearrangements began at the base of the Clitellata (rather than progressing gradually throughout polychaete annelids), that the inter-chromosomal rearrangements continue in several clitellate lineages and that these events have substantially shaped the evolution of the otherwise highly conserved Hox cluster.

The NCBI Comparative Genome Viewer (CGV) is an interactive visualization tool for the analysis of whole-genome eukaryotic alignments

We report a new visualization tool for analysis of whole-genome assembly-assembly alignments, the Comparative Genome Viewer (CGV) (https://ncbi.nlm.nih.gov/genome/cgv/). CGV visualizes pairwise same-species and cross-species alignments provided by National Center for Biotechnology Information (NCBI) using assembly alignment algorithms developed by us and others. Researchers can examine large structural differences spanning chromosomes, such as inversions or translocations. Users can also navigate to regions of interest, where they can detect and analyze smaller-scale deletions and rearrangements within specific chromosome or gene regions. RefSeq or user-provided gene annotation is displayed where available. CGV currently provides approximately 800 alignments from over 350 animal, plant, and fungal species. CGV and related NCBI viewers are undergoing active development to further meet needs of the research community in comparative genome visualization.

Responses to commentaries on "The gene: An appraisal"

The central conclusions of "The Gene: An Appraisal" are that genetic variance does not underpin biological evolution, and, therefore, that genes are not Mendel's units of inheritance. In this response, I will address the criticisms I have received via commentaries on that paper by defending the following statements: 1. Epistasis does not explain the power-law fitness profile of the Long-Term Evolution Experiment (LTEE). The data from the evolution of natural systems displays the power-law form ubiquitously. Epistasis plays no role in evolution. 2. The common characteristics of living things (natural systems) are described by the principle of least action in de Maupertuis's original form, which is synonymous with the 2nd law of thermodynamics and Newton's 2nd law of motion in its complete form, i.e., F = dp/dt. Organisms strive to achieve free energy balance with their environments. 3. Based on an appraisal of the scientific environment between 1880 and 1911, I conclude that Johannsen's genotype conception was perhaps, the only option available to him. 4. The power-law fitness profile of the LTEE falsifies Fisher's Genetical Theory of Natural Selection, Johannsen's genotype conception, and the idea that 'Darwinian evolution' is an exception to the generic thermodynamic process of evolution in natural systems. 5. The use of the technique of genome-wide association to identify the causes and the likelihoods of inherited common diseases and behavioural traits is a 'wild goose chase' because genes are not the units of inheritance.

Extreme genome scrambling in marine planktonic Oikopleura dioica cryptic species

Genome structural variations within species are rare. How selective constraints preserve gene order and chromosome structure is a central question in evolutionary biology that remains unsolved. Our sequencing of several genomes of the appendicularian tunicate Oikopleura dioica around the globe reveals extreme genome scrambling caused by thousands of chromosomal rearrangements, although showing no obvious morphological differences between these animals. The breakpoint accumulation rate is an order of magnitude higher than in ascidian tunicates, nematodes, Drosophila, or mammals. Chromosome arms and sex-specific regions appear to be the primary unit of macrosynteny conservation. At the microsyntenic level, scrambling did not preserve operon structures, suggesting an absence of selective pressure to maintain them. The uncoupling of the genome scrambling with morphological conservation in O. dioica suggests the presence of previously unnoticed cryptic species and provides a new biological system that challenges our previous vision of speciation in which similar animals always share similar genome structures.

Parallel gene size and isoform expansion of ancient neuronal genes

How nervous systems evolved is a central question in biology. A diversity of synaptic proteins is thought to play a central role in the formation of specific synapses leading to nervous system complexity. The largest animal genes, often spanning hundreds of thousands of base pairs, are known to be enriched for expression in neurons at synapses and are frequently mutated or misregulated in neurological disorders and diseases. Although many of these genes have been studied independently in the context of nervous system evolution and disease, general principles underlying their parallel evolution remain unknown. To investigate this, we directly compared orthologous gene sizes across eukaryotes. By comparing relative gene sizes within organisms, we identified a distinct class of large genes with origins predating the diversification of animals and, in many cases, the emergence of neurons as dedicated cell types. We traced this class of ancient large genes through evolution and found orthologs of the large synaptic genes potentially driving the immense complexity of metazoan nervous systems, including in humans and cephalopods. Moreover, we found that while these genes are evolving under strong purifying selection, as demonstrated by low dN/dS ratios, they have simultaneously grown larger and gained the most isoforms in animals. This work provides a new lens through which to view this distinctive class of large and multi-isoform genes and demonstrates how intrinsic genomic properties, such as gene length, can provide flexibility in molecular evolution and allow groups of genes and their host organisms to evolve toward complexity.

Evolutionary history influences the microbiomes of a female symbiotic reproductive organ in cephalopods

Many female squids and cuttlefishes have a symbiotic reproductive organ called the accessory nidamental gland (ANG) that hosts a bacterial consortium involved with egg defense against pathogens and fouling organisms. While the ANG is found in multiple cephalopod families, little is known about the global microbial diversity of these ANG bacterial symbionts. We used 16S rRNA gene community analysis to characterize the ANG microbiome from different cephalopod species and assess the relationship between host and symbiont phylogenies. The ANG microbiome of 11 species of cephalopods from four families (superorder: Decapodiformes) that span seven geographic locations was characterized. Bacteria of class Alphaproteobacteria, Gammaproteobacteria, and Flavobacteriia were found in all species, yet analysis of amplicon sequence variants by multiple distance metrics revealed a significant difference between ANG microbiomes of cephalopod families (weighted/unweighted UniFrac, Bray-Curtis, P = 0.001). Despite being collected from widely disparate geographic locations, members of the family Sepiolidae (bobtail squid) shared many bacterial taxa including (~50%) Opitutae (Verrucomicrobia) and Ruegeria (Alphaproteobacteria) species. Furthermore, we tested for phylosymbiosis and found a positive correlation between host phylogenetic distance and bacterial community dissimilarity (Mantel test r = 0.7). These data suggest that closely related sepiolids select for distinct symbionts from similar bacterial taxa. Overall, the ANGs of different cephalopod species harbor distinct microbiomes and thus offer a diverse symbiont community to explore antimicrobial activity and other functional roles in host fitness.IMPORTANCEMany aquatic organisms recruit microbial symbionts from the environment that provide a variety of functions, including defense from pathogens. Some female cephalopods (squids, bobtail squids, and cuttlefish) have a reproductive organ called the accessory nidamental gland (ANG) that contains a bacterial consortium that protects eggs from pathogens. Despite the wide distribution of these cephalopods, whether they share similar microbiomes is unknown. Here, we studied the microbial diversity of the ANG in 11 species of cephalopods distributed over a broad geographic range and representing 15-120 million years of host divergence. The ANG microbiomes shared some bacterial taxa, but each cephalopod species had unique symbiotic members. Additionally, analysis of host-symbiont phylogenies suggests that the evolutionary histories of the partners have been important in shaping the ANG microbiome. This study advances our knowledge of cephalopod-bacteria relationships and provides a foundation to explore defensive symbionts in other systems.

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.

A chromosome-level genome for the nudibranch gastropod Berghia stephanieae helps parse clade-specific gene expression in novel and conserved phenotypes

BACKGROUND

How novel phenotypes originate from conserved genes, processes, and tissues remains a major question in biology. Research that sets out to answer this question often focuses on the conserved genes and processes involved, an approach that explicitly excludes the impact of genetic elements that may be classified as clade-specific, even though many of these genes are known to be important for many novel, or clade-restricted, phenotypes. This is especially true for understudied phyla such as mollusks, where limited genomic and functional biology resources for members of this phylum have long hindered assessments of genetic homology and function. To address this gap, we constructed a chromosome-level genome for the gastropod Berghia stephanieae (Valdés, 2005) to investigate the expression of clade-specific genes across both novel and conserved tissue types in this species.

RESULTS

The final assembled and filtered Berghia genome is comparable to other high-quality mollusk genomes in terms of size (1.05 Gb) and number of predicted genes (24,960 genes) and is highly contiguous. The proportion of upregulated, clade-specific genes varied across tissues, but with no clear trend between the proportion of clade-specific genes and the novelty of the tissue. However, more complex tissue like the brain had the highest total number of upregulated, clade-specific genes, though the ratio of upregulated clade-specific genes to the total number of upregulated genes was low.

CONCLUSIONS

Our results, when combined with previous research on the impact of novel genes on phenotypic evolution, highlight the fact that the complexity of the novel tissue or behavior, the type of novelty, and the developmental timing of evolutionary modifications will all influence how novel and conserved genes interact to generate diversity.

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.

Fine mapping of the Cepaea nemoralis shell colour and mid-banded loci using a high-density linkage map

Molluscs are a highly speciose phylum that exhibits an astonishing array of colours and patterns, yet relatively little progress has been made in identifying the underlying genes that determine phenotypic variation. One prominent example is the land snail Cepaea nemoralis for which classical genetic studies have shown that around nine loci, several physically linked and inherited together as a 'supergene', control the shell colour and banding polymorphism. As a first step towards identifying the genes involved, we used whole-genome resequencing of individuals from a laboratory cross to construct a high-density linkage map, and then trait mapping to identify 95% confidence intervals for the chromosomal region that contains the supergene, specifically the colour locus (C), and the unlinked mid-banded locus (U). The linkage map is made up of 215,593 markers, ordered into 22 linkage groups, with one large group making up ~27% of the genome. The C locus was mapped to a ~1.3 cM region on linkage group 11, and the U locus was mapped to a ~0.7 cM region on linkage group 15. The linkage map will serve as an important resource for further evolutionary and population genomic studies of C. nemoralis and related species, as well as the identification of candidate genes within the supergene and for the mid-banding phenotype.

Unusual Evolution of Cephalopod Tryptophan Indole-Lyases

Tryptophan indole-lyase (TIL), a pyridoxal-5-phosphate-dependent enzyme, catalyzes the hydrolysis of L-tryptophan (L-Trp) to indole and ammonium pyruvate. TIL is widely distributed among bacteria and bacterial TILs consist of a D2-symmetric homotetramer. On the other hand, TIL genes are also present in several metazoans. Cephalopods have two TILs, TILα and TILβ, which are believed to be derived from a gene duplication that occurred before octopus and squid diverged. However, both TILα and TILβ individually contain disruptive amino acid substitutions for TIL activity, and neither was active when expressed alone. When TILα and TILβ were coexpressed, however, they formed a heterotetramer that exhibited low TIL activity. The loss of TIL activity of the heterotetramer following site-directed mutagenesis strongly suggests that the active heterotetramer contains the TILα/TILβ heterodimer. Metazoan TILs generally have lower kcat values for L-Trp than those of bacterial TILs, but such low TIL activity may be rather suitable for metazoan physiology, where L-Trp is in high demand. Therefore, reduced activity may have been a less likely target for purifying selection in the evolution of cephalopod TILs. Meanwhile, the unusual evolution of cephalopod TILs may indicate the difficulty of post-gene duplication evolution of enzymes with catalytic sites contributed by multiple subunits, such as TIL.

Interactive visualization of whole eukaryote genome alignments using NCBI's Comparative Genome Viewer (CGV)

We report a new visualization tool for analysis of whole genome assembly-assembly alignments, the Comparative Genome Viewer (CGV) (https://ncbi.nlm.nih.gov/genome/cgv/). CGV visualizes pairwise same-species and cross-species alignments provided by NCBI using assembly alignment algorithms developed by us and others. Researchers can examine the alignments between the two assemblies using two alternate views: a chromosome ideogram-based view or a 2D genome dotplot. Whole genome alignment views expose large structural differences spanning chromosomes, such as inversions or translocations. Users can also navigate to regions of interest, where they can detect and analyze smaller-scale deletions and rearrangements within specific chromosome or gene regions. RefSeq or user-provided gene annotation is displayed in the ideogram view where available. CGV currently provides approximately 700 alignments from over 300 animal, plant, and fungal species. CGV and related NCBI viewers are undergoing active development to further meet needs of the research community in comparative genome visualization.

Spiralian genomics and the evolution of animal genome architecture

Recent developments in sequencing technologies have greatly improved our knowledge of phylogenetic relationships and genomic architectures throughout the tree of life. Spiralia, a diverse clade within Protostomia, is essential for understanding the evolutionary history of parasitism, gene conversion, nervous systems and animal body plans. In this review, we focus on the current hypotheses of spiralian phylogeny and investigate the impact of long-read sequencing on the quality of genome assemblies. We examine chromosome-level assemblies to highlight key genomic features that have driven spiralian evolution, including karyotype, synteny and the Hox gene organization. In addition, we show how chromosome rearrangement has influenced spiralian genomic structures. Although spiralian genomes have undergone substantial changes, they exhibit both conserved and lineage-specific features. We recommend increasing sequencing efforts and expanding functional genomics research to deepen insights into spiralian biology.

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.

High-level RNA editing diversifies the coleoid cephalopod brain proteome

Coleoid cephalopods (octopus, squid and cuttlefish) have unusually complex nervous systems. The coleoid nervous system is also the only one currently known to recode the majority of expressed proteins through A-to-I RNA editing. The deamination of adenosine by adenosine deaminase acting on RNA (ADAR) enzymes produces inosine, which is interpreted as guanosine during translation. If this occurs in an open reading frame, which is the case for tens of thousands of editing sites in coleoids, it can recode the encoded protein. Here, we describe recent findings aimed at deciphering the mechanisms underlying high-level recoding and its adaptive potential. We describe the complement of ADAR enzymes in cephalopods, including a recently discovered novel domain in sqADAR1. We further summarize current evidence supporting an adaptive role of high-level RNA recoding in coleoids, and review recent studies showing that a large proportion of recoding sites is temperature-sensitive. Despite these new findings, the mechanisms governing the high level of RNA recoding in coleoid cephalopods remain poorly understood. Recent advances using genome editing in squid may provide useful tools to further study A-to-I RNA editing in these animals.

Neural control of cephalopod camouflage

In Die Another Day, James Bond receives an Aston Martin that can render itself invisible by dynamically reproducing the surroundings on the car's "polymer skin". In what is widely regarded as the worst Bond movie ever, the invisible car scene is cited as the moment the plot plunges into the truly absurd. But what if nature had actually invented such a technology, and did so hundreds of millions of years ago? The coleoid cephalopods - octopus, cuttlefish and squid - are living examples of dynamic camouflage. Their skin is covered with a high-resolution array of 'cellular pixels' (chromatophores) that are controlled by the brain. To disappear into their surroundings, cephalopods recreate an approximation of their environment on their skin by activating different combinations of colored chromatophores. However, unlike the fictional Bond car, whose surface is coated in tiny cameras to detect the environment, cephalopods don't see the world with their skin. Instead, the visual world is detected by the eyes, processed in the brain, and then used to activate motor commands that direct the skin's camouflage pattern. Thus, cephalopod skin patterns are an external manifestation of their internal perception of the world. How do cephalopods approximate the world with their skin? What can this teach us about how brains work? And which neurobiological tools will be needed to uncover the neural basis of camouflage?

Evolution of cephalopod nervous systems

Giant brains have independently evolved twice on this planet, in vertebrates and in cephalopods (Figure 1A). Thus, the brains and nervous systems of cephalopods provide an important counterpoint to vertebrates in the search for generalities of brain organization and function. Their mere existence disproves various hypotheses proposed to explain the evolution of the mind and the human brain, such as cognition and large brains evolved only in long-lived animals with complex social systems and parental care, none of which is true of cephalopods. Therefore, it is worthwhile to review what is known about the evolution of cephalopod nervous systems to consider how it informs our understanding of general principles of brain evolution.

Scaphopoda is the sister taxon to Bivalvia: Evidence of ancient incomplete lineage sorting

The almost simultaneous emergence of major animal phyla during the early Cambrian shaped modern animal biodiversity. Reconstructing evolutionary relationships among such closely spaced branches in the animal tree of life has proven to be a major challenge, hindering understanding of early animal evolution and the fossil record. This is particularly true in the species-rich and highly varied Mollusca where dramatic inconsistency among paleontological, morphological, and molecular evidence has led to a long-standing debate about the group's phylogeny and the nature of dozens of enigmatic fossil taxa. A critical step needed to overcome this issue is to supplement available genomic data, which is plentiful for well-studied lineages, with genomes from rare but key lineages, such as Scaphopoda. Here, by presenting chromosome-level genomes from both extant scaphopod orders and leveraging complete genomes spanning Mollusca, we provide strong support for Scaphopoda as the sister taxon of Bivalvia, revitalizing the morphology-based Diasoma hypothesis originally proposed 50 years ago. Our molecular clock analysis confidently dates the split between Bivalvia and Scaphopoda at ~520 Ma, prompting a reinterpretation of controversial laterally compressed Early Cambrian fossils, including Anabarella, Watsonella, and Mellopegma, as stem diasomes. Moreover, we show that incongruence in the phylogenetic placement of Scaphopoda in previous phylogenomic studies was due to ancient incomplete lineage sorting (ILS) that occurred during the rapid radiation of Conchifera. Our findings highlight the need to consider ILS as a potential source of error in deep phylogeny reconstruction, especially in the context of the unique nature of the Cambrian Explosion.

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.

Ancient origins of complex neuronal genes

How nervous systems evolved is a central question in biology. An increasing diversity of synaptic proteins is thought to play a central role in the formation of specific synapses leading to nervous system complexity. The largest animal genes, often spanning millions of base pairs, are known to be enriched for expression in neurons at synapses and are frequently mutated or misregulated in neurological disorders and diseases. While many of these genes have been studied independently in the context of nervous system evolution and disease, general principles underlying their parallel evolution remain unknown. To investigate this, we directly compared orthologous gene sizes across eukaryotes. By comparing relative gene sizes within organisms, we identified a distinct class of large genes with origins predating the diversification of animals and in many cases the emergence of dedicated neuronal cell types. We traced this class of ancient large genes through evolution and found orthologs of the large synaptic genes driving the immense complexity of metazoan nervous systems, including in humans and cephalopods. Moreover, we found that while these genes are evolving under strong purifying selection as demonstrated by low dN/dS scores, they have simultaneously grown larger and gained the most isoforms in animals. This work provides a new lens through which to view this distinctive class of large and multi-isoform genes and demonstrates how intrinsic genomic properties, such as gene length, can provide flexibility in molecular evolution and allow groups of genes and their host organisms to evolve toward complexity.

Chromosome-level genome assembly of the caenogastropod snail Rapana venosa

The carnivorous gastropod Rapana venosa (Valenciennes, 1846) is one of the most notorious ecological invaders worldwide. Here, we present the first high-quality chromosome-scale reference R. venosa genome obtained via PacBio sequencing, Illumina paired-end sequencing, and high-throughput chromosome conformation capture scaffolding. The assembled genome has a size of 2.30 Gb, with a scaffold N50 length of 64.63 Mb, and is anchored to 35 chromosomes. It contains 29,649 protein-coding genes, 77.22% of which were functionally annotated. Given its high heterozygosity (1.41%) and large proportion of repeat sequences (57.72%), it is one of the most complex genome assemblies. This chromosome-level genome assembly of R. venosa is an important resource for understanding molluscan evolutionary adaption and provides a genetic basis for its biological invasion control.

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.

Creation of an albino squid line by CRISPR-Cas9 and its application for in vivo functional imaging of neural activity

Cephalopods are remarkable among invertebrates for their cognitive abilities, adaptive camouflage, novel structures, and propensity for recoding proteins through RNA editing. Due to the lack of genetically tractable cephalopod models, however, the mechanisms underlying these innovations are poorly understood. Genome editing tools such as CRISPR-Cas9 allow targeted mutations in diverse species to better link genes and function. One emerging cephalopod model, Euprymna berryi, produces large numbers of embryos that can be easily cultured throughout their life cycle and has a sequenced genome. As proof of principle, we used CRISPR-Cas9 in E. berryi to target the gene for tryptophan 2,3 dioxygenase (TDO), an enzyme required for the formation of ommochromes, the pigments present in the eyes and chromatophores of cephalopods. CRISPR-Cas9 ribonucleoproteins targeting tdo were injected into early embryos and then cultured to adulthood. Unexpectedly, the injected specimens were pigmented, despite verification of indels at the targeted sites by sequencing in injected animals (G0s). A homozygote knockout line for TDO, bred through multiple generations, was also pigmented. Surprisingly, a gene encoding indoleamine 2,3, dioxygenase (IDO), an enzyme that catalyzes the same reaction as TDO in vertebrates, was also present in E. berryi. Double knockouts of both tdo and ido with CRISPR-Cas9 produced an albino phenotype. We demonstrate the utility of these albinos for in vivo imaging of Ca²⁺ signaling in the brain using two-photon microscopy. These data show the feasibility of making gene knockout cephalopod lines that can be used for live imaging of neural activity in these behaviorally sophisticated organisms.

RNA recoding in cephalopods tailors microtubule motor protein function

RNA editing is a widespread epigenetic process that can alter the amino acid sequence of proteins, termed "recoding." In cephalopods, most transcripts are recoded, and recoding is hypothesized to be an adaptive strategy to generate phenotypic plasticity. However, how animals use RNA recoding dynamically is largely unexplored. We investigated the function of cephalopod RNA recoding in the microtubule motor proteins kinesin and dynein. We found that squid rapidly employ RNA recoding in response to changes in ocean temperature, and kinesin variants generated in cold seawater displayed enhanced motile properties in single-molecule experiments conducted in the cold. We also identified tissue-specific recoded squid kinesin variants that displayed distinct motile properties. Finally, we showed that cephalopod recoding sites can guide the discovery of functional substitutions in non-cephalopod kinesin and dynein. Thus, RNA recoding is a dynamic mechanism that generates phenotypic plasticity in cephalopods and can inform the characterization of conserved non-cephalopod proteins.

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.

Squid express conserved ADAR orthologs that possess novel features

The coleoid cephalopods display unusually extensive mRNA recoding by adenosine deamination, yet the underlying mechanisms are not well understood. Because the adenosine deaminases that act on RNA (ADAR) enzymes catalyze this form of RNA editing, the structure and function of the cephalopod orthologs may provide clues. Recent genome sequencing projects have provided blueprints for the full complement of coleoid cephalopod ADARs. Previous results from our laboratory have shown that squid express an ADAR2 homolog, with two splice variants named sqADAR2a and sqADAR2b and that these messages are extensively edited. Based on octopus and squid genomes, transcriptomes, and cDNA cloning, we discovered that two additional ADAR homologs are expressed in coleoids. The first is orthologous to vertebrate ADAR1. Unlike other ADAR1s, however, it contains a novel N-terminal domain of 641 aa that is predicted to be disordered, contains 67 phosphorylation motifs, and has an amino acid composition that is unusually high in serines and basic amino acids. mRNAs encoding sqADAR1 are themselves extensively edited. A third ADAR-like enzyme, sqADAR/D-like, which is not orthologous to any of the vertebrate isoforms, is also present. Messages encoding sqADAR/D-like are not edited. Studies using recombinant sqADARs suggest that only sqADAR1 and sqADAR2 are active adenosine deaminases, both on perfect duplex dsRNA and on a squid potassium channel mRNA substrate known to be edited in vivo. sqADAR/D-like shows no activity on these substrates. Overall, these results reveal some unique features in sqADARs that may contribute to the high-level RNA recoding observed in cephalopods.

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.

Differentiation trajectories of the Hydra nervous system reveal transcriptional regulators of neuronal fate

The small freshwater cnidarian polyp Hydra vulgaris uses adult stem cells (interstitial stem cells) to continually replace neurons throughout its life. This feature, combined with the ability to image the entire nervous system (Badhiwala et al., 2021; Dupre & Yuste, 2017) and availability of gene knockdown techniques (Juliano, Reich, et al., 2014; Lohmann et al., 1999; Vogg et al., 2022), makes Hydra a tractable model for studying nervous system development and regeneration at the whole-organism level. In this study, we use single-cell RNA sequencing and trajectory inference to provide a comprehensive molecular description of the adult nervous system. This includes the most detailed transcriptional characterization of the adult Hydra nervous system to date. We identified eleven unique neuron subtypes together with the transcriptional changes that occur as the interstitial stem cells differentiate into each subtype. Towards the goal of building gene regulatory networks to describe Hydra neuron differentiation, we identified 48 transcription factors expressed specifically in the Hydra nervous system, including many that are conserved regulators of neurogenesis in bilaterians. We also performed ATAC-seq on sorted neurons to uncover previously unidentified putative regulatory regions near neuron-specific genes. Finally, we provide evidence to support the existence of transdifferentiation between mature neuron subtypes and we identify previously unknown transition states in these pathways. All together, we provide a comprehensive transcriptional description of an entire adult nervous system, including differentiation and transdifferentiation pathways, which provides a significant advance towards understanding mechanisms that underlie nervous system regeneration.

Sensory specializations drive octopus and squid behaviour

The evolution of new traits enables expansion into new ecological and behavioural niches. Nonetheless, demonstrated connections between divergence in protein structure, function and lineage-specific behaviours remain rare. Here we show that both octopus and squid use cephalopod-specific chemotactile receptors (CRs) to sense their respective marine environments, but structural adaptations in these receptors support the sensation of specific molecules suited to distinct physiological roles. We find that squid express ancient CRs that more closely resemble related nicotinic acetylcholine receptors, whereas octopuses exhibit a more recent expansion in CRs consistent with their elaborated 'taste by touch' sensory system. Using a combination of genetic profiling, physiology and behavioural analyses, we identify the founding member of squid CRs that detects soluble bitter molecules that are relevant in ambush predation. We present the cryo-electron microscopy structure of a squid CR and compare this with octopus CRs1 and nicotinic receptors2. These analyses demonstrate an evolutionary transition from an ancestral aromatic 'cage' that coordinates soluble neurotransmitters or tastants to a more recent octopus CR hydrophobic binding pocket that traps insoluble molecules to mediate contact-dependent chemosensation. Thus, our study provides a foundation for understanding how adaptation of protein structure drives the diversification of organismal traits and behaviour.

Structural basis of sensory receptor evolution in octopus

Chemotactile receptors (CRs) are a cephalopod-specific innovation that allow octopuses to explore the seafloor via 'taste by touch'1. CRs diverged from nicotinic acetylcholine receptors to mediate contact-dependent chemosensation of insoluble molecules that do not readily diffuse in marine environments. Here we exploit octopus CRs to probe the structural basis of sensory receptor evolution. We present the cryo-electron microscopy structure of an octopus CR and compare it with nicotinic receptors to determine features that enable environmental sensation versus neurotransmission. Evolutionary, structural and biophysical analyses show that the channel architecture involved in cation permeation and signal transduction is conserved. By contrast, the orthosteric ligand-binding site is subject to diversifying selection, thereby mediating the detection of new molecules. Serendipitous findings in the cryo-electron microscopy structure reveal that the octopus CR ligand-binding pocket is exceptionally hydrophobic, enabling sensation of greasy compounds versus the small polar molecules detected by canonical neurotransmitter receptors. These discoveries provide a structural framework for understanding connections between evolutionary adaptations at the atomic level and the emergence of new organismal behaviour.

On the origin and evolution of RNA editing in metazoans

Extensive adenosine-to-inosine (A-to-I) editing of nuclear-transcribed mRNAs is the hallmark of metazoan transcriptional regulation. Here, by profiling the RNA editomes of 22 species that cover major groups of Holozoa, we provide substantial evidence supporting A-to-I mRNA editing as a regulatory innovation originating in the last common ancestor of extant metazoans. This ancient biochemistry process is preserved in most extant metazoan phyla and primarily targets endogenous double-stranded RNA (dsRNA) formed by evolutionarily young repeats. We also find intermolecular pairing of sense-antisense transcripts as an important mechanism for forming dsRNA substrates for A-to-I editing in some but not all lineages. Likewise, recoding editing is rarely shared across lineages but preferentially targets genes involved in neural and cytoskeleton systems in bilaterians. We conclude that metazoan A-to-I editing might first emerge as a safeguard mechanism against repeat-derived dsRNA and was later co-opted into diverse biological processes due to its mutagenic nature.

Extensive Recoding of the Neural Proteome in Cephalopods by RNA Editing

The coleoid cephalopods have the largest brains, and display the most complex behaviors, of all invertebrates. The molecular and cellular mechanisms that underlie these remarkable advancements remain largely unexplored. Early molecular cloning studies of squid ion channel transcripts uncovered an unusually large number of A→I RNA editing sites that recoded codons. Further cloning of other neural transcripts showed a similar pattern. The advent of deep-sequencing technologies and the associated bioinformatics allowed the mapping of RNA editing events across the entire neural transcriptomes of various cephalopods. The results were remarkable: They contained orders of magnitude more recoding editing sites than any other taxon. Although RNA editing sites are abundant in most multicellular metazoans, they rarely recode. In cephalopods, the majority of neural transcripts are recoded. Recent studies have focused on whether these events are adaptive, as well as other noncanonical aspects of cephalopod RNA editing.

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.

A Mini-Review on Reflectins, from Biochemical Properties to Bio-Inspired Applications

Some cephalopods (squids, octopuses, and cuttlefishes) produce dynamic structural colors, for camouflage or communication. The key to this remarkable capability is one group of specialized cells called iridocytes, which contain aligned membrane-enclosed platelets of high-reflective reflectins and work as intracellular Bragg reflectors. These reflectins have unusual amino acid compositions and sequential properties, which endows them with functional characteristics: an extremely high reflective index among natural proteins and the ability to answer various environmental stimuli. Based on their unique material composition and responsive self-organization properties, the material community has developed an impressive array of reflectin- or iridocyte-inspired optical systems with distinct tunable reflectance according to a series of internal and external factors. More recently, scientists have made creative attempts to engineer mammalian cells to explore the function potentials of reflectin proteins as well as their working mechanism in the cellular environment. Progress in wide scientific areas (biophysics, genomics, gene editing, etc.) brings in new opportunities to better understand reflectins and new approaches to fully utilize them. The work introduced the composition features, biochemical properties, the latest developments, future considerations of reflectins, and their inspiration applications to give newcomers a comprehensive understanding and mutually exchanged knowledge from different communities (e.g., biology and material).

The role of post-transcriptional modifications during development

While the existence of post-transcriptional modifications of RNA nucleotides has been known for decades, in most RNA species the exact positions of these modifications and their physiological function have been elusive until recently. Technological advances, such as high-throughput next-generation sequencing (NGS) methods and nanopore-based mapping technologies, have made it possible to map the position of these modifications with single nucleotide accuracy, and genetic screens have uncovered the “writer”, “reader” and “eraser” proteins that help to install, interpret and remove such modifications, respectively. These discoveries led to intensive research programmes with the aim of uncovering the roles of these modifications during diverse biological processes. In this review, we assess novel discoveries related to the role of post-transcriptional modifications during animal development, highlighting how these discoveries can affect multiple aspects of development from fertilization to differentiation in many species.

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.

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.