Chromosome‐level genome assembly of the East Asian common octopus ( Octopus sinensis ) using PacBio sequencing and Hi‐C technology

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

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

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

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

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 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.

Better late than never: Optimising the proteomic analysis of field-collected octopus

Proteomics, the temporal study of proteins expressed by an organism, is a powerful technique that can reveal how organisms respond to biological perturbations, such as disease and environmental stress. Yet, the use of proteomics for addressing ecological questions has been limited, partly due to inadequate protocols for the sampling and preparation of animal tissues from the field. Although RNAlater is an ideal alternative to freezing for tissue preservation in transcriptomics studies, its suitability for the field could be more broadly examined. Moreover, existing protocols require samples to be preserved immediately to maintain protein integrity, yet the effects of delays in preservation on proteomic analyses have not been thoroughly tested. Hence, we optimised a proteomic workflow for wild-caught samples. First, we conducted a preliminary in-lab test using SDS-PAGE analysis on aquaria-reared Octopus berrima confirming that RNAlater can effectively preserve proteins up to 6 h after incubation, supporting its use in the field. Subsequently, we collected arm tips from wild-caught Octopus berrima and preserved them in homemade RNAlater immediately, 3 h, and 6 h after euthanasia. Processed tissue samples were analysed by liquid chromatography tandem mass spectrometry to ascertain protein differences between time delay in tissue preservation, as well as the influence of sex, tissue type, and tissue homogenisation methods. Over 3500 proteins were identified from all tissues, with bioinformatic analysis revealing protein abundances were largely consistent regardless of sample treatment. However, nearly 10% additional proteins were detected from tissues homogenised with metal beads compared to liquid nitrogen methods, indicating the beads were more efficient at extracting proteins. Our optimised workflow demonstrates that sampling non-model organisms from remote field sites is achievable and can facilitate extensive proteomic coverage without compromising protein integrity.

Genetic confirmation of Octopus insularis (Leite and Haimovici, 2008) in South Florida, United States using physical features and de novo genome assembly

The distribution of octopuses within the Octopus vulgaris species complex remains inadequately understood. Species determination can be complex and involves characterizing a specimen's physical features and comparing its genetic makeup to other populations. In this study, we present the first genetic confirmation of Octopus insularis (Leite and Haimovici, 2008) inhabiting the coastal waters of the Florida Keys, United States. We employed visual observations to identify species-specific body patterns of three wild-caught octopuses and used de novo genome assembly to confirm their species. All three specimens exhibited a red/white reticulated pattern on their ventral arm surface. Two specimens displayed body pattern components of deimatic display (white eye encircled by a light ring, with darkening around the eye). All visual observations were consistent with distinguishing features of O. insularis. We then compared mitochondrial subunits COI, COIII, and 16S in these specimens across all available annotated octopod sequences, including Sepia apama (Hotaling et al., 2021) as a control outgroup taxon. For species exhibiting intraspecific genomic variation, we included multiple sequences from geographically distinct populations. Laboratory specimens consistently clustered into a single taxonomic node with O. insularis. These findings confirm O. insularis presence in South Florida and suggest a more extensive northern distribution than previously assumed. Whole genome Illumina sequencing of multiple specimens enabled taxonomic identification with well-established DNA barcodes while also generating the first de novo full assembly of O. insularis. Furthermore, constructing and comparing phylogenetic trees for multiple conserved genes is essential for confirming the presence and delineation of cryptic species in the Caribbean.

CaaX-less lamins: Lophotrochozoa provide a glance at the playground of evolution

Nuclear lamins are the main components of the nuclear lamina in many eukaryotes. They are members of the intermediate filament (IF) protein family. Lamins differ from cytoplasmic IF proteins by the presence of a nuclear localisation sequence (NLS) and a C-terminal tetrapeptide, the CaaX motif. The CaaX motif is target of post-translational modifications including isoprenylation, proteolytic processing, and carboxyl-methylation. These modifications, in conjunction with the NLS, direct lamins to the inner nuclear membrane where they assemble into filaments. Lamins lacking a CaaX motif are unable to associate independently with nuclear membranes and remain in the nucleoplasm. So far, three species have been reported to exclusively express CaaX-less lamins. All three belong to the lophotrochozoan lineage. To find out whether they represent rare exceptions, we analysed lamins of representatives of 17 lophotrochozoan phyla. Here we report that all four clades of Rotifera as well as individual taxa of Mollusca and Annelida lack CaaX-lamins, but express lamins with alternative C-termini. Of note, the respective mollusc and annelid groups occupy very different phylogenetic ranks. Most of these alternative C-termini are rich in aromatic residues. A possible function of these residues in membrane association is discussed. Alternative splicing of terebellid lamin transcripts gives rise to two lamin variants, one with a CaaX motif and one with an alternative C-terminus. A similar situation is found in Arenicolidae, Opheliidae, Capitellidae, and Echiura. This points a way, how the switch from lamins carrying a CaaX motif to lamins with alternative C-termini may have occurred.

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.

High-Density Genetic Linkage Map of the Southern Blue-ringed Octopus (Octopodidae: Hapalochlaena maculosa)

Genetic linkage maps provide a useful resource for non-model genomes and can aid in genome reassembly to form more contiguous pseudo-chromosomes. We present the first linkage map of any cephalopod, H. maculosa, composed of 47 linkage groups (LG). A total of 2166 single nucleotide polymorphisms and 2455 presence–absence variant loci were utilised by Lep-Map3 in linkage map construction. The map length spans 2016.62 cM with an average marker distance of 0.85 cM. Integration of the recent H. maculosa genome allowed 1151 scaffolds comprising 34% of the total genomic sequence to be orientated and/or placed using 1278 markers across all 47 LG. The linkage map generated provides a new perspective on HOX gene distribution in octopods. In the H. maculosa linkage map three (SCR, LOX4 and POST1) of six identified HOX genes (HOX1/LAB, SCR, LOX2, LOX4, LOX5, POST1) were located within the same LG (LG 9). The generation of a linkage map for H. maculosa has provided a valuable resource for understanding the evolution of cephalopod genomes and will provide a base for future work.

A high-quality, haplotype-phased genome reconstruction reveals unexpected haplotype diversity in a pearl oyster

Homologous chromosomes in the diploid genome are thought to contain equivalent genetic information, but this common concept has not been fully verified in animal genomes with high heterozygosity. Here we report a near-complete, haplotype-phased, genome assembly of the pearl oyster, Pinctada fucata, using hi-fidelity (HiFi) long reads and chromosome conformation capture data. This assembly includes 14 pairs of long scaffolds (>38 Mb) corresponding to chromosomes (2n = 28). The accuracy of the assembly, as measured by an analysis of k-mers, is estimated to be 99.99997%. Moreover, the haplotypes contain 95.2% and 95.9%, respectively, complete and single-copy BUSCO genes, demonstrating the high quality of the assembly. Transposons comprise 53.3% of the assembly and are a major contributor to structural variations. Despite overall collinearity between haplotypes, one of the chromosomal scaffolds contains megabase-scale non-syntenic regions, which necessarily have never been detected and resolved in conventional haplotype-merged assemblies. These regions encode expanded gene families of NACHT, DZIP3/hRUL138-like HEPN, and immunoglobulin domains, multiplying the immunity gene repertoire, which we hypothesize is important for the innate immune capability of pearl oysters. The pearl oyster genome provides insight into remarkable haplotype diversity in animals.

Cell type diversity in a developing octopus brain

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

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.

Gene Recruitments and Dismissals in the Argonaut Genome Provide Insights into Pelagic Lifestyle Adaptation and Shell-like Eggcase Reacquisition

The paper nautilus or greater argonaut, Argonauta argo, is a species of octopods which is characterized by its pelagic lifestyle and by the presence of a protective spiral-shaped shell-like eggcase in females. To reveal the genomic background of how the species adapted to the pelagic lifestyle and acquired its shell-like eggcase, we sequenced the draft genome of the species. The genome size was 1.1 Gb, which is the smallest among the cephalopods known to date, with the top 215 scaffolds (average length 5,064,479 bp) covering 81% (1.09 Gb) of the total assembly. A total of 26,433 protein-coding genes were predicted from 16,802 assembled scaffolds. From these, we identified nearly intact HOX, Parahox, Wnt clusters, and some gene clusters that could probably be related to the pelagic lifestyle, such as reflectin, tyrosinase, and opsin. The gene models also revealed several homologous genes related to calcified shell formation in Conchiferan mollusks, such as Pif-like, SOD, and TRX. Interestingly, comparative genomics analysis revealed that the homologous genes for such genes were also found in the genome of the shell-less octopus, as well as Nautilus, which has a true outer shell. Therefore, the draft genome sequence of Arg. argo presented here has helped us to gain further insights into the genetic background of the dynamic recruitment and dismissal of genes to form an important, converging extended phenotypic structure such as the shell and the shell-like eggcase. Additionally, it allows us to explore the evolution of from benthic to pelagic lifestyles in cephalopods and octopods.

Genomes of Two Flying Squid Species Provide Novel Sights into Adaptations of Cephalopods to Pelagic Life

Pelagic cephalopods have evolved a series of fascinating traits, such as excellent visual acuity, high-speed agility, and photophores for adaptation to open pelagic oceans. However, the genetic mechanisms underpinning these traits are not well understood. Thus, in this study, we obtained high-quality genomes of two purpleback flying squid species (Sthenoteuthis oualaniensis and Sthenoteuthis sp.), with sizes of 5450 Mb and 5651 Mb, respectively. Comparative genomic analyses revealed that the S-crystallin subfamily SL20-1 associated with visual acuity in the purpleback flying squid lineage was significantly expanded, and the evolution of high-speed agility for the species was accompanied by significant positive selection pressure on genes related to energy metabolism. These molecular signals might have contributed to the evolution of their adaptative predatory and anti-predatory traits. In addition, the transcriptomic analysis provided clear indications of the evolution of the photophores of purpleback flying squids, especially the recruitment of new genes and energy metabolism-related genes which may have played key functional roles in the process.

Optimization of Whole Mount RNA Multiplexed in situ Hybridization Chain Reaction With Immunohistochemistry, Clearing and Imaging to Visualize Octopus Embryonic Neurogenesis

Gene expression analysis has been instrumental to understand the function of key factors during embryonic development of many species. Marker analysis is also used as a tool to investigate organ functioning and disease progression. As these processes happen in three dimensions, the development of technologies that enable detection of gene expression in the whole organ or embryo is essential. Here, we describe an optimized protocol of whole mount multiplexed RNA in situ hybridization chain reaction version 3.0 (HCR v3.0) in combination with immunohistochemistry (IHC), followed by fructose-glycerol clearing and light sheet fluorescence microscopy (LSFM) imaging on Octopus vulgaris embryos. We developed a code to automate probe design which can be applied for designing HCR v3.0 type probe pairs for fluorescent in situ mRNA visualization. As proof of concept, neuronal (Ov-elav) and glial (Ov-apolpp) markers were used for multiplexed HCR v3.0. Neural progenitor (Ov-ascl1) and precursor (Ov-neuroD) markers were combined with immunostaining for phosphorylated-histone H3, a marker for mitosis. After comparing several tissue clearing methods, fructose-glycerol clearing was found optimal in preserving the fluorescent signal of HCR v3.0. The expression that was observed in whole mount octopus embryos matched with the previous expression data gathered from paraffin-embedded transverse sections. Three-dimensional reconstruction revealed additional spatial organization that had not been discovered using two-dimensional methods.

Whole-Genome Sequencing and Genome-Wide Studies of Spiny Head Croaker (Collichthys lucidus) Reveals Potential Insights for Well-Developed Otoliths in the Family Sciaenidae

Spiny head croaker (Collichthys lucidus), belonging to the family Sciaenidae, is a small economic fish with a main distribution in the coastal waters of Northwestern Pacific. Here, we constructed a nonredundant chromosome-level genome assembly of spiny head croaker and also made genome-wide investigations on genome evolution and gene families related to otolith development. A primary genome assembly of 811.23 Mb, with a contig N50 of 74.92 kb, was generated by a combination of 49.12-Gb Illumina clean reads and 35.24 Gb of PacBio long reads. Contigs of this draft assembly were further anchored into chromosomes by integration with additional 185.33-Gb Hi-C data, resulting in a high-quality chromosome-level genome assembly of 817.24 Mb, with an improved scaffold N50 of 26.58 Mb. Based on our phylogenetic analysis, we observed that C. lucidus is much closer to Larimichthys crocea than Miichthys miiuy. We also predicted that many gene families were significantly expanded (p-value <0.05) in spiny head croaker; among them, some are associated with "calcium signaling pathway" and potential "inner ear functions." In addition, we identified some otolith-related genes (such as otol1a that encodes Otolin-1a) with critical deletions or mutations, suggesting possible molecular mechanisms for well-developed otoliths in the family Sciaenidae.

Identification and Characterization of a Rhodopsin Kinase Gene in the Suckers of Octopus vulgaris: Looking around Using Arms?

Simple Summary

Octopus arms are a fascinating and evolutionarily unique sensory organ, with hundreds of motile suckers, each with thousands of sensory cells, lining eight highly flexible arms. Scientifically, there are many open questions regarding the sensory capabilities of the arms and specifically the highly innervated suckers. In our present work, we used a multidisciplinary approach to fully characterize the light-sensing molecule, Ov-GRK1, in the suckers, skin and retina of Octopus vulgaris. We sequenced the O. vulgaris GRK1 gene, defining a phylogenetic tree and performing a 3D structure model prediction. We found differences in the relative expression of mRNA in different sucker types at several locations along the arm, which might indicate a functional difference. Using labeling methods, we localized the expression to the highly sensitive sucker rim. Our findings indicate that octopus suckers, in specific areas of the arm, might have the ability for light sensing. We therefore suggest that suckers are tactile, chemical and light sensors.

Abstract

In their foraging behavior octopuses rely on arm search movements outside the visual field of the eyes. In these movements the environment is explored primarily by the suckers that line the entire length of the octopus arm. In this study, for the first time, we report the complete characterization of a light-sensing molecule, Ov-GRK1, in the suckers, skin and retina of Octopus vulgaris. We sequenced the O. vulgaris GRK1 gene, defining a phylogenetic tree and performing a 3D structure model prediction. Furthermore, we found differences in relative mRNA expression in different sucker types at several arm levels, and localized it through in situ hybridization. Our findings suggest that the suckers in octopus arms are much more multimodal than was previously shown, adding the potential for light sensing to the already known mechanical and chemical sensing abilities.

Leveraging Short-Read Sequencing to Explore the Genomics of Sepiolid Squid

Due to their large size (∼3-5 Gb) and high repetitive content, the study of cephalopod genomes has historically been problematic. However, with the recent sequencing of several cephalopod genomes, including the Hawaiian bobtail squid (Euprymna scolopes), whole-genome studies of these molluscs are now possible. Of particular interest are the sepiolid or bobtail squids, many of which develop photophores in which bioluminescent bacterial symbionts reside. The variable presence of the symbiosis throughout the family allows us to determine regions of the genome that are under selection in symbiotic lineages, potentially providing a mechanism for identifying genes instrumental in the evolution of these mutualistic associations. To this end, we have used high-throughput sequencing to generate sequence from five bobtail squid genomes, four of which maintain symbioses with luminescent bacteria (E. hyllebergi, E. albatrossae, E. scolopes and Rondeletiola minor), and one of which does not (Sepietta neglecta). When we performed K-mer based heterozygosity and genome size estimations, we found that the Euprymna genus has a higher predicted genome size than other bobtail squid (∼ 5 Gb as compared to ∼ 4 Gb) and lower genomic heterozygosity. When we analyzed the repetitive content of the genomes, we found that genomes in the genus Euprymna appear to have recently acquired a significant quantity of LINE elements that are not found in its sister genus Rondeletiola or the closely related Sepietta. Using Abyss-2.0 and then Chromosomer with the published E. scolopes genome as a reference, we generated E. hyllebergi and E. albatrossae genomes of 1.54-1.57 Gb in size, but containing over 78-81% of eukaryotic single-copy othologs. The data we have generated will enable future whole-genome comparisons between these species to determine gene and regulatory content that differs between symbiotic and non-symbiotic lineages, as well as genes associated with symbiosis that are under selection.

Genomic insights into the adaptation and evolution of the nautilus, an ancient but evolving "living fossil"

The nautilus, commonly known as a "living fossil," is endangered and may be at risk of extinction. The lack of genomic information hinders a thorough understanding of its biology and evolution, which can shed light on the conservation of this endangered species. Here, we report the first high-quality chromosome-level genome assembly of Nautilus pompilius. The assembled genome size comprised 785.15 Mb. Comparative genomic analyses indicated that transposable elements (TEs) and large-scale genome reorganizations may have driven lineage-specific evolution in the cephalopods. Remarkably, evolving conserved genes and recent TE insertion activities were identified in N. pompilius, and we speculate that these findings reflect the strong adaptability and long-term survival of the nautilus. We also identified gene families that are potentially responsible for specific adaptation and evolution events. Our study provides unprecedented insights into the specialized biology and evolution of N. pompilius, and the results serve as an important resource for future conservation genomics of the nautilus and closely related species.

African Arowana Genome Provides Insights on Ancient Teleost Evolution

Osteoglossiformes is a basal clade of teleost, evolving since the Jurassic period. The genomes of Osteoglossiformes species would shed light on the evolution and adaptation of teleost. Here, we established a chromosome-level genome of African arowana. Together with the genomes of pirarucu and Asian arowana, we found that they diverged at ∼106.1 million years ago (MYA) and ∼59.2 MYA, respectively, which are coincident with continental separation. Interestingly, we identified a dynamic genome evolution characterized by a fast evolutionary rate and a high pseudogenization rate in African arowana and pirarucu. Additionally, more transposable elements were found in Asian arowana which confer more gene duplications. Moreover, we found the contraction of olfactory receptor and the expansion of UGT in African arowana might be related to its transformation from carnivore to be omnivore. Taken together, we provided valuable genomic resource of Osteoglossidae and revealed the correlation of biogeography and teleost evolution.

•

An evolutionary model of Osteoglossidae along the continental drift is provided

•

A faster evolving rate of African arowana than Asian arowana is revealed

•

The gene duplications of Asian arowana are related to more class I TE insertions

•

A mechanism of African arowana’s feeding habits transition is proposed.