Untangling the roles of microclimate, behaviour and physiological polymorphism in governing vulnerability of intertidal snails to heat stress

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.

Incorporating otolith-isotope inferred field metabolic rate into conservation strategies

Fluctuating ocean conditions are rearranging whole networks of marine communities-from individual-level physiological thresholds to ecosystem function. Physiological studies support predictions from individual-level responses (biochemical, cellular, tissue, respiratory potential) based on laboratory experiments. The otolith-isotope method of recovering field metabolic rate has recently filled a gap for the bony fishes, linking otolith stable isotope composition to in situ oxygen consumption and experienced temperature estimates. Here, we review the otolith-isotope method focusing on the biochemical and physiological processes that yield estimates of field metabolic rate. We identify a multidisciplinary pathway in the application of this method, providing concrete research goals (field, modeling) aimed at linking individual-level physiological data to higher levels of biological organization. We hope that this review will provide researchers with a transdisciplinary 'roadmap', guiding the use of the otolith-isotope method to bridge the gap between individual-level physiology, observational field studies, and modeling efforts, while ensuring that in situ data is central in marine policy-making aimed at mitigating climatic and anthropogenic threats.

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.

The impacts of thermal heterogeneity across microhabitats on post-settlement selection of intertidal mussels

Rapid genetic selection is critical for allowing natural populations to adapt to different thermal environments such as those that occur across intertidal microhabitats with high degrees of thermal heterogeneity. To address the question of how thermal regimes influence selection and adaptation in the intertidal black mussel Mytilisepta virgata, we continuously recorded environmental temperatures in both tidal pools and emergent rock microhabitats and then assessed genetic differentiation, gene expression patterns, RNA editing level, and cardiac performance. Our results showed that the subpopulations in the tidal pool and on emergent rocks had different genetic structures and exhibited different physiological and molecular responses to high-temperature stress. These results indicate that environmental heterogeneity across microhabitats is important for driving genetic differentiation and shed light on the importance of post-settlement selection for adaptively modifying the genetic composition and thermal responses of these intertidal mussels.

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.

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.

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.

Thermodynamic effects drive countergradient responses in the thermal performance of Littorina saxatilis across latitude

Thermal performance curves (TPCs) provide a powerful framework to assess the evolution of thermal sensitivity in populations exposed to divergent selection regimes across latitude. However, there is a lack of consensus regarding the extent to which physiological adjustments that compensate for latitudinal temperature variation (metabolic cold adaptation; MCA) may alter the shape of TPCs, including potential repercussion on upper thermal limits. To address this, we compared TPCs for cardiac activity in latitudinally-separated populations of the intertidal periwinkle Littorina saxatilis. We applied a non-linear TPC modelling approach to explore how different metrics governing the shape of TPCs varied systematically in response to local adaptation and thermal acclimation. Both critical upper limits, and the temperatures at which cardiac performance was maximised, were higher in the northernmost (cold-adapted) population and displayed a countergradient latitudinal trend which was most pronounced following acclimation to low temperatures. We interpret this response as a knock-on consequence of increased standard metabolic rate in high latitude populations, indicating that physiological compensation associated with MCA may indirectly influence variation in upper thermal limits across latitude. Our study highlights the danger of assuming that variation in any one aspect of the TPC is adaptive without appropriate mechanistic and ecological context.

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.

Coastal upwelling generates cryptic temperature refugia

Understanding the effects of climate-mediated environmental variation on the distribution of organisms is critically important in an era of global change. We used wavelet analysis to quantify the spatiotemporal (co)variation in daily water temperature for predicting the distribution of cryptic refugia across 16 intertidal sites that were characterized as 'no', 'weak' or 'strong' upwelling and spanned 2000 km of the European Atlantic Coast. Sites experiencing weak upwelling exhibited high synchrony in temperature but low levels of co-variability at monthly to weekly timescales, whereas the opposite was true for sites experiencing strong upwelling. This suggests upwelling generates temporal thermal refugia that can promote organismal performance by both supplying colder water that mitigates thermal stress during hot Summer months and ensuring high levels of fine-scale variation in temperature that reduce the duration of thermal extremes. Additionally, pairwise correlograms based on the Pearson-product moment correlation coefficient and wavelet coherence revealed scale dependent trends in temperature fluctuations across space, with a rapid decay in strong upwelling sites at monthly and weekly timescales. This suggests upwelling also generates spatial thermal refugia that can 'rescue' populations from unfavorable conditions at local and regional scales. Overall, this study highlights the importance of identifying cryptic spatiotemporal refugia that emerge from fine-scale environmental variation to map potential patterns of organismal performance in a rapidly changing world.

Spatial correlation reverses the compound effect of multiple stressors on rocky shore biofilm

Understanding how multifactorial fluctuating environments affect species and communities remains one of the major challenges in ecology. The spatial configuration of the environment is known to generate complex patterns of correlation among multiple stressors. However, to what extent the spatial correlation between simultaneously fluctuating variables affects ecological assemblages in real‐world conditions remains poorly understood. Here, we use field experiments and simulations to assess the influence of spatial correlation of two relevant climate variables – warming and sediment deposition following heavy precipitation – on the biomass and photosynthetic activity of rocky intertidal biofilm. First, we used a response‐surface design experiment to establish the relation between biofilm, warming, and sediment deposition in the field. Second, we used the response surface to generate predictions of biofilm performance under different scenarios of warming and sediment correlation. Finally, we tested the predicted outcomes by manipulating the degree of correlation between the two climate variables in a second field experiment. Simulations stemming from the experimentally derived response surface showed how the degree and direction (positive or negative) of spatial correlation between warming and sediment deposition ultimately determined the nonlinear response of biofilm biomass (but not photosynthetic activity) to fluctuating levels of the two climate variables. Experimental results corroborated these predictions, probing the buffering effect of negative spatial correlation against extreme levels of warming and sediment deposition. Together, these results indicate that consideration of nonlinear response functions and local‐scale patterns of correlation between climate drivers can improve our understanding and ability to predict ecological responses to multiple processes in heterogeneous environments.

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.

Why are "suboptimal" temperatures preferred in a tropical intertidal ectotherm?

In thermally extreme environments it is challenging for organisms to maximize performance due to risks associated with stochastic variation in temperature and, subsequently, over evolutionary time minimizing the exposure to risk can serve as one of the mechanisms that result in organisms preferring suboptimal temperatures. We tested this hypothesis in a slow-moving intertidal snail on tropical rocky shores, where temperature variability increases with time from 30 min to 20 h when recorded at 30 min intervals (due to short-term environmental autocorrelation where temperatures closer in time are more similar as compared to temperatures over a long period of time). Failure to accommodate temporal variation in thermal stress by selecting cool habitats can result in mortality. Thermal performance curves for different traits (heart rate and locomotion) were measured and compared to the snail's thermal preferences in both the field and laboratory. Predicted performances of the snails were simulated based on thermal performance curves for different traits over multiple time scales and simulated carryover effects. A strong mismatch was found between physiological and behavioural thermal maxima of the snails (physiological thermal maximum being higher by ~ 7 °C), but the snails avoided these maxima and sought temperatures 7 - 14 °C cooler. Such a risk-averse strategy can be explained by their predicted performances where the snails should make decisions about preferred temperatures based on time periods ≥ 5 h to avoid underestimating the temporal variation in body temperature. In extreme and stochastic environments, where the temporal variation in environmental conditions can lead to substantial divergence between instantaneous and time-averaged thermal performances, "cooler is better" and "suboptimal" body temperatures are preferred as they provide sufficient buffer to reduce mortality risk from heat stress.

Linking individual and population patterns of rocky-shore mussels

Individual traits and population parameters can be used as proxies of processes taking place within a range of scales, thus improving the way we can evaluate species response to environmental variability. In intertidal rocky shores, patterns at the within-site scale, i.e., between centimeters to hundreds of meters, are important for understanding the population response into these highly variable environments. Here, we studied a rocky-shore mussel population at the within-site spatial scale (1) to test how intertidal height and orientation of the shore affect individual traits and population parameters, (2) to infer the link between individual and population level features, and (3) to explore the upscaling mechanisms driving population structure and processes. We analyzed the patterns of six population parameters: density, biomass, crowding, median individual size, recruitment and mortality rate, and four individual traits: growth rate, spawning phenology, size and condition index. Crowding was defined as the degree of overlapping of individuals within a given area, for which we created a "crowding index". Mussels were studied along the intertidal height gradient in two rocky shores with contrasted orientation at one site over a full year. Our results showed a significant effect of intertidal height and shore orientation on most of individual traits and population parameters studied. In contrast, biomass contained in a full covered surface did not vary in space nor in time. This pattern likely results from relatively constant crowding and a trade-off between median individuals' size and density. We hypothesize that growth, mortality and recruitment rates may all play roles in the stability of the crowding structure of mussel aggregations. Variation in spawning phenology between the two shores in the study site was also observed, suggesting different temporal dynamics of microclimate conditions. Interestingly, despite the different population size distribution between the two shores, our estimates indicate similar potential reproductive output. We hypothesize that the structure of the patches would tend to maintain or carry a maximum of biomass due to trade-offs between density and size while maintaining and maximizing the reproductive output. The patterns of spatial variability of individual traits and population parameters in our study site suggest that heterogeneous within-site conditions influence variation in individual performance and population processes. These results provide insights about the relationship between individual traits and how these relationships make patterns at the population level emerge. They provide baseline information necessary to improve models of metapopulation with spatially explicit processes.

Physiological Consequences of Oceanic Environmental Variation: Life from a Pelagic Organism's Perspective

To better understand life in the sea, marine scientists must first quantify how individual organisms experience their environment, and then describe how organismal performance depends on that experience. In this review, we first explore marine environmental variation from the perspective of pelagic organisms, the most abundant life forms in the ocean. Generation time, the ability to move relative to the surrounding water (even slowly), and the presence of environmental gradients at all spatial scales play dominant roles in determining the variation experienced by individuals, but this variation remains difficult to quantify. We then use this insight to critically examine current understanding of the environmental physiology of pelagic marine organisms. Physiologists have begun to grapple with the complexity presented by environmental variation, and promising frameworks exist for predicting and/or interpreting the consequences for physiological performance. However, new technology needs to be developed and much difficult empirical work remains, especially in quantifying response times to environmental variation and the interactions among multiple covarying factors. We call on the field of global-change biology to undertake these important challenges.

The Goldilocks Principle: A Unifying Perspective on Biochemical Adaptation to Abiotic Stressors in the Sea

The ability of marine organisms to thrive over wide ranges of environmental stressors that perturb structures of proteins, nucleic acids, and lipids illustrates the effectiveness of adaptation at the biochemical level. A critical role of these adaptations is to achieve a proper balance between structural rigidity, which is necessary for maintaining three-dimensional conformation, and flexibility, which is required to allow changes in conformation during function. The Goldilocks principle refers to this balancing act, wherein structural stability and functional properties are poised at values that are just right for the environment the organism faces. Achieving this balance involves changes in macromolecular sequence and adaptive change in the composition of the aqueous or lipid milieu in which macromolecules function. This article traces the development of the field of biochemical adaptation throughout my career and shows how comparative studies of marine animals from diverse habitats have shed light on fundamental properties of life common to all organisms.

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

Elucidating the physiological mechanisms that underlie thermal stress and discovering how species differ in capacities for phenotypic acclimatization and evolutionary adaptation to this stress is critical for understanding current latitudinal and vertical distribution patterns of species and for predicting their future state in a warming world. Such mechanistic analyses require careful choice of study systems (species and temperature-sensitive traits) and design of laboratory experiments that reflect the complexities of in situ conditions. Here, we critically review a wide range of studies of intertidal molluscs that provide mechanistic accounts of thermal effects across all levels of biological organization - behavioural, organismal, organ level, cellular, molecular, and genomic - and show how temperature-sensitive traits govern distribution patterns and capacities for coping with thermal stress. Comparisons of congeners from different thermal habitats are especially effective means for identifying adaptive variation. We employ these mechanistic analyses to illustrate how species differ in the severity of threats posed by rising temperature. Counterintuitively, we show that some of the most heat-tolerant species may be most threatened by increases in temperatures because of their small thermal safety margins and minimal abilities to acclimatize to higher temperatures. We discuss recent molecular biological and genomic studies that provide critical foundations for understanding the types of evolutionary changes in protein structure, RNA secondary structure, genome content, and gene expression capacities that underlie adaptation to temperature. Duplication of stress-related genes, as found in heat-tolerant molluscs, may provide enhanced capacity for coping with higher temperatures. We propose that the anatomical, behavioural, physiological, and genomic diversity found among intertidal molluscs, which commonly are of critical importance and high abundance in these ecosystems, makes this group of animals a highly appropriate study system for addressing questions about the mechanistic determinants of current and future distribution patterns of intertidal organisms.

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.

Limited plasticity in thermally tolerant ectotherm populations: evidence for a trade-off

Many species face extinction risks owing to climate change, and there is an urgent need to identify which species' populations will be most vulnerable. Plasticity in heat tolerance, which includes acclimation or hardening, occurs when prior exposure to a warmer temperature changes an organism's upper thermal limit. The capacity for thermal acclimation could provide protection against warming, but prior work has found few generalizable patterns to explain variation in this trait. Here, we report the results of, to our knowledge, the first meta-analysis to examine within-species variation in thermal plasticity, using results from 20 studies (19 species) that quantified thermal acclimation capacities across 78 populations. We used meta-regression to evaluate two leading hypotheses. The climate variability hypothesis predicts that populations from more thermally variable habitats will have greater plasticity, while the trade-off hypothesis predicts that populations with the lowest heat tolerance will have the greatest plasticity. Our analysis indicates strong support for the trade-off hypothesis because populations with greater thermal tolerance had reduced plasticity. These results advance our understanding of variation in populations' susceptibility to climate change and imply that populations with the highest thermal tolerance may have limited phenotypic plasticity to adjust to ongoing climate warming.

Where three snail species attach while emersed in relation to heterogenous substrate temperatures underneath intertidal boulders

Mobile intertidal gastropods can employ behavioural thermoregulation to mitigate thermal stress, which may include retreating under boulders when emersed. However, little is known about how gastropod occupancy of under-boulder habitats is associated with any variations in substrate temperature that exist under boulders. Thermal imagery was used to measure the temperature of boulder lower surfaces and investigate how three snail species were associated at low tide with the maximum and average temperatures underneath grey siltstone and quartzite. Lower boulder surfaces had heterogeneous temperatures, with grey siltstone having temperature gradients and quartzite temperature showing mosaics. Temperature differences between the hottest and coolest gradient or mosaic locations were >5 °C; thus there was a range of temperatures that snails could interact with. All three snail species occupied cooler parts of temperature mosaics or gradients, avoiding the hottest areas. Stronger associations were detected on the hotter grey siltstone and for the more-thermally sensitive Nerita atramentosa and Diloma concameratum. Even though snails were associated with cooler areas, some individuals were still exposed to extreme substratum heat (>50 °C). These results suggest that gastropod thermoregulatory behaviour is far more complex than simply retreating underneath boulders at low tide, as there is also a range of under-boulder temperatures that they interact with. Untangling interactions between intertidal gastropods and heterogenous substrate temperatures is important given rocky seashores already represent a thermally-variable and potentially-stressful habitat, which may be exacerbated further given predictions of warming temperatures associated with climate change.

Physiological determinants of biogeography: The importance of metabolic depression to heat tolerance

A quantitative understanding of physiological thermal responses is vital for forecasting species distributional shifts in response to climate change. Many studies have focused on metabolic rate as a global metric for analyzing the sublethal effects of changing environments on physiology. Thermal performance curves (TPCs) have been suggested as a viable analytical framework, but standard TPCs may not fully capture physiological responses, due in part to failure to consider the process of metabolic depression. We derived a model based on the nonlinear regression of biological temperature-dependent rate processes and built a heart rate data set for 26 species of intertidal molluscs distributed from 33°S to ~40°N. We then calculated physiological thermal performance limits with continuous heating using T 1 / 2 H , the temperature at which heart rate is decreased to 50% of the maximal rate, as a more realistic measure of upper thermal limits. Results indicate that heat-induced metabolic depression of cardiac performance is a common adaptive response that allows tolerance of harsh environments. Furthermore, our model accounted for the high inter-individual variability in the shape of cardiac TPCs. We then used these TPCs to calculate physiological thermal safety margins (pTSM), the difference between the maximal operative temperature (95th percentile of field temperatures) and T 1 / 2 H of each individual. Using pTSMs, we developed a physiological species distribution model (pSDM) to forecast future geographic distributions. pSDM results indicate that climate-induced species range shifts are potentially less severe than predicted by a simple correlative SDM. Species with metabolic depression below the optimum temperature will be more thermal resistant at their warm trailing edges. High intraspecific variability further suggests that models based on species-level vulnerability to environmental change may be problematic. This multi-scale, mechanistic understanding that incorporates metabolic depression and inter-individual variability in thermal response enables better predictions about the relationship between thermal stress and species distributions.

Resilience in Greenland intertidal Mytilus: The hidden stress defense

The Arctic is experiencing particularly rapid rates of warming, consequently invasive boreal species are now able to survive the less extreme Arctic winter temperatures. Whilst persistence of intertidal and terrestrial species in the Arctic is primarily determined by their ability to tolerate the freezing winters, air temperatures in the Arctic summer can reach 36 °C in the intertidal, which is beyond the upper thermal limits of many marine species. This is normally lethal for the conspicuous ecosystem engineer Mytilus edulis. Transcriptomic analyses were undertaken on both in situ collected and experimentally warmed animals to understand whether M. edulis is able to tolerate these very high summer temperatures. Surprisingly there was no significant enrichment for Gene Ontology terms (GO) when comparing the inner and outer fjord intertidal animals with outer fjord subtidal (control) animals, representing animals collected at 27 °C, 19 °C and 3 °C respectively. This lack of differentiation indicated a wide acclimation ability in this species. Conversely, significant enrichment for processes such as signal transduction, cytoskeleton and cellular protein modification was identified in the expression profiles of the 22 °C and 32 °C experimentally heated animals. This difference in gene expression between in situ collected and experimentally warmed animals was almost certainly due to the former being acclimated to a fluctuating, but predictable, temperature regime, which has increased their thermal tolerances. Interestingly, there was no evidence for enrichment of the classical cellular stress response in any of the animals sampled. Identification of a massive expansion of the HSPA12 heat shock protein 70 kDa gene family presented the possibility of these genes acting as intertidal regulators underpinning thermal resilience. This expansion has resulted in a modified cellular stress response, as an evolutionary adaptation to the rigour of the invasive intertidal life style. Thus, M. edulis appear to have considerable capacity to withstand the current rates of Arctic warming, and the very large attendant thermal variation.

Temperature adaptations of the thermophilic snail Echinolittorina malaccana: insights from metabolomic analysis

The periwinkle snail Echinolittorina malaccana, for which the upper lethal temperature is near 55°C, is one of the most heat-tolerant eukaryotes known. We conducted a multi-level investigation - including cardiac physiology, enzyme activity, and targeted and untargeted metabolomic analyses - that elucidated a spectrum of adaptations to extreme heat in this organism. All systems examined showed heat intensity-dependent responses. Under moderate heat stress (37-45°C), the snail depressed cardiac activity and entered a state of metabolic depression. The global metabolomic and enzymatic analyses revealed production of metabolites characteristic of oxygen-independent pathways of ATP generation (lactate and succinate) in the depressed metabolic state, which suggests that anaerobic metabolism was the main energy supply pathway under heat stress (37-52°C). The metabolomic analyses also revealed alterations in glycerophospholipid metabolism under extreme heat stress (52°C), which likely reflected adaptive changes to maintain membrane structure. Small-molecular-mass organic osmolytes (glycine betaine, choline and carnitine) showed complex changes in concentration that were consistent with a role of these protein-stabilizing solutes in protection of the proteome under heat stress. This thermophilic species can thus deploy a wide array of adaptive strategies to acclimatize to extremely high temperatures.

Impact of short-term elevated temperature stress on winter-acclimated individuals of the marine gastropod Crepidula fornicata

The intertidal zone is an especially stressful thermal habitat, typically exposing residents to air temperatures for up to 6 h at a time, twice daily. Tolerance to elevated temperatures has been particularly well-studied for a variety of intertidal species, especially with regard to upper thermal limits during summers. However, in recent years, as climates have been changing around the world, temperate zone intertidal organisms have sometimes been exposed to periods of unusually high air temperatures during the winter. The present study sought to examine the impact of elevated temperatures on survival and clearance rates of winter-acclimated intertidal individuals of the sedentary marine suspension-feeding gastropod Crepidula fornicata. Individuals were collected intertidally from Nahant, Massachusetts from late January to early April each year for 5 years, maintained in the laboratory at the acclimation temperature of 6 °C, and exposed in the laboratory for 3 h to temperatures as high as 37 °C in seawater either once or twice, 24 h apart. Although mean clearance rates were substantially reduced for at least the next 12-24 h after individuals were returned to the 6 °C control condition following exposures to elevated temperatures as low as 21-26 °C, we saw little mortality even following two 3 h exposures to 35 °C, or single exposures to 37 °C. Mortality was substantial, however, in one experiment following a double exposure to 37 °C. Smaller individuals (~5-12 mm in shell length) were somewhat more sensitive to the thermal stress than adults were. Intertidal members of C. fornicata in Massachusetts seem well-prepared to deal with the increasing range of winter air temperatures associated with the global climate confusion predicted for future years. Additional studies will be required to understand the physiological and biochemical mechanisms used by winter-acclimated individuals of this species to tolerate such periodic substantial temperature increases of 29-31 °C.

Supratidal existence drives phenotypic divergence, but not speciation, in tropical rocky-shore snails

The mechanisms underlying the evolutionary adaptation of animals that transcend the ecological barrier separating the intertidal and supratidal zones of rocky shores are poorly understood. Different wetting frequencies in these zones in tropical regions (daily vs. seasonally, respectively) impose different physical stressors, which should drive phenotypic variation and ultimately speciation in the animals that inhabit them. We studied morphological, physiological and genetic variation in a tropical high-shore gastropod that transcends these zones [Echinolittorina malaccana (Philippi, 1847)]. Variation in melanization, shell features and evaporative water loss was linked to regular seawater wetting, frequent activity and feeding, and solar exposure in intertidal snails, and to inactivity and prolonged aestivation in the shade during continuous air exposure in supratidal snails. Despite selective pressure for phenotypic divergence, and reproductive isolation of the populations in either zone, their mitochondrial COI gene sequences confirmed that they represent a single species. Speciation in our study case is probably constrained by the limitation on activity, mating and reproduction of supratidal snails, such that their populations can only be sustained through intertidal pelagic larval recruitment. Comparisons with other studies suggest that supratidal speciation and specialization for life in this zone probably require moderation of the abiotic (desiccative) conditions, to facilitate greater activity and interaction of animals during air exposure.

Differences in heat tolerance plasticity between supratidal and intertidal snails indicate complex responses to microhabitat temperature variation

Tropical intertidal gastropods that experience extreme and highly variable daily temperatures have evolved significant and complex heat tolerance plasticity, comprising components that respond to different timescales of temperature variation. An earlier study showed different plasticity attributes in snails from differently-heated coastlines, suggesting lifelong irreversible responses that matched habitat thermal regimes. To determine whether heat tolerance plasticity varied at a finer, within-shore spatial scale, we compared the responses of supratidal (predominantly shade-dwelling) and intertidal (frequently solar-exposed) populations of the tropical thermophilic gastropod, Echinolittorina malaccana. Snails modified lethal temperature (LT₅₀) under warm or cool laboratory acclimation, with the overall variation in LT₅₀ being greater in the supratidal (56.0-58.0 °C) than in the intertidal population (57.1-58.1 °C). Similar maximum LT₅₀s expressed by the populations after warm acclimation suggest a capacity limitation under these temperature conditons. The different minimum LT₅₀s after cool acclimation corresponded with microhabitat temperature and field acclimatization of the snails. Different responses to the same laboratory acclimation treatment imply long-term (and possibly lifelong) thermal acclimatization, which could benefit sedentary organisms that are randomly recruited as larvae from a common thermally-stable aquatic environment to thermally-unpredictable intertidal microhabitats. These findings provide another example of thermal tolerance plasticity operating at microhabitat scales, suggesting the importance of considering microhabitat thermal responses when assessing broad-scale environmental change.

Mussel acclimatization to high, variable temperatures is lost slowly upon transfer to benign conditions

Climate change is increasing the temperature variability animals face, and thermal acclimatization allows animals to adjust adaptively to this variability. Although the rate of heat acclimatization has received some study, little is known about how long these adaptive changes remain without continuing exposure to heat stress. This study explored the rate at which field acclimatization states are lost when temperature variability is minimized during constant submersion. California mussels (Mytilus californianus) with different acclimatization states were collected from high- and low-zone sites (∼12 versus ∼5°C daily temperature ranges, respectively) and then kept submerged at 15°C for 8 weeks. Each week, the cardiac thermal performance of mussels was measured as a metric of acclimatization state: critical (Tcrit) and flatline (Tflat) temperatures were recorded. Over 8 weeks of constant submersion, the mean Tcrit of high-zone mussels decreased by 1.07°C from baseline, but low-zone mussels' mean Tcrit was unchanged. High- and low-zone mussels' mean maximum heart rate (HR) and resting HR decreased ∼12 and 35%, respectively. Tflat was unchanged in both groups. These data suggest that Tcrit and HR are more physiologically plastic in response to the narrowing of an animal's daily temperature range than Tflat is, and that an animal's prior acclimatization state (high versus low) influences the acclimatory capacity of Tcrit Approximately 2 months were required for the cardiac thermal performance of the high-zone mussels to reach that of the low-zone mussels, suggesting that acclimatization to high and variable temperatures may persist long enough to enable these animals to cope with intermittent bouts of heat stress.

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.

Effects of balancing selection and microhabitat temperature variations on heat tolerance of the intertidal black mussel Septifer virgatus

Realistic assessments of the impacts of global warming on population extinction risk are likely to require an integrated analysis of the roles of standing genetic variation, microhabitat thermal complexity, and the inter-individual variation of heat tolerance due to both genetic differences and seasonal acclimatization effects. Here, we examine whether balancing selection and microhabitat temperature heterogeneity can interact to enhance the population persistence to thermal stress for the black mussel Septifer virgatus. We deployed biomimetic data loggers on the shore to measure the microhabitat-specific thermal variation from June 2014 to April 2016. Thermal tolerance of specimens was indexed by measuring effects of temperature on heart rate. Genotyping of specimens was performed using double digestion restriction association RADSeq (ddRADseq). Our results show that inter-individual variations in thermal tolerance correlate significantly with genetic differences at some specific gene loci, and that heterozygotes have higher thermal tolerances than homozygotes. The observed seasonal changes in genotype frequency suggest that these loci are under balancing selection. The ability of thermally resistant heterozygotes to survive in sun-exposed microhabitats acts to balance the loss of homozygotes during summer and enable the persistence of genetic polymorphisms. Population persistence of the mussel is also facilitated by the micro-scale variation in temperature, which provides refugia from thermal stress. Our results emphasize that inter-individual variation in thermal tolerance and in microhabitat heterogeneity in temperature are important for the persistence of populations in rocky shore habitats.

The cellular stress response and temperature: Function, regulation, and evolution

The cellular stress response (CSR) is critical for enabling organisms to cope with thermal damage to proteins, nucleic acids, and membranes. It is a graded response whose properties vary with the degree of cellular damage. Molecular damage has positive, as well as negative, function-perturbing effects. Positive effects include crucial regulatory interactions that orchestrate involvement of the different components of the CSR. Thermally unfolded proteins signal for rapid initiation of transcription of genes encoding heat shock proteins (HSPs), central elements of the heat shock response (HSR). Thermal disruption of messenger RNA (mRNA) secondary structures in untranslated regions leads to the culling of the mRNA pool: thermally labile mRNAs for housekeeping proteins are degraded by exonucleases; heat-resistant mRNAs for stress proteins like HSPs then can monopolize the translational apparatus. Thus, proteins and RNA function as "cellular thermometers," and evolved differences in their thermal stabilities enable rapid initiation of the CSR whenever cell temperature rises significantly above the normal thermal range of a species. Covalent DNA damage, which may result from increased production of reactive oxygen species, is temperature-dependent; its extent may determine cellular survival. High levels of stress that exceed capacities for molecular repair can lead to proteolysis, inhibition of cell division, and programmed cell death (apoptosis). Onset of these processes may occur later in the stress period, after initiation of the HSR, to allow HSPs opportunity to restore protein homeostasis. Delay of these energy costly processes may also result from shortfalls in availability of adenosine triphosphate and reducing power during times of peak stress.

Local cold adaption increases the thermal window of temperate mussels in the Arctic

Species expand towards higher latitudes in response to climate warming, but the pace of this expansion is related to the physiological capacity to resist cold stress. However, few studies exist that have quantified the level of inter-population local adaptation in marine species freeze tolerance, especially in the Arctic. We investigated the importance of cold adaptation and thermal window width towards high latitudes from the temperate to the Arctic region. We measured upper and lower lethal air temperatures (i.e. LT and LT50) in temperate and Arctic populations of blue mussels (Mytilus edulis), and analysed weather data and membrane fatty acid compositions, following emersion simulations. Both populations had similar upper LT (~38 °C), but Arctic mussels survived 4°C colder air temperatures than temperate mussels (-13 vs. -9°C, respectively), corresponding to an 8% increase in their thermal window. There were strong latitudinal relationships between thermal window width and local air temperatures, indicating Arctic mussels are highly adapted to the Arctic environment where the seasonal temperature span exceeds 60°C. Local adaptation and local habitat heterogeneity thus allow leading-edge M. edulis to inhabit high Arctic intertidal zones. This intraspecific pattern provides insight into the importance of accounting for cold adaptation in climate change, conservation and biogeographic studies.

Thermotolerance Divergence Revealed by the Physiological and Molecular Responses in Two Oyster Subspecies of Crassostrea gigas in China

Investigating the physiological mechanisms of closely related species that exhibit distinct geographic distributions and thermal niches is essential for understanding their thermal tolerance capacities and local adaptations in view of climate warming. The variations in upper thermal limits (LT50) under acute heat shock and cardiac activity, standard metabolic rate (SMR), anaerobic metabolite production and molecular responses (expression of molecular chaperones and glycolysis metabolism genes) under increasing temperatures in two oyster subspecies were studied. The populations of two oyster subspecies, Crassostrea gigas gigas and C. gigas angulata, exhibit different latitudinal distributions along the northern and southern coastlines of China, respectively, which experience different environmental conditions. The LT50 was significantly higher, by ∼1°C, in the southern than in the northern oysters. In both subspecies, temperature increases had powerful effects on heart rate, SMR and gene expression. The southern oysters had the highest Arrhenius breakpoint temperatures for heart rate (31.4 ± 0.17°C) and SMR (33.09°C), whereas the heart rate (28.86 ± 0.3°C) and SMR (29.22°C) of the northern oysters were lower. The same patterns were observed for the Q 10 coefficients. More thermal sensitivity was observed in the northern oysters than in their southern counterparts, as the heat-shock proteins (HSPs) in the northern oysters were expressed first and had a higher induction at a lower temperature than those of southern oysters. Furthermore, different expression patterns of energetic metabolism genes (HK, PK, and PEPCK) were observed. In the northern oysters, increasing anaerobic glycolysis genes (PEPCK) and end products (succinate) were found at 36-43°C, indicating a transition from aerobic to anaerobic metabolism and a lower aerobic scope compared with the southern oysters. These two subspecies experience different environmental conditions, and their physiological performances suggested species-specific thermal tolerance windows in which the southern oysters, with mild physiological flexibility, had a higher potential capability to withstand heat stress. Overall, our results indicate that comparing and unifying physiological and molecular mechanisms can provide a framework for understanding the likely effects of global warming on marine ectotherms in intertidal regions.

Mapping physiology: biophysical mechanisms define scales of climate change impacts

The rocky intertidal zone is a highly dynamic and thermally variable ecosystem, where the combined influences of solar radiation, air temperature and topography can lead to differences greater than 15°C over the scale of centimetres during aerial exposure at low tide. For most intertidal organisms this small-scale heterogeneity in microclimates can have enormous influences on survival and physiological performance. However, the potential ecological importance of environmental heterogeneity in determining ecological responses to climate change remains poorly understood. We present a novel framework for generating spatially explicit models of microclimate heterogeneity and patterns of thermal physiology among interacting organisms. We used drone photogrammetry to create a topographic map (digital elevation model) at a resolution of 2 × 2 cm from an intertidal site in Massachusetts, which was then fed into to a model of incident solar radiation based on sky view factor and solar position. These data were in turn used to drive a heat budget model that estimated hourly surface temperatures over the course of a year (2017). Body temperature layers were then converted to thermal performance layers for organisms, using thermal performance curves, creating 'physiological landscapes' that display spatially and temporally explicit patterns of 'microrefugia'. Our framework shows how non-linear interactions between these layers lead to predictions about organismal performance and survivorship that are distinct from those made using any individual layer (e.g. topography, temperature) alone. We propose a new metric for quantifying the 'thermal roughness' of a site (RqT, the root mean square of spatial deviations in temperature), which can be used to quantify spatial and temporal variability in temperature and performance at the site level. These methods facilitate an exploration of the role of micro-topographic variability in driving organismal vulnerability to environmental change using both spatially explicit and frequency-based approaches.

Forecasting species range dynamics with process-explicit models: matching methods to applications

Knowing where species occur is fundamental to many ecological and environmental applications. Species distribution models (SDMs) are typically based on correlations between species occurrence data and environmental predictors, with ecological processes captured only implicitly. However, there is a growing interest in approaches that explicitly model processes such as physiology, dispersal, demography and biotic interactions. These models are believed to offer more robust predictions, particularly when extrapolating to novel conditions. Many process-explicit approaches are now available, but it is not clear how we can best draw on this expanded modelling toolbox to address ecological problems and inform management decisions. Here, we review a range of process-explicit models to determine their strengths and limitations, as well as their current use. Focusing on four common applications of SDMs - regulatory planning, extinction risk, climate refugia and invasive species - we then explore which models best meet management needs. We identify barriers to more widespread and effective use of process-explicit models and outline how these might be overcome. As well as technical and data challenges, there is a pressing need for more thorough evaluation of model predictions to guide investment in method development and ensure the promise of these new approaches is fully realised.

Temperature-induced aerobic scope and Hsp70 expression in the sea cucumber Holothuria scabra

The Aerobic Scope (AS), which reflects the functional capacity for biological fitness, is a highly relevant proxy to determine thermal tolerance in various taxa. Despite the importance of this method, its implementation is often hindered, due to lacking techniques to accurately measure standard- (SMR) and maximal- (MMR) metabolic rates, especially in sluggish marine invertebrates with low oxygen consumption rates, such as sea cucumbers. In this study the AS concept was modified to define a Temperature-induced Aerobic Scope (TAS), based on metabolic rate changes due to temperature adjustments rather than traditionally used physical activity patterns. Consequentially, temperature dependent peak and bottom O₂ consumption rates, defined as Temperature-induced Maximal- (TMMR) and Standard Metabolic Rates (TSMR), respectively, served as MMR and SMR alternatives for the sea cucumber Holothuria scabra. TMMR and TSMR were induced through acute temperature change (2°C per hour; 17–41°C) until critical warm (WT_crit) and cold (CT_crit) temperatures were reached, respectively. In addition, Hsp70 gene expression linked to respiration rates served as synergistic markers to confirm critical threshold temperatures. O₂ consumption of H. scabra peaked distinctly at WT_crit of 38°C (TMMR = 33.2 ± 4.7 μgO₂ g^-1 h^-1). A clear metabolic bottom line was reached at CT_crit of 22°C (TSMR = 2.2 ± 1.4 μgO₂ g^-1 h^-1). Within the thermal window of 22–38°C H. scabra sustained positive aerobic capacity, with assumed optimal performance range between 29–31.5°C (13.85–18.7 μgO₂ g^-1 h^-1). Between 39–41°C H. scabra decreased respiration progressively, while gene expression levels of Hsp70 increased significantly at 41°C, indicating prioritization of heat shock response (HSR) and homeostatic disruption. At the cold end (17–22°C) homeostatic disruption was visible through incrementally increasing energetic expenses to fuel basal maintenance costs, but no Hsp70 overexpression occurred. TMMR, TSMR and TAS proved to be reliable metrics, similar to the traditional energetic key parameters MMR, SMR and AS, to determine a specific aerobic performance window for the sluggish bottom dwelling species H. scabra. In addition, the linkage between respiration physiology and molecular defense mechanisms showed valuable analytical synergies in terms of mechanistic prioritization as response to thermal stress. Overall, this study will help to define lethal temperatures for aquaculture and to predict the effects of environmental stress, such as ocean warming, in H. scabra.

Non-reversible and Reversible Heat Tolerance Plasticity in Tropical Intertidal Animals: Responding to Habitat Temperature Heterogeneity

The theory for thermal plasticity of tropical ectotherms has centered on terrestrial and open-water marine animals which experience reduced variation in diurnal and seasonal temperatures, conditions constraining plasticity selection. Tropical marine intertidal animals, however, experience complex habitat thermal heterogeneity, circumstances encouraging thermal plasticity selection. Using the tropical rocky-intertidal gastropod, Echinolittorina malaccana, we investigated heat tolerance plasticity in terms of laboratory acclimation and natural acclimatization of populations from thermally-dissimilar nearby shorelines. Laboratory treatments yielded similar capacities of snails from either population to acclimate their lethal thermal limit (LT50 variation was ∼2°C). However, the populations differed in the temperature range over which acclimatory adjustments could be made; LT50 plasticity occurred over a higher temperature range in the warm-shore snails compared to the cool-shore snails, giving an overall acclimation capacity for the populations combined of 2.9°C. In addition to confirming significant heat tolerance plasticity in tropical intertidal animals, these findings reveal two plasticity forms, reversible (laboratory acclimation) and non-reversible (population or shoreline specific) plasticity. The plasticity forms should account for different spatiotemporal scales of the environmental temperature variation; reversible plasticity for daily and tidal variations in microhabitat temperature and non-reversible plasticity for lifelong, shoreline temperature conditions. Non-reversible heat tolerance plasticity, likely established after larvae settle on the shore, should be energetically beneficial in preventing heat shock protein overexpression, but also should facilitate widespread colonization of coasts that support thermally-diverse shorelines. This first demonstration of different plasticity forms in benthic intertidal animals supports the hypothesis that habitat heterogeneity (irrespective of latitude) drives thermal plasticity selection. It further suggests that studies not making reference to different spatial scales of thermal heterogeneity, nor seeking how these may drive different thermal plasticity forms, risk misinterpreting ectothermic responses to environmental warming.

High thermal stress responses of Echinolittorina snails at their range edge predict population vulnerability to future warming

Populations at the edge of their species' distribution ranges are typically living at the physiological extreme of the environmental conditions they can tolerate. As a species' response to global change is likely to be largely determined by its physiological performance, subsequent changes in environmental conditions can profoundly influence populations at range edges, resulting in range extensions or retractions. To understand the differential physiological performance among populations at their distribution range edge and center, we measured levels of mRNA for heat shock protein 70 (hsp70) as an indicator of temperature sensitivity in two high-shore littorinid snails, Echinolittorina malaccana and E. radiata, between 1°N to 36°N along the NW Pacific coast. These Echinolittorina snails are extremely heat-tolerant and frequently experience environmental temperatures in excess of 55 °C when emersed. It was assumed that animals exhibiting high temperature sensitivity will synthesize higher levels of mRNA, which will thus lead to higher energetic costs for thermal defense. Populations showed significant geographic variation in temperature sensitivity along their range. Snails at the northern range edge of E. malaccana and southern range edge of E. radiata exhibited higher levels of hsp70 expression than individuals collected from populations at the center of their respective ranges. The high levels of hsp70 mRNA in populations at the edge of a species' distribution range may serve as an adaptive response to locally stressful thermal environments, suggesting populations at the edge of their distribution range are potentially more sensitive to future global warming.

Untangling the roles of microclimate, behaviour and physiological polymorphism in governing vulnerability of intertidal snails to heat stress

Biogeographic distributions are driven by cumulative effects of smaller scale processes. Thus, vulnerability of animals to thermal stress is the result of physiological sensitivities to body temperature (T_b), microclimatic conditions, and behavioural thermoregulation. To understand interactions among these variables, we analysed the thermal tolerances of three species of intertidal snails from different latitudes along the Chinese coast, and estimated potential T_b in different microhabitats at each site. We then empirically determined the temperatures at which heart rate decreased sharply with rising temperature (Arrhenius breakpoint temperature, ABT) and at which it fell to zero (flat line temperature, FLT) to calculate thermal safety margins (TSM). Regular exceedance of FLT in sun-exposed microhabitats, a lethal effect, was predicted for only one mid-latitude site. However, ABTs of some individuals were exceeded at sun-exposed microhabitats in most sites, suggesting physiological impairment for snails with poor behavioural thermoregulation and revealing inter-individual variations (physiological polymorphism) of thermal limits. An autocorrelation analysis of T_b showed that predictability of extreme temperatures was lowest at the hottest sites, indicating that the effectiveness of behavioural thermoregulation is potentially lowest at these sites. These results illustrate the critical roles of mechanistic studies at small spatial scales when predicting effects of climate change.