Thermal reaction norms of growth in Atlantic and Pacific silverside fishes

How organisms may adapt to rising global temperatures is uncertain, but concepts can emerge from studying adaptive physiological trait variations across existing spatial climate gradients. Many ectotherms, particularly fish, have evolved increasing genetic growth capacities with latitude (i.e. countergradient variation (CnGV) in growth), which are thought to be an adaptation primarily to strong gradients in seasonality. In contrast, evolutionary responses to gradients in mean temperature are often assumed to involve an alternative mode, 'thermal adaptation'. We measured thermal growth reaction norms in Pacific silverside populations (Atherinops affinis) occurring across a weak latitudinal temperature gradient with invariant seasonality along the North American Pacific coast. Instead of thermal adaptation, we found novel evidence for CnGV in growth, suggesting that CnGV is a ubiquitous mode of reaction-norm evolution in ectotherms even in response to weak spatial and, by inference, temporal climate gradients. A novel, large-scale comparison between ecologically equivalent Pacific versus Atlantic silversides (Menidia menidia) revealed how closely growth CnGV patterns reflect their respective climate gradients. While steep growth reaction norms and increasing growth plasticity with latitude in M. menidia mimicked the strong, highly seasonal Atlantic coastal gradient, shallow reaction norms and much smaller, latitude-independent growth plasticity in A. affinis resembled the weak Pacific latitudinal temperature gradient.

Theoretical model simulations for variations in thermal reaction norms

Thermal reaction norms for growth rates of six Emiliania huxleyi isolates originating from the central Atlantic (Azores, Portugal) and five isolates from the coastal North Atlantic (Bergen, Norway) were assessed. We used the template mode of variation model to decompose variations in growth rates into modes of biological interest: vertical shift, horizontal shift, and generalist-specialist variation. In line with the actual habitat conditions, isolates from Bergen (Bergen population) grew well at lower temperatures, and isolates from the Azores (Azores population) performed better at higher temperatures. The optimum growth temperature of the Azores population was significantly higher than that of the Bergen population. Neutral genetic differentiation was found between populations by microsatellite analysis. These findings indicate that E. huxleyi populations are adapted to local temperature regimes. Next to between-population variation, we also found variation within populations. Genotype-by-environment interactions resulted in the most pronounced phenotypic differences when isolates were exposed to temperatures outside the range they naturally encounter. Variation in thermal reaction norms between and within populations emphasizes the importance of using more than one isolate when studying the consequences of global change on marine phytoplankton. Phenotypic plasticity and standing genetic variation will be important in determining the potential of natural E. huxleyi populations to cope with global climate change.