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

Potential biomarkers of metal toxicity in deep-sea invertebrates - A critical review of the omics data

Deep-sea mining (DSM) activities are expected to release potentially toxic metal mixtures through the generation of sediment plumes to the marine environment. This may disrupt the normal functioning of biological mechanisms, adversely affecting deep-sea invertebrate organisms. It is thus essential to understand the ecotoxicological effects from these toxic elements in deep-sea organisms and the omics approaches applied to ecotoxicology are seen as promising tools. Here, we provide an overview of the principal biological modifications identified in deep-sea invertebrates when exposed to metals and critically evaluate the current knowledge and discuss which potential biomarkers may be useful after metal exposure. Most of the 50 omics studies on deep-sea invertebrates revised are comparative transcriptomes (n = 41). Forty-three potential biomarker candidates are highlighted from immune system, 46 from cellular metabolism and 29 from oxidative stress. The processes mostly affected by metal toxicity in deep-sea invertebrates are related to innate immune defense; sulfur, chitin, and catabolic metabolism; antioxidation; and detoxification. We acknowledge the current limitations and future perspectives for their uses and emphasize the need to invest in further ecotoxicological studies using the omics approaches.

MolluscDB 2.0: a comprehensive functional and evolutionary genomics database for over 1400 molluscan species

Mollusca represents the second-largest animal phylum but remains less explored genomically. The increase in high-quality genomes and diverse functional genomic data holds great promise for advancing our understanding of molluscan biology and evolution. To address the opportunities and challenges facing the molluscan research community in managing vast multi-omics resources, we developed MolluscDB 2.0 (http://mgbase.qnlm.ac), which integrates extensive functional genomic data and offers user-friendly tools for multilevel integrative and comparative analyses. MolluscDB 2.0 covers 1450 species across all eight molluscan classes and compiles ∼4200 datasets, making it the most comprehensive multi-omics resource for molluscs to date. MolluscDB 2.0 expands the layers of multi-omics data, including genomes, bulk transcriptomes, single-cell transcriptomes, proteomes, epigenomes and metagenomes. MolluscDB 2.0 also more than doubles the number of functional modules and analytical tools, updating 14 original modules and introducing 20 new, specialized modules. Overall, MolluscDB 2.0 provides highly valuable, open-access multi-omics platform for the molluscan research community, expediting scientific discoveries and deepening our understanding of molluscan biology and evolution.

Chromosome-level genome assembly of the ivory shell Babylonia areolata

The ivory shell Babylonia areolata is an economically important marine benthic gastropod known for its rapid growth and high nutritional value. B. areolata is distributed in Southeast Asia and the southeast coastal areas of China. In this study, we constructed a high-quality genome for B. areolata using PacBio, Illumina, and Hi-C sequencing technologies. The genome assembly comprised 35 chromosomal sequences with a total length of 1.65 Gb. The scaffold and contig N50 lengths were 53.17 Mb and 2.64 Mb, respectively, with repeat sequences constituting 64.46% of the genome. Furthermore, 26,130 protein-coding genes and 96.75% of the genome's BUSCOs were identified. This inaugural report of a B. areolata genome provides crucial foundational information for further investigations into the biology, genomics, and genetic improvement of economic traits of this species.

Chromosome-level genome assembly of the sacoglossan sea slug Elysia timida (Risso, 1818)

BACKGROUND

Sequencing and annotating genomes of non-model organisms helps to understand genome architecture, the genetic processes underlying species traits, and how these genes have evolved in closely-related taxa, among many other biological processes. However, many metazoan groups, such as the extremely diverse molluscs, are still underrepresented in the number of sequenced and annotated genomes. Although sequencing techniques have recently improved in quality and quantity, molluscs are still neglected due to difficulties in applying standardized protocols for obtaining genomic data.

RESULTS

In this study, we present the chromosome-level genome assembly and annotation of the sacoglossan sea slug species Elysia timida, known for its ability to store the chloroplasts of its food algae. In particular, by optimizing the long-read and chromosome conformation capture library preparations, the genome assembly was performed using PacBio HiFi and Arima HiC data. The scaffold and contig N50s, at 41.8 Mb and 1.92 Mb, respectively, are approximately 30-fold and fourfold higher compared to other published sacoglossan genome assemblies. Structural annotation resulted in 19,904 protein-coding genes, which are more contiguous and complete compared to publicly available annotations of Sacoglossa with respect to metazoan BUSCOs. We found no evidence for horizontal gene transfer (HGT), i.e. no photosynthetic genes encoded in the sacoglossan nucleus genome. However, we detected genes encoding polyketide synthases in E. timida, indicating that polypropionates are produced. HPLC-MS/MS analysis confirmed the presence of a large number of polypropionates, including known and yet uncharacterised compounds.

CONCLUSIONS

We can show that our methodological approach helps to obtain a high-quality genome assembly even for a "difficult-to-sequence" organism, which may facilitate genome sequencing in molluscs. This will enable a better understanding of complex biological processes in molluscs, such as functional kleptoplasty in Sacoglossa, by significantly improving the quality of genome assemblies and annotations.

Genetic adaptations of marine invertebrates to hydrothermal vent habitats

Hydrothermal vents are unique habitats like an oases of life compared with typical deep-sea, soft-sediment environments. Most animals that live in these habitats are invertebrates, and they have adapted to extreme vent environments that include high temperatures, hypoxia, high sulfide, high metal concentration, and darkness. The advent of next-generation sequencing technology, especially the coming of the new era of omics, allowed more studies to focus on the molecular adaptation of these invertebrates to vent habitats. Many genes linked to hydrothermal adaptation have been studied. We summarize the findings related to these genetic adaptations and discuss which new techniques can facilitate studies in the future.

Tracking the hologenome dynamics in aquatic invertebrates by the holo-2bRAD approach

The "hologenome" concept is an increasingly popular way of thinking about microbiome-host for marine organisms. However, it is challenging to track hologenome dynamics because of the large amount of material, with tracking itself usually resulting in damage or death of the research object. Here we show the simple and efficient holo-2bRAD approach for the tracking of hologenome dynamics in marine invertebrates (i.e., scallop and shrimp) from one holo-2bRAD library. The stable performance of our approach was shown with high genotyping accuracy of 99.91% and a high correlation of r > 0.99 for the species-level profiling of microorganisms. To explore the host-microbe association underlying mass mortality events of bivalve larvae, core microbial species changed with the stages were found, and two potentially associated host SNPs were identified. Overall, our research provides a powerful tool with various advantages (e.g., cost-effective, simple, and applicable for challenging samples) in genetic, ecological, and evolutionary studies.

Environmental specificity of karst cave habitats evidenced by diverse symbiotic bacteria in Opiliones

BACKGROUND

Karst caves serve as natural laboratories, providing organisms with extreme and constant conditions that promote isolation, resulting in a genetic relationship and living environment that is significantly different from those outside the cave. However, research on cave creatures, especially Opiliones, remains scarce, with most studies focused on water, soil, and cave sediments.

RESULTS

The structure of symbiotic bacteria in different caves were compared, revealing significant differences. Based on the alpha and beta diversity, symbiotic bacteria abundance and diversity in the cave were similar, but the structure of symbiotic bacteria differed inside and outside the cave. Microorganisms in the cave play an important role in material cycling and energy flow, particularly in the nitrogen cycle. Although microbial diversity varies inside and outside the cave, Opiliones in Beijing caves and Hainan Island exhibited a strong similarity, indicating that the two environments share commonalities.

CONCLUSIONS

The karst cave environment possesses high microbial diversity and there are noticeable differences among different caves. Different habitats lead to significant differences in the symbiotic bacteria in Opiliones inside and outside the cave, and cave microorganisms have made efforts to adapt to extreme environments. The similarity in symbiotic bacteria community structure suggests a potential similarity in host environments, providing an explanation for the appearance of Sinonychia martensi in caves in the north.

Chromosome-scale Genome Assembly of the Rough Periwinkle Littorina saxatilis

The intertidal gastropod Littorina saxatilis is a model system to study speciation and local adaptation. The repeated occurrence of distinct ecotypes showing different levels of genetic divergence makes L. saxatilis particularly suited to study different stages of the speciation continuum in the same lineage. A major finding is the presence of several large chromosomal inversions associated with the divergence of ecotypes and, specifically, the species offers a system to study the role of inversions in this divergence. The genome of L. saxatilis is 1.35 Gb and composed of 17 chromosomes. The first reference genome of the species was assembled using Illumina data, was highly fragmented (N50 of 44 kb), and was quite incomplete, with a BUSCO completeness of 80.1% on the Metazoan dataset. A linkage map of one full-sibling family enabled the placement of 587 Mbp of the genome into 17 linkage groups corresponding to the haploid number of chromosomes, but the fragmented nature of this reference genome limited the understanding of the interplay between divergent selection and gene flow during ecotype formation. Here, we present a newly generated reference genome that is highly contiguous, with a N50 of 67 Mb and 90.4% of the total assembly length placed in 17 super-scaffolds. It is also highly complete with a BUSCO completeness of 94.1% of the Metazoa dataset. This new reference will allow for investigations into the genomic regions implicated in ecotype formation as well as better characterization of the inversions and their role in speciation.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Genome assembly and population genomic data of a pulmonate snail Ellobium chinense

Ellobium chinense is an airbreathing, pulmonate gastropod species that inhabits saltmarshes in estuaries of the northwestern Pacific. Due to a rapid population decline and their unique ecological niche in estuarine ecosystems, this species has attracted special attention regarding their conservation and the genomic basis of adaptation to frequently changing environments. Here we report a draft genome assembly of E. chinense with a total size of 949.470 Mb and a scaffold N50 of 1.465 Mb. Comparative genomic analysis revealed that the GO terms enriched among four gastropod species are related to signal transduction involved in maintaining electrochemical gradients across the cell membrane. Population genomic analysis using the MSMC model for 14 re-sequenced individuals revealed a drastic decline in Korean and Japanese populations during the last glacial period, while the southern Chinese population retained a much larger effective population size (Ne). These contrasting demographic changes might be attributed to multiple environmental factors during the glacial-interglacial cycles. This study provides valuable genomic resources for understanding adaptation and historical demographic responses to climate change.

Taxonomic Distribution and Molecular Evolution of Mytilectins

R-type lectins are a widespread group of sugar-binding proteins found in nearly all domains of life, characterized by the presence of a carbohydrate-binding domain that adopts a β-trefoil fold. Mytilectins represent a recently described subgroup of β-trefoil lectins, which have been functionally characterized in a few mussel species (Mollusca, Bivalvia) and display attractive properties, which may fuel the development of artificial lectins with different biotechnological applications. The detection of different paralogous genes in mussels, together with the description of orthologous sequences in brachiopods, supports the formal description of mytilectins as a gene family. However, to date, an investigation of the taxonomic distribution of these lectins and their molecular diversification and evolution was still lacking. Here, we provide a comprehensive overview of the evolutionary history of mytilectins, revealing an ancient monophyletic evolutionary origin and a very broad but highly discontinuous taxonomic distribution, ranging from heteroscleromorphan sponges to ophiuroid and crinoid echinoderms. Moreover, the overwhelming majority of mytilectins display a chimera-like architecture, which combines the β-trefoil carbohydrate recognition domain with a C-terminal pore-forming domain, suggesting that the simpler structure of most functionally characterized mytilectins derives from a secondary domain loss.

Mosaic environment-driven evolution of the deep-sea mussel Gigantidas platifrons bacterial endosymbiont

BACKGROUND

The within-species diversity of symbiotic bacteria represents an important genetic resource for their environmental adaptation, especially for horizontally transmitted endosymbionts. Although strain-level intraspecies variation has recently been detected in many deep-sea endosymbionts, their ecological role in environmental adaptation, their genome evolution pattern under heterogeneous geochemical environments, and the underlying molecular forces remain unclear.

RESULTS

Here, we conducted a fine-scale metagenomic analysis of the deep-sea mussel Gigantidas platifrons bacterial endosymbiont collected from distinct habitats: hydrothermal vent and methane seep. Endosymbiont genomes were assembled using a pipeline that distinguishes within-species variation and revealed highly heterogeneous compositions in mussels from different habitats. Phylogenetic analysis separated the assemblies into three distinct environment-linked clades. Their functional differentiation follows a mosaic evolutionary pattern. Core genes, essential for central metabolic function and symbiosis, were conserved across all clades. Clade-specific genes associated with heavy metal resistance, pH homeostasis, and nitrate utilization exhibited signals of accelerated evolution. Notably, transposable elements and plasmids contributed to the genetic reshuffling of the symbiont genomes and likely accelerated adaptive evolution through pseudogenization and the introduction of new genes.

CONCLUSIONS

The current study uncovers the environment-driven evolution of deep-sea symbionts mediated by mobile genetic elements. Its findings highlight a potentially common and critical role of within-species diversity in animal-microbiome symbioses. Video Abstract.

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.

Current knowledge of the Southern Hemisphere marine microbiome in eukaryotic hosts and the Strait of Magellan surface microbiome project

Host-microbe interactions are ubiquitous and play important roles in host biology, ecology, and evolution. Yet, host-microbe research has focused on inland species, whereas marine hosts and their associated microbes remain largely unexplored, especially in developing countries in the Southern Hemisphere. Here, we review the current knowledge of marine host microbiomes in the Southern Hemisphere. Our results revealed important biases in marine host species sampling for studies conducted in the Southern Hemisphere, where sponges and marine mammals have received the greatest attention. Sponge-associated microbes vary greatly across geographic regions and species. Nevertheless, besides taxonomic heterogeneity, sponge microbiomes have functional consistency, whereas geography and aging are important drivers of marine mammal microbiomes. Seabird and macroalgal microbiomes in the Southern Hemisphere were also common. Most seabird microbiome has focused on feces, whereas macroalgal microbiome has focused on the epibiotic community. Important drivers of seabird fecal microbiome are aging, sex, and species-specific factors. In contrast, host-derived deterministic factors drive the macroalgal epibiotic microbiome, in a process known as "microbial gardening". In turn, marine invertebrates (especially crustaceans) and fish microbiomes have received less attention in the Southern Hemisphere. In general, the predominant approach to study host marine microbiomes has been the sequencing of the 16S rRNA gene. Interestingly, there are some marine holobiont studies (i.e., studies that simultaneously analyze host (e.g., genomics, transcriptomics) and microbiome (e.g., 16S rRNA gene, metagenome) traits), but only in some marine invertebrates and macroalgae from Africa and Australia. Finally, we introduce an ongoing project on the surface microbiome of key species in the Strait of Magellan. This is an international project that will provide novel microbiome information of several species in the Strait of Magellan. In the short-term, the project will improve our knowledge about microbial diversity in the region, while long-term potential benefits include the use of these data to assess host-microbial responses to the Anthropocene derived climate change.

Growing in Saltwater: Biotechnological Potential of Novel Methylotuvimicrobium- and Methylomarinum-like Methanotrophic Bacteria

Methanotrophic bacteria that possess a unique ability of using methane as a sole source of carbon and energy have attracted considerable attention as potential producers of a single-cell protein. So far, this biotechnology implied using freshwater methanotrophs, although many regions of the world have limited freshwater resources. This study aimed at searching for novel methanotrophs capable of fast growth in saltwater comparable in composition with seawater. A methane-oxidizing microbial consortium containing Methylomarinum- and Methylotuvimicrobium-like methanotrophs was enriched from sediment from the river Chernavka (water pH 7.5, total salt content 30 g L-1), a tributary river of the hypersaline Lake Elton, southern Russia. This microbial consortium, designated Ch1, demonstrated stable growth on natural gas in a bioreactor in media with a total salt content of 23 to 35.9 g L-1 at a dilution rate of 0.19-0.21 h-1. The highest biomass yield of 5.8 g cell dry weight (CDW)/L with a protein content of 63% was obtained during continuous cultivation of the consortium Ch1 in a medium with a total salt content of 29 g L-1. Isolation attempts resulted in obtaining a pure culture of methanotrophic bacteria, strain Ch1-1. The 16S rRNA gene sequence of strain Ch1-1 displayed 97.09-97.24% similarity to the corresponding gene fragments of characterized representatives of Methylomarinum vadi, methanotrophs isolated from marine habitats. The genome of strain Ch1-1 was 4.8 Mb in size and encoded 3 rRNA operons, and about 4400 proteins. The genome contained the gene cluster coding for ectoine biosynthesis, which explains the ability of strain Ch1-1 to tolerate high salt concentration.

Transcriptome analysis to identify candidate genes related to mammary gland development of Bactrian camel (Camelus bactrianus)

Introduction

The demand for camel milk, which has unique therapeutic properties, is increasing. The mammary gland is the organ in mammals responsible for the production and quality of milk. However, few studies have investigated the genes or pathways related to mammary gland growth and development in Bactrian camels. This study aimed to compare the morphological changes in mammary gland tissue and transcriptome expression profiles between young and adult female Bactrian camels and to explore the potential candidate genes and signaling pathways related to mammary gland development.

Methods

Three 2  years-old female camels and three 5  years-old adult female camels were maintained in the same environment. The parenchyma of the mammary gland tissue was sampled from the camels using percutaneous needle biopsy. Morphological changes were observed using hematoxylin-eosin staining. High-throughput RNA sequencing was performed using the Illumina HiSeq platform to analyze changes in the transcriptome between young and adult camels. Functional enrichment, pathway enrichment, and protein-protein interaction networks were also analyzed. Gene expression was verified using quantitative real-time polymerase chain reaction (qRT-PCR).

Results

Histomorphological analysis showed that the mammary ducts and mammary epithelial cells in adult female camels were greatly developed and differentiated from those in young camels. Transcriptome analysis showed that 2,851 differentially expressed genes were obtained in the adult camel group compared to the young camel group, of which 1,420 were upregulated, 1,431 were downregulated, and 2,419 encoded proteins. Functional enrichment analysis revealed that the upregulated genes were significantly enriched for 24 pathways, including the Hedgehog signaling pathway which is closely related to mammary gland development. The downregulated genes were significantly enriched for seven pathways, among these the Wnt signaling pathway was significantly related to mammary gland development. The protein-protein interaction network sorted the nodes according to the degree of gene interaction and identified nine candidate genes: PRKAB2, PRKAG3, PLCB4, BTRC, GLI1, WIF1, DKK2, FZD3, and WNT4. The expression of fifteen genes randomly detected by qRT-PCR showed results consistent with those of the transcriptome analysis.

Discussion

Preliminary findings indicate that the Hedgehog, Wnt, oxytocin, insulin, and steroid biosynthesis signaling pathways have important effects on mammary gland development in dairy camels. Given the importance of these pathways and the interconnections of the involved genes, the genes in these pathways should be considered potential candidate genes. This study provides a theoretical basis for elucidating the molecular mechanisms associated with mammary gland development and milk production in Bactrian camels.

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.

Geography, not lifestyle, explains the population structure of free-living and host-associated deep-sea hydrothermal vent snail symbionts

BACKGROUND

Marine symbioses are predominantly established through horizontal acquisition of microbial symbionts from the environment. However, genetic and functional comparisons of free-living populations of symbionts to their host-associated counterparts are sparse. Here, we assembled the first genomes of the chemoautotrophic gammaproteobacterial symbionts affiliated with the deep-sea snail Alviniconcha hessleri from two separate hydrothermal vent fields of the Mariana Back-Arc Basin. We used phylogenomic and population genomic methods to assess sequence and gene content variation between free-living and host-associated symbionts.

RESULTS

Our phylogenomic analyses show that the free-living and host-associated symbionts of A. hessleri from both vent fields are populations of monophyletic strains from a single species. Furthermore, genetic structure and gene content analyses indicate that these symbiont populations are differentiated by vent field rather than by lifestyle.

CONCLUSION

Together, this work suggests that, despite the potential influence of host-mediated acquisition and release processes on horizontally transmitted symbionts, geographic isolation and/or adaptation to local habitat conditions are important determinants of symbiont population structure and intra-host composition. Video Abstract.

Duplication and Losses of Opsin Genes in Lophotrochozoan Evolution

Opsins are G-coupled receptors playing a key role in metazoan visual processes. While many studies enriched our understanding of opsin diversity in several animal clades, the opsin evolution in Lophotrochozoa, one of the major metazoan groups, remains poorly understood. Using recently developed phylogenetic approaches, we investigated the opsin evolution in 74 lophotrochozoan genomes. We found that the common ancestor of Lophotrochozoa possessed at least seven opsin paralog groups that underwent divergent evolutionary history in the different phyla. Furthermore, we showed for the first time opsin-related molecules in Bilateria that we named pseudopsins, which may prove critical in uncovering opsin 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.

Gene loss and symbiont switching during adaptation to the deep sea in a globally distributed symbiosis

Chemosynthetic symbioses between bacteria and invertebrates occur worldwide from coastal sediments to the deep sea. Most host groups are restricted to either shallow or deep waters. In contrast, Lucinidae, the most species-rich family of chemosymbiotic invertebrates, has both shallow- and deep-sea representatives. Multiple lucinid species have independently colonized the deep sea, which provides a unique framework for understanding the role microbial symbionts play in evolutionary transitions between shallow and deep waters. Lucinids acquire their symbionts from their surroundings during early development, which may allow them to flexibly acquire symbionts that are adapted to local environments. Via metagenomic analyses of museum and other samples collected over decades, we investigated the biodiversity and metabolic capabilities of the symbionts of 22 mostly deep-water lucinid species. We aimed to test the theory that the symbiont played a role in adaptation to life in deep-sea habitats. We identified 16 symbiont species, mostly within the previously described genus Ca. Thiodiazotropha. Most genomic functions were shared by both shallow-water and deep-sea Ca. Thiodiazotropha, though nitrogen fixation was exclusive to shallow-water species. We discovered multiple cases of symbiont switching near deep-sea hydrothermal vents and cold seeps, where distantly related hosts convergently acquired novel symbionts from a different bacterial order. Finally, analyses of selection revealed consistently stronger purifying selection on symbiont genomes in two extreme habitats - hydrothermal vents and an oxygen-minimum zone. Our findings reveal that shifts in symbiont metabolic capability and, in some cases, acquisition of a novel symbiont accompanied adaptation of lucinids to challenging deep-sea habitats.

The transcriptomic analysis of Planorbarius corneus hemocytes (Gastropoda) naturally infected with Bilharziella polonica (Schistosomatidae)

The study of molluscan innate immunity is essential for understanding the evolution of the immune system. An advance in the knowledge of their immune system can be achieved by increasing the number of model species. Our study focuses on the immunity of Planorbarius corneus, a pulmonate snail widely distributed in Eurasia. These snails are intermediate hosts of many trematodes, including Bilharziella polonica (Schistosomatidae). In this paper we obtained and analyzed transcriptomes of hemocytes of uninfected snails Planorbarius corneus and snails naturally infected with Bilharziella polonica. The transcriptomes were found to contain transcripts encoding all major groups of immune factors previously described for other gastropods. Pathogen-recognition molecules were the most diverse group of immune factors. Comparison of the transcriptomes of the infected and the uninfected molluscs showed that the expression of some genes changed during infection. Our results extend the knowledge of immune responses of pulmonate snails to trematode invasion and promote P. corneus as a new model for the study of molluscan defence reactions.

One-Step Soaking Strategy toward Anti-Swelling Hydrogels with a Stiff "Armor"

Double-network (DN) hydrogels consisting of noncovalent interacting networks are highly desired due to their well-controlled compositions and environmental friendliness, but the low water resistance always impairs their mechanical strength. Here, an anti-swelling hydrogel possessing the core/shell architecture through rational regulation of multiple weak noncovalent interactions is prepared. A composite hydrogel consists of chitosan (CS) and poly(N-acryloyl 2-glycine) (PACG), readily forming the shell-structured DN hydrogel after soaking in a FeCl₃ solution because of in situ formation of chain entanglements, hydrogen bonds, and ionic coordination. The produced DN hydrogels exhibit excellent anti-swelling behaviors and mechanical durability for over half a year, even in some strict situations. Taking the merits of noncovalent bonds in adjustability and reversibility, the swelling property of these hydrogels can be easily customized through control of the ion species and concentrations. A dynamically reversible transition from super-swelling to anti-swelling is realized by breaking up and rebuilding the metal-coordination complexes. This facile but efficient strategy of turning the noncovalent interactions and consequently the mechanics and anti-swelling properties is imperative to achieve the rational design of high-performance hydrogels with specific usage requirements and expand their applicability to a higher stage.

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.

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.

Endosymbiont population genomics sheds light on transmission mode, partner specificity, and stability of the scaly-foot snail holobiont

The scaly-foot snail (Chrysomallon squamiferum) inhabiting deep-sea hydrothermal vents in the Indian Ocean relies on its sulphur-oxidising gammaproteobacterial endosymbionts for nutrition and energy. In this study, we investigate the specificity, transmission mode, and stability of multiple scaly-foot snail populations dwelling in five vent fields with considerably disparate geological, physical and chemical environmental conditions. Results of population genomics analyses reveal an incongruent phylogeny between the endosymbiont and mitochondrial genomes of the scaly-foot snails in the five vent fields sampled, indicating that the hosts obtain endosymbionts via horizontal transmission in each generation. However, the genetic homogeneity of many symbiont populations implies that vertical transmission cannot be ruled out either. Fluorescence in situ hybridisation of ovarian tissue yields symbiont signals around the oocytes, suggesting that vertical transmission co-occurs with horizontal transmission. Results of in situ environmental measurements and gene expression analyses from in situ fixed samples show that the snail host buffers the differences in environmental conditions to provide the endosymbionts with a stable intracellular micro-environment, where the symbionts serve key metabolic functions and benefit from the host’s cushion. The mixed transmission mode, symbiont specificity at the species level, and stable intracellular environment provided by the host support the evolutionary, ecological, and physiological success of scaly-foot snail holobionts in different vents with unique environmental parameters.

Metabolism Interactions Promote the Overall Functioning of the Episymbiotic Chemosynthetic Community of Shinkaia crosnieri of Cold Seeps

Remarkably diverse bacteria have been observed as biofilm aggregates on the surface of deep-sea invertebrates that support the growth of hosts through chemosynthetic carbon fixation. Growing evidence also indicates that community-wide interactions, and especially cooperation among symbionts, contribute to overall community productivity. Here, metagenome-guided metatranscriptomic and metabolic analyses were conducted to investigate the taxonomic composition, functions, and potential interactions of symbionts dwelling on the seta of Shinkaia crosnieri lobsters in a methane cold seep. Methylococcales and Thiotrichales dominated the community, followed by the Campylobacteriales, Nitrosococcales, Flavobacteriales, and Chitinophagales Metabolic interactions may be common among the episymbionts since many separate taxon genomes encoded complementary genes within metabolic pathways. Specifically, Thiotrichales could contribute to detoxification of hydroxylamine that is a metabolic by-product of Methylococcales. Further, Nitrosococcales may rely on methanol leaked from Methylococcales cells that efficiently oxidize methane. Elemental sulfur may also serve as a community good that enhances sulfur utilization that benefits the overall community, as evidenced by confocal Raman microscopy. Stable intermediates may connect symbiont metabolic activities in cyclical oxic-hypoxic fluctuating environments, which then enhance overall community functioning. This hypothesis was partially confirmed via in situ experiments. These results highlight the importance of microbe-microbe interactions in symbiosis and deep-sea adaptation.

IMPORTANCE Symbioses between chemosynthetic bacteria and marine invertebrates are common in deep-sea chemosynthetic ecosystems and are considered critical foundations for deep-sea colonization. Episymbiotic microorganisms tend to form condensed biofilms that may facilitate metabolite sharing among biofilm populations. However, the prevalence of metabolic interactions among deep-sea episymbionts and their contributions to deep-sea adaptations are not well understood due to sampling and cultivation difficulties associated with deep-sea environments. Here, we investigated metabolic interactions among the episymbionts of Shinkaia crosnieri, a dominant chemosynthetic ecosystem lobster species in the Northwest Pacific Ocean. Meta-omics characterizations were conducted alongside in situ experiments to validate interaction hypotheses. Furthermore, imaging analysis was conducted, including electron microscopy, fluorescent in situ hybridization (FISH), and confocal Raman microscopy (CRM), to provide direct evidence of metabolic interactions. The results support the Black Queen Hypothesis, wherein leaked public goods are shared among cohabitating microorganisms to enhance the overall adaptability of the community via cooperation.

Dynamic evolutionary history of spiralian-specific TALE homeobox genes in mollusks

Homeobox genes play essential roles in the early development of many animals. Although the repertoire of most homeobox genes, including three amino acid loop extension (TALE)-type homeobox genes, is conserved in animals, spiralian-TALE (SPILE) genes are a notable exception. In this study, SPILE genes were extracted from the genomic data of 22 mollusk species and classified into four clades (-A/C, -B, -D, and -E) to determine which SPILE genes exhibit dynamic repertoire changes. While SPILE-D and -E duplications were rarely observed, SPILE-B duplication was observed in the bivalve lineage and SPILE-A/C duplication was observed in multiple clades. Conversely, most or all SPILE genes were lost in cephalopods and in some gastropod lineages. SPILE gene expression patterns were also analyzed in multiple mollusk species using publicly available RNA-seq data. The majority of SPILE genes examined, particularly those in the A/C- and B-clades, were specifically expressed during early development, suggesting that most SPILE genes exert specific roles in early development. This comprehensive cataloging and characterization revealed a dynamic evolutionary history, including SPILE-A/C and -B gene duplications and the loss of SPILE genes in several lineages. Furthermore, this study provides a useful resource for studying the molecular mechanism of spiralian early development and the evolution of young and lineage-specific transcription factors.

How many single-copy orthologous genes from whole genomes reveal deep gastropod relationships?

The Gastropoda contains 80% of existing mollusks and is the most diverse animal class second only to the Insecta. However, the deep phylogeny of gastropods has been controversial for a long time. Especially the position of Patellogastropoda is a major uncertainty. Morphology and some mitochondria studies concluded that Patellogastropoda is likely to be sister to all other gastropods (Orthogastropoda hypothesis), while transcriptomic and other mitogenomic studies indicated that Patellogastropoda and Vetigastropoda are sister taxa (Psilogastropoda). With the release of high-quality genomes, orthologous genes can be better identified and serve as powerful candidates for phylogenetic analysis. The question is, given the current limitations on the taxon sampling side, how many markers are needed to provide robust results. Here, we identified single-copy orthologous genes (SOGs) from 14 gastropods species with whole genomes available which cover five main gastropod subclasses. We generated different datasets from 395 to 1610 SOGs by allowing species missing in different levels. We constructed gene trees of each SOG, and inferred species trees from different collections of gene trees. We found as the number of SOGs increased, the inferred topology changed from Patellogastropoda being sister to all other gastropods to Patellogastropoda being sister to Vetigastropoda + Neomphalina (Psilogastropoda s.l.), with considerable support. Our study thus rejects the Orthogastropoda concept showing that the selection of the representative species and use of sufficient informative sites greatly influence the analysis of deep gastropod phylogeny.

Identification of peptidoglycan recognition proteins in hemocytes and kidney of common periwinkle Littorinalittorea

Peptidoglycan Recognition Proteins (PGRPs) are a diverse group of proteins involved in innate immunity. In particular, PGRPs have been shown to participate in immune pattern recognition in various mollusks. However, they have not been described in Caenogastropoda, a large molluscan group comprising sea, freshwater and land snails. In this study, four short PGRPs with molecular weights ranging from 21 to 34 kDa and their isoforms were identified and structurally characterized in the kidney and hemocytic transcriptomes of a caenogastropod mollusk Littorina littorea. All of them (LlPGRP1-4) are secretory, possess a signal peptide and a characteristic N-terminal N-acetylmuramoyl-l-alanine amidase (Ami) domain with conserved Zn²⁺ binding- and amidase catalytic sites. The shortest proteins, LlPGRP1 and LlPGRP2, have no additional conserved motifs on the N-terminus. In longer and most abundantly expressed LlPGRP3 and LlPGRP4 the Ami-domain is combined with an N-terminal SH3-domain and a cysteine-rich motif, respectively. Expression analysis showed that LlPGRPs of the common periwinkle were uninvolved in the immune response to infection with trematode Himasthla elongata though they might act in antibacterial defense.

Duplication of spiralian-specific TALE genes and evolution of the blastomere specification mechanism in the bivalve lineage

Background

Despite the conserved pattern of the cell-fate map among spiralians, bivalves display several modified characteristics during their early development, including early specification of the D blastomere by the cytoplasmic content, as well as the distinctive fate of the 2d blastomere. However, it is unclear what changes in gene regulatory mechanisms led to such changes in cell specification patterns. Spiralian-TALE (SPILE) genes are a group of spiralian-specific transcription factors that play a role in specifying blastomere cell fates during early development in limpets. We hypothesised that the expansion of SPILE gene repertoires influenced the evolution of the specification pattern of blastomere cell fates.

Results

We performed a transcriptome analysis of early development in the purplish bifurcate mussel and identified 13 SPILE genes. Phylogenetic analysis of the SPILE gene in molluscs suggested that duplications of SPILE genes occurred in the bivalve lineage. We examined the expression patterns of the SPILE gene in mussels and found that some SPILE genes were expressed in quartet-specific patterns, as observed in limpets. Furthermore, we found that several SPILE genes that had undergone gene duplication were specifically expressed in the D quadrant, C and D quadrants or the 2d blastomere. These expression patterns were distinct from the expression patterns of SPILE in their limpet counterparts.

Conclusions

These results suggest that, in addition to their ancestral role in quartet specification, certain SPILE genes in mussels contribute to the specification of the C and D quadrants. We suggest that the expansion of SPILE genes in the bivalve lineage contributed to the evolution of a unique cell fate specification pattern in bivalves.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13227-021-00181-2.

Arms race in a cell: genomic, transcriptomic, and proteomic insights into intracellular phage-bacteria interplay in deep-sea snail holobionts

BACKGROUND

Deep-sea animals in hydrothermal vents often form endosymbioses with chemosynthetic bacteria. Endosymbionts serve essential biochemical and ecological functions, but the prokaryotic viruses (phages) that determine their fate are unknown.

RESULTS

We conducted metagenomic analysis of a deep-sea vent snail. We assembled four genome bins for Caudovirales phages that had developed dual endosymbiosis with sulphur-oxidising bacteria (SOB) and methane-oxidising bacteria (MOB). Clustered regularly interspaced short palindromic repeat (CRISPR) spacer mapping, genome comparison, and transcriptomic profiling revealed that phages Bin1, Bin2, and Bin4 infected SOB and MOB. The observation of prophages in the snail endosymbionts and expression of the phage integrase gene suggested the presence of lysogenic infection, and the expression of phage structural protein and lysozyme genes indicated active lytic infection. Furthermore, SOB and MOB appear to employ adaptive CRISPR-Cas systems to target phage DNA. Additional expressed defence systems, such as innate restriction-modification systems and dormancy-inducing toxin-antitoxin systems, may co-function and form multiple lines for anti-viral defence. To counter host defence, phages Bin1, Bin2, and Bin3 appear to have evolved anti-restriction mechanisms and expressed methyltransferase genes that potentially counterbalance host restriction activity. In addition, the high-level expression of the auxiliary metabolic genes narGH, which encode nitrate reductase subunits, may promote ATP production, thereby benefiting phage DNA packaging for replication.

CONCLUSIONS

This study provides new insights into phage-bacteria interplay in intracellular environments of a deep-sea vent snail. Video Abstract.

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.