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

Genome assembly of the chemosynthetic endosymbiont of the hydrothermal vent snail Alviniconcha adamantis from the Mariana Arc

Chemosynthetic animal-microbe symbioses sustain hydrothermal vent communities in the global deep sea. In the Indo-Pacific Ocean, hydrothermal ecosystems are often dominated by gastropod species of the genus Alviniconcha, which live in association with chemosynthetic Gammaproteobacteria or Campylobacteria. While the symbiont genomes of most extant Alviniconcha species have been sequenced, no genome information is currently available for the gammaproteobacterial endosymbiont of Alviniconcha adamantis-a comparatively shallow living species that is thought to be the ancestor to all other present Alviniconcha lineages. Here, we report the first genome sequence for the symbiont of A. adamantis from the Chamorro Seamount at the Mariana Arc. Our phylogenomic analyses show that the A. adamantis symbiont is most closely related to Chromatiaceae endosymbionts of the hydrothermal vent snails Alviniconcha strummeri and Chrysomallon squamiferum, but represents a distinct bacterial species or possibly genus. Overall, the functional capacity of the A. adamantis symbiont appeared to be similar to other chemosynthetic Gammaproteobacteria, though several flagella and chemotaxis genes were detected, which are absent in other gammaproteobacterial Alviniconcha symbionts. These differences might suggest potential contrasts in symbiont transmission dynamics, host recognition, or nutrient transfer. Furthermore, an abundance of genes for ammonia transport and urea usage could indicate adaptations to the oligotrophic waters of the Mariana region, possibly via recycling of host- and environment-derived nitrogenous waste products. This genome assembly adds to the growing genomic resources for chemosynthetic bacteria from hydrothermal vents and will be valuable for future comparative genomic analyses assessing gene content evolution in relation to environment and symbiotic lifestyles.

Host-Microbe Interactions in the Chemosynthetic Riftia pachyptila Symbiosis

All animals are associated with microorganisms; hence, host-microbe interactions are of fundamental importance for life on earth. However, we know little about the molecular basis of these interactions. Therefore, we studied the deep-sea Riftia pachyptila symbiosis, a model association in which the tubeworm host is associated with only one phylotype of endosymbiotic bacteria and completely depends on this sulfur-oxidizing symbiont for nutrition. Using a metaproteomics approach, we identified both metabolic interaction processes, such as substrate transfer between the two partners, and interactions that serve to maintain the symbiotic balance, e.g., host efforts to control the symbiont population or symbiont strategies to modulate these host efforts. We suggest that these interactions are essential principles of mutualistic animal-microbe associations.

The deep-sea tubeworm Riftia pachyptila lacks a digestive system but completely relies on bacterial endosymbionts for nutrition. Although the symbiont has been studied in detail on the molecular level, such analyses were unavailable for the animal host, because sequence information was lacking. To identify host-symbiont interaction mechanisms, we therefore sequenced the Riftia transcriptome, which served as a basis for comparative metaproteomic analyses of symbiont-containing versus symbiont-free tissues, both under energy-rich and energy-limited conditions. Our results suggest that metabolic interactions include nutrient allocation from symbiont to host by symbiont digestion and substrate transfer to the symbiont by abundant host proteins. We furthermore propose that Riftia maintains its symbiont by protecting the bacteria from oxidative damage while also exerting symbiont population control. Eukaryote-like symbiont proteins might facilitate intracellular symbiont persistence. Energy limitation apparently leads to reduced symbiont biomass and increased symbiont digestion. Our study provides unprecedented insights into host-microbe interactions that shape this highly efficient symbiosis.

Metagenomic Signatures of Microbial Communities in Deep-Sea Hydrothermal Sediments of Azores Vent Fields

The organisms inhabiting the deep-seafloor are known to play a crucial role in global biogeochemical cycles. Chemolithoautotrophic prokaryotes, which produce biomass from single carbon molecules, constitute the primary source of nutrition for the higher organisms, being critical for the sustainability of food webs and overall life in the deep-sea hydrothermal ecosystems. The present study investigates the metabolic profiles of chemolithoautotrophs inhabiting the sediments of Menez Gwen and Rainbow deep-sea vent fields, in the Mid-Atlantic Ridge. Differences in the microbial community structure might be reflecting the distinct depth, geology, and distance from vent of the studied sediments. A metagenomic sequencing approach was conducted to characterize the microbiome of the deep-sea hydrothermal sediments and the relevant metabolic pathways used by microbes. Both Menez Gwen and Rainbow metagenomes contained a significant number of genes involved in carbon fixation, revealing the largely autotrophic communities thriving in both sites. Carbon fixation at Menez Gwen site was predicted to occur mainly via the reductive tricarboxylic acid cycle, likely reflecting the dominance of sulfur-oxidizing Epsilonproteobacteria at this site, while different autotrophic pathways were identified at Rainbow site, in particular the Calvin-Benson-Bassham cycle. Chemolithotrophy appeared to be primarily driven by the oxidation of reduced sulfur compounds, whether through the SOX-dependent pathway at Menez Gwen site or through reverse sulfate reduction at Rainbow site. Other energy-yielding processes, such as methane, nitrite, or ammonia oxidation, were also detected but presumably contributing less to chemolithoautotrophy. This work furthers our knowledge of the microbial ecology of deep-sea hydrothermal sediments and represents an important repository of novel genes with potential biotechnological interest.

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.

Sediment Microbial Diversity of Three Deep-Sea Hydrothermal Vents Southwest of the Azores

Menez Gwen, Lucky Strike and Rainbow are the three most visited and well-known deep-sea hydrothermal vent fields in the Azores region, located in the Mid-Atlantic Ridge. Their distinct geological and ecological features allow them to support a diversity of vent communities, which are largely dependent on Bacteria and Archaea capable of anaerobic or microaerophilic metabolism. These communities play important ecological roles through chemoautotrophy, feeding and in establishing symbiotic associations. However, the occurrence and distribution of these microbes remain poorly understood, especially in deep-sea sediments. In this study, we provide for the first time a comparative survey of the sediment-associated microbial communities from these three neighbouring vent fields. Sediment samples collected in the Menez Gwen, Lucky Strike and Rainbow vent fields showed significant differences in trace-metal concentrations and associated microbiomes. The taxonomic profiles of bacterial, archaeal and eukaryotic representatives were assessed by rRNA gene-tag pyrosequencing, identified anaerobic methanogens and microaerobic Epsilonproteobacteria, particularly at the Menez Gwen site, suggesting sediment communities potentially enriched in sub-seafloor microbes rather than from pelagic microbial taxa. Cosmopolitan OTUs were also detected mostly at Lucky Strike and Rainbow sites and affiliated with the bacterial clades JTB255, Sh765B-TzT-29, Rhodospirillaceae and OCS155 marine group and with the archaeal Marine Group I. Some variations in the community composition along the sediment depth were revealed. Elemental contents and hydrothermal influence are suggested as being reflected in the composition of the microbial assemblages in the sediments of the three vent fields. Altogether, these findings represent valuable information for the understanding of the microbial distribution and potential ecological roles in deep-sea hydrothermal fields.

Comparative transcriptome and potential antiviral signaling pathways analysis of the gills in the red swamp crayfish, Procambarus clarkii infected with White Spot Syndrome Virus (WSSV)

Red swamp crayfish is an important model organism for research of the invertebrate innate immunity mechanism. Its excellent disease resistance against bacteria, fungi, and viruses is well-known. However, the antiviral mechanisms of crayfish remain unclear. In this study, we obtained high-quality sequence reads from normal and white spot syndrome virus (WSSV)-challenged crayfish gills. For group normal (GN), 39,390,280 high-quality clean reads were randomly assembled to produce 172,591 contigs; whereas, 34,011,488 high-quality clean reads were randomly assembled to produce 182,176 contigs for group WSSV-challenged (GW). After GO annotations analysis, a total of 35,539 (90.01%), 14,931 (37.82%), 28,221 (71.48%), 25,290 (64.05%), 15,595 (39.50%), and 13,848 (35.07%) unigenes had significant matches with sequences in the Nr, Nt, Swiss-Prot, KEGG, COG and GO databases, respectively. Through the comparative analysis between GN and GW, 12,868 genes were identified as differentially up-regulated DEGs, and 9,194 genes were identified as differentially down-regulated DEGs. Ultimately, these DEGs were mapped into different signaling pathways, including three important signaling pathways related to innate immunity responses. These results could provide new insights into crayfish antiviral immunity mechanism.

Enhanced recovery of Arcobacter spp. using NaCl in culture media and re-assessment of the traits of Arcobacter marinus and Arcobacter halophilus isolated from marine water and shellfish

The genus Arcobacter is a relatively poorly known group of bacteria, and the number of new species and sequences from non-culturable strains has increased considerably in recent years. This study investigates whether using media that contain NaCl might help to improve the recovery of Arcobacter spp. from marine environments. To this aim, 62 water and shellfish samples were analysed in parallel, with both a commonly used culture method (enrichment in Arcobacter-CAT broth followed by culture on Blood Agar) and a new one that supplements the Arcobacter-CAT enrichment broth with 2.5% NaCl (w/v) followed by culturing on Marine Agar. The new method yielded ca. 40% more positive samples and provided a higher diversity of known (11 vs. 7) and unknown (7 vs. 2) Arcobacter species. Among the 11 known species recovered, Arcobacter marinus and Arcobacter halophilus were isolated only by this new method. No more strains of these species have been isolated since their original descriptions, both of which were based only on a single strain. In view of that, the phenotypic characteristics of these species are re-evaluated in the present study, using the new strains. Strains of A. halophilus had the same phenotypic profile as the type strain. However, some strains of A. marinus differed from the type strain in that they did not hydrolyse indoxyl-acetate, becoming, therefore, the first Arcobacter species to show a varying ability to hydrolyse indoxyl-acetate. This study shows to what extent a simple variation to the culture media can have a big influence on positive samples and on the community of species recovered.

In-depth comparative transcriptome analysis of intestines of red swamp crayfish, Procambarus clarkii, infected with WSSV

Crayfish has become one of the most important farmed aquatic species in China due to its excellent disease resistance against bacteria and viruses. However, the antiviral mechanism of crayfish is still not very clear. In the present study, many high-quality sequence reads from crayfish intestine were obtained using Illumina-based transcriptome sequencing. For the normal group (GN), 44,600,142 high-quality clean reads were randomly assembled to produce 125,394 contigs. For the WSSV-challenged group (GW), 47,790,746 high-quality clean reads were randomly assembled to produce 148,983 contigs. After GO annotation, 39,482 unigenes were annotated into three ontologies: biological processes, cellular components, and molecular functions. In addition, 15,959 unigenes were mapped to 25 different COG categories. Moreover, 7,000 DEGs were screened out after a comparative analysis between the GN and GW samples, which were mapped into 250 KEGG pathways. Among these pathways, 36 were obviously changed (P-values < 0.05) and 28 pathways were extremely significantly changed (P-values < 0.01). Finally, five key DEGs involved in the JAK-STAT signaling pathway were chosen for qRT-PCR. The results showed that these five DEGs were obviously up-regulated at 36 h post WSSV infection in crayfish intestine. These results provide new insight into crayfish antiviral immunity mechanisms.

Analysis of synonymous codon usage patterns in sixty-four different bivalve species

Synonymous codon usage bias (CUB) is a defined as the non-random usage of codons encoding the same amino acid across different genomes. This phenomenon is common to all organisms and the real weight of the many factors involved in its shaping still remains to be fully determined. So far, relatively little attention has been put in the analysis of CUB in bivalve mollusks due to the limited genomic data available. Taking advantage of the massive sequence data generated from next generation sequencing projects, we explored codon preferences in 64 different species pertaining to the six major evolutionary lineages in Bivalvia. We detected remarkable differences across species, which are only partially dependent on phylogeny. While the intensity of CUB is mild in most organisms, a heterogeneous group of species (including Arcida and Mytilida, among the others) display higher bias and a strong preference for AT-ending codons. We show that the relative strength and direction of mutational bias, selection for translational efficiency and for translational accuracy contribute to the establishment of synonymous codon usage in bivalves. Although many aspects underlying bivalve CUB still remain obscure, we provide for the first time an overview of this phenomenon in this large, commercially and environmentally important, class of marine invertebrates.

Wine fermentation microbiome: a landscape from different Portuguese wine appellations

Grapes and wine musts harbor a complex microbiome, which plays a crucial role in wine fermentation as it impacts on wine flavour and, consequently, on its final quality and value. Unveiling the microbiome and its dynamics, and understanding the ecological factors that explain such biodiversity, has been a challenge to oenology. In this work, we tackle this using a metagenomics approach to describe the natural microbial communities, both fungal and bacterial microorganisms, associated with spontaneous wine fermentations. For this, the wine microbiome, from six Portuguese wine appellations, was fully characterized as regards to three stages of fermentation - Initial Musts (IM), and Start and End of alcoholic fermentations (SF and EF, respectively). The wine fermentation process revealed a higher impact on fungal populations when compared with bacterial communities, and the fermentation evolution clearly caused a loss of the environmental microorganisms. Furthermore, significant differences (p < 0.05) were found in the fungal populations between IM, SF, and EF, and in the bacterial population between IM and SF. Fungal communities were characterized by either the presence of environmental microorganisms and phytopathogens in the IM, or yeasts associated with alcoholic fermentations in wine must samples as Saccharomyces and non-Saccharomyces yeasts (as Lachancea, Metschnikowia, Hanseniaspora, Hyphopichia, Sporothrix, Candida, and Schizosaccharomyces). Among bacterial communities, the most abundant family was Enterobacteriaceae; though families of species associated with the production of lactic acid (Lactobacillaceae, Leuconostocaceae) and acetic acid (Acetobacteriaceae) were also detected. Interestingly, a biogeographical correlation for both fungal and bacterial communities was identified between wine appellations at IM suggesting that each wine region contains specific and embedded microbial communities which may contribute to the uniqueness of regional wines.

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.

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.

Gene discovery through transcriptome sequencing for the invasive mussel Limnoperna fortunei

The success of the Asian bivalve Limnoperna fortunei as an invader in South America is related to its high acclimation capability. It can inhabit waters with a wide range of temperatures and salinity and handle long-term periods of air exposure. We describe the transcriptome of L. fortunei aiming to give a first insight into the phenotypic plasticity that allows non-native taxa to become established and widespread. We sequenced 95,219 reads from five main tissues of the mussel L. fortunei using Roche’s 454 and assembled them to form a set of 84,063 unigenes (contigs and singletons) representing partial or complete gene sequences. We annotated 24,816 unigenes using a BLAST sequence similarity search against a NCBI nr database. Unigenes were divided into 20 eggNOG functional categories and 292 KEGG metabolic pathways. From the total unigenes, 1,351 represented putative full-length genes of which 73.2% were functionally annotated. We described the first partial and complete gene sequences in order to start understanding bivalve invasiveness. An expansion of the hsp70 gene family, seen also in other bivalves, is present in L. fortunei and could be involved in its adaptation to extreme environments, e.g. during intertidal periods. The presence of toll-like receptors gives a first insight into an immune system that could be more complex than previously assumed and may be involved in the prevention of disease and extinction when population densities are high. Finally, the apparent lack of special adaptations to extremely low O₂ levels is a target worth pursuing for the development of a molecular control approach.

Bivalve omics: state of the art and potential applications for the biomonitoring of harmful marine compounds

The extraordinary progress experienced by sequencing technologies and bioinformatics has made the development of omic studies virtually ubiquitous in all fields of life sciences nowadays. However, scientific attention has been quite unevenly distributed throughout the different branches of the tree of life, leaving molluscs, one of the most diverse animal groups, relatively unexplored and without representation within the narrow collection of well established model organisms. Within this Phylum, bivalve molluscs play a fundamental role in the functioning of the marine ecosystem, constitute very valuable commercial resources in aquaculture, and have been widely used as sentinel organisms in the biomonitoring of marine pollution. Yet, it has only been very recently that this complex group of organisms became a preferential subject for omic studies, posing new challenges for their integrative characterization. The present contribution aims to give a detailed insight into the state of the art of the omic studies and functional information analysis of bivalve molluscs, providing a timely perspective on the available data resources and on the current and prospective applications for the biomonitoring of harmful marine compounds.

Metaproteomics reveals abundant transposase expression in mutualistic endosymbionts

Transposases, enzymes that catalyze the movement of mobile genetic elements, are the most abundant genes in nature. While many bacteria encode an abundance of transposases in their genomes, the current paradigm is that the expression of transposase genes is tightly regulated and generally low due to its severe mutagenic effects. In the current study, we detected the highest number of transposase proteins ever reported in bacteria, in symbionts of the gutless marine worm Olavius algarvensis with metaproteomics. At least 26 different transposases from 12 different families were detected, and genomic and proteomic analyses suggest that many of these are active. This high expression of transposases indicates that the mechanisms for their tight regulation have been disabled or no longer exist.