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