Micro-scale environmental variation amplifies physiological variation among individual mussels

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.

Multidimensional plasticity jointly contributes to rapid acclimation to environmental challenges during biological invasions

Rapid plastic response to environmental changes, which involves extremely complex underlying mechanisms, is crucial for organismal survival during many ecological and evolutionary processes such as those in global change and biological invasions. Gene expression is among the most studied molecular plasticity, while co- or posttranscriptional mechanisms are still largely unexplored. Using a model invasive ascidian Ciona savignyi, we studied multidimensional short-term plasticity in response to hyper- and hyposalinity stresses, covering the physiological adjustment, gene expression, alternative splicing (AS), and alternative polyadenylation (APA) regulations. Our results demonstrated that rapid plastic response varied with environmental context, timescales, and molecular regulatory levels. Gene expression, AS, and APA regulations independently acted on different gene sets and corresponding biological functions, highlighting their nonredundant roles in rapid environmental adaptation. Stress-induced gene expression changes illustrated the use of a strategy of accumulating free amino acids under high salinity and losing/reducing them during low salinity to maintain the osmotic homoeostasis. Genes with more exons were inclined to use AS regulations, and isoform switches in functional genes such as SLC2a5 and Cyb5r3 resulted in enhanced transporting activities by up-regulating the isoforms with more transmembrane regions. The extensive 3'-untranslated region (3'UTR) shortening through APA was induced by both salinity stresses, and APA regulation predominated transcriptomic changes at some stages of stress response. The findings here provide evidence for complex plastic mechanisms to environmental changes, and thereby highlight the importance of systemically integrating different levels of regulatory mechanisms in studying initial plasticity in evolutionary trajectories.

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.

Environment-driven shifts in inter-individual variation and phenotypic integration within subnetworks of the mussel transcriptome and proteome

The environment can alter the magnitude of phenotypic variation among individuals, potentially influencing evolutionary trajectories. However, environmental influences on variation are complex and remain understudied. Populations in heterogeneous environments might exhibit more variation, the amount of variation could differ between benign and stressful conditions, and/or variation might manifest in different ways among stages of the gene‐to‐protein expression cascade or among physiological functions. Here, we explore these three issues by quantifying patterns of inter‐individual variation in both transcript and protein expression levels among California mussels, Mytilus californianus Conrad. Mussels were exposed to five ecologically relevant treatments that varied in the mean and interindividual heterogeneity of body temperature. To target a diverse set of physiological functions, we assessed variation within 19 expression subnetworks, including canonical stress‐response pathways and empirically derived coexpression clusters that represent a diffuse set of cellular processes. Variation in expression was particularly pronounced in the treatments with high mean and heterogeneous body temperatures. However, with few exceptions, environment‐dependent shifts of variation in the transcriptome were not reflected in the proteome. A metric of phenotypic integration provided evidence for a greater degree of constraint on relative expression levels (i.e., stronger correlation) within expression subnetworks in benign, homogeneous environments. Our results suggest that environments that are more stressful on average – and which also tend to be more heterogeneous – can relax these expression constraints and reduce phenotypic integration within biochemical subnetworks. Context‐dependent “unmasking” of functional variation may contribute to interindividual differences in physiological phenotype and performance in stressful environments.

Interactions Between Temperature Variability and Reproductive Physiology Across Traits in an Intertidal Crab

Thermal extremes alter population processes, which can result in part from temperature-induced movement at different spatial and temporal scales. Thermal thresholds for animal movement likely change based on underlying thermal physiology and life-history stage, a topic that requires greater study. The intertidal porcelain crab Petrolisthes cinctipes currently experiences temperatures that can reach near-lethal levels in the high-intertidal zone at low tide. However, the thermal thresholds that trigger migration to cooler microhabitats, and the extent to which crabs move in response to temperature, remain unknown. Moreover, the influence of reproductive status on these thresholds is rarely investigated. We integrated demographic, molecular, behavioral, and physiological measurements to determine if behavioral thermal limits varied due to reproductive state. Demographic data showed a trend for gravid, egg bearing, crabs to appear more often under rocks in the cooler intertidal zone where crab density is highest. In situ expression of 31 genes related to stress, metabolism, and growth in the field differed significantly based on intertidal elevation, with mid-intertidal crabs expressing the gene for the reproductive yolk protein vitellogenin (vg) earlier in the season. Furthermore, VG protein levels were shown to increase with density for female hemolymph. Testing for temperatures that elicit movement revealed that gravid females engage in heat avoidance behavior at lower temperatures (i.e., have a lower voluntary thermal maximum, VTmax) than non-gravid females. VTmax was positively correlated with the temperature of peak firing rate for distal afferent nerve fibers in the walking leg, a physiological relationship that could correspond to the mechanistic underpinning for temperature dependent movement. The vulnerability of marine organisms to global change is predicated by their ability to utilize and integrate physiological and behavioral strategies in response to temperature to maximize survival and reproduction. Interactions between fine-scale temperature variation and reproductive biology can have important consequences for the ecology of species, and is likely to influence how populations respond to ongoing 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.

Acclimation to future climate exposes vulnerability to cold extremes in intertidal sea hares

Highly dynamic environments like estuaries will undergo unpredictable shifts in thermal and salinity regimes with ongoing climate change. These interactive stressors fluctuate predictably and seasonally over historical periods, which has facilitated the evolution of wide environmental tolerance in some estuarine inhabitants. However, physiological and behavioral acclimatization is seasonally based for many estuarine species, meaning that a shift in the unpredictability of climate events and trends will disrupt the effectiveness of evolved tolerance mechanisms. Of particular concern are extreme cold events and high-volume precipitation events, which will acutely and unpredictably alter an estuarine habitat. The eelgrass sea hare, Phyllaplysia taylori, has documented euryhaline and eurythermal tolerance to summer conditions, but the winter environment may pose a greater challenge to seasonally relevant acclimatization scenarios. Here, we characterized lower critical thermal limits, and behavioral responses to stimuli leading up to these limits, in two central California P. taylori populations under four temperature-salinity scenarios in a laboratory acclimation experiment. Acclimation to warmer conditions significantly increased critical thermal minima, while fresher conditions resulted in high mortality. However, the surviving individuals in the fresher conditions were able to respond to stimuli more quickly overall, despite their shortest response time being at a higher temperature than the saltier-acclimated individuals. Within the environmental context of their natural habitats, we find that acclimation to climate change-induced warming will hinder sea hares' ability to weather existing and future cold extremes and precipitation events.

A single heat-stress bout induces rapid and prolonged heat acclimation in the California mussel, Mytilus californianus

Climate change is not only causing steady increases in average global temperatures but also increasing the frequency with which extreme heating events occur. These extreme events may be pivotal in determining the ability of organisms to persist in their current habitats. Thus, it is important to understand how quickly an organism's heat tolerance can be gained and lost relative to the frequency with which extreme heating events occur in the field. We show that the California mussel, Mytilus californianus-a sessile intertidal species that experiences extreme temperature fluctuations and cannot behaviourally thermoregulate-can quickly (in 24-48 h) acquire improved heat tolerance after exposure to a single sublethal heat-stress bout (2 h at 30 or 35°C) and then maintain this improved tolerance for up to three weeks without further exposure to elevated temperatures. This adaptive response improved survival rates by approximately 75% under extreme heat-stress bouts (2 h at 40°C). To interpret these laboratory findings in an ecological context, we evaluated 4 years of mussel body temperatures recorded in the field. The majority (approx. 64%) of consecutive heat-stress bouts were separated by 24-48 h, but several consecutive heat bouts were separated by as much as 22 days. Thus, the ability of M. californianus to maintain improved heat tolerance for up to three weeks after a single sublethal heat-stress bout significantly improves their probability of survival, as approximately 33% of consecutive heat events are separated by 3-22 days. As a sessile animal, mussels likely evolved the capability to rapidly gain and slowly lose heat tolerance to survive the intermittent, and often unpredictable, heat events in the intertidal zone. This adaptive strategy will likely prove beneficial under the extreme heat events predicted with climate change.

Highly dynamic transcriptional reprogramming and shorter isoform shifts under acute stresses during biological invasions

Phenotypic plasticity has been increasingly recognized for its importance in adaptation to novel environments, and initial rapid plastic response to acute stresses usually serves as the stepping stone for future adaptation. Differential gene expression and alternative splicing have been proposed as two underlying mechanisms for rapid plastic response to environmental stresses. Here, we used an invasive model species, Ciona savignyi, to investigate the temporary plastic changes under temperature stresses on gene expression and alternative splicing. Our results revealed rapid and highly dynamic gene expression reprogramming and alternative splicing switch under acute stresses. Distinct transcriptional response profiles were triggered by two types of temperature stresses, showing resilience recovery and increasing divergence under heat and cold challenges, respectively. Interestingly, alternative exons were more inclined to be skipped under both heat and cold stresses, leading to shorter isoforms but with maintained Open Reading Frames (ORFs). Although similar response patterns were observed between differential gene expression and alternative splicing, low overlap between Differentially Expressed Genes (DEGs) and Differentially Alternative Spliced Genes (DASGs) suggests that distinct gene sets and associated functions should be involved in temperature challenges. Thus, alternative splicing should offer an additional layer of plastic response to environmental challenges. Finally, we identified key plastic genes involved in both gene expression regulation and alternative splicing. The results obtained here shed light on adaptation and accommodation mechanisms during biological invasions, particularly for acute environmental changes at early stages of biological invasions such as transport and introduction.

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.

Stress Memory of Recurrent Environmental Challenges in Marine Invasive Species: Ciona robusta as a Case Study

Fluctuating environmental changes impose tremendous stresses on sessile organisms in marine ecosystems, in turn, organisms develop complex response mechanisms to keep adaptive homeostasis for survival. Physiological plasticity is one of the primary lines of defense against environmental challenges, and such defense often relies on the antioxidant defense system (ADS). Hence, it is imperative to understand response mechanisms of ADS to fluctuating environments. Invasive species provide excellent models to study how species cope with environmental stresses, as invasive species encounter sudden, and often recurrent, extensive environmental challenges during the whole invasion process. Here, we studied the roles of ADS on rapid response to recurrent cold challenges in a highly invasive tunicate (Ciona robusta) by simulating cold stresses during its invasion process. We assessed antioxidative indicators, including malondialdehyde (MDA), total antioxidant capacity (T-AOC), superoxide dismutase (SOD), catalase (CAT), and glutathione (GSH), as well as transcriptional changes of ADS-related genes to reveal the physiological plasticity under recurring cold stresses. Our results demonstrated that physiological homeostasis relied on the resilience of ADS, which further accordingly tuned antioxidant activity and gene expression to changing environments. The initial cold stress remodeled baselines of ADS to promote the development of stress memory, and subsequent stress memory largely decreased the physiological response to recurrent environmental challenges. All results here suggest that C. robusta could develop stress memory to maintain physiological homeostasis in changing or harsh environments. The results obtained in this study provide new insights into the mechanism of rapid physiological adaption during biological invasions.

Physiological responses to heat stress in an invasive mussel Mytilus galloprovincialis depend on tidal habitat

Mussels are ecologically important organisms that can survive in subtidal and intertidal zones where they experience thermal stress. We know little about how mussels from different tidal habitats respond to thermal stress. We used the mussel Mytilus galloprovincialis from separate subtidal and intertidal populations to test whether heart rate and indicators of potential aerobic (citrate synthase activity) and anaerobic (cytosolic malate dehydrogenase activity) metabolic capacity are affected by increased temperatures while exposed to air or submerged in water. Subtidal mussels were affected by warming when submerged in water (decreased heart rate) but showed no effect in air. In contrast, intertidal mussels were affected by exposure to air (increased anaerobic capacity) but not by warming. Overall, physiological responses of mussels to thermal stress were dependent on their tidal habitat. These results highlight the importance of considering the natural habitat of mussels when assessing their responses to environmental challenges.

Long-term survival probability, not current habitat quality, predicts dispersal distance in a stream salamander

Dispersal evolves as an adaptive mechanism to optimize individual fitness across the landscape. Specifically, dispersal represents a mechanism to escape fitness costs resulting from changes in environmental conditions. Decades of empirical work suggest that individuals use local habitat cues to make movement decisions, but theory predicts that dispersal can also evolve as a fixed trait, independent of local conditions, in environments characterized by a history of stochastic spatiotemporal variation. Until now, however, both conditional and fixed models of dispersal evolution have primarily been evaluated using emigration data (stay vs. leave), and not dispersal distances: a more comprehensive measure of dispersal. Our goal was to test whether conditional or fixed models of dispersal evolution predict variation in dispersal distance in the stream salamander Gyrinophilus porphyriticus. We quantified variation in habitat conditions using measures of salamander performance from 4 yr of spatially explicit, capture-mark-recapture (CMR) data across three headwater streams in the Hubbard Brook Experimental Forest in central New Hampshire, USA. We used body condition as an index of local habitat quality that individuals may use to make dispersal decisions, and survival probability estimated from multistate CMR models as an index of mortality risk resulting from the long-term history of environmental variation. We found that dispersal distances increased with declining survival probability, indicating that salamanders disperse further in risky environments. Dispersal distances were unrelated to spatial variation in body condition, suggesting that salamanders do not base dispersal distance decisions on local habitat quality. Our study provides the first empirical support for fixed models of dispersal evolution, which predict that dispersal evolves in response to a history of spatiotemporal environmental variation, rather than as a conditional response to current habitat conditions. More broadly, this study underscores the value of assessing alternative scales of environmental variation to gain a more complete and balanced understanding of dispersal evolution.

Hot Rocks and Not-So-Hot Rocks on the Seashore: Patterns and Body-Size Dependent Consequences of Microclimatic Variation in Intertidal Zone Boulder Habitat

Microclimatic variation has emerged as an important driver of many ecological and evolutionary processes. Nonetheless, fine-scale temperature data are still rare in most habitats, limiting our ability to understand the consequences of microclimatic variation under current and future conditions. We measured fine-scale thermal variation in a common, species-rich, but rarely studied habitat with respect to temperature: the airspaces under rocks on intertidal zone boulder shores. The effects of thermal variation were investigated using physiological, behavioral, and demographic responses of the porcelain crab Petrolisthes cinctipes. Habitat temperatures were measured at fine spatial and temporal resolution over 18 months, producing 424,426 temperature records. Microclimatic variation increased with increasing intertidal elevation, particularly with respect to heat extremes. However, mean temperatures were similar across the entire intertidal zone. Overheating risk for P. cinctipes increases with intertidal elevation but is size dependent, as large animals are more heat sensitive than small animals. Still, microclimatic variation high in the intertidal zone provided thermal refugia even under the warmest conditions. Size-dependent thermal responses predicted that large crabs should be rare high in the intertidal zone, which was supported by demographic data. Furthermore, simulations parameterized by our microclimate and organismal data recapitulated demographic patterns. Therefore, interactions between microclimatic variation and size-dependent thermal responses may have significant ecological repercussions that warrant greater attention.

Inter-individual physiological variation in responses to environmental variation and environmental change: Integrating across traits and time

Greater understanding of physiological responses to climate change demands deeper comprehension of the causes and consequences of physiological variation. Increasingly, population trait means are being deconstructed into variable signals at the level of individuals. We advocate for greater consideration of such inter-individual physiological variation and how it both depends on and interacts with environmental variability. First, we review several studies on the intertidal mussel Mytilus californianus to illustrate how the magnitude of inter-individual variation may depend on the environmental context analyzed (i.e., is the mean condition benign or stressful?) and/or on the specific physiological metric investigated. Stressful conditions may reveal or mask variation in disparate ways at different levels of analysis (e.g., transcriptome vs. proteome), but we often lack crucial information regarding the relationships among these different physiological metrics and their consequences for fitness. We then reanalyze several published datasets to ask whether individuals employ divergent strategies over time in response to acute heat stress; such time-dependence would further complicate interpretation of physiological variation. However, definitive conclusions are precluded by limited sample sizes and short timescales in extant datasets. A key remaining challenge is to extend these analytical frameworks to longer periods over which individuals in a population experience repeated, but spatially variable, episodic stress events. We conclude that variation at multiple levels of analysis should be investigated over longer periods and, where possible, within individuals (or genotypes) experiencing repeated environmental challenges. Although difficult in practice, such studies will facilitate improved understanding of potential population-level physiological responses to climate change.

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.

Transcriptome analysis of body wall reveals growth difference between the largest and smallest individuals in the pure and hybrid populations of Apostichopus japonicus

Long-term inbreeding of sea cucumber has resulted in a decrease in its growth rate, which has severely affected yield and economic efficiency. In this study, three Apostichopus japonicus families were constructed and screened into the weight of smallest and largest, which included Russian, Chinese, and their hybrids (RC). We examined the transcriptional profiles of hybrid (RC) and purebred (CC and RR). A total of 49.69 Gb clean reads were obtained, and the Q30 base percentage was above 90.47%. A total of 5191 novel genes were discovered, of which 2592 genes were annotated. Differentially expressed genes (DEGs) were identified, and functional annotation and enrichment analysis were performed. Approximately 1874 DEGs were screened in the Chinese sea cucumber (CC) difference group; 2591 DEGs were obtained in the hybrid sea cucumber difference group (RC), and 3006 DEGs were obtained in the Russian sea cucumber difference group (RR). In Gene Ontology (GO) analysis, highest DEG enrichment was observed for the functional categories of cellular process and metabolic process. In terms of cellular components, DEG enrichment was observed in cell part, cell; for molecular function, DEG enrichment was detected in catalytic activity, binding, hydrolase activity, transferase activity. According to the differential expression analysis, we found that 15 heat shock protein (HSP) genes that have the same expression trends, which were upregulated in the smallest weight of three sea cucumber lines. In addition, COG analysis of defense genes was conducted. All defense genes (ATP-binding cassette transporters (ABCs), multidrug resistance protein (MRPs), and beta-lactamase) showed the same expression trend, which was significantly upregulated in smallest individuals compared to that of largest individuals in RC lines, which implied the smallest individuals are exposed to more pressure during growth. These results may lead to the smallest individuals showing slow growth. Additionally, we selected 12 DEGs to validate the result by qPCR. Those DEGs were included in growth-related and resistance genes. Sequencing of the A. japonicus transcriptome improves our understanding of the transcriptional regulatory apparatus that controls individual development and growth.

Plasticity of thermal tolerance and its relationship with growth rate in juvenile mussels (Mytilus californianus)

Complex life cycles characterized by uncertainty at transitions between larval/juvenile and adult environments could favour irreversible physiological plasticity at such transitions. To assess whether thermal tolerance of intertidal mussels (Mytilus californianus) adjusts to post-settlement environmental conditions, we collected juveniles from their thermally buffered microhabitat from high- and low-shore locations at cool (wave-exposed) and warm (wave-protected) sites. Juveniles were transplanted to unsheltered cages at the two low sites or placed in a common garden. Juveniles transplanted to the warm site for one month in summer had higher thermal tolerance, regardless of origin site. By contrast, common-garden juveniles from all sites had lower tolerance indistinguishable from exposed site transplants. After six months in the field plus a common garden period, there was a trend for higher thermal tolerance at the protected site, while reduced thermal tolerance at both sites indicated seasonal acclimatization. Thermal tolerance and growth rate were inversely related after one but not six months; protected-site transplants were more tolerant but grew more slowly. In contrast to juveniles, adults from low-shore exposed and protected sites retained differences in thermal tolerance after common garden treatment in summer. Both irreversible and reversible forms of plasticity must be considered in organismal responses to changing environments.

Thermal acclimation of fast-growing Japanese Quails (Coturnix japonica) exhibit decreased oxidative stress and increased muscle fiber diameters after acute heat challenges

Many predict dire consequences of increasing temperatures; however, high temperatures in early life may aid animals during extreme thermal events later in life. The underlying physiological mechanisms have not been elucidated. We examined whether developing in warmer temperatures would physiologically benefit adult Japanese Quails (Coturnix japonica Temminck and Schlegel, 1849) by exploring changes in oxidative stress and muscle structure in two quail lines — one selected for control growth and another for fast growth and after acute heat challenges. We used a factorial design to administer four treatment combinations to each line: an initial period of either heat-stress acclimation (3 h every other day to 37 °C) or no acclimation, and after 5 weeks, either an acute heat-stress challenge (8 h at 39 °C) or no challenge. We found that control quails had significantly higher citrate synthase activity than fast-growing quails. Fast-growing quails had higher hydroxyl scavenging capacity than control quails. Peroxyl scavenging capacity decreased in both lines after an acute heat challenge, regardless of acclimation. Finally, fast-growing quails had larger muscle fiber diameters than control quails, and acclimated birds that experienced an acute heat challenge had larger muscle fibers than those that did not experience a heat challenge. Thus, fast-growing quails may physiologically benefit from developing in warmer temperatures.

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.

Intraspecific variation in oxidative stress tolerance in a model cnidarian: Differences in peroxide sensitivity between and within populations of Nematostella vectensis

Nematostella vectensis is a member of the phylum Cnidaria, a lineage that includes anemones, corals, hydras, and jellyfishes. This estuarine anemone is an excellent model system for investigating the evolution of stress tolerance because it is easy to collect in its natural habitat and to culture in the laboratory, and it has a sequenced genome. Additionally, there is evidence of local adaptation to environmental stress in different N. vectensis populations, and abundant protein-coding polymorphisms have been identified, including polymorphisms in proteins that are implicated in stress responses. N. vectensis can tolerate a wide range of environmental parameters, and has recently been shown to have substantial intraspecific variation in temperature preference. We investigated whether different clonal lines of anemones also exhibit differential tolerance to oxidative stress. N. vectensis populations are continually exposed to reactive oxygen species (ROS) generated during cellular metabolism and by other environmental factors. Fifteen clonal lines of N. vectensis collected from four different estuaries were exposed to hydrogen peroxide. Pronounced differences in survival and regeneration were apparent between clonal lines collected from Meadowlands, NJ, Baruch, SC, and Kingsport, NS, as well as among 12 clonal lines collected from a single Cape Cod marsh. To our knowledge, this is the first example of intraspecific variability in oxidative stress resistance in cnidarians or in any marine animal. As oxidative stress often accompanies heat stress in marine organisms, resistance to oxidative stress could strongly influence survival in warming oceans. For example, while elevated temperatures trigger bleaching in corals, oxidative stress is thought to be the proximal trigger of bleaching at the cellular level.

Thermal history and gape of individual Mytilus californianus correlate with oxidative damage and thermoprotective osmolytes

The ability of animals to cope with environmental stress depends – in part – on past experience, yet knowledge of the factors influencing an individual's physiology in nature remains underdeveloped. We used an individual monitoring system to record body temperature and valve gaping behavior of rocky intertidal zone mussels (Mytilus californianus). Thirty individuals were selected from two mussel beds (wave-exposed and wave-protected) that differ in thermal regime. Instrumented mussels were deployed at two intertidal heights (near the lower and upper edges of the mussel zone) and in a continuously submerged tidepool. Following a 23-day monitoring period, measures of oxidative damage to DNA and lipids, antioxidant capacities (catalase activity and peroxyl radical scavenging) and tissue contents of organic osmolytes were obtained from gill tissue of each individual. Univariate and multivariate analyses indicated that inter-individual variation in cumulative thermal stress is a predominant driver of physiological variation. Thermal history over the outplant period was positively correlated with oxidative DNA damage. Thermal history was also positively correlated with tissue contents of taurine, a thermoprotectant osmolyte, and with activity of the antioxidant enzyme catalase. Origin site differences, possibly indicative of developmental plasticity, were only significant for catalase activity. Gaping behavior was positively correlated with tissue contents of two osmolytes. Overall, these results are some of the first to clearly demonstrate relationships between inter-individual variation in recent experience in the field and inter-individual physiological variation, in this case within mussel beds. Such micro-scale, environmentally mediated physiological differences should be considered in attempts to forecast biological responses to a changing environment.

Multimodal in situ datalogging quantifies inter-individual variation in thermal experience and persistent origin effects on gaping behavior among intertidal mussels (Mytilus californianus)

In complex habitats, environmental variation over small spatial scales can equal or exceed larger-scale gradients. This small-scale variation may allow motile organisms to mitigate stressful conditions by choosing benign microhabitats, whereas sessile organisms may rely on other behaviors to cope with environmental stresses in these variable environments. We developed a monitoring system to track body temperature, valve gaping behavior, and posture of individual mussels (Mytilus californianus) in field conditions in the rocky intertidal zone. Neighboring mussels’ body temperatures varied by up to 14°C during low tides. Valve gaping during low tide and postural adjustments, which could theoretically lower body temperature, were not commonly observed. Rather, gaping behavior followed a tidal rhythm at a warm, high intertidal site; this rhythm shifted to a circadian period at a low intertidal site and for mussels continuously submerged in a tidepool. However, individuals within a site varied considerably in time spent gaping when submerged. This behavioral variation could be attributed in part to persistent effects of mussels’ developmental environment. Mussels originating from a wave-protected, warm site gaped more widely, and they remained open for longer periods during high tide than mussels from a wave-exposed, cool site. Variation in behavior was modulated further by recent wave heights and body temperatures during the preceding low tide. These large ranges in body temperatures and durations of valve closure events – which coincide with anaerobic metabolism – support the conclusion that individuals experience “homogeneous” aggregations such as mussel beds in dramatically different fashion, ultimately contributing to physiological variation among neighbors.