Conjugating effects of symbionts and environmental factors on gene expression in deep-sea hydrothermal vent mussels

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.

Genome Study of α-, β-, and γ-Carbonic Anhydrases from the Thermophilic Microbiome of Marine Hydrothermal Vent Ecosystems

Carbonic anhydrases (CAs) are metalloenzymes that can help organisms survive in hydrothermal vents by hydrating carbon dioxide (CO₂). In this study, we focus on alpha (α), beta (β), and gamma (γ) CAs, which are present in the thermophilic microbiome of marine hydrothermal vents. The coding genes of these enzymes can be transferred between hydrothermal-vent organisms via horizontal gene transfer (HGT), which is an important tool in natural biodiversity. We performed big data mining and bioinformatics studies on α-, β-, and γ-CA coding genes from the thermophilic microbiome of marine hydrothermal vents. The results showed a reasonable association between thermostable α-, β-, and γ-CAs in the microbial population of the hydrothermal vents. This relationship could be due to HGT. We found evidence of HGT of α- and β-CAs between Cycloclasticus sp., a symbiont of Bathymodiolus heckerae, and an endosymbiont of Riftia pachyptila via Integrons. Conversely, HGT of β-CA genes from the endosymbiont Tevnia jerichonana to the endosymbiont Riftia pachyptila was detected. In addition, Hydrogenovibrio crunogenus SP-41 contains a β-CA gene on genomic islands (GIs). This gene can be transferred by HGT to Hydrogenovibrio sp. MA2-6, a methanotrophic endosymbiont of Bathymodiolus azoricus, and a methanotrophic endosymbiont of Bathymodiolus puteoserpentis. The endosymbiont of R. pachyptila has a γ-CA gene in the genome. If α- and β-CA coding genes have been derived from other microorganisms, such as endosymbionts of T. jerichonana and Cycloclasticus sp. as the endosymbiont of B. heckerae, through HGT, the theory of the necessity of thermostable CA enzymes for survival in the extreme ecosystem of hydrothermal vents is suggested and helps the conservation of microbiome natural diversity in hydrothermal vents. These harsh ecosystems, with their integral players, such as HGT and endosymbionts, significantly impact the enrichment of life on Earth and the carbon cycle in the ocean.

Identification and characterization of a symbiotic agglutination-related C-type lectin from the hydrothermal vent shrimp Rimicaris exoculata

Rimicaris exoculata (Decapoda: Bresiliidae) is one of the dominant species of hydrothermal vent communities, which inside its gill chamber harbors ectosymbioses with taxonomic invariability while compositional flexibility. Several studies have revealed that the establishment of symbiosis can be initiated and selected by innate immunity-related pattern recognition receptors (PRRs), such as C-type lectins (CTLs). In this research, a CTL was identified in R. exoculata (termed RCTL), which showed high expression at both mRNA and protein levels in the scaphognathite, an organ where the ectosymbionts are attached outside its setae. Linear correlationships were observed between the relative quantities of two major symbionts and the expression of RCTL based on analyzing different shrimp individuals. The recombinant protein of RCTL could recognize and agglutinate the cultivable γ-proteobacterium of Escherichia coli in a Ca²⁺-dependent manner, obeying a dose-dependent and time-cumulative pattern. Unlike conventional crustacean CTLs, the involvement of RCTL could not affect the bacterial growth, which is a key issue for the successful establishment of symbiosis. These results implied that RCTL might play a critical role in symbiotic recognition and attachment to R. exoculata. It also provides insights to understand how R. exoculata adapted to such a chemosynthesis-based environment.

Global host molecular perturbations upon in situ loss of bacterial endosymbionts in the deep-sea mussel Bathymodiolus azoricus assessed using proteomics and transcriptomics

Background

Colonization of deep-sea hydrothermal vents by most invertebrates was made efficient through their adaptation to a symbiotic lifestyle with chemosynthetic bacteria, the primary producers in these ecosystems. Anatomical adaptations such as the establishment of specialized cells or organs have been evidenced in numerous deep-sea invertebrates. However, very few studies detailed global inter-dependencies between host and symbionts in these ecosystems. In this study, we proposed to describe, using a proteo-transcriptomic approach, the effects of symbionts loss on the deep-sea mussel Bathymodiolus azoricus’ molecular biology. We induced an in situ depletion of symbionts and compared the proteo-transcriptome of the gills of mussels in three conditions: symbiotic mussels (natural population), symbiont-depleted mussels and aposymbiotic mussels.

Results

Global proteomic and transcriptomic results evidenced a global disruption of host machinery in aposymbiotic organisms. We observed that the total number of proteins identified decreased from 1118 in symbiotic mussels to 790 in partially depleted mussels and 761 in aposymbiotic mussels. Using microarrays we identified 4300 transcripts differentially expressed between symbiont-depleted and symbiotic mussels. Among these transcripts, 799 were found differentially expressed in aposymbiotic mussels and almost twice as many in symbiont-depleted mussels as compared to symbiotic mussels. Regarding apoptotic and immune system processes – known to be largely involved in symbiotic interactions – an overall up-regulation of associated proteins and transcripts was observed in symbiont-depleted mussels.

Conclusion

Overall, our study showed a global impairment of host machinery and an activation of both the immune and apoptotic system following symbiont-depletion. One of the main assumptions is the involvement of symbiotic bacteria in the inhibition and regulation of immune and apoptotic systems. As such, symbiotic bacteria may increase their lifespan in gill cells while managing the defense of the holobiont against putative pathogens.

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.

Symbiont type and environmental factors affect transcriptome-wide gene expression in the coral Montipora capitata

Reef-building corals may harbor genetically distinct lineages of endosymbiotic dinoflagellates in the genus Symbiodinium, which have been shown to affect important colony properties, including growth rates and resilience against environmental stress. However, the molecular processes underlying these differences are not well understood. In this study, we used whole transcriptome sequencing (RNA-seq) to assess gene expression differences between 27 samples of the coral Montipora capitata predominantly hosting two different Symbiodinium types in clades C and D. The samples were further characterized by their origin from two field sites on Hawai'i Island with contrasting environmental conditions. We found that transcriptome-wide gene expression profiles clearly separated by field site first, and symbiont clade second. With 273 differentially expressed genes (DEGs, 1.3% of all host transcripts), symbiont clade had a measurable effect on host gene expression, but the effect of field site proved almost an order of magnitude higher (1,957 DEGs, 9.6%). According to SNP analysis, we found moderate evidence for host genetic differentiation between field sites (F ST = 0.046) and among corals harboring alternative symbiont clades (F ST = 0.036), suggesting that site-related gene expression differences are likely due to a combination of local adaptation and acclimatization to environmental factors. The correlation between host gene expression and symbiont clade may be due to several factors, including host genotype or microhabitat selecting for alternative clades, host physiology responding to different symbionts, or direct modulation of host gene expression by Symbiodinium. However, the magnitude of these effects at the level of transcription was unexpectedly small considering the contribution of symbiont type to holobiont phenotype.

Metatranscriptomics profile of the gill microbial community during Bathymodiolus azoricus aquarium acclimatization at atmospheric pressure

Background

The deep-sea mussels Bathymodiolus azoricus (Bivalvia: Mytilidae) are the dominant macrofauna subsisting at the hydrothermal vents site Menez Gwen in the Mid-Atlantic Ridge (MAR). Their adaptive success in such challenging environments is largely due to their gill symbiotic association with chemosynthetic bacteria. We examined the response of vent mussels as they adapt to sea-level environmental conditions, through an assessment of the relative abundance of host-symbiont related RNA transcripts to better understand how the gill microbiome may drive host-symbiont interactions in vent mussels during hypothetical venting inactivity.

Results

The metatranscriptome of B. azoricus was sequenced from gill tissues sampled at different time-points during a five-week acclimatization experiment, using Next-Generation-Sequencing. After Illumina sequencing, a total of 181,985,262 paired-end reads of 150 bp were generated with an average of 16,544,115 read per sample. Metatranscriptome analysis confirmed that experimental acclimatization in aquaria accounted for global gill transcript variation. Additionally, the analysis of 16S and 18S rRNA sequences data allowed for a comprehensive characterization of host-symbiont interactions, which included the gradual loss of gill endosymbionts and signaling pathways, associated with stress responses and energy metabolism, under experimental acclimatization. Dominant active transcripts were assigned to the following KEGG categories: “Ribosome”, “Oxidative phosphorylation” and “Chaperones and folding catalysts” suggesting specific metabolic responses to physiological adaptations in aquarium environment.

Conclusions

Gill metagenomics analyses highlighted microbial diversity shifts and a clear pattern of varying mRNA transcript abundancies and expression during acclimatization to aquarium conditions which indicate change in bacterial community activity. This approach holds potential for the discovery of new host-symbiont associations, evidencing new functional transcripts and a clearer picture of methane metabolism during loss of endosymbionts. Towards the end of acclimatization, we observed trends in three major functional subsystems, as evidenced by an increment of transcripts related to genetic information processes; the decrease of chaperone and folding catalysts and oxidative phosphorylation transcripts; but no change in transcripts of gluconeogenesis and co-factors-vitamins.

Physiological impacts of acute Cu exposure on deep-sea vent mussel Bathymodiolus azoricus under a deep-sea mining activity scenario

Over the past years, several studies have been dedicated to understanding the physiological ability of the vent mussel Bathymodiolus azoricus to overcome the high metal concentrations present in their surrounding hydrothermal environment. Potential deep-sea mining activities at Azores Triple junction hydrothermal vent deposits would inevitably lead to the emergence of new fluid sources close to mussel beds, with consequent emission of high metal concentrations and potential resolubilization of Cu from minerals formed during the active phase of the vent field. Copper is an essential metal playing a key role in the activation of metalloenzymes and metalloproteins responsible for important cellular metabolic processes and tissue homeostasis. However, excessive intracellular amounts of reactive Cu ions may cause irreversible damages triggering possible cell apoptosis. In the present study, B. azoricus was exposed to increasing concentrations of Cu for 96h in conditions of temperature and hydrostatic pressure similar to those experienced at the Lucky Strike hydrothermal vent field. Specimens were kept in 1L flasks, exposed to four Cu concentrations: 0μg/L (control), 300, 800 and 1600μg/L and pressurized to 1750bar. We addressed the question of how increased Cu concentration would affect the function of antioxidant defense proteins and expression of antioxidant and immune-related genes in B. azoricus. Both antioxidant enzymatic activities and gene expression were examined in gills, mantle and digestive gland tissues of exposed vent mussels. Our study reveals that stressful short-term Cu exposure has a strong effect on molecular metabolism of the hydrothermal vent mussel, especially in gill tissue. Initially, both the stress caused by unpressurization or by Cu exposure was associated with high antioxidant enzyme activities and tissue-specific transcriptional up-regulation. However, mussels exposed to increased Cu concentrations showed both antioxidant and immune-related gene suppression. Under a mining activity scenario, the release of an excess of dissolved Cu to the vent environment may cause serious changes in cellular defense mechanisms of B. azoricus. This outcome, while adding to our knowledge of Cu toxicity, highlights the potentially deleterious impacts of mining activities on the physiology of deep-sea organisms.

Molecular identification of methane monooxygenase and quantitative analysis of methanotrophic endosymbionts under laboratory maintenance in Bathymodiolus platifrons from the South China Sea

Deep-sea mussels of the genus Bathymodiolus are numerically dominant macrofauna in many cold seep and hydrothermal vent ecosystems worldwide, and they depend on organic carbon produced by symbionts present in the epithelial cells of the gills. Although Bathymodiolus platifrons represents typical methanotrophic endosymbiosis, our understanding of molecular mechanisms of methane oxidization and carbon fixation is still in its infancy. Moreover, the laboratory maintenance of B. platifrons and the symbiont abundance dynamics during maintenance has not been reported. In the present study, we report the first systematic identification and phylogenetic analysis of three subunits of methane monooxygenase (pmoA, pmoB, and pmoC) obtained from the endosymbiotic bacteria found in B. platifrons. The coding sequences (CDS) of the three genes in the B. platifrons endosymbiont were 750, 1,245, and 753 bp, encoding 249, 414, and 250 amino acids, respectively. Sequence alignment and phylogenetic analysis revealed that the symbiont of B. platifrons belongs to the type I methanotrophs. In order to clarify the impact of environmental methane on symbiont abundance, a 34-day laboratory maintenance experiment was conducted in which B. platifrons individuals were acclimatized to methane-present and methane-absent environments. Symbiont abundance was evaluated by calculating the relative DNA content of the methane monooxygenase gene using quantitative real-time PCR. We found that symbiont quantity immediately decreased from its initial level, then continued to gradually decline during maintenance. At 24 and 34 days of maintenance, symbiont abundance in the methane-absent environment had significantly decreased compared to that in the methane-present environment, indicating that the maintenance of symbionts relies on a continuous supply of methane. Our electron microscopy results validated the qPCR analysis. This study enriches our knowledge of the molecular basis and the dynamic changes of the methanotrophic endosymbiosis in B. platifrons, and provides a feasible model biosystem for further investigation of methane oxidization, the carbon fixation process, and environmental adaptations of deep-sea mussels.

Activity of antioxidant enzymes in response to atmospheric pressure induced physiological stress in deep-sea hydrothermal vent mussel Bathymodiolus azoricus

Deep sea hydrothermal Bathymodiolus azoricus mussels from Portuguese EEZ Menez Gwen hydrothermal field possess the remarkable ability to overcome decompression and survive successfully at atmospheric pressure conditions. We investigated the potential use of antioxidant defense enzymes in mussel B. azoricus as biomarkers of oxidative stress induced by long term acclimatization to atmospheric pressure conditions. Mussels collected at Menez Gwen hydrothermal field were acclimatized for two weeks in three distinct conditions suitable of promoting physiological stress, (i) in plain seawater for concomitant endosymbiont bacteria loss, (ii) in plain seawater under metal iron exposure, (iii) constant bubbling methane and pumped sulfide for endosymbiont bacteria survival. The enzymatic activities of superoxide dismutase (SOD), catalase (CAT), and iron storage proteins in addition to electrophoretic profiles were examined in vent mussel gills and digestive gland. Gills showed approximately 3 times more SOD specific activity than digestive glands. On the other hand, digestive glands showed approximately 6 times more CAT specific activity than gills. Iron storage proteins were identified in gill extracts from all experimental conditions mussels. However, in digestive gland extracts only fresh collected mussels and after 2 weeks in FeSO4 showed the presence of iron storage proteins. The differences between SOD, CAT specific activities and the presence of iron storage proteins in the examined tissues reflect dissimilar metabolic and antioxidant activities, as a result of tissue specificities and acclimatization conditions influences on the organism.

Multiple I-Type Lysozymes in the Hydrothermal Vent Mussel Bathymodiolus azoricus and Their Role in Symbiotic Plasticity

The aim of this study was first to identify lysozymes paralogs in the deep sea mussel Bathymodiolus azoricus then to measure their relative expression or activity in different tissue or conditions. B. azoricus is a bivalve that lives close to hydrothermal chimney in the Mid-Atlantic Ridge (MAR). They harbour in specialized gill cells two types of endosymbiont (gram-bacteria): sulphide oxidizing bacteria (SOX) and methanotrophic bacteria (MOX). This association is thought to be ruled by specific mechanism or actors of regulation to deal with the presence of symbiont but these mechanisms are still poorly understood. Here, we focused on the implication of lysozyme, a bactericidal enzyme, in this endosymbiosis. The relative expression of Ba-lysozymes paralogs and the global anti-microbial activity, were measured in natural population (Lucky Strike--1700 m, Mid-Atlantic Ridge), and in in situ experimental conditions. B. azoricus individuals were moved away from the hydrothermal fluid to induce a loss of symbiont. Then after 6 days some mussels were brought back to the mussel bed to induce a re-acquisition of symbiotic bacteria. Results show the presence of 6 paralogs in B. azoricus. In absence of symbionts, 3 paralogs are up-regulated while others are not differentially expressed. Moreover the global activity of lysozyme is increasing with the loss of symbiont. All together these results suggest that lysozyme may play a crucial role in symbiont regulation.

Mercury accumulation in hydrothermal vent mollusks from the southern Tonga Arc, southwestern Pacific Ocean

We provide the mercury (Hg) and monomethylmercury (MMHg) levels of the plume water, sulfide ore, sediment, and mollusks located at the hydrothermal vent fields of the southern Tonga Arc, where active volcanism and intense seismic activity occur frequently. Our objectives were: (1) to address the potential release of Hg from hydrothermal fluids and (2) to examine the distribution of Hg and MMHg levels in hydrothermal mollusks (mussels and snails) harboring chemotrophic bacteria. While high concentrations of Hg in the sediment and Hg, As, and Sb in the sulfide ore indicates that their source is likely hydrothermal fluids, the MMHg concentration in the sediment was orders of magnitude lower than the Hg (<0.001%). It suggests that Hg methylation may have not been favorable in the vent field sediment. In addition, Hg concentrations in the mollusks were much higher (10-100 times) than in other hydrothermal vent environments, indicating that organisms located at the Tonga Arc are exposed to exceedingly high Hg levels. While Hg concentration was higher in the gills and digestive glands than in the mantles and residues of snails and mussels, the MMHg concentrations in the gills and digestive glands were orders of magnitude lower (0.004-0.04%) than Hg concentrations. In summary, our results suggest that the release of Hg from the hydrothermal vent fields of the Tonga Arc and subsequent bioaccumulation are substantial, but not for MMHg.

Post-capture immune gene expression studies in the deep-sea hydrothermal vent mussel Bathymodiolus azoricus acclimatized to atmospheric pressure

Deep-sea hydrothermal vents are extreme habitats that are distributed worldwide in association with volcanic and tectonic events, resulting thus in the establishment of particular environmental conditions, in which high pressure, steep temperature gradients, and potentially toxic concentrations of sulfur, methane and heavy metals constitute driving factors for the foundation of chemosynthetic-based ecosystems. Of all the different macroorganisms found at deep-sea hydrothermal vents, the mussel Bathymodiolus azoricus is the most abundant species inhabiting the vent ecosystems from the Mid-Atlantic Ridge (MAR). In the present study, the effect of long term acclimatization at atmospheric pressure on host-symbiotic associations were studied in light of the ensuing physiological adaptations from which the immune and endosymbiont gene expressions were concomitantly quantified by means of real-time PCR. The expression of immune genes at 0 h, 12 h, 24 h, 36 h, 48 h, 72 h, 1 week and 3 weeks post-capture acclimatization was investigated and their profiles compared across the samples tested. The gene signal distribution for host immune and bacterial genes followed phasic changes in gene expression at 24 h, 1 week and 3 weeks acclimatization when compared to other time points tested during this temporal expression study. Analyses of the bacterial gene expression also suggested that both bacterial density and activity could contribute to shaping the intricate association between endosymbionts and host immune genes whose expression patterns seem to be concomitant at 1 week acclimatization. Fluorescence in situ hybridization was used to assess the distribution and prevalence of endosymbiont bacteria within gill tissues confirming the gradual loss of sulfur-oxidizing (SOX) and methane-oxidizing (MOX) bacteria during acclimatization. The present study addresses the deep-sea vent mussel B. azoricus as a model organism to study how acclimatization in aquaria and the prevalence of symbiotic bacteria are driving the expression of host immune genes. Tight associations, unseen thus far, suggest that host immune and bacterial gene expression patterns reflect distinct physiological responses over the course of acclimatization under aquarium conditions.

Transcriptomic response of the hydrothermal mussel Bathymodiolus azoricus in experimental exposure to heavy metals is modulated by the Pgm genotype and symbiont content

Hydrothermal vent mussels belonging to the genus Bathymodiolus dominate communities at hydrothermal sites of the Mid-Atlantic Ridge. The mussel Bathymodiolus azoricus harbors thiotrophic and methanotrophic symbiotic bacteria in its gills and evolves in naturally highly metal contaminated environments. In the context of investigations on metal tolerance/effect in B. azoricus, we focused our work on the short-term adaptive response (15days) of mussels to different metals exposure at a molecular level using metal concentrations chosen to mimic natural situations at three vents sites. The expression of a set of 38 genes involved in different steps of the metal uptake, detoxication and various metabolisms was analysed by qPCR. Mussels were also genotyped at 10 enzyme loci to explore the relationships among natural genetic variation and gene expression. Relation between symbiont content (both sulfur-oxidizing and methanogen bacteria) and gene expression was also analysed. Our study demonstrated the influence of metal cocktail composition and time exposure on the transcriptome regulation with a specific pattern of regulation observed for the three metal cocktail tested. We also evidenced the significant influence of some specific Pgm genotype on the global gene expression in our experimental populations and a general trend of a higher gene expression in individuals carrying a high symbiont content.

Comparative study of immune responses in the deep-sea hydrothermal vent mussel Bathymodiolus azoricus and the shallow-water mussel Mytilus galloprovincialis challenged with Vibrio bacteria

The deep-sea hydrothermal vent mussel Bathymodiolus azoricus and the continental European coast Mytilus galloprovincialis are two bivalves species living in highly distinct marine habitats. Mussels are filter-feeding animals that may accumulate rapidly bacteria from the environment. Contact with microorganism is thus inevitable during feeding processes where gill tissues assume a strategic importance at the interface between the external milieu and the internal body cavities promoting interactions with potential pathogens during normal filtration and a constant challenge to their immune system. In the present study B. azoricus and M. galloprovincialis were exposed to Vibrio alginolyticus, Vibrio anguillarum and Vibrio splendidus suspensions and to a mixture of these Vibrio suspensions, in order to ascertain the expression level of immune genes in gill samples, from both mussel species. The immune gene expressions were analyzed by means of quantitative-Polymerase Chain Reaction (qPCR). The gene expression results revealed that these bivalve species exhibit significant expression differences between 12 h and 24 h post-challenge times, and between the Vibrio strains used. V. splendidus induced the strongest gene expression level in the two bivalve species whereas the NF-κB and Aggrecan were the most significantly differentially expressed between the two mussel species. When comparing exposure times, both B. azoricus and M. galloprovincialis showed similar percentage of up-regulated genes at 12 h while a marked increased of gene expression was observed at 24 h for the majority of the immune genes in M. galloprovincialis. This contrasts with B. azoricus where the majority of the immune genes were down-regulated at 24 h. The 24 h post-challenge gene expression results clearly bring new evidence supporting time-dependent transcriptional activities resembling acute phase-like responses and different immune responses build-up in these two mussel species when challenged with Vibrio bacteria. High Pressure Liquid Chromatography (HPLC)-Electrospray ionization mass spectrometry (ESI-MS/MS) analyses resulted in different peptide sequences from B. azoricus and M. galloprovincialis gill tissues suggesting that naïve animals present differences, at the protein synthesis level, in their natural environment. B. azoricus proteins sequences, mostly of endosymbiont origin, were related to metabolic, energy production, protein synthesis processes and nutritional demands whereas in M. galloprovincialis putative protein functions were assumed to be related to structural and cellular integrity and signaling functions.

Site-related differences in gene expression and bacterial densities in the mussel Bathymodiolus azoricus from the Menez Gwen and Lucky Strike deep-sea hydrothermal vent sites

The deep-sea hydrothermal vent mussel Bathymodiolus azoricus is a symbiont bearing bivalve that is found in great abundance at the Menez Gwen and Lucky Strike hydrothermal vent sites and in close vicinity of the Azores region near the Mid-Atlantic Ridge (MAR). The physiological relationships that vent mussels have developed with their physical and chemical environments are likely to influence global gene expression profiles providing thus the means to investigate distinct biological markers predicting the origin of Bathymodiolus sp. irrespectively of their geographical localization. Differences found at gene expression levels, and between fluorescence in situ hybridization (FISH) and 16S rRNA amplicon sequencing results provided experimental evidence for the distinction of both Menez Gwen and Lucky Strike vent mussel individuals based on bacterial and vent mussel gene expression signatures and on the constitutive distribution and relative abundance of endosymbiotic bacteria within gill tissues. Our results confirmed the presence of methanotroph endosymbionts in Menez Gwen vent mussels whereas Lucky Strike specimens seem to harbor a different bacterial morphotype when a methane monooxygenase gene specific probe was used. No qualitative differences could be visualized between Menez Gwen and Lucky Strike individuals when tested with a sulfur-oxidizing-related probe. Quantitative PCR (qPCR) studies revealed different gene expression profiles in both Menez Gwen and Lucky Strike mussel gill tissues for the immune genes selected. Genes encoding transcription factors presented noticeably low levels of fold expression whether in Menez Gwen or Lucky Strike animals whereas the genes encoding effector molecules appeared to have higher levels expression in gill tissues from Menez Gwen animals. The peptidoglycan recognition molecule encoding gene, PGRP, presented the highest level of transcriptional activity among the genes analyzed in Menez Gwen mussel gill tissues, seconded by carcinolectin and thus denoting the relevance of immune recognition molecules in early stage of the immune responses onset. Genes regarded as encoding molecules involved in signaling pathways were consistently expressed in both Menez Gwen and Lucky Strike mussel gill tissues. Remarkably, the immunity-related GTPase encoding gene demonstrated, in Lucky Strike samples, the highest level of expression among the signaling molecule encoding genes tested when expressions levels were compared between Menez Gwen and Lucky Strike animals. A differential expression analysis of bacterial genes between Menez Gwen and Lucky Strike mussels indicated a clear expression signature in the latter animal gill tissues. The bacterial community structure ensued from the 16S rRNA sequencing analyses pointed at an unpredicted conservation of endosymbiont bacterial loads between Menez Gwen and Lucky Strike samples. Taken together, our results support the hypothesis that B. azoricus exhibits different transcriptional statuses while living in distinct hydrothermal vent sites may result in distinct gene expressions because of physico-chemical and/or symbiont densities differences.

The transcriptome of Bathymodiolus azoricus gill reveals expression of genes from endosymbionts and free-living deep-sea bacteria

Deep-sea environments are largely unexplored habitats where a surprising number of species may be found in large communities, thriving regardless of the darkness, extreme cold, and high pressure. Their unique geochemical features result in reducing environments rich in methane and sulfides, sustaining complex chemosynthetic ecosystems that represent one of the most surprising findings in oceans in the last 40 years. The deep-sea Lucky Strike hydrothermal vent field, located in the Mid Atlantic Ridge, is home to large vent mussel communities where Bathymodiolus azoricus represents the dominant faunal biomass, owing its survival to symbiotic associations with methylotrophic or methanotrophic and thiotrophic bacteria. The recent transcriptome sequencing and analysis of gill tissues from B. azoricus revealed a number of genes of bacterial origin, hereby analyzed to provide a functional insight into the gill microbial community. The transcripts supported a metabolically active microbiome and a variety of mechanisms and pathways, evidencing also the sulfur and methane metabolisms. Taxonomic affiliation of transcripts and 16S rRNA community profiling revealed a microbial community dominated by thiotrophic and methanotrophic endosymbionts of B. azoricus and the presence of a Sulfurovum-like epsilonbacterium.