Adaptations to the mudflat: Insights from physiological and transcriptional responses to thermal stress in a burrowing bivalve Sinonovacula constricta

Histology, fatty acid composition, antioxidant and glycolipid metabolism, and transcriptome analyses of the acute cold stress response in Phoxinus lagowskii

Water temperature is a crucial environmental factor that significantly affects the physiological and biochemical processes of fish. Due to the occurrence of cold events in aquaculture, it is imperative to investigate how fish respond to cold stress. This study aims to uncover the mechanisms responds to acute cold stress by conducting a comprehensive analysis of the histomorphology, glycolipid metabolic and antioxidant enzymes, fatty acid composition and transcriptome at three temperatures (16 °C, 10 °C and 4 °C) in Phoxinus lagowskii. Our results showed that cold stress not damaged muscle microstructure but caused autophagy (at 10 °C). In addition, serum glucose (Glu) and triglycerides (TG) increased during cold stress. The activities of reactive oxygen species (ROS), superoxide dismutase (SOD), catalase (CAT), fructose phosphokinase (PFK), hexokinase (HK), pyruvate kinase (PK), and malondialdehyde (MDA) content in muscle were measured and analyzed. During cold stress, superoxide dismutase and catalase activities increased, reactive oxygen species content decreased. No significant difference in Glutathione peroxidase (GPx) activity, malondialdehyde and total cholesterol (T-CHO) contents among groups. Phosphokinase and pyruvate kinase activities decreased, and HK activity increased during cold stress. Our study resulted in the identification of a total of 25,400 genes, with 2524 genes showing differential expression across different temperature treatments. Furthermore, KEGG pathway indicated that some pathways upregulated during light cold stress (at 10 °C, including autophagy, and AMP-activated protein kinase (AMPK) signaling pathway. Additionally, circadian rhythm is among the most enriched pathways in genes up-regulated during severe cold stress (at 4 °C). Our findings offer valuable insights into how cold-water fish respond to cold stress.

Behavioural and physiological responses of juvenile geoduck (Panopea zelandica) following acute thermal stress

Climate extremes, such as heatwaves, are expected to become more intense and of longer duration in the near future. These climatic conditions may have a significant impact on the prospects of establishing a new aquaculture industry for the endemic New Zealand geoduck, Panopea zelandica. This study focused on characterising animal behaviour, haemocytes , and heat shock protein (HSP70 & HSP90) mRNA expression following exposure to elevated temperatures, such as those encountered during marine heatwaves around 20 °C and an extreme scenario of 25 °C, contrasted to an ambient temperature of 17 °C. After 24 h of heat challenge, P. zelandica were found to be significantly influenced by the thermal changes, as there were differences recorded in all the responses examined. With increasing temperatures, juvenile geoduck were observed to fully emerge from the sediment a behaviour that has not previously been quantified nor associated with stress in this species. The ability of P. zelandica juveniles to re-bury still warrants further investigation, as adults are unable to do so. Haemocyte analyses revealed an increase in the abundance of granulocytes, cellular aggregations, and size of these aggregations at the highest temperature exposure. Increased expression of the hsp70 gene in the haemolymph after exposure at 25 °C for 24 h was detected and attributed to attempts to mitigate protein denaturation caused by thermal stress. The inter-individual variability in the response of heat shock proteins recorded could aid in future selective breeding programs if it is reflected in net thermotolerance. P. zelandica shows great potential for growing in subtidal habitats around New Zealand, and this study highlights the importance of temperature considerations when selecting potential farm and reseeding locations.

Comparative Transcriptome Analyses Provide New Insights into the Evolution of Divergent Thermal Resistance in Two Eel Gobies

Adaptation to thermal conditions in tidal mudflats always involves tolerating frequent fluctuations and often extreme environmental temperatures. Regulation of gene expression plays a fundamental role in the evolution of these thermal adaptations. To identify the key gene regulatory networks associated with the thermal adaptation, we investigated the capability of cold tolerance, as well as the transcriptomic changes under cold stress in two mudflat inhabitants (Odontamblyopus lacepedii and O. rebecca) with contrasting latitude affinity. Our results revealed a remarkable divergent capacity of cold tolerance (CTmin: 0.61 °C vs. 9.57 °C) between the two gobies. Analysis of transcriptomic changes under cold stress unveiled 193 differentially expressed genes exhibiting similar expression profiles across all tissues and species, including several classic metabolic and circadian rhythm molecules such as ACOD and CIART that may represent the core cold response machinery in eel gobies. Meanwhile, some genes show a unique expression spectrum in the more cold-tolerant O. lacepedii suggesting their roles in the enhanced cold tolerance and hence the extreme thermal adaptations. In addition, a weighted gene co-expression network analysis (WGCNA) revealed a subset of metabolic hub genes including MYH11 and LIPT2 showing distinct down-regulation in O. lacepedii when exposed to cold stress which highlights the role of reduced energy consumption in the enhanced cold tolerance of eel gobies. These findings not only provide new insights into how mudflat teleosts could cope with cold stress and their potential evolutionary strategies for adapting to their thermal environment, but also have important implications for sound management and conservation of their fishery resources in a scenario of global climate warming in the marine realm.

Transcriptomic analysis revealed the dynamic response mechanism to acute ammonia exposure in the ivory shell, Babylonia areolata

The ivory shell (Babylonia areolata) is an economically important shellfish in tropical and subtropical regions, but its intensive culture and biological characteristic of hiding in the sandy substrate make it highly susceptible to ammonia stress. In this study, we investigated the dynamic changes in histopathology, oxidative stress, and transcriptome of the ivory shell at different time points under high concentration (60 mg/L) ammonia exposure. With prolonged exposure to stress, vacuoles appeared in the hepatopancreas while cell volume and intercellular space increased. The activities of superoxide dismutase (SOD) and catalase (CAT) decreased significantly under high concentrations of ammonia-induced stress while malondialdehyde (MDA) levels increased significantly. Integrated analysis of differentially expressed genes (DEGs), weighted gene co-expression network analysis (WGCNA), and quantitative real-time polymerase chain reaction (qRT-PCR) revealed that lipid transport primarily contributed to maintaining cellular homeostasis during the early stage of stress (6 and 12 h). Subsequently, a significant upregulation of oxidation-reduction reactions occurred at the middle stage (24 h), leading to oxidative stress. Finally, during the later stage (48 h), metabolic decomposition provided energy for survival maintenance. Additionally, lysosome and apoptosis were identified as potential key pathways in response to acute ammonia toxicity. Overall, our findings suggest that ivory shells can respond to acute ammonia toxicity via immune and antioxidant defense mechanisms but sustained high concentrations may cause irreversible damage. This study provides valuable insights into the response mechanism of mollusks towards ammonia and serves as a data reference for breeding ammonia-tolerant varieties of ivory shells.

Cardiac performance and heart gene network provide dynamic responses of bay scallop Argopecten irradians irradians exposure to marine heatwaves

The increased frequency of marine heat waves (MHWs) caused by global climate change is predicted to threaten the survival of economic bivalves, therefore having severely adverse effects on local ecological communities and aquaculture production. However, the study of scallops facing MHWs is still scarce, particularly in the scallop Argopecten irradians irradians, which has a significant share of "blue foods" in northern China. In the present study, bay scallop heart was selected to detect its cardiac performance, oxidative impairment and dynamic molecular responses, accompanied by assessing survival variations of individuals in the simulated scenario of MWHs (32 °C) with different time points (0 h, 6 h, 12 h, 24 h, 3 d, 6 d and 10 d). Notably, cardiac indices heart rate (HR), heart amplitude (HA), rate-amplitude product (RAP) and antioxidant enzyme activities superoxide dismutase (SOD) and catalase (CAT) all peaked at 24 h but sharply dropped on 3 d, coinciding with mortality. Transcriptome analysis revealed that the heart actively defended against heat stress at the acute stage (<24 h) via energy supply, misfolded proteins correction and enhanced signal transduction, whereas regulation of the defense response and apoptotic process combined with twice transcription initiation were the dominant responses at the chronic stage (3-10 d). In particular, HSP70 (heat shock protein 70), HSP90 and CALR (calreticulin) in the endoplasmic reticulum were identified as the hub genes (top 5 %) in the HR-associated module via WGCNA (weighted gene co-expression network analysis) trait-module analysis, followed by characterization of their family members and diverse expression patterns under heat exposure. Furthermore, RNAi-mediated knockdown of CALR expression (after 24 h) significantly weakened the thermotolerance of scallops, as evidenced by a drop of 1.31 °C in ABT (Arrhenius break temperature) between the siRNA-injected group and the control group. Our findings elucidated the dynamic molecular responses at the transcriptome level and verified the cardiac functions of CALR in bay scallops confronted with stimulated MHWs.

Impacts of marine heatwaves on byssus production in highly invasive fouling mussels

Marine heatwaves (MHWs) are projected to increase in their frequency, intensity, and duration, causing irreversible and catastrophic consequences for intertidal ecosystems around the world. The highly invasive fouling mussel, Arcuatula senhousia, can cause marked habitat alteration by constructing extremely intense byssal mats, devastating the biodiversity of many intertidal systems, yet very little is known about its fate under conditions of more frequent, hotter and longer MHWs. Here, we assessed impacts of two scenarios of MHWs (low-intensity with 4 °C rise of seawater temperature and high-intensity with 8 °C rise, respectively) on the byssal production of A. senhousia. Mussels exposed to low-intensity MHWs did not show any significant differences in the number, length and diameter of byssal threads, compared with those not thermally stressed. Under high-intensity scenario, the byssus production was significantly depressed, and byssal threads became fewer, shorter and finer, in line with significant decreases in cumulative length and volume. These findings provide a better understanding of responses of invasive fouling mussels such as A. senhousia to MHWs and make a leap forward in linking climate change and biological fouling in marine ecosystems.

Mitochondrial pyruvate carrier influences ganoderic acid biosynthesis in Ganoderma lucidum

Mitochondrial pyruvate carriers (MPCs), located in the inner membrane of mitochondria, are essential carriers for pyruvate to enter mitochondria. MPCs regulate a wide range of intracellular metabolic processes, such as glycolysis, the tricarboxylic acid cycle (TCA cycle), fatty acid metabolism, and amino acid metabolism. However, the metabolic regulation of MPCs in macrofungi is poorly studied. We studied the role of MPCs in Ganoderma lucidum (GlMPC) on ganoderic acid (GA) biosynthesis regulation in G. lucidum. In this study, we found that the mitochondrial/cytoplasmic ratio of pyruvate was downregulated about 75% in GlMPC1- and GlMPC2-silenced transformants compared with wild type (WT). In addition, the GA content was 17.72 mg/g and increased by approximately 50% in GlMPC1- and GlMPC2-silenced transformants compared with WT. By assaying the expression levels of three key enzymes and the enzyme activities of isocitrate dehydrogenase (IDH) and α-ketoglutarate dehydrogenase (α-KGDH) of the TCA cycle in GlMPC1- and GlMPC2-silenced transformants, it was found that the decrease in GlMPCs activity did not significantly downregulate the TCA cycle rate, and the enzyme activity of IDH increased by 44% compared with WT. We then verified that fatty acid β-oxidation (FAO) supplements the TCA cycle by detecting the expression levels of key enzymes involved in FAO. The results showed that compared with WT, the GA content was 1.14 mg/g and reduced by approximately 40% in co-silenced transformants. KEY POINTS: • GlMPCs affects the distribution of pyruvate between mitochondria and the cytoplasm. • Acetyl-CoA produced by FAO maintains the TCA cycle. • Acetyl-CoA produced by FAO promotes the accumulation of GA.

Comparative transcriptomic analysis revealed changes in multiple signaling pathways involved in protein degradation in the digestive gland of Mytilus coruscus during high-temperatures

As a result of global warming, the Mytilus coruscus living attached in the intertidal zone experience extreme and fluctuating changes in temperature, and extreme temperature changes are causing mass mortality of intertidal species. This study explores the transcriptional response of M. coruscus at different temperatures (18 °C, 26 °C, and 33 °C) and different times (0, 12, and 24 h) of action by analyzing the potential temperature of the intertidal zone. In response to high temperatures, several signaling pathways in M. coruscus, ribosome, endocytosis, endoplasmic reticulum stress, protein degradation, and lysosomes, interact to counter the adverse effects of high temperatures on protein homeostasis. Increased expression of key genes, including heat shock proteins (Hsp70, Hsp20, and Hsp110), Lysosome-associated membrane glycoprotein (LAMP), endoplasmic reticulum chaperone (BiP), and baculoviral IAP repeat-containing protein 7 (BIRC7), may further mitigate the effects of heat stress and delay mortality in M. coruscus. These results reveal changes in multiple signaling pathways involved in protein degradation during high-temperature stress, which will contribute to our overall understanding of the molecular mechanisms underlying the response of M. coruscus to high-temperature stress.

Transcriptomic responses reveal impaired physiological performance of the pearl oyster following repeated exposure to marine heatwaves

Marine heatwaves are predicted to become more intense and frequent in the future, possibly threatening the survival of marine organisms and devastating their communities. While recent evidence reveals the adaptability of marine organisms to heatwaves, substantially overlooked is whether they can also adjust to repeated heatwave exposure, which can occur in nature. By analysing transcriptome, we examined the fitness and recoverability of the pearl oyster (Pinctada maxima) after two consecutive heatwaves (24 °C to 32 °C for 3 days; recovery at 24 °C for 4 days). In the first heatwave, 331 differentially expressed genes (DEGs) were found, such as AGE-RAGE, MAPK, JAK-STAT, FoxO and mTOR. Despite the recovery after the first heatwave, 2511 DEGs related to energy metabolism, body defence, cell proliferation and biomineralization were found, where 1655 of them were downregulated, suggesting a strong negative response to the second heatwave. Our findings imply that some marine organisms can indeed tolerate heatwaves by boosting energy metabolism to support molecular defence, cell proliferation and biomineralization, but this capacity can be overwhelmed by repeated exposure to heatwaves. Since recurrence of heatwaves within a short period of time is predicted to be more prevalent in the future, the functioning of marine ecosystems would be disrupted if marine organisms fail to accommodate repeated extreme thermal stress.

Adaptive potential of the Mediterranean mussel Mytilus galloprovincialis to short-term environmental hypoxia

Environmental hypoxia naturally occurs in coastal ecosystems and bivalve mollusks have to frequently face fluctuations of dissolved oxygen concentrations. Exposure to hypoxia is often associated with the change of the antioxidant and functional status in bivalves, and restoration of the normal oxygen supply is considered to induce oxidative stress in tissues of mollusks. The study investigates changes in the activity of two antioxidant enzymes, catalase (CAT) and superoxide dismutase (SOD), as well as the expression level of SOD and CAT genes in gills of the Mediterranean mussel, Mytilus galloprovincialis, under exposure to low dissolved oxygen concentration (2.2 mg L^-1) for 24 h and 72 h, and 24 h reoxygenation period. We also evaluated the intracellular level of reactive oxygen species (ROS), mortality and changes in mitochondrial membrane potential in hemocytes following hypoxia-reoxygenation cycle. 24 h exposure to hypoxia significantly decreased activity of both enzymes, which then recovered up to control levels at the end of 72 h experimental period for SOD and after reoxygenation for CAT. Expression of antioxidant enzyme genes was up-regulated following the 72 h hypoxic exposure period and returned to the basal normoxic level after 24 h reoxygenation. Hypoxia demonstrated a time-dependent effect on the functional state of hemocytes. The 24 h exposure period did not influence aerobic respiration of hemocytes, but prolonged hypoxia (72 h) was associated with a substantial decrease in mitochondrial membrane potential of hemocytes. The intracellular ROS level and mortality of hemocytes did not change under hypoxia. Reoxygenation period was accompanied with a significant decrease of intracellular ROS level. This study indicated that hypoxia did not induce the pronounced oxidative stress in gills and the changes in the antioxidant status were reversible within 24 h of reoxygenation. Hemolymph demonstrated a stable functional state indicating the tolerance of mussels to short-time hypoxia.

Assessing the impact of atmospheric heatwaves on intertidal clams

Heatwaves have become more frequent and intense in the last two decades, resulting in detrimental effects on marine bivalves and ecosystems they sustain. Intertidal clams inhabit the most physiologically challenging habitats in coastal areas and live already near their thermal tolerance limits. However, whether and to what extent atmospheric heatwaves affect intertidal bivalves remain poorly understood. Here, we investigated physiological responses of the Manila clam, Ruditapes philippinarum, to heatwaves at air temperature regimes of 40 °C and 50 °C occurring frequently and occasionally at the present day in the Beibu Gulf, South China Sea. With the increasing intensity of heatwaves and following only two days of aerial exposure, Manila clams suffered 100 % mortality at 50 °C, indicating that they succumb to near future heatwaves, although they survived under various scenarios of moderate heatwaves. The latter is couched in energetic terms across levels of biological organization. Specifically, Manila clams acutely exposed to heatwaves enhanced their standard metabolic rate to fuel essential physiological maintenance, such as increasing activities of SOD, CAT, MDA, and AKP, and expression of HSP70. These strategies occur likely at the expense of fitness-related functions, as best exemplified by significant depressions in activities of enzymes (NKA, CMA, and T-ATP) and expression levels of genes (PT, KHK, CA, CAS, TYR, TNF-BP, and OSER). When heatwaves occurred again, Manila clams can respond and acclimate to thermal stress by implementing a suite of more ATP-efficient and less energy-costly compensatory mechanisms at various levels of biological organization. It is consequently becoming imperative to uncover underlying mechanisms responsible for such positive response and rapid acclimation to recurrent heatwaves.

Gene Co-Expression Network Analysis Reveals the Correlation Patterns Among Genes in Different Temperature Stress Adaptation of Manila Clam

The Manila clam (Ruditapes philippinarum) is one of the most important aquaculture species and widely distributed along the coasts of China, Japan, and Korea. Due to its wide distribution, it can tolerate a wide range of temperature. Studying the gene expression profiles of clam gills had found differentially expressed genes (DEGs) and pathway involved in temperature stress tolerance. A systematic study of cellular response to temperature stress may provide insights into the mechanism of acquired tolerance. Here, weighted gene co-expression network analysis (WGCNA) was carried out using RNA-seq data from gill transcriptome in response to high and low temperature stress. There are a total 32 gene modules, of which 18 gene modules were identified as temperature-related modules. Blue module was one significantly correlated with temperature which was associated with cellular metabolism, apoptosis pathway, ER stress, and others.

Transcriptome analysis reveals acclimation responses of pearl oysters to marine heatwaves

Marine heatwaves (MHWs) are weather-timescale extreme events in the oceans and can have devastating effects on marine bivalves and ecosystems they support, with considerable socio-economic consequences. Yet, the extent to which marine bivalves have the capacity to acclimate and adapt to MHWs remains unknown. Understanding molecular responses to MHWs is imperative to develop strategies for conservation of ecologically and economically important marine organisms. Here, using RNA-Seq, we investigate how various MHWs scenarios elicit molecular changes in threatened and vulnerable pearl oysters, Pinctada maxima (Jameson). Acute exposure of MHWs - mimicked by rapid increases of seawater temperature from 24 °C to 28 °C and 32 °C, respectively - significantly affected the expression levels of metabolic and immune-related genes, with thermal stress-responsive genes especially like HSP20, HSP70 and HSP90 being remarkably up-regulated. Following repeat exposure to MHWs, encouragingly, pearl oysters exhibited evident acclimation responses, as best exemplified by significantly lowered expression levels of key stress-responsive genes involved in metabolism and immunity in comparison to those observed during acute exposure. Findings of the present study provide a better understanding of molecular processes underpinning the acclimation and adaptation of marine bivalves to MHWs in the context of climate change.

Expression of p38MAPK and its regulation of apoptosis under high temperature stress in the razor clam Sinonovacula constricta

p38MAPK is a key branch of the MAPK (mitogen-activated protein kinase) pathway that plays an important role in physiological processes such as apoptosis, cell proliferation and growth. In this experiment, we screened and identified one p38MAPK gene in the razor clam Sinonovacula constricta, which encoded 359 amino acids and was widely expressed in various adult tissues. After 24 h of high temperature stress at 34 °C, the transcript expression of p38MAPK showed significant changes in all tested tissues. In particular in the gill and hepatopancreas tissues, where the expression increased 1.81 and 7.83 times compared with the control group, respectively (P < 0.01). Furthermore, we examined the expression of the apoptosis suppressor gene Bcl-2 and pro-apoptosis gene Bax by knock-down of p38MAPK with dsRNA interference in the gill and hepatopancreas tissues. The obvious up-regulation expression of Bcl-2 and significant suppression of Bax were observed, respectively (P < 0.01). Moreover, the TUNEL staining technique was used to detect apoptosis before and after interference. The degree of apoptosis in the gill and hepatopancreas tissues was reduced after interference with p38MAPK, and the ROS content was significantly reduced (P < 0.01). The results suggested that p38MAPK had a regulatory role in the heat tolerance of razor clams.

Metabolic response of Mercenaria mercenaria under heat and hypoxia stress by widely targeted metabolomic approach

In the context of global climatic changes, marine organisms have been exposed to environmental stressors including heat and hypoxia. This calls for the design of multi-stressors to uncover the impact of oceanic factors on aquatic organisms. So far, little is known about the metabolic response of marine organisms, especially bivalves, to the combined effects of heat and hypoxia. In this study, we employed widely targeted metabolomic analysis to study the metabolic response of gills in hard clam, a heat- and hypoxia-tolerant bivalve. A total of 810 metabolites were identified. Results showed that the heat group (HT) and heat plus hypoxia group (HL) had a higher number of differential metabolites than the hypoxia group (LO). Glycolysis was affected by the heat and heat plus hypoxia stress. Moreover, anaerobic metabolic biomarkers were accumulated marking the onset of anaerobic metabolism. Environmental stresses may affect Tricarboxylic acid (TCA) cycle. Accumulation of carnitine and glycerophospholipid may promote fatty acid β oxidation and maintain cell membrane stability, respectively. The high content of oxidized lipids (i.e., Leukotriene) in HL and HT groups implied that the organisms were under ROS stress. The significantly differential metabolites of organic osmolytes and vitamins might relieve ROS stress. Moreover, accumulation of thermoprotective osmolytes (monosaccharide, Trimethylamine N-oxide (TMAO)) accumulation was helpful to maintain protein homeostasis. This investigation provided new insights into the adaptation mechanisms of hard clam to heat, hypoxia and combined stress at the metabolite level and highlighted the roles of molecules and protectants.

Transcriptomic response of the intertidal limpet Patella vulgata to temperature extremes

Global warming is challenging wild species in land and water. In the intertidal zone, species are already living at their thermal limits, being vulnerable even to small increases in maximum habitat temperatures. Knowledge of the mechanisms by which many intertidal zone species cope with elevated temperatures is limited. We analysed the molecular thermal stress response of the limpet Patella vulgata under slight and frequent (one-day), and extreme and rare (three-day) warming events. Using RNA-seq to assess differential gene expression among treatments, differing molecular responses were obtained in the two treatments, with more changes in gene expression after the three-day event; with one-third of the differentially expressed transcripts being down-regulated. However, across treatments we observed shifts in gene expression for common aspects of the heat stress response including intra-cellular communication, protein chaperoning, proteolysis and cell cycle arrest. Of the 71,675 transcripts obtained, only 259 were differentially expressed after both heating events. From these, 218 defined the core group (i.e. genes induced by thermal stress with similar expression patterns irrespective of the magnitude of the warming event). The core group was composed of already well-studied genes in heat stress responses in intertidal organisms (e.g. heat shock proteins), but also genes from less explored metabolic pathways, e.g. the ubiquitin system, which were also fundamental regardless of the magnitude of the imposed warming. Moreover, we have also identified 41 signaling genes (i.e. a set of genes responding to both events and with expression patterns specific to the intensity of thermal stress), principally including genes involved in the maintenance of extracellular structure that have previously not been identified as part of the response to thermal stress in intertidal zone organisms. These signaling genes will be useful heat stress molecular biomarkers for monitoring heat stress in natural populations.

Gene Expression Profiles in Two Razor Clam Populations: Discerning Drivers of Population Status

With rapidly changing marine ecosystems, shifts in abundance and distribution are being documented for a variety of intertidal species. We examined two adjacent populations of Pacific razor clams (Siliqua patula) in lower Cook Inlet, Alaska. One population (east) supported a sport and personal use fishery, but this has been closed since 2015 due to declines in abundance, and the second population (west) continues to support commercial and sport fisheries. We used gene expression to investigate potential causes of the east side decline, comparing razor clam physiological responses between east and west Cook Inlet. The target gene profile used was developed for razor clam populations in Alaska based on physiological responses to environmental stressors. In this study, we identified no differences of gene expression between east and west populations, leading to two potential conclusions: (1) differences in factors capable of influencing physiology exist between the east and west and are sufficient to influence razor clam populations but are not detected by the genes in our panel, or (2) physiological processes do not account for the differences in abundance, and other factors such as predation or changes in habitat may be impacting the east Cook Inlet population.

Gradual effects of gradient concentrations of polystyrene nanoplastics on metabolic processes of the razor clams

With the widespread occurrence and accumulation of plastic waste in the world, plastic pollution has become a serious threat to ecosystem and ecological security, especially to estuarine and coastal areas. Understanding the impacts of changing nanoplastics concentrations on aquatic organisms living in these areas is essential for revealing the ecological effects caused by plastic pollution. In the present study, we revealed the effects of exposure to gradient concentrations (0.005, 0.05, 0.5 and 50 mg/L) of 75 nm polystyrene nanoplastics (PS-NPs) for 48 h on metabolic processes in muscle tissue of a bivalve, the razor clam Sinonovacula constricta, via metabolomic and transcriptomic analysis. Our results showed that PS-NPs caused dose-dependent adverse effects on energy reserves, membrane lipid metabolism, purine metabolism and lysosomal hydrolases. Exposure to PS-NPs reduced energy reserves, especially lipids. Membrane lipid metabolism was sensitive to PS-NPs with contents of phosphocholines (PC), phosphatidylethanolamines (PE) and phosphatidylserines (PS) increasing and degradation being inhibited in all concentrations. High concentrations of PS-NPs altered the purine metabolism via increasing contents of guanosine triphosphate (GTP) and adenine, which may be needed for DNA repair, and consuming inosine and hypoxanthine. During exposure to low concentrations of PS-NPs, lysosomal hydrolases in S. constricta, especially cathepsins, were inhibited while this influence was improved transitorily in 5 mg/L of PS-NPs. These adverse effects together impacted energy metabolism in S. constricta and disturbed energy homeostasis, which was manifested by the low levels of acetyl-CoA in high concentrations of PS-NPs. Overall, our results revealed the effects of acute exposure to gradient concentrations of PS-NPs on S. constricta, especially its metabolic process, and provide perspectives for understanding the toxicity of dynamic plastic pollution to coastal organisms and ecosystem.

Identification and Characterization of HSP90 Gene Family Reveals Involvement of HSP90, GRP94 and Not TRAP1 in Heat Stress Response in Chlamys farreri

Heat shock proteins 90 (HSP90s) are a class of ubiquitous, highly conserved, and multi-functional molecular chaperones present in all living organisms. They assist protein folding processes to form functional proteins. In the present study, three HSP90 genes, CfHSP90, CfGRP94 and CfTRAP1, were successfully identified in the genome of Chlamys farreri. The length of CfHSP90, CfGRP94 and CfTRAP1 were 7211 bp, 26,457 bp, and 28,699 bp, each containing an open reading frame (ORF) of 2181 bp, 2397 bp, and 2181 bp, and encoding proteins of 726, 798, and 726 amino acids, respectively. A transcriptomic database demonstrated that CfHSP90 and CfGRP94 were the primary functional executors with high expression during larval development and in adult tissues, while CfTRAP1 expression was low. Furthermore, all of the three CfHSP90s showed higher expression in gonads and ganglia as compared with other tissues, which indicated their probable involvement in gametogenesis and nerve signal transmission in C. farreri. In addition, under heat stress, the expressions of CfHSP90 and CfGRP94 were significantly up-regulated in the mantle, gill, and blood, but not in the heart. Nevertheless, the expression of CfTRAP1 did not change significantly in the four tested tissues. Taken together, in coping with heat stress, CfHSP90 and CfGRP94 could help correct protein folding or salvage damaged proteins for cell homeostasis in C. farreri. Collectively, a comprehensive analysis of CfHSP90s in C. farreri was conducted. The study indicates the functional diversity of CfHSP90s in growth, development, and environmental response, and our findings may have implications for the subsequent in-depth exploration of HSP90s in invertebrates.

Membrane lipid metabolism, heat shock response and energy costs mediate the interaction between acclimatization and heat-hardening response in the razor clam Sinonovacula constricta

Thermal plasticity on different time scales, including acclimation/acclimatization and heat-hardening response - a rapid adjustment for thermal tolerance after non-lethal thermal stress, can interact to improve the resilience of organisms to thermal stress. However, little is known about physiological mechanisms mediating this interaction. To investigate the underpinnings of heat-hardening responses after acclimatization in warm seasons, we measured thermal tolerance plasticity, and compared transcriptomic and metabolomic changes after heat hardening at 33 or 37°C followed by recovery of 3 or 24 h in an intertidal bivalve Sinonovacula constricta. Clams showed explicit heat-hardening responses after acclimatization in a warm season. The higher inducing temperature (37°C) caused less effective heat-hardening effects than the inducing temperature that was closer to the seasonal maximum temperature (33°C). Metabolomic analysis highlighted the elevated content of glycerophospholipids in all heat-hardened clams, which may help to maintain the structure and function of the membrane. Heat shock proteins (HSPs) tended to be upregulated after heat hardening at 37°C but not at 33°C, indicating that there was no complete dependency of heat-hardening effects on upregulated HSPs. Enhanced energy metabolism and decreased energy reserves were observed after heat hardening at 37°C, suggesting more energy costs during exposure to a higher inducing temperature, which may restrict heat-hardening effects. These results highlight the mediating role of membrane lipid metabolism, heat shock responses and energy costs in the interaction between heat-hardening response and seasonal acclimatization, and contribute to the mechanistic understanding of evolutionary change and thermal plasticity during global climate change.

Mitochondrial responses towards intermittent heat shocks in the eastern oyster, Crassostrea virginica

Frequent heat waves caused by climate change can give rise to physiological stress in many animals, particularly in sessile ectotherms such as bivalves. Most studies characterizing thermal stress in bivalves focus on evaluating the responses to a single stress event. This does not accurately reflect the reality faced by bivalves, which are often subject to intermittent heat waves. Here, we investigated the effect of intermittent heat stress on mitochondrial functions of the eastern oyster, Crassostrea virginica, which play a key role in setting the thermal tolerance of ectotherms. Specifically, we measured changes in mitochondrial oxygen consumption and H2O2 emission rates before, during and after intermittent 7.5°C heat shocks in oysters acclimated to 15 and 22.5°C. Our results showed that oxygen consumption was impaired following the first heat shock at both acclimation temperatures. After the second heat shock, results for oysters acclimated to 15°C indicated a return to normal. However, oysters acclimated to 22.5°C struggled more with the compounding effects of intermittent heat shocks as denoted by an increased contribution of FAD-linked substrates to mitochondrial respiration as well as high levels of H2O2 emission rates. However, both acclimated populations showed signs of potential recovery 10 days after the second heat shock, reflecting a surprising resilience to heat waves by C. virginica. Thus, this study highlights the important role of acclimation in the oyster's capacity to weather intermittent heat shock.

Identification, Characterization, and Expression Analysis Reveal Diverse Regulated Roles of Three MAPK Genes in Chlamys farreri Under Heat Stress

Mitogen-activated protein kinase (MAPK) cascades are fundamental signal transduction modules in all eukaryotic organisms, participating growth and development, as well as stress response. In the present study, three MAPK genes were successfully identified from the genome of Chlamys farreri, respectively, named CfERK1/2, CfJNK, and Cfp38, and only one copy of ERK, JNK, and p38 were detected. Domain analysis indicated that CfMAPKs possessed the typical domains, including S_TKc, Pkinase, and PKc_like domain. Phylogenetic analysis showed that three CfMAPKs of MAPK subfamilies exists in the common ancestor of vertebrates and invertebrates. All CfMAPKs specifically expressed during larval development and in adult tissues, and the expression level of CfERK1/2 and Cfp38 was apparently higher than that of CfJNK. Under heat stress, the expression of CfERK1/2 and Cfp38 were significantly downregulated and then upregulated in four tissues, while the expression of CfJNK increased in all tissues; these different expression patterns suggested a different molecular mechanism of CfMAPKs for bivalves to adapt to temperature changes. The diversity of CfMAPKs and their specific expression patterns provide valuable information for better understanding of the functions of MAPK cascades in bivalves.

Constructing a de novo transcriptome and a reference proteome for the bivalve Scrobicularia plana: Comparative analysis of different assembly strategies and proteomic analysis

Scrobicularia plana is a coastal and estuarine bivalve widely used in ecotoxicological studies. However, the underlying molecular mechanisms for S. plana pollutant responses are hardly known due to the lack of molecular databases. Thus, in this study we present a holistic approach to assess a robust reference transcriptome and proteome of this clam. A mixture of control and metal-exposed individuals was used for mRNA isolation. Four sets of high quality filtered preprocessed reads were generated (two quality scores and two sequenced lengths) and assembled with Mira, Ray and Trinity algorithms. The sixty-four generated assemblies were refined, filtered and evaluated for their proteomic quality. Eight assemblies presented top Detonate scores but one was selected due to its compactness and biological representation, which was generated: (i) from the highest quality dataset (Q20L100), (ii) using Trinity algorithm with all k-mers (AtKa), (iii) removing redundancy by CD-HIT (RR80), and (iv) filtering out poor contigs (F), that was subsequently named Q20L100AtKaRR80F. S. plana proteomic analysis revealed 10,017 peptide groups that corresponded to 2066 proteins with a wide coverage of molecular functions and biological processes, confirming the strength of the database generated.

Contrasting responsiveness of four ecologically and economically important bivalves to simulated heat waves

Heat waves are expected to increase in duration and frequency, impacting coastal ecosystems, especially intertidal organisms living near their thermal tolerance limits. Sedentary infaunal species are limited to some extent in escapes from sudden temperature changes, rather modifications to their physiology and behaviour are expected. This may lead to strong ecological and economic impacts on commercial bivalve species, such as Venerupis corrugata, Ruditapes decussatus, the introduced Ruditapes philippinarum and Cerastoderma edule, the most relevant in NW Spain. We investigated lethal and sublethal effects of heat during low tide on these species in the laboratory. Summer temperatures experienced within field, shallow sediments at approximately 2 cm depth i.e. 20 °C (control), 27 °C, 32 °C, and 37 °C, were replicated during four consecutive days and the diffusion of heat at the burrowing depth of each species was estimated; temperature exposure was expressed as degree hours above 22 °C. After two days of tidal exposure, C. edule and V. corrugata suffered significant mortalities, and also the most dramatic decrease in scope for growth (SFG) as well as reduction in burrowing activity. After four days under stress, all species had negative SFG. On recovery, species showed compensation at longer exposures, particularly C. edule. These effects of temperature on mortality, growth potential and burrowing ability may increase the time to achieve commercial size and exposure to predation. Particularly, V. corrugata, with a center of distribution lower in the intertidal and subtidal, and C. edule, shallower in the sediment, may be the most affected. Clearly the intensity and frequency of heat waves will affect these key species in the intertidal sediment flats changing ecosystem functioning and fisheries management strategies.

Living on the upper intertidal mudflat: Different behavioral and physiological responses to high temperature between two sympatric cerithidea snails with divergent habitat-use strategies

Temperature plays a major role in controlling species' distributions, and small-scale variation in the thermal environment are potentially an important factor that governs distributions on a local scale. For untangling the roles of behavioral and physiological adaptations on species' distribution at a small-scale level, we carried out a comparative study of two mudflat snails (genus Cerithidea) by determining these congeners' burying behavior, lethal temperature, cardiac performance and heat-shock protein (hsp70) gene expression. These two sympatric snails occupy different microhabitats on the upper intertidal mudflat. During periods of emersion, C. cingulata inhabits the open mudflat and C. largillierti usually aggregates around small rocks on the upper intertidal mudflat. Our results indicate that the two Cerithidea congeners show different behavioral and physiological responses to high temperature. Compared to C. largillierti, C. cingulata prefers to bury into the mud, has a higher thermal limit and a higher level of inducible expression of hsp70 mRNA, implying important roles of behavioral and physiological adaptations to the harsh thermal environment on the open mudflat. Furthermore, results of generalized additive modelling (GAM) analysis of cardiac performance and coefficient of variation (CV) of hsp70 mRNA expression showed high inter-individual variation in C. cingulata. These results highlight the importance of behavioral and physiological adaptions in sympatric species' distributions on the mudflat and help to shed light on the mechanisms of how small-scale differences in the thermal environment shape sympatric species' distributions.