The flexibility afforded to genotypes in different environments by phenotypic plasticity is of interest to biologists studying thermal adaptation because of the thermal lability of many traits. Differences in thermal performance and reaction norms can provide insight into the evolution of thermal adaptation to explore broader questions such as species distributions and persistence under climate change. One approach is to study the effects of temperature on fitness, morphological and more recently gene expression traits in populations from different climatic origins. The diverse climatic conditions experienced by Drosophila melanogaster along the eastern Australian temperate-tropical gradient are ideal given the high degree of continuous trait differentiation, but reaction norm variation has not been well studied in this system. Here, we reared a tropical and temperate population from the ends of the gradient over six developmental temperatures and examined reaction norm variation for five quantitative traits including thermal performance for fecundity, and reaction norms for thermotolerance, body size, viability and 23 transcript-level traits. Despite genetic variation for some quantitative traits, we found no differentiation between the populations for fecundity thermal optima and breadth, and the reaction norms for the other traits were largely parallel, supporting previous work suggesting that thermal evolution occurs by changes in trait means rather than by reaction norm shifts. We examined reaction norm variation in our expanded thermal regime for a gene set shown to previously exhibit GxE for expression plasticity in east Australian flies, as well as key heat-shock genes. Although there were differences in curvature between the populations suggesting a higher degree of thermal plasticity in expression patterns than for the quantitative traits, we found little evidence to support a role for genetic variation in maintaining expression plasticity.

Thermal performance curves (TPCs) are intended to approximate the relationship between temperature and fitness, and are commonly integrated into species distributional models for understanding climate change responses. However, TPCs may vary across traits because selection and environmental sensitivity (plasticity) differ across traits or because the timing and duration of the temperature exposure, here termed time scale, may alter trait variation. Yet, the extent to which TPCs vary temporally and across traits is rarely considered in assessments of climate change responses. Using a common garden approach, we estimated TPCs for standard metabolic rate (SMR), and activity in Drosophila melanogaster at three test temperatures (16, 25 and 30°C), using flies from each of six developmental temperatures (16, 18, 20, 25, 28 and 30°C). We examined the effects of time scale of temperature exposure (minutes/ hours versus days/weeks) in altering TPC shape and position, and commonly used descriptors of the TPC: thermal optimum (T opt), thermal limits (T min and T max) and thermal breadth (T br). In addition, we collated previously published estimates of TPCs for fecundity and egg-to-adult viability in D. melanogaster. We found that the descriptors of the TPCs varied across traits (egg-to-adult viability, SMR, activity and fecundity), but variation in TPCs within these traits was small across studies when measured at the same time scales. The time scale at which traits were measured contributed to greater variation in TPCs than the observed variance across traits, although the relative importance of time scale differed depending on the trait (activity versus fecundity). Variation in the TPC across traits and time scales suggests that TPCs using single traits may not be an accurate predictor of fitness and thermal adaptation across environments.

The cellular stress response has long been the primary model for studying the molecular basis of thermal adaptation, yet the link between gene expression, RNA metabolism and physiological responses to thermal stress remains largely unexplored. We address this by comparing the transcriptional and physiological responses of three geographically distinct populations of Drosophila melanogaster from eastern Australia in response to, and recovery from, a severe heat stress with and without a prestress hardening treatment. We focus on starvin (stv), recently identified as an important thermally responsive gene. Intriguingly, stv encodes seven transcripts from alternative transcription sites and alternative splicing, yet appears to be rapidly heat inducible. First, we show genetic differences in upper thermal limits of the populations tested. We then demonstrate that the stv locus does not ubiquitously respond to thermal stress but is expressed as three distinct thermal and temporal RNA phenotypes (isoforms). The shorter transcript isoforms are rapidly upregulated under stress in all populations and show similar molecular signatures to heat-shock proteins. Multiple stress exposures seem to generate a reserve of pre-mRNAs, effectively 'priming' the cells for subsequent stress. Remarkably, we demonstrate a bypass in the splicing blockade in these isoforms, suggesting an essential role for these transcripts under heat stress. Temporal profiles for the weakly heat responsive stv isoform subset show opposing patterns in the two most divergent populations. Innate and induced transcriptome responses to hyperthermia are complex, and warrant moving beyond gene-level analyses.

No abstract