Insights into deep-sea adaptations and host-symbiont interactions: A comparative transcriptome study on Bathymodiolus mussels and their coastal relatives

Distinct interspecies thermal resistance strategies exhibited by euplanktonic, tychoplanktonic and benthic diatoms under marine heatwaves

Extreme climate events, such as marine heatwaves (MHWs), are expected to occur more frequently and intensely in the future, resulting in a substantial impact on marine life. The way that diatoms respond to MHWs may have crucial effects on global primary production and biogeochemical cycles. Euplanktonic diatoms appear to benefit from MHWs directly, but this phenomenon needs an explanation. As concerns tychoplanktonic and benthic diatoms, no studies have been addressed on their thermal response strategies. To address this, we investigated the responses and underlying mechanisms of three typical growth forms of diatoms, Pseudo-nitzschia multiseries (euplanktonic), Paralia guyana (tychoplanktonic) and Navicula avium (benthic), under heat stress by combining a growth experiment with transcriptomic analysis. Our results showed that the physiological responses of diatoms to MHWs and underlying molecular mechanisms are largely related to their growth forms. The euplanktonic diatom was first depressed, but then had a distinct increase in the growth rate accompanied by inducing zeatin and unsaturated fatty acid biosynthesis and repressing substance assimilation and energy metabolism. Contrarily, the benthic diatom showed elevated substance and energy demands for macromolecules accumulation by reducing cell division and increasing photosynthesis and nitrogen assimilation. The tychoplanktonic diatom exhibited higher physiological plasticity to maintain growth and cellular homeostasis. Our results indicate the increased rate of cell division in euplanktonic diatoms under heat stress is likely an emergency response strategy promoting diatom dispersal for survival, but at the cost of disturbances of metabolic balance.

Shallow-water mussels (Mytilus galloprovincialis) adapt to deep-sea environment through transcriptomic and metagenomic insights

Recent studies have unveiled the deep sea as a rich biosphere, populated by species descended from shallow-water ancestors post-mass extinctions. Research on genomic evolution and microbial symbiosis has shed light on how these species thrive in extreme deep-sea conditions. However, early adaptation stages, particularly the roles of conserved genes and symbiotic microbes, remain inadequately understood. This study examined transcriptomic and microbiome changes in shallow-water mussels Mytilus galloprovincialis exposed to deep-sea conditions at the Site-F cold seep in the South China Sea. Results reveal complex gene expression adjustments in stress response, immune defense, homeostasis, and energy metabolism pathways during adaptation. After 10 days of deep-sea exposure, shallow-water mussels and their microbial communities closely resembled those of native deep-sea mussels, demonstrating host and microbiome convergence in response to adaptive shifts. Notably, methanotrophic bacteria, key symbionts in native deep-sea mussels, emerged as a dominant group in the exposed mussels. Host genes involved in immune recognition and endocytosis correlated significantly with the abundance of these bacteria. Overall, our analyses provide insights into adaptive transcriptional regulation and microbiome dynamics of mussels in deep-sea environments, highlighting the roles of conserved genes and microbial community shifts in adapting to extreme environments.

Transcriptomic responses and evolutionary insights of deep-sea and shallow-water mussels under high hydrostatic pressure condition

Marine mussels inhabit a wide range of ocean depths, necessitating unique adaptations to cope with varying hydrostatic pressures. This study investigates the transcriptomic responses and evolutionary adaptations of the deep-sea mussel Gigantidas platifrons and the shallow-water mussel Mytilus galloprovincialis to high hydrostatic pressure (HHP) conditions. By exposing atmospheric pressure (AP) acclimated G. platifrons and M. galloprovincialis to HHP, we aim to simulate extreme environmental challenges and assess their adaptive mechanisms. Through comparative transcriptomic analysis, we identified both conserved and species-specific mechanisms of adaptation, with a notable change in gene expression associated with immune system, substance transport, protein ubiquitination, apoptosis, lipid metabolism and antioxidant processes in both species. G. platifrons demonstrated an augmented lipid metabolism, whereas M. galloprovincialis exhibited a dampened immune function. Additionally, the expressed pattern of deep-sea mussel G. platifrons were more consistent than shallow-water mussel M. galloprovincialis under hydrostatic pressures changed conditions which corresponding the long-term living stable deep-sea environment. Moreover, evolutionary analysis pinpointed positively selected genes in G. platifrons that are linked to transmembrane transporters, DNA repair and replication, apoptosis, ubiquitination which are important to cell structural integrity, substances transport, and cellular growth regulation. This indicates a specialized adaptation strategy in G. platifrons to cope with the persistent HHP conditions of the deep sea. These results offer significant insights into the molecular underpinnings of mussel adaptation to varied hydrostatic conditions and enhance our comprehension of the evolutionary forces driving their depth-specific adaptations.

Chromosome-level genome assembly of a deep-sea Venus flytrap sea anemone sheds light upon adaptations to an extremely oligotrophic environment

The Venus flytrap sea anemone Actinoscyphia liui inhabits the nutrient-limited deep ocean in the tropical western Pacific. Compared with most other sea anemones, it has undergone a distinct modification of body shape similar to that of the botanic flytrap. However, the molecular mechanism by which such a peculiar sea anemone adapts to a deep-sea oligotrophic environment is unknown. Here, we report the chromosomal-level genome of A. liui constructed from PacBio and Hi-C data. The assembled genome is 522 Mb in size and exhibits a continuous scaffold N50 of 58.4 Mb. Different from most other sea anemones, which typically possess 14-18 chromosomes per haplotype, A. liui has only 11. The reduced number of chromosomes is associated with chromosome fusion, which likely represents an adaptive strategy to economize energy in oligotrophic deep-sea environments. Comparative analysis with other deep-sea sea anemones revealed adaptive evolution in genes related to cellular autophagy (TMBIM6, SESN1, SCOCB and RPTOR) and mitochondrial energy metabolism (MDH1B and KAD2), which may aid in A. liui coping with severe food scarcity. Meanwhile, the genome has undergone at least two rounds of expansion in gene families associated with fast synaptic transmission, facilitating rapid responses to water currents and prey. Positive selection was detected on putative phosphorylation sites of muscle contraction-related proteins, possibly further improving feeding efficiency. Overall, the present study provides insights into the molecular adaptation to deep-sea oligotrophic environments and sheds light upon potential effects of a novel morphology on the evolution of Cnidaria.

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.

Different parts of the mussel Gigantidas haimaensis holobiont responded differently to deep-sea sampling stress

Acute environmental changes cause stress during conventional deep-sea biological sampling without in situ fixation and affect gene expressions of samples collected. However, the degree of influence and underlying mechanisms are hardly investigated. Here, we conducted comparative transcriptomic analyses between in situ and onboard fixed gills and between in situ and onboard fixed mantles of deep-sea mussel Gigantidas haimaensis to assess the effects of incidental sampling stress. Results showed that transcription, translation, and energy metabolism were upregulated in onboard fixed gills and mantles, thereby mobilizing rapid gene expression to tackle the stress. Autophagy and phagocytosis that related to symbiotic interactions between the host and endosymbiont were downregulated in the onboard fixed gills. These findings demonstrated that symbiotic gill and nonsymbiotic mantle responded differently to sampling stress, and symbiosis in the gill was perturbed. Further comparative metatranscriptomic analysis between in situ and onboard fixed gills revealed that stress response genes, peptidoglycan biosynthesis, and methane fixation were upregulated in the onboard fixed endosymbiotic Gammaproteobacteria inside the gills, implying that energy metabolism of the endosymbiont was increased to cope with sampling stress. Furthermore, comparative analysis between the mussel G. haimaensis and the limpet Bathyacmaea lactea transcriptomes resultedidentified six transcription factor orthologs upregulated in both onboard fixed mussel mantles and limpets, including sharply increased early growth response protein 1 and Kruppel-like factor 5. They potentially play key roles in initiating the response of sampled deep-sea macrobenthos to sampling stress. Our results clearly show that in situ fixed biological samples are vital for studying deep-sea environmental adaptation.

Toxicology assessment of deep-sea mining impacts on Gigantidas platifrons: A comparative in situ and laboratory metal exposure study

Deep-sea toxicology is essential for deep-sea environmental impact assessment. Yet most toxicology experiments are conducted solely in laboratory settings, overlooking the complexities of the deep-sea environment. Here we carried out metal exposure experiments in both the laboratory and in situ, to compare and evaluate the response patterns of Gigantidas platifrons to metal exposure (copper [Cu] or cadmium [Cd] at 100 μg/L for 48 h). Metal concentrations, traditional biochemical parameters, and fatty acid composition were assessed in deep-sea mussel gills. The results revealed significant metal accumulation in deep-sea mussel gills in both laboratory and in situ experiments. Metal exposure could induce oxidative stress, neurotoxicity, an immune response, altered energy metabolism, and changes to fatty acid composition in mussel gills. Interestingly, the metal accumulating capability, biochemical response patterns, and fatty acid composition each varied under differing experimental systems. In the laboratory setting, Cd-exposed mussels exhibited a higher value for integrated biomarker response (IBR) while in situ the Cu-exposed mussels instead displayed a higher IBR value. This study emphasizes the importance of performing deep-sea toxicology experiments in situ and contributes valuable data to a standardized workflow for deep-sea toxicology assessment.

Comparative study of integrated bio-responses in deep-sea and nearshore mussels upon abiotic condition changes: Insight into distinct regulation and adaptation

Deep-sea mussels, one of the dominant species in most deep-sea ecosystems, have long been used as model organisms to investigate the adaptations and symbiotic relationships of deep-sea macrofauna under laboratory conditions due to their ability to survive under atmospheric pressure. However, the impact of additional abiotic conditions beyond pressure, such as temperature and light, on their physiological characteristics remains unknown. In this study, deep-sea mussels (Gigantidas platifrons) from cold seep of the South China Sea, along with nearshore mussels (Mytilus coruscus) from the East China Sea, were reared in unfavorable abiotic conditions for up to 8 days. Integrated biochemical indexes including antioxidant defense, immune ability and energy metabolism were investigated in the gill and digestive gland, while cytotoxicity was determined in hemocytes of both types of mussels. The results revealed mild bio-responses in two types of mussels in the laboratory, represented by the effective antioxidant defense with constant total antioxidant capability level and malondialdehyde content. There were also disparate adaptations in deep-sea and nearshore mussels. In deep-sea mussels, significantly increased immune response and energy reservation were observed in gills, together with the elevated cytotoxicity in hemocytes, implying the more severe biological adaptation was required, mainly due to the symbiotic bacteria loss under laboratory conditions. On the contrary, insignificant biological responses were exhibited in nearshore mussels except for the increased energy consumption, indicating the trade-off strategy to use more energy to deal with potential stress. Overall, this comparative study highlights the basal bio-responses of deep-sea and nearshore mussels out of their native environments, providing evidence that short-term culture of both mussels under easily achievable laboratory conditions would not dramatically alter their biological status. This finding will assist in broadening the application of deep-sea mussels as model organism in future research regardless of the specialized research equipment.

To live or die: "Fine-tuning" adaptation revealed by systemic analyses in symbiotic bathymodiolin mussels from diverse deep-sea extreme ecosystems

Hydrothermal vents (HVs) and cold seeps (CSs) are typical deep-sea extreme ecosystems with their own geochemical characteristics to supply the unique living conditions for local communities. Once HVs or CSs stop emission, the dramatic environmental change would pose survival risks to deep-sea organisms. Up to now, limited knowledge has been available to understand the biological responses and adaptive strategy to the extreme environments and their transition from active to extinct stage, mainly due to the technical difficulties and lack of representative organisms. In this study, bathymodiolin mussels, the dominant and successful species surviving in diverse deep-sea extreme ecosystems, were collected from active and extinct HVs (Southwest Indian Ocean) or CSs (South China Sea) via two individual cruises. The transcriptomic analysis and determination of multiple biological indexes in stress defense and metabolic systems were conducted in both gills and digestive glands of mussels, together with the metagenomic analysis of symbionts in mussels. The results revealed the ecosystem- and tissue-specific transcriptional regulation in mussels, addressing the autologous adaptations in antioxidant defense, energy utilization and key compounds (i.e. sulfur) metabolism. In detail, the successful antioxidant defense contributed to conquering the oxidative stress induced during the unavoidable metabolism of xenobiotics commonly existing in the extreme ecosystems; changes in metabolic rate functioned to handle toxic matters in different surroundings; upregulated gene expression of sulfide:quinone oxidoreductase indicated an active sulfide detoxification in mussels from HVs and active stage of HVs & CSs. Coordinately, a heterologous adaptation, characterized by the functional compensation between symbionts and mussels in energy utilization, sulfur and carbon metabolism, was also evidenced by the bacterial metagenomic analysis. Taken together, a new insight was proposed that symbiotic bathymodiolin mussels would develop a "finetuning" strategy combining the autologous and heterologous regulations to fulfill the efficient and effective adaptations for successful survival.

mTORC1 regulates phagosome digestion of symbiotic bacteria for intracellular nutritional symbiosis in a deep-sea mussel

Phagocytosis is one of the methods used to acquire symbiotic bacteria to establish intracellular symbiosis. A deep-sea mussel, Bathymodiolus japonicus, acquires its symbiont from the environment by phagocytosis of gill epithelial cells and receives nutrients from them. However, the manner by which mussels retain the symbiont without phagosome digestion remains unknown. Here, we show that controlling the mechanistic target of rapamycin complex 1 (mTORC1) in mussels leads to retaining symbionts in gill cells. The symbiont is essential for the host mussel nutrition; however, depleting the symbiont's energy source triggers the phagosome digestion of symbionts. Meanwhile, the inhibition of mTORC1 by rapamycin prevented the digestion of the resident symbionts and of the engulfed exogenous dead symbionts in gill cells. This indicates that mTORC1 promotes phagosome digestion of symbionts under reduced nutrient supply from the symbiont. The regulation mechanism of phagosome digestion by mTORC1 through nutrient signaling with symbionts is key for maintaining animal-microbe intracellular nutritional symbiosis.

Metabolic handoffs between multiple symbionts may benefit the deep-sea bathymodioline mussels

Bathymodioline mussels rely on thiotrophic and/or methanotrophic chemosynthetic symbionts for nutrition, yet, secondary heterotrophic symbionts are often present and play an unknown role in the fitness of the organism. The bathymodioline Idas mussels that thrive in gas seeps and on sunken wood in the Mediterranean Sea and the Atlantic Ocean, host at least six symbiont lineages that often co-occur. These lineages include the primary symbionts chemosynthetic methane- and sulfur-oxidizing gammaproteobacteria, and the secondary symbionts, Methylophagaceae, Nitrincolaceae and Flavobacteriaceae, whose physiology and metabolism are obscure. Little is known about if and how these symbionts interact or exchange metabolites. Here we curated metagenome-assembled genomes of Idas modiolaeformis symbionts and used genome-centered metatranscriptomics and metaproteomics to assess key symbiont functions. The Methylophagaceae symbiont is a methylotrophic autotroph, as it encoded and expressed the ribulose monophosphate and Calvin-Benson-Bassham cycle enzymes, particularly RuBisCO. The Nitrincolaceae ASP10-02a symbiont likely fuels its metabolism with nitrogen-rich macromolecules and may provide the holobiont with vitamin B12. The Urechidicola (Flavobacteriaceae) symbionts likely degrade glycans and may remove NO. Our findings indicate that these flexible associations allow for expanding the range of substrates and environmental niches, via new metabolic functions and handoffs.

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.

Interactions among deep-sea mussels and their epibiotic and endosymbiotic chemoautotrophic bacteria: Insights from multi-omics analysis

Endosymbiosis with Gammaproteobacteria is fundamental for the success of bathymodioline mussels in deep-sea chemosynthesis-based ecosystems. However, the recent discovery of Campylobacteria on the gill surfaces of these mussels suggests that these host-bacterial relationships may be more complex than previously thought. Using the cold-seep mussel ( Gigantidas haimaensis) as a model, we explored this host-bacterial system by assembling the host transcriptome and genomes of its epibiotic Campylobacteria and endosymbiotic Gammaproteobacteria and quantifying their gene and protein expression levels. We found that the epibiont applies a sulfur oxidizing (SOX) multienzyme complex with the acquisition of soxB from Gammaproteobacteria for energy production and switched from a reductive tricarboxylic acid (rTCA) cycle to a Calvin-Benson-Bassham (CBB) cycle for carbon assimilation. The host provides metabolic intermediates, inorganic carbon, and thiosulfate to satisfy the materials and energy requirements of the epibiont, but whether the epibiont benefits the host is unclear. The endosymbiont adopts methane oxidation and the ribulose monophosphate pathway (RuMP) for energy production, providing the major source of energy for itself and the host. The host obtains most of its nutrients, such as lysine, glutamine, valine, isoleucine, leucine, histidine, and folate, from the endosymbiont. In addition, host pattern recognition receptors, including toll-like receptors, peptidoglycan recognition proteins, and C-type lectins, may participate in bacterial infection, maintenance, and population regulation. Overall, this study provides insights into the complex host-bacterial relationships that have enabled mussels and bacteria to thrive in deep-sea chemosynthetic ecosystems.

Rhythms and Clocks in Marine Organisms

The regular movements of waves and tides are obvious representations of the oceans' rhythmicity. But the rhythms of marine life span across ecological niches and timescales, including short (in the range of hours) and long (in the range of days and months) periods. These rhythms regulate the physiology and behavior of individuals, as well as their interactions with each other and with the environment. This review highlights examples of rhythmicity in marine animals and algae that represent important groups of marine life across different habitats. The examples cover ecologically highly relevant species and a growing number of laboratory model systems that are used to disentangle key mechanistic principles. The review introduces fundamental concepts of chronobiology, such as the distinction between rhythmic and endogenous oscillator-driven processes. It also addresses the relevance of studying diverse rhythms and oscillators, as well as their interconnection, for making better predictions of how species will respond to environmental perturbations, including climate change. As the review aims to address scientists from the diverse fields of marine biology, ecology, and molecular chronobiology, all of which have their own scientific terms, we provide definitions of key terms throughout the article.

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.

Insights into symbiotic interactions from metatranscriptome analysis of deep-sea mussel Gigantidas platifrons under long-term laboratory maintenance

Symbioses between invertebrates and chemosynthetic bacteria are of fundamental importance in deep-sea ecosystems, but the mechanisms that enable their symbiont associations are still largely undescribed, owing to the culturable difficulties of deep-sea lives. Bathymodiolinae mussels are remarkable in their ability to overcome decompression and can be maintained successfully for an extended period under atmospheric pressure, thus providing a model for investigating the molecular basis of symbiotic interactions. Herein, we conducted metatranscriptome sequencing and gene co-expression network analysis of Gigantidas platifrons under laboratory maintenance with gradual loss of symbionts. The results revealed that one-day short-term maintenance triggered global transcriptional perturbation in symbionts, but little gene expression changes in mussel hosts, which were mainly involved in responses to environmental changes. Long-term maintenance with depleted symbionts induced a metabolic shift in the mussel host. The most notable changes were the suppression of sterol biosynthesis and the complementary activation of terpenoid backbone synthesis in response to the reduction of bacteria-derived terpenoid sources. In addition, we detected the upregulation of host proteasomes responsible for amino acid deprivation caused by symbiont depletion. Additionally, a significant correlation between host microtubule motor activity and symbiont abundance was revealed, suggesting the possible function of microtubule-based intracellular trafficking in the nutritional interaction of symbiosis. Overall, by analyzing the dynamic transcriptomic changes during the loss of symbionts, our study highlights the nutritional importance of symbionts in supplementing terpenoid compounds and essential amino acids and provides insight into the molecular mechanisms and strategies underlying the symbiotic interactions in deep-sea ecosystems.

Global patterns in symbiont selection and transmission strategies in sponges

Sponges host dense and diverse communities of microbes (known as the microbiome) beneficial for the host nutrition and defense. Symbionts in turn receive shelter and metabolites from the sponge host, making their relationship beneficial for both partners. Given that sponge-microbes associations are fundamental for the survival of both, especially the sponge, such relationship is maintained through their life and even passed on to the future generations. In many organisms, the microbiome has profound effects on the development of the host, but the influence of the microbiome on the reproductive and developmental pathways of the sponges are less understood. In sponges, microbes are passed on to oocytes, sperm, embryos, and larvae (known as vertical transmission), using a variety of methods that include direct uptake from the mesohyl through phagocytosis by oocytes to indirect transmission to the oocyte by nurse cells. Such microbes can remain in the reproductive elements untouched, for transfer to offspring, or can be digested to make the yolky nutrient reserves of oocytes and larvae. When and how those decisions are made are fundamentally unanswered questions in sponge reproduction. Here we review the diversity of vertical transmission modes existent in the entire phylum Porifera through detailed imaging using electron microscopy, available metabarcoding data from reproductive elements, and macroevolutionary patterns associated to phylogenetic constraints. Additionally, we examine the fidelity of this vertical transmission and possible reasons for the observed variability in some developmental stages. Our current understanding in marine sponges, however, is that the adult microbial community is established by a combination of both vertical and horizontal (acquisition from the surrounding environment in each new generation) transmission processes, although the extent in which each mode shapes the adult microbiome still remains to be determined. We also assessed the fundamental role of filtration, the cellular structures for acquiring external microbes, and the role of the host immune system, that ultimately shapes the stable communities of prokaryotes observed in adult sponges.

Adaptation and evolution of the sea anemone Alvinactis sp. to deep-sea hydrothermal vents: A comparison using transcriptomes

Sea anemones are diverse and ecologically successful members of Anthozoa. They are often found in intertidal and shallow waters, although a few of them inhabit harsher living conditions, such as deep-sea hydrothermal vents. Here, we sequenced the transcriptome of the vent sea anemone Alvinactis sp., which was collected from Edmond vent along the central Indian Ocean ridge at a depth of 3275 m, to explore the molecular mechanisms related to adaptation to vents. Compared with another deep-sea anemone (Paraphelliactis xishaensis) and five shallow water sea anemones, a total of 117 positively selected genes and 46 significantly expanded gene families were found in Alvinactis sp. specifically that may be related to its vent-specific aspect of adaptation. In addition, 127 positively selected genes and 23 significantly expanded gene families that were found in both Alvinactis sp. and P. xishaensis. Among these, vent-specific adaptations of Alvinactis sp. may involve genetic alterations in peroxisome, ubiquitin-mediated protein degradation, oxidative phosphorylation, and endocytosis, and its deep-sea adaptation may involve changes in genetic information processing. Differentially expressed genes between Alvinactis sp. and the deep-sea anemone P. xishaensis were enriched in a variety of pathways related to adaptation, such as energy metabolism, genetic information processing, endocytosis, and peroxisomes. Overall, we provided the first transcriptome of sea anemones that inhabit vents, which enriches our knowledge of deep-sea hydrothermal vent adaptation and the diversity of sea anemones.

Full-Length Transcriptome Comparison Provides Novel Insights into the Molecular Basis of Adaptation to Different Ecological Niches of the Deep-Sea Hydrothermal Vent in Alvinocaridid Shrimps

The deep-sea hydrothermal vent ecosystem is one of the extreme chemoautotrophic environments. Shinkaicaris leurokolos Kikuchi and Hashimoto, 2000, and Alvinocaris longirostris Kikuchi and Ohta, 1995, are typically co-distributed and closely related alvinocaridid shrimps in hydrothermal vent areas with different ecological niches, providing an excellent model for studying the adaptive evolution mechanism of animals in the extreme deep-sea hydrothermal vent environment. The shrimp S. leurokolos lives in close proximity to the chimney vent discharging high-temperature fluid, while A. longirostris inhabits the peripheral areas of hydrothermal vents. In this study, full-length transcriptomes of S. leurokolos and A. longirostris were generated using a combination of single-molecule real-time (SMRT) and Illumina RNA-seq technology. Expression analyses of the transcriptomes showed that among the top 30% of highly expressed genes of each species, more genes related to sulfide and heavy metal metabolism (sulfide: quinone oxidoreductase, SQR; persulfide dioxygenase, ETHE1; thiosulfate sulfurtransferase, TST, and ferritin, FRI) were specifically highly expressed in S. leurokolos, while genes involved in maintaining epibiotic bacteria or pathogen resistance (beta-1,3-glucan-binding protein, BGBP; endochitinase, CHIT; acidic mammalian chitinase, CHIA, and anti-lipopolysaccharide factors, ALPS) were highly expressed in A. longirostris. Gene family expansion analysis revealed that genes related to anti-oxidant metabolism (cytosolic manganese superoxide dismutase, SODM; glutathione S-transferase, GST, and glutathione peroxidase, GPX) and heat stress (heat shock cognate 70 kDa protein, HSP70 and heat shock 70 kDa protein cognate 4, HSP7D) underwent significant expansion in S. leurokolos, while CHIA and CHIT involved in pathogen resistance significantly expanded in A. longirostris. Finally, 66 positively selected genes (PSGs) were identified in the vent shrimp S. leurokolos. Most of the PSGs were involved in DNA repair, antioxidation, immune defense, and heat stress response, suggesting their function in the adaptive evolution of species inhabiting the extreme vent microhabitat. This study provides abundant genetic resources for deep-sea invertebrates, and is expected to lay the foundation for deep decipherment of the adaptive evolution mechanism of shrimps in a deep-sea chemosynthetic ecosystem based on further whole-genome comparison.

Gene expression profiles provide insights into the survival strategies in deep-sea mussel (Bathymodiolus platifrons) of different developmental stages

Background

Deep-sea mussels living in the cold seeps with enormous biomass act as the primary consumers. They are well adapted to the extreme environment where light is absent, and hydrogen sulfide, methane, and other hydrocarbon-rich fluid seepage occur. Despite previous studies on diversity, role, evolution, and symbiosis, the changing adaptation patterns during different developmental stages of the deep-sea mussels remain largely unknown.

Results

The deep-sea mussels (Bathymodiolus platifrons) of two developmental stages were collected from the cold seep during the ocean voyage. The gills, mantles, and adductor muscles of these mussels were used for the Illumina sequencing. A total of 135 Gb data were obtained, and subsequently, 46,376 unigenes were generated using de-novo assembly strategy. According to the gene expression analysis, amounts of genes were most actively expressed in the gills, especially genes involved in environmental information processing. Genes encoding Toll-like receptors and sulfate transporters were up-regulated in gills, indicating that the gill acts as both intermedium and protective screen in the deep-sea mussel. Lysosomal enzymes and solute carrier responsible for nutrients absorption were up-regulated in the older mussel, while genes related to toxin resistance and autophagy were up-regulated in the younger one, suggesting that the older mussel might be in a vigorous stage while the younger mussel was still paying efforts in survival and adaptation.

Conclusions

In general, our study suggested that the adaptation capacity might be formed gradually during the development of deep-sea mussels, in which the gill and the symbionts play essential roles.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08505-9.

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

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

microRNAs facilitate comprehensive responses of Bathymodiolinae mussel against symbiotic and nonsymbiotic bacteria stimulation

Bathymodiolinae mussels are dominant species in cold seeps and hydrothermal vents and could harbor endosymbionts in gill bacteriocytes. However, mechanisms underlying the symbiosis have remained largely undisclosed for years. In the present study, the global expression pattern of immune-related genes and miRNAs were surveyed in Gigantidas platifrons during bacterial challenges using enriched symbiotic methane oxidation bacteria MOBs or nonsymbiotic Vibrio. As a result, multiple pattern recognition receptors were found differentially expressed at 12 h and 24 h post bacteria challenges and distinctly clustered between stimulations. Dozens of immune effectors along with signal transducers were also modulated simultaneously during MOB or Vibrio challenge. A total of 459 miRNAs were identified in the gill while some were differentially expressed post MOB or nonsymbiotic bacteria challenge. A variety of immune-related genes were annotated as target genes of aforesaid differentially expressed miRNAs. As a result, biological processes including the immune recognition, lysosome activity and bacteria engulfment were suggested to be dynamically modulated by miRNAs in either symbiotic or nonsymbiotic bacteria challenge. It was suggested that G. platifrons mussels could maintain a robust immune response against invading pathogens while establishing symbiosis with chemosynthetic bacteria with the orchestra of immune-related genes and miRNAs.

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.

A Toll-like receptor identified in Gigantidas platifrons and its potential role in the immune recognition of endosymbiotic methane oxidation bacteria

Symbiosis with chemosynthetic bacteria is an important ecological strategy for the deep-sea megafaunas including mollusks, tubeworms and crustacean to obtain nutrients in hydrothermal vents and cold seeps. How the megafaunas recognize symbionts and establish the symbiosis has attracted much attention. Bathymodiolinae mussels are endemic species in both hydrothermal vents and cold seeps while the immune recognition mechanism underlying the symbiosis is not well understood due to the nonculturable symbionts. In previous study, a lipopolysaccharide (LPS) pull-down assay was conducted in Gigantidas platifrons to screen the pattern recognition receptors potentially involved in the recognition of symbiotic methane-oxidizing bacteria (MOB). Consequently, a total of 208 proteins including GpTLR13 were identified. Here the molecular structure, expression pattern and immune function of GpTLR13 were further analyzed. It was found that GpTLR13 could bind intensively with the lipid A structure of LPS through surface plasmon resonance analysis. The expression alternations of GpTLR13 transcripts during a 28-day of symbiont-depletion assay were investigated by real-time qPCR. As a result, a robust decrease of GpTLR13 transcripts was observed accompanying with the loss of symbionts, implying its participation in symbiosis. In addition, GpTLR13 transcripts were found expressed exclusively in the bacteriocytes of gills of G. platifrons by in situ hybridization. It was therefore speculated that GpTLR13 may be involved in the immune recognition of symbiotic methane-oxidizing bacteria by specifically recognizing the lipid A structure of LPS. However, the interaction between GpTLR13 and symbiotic MOB was failed to be addressed due to the nonculturable symbionts. Nevertheless, the present result has provided with a promising candidate as well as a new approach for the identification of symbiont-related genes in Bathymodiolinae mussels.

A glimpse of deep-sea adaptation in chemosynthetic holobionts: Depressurization causes DNA fragmentation and cell death of methanotrophic endosymbionts rather than their deep-sea Bathymodiolinae host

Bathymodiolinae mussels are typical species in deep-sea cold seeps and hydrothermal vents and an ideal model for investigating chemosynthetic symbiosis and the influence of high hydrostatic pressure on deep-sea organisms. Herein, the potential influence of depressurization on DNA fragmentation and cell death in Bathymodiolinae hosts and their methanotrophic symbionts were surveyed using isobaric and unpressurized samples. As a hallmark of cell death, massive DNA fragmentation was observed in methanotrophic symbionts from unpressurized Bathymodiolinae while several endonucleases and restriction enzymes were upregulated. Additionally, genes involved in DNA repair, glucose/methane metabolism as well as two-component regulatory system were also differentially expressed in depressurized symbionts. DNA fragmentation and programmed cell death, however, were rarely detected in the host bacteriocytes owing to the orchestrated upregulation of inhibitor of apoptosis genes and downregulation of caspase genes. Meanwhile, diverse host immune recognition receptors were promoted during depressurization, probably enabling the regain of symbionts. When the holobionts were subjected to a prolonged acclimation at atmospheric pressure, alternations in both the DNA fragmentation and the expression atlas of aforesaid genes were continuously observed in symbionts, demonstrating the persistent influence of depressurization. Contrarily, the host cells demonstrated certain tolerance against depressurization stress as expression level of some immune-related genes returned to the basal level in isobaric samples. Altogether, the present study illustrates the distinct stress responses of Bathymodiolinae hosts and their methanotrophic symbionts against depressurization, which could provide further insight into the deep-sea adaptation of Bathymodiolinae holobionts while highlighting the necessity of using isobaric sampling methods in deep-sea research.

Molecular analyses of the gill symbiosis of the bathymodiolin mussel Gigantidas platifrons

Although the deep-sea bathymodiolin mussels have been intensively studied as a model of animal-bacteria symbiosis, it remains challenging to assess the host-symbiont interactions due to the complexity of the symbiotic tissue-the gill. Using cold-seep mussel Gigantidas platifrons as a model, we isolated the symbiont harboring bacteriocytes and profiled the transcriptomes of the three major parts of the symbiosis-the gill, the bacteriocyte, and the symbiont. This breakdown of the complex symbiotic tissue allowed us to characterize the host-symbiont interactions further. Our data showed that the gill's non-symbiotic parts play crucial roles in maintaining and protecting the symbiosis; the bacteriocytes supply the symbiont with metabolites, control symbiont population, and shelter the symbiont from phage infection; the symbiont dedicates to the methane oxidation and energy production. This study demonstrates that the bathymodiolin symbiosis interacts at the tissue, cellular, and molecular level, maintaining high efficiency and harmonic chemosynthetic micro niche.

The cytokine MIF controls daily rhythms of symbiont nutrition in an animal-bacterial association

The recent recognition that many symbioses exhibit daily rhythms has encouraged research into the partner dialogue that drives these biological oscillations. Here we characterized the pivotal role of the versatile cytokine macrophage migration inhibitory factor (MIF) in regulating a metabolic rhythm in the model light-organ symbiosis between Euprymna scolopes and Vibrio fischeri As the juvenile host matures, it develops complex daily rhythms characterized by profound changes in the association, from gene expression to behavior. One such rhythm is a diurnal shift in symbiont metabolism triggered by the periodic provision of a specific nutrient by the mature host: each night the symbionts catabolize chitin released from hemocytes (phagocytic immune cells) that traffic into the light-organ crypts, where the population of V. fischeri cells resides. Nocturnal migration of these macrophage-like cells, together with identification of an E. scolopes MIF (EsMIF) in the light-organ transcriptome, led us to ask whether EsMIF might be the gatekeeper controlling the periodic movement of the hemocytes. Western blots, ELISAs, and confocal immunocytochemistry showed EsMIF was at highest abundance in the light organ. Its concentration there was lowest at night, when hemocytes entered the crypts. EsMIF inhibited migration of isolated hemocytes, whereas exported bacterial products, including peptidoglycan derivatives and secreted chitin catabolites, induced migration. These results provide evidence that the nocturnal decrease in EsMIF concentration permits the hemocytes to be drawn into the crypts, delivering chitin. This nutritional function for a cytokine offers the basis for the diurnal rhythms underlying a dynamic symbiotic conversation.

The Diversity, Composition, and Putative Functions of Gill-Associated Bacteria of Bathymodiolin Mussel and Vesicomyid Clam from Haima Cold Seep, South China Sea

The Haima cold seep, which is one of the two active cold seeps in the South China Sea, is known for its great ecological importance. The seep bivalves are assumed to depend mainly on their bacterial symbiosis for survival and growth. However, information on the bacterial diversity, composition, and putative function of gill-associated of dominant dwelling animals in Haima cold seep remain elusive. Herein, we adopted a high-throughput sequencing of 16S rRNA gene amplicons, and function prediction methods (Functional Annotation of Prokaryotic Taxa (FAPROTAX) and Phylogenetic Investigation of Communities by Reconstruction of Unobserved States (PICURUSTs)) to purposely illustrate the taxonomic and phylogenetic diversity, composition, and putative functions of the symbionts in bathymodiolin mussel Gigantidas haimaensis (Bivalvia: Mytilidae: Gigantidas) and vesicomyid clam Archivesica marissinica (Bivalvia: Glossoidea: Vesicomyidae). The predominant microbes of both species were Proteobacteria and Gammaproteobacteria on the phylum and class level, respectively. The taxonomic and phylogenetic diversity of gill microbial communities in G. haimaensis were significantly different from those in A. marissinica (p < 0.05). Nine functional groups, including seven carbon-related biogeochemical groups, were identified through the FAPROTAX analysis. However, the most dominant groups for G. haimaensis and A. marissinica were both chemoheterotrophic. G. haimaensis and A. marissinica shared many pathways, however, 16 obtained Kyoto Encyclopedia of Genes and Genomes (KEGG) orthologous groups (42.11%) significantly differed between the two species (p < 0.05). These findings would provide insight into the functions of microbes in the element cycling and energy flow as well as the host-symbiont relationship of bivalves in the Haima cold seep environment.

Molecular Response to High Hydrostatic Pressure: Time-Series Transcriptomic Analysis of Shallow-Water Sea Cucumber Apostichopus japonicus

Hydrostatic pressure is a key environmental factor constraining the benthic migration of shallow-water invertebrates. Although many studies have examined the physiological effects of high hydrostatic pressure on shallow-water invertebrates, the molecular response to high pressure is not fully understood. This question has received increasing attention because ocean warming is forcing the bathymetric migrations of shallow-water invertebrates. Here, we applied time-series transcriptomic analysis to high-pressure incubated and atmospheric pressure-recovered shallow-water sea cucumber (Apostichopus japonicus) to address this question. A total of 44 samples from 15 experimental groups were sequenced. Our results showed that most genes responded to pressure stress at the beginning when pressure was changed, but significant differences of gene expression appeared after 4 to 6 h. Transcription was the most sensitive biological process responding to high-pressure exposure, which was enriched among up-regulated genes after 2 h, followed by ubiquitination (4 h), endocytosis (6 h), stress response (6 h), methylation regulation (24 h), and transmembrane transportation (24 h). After high-pressure incubation, all these biological processes remained up-regulated within 4–6 h at atmospheric pressure. Overall, our results revealed the dynamic transcriptional response of A. japonicus to high-pressure exposure. Additionally, few quantitative or functional responses related to A. japonicus on transcriptional level were introduced by hydrostatic pressure changes after 1 h, and main biological responses were introduced after 4 h, suggesting that, when hydrostatic pressure is the mainly changed environmental factor, it will be better to fix sea cucumber samples for transcriptomic analysis within 1 h, but 4 h will be also acceptable.

Comparative transcriptomic analysis of deep- and shallow-water barnacle species (Cirripedia, Poecilasmatidae) provides insights into deep-sea adaptation of sessile crustaceans

BACKGROUND

Barnacles are specialized marine organisms that differ from other crustaceans in possession of a calcareous shell, which is attached to submerged surfaces. Barnacles have a wide distribution, mostly in the intertidal zone and shallow waters, but a few species inhabit the deep-sea floor. It is of interest to investigate how such sessile crustaceans became adapted to extreme deep-sea environments. We sequenced the transcriptomes of a deep-sea barnacle, Glyptelasma gigas collected at a depth of 731 m from the northern area of the Zhongjiannan Basin, and a shallow-water coordinal relative, Octolasmis warwicki. The purpose of this study was to provide genetic resources for investigating adaptation mechanisms of deep-sea barnacles.

RESULTS

Totals of 62,470 and 51,585 unigenes were assembled for G. gigas and O. warwicki, respectively, and functional annotation of these unigenes was made using public databases. Comparison of the protein-coding genes between the deep- and shallow-water barnacles, and with those of four other shallow-water crustaceans, revealed 26 gene families that had experienced significant expansion in G. gigas. Functional annotation showed that these expanded genes were predominately related to DNA repair, signal transduction and carbohydrate metabolism. Base substitution analysis on the 11,611 single-copy orthologs between G. gigas and O. warwicki indicated that 25 of them were distinctly positive selected in the deep-sea barnacle, including genes related to transcription, DNA repair, ligand binding, ion channels and energy metabolism, potentially indicating their importance for survival of G. gigas in the deep-sea environment.

CONCLUSIONS

The barnacle G. gigas has adopted strategies of expansion of specific gene families and of positive selection of key genes to counteract the negative effects of high hydrostatic pressure, hypoxia, low temperature and food limitation on the deep-sea floor. These expanded gene families and genes under positive selection would tend to enhance the capacities of G. gigas for signal transduction, genetic information processing and energy metabolism, and facilitate networks for perceiving and responding physiologically to the environmental conditions in deep-sea habitats. In short, our results provide genomic evidence relating to deep-sea adaptation of G. gigas, which provide a basis for further biological studies of sessile crustaceans in the deep sea.

Metal adaptation strategies of deep-sea Bathymodiolus mussels from a cold seep and three hydrothermal vents in the West Pacific

Deep-sea Bathymodiolus mussels are ubiquitous in most cold seeps and hydrothermal fields, where they have adapted to various toxic environments including high metal exposure. However, there is scarce knowledge of metal accumulation and metal-related biomarkers in B. mussels. Here, we present data for metal concentrations (Ag, Cd, Cr, Cu, Fe, Mn, Pb, and Zn) and metal related biomarkers (superoxide dismutase-SOD, catalase-CAT, glutathione peroxidase-GPX, glutathione-GSH, metallothioneins-MTs, and lipid peroxidation-LPO) in different tissues of B. mussels from four different deep-sea geochemical settings: one cold seep and three vent fields in the West Pacific Ocean. Results showed that mussel gills generally exhibited higher metal enrichment than the mantle. Mussels from hydrothermal vents usually had higher metal concentrations (Fe, Cr, Cd, and Pb) than those from cold seep, which could be related to their higher contents in fluids or sediments. However, despite quite different metals loads among the geochemical environment settings, Mn, Zn, and Cu concentrations varied over a smaller range across the sampling sites, implying biological regulation by deep-sea mussels for these elements. Several statistically significant correlations were observed between SOD, CAT, GSH, MTs, and metal levels in analyzed tissues. Although the vent ecosystem is harsher than the cold seep ecosystem, according to our results their mussels' biomarker levels were not so different. This finding suggests that some adaptive or compensatory mechanisms may occur in chronically polluted deep-sea mussels. Principal component analysis allowed for distinguishing different deep-sea settings, indicating that B. mussels are robust indicators of their living environments. We also compared our results with those reported for coastal mussels. To our best knowledge, this is the first integrated study to report metal accumulation and metal-related biomarkers in the deep-sea B. mussels from the West Pacific.

Phylogeny and highland adaptation of Chinese species in Allium section Daghestanica (Amaryllidaceae) revealed by transcriptome sequencing

Allium L. is one of the largest monocotyledonous genera with extensive distribution in the Northern Hemisphere. The fundamental phylogenies of Allium have been investigated using many morphological and molecular characters. However, the morphological characters may not agree with the molecular results in some Allium groups or sections (such as the Chinese Allium section Daghestanica), which may result in ambiguous species relationships and hinder further evolutionary and adaptive researches. Here, transcriptome sequences of the six Chinese endemics from Allium section Daghestanica were collected, with their single-copy genes (SCGs) were extracted. The interspecies relationships were analyzed using concatenation and coalescent methods. The branch-site model (BSM) was conducted to detect the positively selected genes (PSGs) in five highland species of this section. Based on 1644, 1281 and 1580 SCGs in flowers, leaves, and flowers-leaves combination respectively, a robust consistent and well-resolved phylogeny was generated from the concatenation method. Strong conflicts among individual gene trees were detected in the coalescent method, and morphological characters were incongruent with molecular relationships to some degree. Many PSGs were involved in responses of various stresses and stimuli (e.g. hypoxia, low temperature, aridity), DNA repair, metabolism, nutrient or energy intake, photosynthesis, and signal transduction. Our study revealed a clear interspecies relationship of Chinese endemics in Allium section Daghestanica and suggested that the discordance between morphological characters and molecular relationships might result from that the former are more susceptible to convergence compared with the latter. PSGs detected in our study may provide some insights into highland adaptation in Allium species.

Insights into high-pressure acclimation: comparative transcriptome analysis of sea cucumber Apostichopus japonicus at different hydrostatic pressure exposures

Background

Global climate change is predicted to force the bathymetric migrations of shallow-water marine invertebrates. Hydrostatic pressure is proposed to be one of the major environmental factors limiting the vertical distribution of extant marine invertebrates. However, the high-pressure acclimation mechanisms are not yet fully understood.

Results

In this study, the shallow-water sea cucumber Apostichopus japonicus was incubated at 15 and 25 MPa at 15 °C for 24 h, and subjected to comparative transcriptome analysis. Nine samples were sequenced and assembled into 553,507 unigenes with a N50 length of 1204 bp. Three groups of differentially expressed genes (DEGs) were identified according to their gene expression patterns, including 38 linearly related DEGs whose expression patterns were linearly correlated with hydrostatic pressure, 244 pressure-sensitive DEGs which were up-regulated at both 15 and 25 MPa, and 257 high-pressure-induced DEGs which were up-regulated at 25 MPa but not up-regulated at 15 MPa.

Conclusions

Our results indicated that the genes and biological processes involving high-pressure acclimation are similar to those related to deep-sea adaptation. In addition to representative biological processes involving deep-sea adaptation (such as antioxidation, immune response, genetic information processing, and DNA repair), two biological processes, namely, ubiquitination and endocytosis, which can collaborate with each other and regulate the elimination of misfolded proteins, also responded to high-pressure exposure in our study. The up-regulation of these two processes suggested that high hydrostatic pressure would lead to the increase of misfolded protein synthesis, and this may result in the death of shallow-water sea cucumber under high-pressure exposure.

Identification and characterization of endosymbiosis-related immune genes in deep-sea mussels Gigantidas platifrons

Deep-sea mussels of the subfamily Bathymodiolinae are common and numerically dominant species widely distributed in cold seeps and hydrothermal vents. During long-time evolution, deep-sea mussels have evolved to be well adapted to the local environment of cold seeps and hydrothermal vents by various ways, especially by establishing endosymbiosis with chemotrophic bacteria. However, biological processes underlying the establishment and maintenance of symbiosis between host mussels and symbionts are largely unclear. In the present study, Gigantidas platifrons genes possibly involved in the symbiosis with methane oxidation symbionts were identified and characterized by Lipopolysaccharide (LPS) pull-down and in situ hybridization. Five immune related proteins including Toll-like receptor 2 (TLR2), integrin, vacuolar sorting protein (VSP), matrix metalloproteinase 1 (MMP1), and leucine-rich repeat (LRR-1) were identified by LPS pull-down assay. These five proteins were all conserved in either molecular sequences or functional domains and known to be key molecules in host immune recognition, phagocytosis, and lysosome-mediated digestion. Furthermore, in situ hybridization of LRR-1, TLR2 and VSP genes was conducted to investigate their expression patterns in gill tissues of G. platifrons. Consequently, LRR-1, TLR2, and VSP genes were found expressed exclusively in the bacteriocytes of G. platifrons. Therefore, it was suggested that TLR2, integrin, VSP, MMP1, and LRR-1 might be crucial molecules in the symbiosis between G. platifrons and methane oxidation bacteria by participating in symbiosis-related immune processes.

Organ transcriptomes of the lucinid clam Loripes orbiculatus (Poli, 1791) provide insights into their specialised roles in the biology of a chemosymbiotic bivalve

BACKGROUND

The lucinid clam Loripes orbiculatus lives in a nutritional symbiosis with sulphur-oxidizing bacteria housed in its gills. Although our understanding of the lucinid endosymbiont physiology and metabolism has made significant progress, relatively little is known about how the host regulates the symbiosis at the genetic and molecular levels. We generated transcriptomes from four L. orbiculatus organs (gills, foot, visceral mass, and mantle) for differential expression analyses, to better understand this clam's physiological adaptations to a chemosymbiotic lifestyle, and how it regulates nutritional and immune interactions with its symbionts.

RESULTS

The transcriptome profile of the symbiont-housing gill suggests the regulation of apoptosis and innate immunity are important processes in this organ. We also identified many transcripts encoding ion transporters from the solute carrier family that possibly allow metabolite exchange between host and symbiont. Despite the clam holobiont's clear reliance on chemosynthesis, the clam's visceral mass, which contains the digestive tract, is characterised by enzymes involved in digestion, carbohydrate recognition and metabolism, suggesting that L. orbiculatus has a mixotrophic diet. The foot transcriptome is dominated by the biosynthesis of glycoproteins for the construction of mucus tubes, and receptors that mediate the detection of chemical cues in the environment.

CONCLUSIONS

The transcriptome profiles of gills, mantle, foot and visceral mass provide insights into the molecular basis underlying the functional specialisation of bivalve organs adapted to a chemosymbiotic lifestyle.

An LRR-domain containing protein identified in Bathymodiolus platifrons serves as intracellular recognition receptor for the endosymbiotic methane-oxidation bacteria

As domain species in seep and vent ecosystem, Bathymodioline mussels has been regarded as a model organism in investigating deep sea chemosymbiosis. However, mechanisms underlying their symbiosis with chemosynthetic bacteria, especially how the host recognizes symbionts, have remained largely unsolved. In the present study, a modified pull-down assay was conducted using enriched symbiotic methane-oxidation bacteria as bait and gill proteins of Bathymodiolus platifrons as a target to isolate pattern recognition receptors involved in the immune recognition of symbionts. As a result, a total of 47 proteins including BpLRR-1 were identified from the pull-down assay. It was found that complete cDNA sequence of BpLRR-1 contained an open reading frame of 1479 bp and could encode a protein of 492 amino acid residues with no signal peptide or transmembrane region but eight LRR motif and two EFh motif. The binding patterns of BpLRR-1 against microbial associated molecular patterns were subsequently investigated by surface plasmon resonance analysis and LPS pull-down assay. Consequently, BpLRR-1 was found with high binding affinity with LPS and suggested as a key molecule in recognizing symbionts. Besides, transcripts of BpLRR-1 were found decreased significantly during symbiont depletion assay yet increased rigorously during symbionts or nonsymbiotic Vibrio alginolyticus challenge, further demonstrating its participation in the chemosynthetic symbiosis. Collectively, these results suggest that BpLRR-1 could serve as an intracellular recognition receptor for the endosymbionts, providing new hints for understanding the immune recognition in symbiosis of B. platifrons.

The complete mitochondrial genome of Calyptogena marissinica (Heterodonta: Veneroida: Vesicomyidae): Insight into the deep-sea adaptive evolution of vesicomyids

The deep-sea chemosynthetic environment is one of the most extreme environments on the Earth, with low oxygen, high hydrostatic pressure and high levels of toxic substances. Species of the family Vesicomyidae are among the dominant chemosymbiotic bivalves found in this harsh habitat. Mitochondria play a vital role in oxygen usage and energy metabolism; thus, they may be under selection during the adaptive evolution of deep-sea vesicomyids. In this study, the mitochondrial genome (mitogenome) of the vesicomyid bivalve Calyptogena marissinica was sequenced with Illumina sequencing. The mitogenome of C. marissinica is 17,374 bp in length and contains 13 protein-coding genes, 2 ribosomal RNA genes (rrnS and rrnL) and 22 transfer RNA genes. All of these genes are encoded on the heavy strand. Some special elements, such as tandem repeat sequences, “G(A)_nT” motifs and AT-rich sequences, were observed in the control region of the C. marissinica mitogenome, which is involved in the regulation of replication and transcription of the mitogenome and may be helpful in adjusting the mitochondrial energy metabolism of organisms to adapt to the deep-sea chemosynthetic environment. The gene arrangement of protein-coding genes was identical to that of other sequenced vesicomyids. Phylogenetic analyses clustered C. marissinica with previously reported vesicomyid bivalves with high support values. Positive selection analysis revealed evidence of adaptive change in the mitogenome of Vesicomyidae. Ten potentially important adaptive residues were identified, which were located in cox1, cox3, cob, nad2, nad4 and nad5. Overall, this study sheds light on the mitogenomic adaptation of vesicomyid bivalves that inhabit the deep-sea chemosynthetic environment.

Transcriptomic analysis reveals insights into deep-sea adaptations of the dominant species, Shinkaia crosnieri (Crustacea: Decapoda: Anomura), inhabiting both hydrothermal vents and cold seeps

Background

Hydrothermal vents and cold seeps are typical deep-sea chemosynthetically-driven ecosystems that allow high abundance of specialized macro-benthos. To gather knowledge about the genetic basis of adaptation to these extreme environments, species shared between different habitats, especially for the dominant species, are of particular interest. The galatheid squat lobster, Shinkaia crosnieri Baba and Williams, 1998, is one of the few dominant species inhabiting both deep-sea hydrothermal vents and cold seeps. In this study, we performed transcriptome analyses of S. crosnieri collected from the Iheya North hydrothermal vent (HV) and a cold seep in the South China Sea (CS) to provide insights into how this species has evolved to thrive in different deep-sea chemosynthetic ecosystems.

Results

We analyzed 5347 orthologs between HV and CS to identify genes under positive selection through the maximum likelihood approach. A total of 82 genes were identified to be positively selected and covered diverse functional categories, potentially indicating their importance for S. crosnieri to cope with environmental heterogeneity between deep-sea vents and seeps. Among 39,806 annotated unigenes, a large number of differentially expressed genes (DEGs) were identified between HV and CS, including 339 and 206 genes significantly up-regulated in HV and CS, respectively. Most of the DEGs associated with stress response and immunity were up-regulated in HV, possibly allowing S. crosnieri to increase its capability to manage more environmental stresses in the hydrothermal vents.

Conclusions

We provide the first comprehensive transcriptomic resource for the deep-sea squat lobster, S. crosnieri, inhabiting both hydrothermal vents and cold seeps. A number of stress response and immune-related genes were positively selected and/or differentially expressed, potentially indicating their important roles for S. crosnieri to thrive in both deep-sea vents and cold seeps. Our results indicated that genetic adaptation of S. crosnieri to different deep-sea chemosynthetic environments might be mediated by adaptive evolution of functional genes related to stress response and immunity, and alterations in their gene expression that lead to different stress resistance. However, further work is required to test these proposed hypotheses. All results can constitute important baseline data for further studies towards elucidating the adaptive mechanisms in deep-sea crustaceans.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5753-7) contains supplementary material, which is available to authorized users.

Comparative transcriptome analysis of Eogammarus possjeticus at different hydrostatic pressure and temperature exposures

Hydrostatic pressure is an important environmental factor affecting the vertical distribution of marine organisms. Laboratory-based studies have shown that many extant shallow-water marine benthic invertebrates can tolerate hydrostatic pressure outside their known natural distributions. However, only a few studies have focused on the molecular mechanisms of pressure acclimatisation. In the present work, we examined the pressure tolerance of the shallow-water amphipod Eogammarus possjeticus at various temperatures (5, 10, 15, and 20 °C) and hydrostatic pressures (0.1–30 MPa) for 16 h. Six of these experimental groups were used for transcriptome analysis. We found that 100% of E. possjeticus survived under 20 MPa at all temperature conditions for 16 h. Sequence assembly resulted in 138, 304 unigenes. Results of differential expression analysis revealed that 94 well-annotated genes were up-regulated under high pressure. All these findings indicated that the pressure tolerance of E. possjeticus was related to temperature. Several biological processes including energy metabolism, antioxidation, immunity, lipid metabolism, membrane-related process, genetic information processing, and DNA repair are probably involved in the acclimatisation in deep-sea environments.

Gill symbionts of the cold-seep mussel Bathymodiolus platifrons: Composition, environmental dependency and immune control

Deep-sea Bathymodiolus mussels depend on the organic carbon supplied by symbionts inside their gills. In this study, optimized methods of quantitative real-time PCR and fluorescence in situ hybridization targeted to both mRNA and 16S rRNA were used to investigate the gill symbionts of the cold-seep mussel Bathymodiolus platifrons, including species composition, environmental dependency and immune control by the host. Our results showed that methanotrophs were the major symbiotic bacteria in the gills of B. platifrons, while thiotrophs were scarce. In the mussels freshly collected from the deep sea, methanotrophs were housed in bacteriocytes in a unique circular pattern, and a lysosome-related gene (VAMP) encoding a vesicle-associated membrane protein was expressed at a high level and presented exactly where the methanotrophs occurred. After the mussels were reared for three months in aquaria without methane supply, the abundance of methanotrophs decreased significantly and their circle-shaped distribution pattern disappeared; in addition, the expression of VAMP decreased significantly. These results suggest that the symbiosis between B. platifrons and methanotrophs is influenced by the environment and that the lysosomal system plays an important immune role in controlling the abundance of endosymbionts in host. This study provides a reliable method for investigating symbionts in deep-sea mussels and enriches the knowledge about symbionts in B. platifrons.

High rates of apoptosis visualized in the symbiont-bearing gills of deep-sea Bathymodiolus mussels

Symbiosis between Bathymodiolus and Gammaproteobacteria allows these deep-sea mussels to live in toxic environments such as hydrothermal vents and cold seeps. The quantity of endosymbionts within the gill-bacteriocytes appears to vary according to the hosts environment; however, the mechanisms of endosymbiont population size regulation remain obscure. We investigated the possibility of a control of endosymbiont density by apoptosis, a programmed cell death, in three mussel species. Fluorometric TUNEL and active Caspase-3-targeting antibodies were used to visualize and quantify apoptotic cells in mussel gills. To control for potential artefacts due to depressurization upon specimen recovery from the deep-sea, the apoptotic rates between mussels recovered unpressurised, versus mussels recovered in a pressure-maintaining device, were compared in two species from hydrothermal vents on the Mid-Atlantic Ridge: Bathymodiolus azoricus and B. puteoserpentis. Results show that pressurized recovery had no significant effect on the apoptotic rate in the gill filaments. Apoptotic levels were highest in the ciliated zone and in the circulating hemocytes, compared to the bacteriocyte zone. Apoptotic gill-cells in B. aff. boomerang from cold seeps off the Gulf of Guinea show similar distribution patterns. Deep-sea symbiotic mussels have much higher rates of apoptosis in their gills than the coastal mussel Mytilus edulis, which lacks chemolithoautotrophic symbionts. We discuss how apoptosis might be one of the mechanisms that contribute to the adaptation of deep-sea mussels to toxic environments and/or to symbiosis.

Gonadal Transcriptome Analysis of Pacific Abalone Haliotis discus discus: Identification of Genes Involved in Germ Cell Development

Little is known about the molecular mechanisms governing gonadal developmental processes in abalones. Here, we conducted transcriptome analysis of Pacific abalone Haliotis discus discus for gene discovery in the brain, ovary, testis, and unfertilized eggs. Among the annotated unigenes, 48.6% of unigenes were identified by Venn diagram analysis as having universal or tissue-specific expression. Twenty-three genes with gonad-biased gene ontology (GO) terms were first obtained. Secondly, 36 genes were found by screening known gene names related to germ cell development. Finally, 17 genes were obtained by querying the annotated unigene database for zygotically expressed gonadal genes (ovary and testis) and maternally expressed gonadal genes (ovary, testis, and unfertilized eggs) using keywords related to reproduction. To further verify tissue distribution pattern and subcellular localization of these genes, RT-PCR and in situ hybridization were performed using a unigene encoding a germ cell marker, vasa, as control. The results showed that vasa was expressed mainly in the early developmental stages of germ cells in both sexes. One of the candidate genes, vitelline envelope zona pellucida domain protein 12 (ZP12), was expressed in the primordial germ cells of immature gonad and early developmental stages of germ cells of the adult female. The results obtained from the present study suggest that vasa and ZP12 are involved in germ cell development of Pacific abalone and that ZP12 is an especially useful germ cell-specific marker in immature adults. The current gonadal transcriptome profile is an extensive resource for future reproductive molecular biology studies of this species.

First comprehensive analysis of lysine acetylation in Alvinocaris longirostris from the deep-sea hydrothermal vents

Background

Deep-sea hydrothermal vents are unique chemoautotrophic ecosystems with harsh conditions. Alvinocaris longirostris is one of the dominant crustacean species inhabiting in these extreme environments. It is significant to clarify mechanisms in their adaptation to the vents. Lysine acetylation has been known to play critical roles in the regulation of many cellular processes. However, its function in A. longirostris and even marine invertebrates remains elusive. Our study is the first, to our knowledge, to comprehensively investigate lysine acetylome in A. longirostris.

Results

In total, 501 unique acetylation sites from 206 proteins were identified by combination of affinity enrichment and high-sensitive-massspectrometer. It was revealed that Arg, His and Lys occurred most frequently at the + 1 position downstream of the acetylation sites, which were all alkaline amino acids and positively charged. Functional analysis revealed that the protein acetylation was involved in diverse cellular processes, such as biosynthesis of amino acids, citrate cycle, fatty acid degradation and oxidative phosphorylation. Acetylated proteins were found enriched in mitochondrion and peroxisome, and many stress response related proteins were also discovered to be acetylated, like arginine kinases, heat shock protein 70, and hemocyanins. In the two hemocyanins, nine acetylation sites were identified, among which one acetylation site was unique in A. longirostris when compared with other shallow water shrimps. Further studies are warranted to verify its function.

Conclusion

The lysine acetylome of A. longirostris is investigated for the first time and brings new insights into the regulation function of the lysine acetylation. The results supply abundant resources for exploring the functions of acetylation in A. longirostris and other shrimps.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4745-3) contains supplementary material, which is available to authorized users.

Bias and Correction in RNA-seq Data for Marine Species

RNA-seq is a recently developed approach widely used for transcriptome profiling in biological analyses that use next-generation sequencing technologies. Accurate estimation of gene expression levels is critical for answering biological questions. Here, we show that the commonly used measure of gene expression levels, fragments per kilobase of transcript per million mapped reads (FPKM), is biased in transcript length, GC content, and dinucleotide frequencies in the RNA-seq analysis of marine species. We used a generalized linear model to correct the observed biases of FPKM. We used RNA-seq data sets from eight species obtained by different sequencing methods to evaluate the correction methods. Our work contributes to the understanding of potential technical artifacts in RNA-seq experiments for marine species, and presents a means by which more accurate gene expression measures can be obtained.