Combined metabolome and proteome analysis of the mantle tissue from Pacific oyster Crassostrea gigas exposed to elevated pCO2

The molecular response of Mytilus coruscus mantle to shell damage under acute acidified sea water revealed by iTRAQ based quantitative proteomic analysis

Mytilus coruscus is an economically important marine bivalve that lives in estuarine sea areas with seasonal coastal acidification and frequently suffers shell injury in the natural environment. However, the molecular responses and biochemical properties of Mytilus under these conditions are not fully understood. In the present study, we employed tandem mass spectrometry combined with isobaric tagging to identify differentially expressed proteins in the mantle tissue of M. coruscus under different short-term treatments, including shell-complete mussels raised in normal seawater (pH 8.1), shell-damaged mussels raised in normal seawater (pH 8.1), and acidified seawater (pH 7.4). A total of 2694 proteins were identified in the mantle, and analysis of their relative abundance from the three different treatments revealed alterations in the proteins involved in immune regulation, oxidation-reduction processes, protein folding and processing, energy provision, and cytoskeleton. The results obtained by quantitative proteomic analysis of the mantle allowed us to delineate the molecular strategies adopted by M. coruscus in the shell repair process in acidified environments, including an increase in proteins involved in oxidation-reduction processes, protein processing, and cell growth at the expense of proteins involved in immune capacity and energy metabolism. SIGNIFICANCE: The impact of global ocean acidification on calcifying organisms has become a major ecological and environmental problem in the world. Mytilus coruscus is an economically important marine bivalve living in estuary sea area with seasonal coastal acidification, and frequently suffering shell injury in natural environment. Molecular responses of M coruscus under the shell damage and acute acidification is still largely unknown. For this reason, iTRAQ based quantitative proteomic and histological analysis of the mantle from M. coruscus under shell damage and acute acidification were performed, for revealing the proteomic response and possible adaptation mechanism of Mytilus under combined shell damage and acidified sea water, and understanding how the mussel mantle implement a shell-repair process under acidified sea water. Our study provides important data for understanding the shell repair process and proteomic response of Mytilus under ocean acidification, and providing insights into potential adaptation of mussels to future global change.

Molecular responses to ocean acidification in an Antarctic bivalve and an ascidian

Southern Ocean organisms are considered particularly vulnerable to Ocean acidification (OA), as they inhabit cold waters where calcite-aragonite saturation states are naturally low. It is also generally assumed that OA would affect calcifying animals more than non-calcifying animals. In this context, we aimed to study the impact of reduced pH on both types of species: the ascidian Cnemidocarpa verrucosa sp. A, and the bivalve Aequiyoldia eightsii, from an Antarctic fjord. We used gene expression profiling and enzyme activity to study the responses of these two Antarctic benthic species to OA. We report the results of an experiment lasting 66 days, comparing the molecular mechanisms underlying responses under two pCO2 treatments (ambient and elevated pCO2). We observed 224 up-regulated and 111 down-regulated genes (FC ≥ 2; p-value ≤ 0.05) in the ascidian. In particular, the decrease in pH caused an upregulation of genes involved in the immune system and antioxidant response. While fewer differentially expressed (DE) genes were observed in the infaunal bivalve, 34 genes were up-regulated, and 69 genes were downregulated (FC ≥ 2; p-value ≤ 0.05) in response to OA. We found downregulated genes involved in the oxidoreductase pathway (such as glucose dehydrogenase and trimethyl lysine dioxygenase), while the heat shock protein 70 was up-regulated. This work addresses the effect of OA in two common, widely distributed Antarctic species, showing striking results. Our major finding highlights the impact of OA on the non-calcifying species, a result that differ from the general trend, which describes a higher impact on calcifying species. This calls for discussion of potential effects on non-calcifying species, such as ascidians, a diverse and abundant group that form extended three-dimensional clusters in shallow waters and shelf areas in the Southern Ocean.

Comparative Single-Cell Transcriptomics Reveals Novel Genes Involved in Bivalve Embryonic Shell Formation and Questions Ontogenetic Homology of Molluscan Shell Types

Mollusks are known for their highly diverse repertoire of body plans that often includes external armor in form of mineralized hardparts. Representatives of the Conchifera, one of the two major lineages that comprises taxa which originated from a uni-shelled ancestor (Monoplacophora, Gastropoda, Cephalopoda, Scaphopoda, Bivalvia), are particularly relevant regarding the evolution of mollusk shells. Previous studies have found that the shell matrix of the adult shell (teleoconch) is rapidly evolving and that the gene set involved in shell formation is highly taxon-specific. However, detailed annotation of genes expressed in tissues involved in the formation of the embryonic shell (protoconch I) or the larval shell (protoconch II) are currently lacking. Here, we analyzed the genetic toolbox involved in embryonic and larval shell formation in the quagga mussel Dreissena rostriformis using single cell RNA sequencing. We found significant differences in genes expressed during embryonic and larval shell secretion, calling into question ontogenetic homology of these transitory bivalve shell types. Further ortholog comparisons throughout Metazoa indicates that a common genetic biomineralization toolbox, that was secondarily co-opted into molluscan shell formation, was already present in the last common metazoan ancestor. Genes included are engrailed, carbonic anhydrase, and tyrosinase homologs. However, we found that 25% of the genes expressed in the embryonic shell field of D. rostriformis lack an ortholog match with any other metazoan. This indicates that not only adult but also embryonic mollusk shells may be fast-evolving structures. We raise the question as to what degree, and on which taxonomic level, the gene complement involved in conchiferan protoconch formation may be lineage-specific or conserved across taxa.

Effects of tralopyril on histological, biochemical and molecular impacts in Pacific oyster, Crassostrea gigas

Recently, the toxic effects of tralopyril, as a new antifouling biocide, on aquatic organisms have aroused widespread attention about the potential toxicity. However, the mechanism of tralopyril on marine mollusks has not been elaborated clearly. In this study, the histological, biochemical and molecular impacts of tralopyril on adult Crassostrea gigas were investigated. The results indicated that the 96 h LC50 of tralopyril to adult Crassostrea gigas was 911 μg/L. After exposure to tralopyril (0, 40, 80 and 160 μg/L) for 6 days, the mantle mucus secretion coverage ratio of Crassostrea gigas was increased with a dose-dependent pattern. Catalase (CAT) activity was significantly increased, amylase (AMS) activity, acid phosphatase (ACP) activity and calcium ion (Ca²⁺) concentration significantly decreased. Meanwhile, integrated biomarker responses (IBR) index suggested that higher concentrations of tralopyril caused severer damage to Crassostrea gigas. In addition, the mRNA expression levels of biomineralization related genes in the mantle were significantly upregulated. Collectively, this study firstly revealed the histological, biochemical and molecular impacts of tralopyril exposure on adult Crassostrea gigas, which provided new insights for understanding the toxicity of tralopyril in marine mollusks.

Repeated exposure to simulated marine heatwaves enhances the thermal tolerance in pearl oysters

In an era of unprecedented climate change, marine heatwaves (MHWs) are projected to increase in frequency, intensity, and duration, severely affecting marine organisms and fisheries and causing profound ecological and socioeconomic impacts. However, very little is known about effects of MHWs on ecologically and economically important bivalve species. Here, we investigate how pearl oysters, Pinctada maxima (Jameson), respond to MHWs under various simulated scenarios. Acute 3-day exposure to MHWs, mimicked by increasing the ambient seawater temperature from 24°C to 28°C, 32°C, and 36°C, respectively, induced significant changes (initially sustained increase and late decrease) in activities of antioxidant enzymes (GSH-Px, SOD, CAT, MDA, and T-AOC) and biomineralizaiton-related enzymes (AKP and ACP). Likewise, energy-metabolizing enzymes (NKA, CMA, and T-ATP) showed remarkable acute responses, with significant increases occurring at the beginning and end of thermal exposure. Following repeated exposure to MHWs, without exception, acute responses of assayed enzymes to MHWs were significantly alleviated, implying that pearl oysters have the ability to implement more efficient and less costly energy-utilizing strategies to compensate for thermal stress induced physiological interferences. These findings indicate that marine bivalves can respond plastically and acclimate rapidly to MHWs and pave the way for predicting the fate of this important taxonomic groups in rapidly changing oceans.

A proteomic analysis of the effect of ocean acidification on the haemocyte proteome of the South African abalone Haliotis midae

As a result of increasing CO₂ emissions and the prevalence of global climate change, ocean acidification (OA) is becoming more pervasive, affecting many trophic levels, particularly those that rely on succinctly balanced ocean chemistry. This ultimately threatens community structures, as well as the future sustainability of the fishing/aquaculture industry. Understanding the molecular stress response of key organisms will aid in predicting their future survivability under changing environmental conditions. This study sought to elucidate the molecular stress response of the South African abalone, Haliotis midae, an understudied organism with high economic value, utilising a high throughput iTRAQ-based proteomics methodology. Adult abalone were exposed to control (pH 7.9) and experimental (pH 7.5) conditions for 12, 72 and 168 h, following which protein was isolated from sampled haemocytes and subsequently processed. iTRAQ-labelled peptides were analysed using mass spectrometry, while an array of bioinformatics tools was utilised for analysing the proteomic data. COG analysis identified "Cytoskeleton", "Translation, ribosomal structure and biogenesis", "Post-translational modification, protein turnover, chaperones", and "Intracellular trafficking, secretion and vesicular transport" to be the most enriched functional classes, while statistical analysis identified a total of 33 up-regulated and 23 down-regulated effectors of OA stress in abalone. Several of the up-regulated proteins that were identified function in central metabolism (ENO1, PGK, DUOX1, GPD2), the stress/immune response (CAMKI, HSPA5/GRP78, MAPKI), and cytoskeleton, protein sorting and signal transduction (IQGAP1, MYO9B, TLN1, RDX, TCP-1/CCT, SNX6, CHMP1a, VPS13a). Protein-protein interactions were predicted using STRING DB, Cytoscape and Ingenuity Pathway Analysis, providing a model of the effects of OA on the H. midae haemocyte proteome. The data indicated that H. midae underwent a metabolic shift under OA conditions to utilize more energy-efficient mechanisms of ATP generation, while attempts at restoring haemocyte stabilisation and homeostasis were reflected by up-regulation of oxidative stress and cytoskeletal proteins. Our results support other molluscan studies that report a complex array of overlapping functions of both the stress and immune response systems. This interplay of the mounted stress and immune response is maintained and observed through the up-regulation of proteins involved in protein synthesis and turnover, as well as intracellular signalling and transport. The data presented in this study highlight the value of employing sensitive and robust -omics technologies for assessing the effects of changing environmental conditions on marine organisms.

Temperature but not ocean acidification affects energy metabolism and enzyme activities in the blue mussel, Mytilus edulis

In mosaic marine habitats, such as intertidal zones, ocean acidification (OA) is exacerbated by high variability of pH, temperature, and biological CO2 production. The nonlinear interactions among these drivers can be context-specific and their effect on organisms in these habitats remains largely unknown, warranting further investigation.We were particularly interested in Mytilus edulis (the blue mussel) from intertidal zones of the Gulf of Maine (GOM), USA, for this study. GOM is a hot spot of global climate change (average sea surface temperature (SST) increasing by >0.2°C/year) with >60% decline in mussel population over the past 40 years.Here, we utilize bioenergetic underpinnings to identify limits of stress tolerance in M. edulis from GOM exposed to warming and OA. We have measured whole-organism oxygen consumption rates and metabolic biomarkers in mussels exposed to control and elevated temperatures (10 vs. 15°C, respectively) and current and moderately elevated P CO2 levels (~400 vs. 800 µatm, respectively).Our study demonstrates that adult M. edulis from GOM are metabolically resilient to the moderate OA scenario but responsive to warming as seen in changes in metabolic rate, energy reserves (total lipids), metabolite profiles (glucose and osmolyte dimethyl amine), and enzyme activities (carbonic anhydrase and calcium ATPase).Our results are in agreement with recent literature that OA scenarios for the next 100-300 years do not affect this species, possibly as a consequence of maintaining its in vivo acid-base balance.

Metabolic response of Scapharca subcrenata to heat stress using GC/MS-based metabolomics

Marine mollusks are commonly subjected to heat stress. To evaluate the effects of heat stress on the physiological metabolism of the ark shell Scapharca subcrenata, clams were exposed to different high temperatures (24, 28 and 32 °C) for 72 h. The oxygen consumption and ammonia excretion rates were measured at 2, 12, 24, 48 and 72 h. The results indicated that the metabolic rates of the ark shell significantly increased with increasing heat stress, accompanied by mortalities in response to prolonged exposure. A metabolomics approach based on gas chromatography coupled with mass spectrometry was further applied to assess the changes of metabolites in the mantle of the ark shell at 32 °C. Moreover, multivariate and pathway analyses were conducted for the different metabolites. The results showed that the heat stress caused changes in energy metabolism, amino acid metabolism, osmotic regulation, carbohydrate metabolism and lipid metabolism through different metabolic pathways. These results are consistent with the significant changes of oxygen consumption rate and ammonia excretion rate. The present study contributes to the understanding of the impacts of heat stress on intertidal bivalves and elucidates the relationship between individual-level responses and underlying molecular metabolic dynamics.

Effects of acidification on the proteome during early development of Babylonia areolata

Increases in atmospheric CO2 partial pressure have lowered seawater pH in marine ecosystems, a process called ocean acidification (OA). The effects of OA during the critical stages of larval development may have disastrous consequences for some marine species, including Babylonia areolata (Link 1807), a commercially important sea snail in China and South East Asia. To investigate how OA affects the proteome of Babylonia areolata, here we used label-free proteomics to study protein changes in response to acidified (pH 7.6) or ambient seawater (pH 8.1) during three larvae developmental stages of B. areolata, namely, the veliger larvae before attachment (E1), veliger larvae after attachment (E2), and carnivorous juvenile snail (E3). In total, we identified 720 proteins. This result suggested that acidification seriously affects late veliger stage after attachment (E2). Further examination of the roles of differentially expressed proteins, which include glutaredoxin, heat-shock protein 70, thioredoxin, catalase, cytochrome-c-oxidase, peroxiredoxin 6, troponin T, CaM kinase II alpha, proteasome subunit N3 and cathepsin L, will be important for understanding the molecular mechanisms underlying pH reduction.

Metabolic responses to elevated pCO2 in the gills of the Pacific oyster (Crassostrea gigas) using a GC-TOF-MS-based metabolomics approach

Rising atmospheric carbon dioxide (CO₂), primarily from anthropogenic emissions, are resulting in increasing absorption of CO₂ by the oceans, leading to a decline in oceanic pH in a process known as ocean acidification (OA). There is a growing body of evidence demonstrating the potential effect of OA on the energetics/physiology and consequently life-history traits of commensally important marine organisms. However, despite this little is known of how fundamental metabolic pathways that underpin changes in organismal physiology are affected by OA. Consequently, a gas chromatography time-of-flight mass spectrometry (GC-TOF-MS) based metabolic profiling approach was applied to examine the metabolic responses of Crassostrea gigas to elevated pCO₂ levels, under otherwise natural field conditions. Oysters were exposed natural environmental pCO₂ (~625.40 μatm) and elevated pCO₂ (~1432.94 μatm) levels for 30 days. Results indicated that 36 differential metabolites were identified. Differential metabolites were mapped in the Kyoto Encyclopedia of Genes and Genomes (KEGG) database to search for the related metabolic pathways. Pathway enrichment analysis indicates that alanine, aspartate and glutamate metabolism and glycine, serine and threonine metabolism were the most statistically enriched pathways. Further analysis suggested that elevated pCO₂ disturb the TCA cycle via succinate accumulation and C. gigas most likely adjust their energy metabolic via alanine and GABA accumulation accordingly to cope with elevated pCO₂. These findings provide an understanding of the molecular mechanisms involved in modulating C. gigas metabolism under elevated pCO₂.

Metabolome responses of the sea cucumber Apostichopus japonicus to multiple environmental stresses: Heat and hypoxia

Economically important marine organisms face severe environmental challenges, such as high temperature and low dissolved oxygen, from global climate change. Adverse environmental factors impact the survival and growth of economically important marine organisms, thereby negatively influencing the aquaculture industry. However, little is known about the responses of sea cucumbers to combined environmental co-stressors till now. In this study, ultra-performance liquid chromatography (UPLC) was utilized to obtain metabolic profiles of sea cucumbers. Changes in the concentrations of 84, 68, and 417 metabolites related to the responses of sea cucumbers to heat (26 °C), hypoxia (2 mg/L) and the combined stress, respectively, were observed and analyzed. Representative biomarkers were discussed in detail, including deltaline, fusarin C, halichondrin B and rapanone. The concentration of metabolites involved in the regulation of energy metabolism, including amino acid, carbohydrate and lipid metabolism were significantly changed, and the tricarboxylic acid (TCA)-cycle was significantly altered under heat plus hypoxia. We interpreted these changes partly as an adaptation mechanism in response to environmental stress. Based on the decreased accumulation of glutamine, we hypothesized that heat stress is the main factor that interferes with the process of glutamic acid-glutamine metabolism. The present study showed that combined environmental stressors have a more extensive impact on the metabolites of the respiratory tree in sea cucumbers than single stress. These results would facilitate further development of the sea cucumber as an echinoderm model to study mechanisms of response to adverse environments, as well as to help advance knowledge of the adaptation of marine organisms to global climate change.

Gender-specific metabolic responses in gonad of mussel Perna viridis to triazophos

Triazophos, as a lipophilic organophosphate pesticide, displays higher bioaccumulation in the gonads of shellfish. To study the reproductive toxicity of triazophos, we applied metabolomics to characterize the gender-specific metabolic responses in mussel Perna viridis exposed to triazophos. Metabolites were differently altered by triazophos in ovaries of mussel at different concentrations and time intervals, while basically similar metabolic response patterns were observed in male mussels at the two tested concentrations after exposure for 24 and 48h. The significant changes of metabolites in ovaries of mussel exhibited the disturbances in energy metabolism and osmotic regulation, while in male samples triazophos only affected the energy metabolism. Moreover, glycine, sn-glycero-3-phosphocholine, ethanol, aspartate, etc. exhibited consistent variation tendency in both male and female individuals. While the changes of homarine, betaine, taurine, hypotaurine, malonate, β-alanine, succinate, and choline showed obviously gender-specific responses. Overall, this study confirmed the gender-specific responses in gonad of P. viridis to triazophos exposure.

Transcriptomic profiling of adaptive responses to ocean acidification

Some populations of marine organisms appear to have inherent tolerance or the capacity for acclimation to stressful environmental conditions, including those associated with climate change. Sydney rock oysters from the B2 breeding line exhibit resilience to ocean acidification (OA) at the physiological level. To understand the molecular basis of this physiological resilience, we analysed the gill transcriptome of B2 oysters that had been exposed to near-future projected ocean pH over two consecutive generations. Our results suggest that the distinctive performance of B2 oysters in the face of OA is mediated by the selective expression of genes involved in multiple cellular processes. Subsequent high-throughput qPCR revealed that some of these transcriptional changes are exclusive to B2 oysters and so may be associated with their resilience to OA. The intracellular processes mediated by the differentially abundant genes primarily involve control of the cell cycle and maintenance of cellular homeostasis. These changes may enable B2 oysters to prevent apoptosis resulting from oxidative damage or to alleviate the effects of apoptosis through regulation of the cell cycle. Comparative analysis of the OA conditioning effects across sequential generations supported the contention that B2 and wild-type oysters have different trajectories of changing gene expression and responding to OA. Our findings reveal the broad set of molecular processes underlying transgenerational conditioning and potential resilience to OA in a marine calcifier. Identifying the mechanisms of stress resilience can uncover the intracellular basis for these organisms to survive and thrive in a rapidly changing ocean.

A Preliminary Study on the Pattern, the Physiological Bases and the Molecular Mechanism of the Adductor Muscle Scar Pigmentation in Pacific Oyster Crassostrea gigas

The melanin pigmentation of the adductor muscle scar and the outer surface of the shell are among attractive features and their pigmentation patterns and mechanism still remains unknown in the Pacific oyster Crassostrea gigas. To study these pigmentation patterns, the colors of the adductor muscle scar vs. the outer surface of the shell on the same side were compared. No relevance was found between the colors of the adductor muscle scars and the corresponding outer surface of the shells, suggesting that their pigmentation processes were independent. Interestingly, a relationship between the color of the adductor muscle scars and the dried soft-body weight of Pacific oysters was found, which could be explained by the high hydroxyl free radical scavenging capacity of the muscle attached to the black adductor muscle scar. After the transcriptomes of pigmented and unpigmented adductor muscles and mantles were studied by RNAseq and compared, it was found that the retinol metabolism pathway were likely to be involved in melanin deposition on the adductor muscle scar and the outer surface of the shell, and that the different members of the tyrosinase or Cytochrome P450 gene families could play a role in the independent pigmentation of different organs.

Contrasting impacts of ocean acidification and warming on the molecular responses of CO2-resilient oysters

BACKGROUND

This study characterises the molecular processes altered by both elevated CO₂ and increasing temperature in oysters. Differences in resilience of marine organisms against the environmental stressors associated with climate change will have significant implications for the sustainability of coastal ecosystems worldwide. Some evidence suggests that climate change resilience can differ between populations within a species. B2 oysters represent a unique genetic resource because of their capacity to better withstand the impacts of elevated CO₂ at the physiological level, compared to non-selected oysters from the same species (Saccostrea glomerata). Here, we used proteomic and transcriptomic analysis of gill tissue to evaluate whether the differential response of B2 oysters to elevated CO₂ also extends to increased temperature.

RESULTS

Substantial and distinctive effects on protein concentrations and gene expression were evident among B2 oysters responding to elevated CO₂ or elevated temperature. The combination of both stressors also altered oyster gill proteomes and gene expression. However, the impacts of elevated CO₂ and temperature were not additive or synergistic, and may be antagonistic.

CONCLUSIONS

The data suggest that the simultaneous exposure of CO₂-resilient oysters to near-future projected ocean pH and temperature results in complex changes in molecular processes in order to prevent stress-induced cellular damage. The differential response of B2 oysters to the combined stressors also indicates that the addition of thermal stress may impair the resilience of these oysters to decreased pH. Overall, this study reveals the intracellular mechanisms that might enable marine calcifiers to endure the emergent, adverse seawater conditions resulting from climate change.

Effects of elevated carbon dioxide on contraction force and proteome composition of sea urchin tube feet

This study examined how contraction force and protein profiles of the tube feet of the sea urchin (Pseudocentrotus depressus) were affected when acclimated to 400 (control), 2000 and 10,000μatm CO₂ for 48days. Acclimation to higher CO₂ conditions significantly reduced contraction force of the tube feet. Two-dimensional gel electrophoresis showed that eight spots changed in protein volume: six up-regulated and two down-regulated. Using matrix-assisted laser desorption/ionization-quadrupole ion trap-time of flight mass spectrometry, three up-regulated spots (tubulin beta chain, tropomyosin fragment, and actin N-terminal fragment) and two down-regulated spots (actin C-terminal fragment and myosin light chain) were identified. One possible interpretation of the results is that elevated CO₂ weakened contraction of the tube feet muscle through an alteration of proteome composition, mainly associated with post-translational processing/proteolysis of muscle-related proteins.

Triomics Analysis of Imatinib-Treated Myeloma Cells Connects Kinase Inhibition to RNA Processing and Decreased Lipid Biosynthesis

The combination of metabolomics, lipidomics, and phosphoproteomics that incorporates triple stable isotope labeling by amino acids in cell culture (SILAC) protein labeling, as well as (13)C in vivo metabolite labeling, was demonstrated on BCR-ABL-positive H929 multiple myeloma cells. From 11 880 phosphorylation sites, we confirm that H929 cells are primarily signaling through the BCR-ABL-ERK pathway, and we show that imatinib treatment not only downregulates phosphosites in this pathway but also upregulates phosphosites on proteins involved in RNA expression. Metabolomics analyses reveal that BCR-ABL-ERK signaling in H929 cells drives the pentose phosphate pathway (PPP) and RNA biosynthesis, where pathway inhibition via imatinib results in marked PPP impairment and an accumulation of RNA nucleotides and negative regulation of mRNA. Lipidomics data also show an overall reduction in lipid biosynthesis and fatty acid incorporation with a significant decrease in lysophospholipids. RNA immunoprecipitation studies confirm that RNA degradation is inhibited with short imatinib treatment and transcription is inhibited upon long imatinib treatment, validating the triomics results. These data show the utility of combining mass spectrometry-based "-omics" technologies and reveals that kinase inhibitors may not only downregulate phosphorylation of their targets but also induce metabolic events via increased phosphorylation of other cellular components.