De Novo Genome Assembly of Limpet Bathyacmaea lactea (Gastropoda: Pectinodontidae): The First Reference Genome of a Deep-Sea Gastropod Endemic to Cold Seeps

Photosymbiosis shaped animal genome architecture and gene evolution as revealed in giant clams

Symbioses are major drivers of organismal diversification and phenotypic innovation. However, how long-term symbioses shape whole genome evolution in metazoans is still underexplored. Here, we use a giant clam (Tridacna maxima) genome to demonstrate how symbiosis has left complex signatures in an animal's genome. Giant clams thrive in oligotrophic waters by forming a remarkable association with photosymbiotic dinoflagellate algae. Genome-based demographic inferences uncover a tight correlation between T. maxima global population change and major paleoclimate and habitat shifts, revealing how abiotic and biotic factors may dictate T. maxima microevolution. Comparative analyses reveal genomic features that may be symbiosis-driven, including expansion and contraction of immunity-related gene families and a large proportion of lineage-specific genes. Strikingly, about 70% of the genome is composed of repetitive elements, especially transposable elements, most likely resulting from a symbiosis-adapted immune system. This work greatly enhances our understanding of genomic drivers of symbiosis that underlie metazoan evolution and diversification.

Metals and metallothionein evolution in snails: a contribution to the concept of metal-specific functionality from an animal model group

This is a critical review of what we know so far about the evolution of metallothioneins (MTs) in Gastropoda (snails, whelks, limpets and slugs), an important class of molluscs with over 90,000 known species. Particular attention will be paid to the evolution of snail MTs in relation to the role of some metallic trace elements (cadmium, zinc and copper) and their interaction with MTs, also compared to MTs from other animal phyla. The article also highlights the important distinction, yet close relationship, between the structural and metal-selective binding properties of gastropod MTs and their physiological functionality in the living organism. It appears that in the course of the evolution of Gastropoda, the trace metal cadmium (Cd) must have played an essential role in the development of Cd-selective MT variants. It is shown how the structures and Cd-selective binding properties in the basal gastropod clades have evolved by testing and optimizing different combinations of ancestral and novel MT domains, and how some of these domains have become established in modern and recent gastropod clades. In this context, the question of how adaptation to new habitats and lifestyles has affected the original MT traits in different gastropod lineages will also be addressed. The 3D structures and their metal binding preferences will be highlighted exemplarily in MTs of modern littorinid and helicid snails. Finally, the importance of the different metal requirements and pathways in snail tissues and cells for the shaping and functionality of the respective MT isoforms will be shown.

Deciphering the population structure and genetic basis of growth traits from whole-genome resequencing of the leopard coral grouper ( Plectropomus leopardus)

The leopard coral grouper ( Plectropomus leopardus) is a species of significant economic importance. Although artificial cultivation of P. leopardus has thrived in recent decades, the advancement of selective breeding has been hindered by the lack of comprehensive population genomic data. In this study, we identified over 8.73 million single nucleotide polymorphisms (SNPs) through whole-genome resequencing of 326 individuals spanning six distinct groups. Furthermore, we categorized 226 individuals with high-coverage sequencing depth (≥14×) into eight clusters based on their genetic profiles and phylogenetic relationships. Notably, four of these clusters exhibited pronounced genetic differentiation compared with the other populations. To identify potentially advantageous loci for P. leopardus, we examined genomic regions exhibiting selective sweeps by analyzing the nucleotide diversity ( θπ) and fixation index ( F ST) in these four clusters. Using these high-coverage resequencing data, we successfully constructed the first haplotype reference panel specific to P. leopardus. This achievement holds promise for enabling high-quality, cost-effective imputation methods. Additionally, we combined low-coverage sequencing data with imputation techniques for a genome-wide association study, aiming to identify candidate SNP loci and genes associated with growth traits. A significant concentration of these genes was observed on chromosome 17, which is primarily involved in skeletal muscle and embryonic development and cell proliferation. Notably, our detailed investigation of growth-related SNPs across the eight clusters revealed that cluster 5 harbored the most promising candidate SNPs, showing potential for genetic selective breeding efforts. These findings provide a robust toolkit and valuable insights into the management of germplasm resources and genome-driven breeding initiatives targeting P. leopardus.

Insights into the genome of the 'Loco' Concholepas concholepas (Gastropoda: Muricidae) from low-coverage short-read sequencing: genome size, ploidy, transposable elements, nuclear RNA gene operon, mitochondrial genome, and phylogenetic placement in the family Muricidae

BACKGROUND

The Peruvian 'chanque' or Chilean 'loco' Concholepas concholepas is an economically, ecologically, and culturally important muricid gastropod heavily exploited by artisanal fisheries in the temperate southeastern Pacific Ocean. In this study, we have profited from a set of bioinformatics tools to recover important biological information of C. concholepas from low-coverage short-read NGS datasets. Specifically, we calculated the size of the nuclear genome, ploidy, and estimated transposable elements content using an in silico k-mer approach, we discovered, annotated, and quantified those transposable elements, we assembled and annotated the 45S rDNA RNA operon and mitochondrial genome, and we confirmed the phylogenetic position of C. concholepas within the muricid subfamily Rapaninae based on translated protein coding genes.

RESULTS

Using a k-mer approach, the haploid genome size estimated for the predicted diploid genome of C. concholepas varied between 1.83 Gbp (with kmer = 24) and 2.32 Gbp (with kmer = 36). Between half and two thirds of the nuclear genome of C. concholepas was composed of transposable elements. The most common transposable elements were classified as Long Interspersed Nuclear Elements and Short Interspersed Nuclear Elements, which were more abundant than DNA transposons, simple repeats, and Long Terminal Repeats. Less abundant repeat elements included Helitron mobile elements, 45S rRNA DNA, and Satellite DNA, among a few others.The 45S rRNA DNA operon of C. concholepas that encodes for the ssrRNA, 5.8S rRNA, and lsrRNA genes was assembled into a single contig 8,090 bp long. The assembled mitochondrial genome of C. concholepas is 15,449 bp long and encodes 13 protein coding genes, two ribosomal genes, and 22 transfer RNAs.

CONCLUSION

The information gained by this study will inform the assembly of a high quality nuclear genome for C. concholepas and will support bioprospecting and biomonitoring using environmental DNA to advance development of conservation and management plans in this overexploited marine snail.

The genome of the deep-sea anemone Actinernus sp. contains a mega-array of ANTP-class homeobox genes

Members of the phylum Cnidaria include sea anemones, corals and jellyfish, and have successfully colonized both marine and freshwater habitats throughout the world. The understanding of how cnidarians adapt to extreme environments such as the dark, high-pressure deep-sea habitat has been hindered by the lack of genomic information. Here, we report the first chromosome-level deep-sea cnidarian genome, of the anemone Actinernus sp., which was 1.39 Gbp in length and contained 44 970 gene models including 14 806 tRNA genes and 30 164 protein-coding genes. Analyses of homeobox genes revealed the longest chromosome hosts a mega-array of Hox cluster, HoxL, NK cluster and NKL homeobox genes; until now, such an array has only been hypothesized to have existed in ancient ancestral genomes. In addition to this striking arrangement of homeobox genes, analyses of microRNAs revealed cnidarian-specific complements that are distinctive for nested clades of these animals, presumably reflecting the progressive evolution of the gene regulatory networks in which they are embedded. Also, compared with other sea anemones, circadian rhythm genes were lost in Actinernus sp., which likely reflects adaptation to living in the dark. This high-quality genome of a deep-sea cnidarian thus reveals some of the likely molecular adaptations of this ecologically important group of metazoans to the extreme deep-sea environment. It also deepens our understanding of the evolution of genome content and organization of animals in general and cnidarians in particular, specifically from the viewpoint of key developmental control genes like the homeobox-encoding genes, where we find an array of genes that until now has only been hypothesized to have existed in the ancient ancestor that pre-dated both the cnidarians and bilaterians.

The First Genome Survey of the Snail Provanna glabra Inhabiting Deep-Sea Hydrothermal Vents

The snail P. glabra is an endemic species in deep-sea chemosynthetic ecosystems of the Northwest Pacific Ocean. To obtain more genetic information on this species and provide the basis for subsequent whole-genome map construction, a genome survey was performed on this snail from the hydrothermal vent of Okinawa Trough. The genomic size of P. glabra was estimated to be 1.44 Gb, with a heterozygosity of 1.91% and a repeated sequence content of 69.80%. Based on the sequencing data, a draft genome of 1.32 Gb was assembled. Transposal elements (TEs) accounted for 40.17% of the entire genome, with DNA transposons taking the highest proportion. It was found that most TEs were inserted in the genome recently. In the simple sequence repeats, the dinucleotide motif was the most enriched microsatellite type, accounting for 53% of microsatellites. A complete mitochondrial genome of P. glabra with a total length of 16,268 bp was assembled from the sequencing data. After comparison with the published mitochondrial genome of Provanna sp. from a methane seep, 331 potential single nucleotide polymorphism (SNP) sites were identified in protein-coding genes (PCGs). Except for the cox1 gene, nad2, nad4, nad5, and cob genes are expected to be candidate markers for population genetic and phylogenetic studies of P. glabra and other deep-sea snails. Compared with shallow-water species, three mitochondrial genes of deep-sea gastropods exhibited a higher evolutionary rate, indicating strong selection operating on mitochondria of deep-sea species. This study provides insights into the genome characteristics of P. glabra and supplies genomic resources for further studies on the adaptive evolution of the snail in extreme deep-sea chemosynthetic environments.

A chromosome-level genome assembly of a deep-sea starfish (Zoroaster cf. ophiactis)

Understanding of adaptation and evolution of organisms in the deep sea requires more genomic resources. Zoroaster cf. ophiactis is a sea star in the family Zoroasteridae occurring exclusively in the deep sea. In this study, a chromosome-level genome assembly for Z. cf. ophiactis was generated by combining Nanopore long-read, Illumina short-read, and Hi-C sequencing data. The final assembly was 1,002.0 Mb in length, with a contig N50 of 376 Kb and a scaffold N50 of 40.4 Mb, and included 22 pseudo-chromosomes, covering 92.3% of the assembly. Completeness analysis evaluated with BUSCO revealed that 95.91% of the metazoan conserved genes were complete. Additionally, 39,426 protein-coding genes were annotated for this assembly. This chromosome-level genome assembly represents the first high-quality genome for the deep-sea Asteroidea, and will provide a valuable resource for future studies on evolution and adaptation of deep-sea echinoderms.

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.

Transcriptome analysis of gene expression profiling from the deep sea in situ to the laboratory for the cold seep mussel Gigantidas haimaensis

Background

The deep-sea mussel Gigantidas haimaensis is a representative species from the Haima cold seep ecosystem in the South China Sea that establishes endosymbiosis with chemotrophic bacteria. During long-term evolution, G. haimaensis has adapted well to the local environment of cold seeps. Until now, adaptive mechanisms responding to environmental stresses have remained poorly understood.

Results

In this study, transcriptomic analysis was performed for muscle tissue of G. haimaensis in the in situ environment (MH) and laboratory environment for 0 h (M0), 3 h (M3) and 9 h (M9), and 187,368 transcript sequences and 22,924 annotated differentially expressed genes (DEGs) were generated. Based on Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, these DEGs were enriched with a broad spectrum of biological processes and pathways, including those associated with antioxidants, apoptosis, chaperones, immunity and metabolism. Among these significantly enriched pathways, protein processing in the endoplasmic reticulum and metabolism were the most affected metabolic pathways. These results may imply that G. haimaensis struggles to support the life response to environmental change by changing gene expression profiles.

Conclusion

The present study provides a better understanding of the biological responses and survival strategies of the mussel G. haimaensis from deep sea in situ to the laboratory environment.

Genome of a giant isopod, Bathynomus jamesi, provides insights into body size evolution and adaptation to deep-sea environment

BACKGROUND

The deep-sea may be regarded as a hostile living environment, due to low temperature, high hydrostatic pressure, and limited food and light. Isopods, a species-rich group of crustaceans, are widely distributed across different environments including the deep sea and as such are a useful model for studying adaptation, migration, and speciation. Similar to other deep-sea organisms, giant isopods have larger body size than their shallow water relatives and have large stomachs and fat bodies presumably to store organic reserves. In order to shed light on the genetic basis of these large crustaceans adapting to the oligotrophic environment of deep-sea, the high-quality genome of a deep-sea giant isopod Bathynomus jamesi was sequenced and assembled.

RESULTS

B. jamesi has a large genome of 5.89 Gb, representing the largest sequenced crustacean genome to date. Its large genome size is mainly attributable to the remarkable proliferation of transposable elements (84%), which may enable high genome plasticity for adaptive evolution. Unlike its relatives with small body size, B. jamesi has expanded gene families related to pathways of thyroid and insulin hormone signaling that potentially contribute to its large body size. Transcriptomic analysis showed that some expanded gene families related to glycolysis and vesicular transport were specifically expressed in its digestive organs. In addition, comparative genomics and gene expression analyses in six tissues suggested that B. jamesi has inefficient lipid degradation, low basal metabolic rate, and bulk food storage, suggesting giant isopods adopt a more efficient mechanism of nutrient absorption, storage, and utilization to provide sustained energy supply for their large body size.

CONCLUSIONS

Taken together, the giant isopod genome may provide a valuable resource for understanding body size evolution and adaptation mechanisms of macrobenthic organisms to deep-sea environments.

Comparative transcriptomic analysis of in situ and onboard fixed deep-sea limpets reveals sample preparation-related differences

Precise gene expression reflects the molecular response of deep-sea organisms to their harsh living environments. However, changes in environmental factors during lifting samples from the deep sea to a research vessel can also affect gene expression. By using the transcriptomic approach, we compared the gene expression profiles of the onboard fixed with the in situ fixed samples of the deep-sea limpet Bathyacmaea lactea. Our results revealed that the concomitant stress during conventional deep-sea sampling without RNA in situ fixation greatly influenced the gene expression. Various biological activities, such as cell and tissue structure, lysosomal activity, fluid balance, and unsaturated fatty acid metabolism, were perturbed, suggesting that the sampling stress has exerted systemic impacts on the life of the limpets. These findings clearly illustrate that deep-sea samples without RNA in situ fixation can easily lead to biased results in gene expression analysis, which requires to be appropriately addressed in future studies.

The genome of an apodid holothuroid (Chiridota heheva) provides insights into its adaptation to a deep-sea reducing environment

Cold seeps and hydrothermal vents are deep-sea reducing environments that are characterized by lacking oxygen and photosynthesis-derived nutrients. Most animals acquire nutrition in cold seeps or hydrothermal vents by maintaining epi- or endosymbiotic relationship with chemoautotrophic microorganisms. Although several seep- and vent-dwelling animals hosting symbiotic microbes have been well-studied, the genomic basis of adaptation to deep-sea reducing environment in nonsymbiotic animals is still lacking. Here, we report a high-quality genome of Chiridota heheva Pawson & Vance, 2004, which thrives by extracting organic components from sediment detritus and suspended material, as a reference for nonsymbiotic animal's adaptation to deep-sea reducing environments. The expansion of the aerolysin-like protein family in C. heheva compared with other echinoderms might be involved in the disintegration of microbes during digestion. Moreover, several hypoxia-related genes (Pyruvate Kinase M2, PKM2; Phospholysine Phosphohistidine Inorganic Pyrophosphate Phosphatase, LHPP; Poly(A)-specific Ribonuclease Subunit PAN2, PAN2; and Ribosomal RNA Processing 9, RRP9) were subject to positive selection in the genome of C. heheva, which contributes to their adaptation to hypoxic environments.

The First Draft Genome of a Cold-Water Coral Trachythela sp. (Alcyonacea: Stolonifera: Clavulariidae)

Cold-water corals (CWCs) are important habitats for creatures in the deep-sea environment, but they have been degraded by anthropogenic activity. So far, no genome for any CWC has been reported. Here, we report a draft genome of Trachythela sp., which represents the first genome of CWCs to date. In total, 56 and 65 Gb of raw reads were generated from Illumina and Nanopore sequencing platforms, respectively. The final assembled genome was 578.26 Mb, which consisted of 396 contigs with a contig N50 of 3.56 Mb, and the genome captured 90.1% of the metazoan Benchmarking Universal Single-Copy Orthologs. We identified 335 Mb (57.88% of the genome) of repetitive elements, which is a higher proportion compared with others in the Cnidarians, along with 35,305 protein-coding genes. We also detected 483 expanded and 51 contracted gene families, and many of them were associated with longevity, ion transposase, heme-binding nicotinamide adenine dinucleotide, and metabolic regulators of transcription. Overall, we believe this genome will serve as an important resource for studies on community protection for CWCs.

Convergent Evolution and Structural Adaptation to the Deep Ocean in the Protein-Folding Chaperonin CCTα

The deep ocean is the largest biome on Earth and yet it is among the least studied environments of our planet. Life at great depths requires several specific adaptations; however, their molecular mechanisms remain understudied. We examined patterns of positive selection in 416 genes from four brittle star (Ophiuroidea) families displaying replicated events of deep-sea colonization (288 individuals from 216 species). We found consistent signatures of molecular convergence in functions related to protein biogenesis, including protein folding and translation. Five genes were recurrently positively selected, including chaperonin-containing TCP-1 subunit α (CCTα), which is essential for protein folding. Molecular convergence was detected at the functional and gene levels but not at the amino-acid level. Pressure-adapted proteins are expected to display higher stability to counteract the effects of denaturation. We thus examined in silico local protein stability of CCTα across the ophiuroid tree of life (967 individuals from 725 species) in a phylogenetically corrected context and found that deep-sea-adapted proteins display higher stability within and next to the substrate-binding region, which was confirmed by in silico global protein stability analyses. This suggests that CCTα displays not only structural but also functional adaptations to deep-water conditions. The CCT complex is involved in the folding of ∼10% of newly synthesized proteins and has previously been categorized as a "cold-shock" protein in numerous eukaryotes. We thus propose that adaptation mechanisms to cold and deep-sea environments may be linked and highlight that efficient protein biogenesis, including protein folding and translation, is a key metabolic deep-sea adaptation.