Genomic adaptations to chemosymbiosis in the deep-sea seep-dwelling tubeworm Lamellibrachia luymesi

Chemosynthetic alphaproteobacterial diazotrophs reside in deep-sea cold-seep bottom waters

Nitrogen (N)-fixing organisms, also known as diazotrophs, play a crucial role in N-limited ecosystems by controlling the production of bioavailable N. The carbon-dominated cold-seep ecosystems are inherently N-limited, making them hotspots of N fixation. However, the knowledge of diazotrophs in cold-seep ecosystems is limited compared to other marine ecosystems. In this study, we used multi-omics to investigate the diversity and catabolism of diazotrophs in deep-sea cold-seep bottom waters. Our findings showed that the relative abundance of diazotrophs in the bacterial community reached its highest level in the cold-seep bottom waters compared to the cold-seep upper waters and non-seep bottom waters. Remarkably, more than 98% of metatranscriptomic reads aligned on diazotrophs in cold-seep bottom waters belonged to the genus Sagittula, an alphaproteobacterium. Its metagenome-assembled genome, named Seep-BW-D1, contained catalytic genes (nifHDK) for nitrogen fixation, and the nifH gene was actively transcribed in situ. Seep-BW-D1 also exhibited chemosynthetic capability to oxidize C1 compounds (methanol, formaldehyde, and formate) and thiosulfate (S2O32-). In addition, we observed abundant transcripts mapped to genes involved in the transport systems for acetate, spermidine/putrescine, and pectin oligomers, suggesting that Seep-BW-D1 can utilize organics from the intermediates synthesized by methane-oxidizing microorganisms, decaying tissues from cold-seep benthic animals, and refractory pectin derived from upper photosynthetic ecosystems. Overall, our study corroborates that carbon-dominated cold-seep bottom waters select for diazotrophs and reveals the catabolism of a novel chemosynthetic alphaproteobacterial diazotroph in cold-seep bottom waters.

IMPORTANCE

Bioavailable nitrogen (N) is a crucial element for cellular growth and division, and its production is controlled by diazotrophs. Marine diazotrophs contribute to nearly half of the global fixed N and perform N fixation in various marine ecosystems. While previous studies mainly focused on diazotrophs in the sunlit ocean and oxygen minimum zones, recent research has recognized cold-seep ecosystems as overlooked N-fixing hotspots because the seeping fluids in cold-seep ecosystems introduce abundant bioavailable carbon but little bioavailable N, making most cold seeps inherently N-limited. With thousands of cold-seep ecosystems detected at continental margins worldwide in the past decades, the significant role of cold seeps in marine N biogeochemical cycling is emphasized. However, the diazotrophs in cold-seep bottom waters remain poorly understood. Through multi-omics, this study identified a novel alphaproteobacterial chemoheterotroph belonging to Sagittula as one of the most active diazotrophs residing in cold-seep bottom waters and revealed its catabolism.

Deciphering deep-sea chemosynthetic symbiosis by single-nucleus RNA-sequencing

Bathymodioline mussels dominate deep-sea methane seep and hydrothermal vent habitats and obtain nutrients and energy primarily through chemosynthetic endosymbiotic bacteria in the bacteriocytes of their gill. However, the molecular mechanisms that orchestrate mussel host-symbiont interactions remain unclear. Here, we constructed a comprehensive cell atlas of the gill in the mussel Gigantidas platifrons from the South China Sea methane seeps (1100 m depth) using single-nucleus RNA-sequencing (snRNA-seq) and whole-mount in situ hybridisation. We identified 13 types of cells, including three previously unknown ones, and uncovered unknown tissue heterogeneity. Every cell type has a designated function in supporting the gill's structure and function, creating an optimal environment for chemosynthesis, and effectively acquiring nutrients from the endosymbiotic bacteria. Analysis of snRNA-seq of in situ transplanted mussels clearly showed the shifts in cell state in response to environmental oscillations. Our findings provide insight into the principles of host-symbiont interaction and the bivalves' environmental adaption mechanisms.

Single-cell RNA-seq reveals distinct metabolic "microniches" and close host-symbiont interactions in deep-sea chemosynthetic tubeworm

Vestimentiferan tubeworms that thrive in deep-sea chemosynthetic ecosystems rely on a single species of sulfide-oxidizing gammaproteobacterial endosymbionts housed in a specialized symbiotic organ called trophosome as their primary carbon source. While this simple symbiosis is remarkably productive, the host-symbiont molecular interactions remain unelucidated. Here, we applied an approach for deep-sea in situ single-cell fixation in a cold-seep tubeworm, Paraescarpia echinospica. Single-cell RNA sequencing analysis and further molecular characterizations of both the trophosome and endosymbiont indicate that the tubeworm maintains two distinct metabolic "microniches" in the trophosome by controlling the availability of chemosynthetic gases and metabolites, resulting in oxygenated and hypoxic conditions. The endosymbionts in the oxygenated niche actively conduct autotrophic carbon fixation and are digested for nutrients, while those in the hypoxic niche conduct anaerobic denitrification, which helps the host remove ammonia waste. Our study provides insights into the molecular interactions between animals and their symbiotic microbes.

Genetic adaptations of sea anemone to hydrothermal environment

Hydrothermal vent habitats are characterized by high hydrostatic pressure, darkness, and the continuous release of toxic metal ions into the surrounding environment where sea anemones and other invertebrates thrive. Nevertheless, the understanding of metazoan metal ion tolerances and environmental adaptations remains limited. We assembled a chromosome-level genome for the vent sea anemone, Alvinactis idsseensis sp. nov. Comparative genomic analyses revealed gene family expansions and gene innovations in A. idsseensis sp. nov. as a response to high concentrations of metal ions. Impressively, the metal tolerance proteins MTPs is a unique evolutionary response to the high concentrations of Fe2+ and Mn2+ present in the environments of these anemones. We also found genes associated with high concentrations of polyunsaturated fatty acids that may respond to high hydrostatic pressure and found sensory and circadian rhythm-regulated genes that were essential for adaptations to darkness. Overall, our results provide insights into metazoan adaptation to metal ions, high pressure, and darkness in hydrothermal vents.

Third-Generation Sequencing Reveals the Adaptive Role of the Epigenome in Three Deep-Sea Polychaetes

The roles of DNA methylation in invertebrates are poorly characterized, and critical data are missing for the phylum Annelida. We fill this knowledge gap by conducting the first genome-wide survey of DNA methylation in the deep-sea polychaetes dominant in deep-sea vents and seeps: Paraescarpia echinospica, Ridgeia piscesae, and Paralvinella palmiformis. DNA methylation calls were inferred from Oxford Nanopore sequencing after assembling high-quality genomes of these animals. The genomes of these worms encode all the key enzymes of the DNA methylation metabolism and possess a mosaic methylome similar to that of other invertebrates. Transcriptomic data of these polychaetes support the hypotheses that gene body methylation strengthens the expression of housekeeping genes and that promoter methylation acts as a silencing mechanism but not the hypothesis that DNA methylation suppresses the activity of transposable elements. The conserved epigenetic profiles of genes responsible for maintaining homeostasis under extreme hydrostatic pressure suggest DNA methylation plays an important adaptive role in these worms.

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.

Distinct genomic routes underlie transitions to specialised symbiotic lifestyles in deep-sea annelid worms

Bacterial symbioses allow annelids to colonise extreme ecological niches, such as hydrothermal vents and whale falls. Yet, the genetic principles sustaining these symbioses remain unclear. Here, we show that different genomic adaptations underpin the symbioses of phylogenetically related annelids with distinct nutritional strategies. Genome compaction and extensive gene losses distinguish the heterotrophic symbiosis of the bone-eating worm Osedax frankpressi from the chemoautotrophic symbiosis of deep-sea Vestimentifera. Osedax's endosymbionts complement many of the host's metabolic deficiencies, including the loss of pathways to recycle nitrogen and synthesise some amino acids. Osedax's endosymbionts possess the glyoxylate cycle, which could allow more efficient catabolism of bone-derived nutrients and the production of carbohydrates from fatty acids. Unlike in most Vestimentifera, innate immunity genes are reduced in O. frankpressi, which, however, has an expansion of matrix metalloproteases to digest collagen. Our study supports that distinct nutritional interactions influence host genome evolution differently in highly specialised symbioses.

Sulfur, sterol and trehalose metabolism in the deep-sea hydrocarbon seep tubeworm Lamellibrachia luymesi

BACKGROUND

Lamellibrachia luymesi dominates cold sulfide-hydrocarbon seeps and is known for its ability to consume bacteria for energy. The symbiotic relationship between tubeworms and bacteria with particular adaptations to chemosynthetic environments has received attention. However, metabolic studies have primarily focused on the mechanisms and pathways of the bacterial symbionts, while studies on the animal hosts are limited.

RESULTS

Here, we sequenced the transcriptome of L. luymesi and generated a transcriptomic database containing 79,464 transcript sequences. Based on GO and KEGG annotations, we identified transcripts related to sulfur metabolism, sterol biosynthesis, trehalose synthesis, and hydrolysis. Our in-depth analysis identified sulfation pathways in L. luymesi, and sulfate activation might be an important detoxification pathway for promoting sulfur cycling, reducing byproducts of sulfide metabolism, and converting sulfur compounds to sulfur-containing organics, which are essential for symbiotic survival. Moreover, sulfide can serve directly as a sulfur source for cysteine synthesis in L. luymesi. The existence of two pathways for cysteine synthesis might ensure its participation in the formation of proteins, heavy metal detoxification, and the sulfide-binding function of haemoglobin. Furthermore, our data suggested that cold-seep tubeworm is capable of de novo sterol biosynthesis, as well as incorporation and transformation of cycloartenol and lanosterol into unconventional sterols, and the critical enzyme involved in this process might have properties similar to those in the enzymes from plants or fungi. Finally, trehalose synthesis in L. luymesi occurs via the trehalose-6-phosphate synthase (TPS) and trehalose-6-phosphate phosphatase (TPP) pathways. The TPP gene has not been identified, whereas the TPS gene encodes a protein harbouring conserved TPS/OtsA and TPP/OtsB domains. The presence of multiple trehalases that catalyse trehalose hydrolysis could indicate the different roles of trehalase in cold-seep tubeworms.

CONCLUSIONS

We elucidated several molecular pathways of sulfate activation, cysteine and cholesterol synthesis, and trehalose metabolism. Contrary to the previous analysis, two pathways for cysteine synthesis and the cycloartenol-C-24-methyltransferase gene were identified in animals for the first time. The present study provides new insights into particular adaptations to chemosynthetic environments in L. luymesi and can serve as the basis for future molecular studies on host-symbiont interactions and biological evolution.

Hologenome analysis reveals independent evolution to chemosymbiosis by deep-sea bivalves

BACKGROUND

Bivalves have independently evolved a variety of symbiotic relationships with chemosynthetic bacteria. These relationships range from endo- to extracellular interactions, making them ideal for studies on symbiosis-related evolution. It is still unclear whether there are universal patterns to symbiosis across bivalves. Here, we investigate the hologenome of an extracellular symbiotic thyasirid clam that represents the early stages of symbiosis evolution.

RESULTS

We present a hologenome of Conchocele bisecta (Bivalvia: Thyasiridae) collected from deep-sea hydrothermal vents with extracellular symbionts, along with related ultrastructural evidence and expression data. Based on ultrastructural and sequencing evidence, only one dominant Thioglobaceae bacteria was densely aggregated in the large bacterial chambers of C. bisecta, and the bacterial genome shows nutritional complementarity and immune interactions with the host. Overall, gene family expansions may contribute to the symbiosis-related phenotypic variations in different bivalves. For instance, convergent expansions of gaseous substrate transport families in the endosymbiotic bivalves are absent in C. bisecta. Compared to endosymbiotic relatives, the thyasirid genome exhibits large-scale expansion in phagocytosis, which may facilitate symbiont digestion and account for extracellular symbiotic phenotypes. We also reveal that distinct immune system evolution, including expansion in lipopolysaccharide scavenging and contraction of IAP (inhibitor of apoptosis protein), may contribute to the different manners of bacterial virulence resistance in C. bisecta.

CONCLUSIONS

Thus, bivalves employ different pathways to adapt to the long-term co-existence with their bacterial symbionts, further highlighting the contribution of stochastic evolution to the independent gain of a symbiotic lifestyle in the lineage.

The genome of a vestimentiferan tubeworm (Ridgeia piscesae) provides insights into its adaptation to a deep-sea environment

BACKGROUND

Vestimentifera (Polychaeta, Siboglinidae) is a taxon of deep-sea worm-like animals living in deep-sea hydrothermal vents, cold seeps, and organic falls. The morphology and lifespan of Ridgeia piscesae, which is the only vestimentiferan tubeworm species found in the hydrothermal vents on the Juan de Fuca Ridge, vary greatly according to endemic environment. Recent analyses have revealed the genomic basis of adaptation in three vent- and seep-dwelling vestimentiferan tubeworms. However, the evolutionary history and mechanism of adaptation in R. piscesae, a unique species in the family Siboglinidae, remain to be investigated.

RESULT

We assembled a draft genome of R. piscesae collected at the Cathedral vent of the Juan de Fuca Ridge. Comparative genomic analysis showed that vent-dwelling tubeworms with a higher growth rate had smaller genome sizes than seep-dwelling tubeworms that grew much slower. A strong positive correlation between repeat content and genome size but not intron size and the number of protein-coding genes was identified in these deep-sea tubeworm species. Evolutionary analysis revealed that Ridgeia pachyptila and R. piscesae, the two tubeworm species that are endemic to hydrothermal vents of the eastern Pacific, started to diverge between 28.5 and 35 million years ago. Four genes involved in cell proliferation were found to be subject to positive selection in the genome of R. piscesae.

CONCLUSION

Ridgeia pachyptila and R. piscesae started to diverge after the formation of the Gorda/Juan de Fuca/Explorer ridge systems and the East Pacific Rise. The high growth rates of vent-dwelling tubeworms might be derived from their small genome sizes. Cell proliferation is important for regulating the growth rate in R. piscesae.

Metal bioaccumulation and transfer in benthic species based on isotopic signatures of individual amino acids in South China Sea cold seep environments

Cold seeps are deep-sea 'oases' with dense and dominant coexisting populations of large mussels and tubeworms under extreme environments. Under such natural source of high metal concentrations, the present study investigated the metal bioaccumulation and transfer with trophic positions in six benthic species by the isotopic δ¹⁵N and δ¹³C signatures in the active Haima cold seep, South China Sea. Comparing the isotopic signatures of bulk-tissue and amino acids by compound-specific isotopic analysis (CSIA-AA), we found that the bulk trophic (TP_B) values in the benthos except mussels were significantly higher than those of CSIA-based TP_Glu-Phe values. The estimated CSIA-based TP_Glu-Phe values showed a relatively compressed food chain with much changeable and unique amino acid isotopic heterogeneity, followed slim tubeworms (1.20)<mussels (1.38)<clams (1.52)<brittle stars (1.82)<giant tubeworms (2.16)<shrimps (2.31). All species accumulated relatively high concentrations of Fe, Zn, Cu, and Cr, especially for Zn in clams. Pearson correlation analysis showed that most metals had no significant relationship between their bioaccumulation and trophic positions, whereas Hg showed a significantly positive bioaccumulation through trophic transfer in such a compressed food chain. Water exposure was a major metal source rather than bacterial assimilation for most metals in the cold seep higher consumers. Hyperaccumulation of specific metals in some tissues of different benthos indicated different metal overflows in the Haima cold seep (As and Ni for tubeworms, Zn and Cd for clam gills, Ag and Cu for mussel gills). This study demonstrated high metal adaptations in different species and stable isotopic characteristics of amino acid metabolism in a natural high metal source of an active deep-sea cold seep, which is important for deep-sea development and environmental protection.

Immune diversity in lophotrochozoans, with a focus on recognition and effector systems

Lophotrochozoa is one of the most species-rich but immunologically poorly explored phyla. Although lack of acquired response in a narrow sense, lophotrochozoans possess various genetic mechanisms that enhance the diversity and specificity of innate immune system. Here, we review the recent advances of comparative immunology studies in lophotrochozoans with focus on immune recognition and effector systems. Haemocytes and coelomocytes are general important yet understudied player. Comparative genomics studies suggest expansion and functional divergence of lophotrochozoan immune reorganization systems is not as "homogeneous and simple" as we thought including the large-scale expansion and molecular divergence of pattern recognition receptors (PRRs) (TLRs, RLRs, lectins, etc.) and signaling adapters (MyD88s etc.), significant domain recombination of immune receptors (RLR, NLRs, lectins, etc.), extensive somatic recombination of fibrinogenrelated proteins (FREPs) in snails. Furthermore, there are repeatedly identified molecular mechanisms that generate immune effector diversity, including high polymorphism of antimicrobial peptides and proteins (AMPs), reactive oxygen and nitrogen species (RONS) and cytokines. Finally, we argue that the next generation omics tools and the recently emerged genome editing technicism will revolutionize our understanding of innate immune system in a comparative immunology perspective.

The Fox Gene Repertoire in the Annelid Owenia fusiformis Reveals Multiple Expansions of the foxQ2 Class in Spiralia

Fox genes are a large and conserved family of transcription factors involved in many key biological processes, including embryogenesis and body patterning. Although the role of Fox genes has been studied in an array of model systems, comprehensive comparative studies in Spiralia-a large clade of invertebrate animals including molluscs and annelids-are scarce but much needed to better understand the evolutionary history of this gene family. Here, we reconstruct and functionally characterize the Fox gene complement in the annelid Owenia fusiformis, a slow evolving species and member of the sister group to all remaining annelids. The genome of O. fusiformis contains at least a single ortholog for 20 of the 22 Fox gene classes that are ancestral to Bilateria, including an ortholog of the recently discovered foxT class. Temporal and spatial expression dynamics reveal a conserved role of Fox genes in gut formation, mesoderm patterning, and apical organ and cilia formation in Annelida and Spiralia. Moreover, we uncover an ancestral expansion of foxQ2 genes in Spiralia, represented by 11 paralogs in O. fusiformis. Notably, although all foxQ2 copies have apical expression in O. fusiformis, they show variable spatial domains and staggered temporal activation, which suggest cooperation and sub-functionalization among foxQ2 genes for the development of apical fates in this annelid. Altogether, our study informs the evolution and developmental roles of Fox genes in Annelida and Spiralia generally, providing the basis to explore how regulatory changes in Fox gene expression might have contributed to developmental and morphological diversification in Spiralia.

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.

Novel insights on obligate symbiont lifestyle and adaptation to chemosynthetic environment as revealed by the giant tubeworm genome

The mutualism between the giant tubeworm Riftia pachyptila and its endosymbiont Candidatus Endoriftia persephone has been extensively researched over the past 40 years. However, the lack of the host whole genome information has impeded the full comprehension of the genotype/phenotype interface in Riftia. Here we described the high-quality draft genome of Riftia, its complete mitogenome, and tissue-specific transcriptomic data. The Riftia genome presents signs of reductive evolution, with gene family contractions exceeding expansions. Expanded gene families are related to sulphur metabolism, detoxification, anti-oxidative stress, oxygen transport, immune system, and lysosomal digestion, reflecting evolutionary adaptations to the vent environment and endosymbiosis. Despite the derived body plan, the developmental gene repertoire in the gutless tubeworm is extremely conserved with the presence of a near intact and complete Hox cluster. Gene expression analyses establishes that the trophosome is a multi-functional organ marked by intracellular digestion of endosymbionts, storage of excretory products and haematopoietic functions. Overall, the plume and gonad tissues both in contact to the environment harbour highly expressed genes involved with cell cycle, programmed cell death, and immunity indicating a high cell turnover and defence mechanisms against pathogens. We posit that the innate immune system plays a more prominent role into the establishment of the symbiosis during the infection in the larval stage, rather than maintaining the symbiostasis in the trophosome. This genome bridges four decades of physiological research in Riftia, whilst simultaneously provides new insights into the development, whole organism functions and evolution in the giant tubeworm.

Comparative genomic and transcriptomic analyses of transposable elements in polychaetous annelids highlight LTR retrotransposon diversity and evolution

Background

With the expansion of high throughput sequencing, we now have access to a larger number of genome-wide studies analyzing the Transposable elements (TEs) composition in a wide variety of organisms. However, genomic analyses often remain too limited in number and diversity of species investigated to study in depth the dynamics and evolutionary success of the different types of TEs among metazoans. Therefore, we chose to investigate the use of transcriptomes to describe the diversity of TEs in phylogenetically related species by conducting the first comparative analysis of TEs in two groups of polychaetes and evaluate the diversity of TEs that might impact genomic evolution as a result of their mobility.

Results

We present a detailed analysis of TEs distribution in transcriptomes extracted from 15 polychaetes depending on the number of reads used during assembly, and also compare these results with additional TE scans on associated low-coverage genomes. We then characterized the clades defined by 1021 LTR-retrotransposon families identified in 26 species. Clade richness was highly dependent on the considered superfamily. Copia elements appear rare and are equally distributed in only three clades, GalEa, Hydra and CoMol. Among the eight BEL/Pao clades identified in annelids, two small clades within the Sailor lineage are new for science. We characterized 17 Gypsy clades of which only 4 are new; the C-clade largely dominates with a quarter of the families. Finally, all species also expressed for the majority two distinct transcripts encoding PIWI proteins, known to be involved in control of TEs mobilities.

Conclusions

This study shows that the use of transcriptomes assembled from 40 million reads was sufficient to access to the diversity and proportion of the transposable elements compared to those obtained by low coverage sequencing. Among LTR-retrotransposons Gypsy elements were unequivocally dominant but results suggest that the number of Gypsy clades, although high, may be more limited than previously thought in metazoans. For BEL/Pao elements, the organization of clades within the Sailor lineage appears more difficult to establish clearly. The Copia elements remain rare and result from the evolutionary consistent success of the same three clades.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13100-021-00252-0.

Genomic signatures supporting the symbiosis and formation of chitinous tube in the deep-sea tubeworm Paraescarpia echinospica

Vestimentiferan tubeworms are iconic animals that present as large habitat-forming chitinized tube bushes in deep-sea chemosynthetic ecosystems. They are gutless and depend entirely on their endosymbiotic sulfide-oxidizing chemoautotrophic bacteria for nutrition. Information on the genomes of several siboglinid endosymbionts has improved our understanding of their nutritional supplies. However, the interactions between tubeworms and their endosymbionts remain largely unclear due to a paucity of host genomes. Here, we report the chromosome-level genome of the vestimentiferan tubeworm Paraescarpia echinospica. We found that the genome has been remodeled to facilitate symbiosis through the expansion of gene families related to substrate transfer and innate immunity, suppression of apoptosis, regulation of lysosomal digestion, and protection against oxidative stress. Furthermore, the genome encodes a programmed cell death pathway that potentially controls the endosymbiont population. Our integrated genomic, transcriptomic, and proteomic analyses uncovered matrix proteins required for the formation of the chitinous tube and revealed gene family expansion and co-option as evolutionary mechanisms driving the acquisition of this unique supporting structure for deep-sea tubeworms. Overall, our study provides novel insights into the host’s support system that has enabled tubeworms to establish symbiosis, thrive in deep-sea hot vents and cold seeps, and produce the unique chitinous tubes in the deep sea.

Genome-wide characterization of LTR retrotransposons in the non-model deep-sea annelid Lamellibrachia luymesi

Background

Long Terminal Repeat retrotransposons (LTR retrotransposons) are mobile genetic elements composed of a few genes between terminal repeats and, in some cases, can comprise over half of a genome’s content. Available data on LTR retrotransposons have facilitated comparative studies and provided insight on genome evolution. However, data are biased to model systems and marine organisms, including annelids, have been underrepresented in transposable elements studies. Here, we focus on genome of Lamellibrachia luymesi, a vestimentiferan tubeworm from deep-sea hydrocarbon seeps, to gain knowledge of LTR retrotransposons in a deep-sea annelid.

Results

We characterized LTR retrotransposons present in the genome of L. luymesi using bioinformatic approaches and found that intact LTR retrotransposons makes up about 0.1% of L. luymesi genome. Previous characterization of the genome has shown that this tubeworm hosts several known LTR-retrotransposons. Here we describe and classify LTR retrotransposons in L. luymesi as within the Gypsy, Copia and Bel-pao superfamilies. Although, many elements fell within already recognized families (e.g., Mag, CSRN1), others formed clades distinct from previously recognized families within these superfamilies. However, approximately 19% (41) of recovered elements could not be classified. Gypsy elements were the most abundant while only 2 Copia and 2 Bel-pao elements were present. In addition, analysis of insertion times indicated that several LTR-retrotransposons were recently transposed into the genome of L. luymesi, these elements had identical LTR’s raising possibility of recent or ongoing retrotransposon activity.

Conclusions

Our analysis contributes to knowledge on diversity of LTR-retrotransposons in marine settings and also serves as an important step to assist our understanding of the potential role of retroelements in marine organisms. We find that many LTR retrotransposons, which have been inserted in the last few million years, are similar to those found in terrestrial model species. However, several new groups of LTR retrotransposons were discovered suggesting that the representation of LTR retrotransposons may be different in marine settings. Further study would improve understanding of the diversity of retrotransposons across animal groups and environments.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07749-1.

Contrasting influences on bacterial symbiont specificity by co-occurring deep-sea mussels and tubeworms

Relationships fueled by sulfide between deep-sea invertebrates and bacterial symbionts are well known, yet the diverse overlapping factors influencing symbiont specificity are complex. For animals that obtain their symbionts from the environment, both host identity and geographic location can impact the ultimate symbiont partner. Bacterial symbionts were analysed for three co-occurring species each of Bathymodiolus mussels and vestimentiferan tubeworms, from three deep methane seeps off the west coast of Costa Rica. The bacterial internal transcribed spacer gene was analysed via direct and barcoded amplicon sequencing to reveal fine-scale symbiont diversity. Each of the three mussel species (B. earlougheri, B. billschneideri and B. nancyschneideri) hosted genetically distinct thiotrophic endosymbionts, despite living nearly side-by-side in their habitat, suggesting that host identity is crucial in driving symbiont specificity. The dominant thiotrophic symbiont of co-occurring tubeworms Escarpia spicata and Lamellibrachia (L. barhami and L. donwalshi), on the other hand, was identical regardless of host species or sample location, suggesting lack of influence by either factor on symbiont selectivity in this group of animals. These findings highlight the specific, yet distinct, influences on the environmental acquisition of symbionts in two foundational invertebrates with similar lifestyles, and provide a rapid, precise method of examining symbiont identities.

Hologenome analysis reveals dual symbiosis in the deep-sea hydrothermal vent snail Gigantopelta aegis

Animals endemic to deep-sea hydrothermal vents often form obligatory symbioses with bacteria, maintained by intricate host-symbiont interactions. Most genomic studies on holobionts have not investigated both sides to similar depths. Here, we report dual symbiosis in the peltospirid snail Gigantopelta aegis with two gammaproteobacterial endosymbionts: a sulfur oxidiser and a methane oxidiser. We assemble high-quality genomes for all three parties, including a chromosome-level host genome. Hologenomic analyses reveal mutualism with nutritional complementarity and metabolic co-dependency, highly versatile in transporting and using chemical energy. Gigantopelta aegis likely remodels its immune system to facilitate dual symbiosis. Comparisons with Chrysomallon squamiferum, a confamilial snail with a single sulfur-oxidising gammaproteobacterial endosymbiont, show that their sulfur-oxidising endosymbionts are phylogenetically distant. This is consistent with previous findings that they evolved endosymbiosis convergently. Notably, the two sulfur-oxidisers share the same capabilities in biosynthesising nutrients lacking in the host genomes, potentially a key criterion in symbiont selection.

The complete mitochondrial genome of a cold seep gastropod Phymorhynchus buccinoides (Neogastropoda: Conoidea: Raphitomidae)

Phymorhynchus is a genus of deep-sea snails that are most distributed in hydrothermal vent or cold seep environments. In this study, we presented the complete mitochondrial genome of P. buccinoides, a cold seep snail from the South China Sea. It is the first mitochondrial genome of a cold seep member of the superfamily Conoidea. The mitochondrial genome is 15,764 bp in length, and contains 13 protein-coding genes (PCGs), 2 rRNA genes, and 22 tRNA genes. These genes are encoded on the positive strand, except for 8 tRNA genes that are encoded on the negative strand. The start codon ATG and 3 types of stop codons, TAA, TAG and the truncated termination codon T, are used in the 13 PCGs. All 13 PCGs in the 26 species of Conoidea share the same gene order, while several tRNA genes have been translocated. Phylogenetic analysis revealed that P. buccinoides clustered with Typhlosyrinx sp., Eubela sp., and Phymorhynchus sp., forming the Raphitomidae clade, with high support values. Positive selection analysis showed that a residue located in atp6 (18 S) was identified as the positively selected site with high posterior probabilities, suggesting potential adaption to the cold seep environment. Overall, our data will provide a useful resource on the evolutionary adaptation of cold seep snails for future studies.

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

Cold seeps, characterized by the methane, hydrogen sulfide, and other hydrocarbon chemicals, foster one of the most widespread chemosynthetic ecosystems in deep sea that are densely populated by specialized benthos. However, scarce genomic resources severely limit our knowledge about the origin and adaptation of life in this unique ecosystem. Here, we present a genome of a deep-sea limpet Bathyacmaea lactea, a common species associated with the dominant mussel beds in cold seeps. We yielded 54.6 gigabases (Gb) of Nanopore reads and 77.9-Gb BGI-seq raw reads, respectively. Assembly harvested a 754.3-Mb genome for B. lactea, with 3,720 contigs and a contig N50 of 1.57 Mb, covering 94.3% of metazoan Benchmarking Universal Single-Copy Orthologs. In total, 23,574 protein-coding genes and 463.4 Mb of repetitive elements were identified. We analyzed the phylogenetic position, substitution rate, demographic history, and TE activity of B. lactea. We also identified 80 expanded gene families and 87 rapidly evolving Gene Ontology categories in the B. lactea genome. Many of these genes were associated with heterocyclic compound metabolism, membrane-bounded organelle, metal ion binding, and nitrogen and phosphorus metabolism. The high-quality assembly and in-depth characterization suggest the B. lactea genome will serve as an essential resource for understanding the origin and adaptation of life in the cold seeps.

Genomic adaptations to chemosymbiosis in the deep-sea seep-dwelling tubeworm Lamellibrachia luymesi

BACKGROUND

Symbiotic relationships between microbes and their hosts are widespread and diverse, often providing protection or nutrients, and may be either obligate or facultative. However, the genetic mechanisms allowing organisms to maintain host-symbiont associations at the molecular level are still mostly unknown, and in the case of bacterial-animal associations, most genetic studies have focused on adaptations and mechanisms of the bacterial partner. The gutless tubeworms (Siboglinidae, Annelida) are obligate hosts of chemoautotrophic endosymbionts (except for Osedax which houses heterotrophic Oceanospirillales), which rely on the sulfide-oxidizing symbionts for nutrition and growth. Whereas several siboglinid endosymbiont genomes have been characterized, genomes of hosts and their adaptations to this symbiosis remain unexplored.

RESULTS

Here, we present and characterize adaptations of the cold seep-dwelling tubeworm Lamellibrachia luymesi, one of the longest-lived solitary invertebrates. We sequenced the worm's ~ 688-Mb haploid genome with an overall completeness of ~ 95% and discovered that L. luymesi lacks many genes essential in amino acid biosynthesis, obligating them to products provided by symbionts. Interestingly, the host is known to carry hydrogen sulfide to thiotrophic endosymbionts using hemoglobin. We also found an expansion of hemoglobin B1 genes, many of which possess a free cysteine residue which is hypothesized to function in sulfide binding. Contrary to previous analyses, the sulfide binding mediated by zinc ions is not conserved across tubeworms. Thus, the sulfide-binding mechanisms in sibgolinids need to be further explored, and B1 globins might play a more important role than previously thought. Our comparative analyses also suggest the Toll-like receptor pathway may be essential for tolerance/sensitivity to symbionts and pathogens. Several genes related to the worm's unique life history which are known to play important roles in apoptosis, cell proliferation, and aging were also identified. Last, molecular clock analyses based on phylogenomic data suggest modern siboglinid diversity originated in 267 mya (± 70 my) support previous hypotheses indicating a Late Mesozoic or Cenozoic origins of approximately 50-126 mya for vestimentiferans.

CONCLUSIONS

Here, we elucidate several specific adaptations along various molecular pathways that link phenome to genome to improve understanding of holobiont evolution. Our findings of adaptation in genomic mechanisms to reducing environments likely extend to other chemosynthetic symbiotic systems.