Molecular identification of differentially regulated genes in the hydrothermal-vent species Bathymodiolus thermophilus and Paralvinella pandorae in response to temperature

Genetic adaptations of marine invertebrates to hydrothermal vent habitats

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

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.

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.

A dynamic history of admixture from Mediterranean and Carpathian glacial refugia drives genomic diversity in the bank vole

Understanding the historical contributions of differing glacial refugia is key to evaluating the roles of microevolutionary forces, such as isolation, introgression, and selection in shaping genomic diversity in present-day populations. In Europe, where both Mediterranean and extra-Mediterranean (e.g., Carpathian) refugia of the bank vole (Clethrionomys glareolus) have been identified, mtDNA indicates that extra-Mediterranean refugia were the main source of colonization across the species range, while Mediterranean peninsulas harbor isolated, endemic lineages. Here, we critically evaluate this hypothesis using previously generated genomic data (>6,000 SNPs) for over 800 voles, focusing on genomic contributions to bank voles in central Europe, a key geographic area in considering range-wide colonization. The results provide clear evidence that both extra-Mediterranean (Carpathian) and Mediterranean (Spanish, Calabrian, and Balkan) refugia contributed to the ancestry and genomic diversity of bank vole populations across Europe. Few strong barriers to dispersal and frequent admixture events in central Europe have led to a prominent mid-latitude peak in genomic diversity. Although the genomic contribution of the centrally located Carpathian refugium predominates, populations in different parts of Europe have admixed origins from Mediterranean (28%-47%) and the Carpathian (53%-72%) sources. We suggest that the admixture from Mediterranean refugia may have provisioned adaptive southern alleles to more northern populations, facilitating the end-glacial spread of the admixed populations and contributing to increased bank vole diversity in central Europe. This study adds critical details to the complex end-glacial colonization history of this well-studied organism and underscores the importance of genomic data in phylogeographic interpretation.

Differential gene expression indicates modulated responses to chronic and intermittent hypoxia in corallivorous fireworms (Hermodice carunculata)

Climate models predict an increase in extent, frequency, and duration of marine hypoxia events in the twenty first century. A better understanding of organismal responses to hypoxia in individual species is a crucial step for predicting ecosystem responses. We experimentally subjected a common invertebrate, the bearded fireworm (Hermodice carunculata) to two levels of chronic hypoxia and, in a separate experiment, to intermittent hypoxia. We found components of the conserved hypoxia-inducible factor (HIF) pathway and show a modulated response to hypoxia depending on the severity of hypoxic stress: under mild hypoxia, only the HIF-1α subunit is upregulated, while expression of the other subunit, aryl hydrocarbon nuclear translator, only increases significantly at more severe hypoxia levels. The chronic trials revealed down-regulation of genes related to cell adhesion, transport, development and heme-binding, and up-regulation of genes related to glycolysis, oxygen binding, cell differentiation, digestive and reproductive function. The intermittent hypoxia trials revealed an upregulation of heme transporter activity during hypoxia, and our time series analysis characterized nine clusters of genes with similar expression patterns. Our findings suggest that H. carunculata is likely to tolerate, and be resilient to, predicted future hypoxia conditions.

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.

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

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

Proteomic responses to metal-induced oxidative stress in hydrothermal vent-living mussels, Bathymodiolus sp., on the Southwest Indian Ridge

Bathymodiolin mussels are amongst the dominant fauna occupying hydrothermal vent ecosystems throughout the World's oceans. This subfamily inhabits a highly ephemeral and variable environment, where exceptionally high concentrations of reduced sulphur species and heavy metals necessitate adaptation of specialised detoxification mechanisms. Whilst cellular responses to common anthropogenic pollutants are well-studied in shallow-water species, they remain limited in deep-sea vent fauna. Bathymodiolus sp. were sampled from two newly-discovered vent sites on the Southwest Indian Ridge (Tiamat and Knuckers Gaff) by the remotely operated vehicle (ROV) Kiel 6000 during the RRS James Cook cruise, JC 067 in November 2011. Here, we use redox proteomics to investigate the effects of tissue metal accumulation on protein expression and thiol oxidation in gill. Following 2D PAGE, we demonstrate a significant difference in intensity in 30 protein spots in this organ between the two vent sites out of 205 matched spots. We also see significant variations in thiol oxidation in 15 spots, out of 143 matched. At Tiamat, 23 protein spots are up-regulated compared to Knuckers Gaff and we identify 5 of these with important roles in metabolism, cell structure, stress response, and redox homeostasis. We suggest that increased metal exposure triggers changes in the proteome, regulating tissue uptake. This is evident both between vent sites and across a chemical gradient within the Knuckers Gaff vent site. Our findings highlight the importance of proteomic plasticity in successful adaptation to the spatially and temporally fluctuating chemical environments that are characteristic of hydrothermal vent habitats.

First evidence of protein profile alteration due to the main cocaine metabolite (benzoylecgonine) in a freshwater biological model

Illicit drugs represent not only a great social problem but are also considered an environmental problem because their use and, often, abuse release large amounts of parent compounds, and especially their metabolites, into freshwaters. One of the most commonly used drugs is cocaine, which is the second most prevalent drug in Europe (accounting for almost 30% of all cocaine users worldwide). Cocaine is rapidly metabolised in humans to benzoylecgonine (35-54%), ecgonine methyl ester (32-49%) and norcocaine (5%), which are eliminated in the urine and are only partially removed by wastewater treatment plants (WWTPs). Because no studies have previously been carried out to evaluate the possible risks due to cocaine and its metabolites in non-target organisms, we applied a multi-disciplinary approach to investigate the possible environmental risk related to benzoylecgonine (BE), the main metabolite of cocaine. Previous studies carried out by means of a biomarker suite and the redox-proteomic approach showed an imbalance of anti-oxidant enzyme activities and several genotoxic effects to be caused by environmental BE concentrations in the freshwater bivalve Zebra mussel (Dreissena polymorpha). This report presents the results obtained in the last step of this study, based on a proteomics analysis. We analysed the protein expression profile in the gills of Zebra mussels exposed to two different concentrations (0.5 and 1 μg/L) of BE for 14 days through 2-DE and mass spectrometry analysis (RP-UPLC ESI-LTQ-Orbitrap). Our results highlight significant changes in some proteins in gill cells whose functions are crucial for overall metabolism. In particular, we detected a probable effect of BE on calcium homeostasis and a consequent imbalance of oxidative stress, as verified for vertebrates.

Seasonal variability of antioxidant biomarkers in mud crabs (Scylla serrata)

Studies on oxidative stress (OS) in crustacea are widely used as ecotoxicological indices to assess the environment risk produced by the impact of several stressor and pollutants. In the present study, effects of seasonality on OS physiology markers such as antioxidant enzymes (superoxide dismutase, catalase, glutathione peroxidase and glutathione reductase), small antioxidant molecules (ascorbic acid and reduced glutathione), oxidative stress indices (lipid peroxidation, protein carbonylation and hydrogen peroxide) and total antioxidant capacity in hepatopancreas, gills and abdominal muscle of adult mud crab Scylla serrata, sampled from Chilika lagoon of India, were determined in winter, summer and rainy seasons. Results indicate that variations in enzymatic and non-enzymatic antioxidants with relation to season were not only tissue specific but also were gender specific. The levels of OS parameters were higher in hepatopancreas in comparison to gills and abdominal muscle of the crabs in all seasons. OS indices in tissues of the crabs were mainly higher in summer season when temperature and salinity of the lagoon were high with low oxygen content. Although OS was lower in winter season and moderate in rainy season in tissues of male crabs, it was higher in gills and hepatopancreas of females in rainy season. Correlation analyses between hydrological parameters of the lagoon (temperature, salinity and dissolved oxygen content) and OS physiology parameters in tissues of crabs suggest that abiotic factors influence the levels of antioxidant enzymes and, thereby the OS status in a tissue and sex specific manner. Collectively, the results of the present work suggest that further investigation is warranted before using OS parameters in S. serrata as biomarkers to monitor estuarine environment as these are influenced by gender, tissue and season.

Cloning and characterization of a Δ9-desaturase gene of the Antarctic fish Chionodraco hamatus and Trematomus bernacchii

Chionodraco hamatus and Trematomus bernacchii are perciforms, members of the fish suborder Notothenioidei that live in the Antarctic Ocean and experience very cold and persistent environmental temperature. These fish have biochemical and molecular features that allow them to live at these extreme cold temperatures. Fine tuning of the level of unsaturated fatty acids content in membrane is a key mechanism of living organisms to adapt to cold and high temperatures. Desaturases are key enzymes that synthesize unsaturated fatty acyl-CoAs from saturated fatty acids. We cloned and sequenced a Δ(9)-desaturase gene and its cDNA of C. hamatus, and the cDNA of T. bernacchii. The coded proteins are virtually identical and share homology to other Δ(9)-desaturase fish sequences. These proteins contain, in the first trans-membrane domain, two cysteine residues that may form a disulfur bond present in the corresponding membrane region of Δ(9)-desaturase proteins of other Antarctic fish but not in Eleginops maclovinus that experiences higher environmental temperatures and in all other Δ(9)-desaturase genes of mammals present in data bases. C. hamatus Δ(9)-desaturase gene complements a Saccharomyces cerevisiae mutant lacking Δ(9)-desaturase (Ole1) gene. Analysis of sequence homology of the trans-membrane domains of Δ(9)-desaturase and the cytoplasmic region of the same proteins of Antarctic fish, non-Antarctic fish and mammals suggest that the significant differences found in the homologous sequences of the first trans-membrane domain may be due to the specific lipid content of their membrane.

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

Background

The deep-sea hydrothermal vent mussel Bathymodiolus azoricus harbors thiotrophic and methanotrophic symbiotic bacteria in its gills. While the symbiotic relationship between this hydrothermal mussel and these chemoautotrophic bacteria has been described, the molecular processes involved in the cross-talking between symbionts and host, in the maintenance of the symbiois, in the influence of environmental parameters on gene expression, and in transcriptome variation across individuals remain poorly understood. In an attempt to understand how, and to what extent, this double symbiosis affects host gene expression, we used a transcriptomic approach to identify genes potentially regulated by symbiont characteristics, environmental conditions or both. This study was done on mussels from two contrasting populations.

Results

Subtractive libraries allowed the identification of about 1000 genes putatively regulated by symbiosis and/or environmental factors. Microarray analysis showed that 120 genes (3.5% of all genes) were differentially expressed between the Menez Gwen (MG) and Rainbow (Rb) vent fields. The total number of regulated genes in mussels harboring a high versus a low symbiont content did not differ significantly. With regard to the impact of symbiont content, only 1% of all genes were regulated by thiotrophic (SOX) and methanotrophic (MOX) bacteria content in MG mussels whereas 5.6% were regulated in mussels collected at Rb. MOX symbionts also impacted a higher proportion of genes than SOX in both vent fields. When host transcriptome expression was analyzed with respect to symbiont gene expression, it was related to symbiont quantity in each field.

Conclusions

Our study has produced a preliminary description of a transcriptomic response in a hydrothermal vent mussel host of both thiotrophic and methanotrophic symbiotic bacteria. This model can help to identify genes involved in the maintenance of symbiosis or regulated by environmental parameters. Our results provide evidence of symbiont effect on transcriptome regulation, with differences related to type of symbiont, even though the relative percentage of genes involved remains limited. Differences observed between the vent site indicate that environment strongly influences transcriptome regulation and impacts both activity and relative abundance of each symbiont. Among all these genes, those participating in recognition, the immune system, oxidative stress, and energy metabolism constitute new promising targets for extended studies on symbiosis and the effect of environmental parameters on the symbiotic relationships in B. azoricus.

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

Background

Bathymodiolus azoricus is a deep-sea hydrothermal vent mussel found in association with large faunal communities living in chemosynthetic environments at the bottom of the sea floor near the Azores Islands. Investigation of the exceptional physiological reactions that vent mussels have adopted in their habitat, including responses to environmental microbes, remains a difficult challenge for deep-sea biologists. In an attempt to reveal genes potentially involved in the deep-sea mussel innate immunity we carried out a high-throughput sequence analysis of freshly collected B. azoricus transcriptome using gills tissues as the primary source of immune transcripts given its strategic role in filtering the surrounding waterborne potentially infectious microorganisms. Additionally, a substantial EST data set was produced and from which a comprehensive collection of genes coding for putative proteins was organized in a dedicated database, "DeepSeaVent" the first deep-sea vent animal transcriptome database based on the 454 pyrosequencing technology.

Results

A normalized cDNA library from gills tissue was sequenced in a full 454 GS-FLX run, producing 778,996 sequencing reads. Assembly of the high quality reads resulted in 75,407 contigs of which 3,071 were singletons. A total of 39,425 transcripts were conceptually translated into amino-sequences of which 22,023 matched known proteins in the NCBI non-redundant protein database, 15,839 revealed conserved protein domains through InterPro functional classification and 9,584 were assigned with Gene Ontology terms. Queries conducted within the database enabled the identification of genes putatively involved in immune and inflammatory reactions which had not been previously evidenced in the vent mussel. Their physical counterpart was confirmed by semi-quantitative quantitative Reverse-Transcription-Polymerase Chain Reactions (RT-PCR) and their RNA transcription level by quantitative PCR (qPCR) experiments.

Conclusions

We have established the first tissue transcriptional analysis of a deep-sea hydrothermal vent animal and generated a searchable catalog of genes that provides a direct method of identifying and retrieving vast numbers of novel coding sequences which can be applied in gene expression profiling experiments from a non-conventional model organism. This provides the most comprehensive sequence resource for identifying novel genes currently available for a deep-sea vent organism, in particular, genes putatively involved in immune and inflammatory reactions in vent mussels.

The characterization of the B. azoricus transcriptome will facilitate research into biological processes underlying physiological adaptations to hydrothermal vent environments and will provide a basis for expanding our understanding of genes putatively involved in adaptations processes during post-capture long term acclimatization experiments, at "sea-level" conditions, using B. azoricus as a model organism.

Identification of differentially expressed genes in the hydrothermal vent shrimp Rimicaris exoculata exposed to heat stress

The deep-sea vent shrimp Rimicaris exoculata dominates the vagile megafauna at most vent sites along the Mid-Atlantic Ridge. This shrimp swarms around the hot end of the hydrothermal biotope where temperature can exceed its critical maximal temperature (33-38.5 ± 2°C). It may therefore be subjected to a thermal regime that is assumed to be stressful for animals. In this study, we used a global transcriptomic approach by constructing suppression subtractive hybridization cDNA libraries in order to identify specific up- and down-regulated genes in R. exoculata exposed to a severe heat stress (1h at 30°C). A total of 218 sequences representing potentially highly expressed genes in thermally stressed shrimp were obtained. Expression of 11 genes involved in various cell functions was quantified in control and heat shocked specimens using real-time PCR. Differential expression was observed for some specific genes such as mannose receptor C1, metalloprotease, histone H1, and hemocyanin with a strong up-regulation of several genes encoding heat shock proteins. These results suggest that R. exoculata is affected at both cellular and molecular levels by sustained exposure at 30°C. The sequenced ESTs presented here will provide an excellent basis for future thermal stress studies on deep-sea vent fauna.

Global depression in gene expression as a response to rapid thermal changes in vent mussels

Hydrothermal vent mussels belonging to the genus Bathymodiolus are distributed worldwide and dominate communities at shallow Atlantic hydrothermal sites. While organisms inhabiting coastal ecosystems are subjected to predictable oscillations of physical and chemical variables owing to tidal cycles, the vent mussels sustain pronounced temperature changes over short periods of time, correlated to the alternation of oxic/anoxic phases. In this context, we focused on the short-term adaptive response of mussels to temperature change at a molecular level. The mRNA expression of 23 genes involved in various cell functions of the vent mussel Bathymodiolus azoricus was followed after heat shocks for either 30 or 120 min, at 25 and 30 degrees C over a 48 h recovery period at 5 degrees C. Mussels were genotyped at 10 enzyme loci to explore a relationship between natural genetic variation, gene expression and temperature adaptation. Results indicate that the mussel response to increasing temperature is a depression in gene expression, such a response being genotypically correlated at least for the Pgm-1 locus. This suggests that an increase in temperature could be a signal triggering anaerobiosis for B. azoricus or this latter alternatively behaves more like a 'cold' stenotherm species, an attribute more related to its phylogenetic history, a cold seeps/wood fall origin.