Endosymbiosis with chemosynthetic bacteria has enabled many deep-sea invertebrates to thrive at hydrothermal vents and cold seeps, but most previous studies on this mutualism have focused on the bacteria only. Vesicomyid clams dominate global deep-sea chemosynthesis-based ecosystems. They differ from most deep-sea symbiotic animals in passing their symbionts from parent to offspring, enabling intricate co-evolution between the host and the symbiont. Here, we sequenced the genomes of the clam Archivesica marissinica (Bivalvia: Vesicomyidae) and its bacterial symbiont to understand the genomic/metabolic integration behind this symbiosis. At 1.52 gigabases, the clam genome encodes 28 genes horizontally transferred from bacteria, a large number of pseudogenes and transposable elements whose massive expansion corresponded to the timing of the rise and subsequent divergence of symbiont-bearing vesicomyids. The genome exhibits gene family expansion in cellular processes that likely facilitate chemoautotrophy, including gas delivery to support energy and carbon production, metabolite exchange with the symbiont, and regulation of the bacteriocyte population. Contraction in cellulase genes is likely adaptive to the shift from phytoplankton-derived to bacteria-based food. It also shows contraction in bacterial recognition gene familie, indicative of suppressed immune response to the endosymbiont. The gammaproteobacterium endosymbiont has a reduced genome of 1.03 megabases but retains complete pathways for sulfur oxidation, carbon fixation, and biosynthesis of 20 common amino acids, indicating the host’s high dependence on the symbiont for nutrition. Overall, the host-symbiont genomes show not only tight metabolic complementarity, but also distinct signatures of co-evolution allowing the vesicomyids to thrive in chemosynthesis-based ecosystems.

No abstract

Article impact statement: IUCN Red List assessments robustly compare extinction risk from deep-sea mining at understudied hydrothermal vent habitats.

Choosing the optimum assembly approach is essential to achieving a high-quality genome assembly suitable for comparative and evolutionary genomic investigations. Significant recent progress in long-read sequencing technologies such as PacBio and Oxford Nanopore Technologies (ONT) has also brought about a large variety of assemblers. Although these have been extensively tested on model species such as Homo sapiens and Drosophila melanogaster , such benchmarking has not been done in Mollusca, which lacks widely adopted model species. Molluscan genomes are notoriously rich in repeats and are often highly heterozygous, making their assembly challenging. Here, we benchmarked 10 assemblers based on ONT raw reads from two published molluscan genomes of differing properties, the gastropod Chrysomallon squamiferum (356.6 Mb, 1.59% heterozygosity) and the bivalve Mytilus coruscus (1593 Mb, 1.94% heterozygosity). By optimizing the assembly pipeline, we greatly improved both genomes from previously published versions. Our results suggested that 40–50X of ONT reads are sufficient for high-quality genomes, with Flye being the recommended assembler for compact and less heterozygous genomes exemplified by C. squamiferum , while NextDenovo excelled for more repetitive and heterozygous molluscan genomes exemplified by M. coruscus . A phylogenomic analysis using the two updated genomes with 32 other published high-quality lophotrochozoan genomes resulted in maximum support across all nodes, and we show that improved genome quality also leads to more complete matrices for phylogenomic inferences. Our benchmarking will ensure efficiency in future assemblies for molluscs and perhaps also for other marine phyla with few genomes available. This article is part of the Theo Murphy meeting issue ‘Molluscan genomics: broad insights and future directions for a neglected phylum’.

The distribution of species among spatially isolated habitat patches supports regional biodiversity and stability, so understanding the underlying processes and structure is a key target of conservation. Although multivariate statistics can infer the connectivity processes driving species distribution, such as dispersal and habitat suitability, they rarely explore the structure. Methods from graph theory, applied to distribution data, give insights into both connectivity pathways and processes by intuitively formatting the data as a network of habitat patches. We apply these methods to empirical data from the hydrothermal vent habitats of the Northwest Pacific. Hydrothermal vents are “oases” of biological productivity and endemicity on the seafloor that are imminently threatened by anthropogenic disturbances with unknown consequences to biodiversity. Here, we describe the structure of species assemblage networks at hydrothermal vents, how local and regional parameters affect their structure, and the implications for conservation. Two complementary networks were formed from an extensive species assemblage dataset: a similarity network of vent site nodes linked by weighted edges based on their pairwise assemblage similarity and a bipartite network of species nodes linked to vent site nodes at which they are present. Using these networks, we assessed the role of individual vent sites in maintaining network connectivity and identified biogeographic sub‐regions. The three sub‐regions and two outlying sites are separated by their spatial arrangement and local environmental filters. Both networks detected vent sites that play a disproportionately important role in regional pathways, while the bipartite network also identified key vent sites maintaining the distinct species assemblages of their sub‐regions. These regional connectivity pathways provide insights into historical colonization routes, while sub‐regional connectivity pathways are of value when selecting sites for conservation and/or estimating the multivent impacts from proposed deep‐sea mining.