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.

Tracing the evolution of tissue inhibitor of metalloproteinases in Metazoa with the Pteria penguin genome.

Tissue inhibitors of metalloproteinase (TIMPs) play a pivotal role in regulating extracellular matrix (ECM) dynamics and have been extensively studied in vertebrates. However, understanding their evolution across invertebrate phyla is limited. Utilizing the high-quality Pteria penguin genome, we conducted phylogenomic orthology analyses across metazoans, revealing the emergence and distribution of the TIMP gene family. Our findings show that TIMP repertoires originated during eumetazoan radiation, experiencing independent duplication events in different clades, resulting in varied family sizes. Particularly, Pteriomorphia bivalves within Mollusca exhibited the most significant expansion and displayed the most diverse TIMP repertoires among metazoans. These expansions were attributed to multiple gene duplication events, potentially driven by the demands for functional diversification related to multiple adaptive traits, contributing to the adaptation of Pteriomorphia bivalves as stationary filter feeders. In this context, Pteriomorphia bivalves offer a promising model for studying invertebrate TIMP evolution.

Genome assembly and population genomic data of a pulmonate snail Ellobium chinense.

Ellobium chinense is an airbreathing, pulmonate gastropod species that inhabits saltmarshes in estuaries of the northwestern Pacific. Due to a rapid population decline and their unique ecological niche in estuarine ecosystems, this species has attracted special attention regarding their conservation and the genomic basis of adaptation to frequently changing environments. Here we report a draft genome assembly of E. chinense with a total size of 949.470 Mb and a scaffold N50 of 1.465 Mb. Comparative genomic analysis revealed that the GO terms enriched among four gastropod species are related to signal transduction involved in maintaining electrochemical gradients across the cell membrane. Population genomic analysis using the MSMC model for 14 re-sequenced individuals revealed a drastic decline in Korean and Japanese populations during the last glacial period, while the southern Chinese population retained a much larger effective population size (Ne). These contrasting demographic changes might be attributed to multiple environmental factors during the glacial-interglacial cycles. This study provides valuable genomic resources for understanding adaptation and historical demographic responses to climate change.

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.

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.

Evaluation of species thermal sensitivity with individual-based physiological performance.

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

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

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

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.

Thermal heterogeneity is an important factor for maintaining the genetic differentiation pattern of the pelagic barnacle Lepas anatifera in the northwest Pacific.

Macrobenthic invertebrates are ubiquitously distributed in the epipelagic zone of the open ocean. Yet, our understanding of their genetic structure patterns remains poorly understood. Investigating the genetic differentiation patterns of pelagic Lepas anatifera and clarifying the potential roles of temperature maintaining this pattern are crucial for our understanding of the distribution and biodiversity of pelagic macrobenthos. In the present study, mitochondrial cytochrome oxidase subunit I (mtDNA COI) from three South China Sea (SCS) populations and six Kuroshio Extension (KE) region populations of L. anatifera sampled from fixed buoys and genome-wide SNPs from a subset of populations (two SCS populations and four KE region populations) were sequenced and analyzed for investigating the genetic pattern of the pelagic barnacle. Water temperature was different among sampling sites; in other words, the water temperature decreased with latitude increases, and the water temperature on the surface was higher than in the subsurface. Our result showed that three lineages with clear genetic differentiation were found in different geographical locations and depths based on mtDNA COI, all SNPs, neutral SNPs, and outlier SNPs. Lineage 1 and lineage 2 were dominant in the subsurface populations and surface populations from the KE region, respectively. Lineage 3 was dominant in the SCS populations. Historical events during the Pliocene epoch shaped the differentiation of the three lineages, while, nowadays, temperature heterogeneity maintains the current genetic pattern of L. anatifera in the northwest Pacific. The subsurface populations were genetically isolated from the surface populations in the Kuroshio Extension (KE) region, implying small-scale vertical thermal heterogeneity was also an important factor maintaining the genetic differentiation pattern of the pelagic species.

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.

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

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

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.

Aerobic scope in fishes with different lifestyles and across habitats: Trade-offs among hypoxia tolerance, swimming performance and digestion.

Exercise and aerobic scope in fishes have attracted scientists' attention for several decades. While it has been suggested that aerobic scope may limit behavioral expression and tolerance to environmental stressors in fishes, the exact importance of aerobic scope in an ecological context remains poorly understood. In this review, we examine the ecological relevance of aerobic scope by reconsidering and reanalyzing the existing literature on Chinese freshwater fishes across a wide-range of habitats and lifestyles. The available evidence suggests that natural selection in fast-flowing aquatic habitats may favor species with a high aerobic scope and anaerobic capacity for locomotion, whereas in relatively slow-flowing habitats, hypoxia tolerance may be favored at the cost of reduced locomotor capacity. In addition, while physical activity can usually cause fishes from fast-flowing habitats to reach their aerobic metabolic ceiling (i.e., maximum metabolic rate), possibly due to selection pressure on locomotion, most species from slow-flowing habitats can only reach their metabolic ceiling during digestion, either alone or in combination with physical activity. Overall, we suggest that fish exhibit a continuum of metabolic types, from a 'visceral metabolic type' with a higher digestive performance to a 'locomotion metabolic type' which appears to have reduced capacity for digestion but enhanced locomotor performance. Generally, locomotor-type species can either satisfy the demands of their high swimming capacity with a high oxygen uptake capacity or sacrifice digestion while swimming. In contrast, most visceral-type species show a pronounced decrease in swimming performance while digesting, probably owing to conflicts within their aerobic scope. In conclusion, the ecological relevance of aerobic scope and the consequent effects on other physiological functions are closely related to habitat and the lifestyle of a given species. These results suggest that swimming performance, digestion and hypoxia tolerance might coevolve due to dependence on metabolic traits such as aerobic scope.

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 biogeographical break of the Yangtze River Estuary in recent decades. Data from neutral single nucleotide polymorphisms (SNPs) 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 for 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 biogeographical pattern of N. yoldii has formed with the occurrence of demographic, physiological and genetic changes, and emphasize the roles of adaptation of marginal populations during range expansion.

Two complete mitochondrial genomes of the barnacle Lepas anatifera Linnaeus, 1758 (Scalpellomorpha, Lepadidae) implying the possibility of cryptic speciation.

The barnacle Lepas anatifera Linnaeus, 1758 (Scalpellomorpha, Lepadidae) is a worldwide distributed species. For investigating its genetic diversity in the northwest Pacific, two complete mitochondrial genomes were determined and analyzed. The lengths of the two complete mitogenomes were 15,708 bp and 15,703 bp, respectively. Both of them contained typical 37 genes with an identical order to L. anserifera Linnaeus, 1767 and L. australis Darwin, 1851 mitogenome. Except for ND1 and ND2, 11 protein-coding genes (PCGs) started with an ATN initiation codon (ATA, ATG, ATC, and ATT). Twelve PCGs were terminated with TAA or TAG stop codon, whereas ND1 possessed an incomplete termination codon (T-). Phylogenetic analysis revealed that L. australis and L. anserifera clustered together, and then with L. anatifera. The distinct genetic distances (0.17) based on concatenated sequence of 13 PCGs between the two mitogenomes of L. anatifera suggest the existence of cryptic speciation. Additional samples from multiple localities should be collected and analyzed to deepen the understanding of cryptic diversity within the northwest Pacific.

Mapping the potential for offshore aquaculture of salmonids in the Yellow Sea.

Mariculture has been one of the fastest-growing global food production sectors over the past three decades. With the congestion of space and deterioration of the environment in coastal regions, offshore aquaculture has gained increasing attention. Atlantic salmon (Salmo salar) and rainbow trout (Oncorhynchus mykiss) are two important aquaculture species and contribute to 6.1% of world aquaculture production of finfish. In the present study, we established species distribution models (SDMs) to identify the potential areas for offshore aquaculture of these two cold-water fish species considering the mesoscale spatio-temporal thermal heterogeneity of the Yellow Sea. The values of the area under the curve (AUC) and the true skill statistic (TSS) showed good model performance. The suitability index (SI), which was used in this study to quantitatively assess potential offshore aquaculture sites, was highly dynamic at the surface water layer. However, high SI values occurred throughout the year at deeper water layers. The potential aquaculture areas for S. salar and O. mykiss in the Yellow Sea were estimated as 52,270 ± 3275 (95% confidence interval, CI) and 146,831 ± 15,023 km2, respectively. Our results highlighted the use of SDMs in identifying potential aquaculture areas based on environmental variables. Considering the thermal heterogeneity of the environment, this study suggested that offshore aquaculture for Atlantic salmon and rainbow trout was feasible in the Yellow Sea by adopting new technologies (e.g., sinking cages into deep water) to avoid damage from high temperatures in summer. Supplementary Information: The online version contains supplementary material available at 10.1007/s42995-022-00141-2.

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.

Thermal adaptation of mRNA secondary structure: stability versus lability.

Macromolecular function commonly involves rapidly reversible alterations in three-dimensional structure (conformation). To allow these essential conformational changes, macromolecules must possess higher order structures that are appropriately balanced between rigidity and flexibility. Because of the low stabilization free energies (marginal stabilities) of macromolecule conformations, temperature changes have strong effects on conformation and, thereby, on function. As is well known for proteins, during evolution, temperature-adaptive changes in sequence foster retention of optimal marginal stability at a species' normal physiological temperatures. Here, we extend this type of analysis to messenger RNAs (mRNAs), a class of macromolecules for which the stability-lability balance has not been elucidated. We employ in silico methods to determine secondary structures and estimate changes in free energy of folding (ΔGfold) for 25 orthologous mRNAs that encode the enzyme cytosolic malate dehydrogenase in marine mollusks with adaptation temperatures spanning an almost 60 °C range. The change in free energy that occurs during formation of the ensemble of mRNA secondary structures is significantly correlated with adaptation temperature: ΔGfold values are all negative and their absolute values increase with adaptation temperature. A principal mechanism underlying these adaptations is a significant increase in synonymous guanine + cytosine substitutions with increasing temperature. These findings open up an avenue of exploration in molecular evolution and raise interesting questions about the interaction between temperature-adaptive changes in mRNA sequence and in the proteins they encode.