High-throughput sequencing and analysis of the gill tissue transcriptome from the deep-sea hydrothermal vent mussel Bathymodiolus azoricus

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.

Transcriptome analysis of gene expression profiling from the deep sea in situ to the laboratory for the cold seep mussel Gigantidas haimaensis

Background

The deep-sea mussel Gigantidas haimaensis is a representative species from the Haima cold seep ecosystem in the South China Sea that establishes endosymbiosis with chemotrophic bacteria. During long-term evolution, G. haimaensis has adapted well to the local environment of cold seeps. Until now, adaptive mechanisms responding to environmental stresses have remained poorly understood.

Results

In this study, transcriptomic analysis was performed for muscle tissue of G. haimaensis in the in situ environment (MH) and laboratory environment for 0 h (M0), 3 h (M3) and 9 h (M9), and 187,368 transcript sequences and 22,924 annotated differentially expressed genes (DEGs) were generated. Based on Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, these DEGs were enriched with a broad spectrum of biological processes and pathways, including those associated with antioxidants, apoptosis, chaperones, immunity and metabolism. Among these significantly enriched pathways, protein processing in the endoplasmic reticulum and metabolism were the most affected metabolic pathways. These results may imply that G. haimaensis struggles to support the life response to environmental change by changing gene expression profiles.

Conclusion

The present study provides a better understanding of the biological responses and survival strategies of the mussel G. haimaensis from deep sea in situ to the laboratory environment.

Molecular characterization of peptidoglycan recognition proteins from Mytilus coruscus

Mytilus shows great immune resistance to various bacteria from the living waters, indicating a complex immune recognition mechanism against various microbes. Peptidoglycan recognition proteins (PGRPs) play an important role in the defense against invading microbes via the recognition of the immunogenic substance peptidoglycan (PGN). Therefore, eight PGRPs were identified from the gill transcriptome of Mytilus coruscus. The sequence features, expression pattern in various organs and larval development stages, and microbes induced expression profiles of these Mytilus PGRPs were determined. Our data revealed the constitutive expression of PGRPs in various organs with relative higher expression level in immune-related organs. The expression of PGRPs is developmentally regulated, and most PGRPs are undetectable in larvae stages. The expression level of most PGRPs was significantly increased with in vivo microbial challenges, showing strong response to Gram-positive strain in gill and digestive gland, strong response to Gram-negative strain in hemocytes, and relative weaker response to fungus in the three tested organs. In addition, the function analysis of the representative recombinant expressed PGRP (rMcPGRP-2) confirmed the antimicrobial and agglutination activities, showing the immune-related importance of PGRP in Mytilus. Our work suggests that Mytilus PGRPs can act as pattern recognition receptors to recognize the invading microorganisms and the antimicrobial effectors during the innate immune response of Mytilus.

Proteome and Transcriptome Analysis of Gonads Reveals Intersex in Gigantidas haimaensis

Sex has proven to be one of the most intriguing areas of research across evolution, development, and ecology. Intersex or sex change occurs frequently in molluscs. The deep-sea mussel Gigantidas haimaensis often dominates within Haima cold seep ecosystems, but details of their reproduction remain unknown. Herein, we conducted a combined proteomic and transcriptomic analysis of G. haimaensis gonads to provide a systematic understanding of sexual development in deep-sea bivalves. A total of 2,452 out of 42,238 genes (5.81%) and 288 out of 7,089 proteins (4.06%) were significantly differentially expressed between ovaries and testes with a false discovery rate (FDR) <0.05. Candidate genes involved in sexual development were identified; among 12 differentially expressed genes between sexes, four ovary-biased genes (β-catenin, fem-1, forkhead box L2 and membrane progestin receptor α) were expressed significantly higher in males than females. Combining histological characteristics, we speculate that the males maybe intersex undergoing sex change, and implied that these genes may be involved in the process of male testis converting into female gonads in G. haimaensis. The results suggest that this adaptation may be based on local environmental factors, sedentary lifestyles, and patchy distribution, and sex change may facilitate adaptation to a changing environment and expansion of the population. The findings provide a valuable genetic resource to better understand the mechanisms of sex change and survival strategies in deep-sea bivalves.

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.

Comparison of lipids and fatty acids among tissues of two semiterrestrial crabs reveals ecophysiological adaptations in changing coastal environments

Decapods have successfully colonized changing coastal habitats throughout the world by adapting their behavior, physiology, and biochemistry. Biochemical reserves, such as lipids and fatty acids (FAs), play fundamental roles in this adaptation process. These energy reserves are key for the development of decapods and their composition mainly depends on the type and quality of food available in their habitats. This study evaluated the lipid content and FA composition of three tissues (hepatopancreas, gills, and muscle) in two widely distributed, semi-terrestrial coastal crab species in Chile, Cyclograpsus cinereus from the upper intertidal and Hemigrapsus crenulatus from estuaries. This evaluation aimed to assess the physiological role of the bioenergetic reserves of these crabs, which tolerate fluctuating environmental conditions. Our results showed that both species had a higher lipid content in the hepatopancreas and a lower lipid content in its gills and muscle. All three of the evaluated tissues in C. cinereus showed high contents of saturated fatty acids (SFAs), and its hepatopancreas displayed the highest contents of monounsaturated (MUFAs) and polyunsaturated fatty acids (PUFAs). In turn, H. crenulatus had the highest contents of MUFAs and PUFAs in its gills and muscle tissues, including an important amount of eicosapentaenoic acid (EPA). The FA content of C. cinereus may indicate an adaptive physiological response aimed at maintaining its cellular fluid balance during periods of desiccation in the upper intertidal zone. In contrast, the FAs found in H. crenulatus may be linked to the high activity of the sodium‑potassium pump in its gills, in order to maintain osmoregulation in estuaries.

RNA-seq analysis of the mantle transcriptome from Mytilus edulis during a seasonal spawning event in deep and shallow water culture sites on the northeast coast of Newfoundland, Canada

The blue mussel (Mytilus edulis) has global commercial and ecological importance both in wild and cultured conditions. However there is a qualitative and quantitative lack of knowledge of the molecular mechanisms associated with its reproductive physiology, especially with reference to environmental interactions. Here we initiated a transcriptomic analysis (RNA-sequencing (RNA-seq)) of the mantle from both sexes sampled during a seasonal spawning event and from two culture depths (shallow-5 m; deep- 15 m). Mantle libraries were produced from 3 males and 3 females sampled from each of two shallow sites and two deep sites for a total of 12 replicate male and 12 replicate female libraries (24 total libraries). Overall a total of 2.3 billion raw 100 base reads with an average of 96.5 million reads/library were obtained and assembled into 296,118 transcripts with an average length of 568 bp. Overall, 315 transcripts from male libraries and 25 from female libraries were found to be upregulated in deep water as compared to shallow (edgeR adjusted p value ≤ 0.05). Conversely, 126 transcripts from male libraries and 135 from female libraries were found to be significantly downregulated at the same depth. Thirteen transcripts were selected for qPCR validation based on importance in reproduction, antimicrobial defense and metabolism. Of these, 9 RNA-seq identified transcripts were shown by qPCR to be differentially expressed between groups: 2 were upregulated in deep compared with shallow water (dhx38, mt-co1), 2 were upregulated for female compared with male mantle (pias2, mapkap1) and 6 genes (fndc3a, acbd3, klhl10, ccnb3, armc4, mt-co1) showed to be upregulated in males compared to females. The majority of qPCR studied transcripts were identified as involved in gamete development based on the UniProt database. This study further characterizes the importance of the mantle transcriptome during reproductive activities of M. edulis.

Phylogenetic Relationships and Adaptation in Deep-Sea Mussels: Insights from Mitochondrial Genomes

Mitochondrial genomes (mitogenomes) are an excellent source of information for phylogenetic and evolutionary studies, but their application in marine invertebrates is limited. In the present study, we utilized mitogenomes to elucidate the phylogeny and environmental adaptation in deep-sea mussels (Mytilidae: Bathymodiolinae). We sequenced and assembled seven bathymodioline mitogenomes. A phylogenetic analysis integrating the seven newly assembled and six previously reported bathymodioline mitogenomes revealed that these bathymodiolines are divided into three well-supported clades represented by five Gigantidas species, six Bathymodiolus species, and two "Bathymodiolus" species, respectively. A Common interval Rearrangement Explorer (CREx) analysis revealed a gene order rearrangement in bathymodiolines that is distinct from that in other shallow-water mytilids. The CREx analysis also suggested that reversal, transposition, and tandem duplications with subsequent random gene loss (TDRL) may have been responsible for the evolution of mitochondrial gene orders in bathymodiolines. Moreover, a comparison of the mitogenomes of shallow-water and deep-sea mussels revealed that the latter lineage has experienced relaxed purifying selection, but 16 residues of the atp6, nad4, nad2, cob, nad5, and cox2 genes have underwent positive selection. Overall, this study provides new insights into the phylogenetic relationships and mitogenomic adaptations of deep-sea mussels.

Transcriptomic Response of Mussel Gills After a Vibrio splendidus Infection Demonstrates Their Role in the Immune Response

Mussels (Mytilus galloprovincialis) are filter feeder bivalves that are constantly in contact with a wide range of microorganisms, some of which are potentially pathogenic. How mussels recognize and respond to pathogens has not been fully elucidated to date; therefore, we investigated the immune mechanisms that these animals employ in response to a bacterial bath infection from the surrounding water, mimicking the response that mussels mount under natural conditions. After the bath infection, mussels were able to remove the bacteria from their bodies and from the water tank. Accordingly, antibacterial activity was detected in gill extracts, demonstrating that this tissue plays a central role in removing and clearing potential pathogens. A transcriptomic study performed after a bath infection with Vibrio splendidus identified a total of 1,156 differentially expressed genes. The expression levels of genes contributing to a number of biological processes, such as immune response activation pathways and their regulation with cytokines, cell recognition, adhesion and apoptosis, were significantly modulated after infection, suggesting that the gills play important roles in pathogen recognition, as well as being activators and regulators of the mussel innate immune response. In addition to RNA-seq analysis, long non-coding RNAs and their neighboring genes were also analyzed and exhibited modulation after the bacterial challenge. The response of gills against bath infection was compared with the findings of a previous transcriptomic study on hemocytes responding to systemic infection, demonstrating the different and specific functions of gills. The results of this study indicate that recognition processes occur in the gill, thereby activating the effector agents of the immune response to overcome bacterial infection.

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.

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.

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.

De novo transcriptome analysis of stressed blood clam (Anadara broughtonii) and identification of genes associated with hemoglobin

BACKGROUND

Blood clam (Anadara broughtonii) is a commercially important marine bivalve characterised by the red blood. Recently, the clams have been subjected to severe resource recession. Multiple environmental stressors are indispensable for the recession.

OBJECTIVE

We aimed to investigate the transcriptome information of blood clam under environmental stressors.

METHODS

Paired-end Illumina HiSeq™ 2500 sequencing technology was employed for cDNA library construction and Illumina sequencing. Several public databases were introduced for gene annotation, and Kyoto Encyclopedia of Genes and Genomes (KEGG) database was used for pathways analyses. The Open Reading Frame of annotated hemoglobin (Hb) was predicted and validated by DNAMAN 6.0 and NCBI BLASTx, respectively.

RESULTS

A total of 242,919 transcripts were generated, 116,264 unigenes were subsequently assembled with an average length of 747 base pairs, and 33,776 unigenes were successfully annotated. Gene Ontology (GO) categories indicated that the terms of cellular processes, metabolic processes, cell, cell part, binding, and catalytic activity were dominant. KEGG pathway analyses suggested ribosome, oxytocin, focal adhesion, Ras, and PI3K-Akt were the largest signaling pathway groups, and many presented pathways (Ras, Rap1, and MAPK, etc.) were related to apoptosis, immune and stress response. In addition, a total of 19,306 potential simple sequence repeats (SSRs) were detected in 15,852 sequences. Six hemoglobin-related genes with complete conserved domain sequences were identified and 3 of them were predicted as HbI, HbIIα, and HbIIβ.

CONCLUSION

This study provides transcriptome responses to multiple environmental stressors in blood clams and would provide interesting hints for further studies.

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.

Divergence history and hydrothermal vent adaptation of decapod crustaceans: A mitogenomic perspective

Decapod crustaceans, such as alvinocaridid shrimps, bythograeid crabs and galatheid squat lobsters are important fauna in the hydrothermal vents and have well adapted to hydrothermal vent environments. In this study, eighteen mitochondrial genomes (mitogenomes) of hydrothermal vent decapods were used to explore the evolutionary history and their adaptation to the hydrothermal vent habitats. BI and ML algorithms produced consistent phylogeny for Decapoda. The phylogenetic relationship revealed more evolved positions for all the hydrothermal vent groups, indicating they migrated from non-vent environments, instead of the remnants of ancient hydrothermal vent species, which support the extinction/repopulation hypothesis. The divergence time estimation on the Alvinocarididae, Bythograeidae and Galatheoidea nodes are located at 75.20, 56.44 and 47.41–50.43 Ma, respectively, which refers to the Late Cretaceous origin of alvinocaridid shrimps and the Early Tertiary origin of bythograeid crabs and galatheid squat lobsters. These origin stories are thought to associate with the global deep-water anoxic/dysoxic events. Total eleven positively selected sites were detected in the mitochondrial OXPHOS genes of three lineages of hydrothermal vent decapods, suggesting a link between hydrothermal vent adaption and OXPHOS molecular biology in decapods. This study adds to the understanding of the link between mitogenome evolution and ecological adaptation to hydrothermal vent habitats in decapods.

Gene expression study in Bathymodiolus azoricus populations from three North Atlantic hydrothermal vent sites

The deep-sea hydrothermal vents are known as harsh environments, abundant in animal diversity surrounded by fluids with specific physiological and chemical composition. Bathymodiolus azoricus mussels are endemic species dwelling at hydrothermal vent sites and at distinct depth ranges. Mussels from Menez Gwen (MG), Lucky Strike (LS), Rainbow (Rb) were collected at 800 m, 1730 m and 2310 m depths respectively, along the Mid-Atlantic Ridge. Five different tissues including gill, digestive gland, mantle, adductor muscle and foot from MG, LS and Rb mussels were selected for gene expression analyses by qPCR. 30 genes were tested to investigate the level of immune and apoptotic gene expression among B. azoricus populations. Statistical analyses confirmed tissue-specific gene expression differences among the five tissues. The digestive gland tissue showed a higher transcriptional activity characterized by an up-regulation of gene activities, contrary to what was assessed in the adductor muscle tissue. Five categories included recognition, signaling, transcription, effector and apoptotic genes were analyzed in this study. The majority of genes differed in levels of expression between MG/LS and LS/Rb in the digestive gland. Our findings suggest that gene expression profiles are inherent to the tissue analyzed, thus implying an immune tissue-specificity controlling defense responses across B. azoricus mussel body as a whole.

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.

Protein expression profiles in Bathymodiolus azoricus exposed to cadmium

Proteomic changes in the "gill-bacteria complex" of the hydrothermal vent mussel B. azoricus exposed to cadmium in pressurized chambers ((Incubateurs Pressurises pour l'Observation en Culture d'Animaux Marins Profonds - IPOCAMP) were analyzed and compared with the non-exposed control group. 2-D Fluorescence Difference Gel Electrophoresis (2D-DIGE) showed that less than 1.5% of the proteome of mussels and symbiotic bacteria were affected by a short-term (24 h) Cd exposure. Twelve proteins of the more abundant differentially expressed proteins of which six were up-regulated and six were down-regulated were excised, digested and identified by mass spectrometry. The identified proteins included structural proteins (actin/actin like proteins), metabolic proteins (calreticulin/calnexin, peptidyl-prolyl cis-trans isomerase, aminotransferase class-III, electron transfer flavoprotein, proteasome, alpha-subunit and carbonic anhydrase) and stress response proteins (chaperone protein htpG, selenium-binding protein and glutathione transferases). All differently expressed proteins are tightly connected to Cd exposure and are affected by oxidative stress. It was also demonstrated that B. azoricus was well adapted to Cd contamination therefore B. azoricus from hydrothermal vent areas may be considered a good bioindicator.

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.

RNA-Seq Analysis for Assessing the Early Response to DSP Toxins in Mytilus galloprovincialis Digestive Gland and Gill

The harmful effects of diarrhetic shellfish poisoning (DSP) toxins on mammalian cell lines have been widely assessed. Studies in bivalves suggest that mussels display a resistance to the cytogenotoxic effects of DSP toxins. Further, it seems that the bigger the exposure, the more resistant mussels become. To elucidate the early genetic response of mussels against these toxins, the digestive gland and the gill transcriptomes of Mytilus galloprovincialis after Prorocentrum lima exposure (100,000 cells/L, 48 h) were de novo assembled based on the sequencing of 8 cDNA libraries obtained using an Illumina HiSeq 2000 platform. The assembly provided 95,702 contigs. A total of 2286 and 4523 differentially expressed transcripts were obtained in the digestive gland and the gill, respectively, indicating tissue-specific transcriptome responses. These transcripts were annotated and functionally enriched, showing 44 and 60 significant Pfam families in the digestive gland and the gill, respectively. Quantitative PCR (qPCR) was performed to validate the differential expression patterns of several genes related to lipid and carbohydrate metabolism, energy production, genome integrity and defense, suggesting their participation in the protective mechanism. This work provides knowledge of the early response against DSP toxins in the mussel M. galloprovincialis and useful information for further research on the molecular mechanisms of the bivalve resistance to these toxins.

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.

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.

Identification of distinct haemocyte populations from the freshwater bivalves swan mussel (Anodonta cygnea) and duck mussel (Anodonta anatina) using wheat-germ agglutinin

Haemocytes play a major role in molluscs immunity. Functional studies are, however, impaired by limited available experimental tools to identify and sort distinct haemocyte populations. Therefore, using nonlethal methods, we aimed at evaluating whether lectin staining combined with flow cytometry could be used to distinguish circulating haemocyte populations from two freshwater bivalves of the family Unionidae, the duck mussel (Anodonta anatina (L., 1758)) and the swan mussel (Anodonta cygnea (L., 1758)). Based on classical classification, haemocytes were distinguished as granulocytes and hyalinocytes and cytological features were visualized using transmission microscopy and staining techniques. Size, granularity, viability, and surface staining using lectins as specific probes were analysed by flow cytometry and fluorescence microscopy. The microscopic proportions of granulocytes and hyalinocytes significantly differed, being of 70% and 30% for A. cygnea and of 85% and 15% for A. anatina, respectively. Two haemocyte populations were sorted by flow cytometry based on size and granularity and confirmed as granulocytes and hyalinocytes. Interestingly, two different granulocyte populations could be further discriminated in A. cygnea according to their binding affinity to wheat-germ agglutinin (WGA), whereas granulocytes of A. anatina all stained similarly. Our results show that WGA labelling combined with flow cytometry can be used to better discriminate Anodonta haemocyte populations and obtain purified populations for functional studies.

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.

Comparative immunogenomics of molluscs

Comparative immunology, studying both vertebrates and invertebrates, provided the earliest descriptions of phagocytosis as a general immune mechanism. However, the large scale of animal diversity challenges all-inclusive investigations and the field of immunology has developed by mostly emphasizing study of a few vertebrate species. In addressing the lack of comprehensive understanding of animal immunity, especially that of invertebrates, comparative immunology helps toward management of invertebrates that are food sources, agricultural pests, pathogens, or transmit diseases, and helps interpret the evolution of animal immunity. Initial studies showed that the Mollusca (second largest animal phylum), and invertebrates in general, possess innate defenses but lack the lymphocytic immune system that characterizes vertebrate immunology. Recognizing the reality of both common and taxon-specific immune features, and applying up-to-date cell and molecular research capabilities, in-depth studies of a select number of bivalve and gastropod species continue to reveal novel aspects of molluscan immunity. The genomics era heralded a new stage of comparative immunology; large-scale efforts yielded an initial set of full molluscan genome sequences that is available for analyses of full complements of immune genes and regulatory sequences. Next-generation sequencing (NGS), due to lower cost and effort required, allows individual researchers to generate large sequence datasets for growing numbers of molluscs. RNAseq provides expression profiles that enable discovery of immune genes and genome sequences reveal distribution and diversity of immune factors across molluscan phylogeny. Although computational de novo sequence assembly will benefit from continued development and automated annotation may require some experimental validation, NGS is a powerful tool for comparative immunology, especially increasing coverage of the extensive molluscan diversity. To date, immunogenomics revealed new levels of complexity of molluscan defense by indicating sequence heterogeneity in individual snails and bivalves, and members of expanded immune gene families are expressed differentially to generate pathogen-specific defense responses.

Genome sequence of the sulfur-oxidizing Bathymodiolus thermophilus gill endosymbiont

Bathymodiolus thermophilus, a mytilid mussel inhabiting the deep-sea hydrothermal vents of the East Pacific Rise, lives in symbiosis with chemosynthetic Gammaproteobacteria within its gills. The intracellular symbiont population synthesizes nutrients for the bivalve host using the reduced sulfur compounds emanating from the vents as energy source. As the symbiont is uncultured, comprehensive and detailed insights into its metabolism and its interactions with the host can only be obtained from culture-independent approaches such as genomics and proteomics. In this study, we report the first draft genome sequence of the sulfur-oxidizing symbiont of B. thermophilus, here tentatively named Candidatus Thioglobus thermophilus. The draft genome (3.1 Mb) harbors 3045 protein-coding genes. It revealed pathways for the use of sulfide and thiosulfate as energy sources and encodes the Calvin-Benson-Bassham cycle for CO₂ fixation. Enzymes required for the synthesis of the tricarboxylic acid cycle intermediates oxaloacetate and succinate were absent, suggesting that these intermediates may be substituted by metabolites from external sources. We also detected a repertoire of genes associated with cell surface adhesion, bacteriotoxicity and phage immunity, which may perform symbiosis-specific roles in the B. thermophilus symbiosis.

Short-chain alkanes fuel mussel and sponge Cycloclasticus symbionts from deep-sea gas and oil seeps

Cycloclasticus bacteria are ubiquitous in oil-rich regions of the ocean and are known for their ability to degrade polycyclic aromatic hydrocarbons (PAHs). In this study, we describe Cycloclasticus that have established a symbiosis with Bathymodiolus heckerae mussels and poecilosclerid sponges from asphalt-rich, deep-sea oil seeps at Campeche Knolls in the southern Gulf of Mexico. Genomic and transcriptomic analyses revealed that in contrast to all known Cycloclasticus, the symbiotic Cycloclasticus appeared to lack the genes needed for PAH degradation. Instead, these symbionts use propane and other short-chain alkanes such as ethane and butane as carbon and energy sources, thus expanding the limited range of substrates known to power chemosynthetic symbioses. Analyses of short-chain alkanes in the environment of the Campeche Knolls symbioses revealed that these are present at high concentrations (in the µM to mM range). Comparative genomic analyses revealed high similarities between the genes used by the symbiotic Cycloclasticus to degrade short-chain alkanes and those of free-living Cycloclasticus that bloomed during the Deepwater Horizon (DWH) oil spill. Our results indicate that the metabolic versatility of bacteria within the Cycloclasticus clade is higher than previously assumed, and highlight the expanded role of these keystone species in the degradation of marine hydrocarbons.

Exploring the molecular basis of adaptive evolution in hydrothermal vent crab Austinograea alayseae by transcriptome analysis

Elucidating the genetic mechanisms of adaptation to the hydrothermal vent in organisms at genomic level is significant for understanding the adaptive evolution process in the extreme environment. We performed RNA-seq on four different tissues of a vent crab species, Austinograea alayseae, producing 725,461 unigenes and 134,489 annotated genes. Genes related to sensory, circadian rhythm, hormone, hypoxia stress, metal detoxification and immunity were identified. It was noted that in the degenerated eyestalk, transcription of phototransduction related genes which are important for retinal function was greatly reduced; three crucial neuropeptide hormones, one molt-inhibiting and two crustacean hyperglycemic hormone precursors were characterized with conserved domains; hypoxia-inducible factor 1 and two novel isoforms of metallothioneins in the vent crabs were discovered. An analysis of 6,932 orthologs among three crabs A. alayseae, Portunus trituberculutus and Eriocheir sinensis revealed 19 positive selected genes (PSGs). Most of the PSGs were involved in immune responses, such as crustins and anti-lipopolysaccharide factor, suggesting their function in the adaptation to environment. The characterization of the first vent crab transcriptome provides abundant resources for genetic and evolutionary studies of this species, and paves the way for further investigation of vent adaptation process in crabs.

Adaptation and evolution of deep-sea scale worms (Annelida: Polynoidae): insights from transcriptome comparison with a shallow-water species

Polynoid scale worms (Polynoidae, Annelida) invaded deep-sea chemosynthesis-based ecosystems approximately 60 million years ago, but little is known about their genetic adaptation to the extreme deep-sea environment. In this study, we reported the first two transcriptomes of deep-sea polynoids (Branchipolynoe pettiboneae, Lepidonotopodium sp.) and compared them with the transcriptome of a shallow-water polynoid (Harmothoe imbricata). We determined codon and amino acid usage, positive selected genes, highly expressed genes and putative duplicated genes. Transcriptome assembly produced 98,806 to 225,709 contigs in the three species. There were more positively charged amino acids (i.e., histidine and arginine) and less negatively charged amino acids (i.e., aspartic acid and glutamic acid) in the deep-sea species. There were 120 genes showing clear evidence of positive selection. Among the 10% most highly expressed genes, there were more hemoglobin genes with high expression levels in both deep-sea species. The duplicated genes related to DNA recombination and metabolism, and gene expression were only enriched in deep-sea species. Deep-sea scale worms adopted two strategies of adaptation to hypoxia in the chemosynthesis-based habitats (i.e., rapid evolution of tetra-domain hemoglobin in Branchipolynoe or high expression of single-domain hemoglobin in Lepidonotopodium sp.).

Adaptation to deep-sea chemosynthetic environments as revealed by mussel genomes

Hydrothermal vents and methane seeps are extreme deep-sea ecosystems that support dense populations of specialized macro-benthos such as mussels. But the lack of genome information hinders the understanding of the adaptation of these animals to such inhospitable environments. Here we report the genomes of a deep-sea vent/seep mussel (Bathymodiolus platifrons) and a shallow-water mussel (Modiolus philippinarum). Phylogenetic analysis shows that these mussel species diverged approximately 110.4 million years ago. Many gene families, especially those for stabilizing protein structures and removing toxic substances from cells, are highly expanded in B. platifrons, indicating adaptation to extreme environmental conditions. The innate immune system of B. platifrons is considerably more complex than that of other lophotrochozoan species, including M. philippinarum, with substantial expansion and high expression levels of gene families that are related to immune recognition, endocytosis and caspase-mediated apoptosis in the gill, revealing presumed genetic adaptation of the deep-sea mussel to the presence of its chemoautotrophic endosymbionts. A follow-up metaproteomic analysis of the gill of B. platifrons shows methanotrophy, assimilatory sulfate reduction and ammonia metabolic pathways in the symbionts, providing energy and nutrients, which allow the host to thrive. Our study of the genomic composition allowing symbiosis in extremophile molluscs gives wider insights into the mechanisms of symbiosis in other organisms such as deep-sea tubeworms and giant clams.

Metabolic and physiological interdependencies in the Bathymodiolus azoricus symbiosis

The hydrothermal vent mussel Bathymodiolus azoricus lives in an intimate symbiosis with two types of chemosynthetic Gammaproteobacteria in its gills: a sulfur oxidizer and a methane oxidizer. Despite numerous investigations over the last decades, the degree of interdependence between the three symbiotic partners, their individual metabolic contributions, as well as the mechanism of carbon transfer from the symbionts to the host are poorly understood. We used a combination of proteomics and genomics to investigate the physiology and metabolism of the individual symbiotic partners. Our study revealed that key metabolic functions are most likely accomplished jointly by B. azoricus and its symbionts: (1) CO₂ is pre-concentrated by the host for carbon fixation by the sulfur-oxidizing symbiont, and (2) the host replenishes essential biosynthetic TCA cycle intermediates for the sulfur-oxidizing symbiont. In return (3), the sulfur oxidizer may compensate for the host's putative deficiency in amino acid and cofactor biosynthesis. We also identified numerous ‘symbiosis-specific' host proteins by comparing symbiont-containing and symbiont-free host tissues and symbiont fractions. These proteins included a large complement of host digestive enzymes in the gill that are likely involved in symbiont digestion and carbon transfer from the symbionts to the host.

Differential Gene Expression during Larval Metamorphic Development in the Pearl Oyster, Pinctada fucata, Based on Transcriptome Analysis

P. fucata experiences a series of transformations in appearance, from swimming larvae to sessile juveniles, during which significant changes in gene expression likely occur. Thus, P. fucata could be an ideal model in which to study the molecular mechanisms of larval metamorphosis during development in invertebrates. To study the molecular driving force behind metamorphic development in larvae of P. fucata, transcriptomes of five larval stages (trochophore, D-shape, umbonal, eyespots, and spats) were sequenced using an Illumina HiSeq™ 2000 system and assembled and characterized with the transcripts of six tissues. As a result, a total of 174,126 unique transcripts were assembled and 60,999 were annotated. The number of unigenes varied among the five larval stages. Expression profiles were distinctly different between trochophore, D-shape, umbonal, eyespots, and spats larvae. As a result, 29 expression trends were sorted, of which eight were significant. Among others, 80 development-related, differentially expressed unigenes (DEGs) were identified, of which the majority were homeobox-containing genes. Most DEGs occurred among trochophore, D-shaped, and UES (umbonal, eyespots, and spats) larvae as verified by qPCR. Principal component analysis (PCA) also revealed significant differences in expression among trochophore, D-shaped, and UES larvae with ten transcripts identified but no matching annotations.

Transcriptome Analysis of the Sydney Rock Oyster, Saccostrea glomerata: Insights into Molluscan Immunity

BACKGROUND

Oysters have important ecological functions in their natural environment, acting as global carbon sinks and improving water quality by removing excess nutrients from the water column. During their life-time oysters are exposed to a variety of pathogens that can cause severe mortality in a range of oyster species. Environmental stressors encountered in their habitat can increase the susceptibility of oysters to these pathogens and in general have been shown to impact on oyster immunity, making immune parameters expressed in these marine animals an important research topic.

RESULTS

Paired-end Illumina high throughput sequencing of six S. glomerata tissues exposed to different environmental stressors resulted in a total of 484,121,702 paired-end reads. When reads and assembled transcripts were compared to the C. gigas genome, an overall low level of similarity at the nucleotide level, but a relatively high similarity at the protein level was observed. Examination of the tissue expression pattern showed that some transcripts coding for cathepsins, heat shock proteins and antioxidant proteins were exclusively expressed in the haemolymph of S. glomerata, suggesting a role in innate immunity. Furthermore, analysis of the S. glomerata ORFs showed a wide range of genes potentially involved in innate immunity, from pattern recognition receptors, components of the Toll-like signalling and apoptosis pathways to a complex antioxidant defence mechanism.

CONCLUSIONS

This is the first large scale RNA-Seq study carried out in S. glomerata, showing the complex network of innate immune components that exist in this species. The results confirmed that many of the innate immune system components observed in mammals are also conserved in oysters; however, some, such as the TLR adaptors MAL, TRIF and TRAM are either missing or have been modified significantly. The components identified in this study could help explain the oysters' natural resilience against pathogenic microorganisms encountered in their natural environment.

The Shell of the Invasive Bivalve Species Dreissena polymorpha: Biochemical, Elemental and Textural Investigations

The zebra mussel Dreissena polymorpha is a well-established invasive model organism. Although extensively used in environmental sciences, virtually nothing is known of the molecular process of its shell calcification. By describing the microstructure, geochemistry and biochemistry/proteomics of the shell, the present study aims at promoting this species as a model organism in biomineralization studies, in order to establish a bridge with ecotoxicology, while sketching evolutionary conclusions. The shell of D. polymorpha exhibits the classical crossed-lamellar/complex crossed lamellar combination found in several heterodont bivalves, in addition to an external thin layer, the characteristics of which differ from what was described in earlier publication. We show that the shell selectively concentrates some heavy metals, in particular uranium, which predisposes D. polymorpha to local bioremediation of this pollutant. We establish the biochemical signature of the shell matrix, demonstrating that it interacts with the in vitro precipitation of calcium carbonate and inhibits calcium carbonate crystal formation, but these two properties are not strongly expressed. This matrix, although overall weakly glycosylated, contains a set of putatively calcium-binding proteins and a set of acidic sulphated proteins. 2D-gels reveal more than fifty proteins, twenty of which we identify by MS-MS analysis. We tentatively link the shell protein profile of D. polymorpha and the peculiar recent evolution of this invasive species of Ponto-Caspian origin, which has spread all across Europe in the last three centuries.

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.

Shotgun proteomics to unravel marine mussel (Mytilus edulis) response to long-term exposure to low salinity and propranolol in a Baltic Sea microcosm

Pharmaceuticals, among them the β-adrenoceptor blocker propranolol, are an important group of environmental contaminants reported in European waters. Laboratory exposure to pharmaceuticals on marine species has been performed without considering the input of the ecosystem flow. To unravel the ecosystem response to long-term exposure to propranolol we have performed long-term exposure to propranolol and low salinity in microcosms. We applied shotgun proteomic analysis to gills of Mytilus edulis from those Baltic Sea microcosms and identified 2071 proteins with a proteogenomic strategy. The proteome profiling patterns from the 587 highly reproductive proteins among groups define salinity as a key factor in the mussel's response to propranolol. Exposure at low salinity drives molecular mechanisms of adaptation based on a decrease in the abundance of several cytoskeletal proteins, signalling and intracellular membrane trafficking pathway combined with a response towards the maintenance of transcription and translation. The exposure to propranolol combined with low salinity modulates the expression of structural proteins including cilia functions and decreases the expression of membrane protein transporters. This study reinforces the environment concerns of the impact of low salinity in combination with anthropogenic pollutants and anticipates critical physiological conditions for the survival of the blue mussel in the northern areas.

BIOLOGICAL SIGNIFICANCE

Applying shotgun proteomic analysis to M. edulis gills samples from a long-term microcosm exposure to propranolol and following a proteogenomic identification strategy, we have identified 2071 proteins. The proteomic analysis unrevealed which molecular mechanisms drive the adaptation to low salinity stress and how salinity modulates the effects of exposure to propranolol. These results reinforce the idea of the impact of low salinity in combination with anthropogenic pollutants and anticipate critical physiological condition.

Comparative physiological, biochemical, and molecular thermal stress response profiles for two Unionid freshwater mussel species

Freshwater mussels, aquatic keystone species, are in global decline. Long life spans, sedentary lifestyles, and unique reproductive strategies involving obligate parasitic stages make unionid freshwater mussels particularly sensitive to environmental perturbations resulting from global climate change. A greater understanding of the mechanisms by which closely related species differ in their response to thermal challenge is critical for successful conservation and management practices. As such, both an acute heat shock and a chronic warming simulation were conducted in order to evaluate responses between hypothesized thermally tolerant, Villosa lienosa, and thermally sensitive, Villosa nebulosa, freshwater mussels in response to predicted thermal warming. Multiple biological responses were quantified, including mortality, condition index, growth rates, glycogen and triglyceride content, and candidate gene expression. During acute heat shock, both species upregulated HSP90 and HSP70, though V. lienosa showed consistently greater transcript levels during upregulation. This pattern was consistent during the chronic warming simulation, with V. nebulosa showing greater induction of HSP60. Chronic warming stimulated increases in condition index for V. nebulosa, however declines in growth rates during a recovery period were observed with no concurrent tissue glycogen levels changes. This contrasts with V. lienosa, where tissue glycogen significantly increased during chronic warming, though no response was observed for condition index or growth rates. These differences might indicate disparate thermal stress response mechanisms correlated with metabolic demands and resource utilization. These biological differences could thus be a factor influencing current ranges and these two species‘ future ability to cope with persistent warming in their native habitats.

Allatostatin-type A, kisspeptin and galanin GPCRs and putative ligands as candidate regulatory factors of mantle function

Allatostatin-type A (AST-A), kisspeptin (KISS) and galanin (GAL) G-protein coupled receptor (GPCR) systems share a common ancestral origin in arthropods and the vertebrates where they regulate metabolism and reproduction. The molluscs are the second most diverse phylum in the animal kingdom, they occupy an important phylogenetic position, and their genome is more similar to deuterostomes than the arthropods and nematodes and thus they are good models for studies of gene family evolution and function. This mini-review intends to extend the current knowledge about AST-A, KISS and GAL GPCR system evolution and their putative function in the mollusc mantle. Comparative evolutionary analysis of the target GPCR systems was established by identifying homologues in genomes and tissue transcriptome datasets available for molluscs and comparing them to those of other metazoan systems. Studies in arthropods have revealed the existence of the AST-A system but the loss of homologues of the KISS and GAL systems. Homologues of the insect AST-AR and vertebrate KISSR genes were found in molluscs but putative GALR genes were absent. Receptor gene number suggested that members of this family have suffered lineage specific evolution during the molluscan radiation. In molluscs, orthologues of the insect AST-A peptides were not identified but buccalin peptides that are structurally related were identified and are putative receptor agonists. The identification of AST-AR and KISSR genes in molluscs strengthens the hypotheses that in metazoans members of the AST-AR subfamily share evolutionary proximity with KISSRs. The variable number of receptors and large repertoire of buccalin peptides may be indicative of the functional diversity of the AST-AR/KISSR systems in molluscs. The identification of AST-A and KISS receptors and ligands in the mantle transcriptome indicates that in molluscs they may have acquired a novel function and may play a role in shell development or sensory detection in the mantle.

Vibrio diabolicus challenge in Bathymodiolus azoricus populations from Menez Gwen and Lucky Strike hydrothermal vent sites

Menez Gwen (MG) and Lucky Strike (LS) deep-sea hydrothermal vents are located at 850 m and 1730 m depths respectively and support chemosynthesis-based ecosystems partially differing in heavy metal concentration, temperature range, and faunistic composition. The successfully adapted deep-sea vent mussel Bathymodiolus azoricus is found at both vent locations. In such inhospitable environments survival strategies rely on the establishment of bacteria-vent animal symbiosis In spite of the toxic nature of deep-sea vents, the problem of microbial threat and the need for immunity exist in B. azoricus. This study aims at investigating the immune system of B. azoricus from MG and LS populations by comparing immune gene expressions profiles using the deep-sea vent-related Vibrio diabolicus. Expression of nineteen immune genes was analyzed from gill, digestive gland and mantle tissues upon 3 h, 12 h and 24 h V. diabolicus challenges. Based on quantitative-Polymerase Chain Reaction (qPCR) significant gene expression differences were found among MG and LS populations and challenge times MG mussels revealed that gill and digestive gland gene expression levels were remarkably higher than those from LS mussels. Expression of Carcinolectin, Serpin-2, SRCR, IRGs, RTK, TLR2, NF-κB, HSP70 and Ferritin genes was greater in MG than LS mussels. In contrast, mantle tissue from LS mussels revealed the highest peak of expression at 24 h for most genes analyzed. The activation of immune signaling pathways demonstrated that gene expression profiles are distinct between the two mussel populations. These differences may possibly ensue from intrinsic immune transcriptional activities upon which host responses are modulated in presence of V. diabolicus. mRNA transcript variations were assessed during 24 h acclimatization taking into account the partial depuration to which mussels were subjected to. Additionally, gene expression differences may reflect still accountable effects from the presence of vent remaining microfluidic environments within the tissues analyzed.

High-throughput transcriptome sequencing of the cold seep mussel Bathymodiolus platifrons

Bathymodiolid mussels dominate hydrothermal vents, cold methane/sulfide-hydrocarbon seeps, and other sites of organic enrichment. Here, we aimed to explore the innate immune system and detoxification mechanism of the deep sea mussel Bathymodiolus platifrons collected from a methane seep in the South China Sea. We sequenced the transcriptome of the mussels’ gill, foot and mantle tissues and generated a transcriptomic database containing 96,683 transcript sequences. Based on GO and KEGG annotations, we reported transcripts that were related to the innate immune system, heavy metal detoxification and sulfide metabolic genes. Our in-depth analysis on the isoforms of peptidoglycan recognition protein (PGRP) that have different cellular location and potentially differential selectivity towards peptidoglycan (PGN) from gram-positive and gram-negative bacteria were differentially expressed in different tissues. We also reported a potentially novel form of metallothionein and the production of phytochelatin in B. platifrons, which has not been reported in any of its coastal relative Mytilus mussel species. Overall, the present study provided new insights into heavy metal and sulfide metabolism in B. platifrons and can be served as the basis for future molecular studies on host-symbiont interactions in cold seep mussels.

RNA-Seq in Mytilus galloprovincialis: comparative transcriptomics and expression profiles among different tissues

Background

The Mediterranean mussel (Mytilus galloprovincialis) is a cosmopolitan, cultured bivalve with worldwide commercial and ecological importance. However, there is a qualitative and quantitative lack of knowledge of the molecular mechanisms involved in the physiology and immune response of this mollusc. In order to start filling this gap, we have studied the transcriptome of mantle, muscle and gills from naïve Mediterranean mussels and hemocytes exposed to distinct stimuli.

Results

A total of 393,316 million raw RNA-Seq reads were obtained and assembled into 151,320 non-redundant transcripts with an average length of 570 bp. Only 55 % of the transcripts were shared across all tissues. Hemocyte and gill transcriptomes shared 60 % of the transcripts while mantle and muscle transcriptomes were most similar, with 77 % shared transcripts. Stimulated hemocytes showed abundant defense and immune-related proteins, in particular, an extremely high amount of antimicrobial peptides. Gills expressed many transcripts assigned to both structure and recognition of non-self patterns, while in mantle many transcripts were related to reproduction and shell formation. Moreover, this tissue presented additional and interesting hematopoietic, antifungal and sensorial functions. Finally, muscle expressed many myofibril and calcium-related proteins and was found to be unexpectedly associated with defense functions. In addition, many metabolic routes related to cancer were represented.

Conclusions

Our analyses indicate that whereas the transcriptomes of these four tissues have characteristic expression profiles in agreement with their biological structures and expected functions, tissue-specific transcriptomes reveal a complex and specialized functions.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1817-5) contains supplementary material, which is available to authorized users.

An updated molecular basis for mussel immunity

Non-self recognition with the consequent tolerance or immune reaction is a crucial process to succeed as living organisms. At the same time the interactions between host species and their microbiome, including potential pathogens and parasites, significantly contribute to animal life diversity. Marine filter-feeding bivalves, mussels in particular, can survive also in heavily anthropized coastal waters despite being constantly surrounded by microorganisms. Based on the first outline of the Mytilus galloprovincialis immunome dated 2011, the continuously growing transcript data and the recent release of a draft mussel genome, we explored the available sequence data and scientific literature to reinforce our knowledge on the main gene-encoded elements of the mussel immune responses, from the pathogen recognition to its clearance. We carefully investigated molecules specialized in the sensing and targeting of potential aggressors, expected to show greater molecular diversification, and outlined, whenever relevant, the interconnected cascades of the intracellular signal transduction. Aiming to explore the diversity of extracellular, membrane-bound and intracellular pattern recognition receptors in mussel, we updated a highly complex immune system, comprising molecules which are described here in detail for the first time (e.g. NOD-like receptors) or which had only been partially characterized in bivalves (e.g. RIG-like receptors). Overall, our comparative sequence analysis supported the identification of over 70 novel full-length immunity-related transcripts in M. galloprovincialis. Nevertheless, the multiplicity of gene functions relevant to immunity, the involvement of part of them in other vital processes, and also the lack of a refined mussel genome make this work still not-exhaustive and support the development of more specific studies.

Key metabolic pathways involved in xenobiotic biotransformation and stress responses revealed by transcriptomics of the mangrove oyster Crassostrea brasiliana

The Brazilian oyster Crassostrea brasiliana was challenged to three common environmental contaminants: phenanthrene, diesel fuel water-accommodated fraction (WAF) and domestic sewage. Total RNA was extracted from the gill and digestive gland, and cDNA libraries were sequenced using the 454 FLX platform. The assembled transcriptome resulted in ̃20,000 contigs, which were annotated to produce the first de novo transcriptome for C. brasiliana. Sequences were screened to identify genes potentially involved in the biotransformation of xenobiotics and associated antioxidant defence mechanisms. These gene families included those of the cytochrome P450 (CYP450), 70kDa heat shock, antioxidants, such as glutathione S-transferase, superoxide dismutase, catalase and also multi-drug resistance proteins. Analysis showed that the massive expansion of the CYP450 and HSP70 family due to gene duplication identified in the Crassostrea gigas genome also occurred in C. brasiliana, suggesting these processes form the base of the Crassostrea lineage. Preliminary expression analyses revealed several candidates biomarker genes that were up-regulated during each of the three treatments, suggesting the potential for environmental monitoring.