An integrated, multi-level analysis of thermal effects on intertidal molluscs for understanding species distribution patterns

Pulsative venous return from the branchial vessels to the heart of the bivalve Mytilus galloprovincialis supports the constant-volume mechanism

In bivalves and gastropods, ventricle contraction causes a negative pressure in the auricles and increases venous return from the afferent oblique vein (AOV): the constant-volume (CV) mechanism. The flow in the AOV should be a pulsative flow synchronized with the ventricular contraction. The flow in the heart and adjacent vessels of Mytilus galloprovincialis were measured by magnetic resonance imaging to confirm this hypothesis. Under a regular heartbeat, pulsative flows in the AOV and branchial vessels (BVs) were almost completely synchronized with the flow in the aorta, while filling of the ventricle was in the opposite phase. Flows in the BVs were directed to the posterior direction, and a pair of BVs in the gill axes (the efferent BVs) were connected to the AOV. Based on the images of the whole pathway of the AOV in an oblique slice, we confirmed that haemolymph flow was evoked from the efferent BVs and flow into the ventricle via the auricle was completed in a single heartbeat. Therefore, the walls of the AOV and BVs could resist negative transmural pressure caused by the ventricular contraction. In conclusion, the auricle, the AOV and the BVs, including the gill filaments, act as a suction pump. The pulsative venous return is driven by the negative pressure of the AOV as in the CV mechanism, and the negative pressure in the efferent BVs could draw haemolymph from the sinus via the gill and the afferent BVs. Therefore, Mytilus can start and stop its heartbeat as necessary.

Environmental and social framework to protect marine bivalves under extreme weather events

Rising ocean temperatures, a consequence of anthropogenic climate change, are increasing the frequency, intensity, and magnitude of extreme marine heatwaves (MHWs). These persistent anomalous warming events can have severe ecological and socioeconomic impacts, threatening ecologically and economically vital organisms such as bivalves and the ecosystems they support. Developing robust environmental and social frameworks to enhance the resilience and adaptability of bivalve aquaculture is critical to ensuring the sustainability of this crucial food source. This review synthesizes the current understanding of the physiological and ecological impacts of MHWs on commercially important bivalve species farmed globally. We propose an integrated risk assessment framework that encompasses environmental monitoring, farm-level preparedness planning, and community-level social support systems to safeguard bivalve aquaculture. Specifically, we examine heatwave prediction models, local mitigation strategies, and social programs that could mitigate the impacts on bivalve farms and vulnerable coastal communities economically dependent on this fishery. At the farm level, adaptation strategies such as selective breeding for heat-tolerant strains, optimized site selection, and adjustments to culture practices can improve survival outcomes during MHWs. Robust disease surveillance and management programs are essential for early detection and rapid response. Furthermore, we highlight the importance of stakeholder engagement, knowledge exchange, and collaborative governance in developing context-specific, inclusive, and equitable safeguard systems. Proactive measures, such as advanced forecasting tools like the California Current Marine Heat Wave Tracker developed by NOAA's Southwest Fisheries Science Center, enable preemptive action before losses occur. Coordinated preparation and response, underpinned by continuous monitoring and adaptive management, promise to protect these climate-vulnerable food systems and coastal communities. However, sustained research, innovation, and cross-sector collaboration are imperative to navigate the challenges posed by our rapidly changing oceans.

The role of species thermal plasticity for alien species invasibility in a changing climate: A case study of Lophocladia trichoclados

The Mediterranean Sea provides fertile ground for understanding the complex interplay between invasive species and native habitats, particularly within the context of climate change. This thermal tolerance study reveals the remarkable ability of Lophocladia trichoclados, a red algae species that has proven highly invasive, to adapt to varying temperatures, particularly thriving in colder Mediterranean waters, where it can withstand temperatures as low as 14 °C, a trait not observed in its native habitat. This rapid acclimation, occurring in less than a century, might entail a trade-off with high temperature resistance. Additionally, all sampled populations in the Mediterranean share the same haplotype, suggesting a common origin and the possibility that we might be facing an exceptionally acclimatable and invasive strain. This high degree of acclimatability could determine the future spread capacity in a changing scenario, highlighting the importance of considering both acclimation and adaptation in understanding the expansion of invasive species' ranges.

Flight or fight: different strategies of intertidal periwinkle Littoraria sinensis coping with high temperature across populations

Intertidal organisms usually live near their upper thermal limits, and are vulnerable to future global warming. As a vital response to thermal stress, thermoregulatory strategy in physiological and behavioral performance is essential for organisms coping with thermal stress and surviving the changing world. To investigate the relationship between the thermoregulatory strategy and habitat temperature, in the present study, we comparatively investigated the thermal responsive strategy among different geographic populations of the supralittoral snail Littoraria sinensis by determining snails' cardiac function and behavioral performance. Our results indicated that populations inhabiting high ambient temperatures had higher sublethal temperatures (i.e. Arrhenius breakpoint temperatures, ABTs, the temperature at which the heart rate shapely decreases with further heating) and lethal temperatures (i.e. Flatline temperatures, FLTs, the temperature at which heart rate ceases), and behaved less actively (e.g. shorter moving distances and shorter moving time) in the face of high and rising temperatures-a physiological fight strategy. On the other hand, populations at relatively low ambient temperatures had relatively lower physiological upper thermal limits with lower ABTs and FLTs and moved more actively in the face of high and rising temperatures-a behavioral flight strategy. These results demonstrate that the thermoregulatory strategies of the snails are closely related to their habitat temperatures and are different among populations surviving divergent thermal environments.

Genome-wide sequencing reveals geographical variations in the thermal adaptation of an aquaculture species with frequent seedling introductions

Climate change and human activity are essential factors affecting marine biodiversity and aquaculture, and understanding the impacts of human activities on the genetic structure to increasing high temperatures is crucial for sustainable aquaculture and marine biodiversity conservation. As a commercially important bivalve, the Manila clam Ruditapes philippinarum is widely distributed along the coast of China, and it has been frequently introduced from Fujian Province, China, to other regions for aquaculture. In this study, we collected four populations of Manila clams from different areas to evaluate their thermal tolerance by measuring cardiac performance and genetic variations using whole-genome resequencing. The upper thermal limits of the clams showed high variations within and among populations. Different populations displayed divergent genetic compositions, and the admixed population was partly derived from the Zhangzhou population in Fujian Province, implying a complex genomic landscape under the influence of local genetic sources and human introductions. Multiple single nucleotide polymorphisms (SNPs) were associated with the cardiac functional traits, and some of these SNPs can affect the codon usage and the structural stability of the resulting protein. This study shed light on the importance of establishing long-term ecological and genetic monitoring programs at the local level to enhance resilience to future climate change.

Multiple genetic sources facilitate the northward range expansion of an intertidal oyster along China's coast

Coastal artificial structures on the former mudflats provide available habitats for the rocky intertidal species which can establish new populations in these emerging habitats over their former distribution range limits. As a former southern species, the oyster Crassostrea sikamea has become a pioneer and rapidly invaded the artificial shorelines in northern China. We used a seascape genomics approach to investigate the population structure and genetic sources of C. sikamea on the coastal artificial structures, which is crucial for understanding the genetic mechanisms driving species distribution range expansion and invasion pathway of intertidal species. Five C. sikamea populations, including two artificial substrate populations (WGZ and ZAP), one oyster reef population (LS), and two natural rocky shore populations (ZS and XM), were measured using single nucleotide polymorphism (SNPs) obtained from double digest restriction-site associated DNA sequencing (ddRAD-Seq). Redundancy analyses (RDA) were implemented for investigating the relationship between local temperature variables and the temperature adaptability of C. sikamea. Genetic diversity, direction and strength of gene flow, and population structure all revealed that the LS and ZS populations were the genetic sources for the oyster populations on the emerging northern coastal artificial structures. Results of RDA showed that there were different adaptive potentials for northern and southern populations to local temperature variables and the oyster reef population which frequently suffers from heat stress owned high heat adaptability. The ZS population as a genetic source nearby the Yangtze River estuary provided mass larvae for the northern populations, and the other genetic source, the heat-tolerant LS population, in the oyster reef played an important role in the post-settlement success by providing preadapted genotypes. These results highlight the importance of multiple sources with divergent adaptative capabilities for biological invasion, and also emphasize the importance of the oyster reef in coastal biodiversity and conservation.

Not just range limits: Warming rate and thermal sensitivity shape climate change vulnerability in a species range center

Climate change manifests unevenly across space and time and produces complex patterns of stress for ecological systems. Species can also show substantial among-population variability in response to environmental change across their geographic range due to evolutionary processes. Explanatory factors or their proxies, such as temperature and latitude, help parse these sources of environmental and intraspecific variability; however, overemphasizing latitudinal trends can obscure the role of local environmental conditions in shaping population vulnerability to climate change. Focusing on the geographic center of a species range to disentangle latitude, we test the hypothesis that populations from warmer regions of a species range are more vulnerable to ocean warming. We conducted a mesocosm experiment and field reciprocal transplant with four populations of a marine snail, Nucella lamellosa, from two regions in British Columbia, Canada, that differ in thermal characteristics: the Central Coast, a cool region, and the Strait of Georgia, one of the warmest regions of this species' range and one that is warming faster than the Central Coast. Populations from the Strait of Georgia experienced growth reductions at contemporary summertime seawater temperatures in the laboratory and showed stark reductions in survival and growth under future seawater conditions and when outplanted at their native transplant sites. This indicates a high vulnerability to ocean warming, especially given the faster rate of ocean warming in this region. In contrast, populations from the cooler Central Coast demonstrated high performance at contemporary seawater temperatures and high growth and survival in projected future seawater temperatures and at their native outplant sites. Given their position within the geographic center of N. lamellosa's range, extirpation events in the vulnerable Strait of Georgia populations could compromise connectivity within the metapopulation and lead to gaps across this species' range. Overall, our study supports predictions that populations from warm regions of species ranges are more vulnerable to environmental warming, suggests that the Strait of Georgia and other inland or coastal seas could be focal points for climate change effects and ecological transformation, and emphasizes the importance of analyzing climate change vulnerability in the context of regional environmental data and throughout a species' range.

Multi-omics reveals largely distinct transcript- and protein-level responses to the environment in an intertidal mussel

Organismal responses to stressful environments are influenced by numerous transcript- and protein-level mechanisms, and the relationships between expression changes at these levels are not always straightforward. Here, we used paired transcriptomic and proteomic datasets from two previous studies from gill of the California mussel, Mytilus californianus, to explore how simultaneous transcript and protein abundance patterns may diverge under different environmental scenarios. Field-acclimatized mussels were sampled from two disparate intertidal sites; individuals from one site were subjected to three further treatments (common garden, low-intertidal or high-intertidal outplant) that vary in temperature and feeding time. Assessing 1519 genes shared between the two datasets revealed that both transcript and protein expression patterns differentiated the treatments at a global level, despite numerous underlying discrepancies. There were far more instances of differential expression between treatments in transcript only (1451) or protein only (226) than of the two levels shifting expression concordantly (68 instances). Upregulated expression of cilium-associated transcripts (likely related to feeding) was associated with relatively benign field treatments. In the most stressful treatment, transcripts, but not proteins, for several molecular chaperones (including heat shock proteins and endoplasmic reticulum chaperones) were more abundant, consistent with a threshold model for induction of translation of constitutively available mRNAs. Overall, these results suggest that the relative importance of transcript- and protein-level regulation (translation and/or turnover) differs among cellular functions and across specific microhabitats or environmental contexts. Furthermore, the degree of concordance between transcript and protein expression can vary across benign versus acutely stressful environmental conditions.

Roles of multi-level temperature-adaptive responses and microhabitat variation in establishing distributions of intertidal species

How intertidal species survive their harsh environment and how best to evaluate and forecast range shifts in species distribution are two important and closely related questions for intertidal ecologists and global change biologists. Adaptive variation in responses of organisms to environmental change across all levels of biological organization - from behavior to molecular systems - is of key importance in setting distribution patterns, yet studies often neglect the interactions of diverse types of biological variation (e.g. differences in thermal optima owing to genetic and acclimation-induced effects) with environmental variation, notably at the scale of microhabitats. Intertidal species have to cope with extreme and frequently changing thermal stress, and have shown high variation in thermal sensitivities and adaptive responses at different levels of biological organization. Here, I review the physiological and biochemical adaptations of intertidal species to environmental temperature on multiple spatial and temporal scales. With fine-scale datasets for the thermal limits of individuals and for environmental temperature variation at the microhabitat scale, we can map the thermal sensitivity for each individual in different microhabitats, and then scale up the thermal sensitivity analysis to the population level and, finally, to the species level by incorporating physiological traits into species distribution models. These more refined mechanistic models that include consideration of physiological variations have higher predictive power than models that neglect these variations, and they will be crucial to answering the questions posed above concerning adaptive mechanisms and the roles they play in governing distribution patterns in a rapidly changing world.

Evaluation of species thermal sensitivity with individual-based physiological performance

Ignoring intraspecific variations can prevent us from accurately assessing species' thermal sensitivity to global warming. Individual-based physiological performance provides a feasible solution to depict species' thermal sensitivity using a bottom-up approach. We measured the cardiac performance of intertidal bivalves (1159 individuals from multiple populations of six bivalves), determined the upper thermal limit of each individual, calculated the proportions of individuals suffering sublethal/lethal heat stress, and mapped sensitive regions to high temperatures. Results showed that high inter-individual variations of physiological performance existed in levels of populations and species, and species' thermal sensitivity was positively related to the intraspecific variations of heat tolerance. This bottom-up approach scaled up from individual, population to species emphasizes the importance of individual-based physiology performance in assessing thermal sensitivity across different hierarchical levels and enables better evaluating and forecasting of species responses to global warming.

Exploring the adaptability of the secondary structure of mRNA to temperature in intertidal snails based on SHAPE experiments

RNA-based thermal regulation is an important strategy for organisms to cope with temperature changes. Inhabiting the intertidal rocky shore, a key interface of the ocean, atmosphere and terrestrial environments, intertidal species have developed variable thermal adaptation mechanisms; however, adaptions at the RNA level remain largely uninvestigated. To examine the relationship between mRNA structural stability and species distribution, in the present study, the secondary structure of cytosolic malate dehydrogenase (cMDH) mRNA of Echinolittorina malaccana, Echinolittorina radiata and Littorina brevicula was determined using selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE), and the change in folding free energy of formation (ΔGfold) was calculated. The results showed that ΔGfold increased as the temperature increased. The difference in ΔGfold (ΔΔGfold) between two specific temperatures (25 versus 0°C, 37 versus 0°C and 57 versus 0°C) differed among the three species, and the ΔΔGfold value of E. malaccana was significantly lower than those of E. radiata and L. brevicula. The number of stems of cMDH mRNA of the snails decreased with increasing temperature, and the breakpoint temperature of E. malaccana was the highest among these. The number of loops was also reduced with increasing temperature, while the length of the loop structure increased accordingly. Consequently, these structural changes can potentially affect the translational efficiency of mRNA. These results imply that there were interspecific differences in the thermal stability of RNA secondary structures in intertidal snails, and these differences may be related to snail distribution.

Contribution of HIF-1α to Heat Shock Response by Transcriptional Regulation of HSF1/HSP70 Signaling Pathway in Pacific Oyster, Crassostrea gigas

Ocean temperature rising drastically threatens the adaptation and survival of marine organisms, causing serious ecological impacts and economic losses. It is crucial to understand the adaptive mechanisms of marine organisms in response to high temperature. In this study, a novel regulatory mechanism that is mediated by hypoxia-inducible factor-1α (HIF-1α) was revealed in Pacific oyster (Crassostrea gigas) in response to heat stress. We identified a total of six HIF-1α genes in the C. gigas genome, of which HIF-1α and HIF-1α-like5 were highly induced under heat stress. We found that the HIF-1α and HIF-1α-like5 genes played critical roles in the heat shock response (HSR) through upregulating the expression of heat shock protein (HSP). Knocking down of HIF-1α via RNA interference (RNAi) inhibited the expression of heat shock factor 1 (HSF1) and HSP70 genes in C. gigas under heat stress. Both HIF-1α and HIF-1α-like5 promoted the transcriptional activity of HSF1 by binding to hypoxia response elements (HREs) within the promoter region. Furthermore, the survival of C. gigas under heat stress was significantly decreased after knocking down of HIF-1α. This work for the first time revealed the involvement of HIF-1α/HSF1/HSP70 pathway in response to heat stress in the oyster and provided an insight into adaptive mechanism of bivalves in the face of ocean warming.

Surviving hot summer: Roles of phenotypic plasticity of intertidal mobile species considering microhabitat environmental heterogeneity

For species inhabiting warming and variable thermal environment, coordinated changes in heat tolerance to temperature fluctuations, which largely depend on phenotypic plasticity, are pivotal in buffering high temperatures. Determining the roles of phenotypic plasticity in wild populations and common garden experiments help us understand how organisms survive hot summer and the warming world. We thus monitored the operative temperature of the intertidal limpets Cellana toreuma in both emergent rock and tidal pool microhabitats from June to October 2021, determined the variations of upper thermal limits of short-term acclimated and long-term acclimated limpets from different microhabitats (emergent rock and tidal pool), and further calculated the relationship between the upper thermal limits and acclimation capacity. Our results indicated that living on the emergent rock, limpets encountered more extreme events in summer. For the short-term acclimated samples, limpets on the emergent rock exhibited obvious variations of sublethal thermal limit (i.e., Arrhenius Break Point of cardiac performance, ABT) during summer months, however, this variation of ABT was absent in the limpets in the tidal pool. After the laboratory long-term acclimation, the ABTs and FLTs (Flat Line Temperature of cardiac performance, as an indicator of lethal temperature) of limpets both on the rock and in the tidal pool increased significantly in October, implying the potential existence of selection during the hot summer. Our results further showed that environmental temperature was an important driver of phenotypic plasticity. This study highlighted the changes in the thermal tolerance of intertidal limpets during summer in different microhabitats.

Intertidal limits shape covariation between metabolic plasticity, oxidative stress and telomere dynamics in Pacific oyster (Crassostrea gigas)

In intertidal zones, species such as sessile shellfish exhibit extended phenotypic plasticity to face rapid environmental changes, but whether frequent exposure to intertidal limits of the distribution range impose physiological costs for the animal remains elusive. Here, we explored how phenotypic plasticity varied along foreshore range at multiple organization levels, from molecular to cellular and whole organism acclimatization, in the Pacific oyster (Crassostrea gigas). We exposed 7-month-old individuals for up to 16 months to three foreshore levels covering the vertical range for this species, representing 20, 50 and 80% of the time spent submerged monthly. Individuals at the upper range limit produced energy more efficiently, as seen by steeper metabolic reactive norms and unaltered ATP levels despite reduced mitochondrial density. By spending most of their time emerged, oysters mounted an antioxidant shielding concomitant with lower levels of pro-oxidant proteins and postponed age-related telomere attrition. Instead, individuals exposed at the lower limit range near subtidal conditions showed lower energy efficiencies, greater oxidative stress and shorter telomere length. These results unraveled the extended acclimatization strategies and the physiological costs of living too fast in subtidal conditions for an intertidal species.

Genetic Diversity, Population Structure, and Environmental Adaptation Signatures of Chinese Coastal Hard-Shell Mussel Mytilus coruscus Revealed by Whole-Genome Sequencing

The hard-shell mussel (Mytilus coruscus) is widespread in the temperate coastal areas of the northwest Pacific and holds a significant position in the shellfish aquaculture market in China. However, the natural resources of this species have been declining, and population genetic studies of M. coruscus are also lacking. In this study, we conducted whole-genome resequencing (WGR) of M. coruscus from eight different latitudes along the Chinese coast and identified a total of 25,859,986 single nucleotide polymorphism (SNP) markers. Our findings indicated that the genetic diversity of M. coruscus from the Zhoushan region was lower compared with populations from other regions. Furthermore, we observed that the evolutionary tree clustered into two primary branches, and the Zhangzhou (ZZ) population was in a separate branch. The ZZ population was partly isolated from populations in other regions, but the distribution of branches was not geographically homogeneous, and a nested pattern emerged, consistent with the population differentiation index (FST) results. To investigate the selection characteristics, we utilized the northern M. coruscus populations (Dalian and Qingdao) and the central populations (Zhoushan and Xiangshan) as reference populations and the southern ZZ population as the target population. Our selection scan analysis identified several genes associated with thermal responses, including Hsp70 and CYP450. These genes may play important roles in the adaptation of M. coruscus to different living environments. Overall, our study provides a comprehensive understanding of the genomic diversity of coastal M. coruscus in China and is a valuable resource for future studies on genetic breeding and the evolutionary adaptation of this species.

The environmental cellular stress response: the intertidal as a multistressor model

The wild poses a multifaceted challenge to the maintenance of cellular function. Therefore, a multistressor approach is essential to predict the cellular mechanisms which promote homeostasis and underpin whole-organism tolerance. The intertidal zone is particularly dynamic, and thus, its inhabitants provide excellent models to assess mechanisms underpinning multistressor tolerance. Here, we critically review our current understanding of the regulation of the cellular stress response (CSR) under multiple abiotic stressors in intertidal organisms and consider to what extent a multistressor approach brings us closer to understanding responses in the wild. The function of the CSR has been well documented in laboratory and field exposures with a view to understanding single-stressor thermal effects. Multistressor studies still remain relatively limited in comparison but have applied three main approaches: (i) laboratory application of multiple stressors in isolation, (ii) multiple stressors applied in combination, and (iii) field-based correlation of multiple stressors against the CSR. The application of multiple stressors in isolation has allowed the identification of putative, shared stress pathways but overlooks non-additive stressor interactions on the CSR. Combined stressor studies are relatively limited in number but already highlight variable effects on the CSR dependent upon stressor type, timing, and magnitude. Field studies have allowed the identification of responsive components of the CSR to various stressors in situ but are correlative, not causative. A combined approach involving laboratory multistressor studies linking the CSR to whole-organism tolerance as well as field studies is required if we are to understand the role of the CSR in the natural environment.

Coping with harsh heat environments: molecular adaptation of metabolic depression in the intertidal snail Echinolittorina radiata

Harsh thermal environments in the rocky intertidal zone pose serious physiological and molecular challenges to the inhabitants. Metabolic depression is regarded as an energy-conserving feature of intertidal species. To understand the molecular mechanism of metabolic depression, we investigated physiological and transcriptomic responses in the intertidal snail Echinolittorina radiata. The metabolic rate and expression of most genes were insensitive to temperatures ranging from 33 to 45 °C and then increased with further heating to 52 °C. Different from other genes, the genes involved in heat shock response (HSR) and oxidative stress response (OSR) (e.g., genes encoding heat shock protein 70 (HSP70) and cytochrome P450 protein (CYP450)) kept upregulating during metabolic depression. These high levels of HSR and OSR genes should be important for surviving the harsh thermal environments on the rocky shore. In the population experiencing more frequent moderate heat events, the depression breadth was larger, and the change in magnitude of upregulation was insensitive for HSR genes (e.g., HSP70s) but heat-sensitive for OSR genes (e.g., CYP450s) at the temperature of 37 to 45 °C. These findings indicate that both the thermal sensitivity of HSR and OSR genes and the insensitivity of metabolic genes are crucial for surviving extreme intertidal environments, and different populations of the same species rely on various physiological mechanisms to differing extents to deal with heat stress. The cellular stress response is not a "one size fits all" response across populations largely depending on local thermal regimes.

Effects of heating rate and shell colour on the cardiac thermal performance in a polymorphic gastropod Batillaria attramentaria

Heating rate has gained extensive attention in mechanistic understanding of physiological responses to changing thermal conditions in the context of climate change. In polymorphic gastropods, differences in the absorption of solar energy between dark- and light-coloured individuals lead to supposable differences in their heating rates and body temperatures in sunshine. In the present study, we examined the effect of heating rate on heart rate (HR) in a polymorphic gastropod Batillaria attramentaria. By using biomimetic models, we found that daily maximum temperature of snails with a dark unbanded shell (D-type morph) was higher than snails with a white line on the upper side of each whorl (UL-type morph) by 0.6 °C when exposed to sunlight, but there was no apparent difference in heating rates between D- and UL-type models. We measured HR of snails at various heating rates from 3.0 to 9.0 °C h^-1. Faster heating rates significantly increased maximum thermal tolerance in both D- and UL-type snails, highlighting the importance to have thorough knowledge on the heating rate in the field to obtain accurate maximum thermal limit of gastropods. Critical temperature at which HR precipitously declines was higher in D-type snails than UL-type snails. Our results suggested that the impacts of heating rate as well as the shell colour should be considered to gain a mechanistic understanding of the population dynamics of polymorphic gastropods.

Temperature variability and metabolic adaptation in terrestrial and aquatic ectotherms

Thermodynamics is a major factor determining rates of energy expenditure, rates of biochemical dynamics, and ultimately the biological and ecological processes linked with resilience to global warming in ectothermic organisms. Nonetheless, whether ectothermic organisms exhibit general adaptive metabolic responses to cope with worldwide variation in thermal conditions has remained as an open question. Here we combine a model comparison approach with a global dataset of standard metabolic rates (SMR), including 1,160 measurements across 788 species of aquatic invertebrates, insects, fishes, amphibians and reptiles, to investigate the association between metabolic rates and environmental temperatures in their respective habitats. Our analyses suggest that variation in SMR after removing allometric and thermodynamic effects is best explained by the temperature range encountered across seasons, which always provided a better fit than the average temperature for the hottest and coldest month and mean annual temperatures. This pattern was consistent across taxonomic groups and robust to sensitivity analyses. Nonetheless, aquatic and terrestrial lineages responded differently to seasonality, with SMR declining - 6.8% °C^-1 of thermal range across seasons in aquatic organisms and increasing 2.8% °C^-1 in terrestrial organisms. These responses may reflect alternative strategies to mitigate the impact of increments in warmer temperatures on energy expenditure, either by means of metabolic reduction in thermally homogeneous water bodies or effective behavioral thermoregulation to exploit temperature heterogeneity on land.

Proteostasis in ice: the role of heat shock proteins and ubiquitin in the freeze tolerance of the intertidal mussel, Mytilus trossulus

The bay mussel, Mytilus trossulus, is an animal that can survive extracellular ice formation. Depending on air and ocean temperatures, freeze tolerant intertidal organisms, like M. trossulus, may freeze and thaw many times during the winter. Freezing can cause protein denaturation, leading to an induction of the heat shock response with expression of chaperone proteins like the 70 kDa heat shock protein (HSP70), and an increase in ubiquitin-conjugated proteins. There has been little work on the mechanisms of freeze tolerance in intertidal species, limiting our understanding of this survival strategy. Additionally, this limited research has focused solely on the effects of single freezing events, but the act of repeatedly crossing the freezing threshold may present novel physiological or biochemical stressors that have yet to be discovered. Mytilus are important ecosystem engineers and provide habitat for other intertidal species, thus understanding their physiology under thermal extremes is important for preserving shoreline health. We predicted that repeated freeze exposures would increase mortality, upregulate HSP70 expression, and increase ubiquitin conjugates in mussels, relative to single, prolonged freeze exposures. Mytilus trossulus from Vancouver, Canada were repeatedly frozen for a combination of 1 × 8 h, 2 × 4 h, or 4 × 2 h. We then compared mortality, HSP70 expression, and the quantity of ubiquitin-conjugated proteins across experimental groups. We found a single 8-h freeze caused significantly more mortality than repeated freeze-thaw cycles. We also found that HSP70 and ubiquitinated protein was upregulated exclusively after freeze-thaw cycles, suggesting that freeze-thaw cycles offer a period of damage repair between freezes. This indicates that freeze-thaw cycles, which happen naturally in the intertidal, are crucial for M. trossulus survival in sub-zero temperatures.

Variations of Supercooling Capacity in Intertidal Gastropods

Winter low-temperature confines species distribution. Intertidal gastropods are distributed from tropical to polar zones, facing variable intensities and durations of low temperatures. They usually set their supercooling points (SCPs) at high subzero temperatures to stimulate freezing. However, the variations in SCP in intertidal gastropods at intraspecific and interspecific levels remain poorly understood. Here, we measured the body size, cooling rate, and SCP of nine intertidal gastropod species in China. These species were distributed in high or middle intertidal zone with different geographic distributions. The average SCPs (-4.27~-7.10 °C) and the coefficients of variation of SCP (22.6%~45.9%) were high in all species. At the intraspecific level, the supercooling capacity was positively correlated with the cooling rate. Interspecifically, the supercooling capacity was closely related to the cooling rate, and also to the species' geographical distribution. Northern high-shore species showed lower SCPs, while southern high-shore species had higher SCPs. There was no difference in SCP between widespread high- and mid-shore species. Our results indicated that the supercooling capability is potentially an adaptative response to the local winter temperatures, and the cooling rate is a factor in determining the difference in SCP at the intraspecific and interspecific levels.

Microhabitat-specific diurnal metabolomic responses of the intertidal limpet Cellana toreuma to winter low temperature

High-throughput determination of circadian rhythms in metabolic response and their divergent patterns in various microhabitats are crucial for understanding how organisms respond to environmental stresses. A mid-intertidal limpet Cellana toreuma was collected at various time points across both daytime and nighttime in winter during low tide for investigating the diurnal metabolomic responses to cold stress and elucidating the divergent metabolic responses to temperature variations across microhabitats. Temperatures of emergent rock microhabitats were lower than the tidal pool and even aggravated at night. A series of metabolomic responses exhibited coordinated diurnal changes in winter. Metabolic responses which were associated with cellular stress responses and energy metabolism of emergent rock microhabitat individuals were highly induced compared to the tidal pool ones. This study shed light on the diurnal patterns of metabolomic responses of intertidal molluscs in the field and emphasized the variations in metabolic responses between microhabitats.

Genome-wide sequencing identifies a thermal-tolerance related synonymous mutation in the mussel, Mytilisepta virgata

The roles of synonymous mutations for adapting to stressful thermal environments are of fundamental biological and ecological interests but poorly understood. To study whether synonymous mutations influence thermal adaptation at specific microhabitats, a genome-wide genotype-phenotype association analysis is carried out in the black mussels Mytilisepta virgata. A synonymous mutation of Ubiquitin-specific Peptidase 15 (MvUSP15) is significantly associated with the physiological upper thermal limit. The individuals carrying GG genotype (the G-type) at the mutant locus possess significantly lower heat tolerance compared to the individuals carrying GA and AA genotypes (the A-type). When heated to sublethal temperature, the G-type exhibit higher inter-individual variations in MvUSP15 expression, especially for the mussels on the sun-exposed microhabitats. Taken together, a synonymous mutation in MvUSP15 can affect the gene expression profile and interact with microhabitat heterogeneity to influence thermal resistance. This integrative study sheds light on the ecological importance of adaptive synonymous mutations as an underappreciated genetic buffer against heat stress and emphasizes the importance of integrative studies at a microhabitat scale for evaluating and predicting the impacts of climate change.

The effects of temperature and pH on the reproductive ecology of sand dollars and sea urchins: Impacts on sperm swimming and fertilization

In an era of climate change, impacts on the marine environment include warming and ocean acidification. These effects can be amplified in shallow coastal regions where conditions often fluctuate widely. This type of environmental variation is potentially important for many nearshore species that are broadcast spawners, releasing eggs and sperm into the water column for fertilization. We conducted two experiments to investigate: 1) the impact of water temperature on sperm swimming characteristics and fertilization rate in sand dollars (Dendraster excentricus; temperatures 8-38°C) and sea urchins (Mesocentrotus franciscanus; temperatures 8-28°C) and; 2) the combined effects of multiple stressors (water temperature and pH) on these traits in sand dollars. We quantify thermal performance curves showing that sand dollar fertilization rates, sperm swimming velocities, and sperm motility display remarkably wide thermal breadths relative to red urchins, perhaps reflecting the wider range of water temperatures experienced by sand dollars at our field sites. For sand dollars, both temperature (8, 16, 24°C) and pH (7.1, 7.5, 7.9) affected fertilization but only temperature influenced sperm swimming velocity and motility. Although sperm velocities and fertilization were positively correlated, our fertilization kinetics model dramatically overestimated measured rates and this discrepancy was most pronounced under extreme temperature and pH conditions. Our results suggest that environmental stressors like temperature and pH likely impair aspects of the reproductive process beyond simple sperm swimming behavior.

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.

Solutions: how adaptive changes in cellular fluids enable marine life to cope with abiotic stressors

The seas confront organisms with a suite of abiotic stressors that pose challenges for physiological activity. Variations in temperature, hydrostatic pressure, and salinity have potential to disrupt structures, and functions of all molecular systems on which life depends. During evolution, sequences of nucleic acids and proteins are adaptively modified to "fit" these macromolecules for function under the particular abiotic conditions of the habitat. Complementing these macromolecular adaptations are alterations in compositions of solutions that bathe macromolecules and affect stabilities of their higher order structures. A primary result of these "micromolecular" adaptations is preservation of optimal balances between conformational rigidity and flexibility of macromolecules. Micromolecular adaptations involve several families of organic osmolytes, with varying effects on macromolecular stability. A given type of osmolyte generally has similar effects on DNA, RNA, proteins and membranes; thus, adaptive regulation of cellular osmolyte pools has a global effect on macromolecules. These effects are mediated largely through influences of osmolytes and macromolecules on water structure and activity. Acclimatory micromolecular responses are often critical in enabling organisms to cope with environmental changes during their lifetimes, for example, during vertical migration in the water column. A species' breadth of environmental tolerance may depend on how effectively it can vary the osmolyte composition of its cellular fluids in the face of stress. Micromolecular adaptations remain an under-appreciated aspect of evolution and acclimatization. Further study can lead to a better understanding of determinants of environmental tolerance ranges and to biotechnological advances in designing improved stabilizers for biological materials.

Demographic, physiological, and genetic factors linked to the poleward range expansion of the snail Nerita yoldii along the shoreline of China

Species range shift is one of the most significant consequences of climate change in the Anthropocene. A comprehensive study, including demographic, physiological, and genetic factors linked to poleward range expansion, is crucial for understanding how the expanding population occupies the new habitat. In the present study, we investigated the demographic, physiological, and genetic features of the intertidal gastropod Nerita yoldii, which has extended its northern limit by ~200 km over the former biogeographic break of the Yangtze River Estuary during recent decades. The neutral SNPs data showed that the new marginal populations formed a distinct cluster established by a few founders. Demographic modelling analysis revealed that the new marginal populations experienced a strong genetic bottleneck followed by recent demographic expansion. Successful expansion that overcame the founder effect might be attributed to its high capacity of rapid population growth and multiple introductions. According to the non-neutral SNPs under diversifying selection, there were high levels of heterozygosity in the new marginal populations, which might be beneficial for adapting to the novel thermal conditions. The common garden experiment showed that the new marginal populations have evolved divergent transcriptomic and physiological responses to heat stress, allowing them to occupy and survive in the novel environment. Lower transcriptional plasticity was observed in the new marginal populations. These results suggest a new biogeographic pattern of N. yoldii has formed with the occurrence of demographic, physiologic, and genetic changes, and emphasize the roles of adaptation of marginal populations during range expansion.

Effects of heat acclimation on cardiac function in the intertidal mussel Mytilus californianus: can laboratory-based indices predict survival in the field?

Thermal performance curves are commonly used to investigate the effects of heat acclimation on thermal tolerance and physiological performance. However, recent work indicates that the metrics of these curves heavily depend on experimental design and may be poor predictors of animal survival during heat events in the field. In intertidal mussels, cardiac thermal performance (CTP) tests have been widely used as indicators of animals' acclimation or acclimatization state, providing two indices of thermal responses: critical temperature (Tcrit; the temperature above which heart rate abruptly declines) and flatline temperature (Tflat; the temperature where heart rate ceases). Despite the wide use of CTP tests, it remains largely unknown how Tcrit and Tflat change within a single individual after heat acclimation, and whether changes in these indices can predict altered survival in the field. Here, we addressed these issues by evaluating changes in CTP indices in the same individuals before and after heat acclimation. For control mussels, merely reaching Tcrit was not lethal, whereas remaining at Tcrit for ≥10 min was lethal. Heat acclimation significantly increased Tcrit only in mussels with an initially low Tcrit (<35°C), but improved their survival time above Tcrit by 20 min on average. Tflat increased by ∼1.6°C with heat acclimation, but it is unlikely that increased Tflat improves survival in the field. In summary, Tcrit and Tflat per se may fall short of providing quantitative indices of thermal tolerance in mussels; instead, a combination of Tcrit and tolerance time at temperatures ≥Tcrit better defines changes in thermal tolerance with heat acclimation.