Genetic constraints for thermal coadaptation in Drosophila subobscura

Background Behaviour has been traditionally viewed as a driver of subsequent evolution because behavioural adjustments expose organisms to novel environments, which may result in a correlated evolution on other traits. In Drosophila subobscura, thermal preference and heat tolerance are linked to chromosomal inversion polymorphisms that show parallel latitudinal clines worldwide, such that "cold-climate" ("warm-climate") chromosome arrangements collectively favour a coherent response to colder (warmer) settings as flies carrying them prefer colder (warmer) conditions and have lower (higher) knock out temperatures. Yet, it is not clear whether a genetic correlation between thermal preference and heat tolerance can partially underlie such response. Results We have analyzed the genetic basis of thermal preference and heat tolerance using isochromosomal lines in D. subobscura. Chromosome arrangements on the O chromosome were known to have a biometrical effect on thermal preference in a laboratory temperature gradient, and also harbour several genes involved in the heat shock response; in particular, the genes Hsp68 and Hsp70. Our results corroborate that arrangements on chromosome O affect adult thermal preference in a laboratory temperature gradient, with cold-climate Ost carriers displaying a lower thermal preference than their warm-climate O3+4 and O3+4+8 counterparts. However, these chromosome arrangements did not have any effect on adult heat tolerance and, hence, we putatively discard a genetic covariance between both traits arising from linkage disequilibrium between genes affecting thermal preference and candidate genes for heat shock resistance. Nonetheless, a possible association of juvenile thermal preference and heat resistance warrants further analysis. Conclusions Thermal preference and heat tolerance in the isochromosomal lines of D. subobscura appear to be genetically independent, which might potentially prevent a coherent response of behaviour and physiology (i.e., coadaptation) to thermal selection. If this pattern is general to all chromosomes, then any correlation between thermal preference and heat resistance across latitudinal gradients would likely reflect a pattern of correlated selection rather than genetic correlation.


Background
Ectotherms exhibit a suite of behavioural and physiological strategies to cope with spatiotemporal variation in ambient temperature [1]. For instance, behavioural adjustments (e.g. modifying daily activity patterns and selecting favourable microclimates; [2]) can buffer the impact of sub-optimal temperatures, and are the main means of thermoregulation in small insects [3][4][5]. Although such adjustments can enable ectotherms to maintain relatively constant body temperatures (T b ) at different seasons and/or latitudes [2,6], the observation of cyclical seasonal changes in genetic markers putatively related to thermal adaptation [7,8] and the clinal variation in thermal stress tolerance in some Drosophila species [9][10][11] suggest that behavioural thermoregulation may be insufficient to fully compensate shifts in environmental temperature [12].
If behavioural thermoregulation is not fully compensatory and climate variation influences the actual T b and physiological performance of organisms distributed over broad latitudinal ranges (i.e., performance falls below its optimum during cooling and warming), then temperature is more than just a key environmental factor that affects development, growth, and survival of individuals [13,14]: it is likely the main selective agent that drives -directly or indirectly -the evolution of clinal patterns in genetic, phenotypic, and life history traits. Furthermore, the divergence of thermal optima in the different subpopulations according to the T b experienced by the organism is expected to bolster a covariance between behavioural shifts (thermal preference) and performance [6,15]. This is related to the idea of "coadaptation" [16], where natural selection is supposed to favour the harmonious adjustment among the suite of (co-)evolving traits [7,17]. Parallel clines on different continents or along independent temperature gradients can thus offer an invaluable opportunity to study thermal coadaptation since the role of temperature in driving those clines is quite compelling.
Some widespread latitudinal clines in Drosophila also provide an additional advantage for studies of thermal coadaptation: there is a relatively well-known historical record following the invasion of a new geographical region (e.g. [18,19]). Perhaps the best example is that of Drosophila subobscura, a native Palaearctic species that invaded the Americas about 30 years ago, and spread rapidly on both South and North America. Clinal patterns for phenotypic traits and genetic polymorphisms emerged very rapidly during these two independent colonization events [20][21][22]. For instance, North American populations soon evolved decreased desiccation resistance with increasing latitude as expected, which matches the pattern found in Old World populations and suggests that strong selection for thermal-related traits along latitudinal gradients is taking place. On the other hand, in South America this trait shows the opposite pattern: higher desiccation tolerance is observed in colder areas [23]. Contrasting outcomes were also observed for other clinally varying traits -wing cell size and cell number [24], and wing shape [20,25] -where the role of temperature remains elusive, which apparently suggests that selective pressures vary in the different clines. An alternative explanation, however, is that evolution can sometimes be constrained by antagonistic genetic correlations (i.e., genetic correlations among traits that are not in accord with the direction of selection [26,27]) arising from linkage disequilibrium between alleles at different loci, and patterns of linkage disequilibrium can vary among populations or seasons [28,29]. In this context, we now know that contrasting wing shape clines in D. subobscura came out as a correlated response of the world-wide parallel inversion clines [21] because inversion-shape relationships in native and colonizing populations are opposite (presumably due to the different associations between inversions and particular alleles which influence the trait), probably as a result of the bottleneck effect that occurred during the colonization of America [30]. Besides, different patterns of linkage disequilibrium could result from variability in migration rates between genetically differentiated populations in the various latitudinal clines [31]. In summary, conflicting outcomes between old and rapidly evolving new clines should probably not be viewed as a nuisance, but as reminder that an appropriate knowledge of the underlying genetic architecture is required to further understand why (or why not) these inconsistencies arise. More specifically, if behaviour "drives" the subsequent parallel evolution in morphology and physiology as predicted ( [6]; but see [32]), it is essential to analyze the genetic basis of thermal preference and temperaturerelated traits to see whether or not thermal coadaptation can happen along a cline.
We have recently undertaken a within-population large-scale study to analyze the association between chromosomal inversion polymorphisms that show parallel latitudinal clines in native and colonizing populations of D. subobscura, with the thermal preferences (T p : the preferred body temperature in a laboratory thermal gradient, which we expect to correlate with the thermal optimum for performance; [33]) and knock out temperatures (T ko : the temperature required to knock out a fly in a waterbath) of their carriers [34]. The main results can be summarized as follows: (i) flies carrying "cold-adapted" or "cold-climate" chromosome arrangements (i.e., those chromosome arrangements in all five major acrocentric chromosomes that show a negative correlation coefficient with maximum temperatures along the cline, or a positive correlation coefficient with latitude in Palaearctic populations; [35,36]) prefer a lower T p and had a lower T ko , in accordance with the natural patterns; (ii) different chromosomes were responsible for the bulk of the genetic variation in T p (chromosomes A and O) and T ko (chromosome E); and (iii) T p and T ko were phenotypically uncorrelated, which agrees with the observation that different independently segregating chromosomes were mainly responsible for the corresponding associations. Taken at a face value, behavioural thermoregulation and performance were indeed "coadapted" in the sense that cold-climate (warm-climate) chromosome arrangements collectively favour a coherent response to colder (warmer) environments, but this was not due to a genetic covariance of behaviour and physiology. There were, however, two potential limitations in the study. First, each individual fly was scored for only one chromosome of its diploid set and, hence, dominance effects (if any) where hidden in the analysis. Second, both intraand interchromosomal contributions were mixed because the assayed flies had the genetic background from the sampled wild population. Although it might be argued that this protocol is somehow closer to what happens in nature, these uncontrolled factors might have precluded a better characterization of the underlying genetic effects. Accordingly, although the amount of genetic variation on T p and T ko explained by the combined effect of all chromosomes carrying at least one cold-climate gene arrangement was statistically significant, it only accounted for 1% of the total phenotypic variation [34].
Here we examine if T p and T ko are genetically correlated and might evolve in a coherent fashion in response to selection; i.e., whether behaviour and physiology are coadapted at the genetic level. We take advantage of the fact that the polymorphic inversions on chromosome O appear to be associated with behavioural thermoregulation in D. subobscura [34], and that this is the only chromosome that can be used to measure the expression of associated traits in replicated inbred and outbred genotypes. Namely, chromosome O is the only one for which a balancer stock (Va/Ba: Varicose/Bare; [37]) is available (a balancer is a specially constructed chromosome that carries a dominant morphological marker that is homozygous lethal and multiple inversions to suppress recombination). This is the longest chromosome in D. subobscura (190 cM which correspond to approximately 31 Mb [38]), and is homologous to arm 3R in D. melanogaster [39,40]. Some chromosome arrangements (O st and O 3+4 ) show conspicuous northwest-southwest latitudinal clines in Palaearctic populations ( Figure 1a). Chromosome O harbours several genes involved in the heat shock response [41]; in particular, gene Hsp68 (located in section O(89A) [42,43] and relatively close to the proximal breakpoint of inversion O 8 [44]), and gene Hsp70 (located in section O(94A) [42,43] and included inside the warm-climate chromosomal arrangement O 3+4 , and close to the distal breakpoint of inversion O 8 [44]) ( Figure 1b). Hsp70 appears to be the primary protein involved in thermotolerance in D. melanogaster [45] though apparently not in other Drosophila species [46] -, and Hsp70 allele frequencies show latitudinal clines and change in response to thermal evolution in the laboratory [47]. In addition, correlated responses to selection for knock down resistance at 39°C have also been found for Hsp68 in D. melanogaster [48].
Previous work also showed that D. subobscura flies carrying O chromosomes derived from replicated thermal lines [49,50] that had evolved in the laboratory at warm temperatures (22°C) had a higher total net fitness than its cold-adapted (13°C) counterparts; that is, a significant shift in thermal optima was observed [51]. All in all, it seems that there is indeed room for the coevolution of behaviour and physiological tolerance in D. subobscura. However, we show here that thermal preference and heat tolerance appear to be genetically independent. Therefore, any latitudinal correlation between both traits would likely reflect a pattern of correlated selection across populations rather than within-population genetic correlations.

Experimental settings
In south-western European populations, the most frequent chromosome arrangements for chromosome O are O st , O 3+4 , O 3+4+7 , and O 3+4+8 [52]. The first two arrangements show a clear contrasting clinal pattern in original Palaearctic populations, with O st increasing and O 3+4 decreasing in frequency with increasing latitude [35,56] (Figure 1a). Arrangement O 3+4+8 is also interesting because in historical times it was mainly restricted to the Mediterranean region, being the most abundant chromosomal arrangement in northern Africa [53]. However, in the last decades its distribution has changed dramatically and recent surveys revealed frequencies as high as 22.6% in Groningen, Netherlands, where it was previously absent [22,54]. Six independent isochromosomal lines for each of these three arrangements (i.e., has been detected among these arrangements [55]. In other words, there are compelling reasons to think that the chromosome arrangements used in this work are genetically differentiated for Hsp70, and probably also for Hsp68 since inversion effects can extend as far as 1000 kilobases outside from breakpoints [56,57]. Following Santos et al. [58] the experimental flies were obtained from 54 crosses, which will be referred to as ,..., with 18 crosses in total), or as outbred including both structural ,..., with 18 cyclically permuted reciprocal crosses in total) and het- ,..., ; j ≠ k; with 18 reciprocal crosses in total). Two developmental temperatures were used in the experiment to study potentially important effects of phenotypic plasticity: 18°C and 22°C. The reason for this was the huge difference (about 7°C-8°C) between our previous estimate of T p (pooled average 16.6°C; [34]) in D. subobscura flies raised at 18°C, and that obtained by Huey and Pascual (23.7°C; [12]) where flies were raised at 22°C. Even though the flies assayed came from different sourcessouth-western Europe in Rego et al. [34], and North America in Huey and Pascual [12] -, which could account for the observed difference because thermal responses can vary between populations [59], it remains to be seen whether developmental plasticity can affect estimates of thermal preference and heat tolerance.

Association between thermal preference and knock out temperature
The phenotypic correlation between T p and T ko was assessed from their partial correlation coefficient, holding constant the variables developmental temperature, sex, plate hour, and water bath (see Methods). In no case were the partial correlations statistically significant: , t = 0.004, df = 1054, P = 0.997) and, hence, they will be analyzed separately in what follows.

Consanguinity and developmental effects a) Thermal preference
Inbreeding and developmental temperature effects on T p were simultaneously analyzed by contrasting isogenic vs. outbred homokaryotypic flies reared at both experimental temperatures ( Figure 2). The factorial analysis of covariance (ANCOVA) only detected statistically significant differences for karyotypes, karyotype × inbreeding interaction, and karyotype × developmental temperature interaction effects (Table 1) sexes (females: 19.0°C ± 3.6°C; males: 18.5°C ± 3.6°C), although in this last case the effect was marginally nonsignificant (P = 0.053). Permutation tests (see Methods) corroborated that the three assayed karyotypes differ in T p (P = 0.001).
Scheffé post hoc tests using the mean square of the nested "cross" effect as the error term showed that the thermal preference Males Figure 2 Inbreeding and temperature effects on thermal preference. Homokaryotipic averages for T p (in°C with 95% confidence intervals) in inbred (left panels) and outbred (right panels) crosses according to sex and developmental temperature.
them. The difference is consistent for both isogenic and outbred flies ( Figure 2). From the present data we can conclude that the preferred temperature ranges or "set point" (T set ) ranges (central 50% of preferred body temperatures; [60] Gene arrangement effects in the outbred lines a) Thermal preference The genetic and environmental (developmental temperature) contributions of chromosome O to T p (and T ko ; below) was assessed from the outbred crosses including all possible karyotypes. Outbred crosses are obviously more relevant to the real situation because inbred genotypes are homozygous for deleterious alleles, and also for alleles that might display heterozygote advantage in the original outbred population. The only statistically significant effects detected by the ANCOVA model (Table 3) were those arising from genetic differences among karyotypes (permutation tests corroborated that the three assayed karyotypes differ in T p ; P = 0.0018) and sexes, with females having a higher T p (mean ± SD: 18.7°C ± 3.6°C) than males (18.0°C ± 3.6°C). As above, average T p was slightly lower for flies reared at 18°C (18.1°C ± 4.0°C) than at 22°C (18.6°C ± 3.2°C), but the difference was marginally non-significant (P = 0.069).

b) Knock out temperature
The ANCOVA for T ko (Table 4) did not detect any difference among karyotypes, in accordance with the previous findings for the inbred crosses. Similarly, the main differences arose between developmental temperature (flies reared at 18°C: 33.6°C ± 1.9°C; flies reared at 22°C: 32.8°C ± 2.3°C) and sex (females: 33.7°C ± 1.8°C; males: 32.7°C ± 2.4°C).
The genetic correlation between T p and T ko after pooling O 3+4 and O 3+4+8 was r p = -0.130 (P = 0.917). Again, the conclusion is that these two traits are uncorrelated. Figure 4 plots the genotypic values in the additive-dominance scales for T p and T ko , together with their statistical significance obtained from the appropriate contrasts (Table 3, 4).

c) Average effects on thermal preference
Our experiment only provides an estimation of the gene (chromosome O) action on T p and does not allow inferences to the base population. It is possible, however, to obtain estimates of the average effects, or "statistically additive effects", by taking into account the gene action and allelic (chromosome arrangement) frequencies in the natural populations [61]. Assuming that the chromosome arrangement effects are roughly the same along the cline (for a measure of climatic temperatures along the Palaearctic cline see Figure 1 in [62]), Table 5  Males Figure 3 Inbreeding and temperature effects on knock out temperature. Homokaryotipic averages for T ko (in°C with 95% confidence intervals) in inbred (left panels) and outbred (right panels) crosses according to sex and developmental temperature.  [52,54]. The interpretation is that flies inheriting a O st chromosome will choose a temperature ranging from around 0.31°C -0.45°C below the average temperature chosen by the population (conversely, flies carrying warm-climate chromosome arrangements will choose a temperature ranging from around 0.03°C -0.52°C above the average). Combined with our previous results with chromosome A (which is the sex chromosome and additive values can be estimated using males' T p ; [34]), where gene arrangement A st exhibits a similar latitudinal pattern than O st and flies carrying A st also display a laboratory thermal preference towards colder temperature, the conclusion is that flies inheriting simultaneously A st and O st will choose temperatures ranging from approximately 0.5°C -1.0°C below the average (these estimates assume perfect additivity).

Discussion
The present results with isogenic lines and their crosses corroborate and extend our previous work with wild flies from south-western Europe [34]. They confirm that arrangements on chromosome O have a biometrical effect on thermal preference in a laboratory temperature gradient, with cold-climate O st carriers displaying a lower T p than their warm-climate O 3+4 and O 3+4+8 counterparts. In addition, T p and T ko were again found to be uncorrelated, and we can now discard a potential genetic covariance between both traits arising from linkage disequilibrium between genes affecting thermal preference and candidate genes for heat shock resistance (i.e., Hsp68 and Hsp70; [42,43]) located inside of, or close to, the chromosome regions covered by the inversions analyzed here (see Background). In other words, we conclude that variation on O chromosome arrangements does not have any effect on knock out temperature (but see below). Note, however, that this does not imply that genes on chromosome O have no effect on   T ko (actually, statistically significant differences were detected among crosses within karyotypes; Table 4); it simply indicates that any allelic variation of putative genes influencing this trait is not in linkage disequilibrium with inversions on this chromosome.
The new findings were: (i) a lack of inbreeding depression for both T p and T ko ; (ii) a lack of phenotypic plasticity for T p according to the temperature at which the flies were raised (18°C and 22°C); and (iii) a substantial effect of developmental temperature on T ko .   thermal limits is the heat shock response, and in D. melanogaster the major inducible heat shock protein Hsp70 appears to be the primary protein involved in thermotolerance [45,63]. Recent work, however, questions the pervasive role of Hsp70 in the mediation of the heat stress response and suggests that it may be life-stage specific, being important in larvae but not in adults [64]. Our results are apparently consistent with the lack of association between Hsp70 and adult heat resistance (but see further discussion below), although also raise a caveat to the conclusion that there is no covariance between T p and T ko . Thus, it could be the case that Hsp70 variation across karyotypes is associated with juvenile tolerance to heat stress, an important trait in Drosophila particularly in summer when larval feeding patches can become lethally hot [65]. This possibility warrants further analysis.
An important concern here is that Hsp70 production might not be inducible in the dynamic experimental protocol we used to estimate upper thermal tolerance, where temperature increased 0.1°C min -1 . One apparently compelling reason for this is that the estimated maximum thermal limits that D. melanogaster can tolerate decrease from approximately 39.9°C with heating rate 0.5°C min -1 to 38.7°C with heating rate 0.1°C min -1 [66], a puzzling result because slower heating rates should allow individuals to acclimatize to new temperatures and also because slow heating rates pre-exposes individuals to non-lethal high temperatures ("hardening"), which increases heat shock resistance [10]. We have recently discussed why these conflicting outcomes arise, and suggest that the contribution of other stressors (e.g. higher desiccation in long thermal tolerance assays associated with slow warming rates) can potentially overshadow thermal acclimation effects in dynamic assays with varying heating rates [67]. In other words, we challenge the idea that induced thermotolerance does not occur in dynamic assays with slow heating rates. At this stage this is just speculative because Hsp70 production was not measured in our flies, but the problem is important because Drosophila adults are likely to experience slow heating rates in nature of 0.06 -0.1°C min -1 [66,68] and further empirical studies are required to explain the apparently inconsistent findings.
The pooled average T p here was (mean ± SD) 18.4°C ± 3.6°C (T set : 15.4°C -21.2°C; these figures include only outbred lines) and about the same at both rearing temperatures. The difference with our previous estimate for wild-flies from Adraga (16.6°C, T set :12.4°C -20.4°C; [34]) does not seem to be overreached, and could be partially explained by the fact that the present flies were genetically homogeneous for all chromosomes from the ch-cu marker strain but chromosome O (recall that the sex chromosome A also had a significant effect on T p ; [34]). This strain has a long history of maintenance at 18°C in the laboratory. In any case, our estimates remain substantially lower than that from Huey and Pascual (23.7°C, T set : 21.2°C -25.9°C; [12]), and the difference cannot be accounted by flies' rearing temperature. No reasonable explanation for the discrepancy can be offered at this moment, but the additional result that developmental temperature substantially affected T ko makes us confidently conclude that our estimates are indeed closer to the actual T p of the species. Flies reared at 22°C showed lower heat resistance than their counterparts reared at 18°C (32.8°C vs. 33.6°C; outbred lines), which could be a consequence of their smaller size due to the inverse relationship between body size and developmental temperature [69,70]. However, resistance to heat does not seem to be associated with body size [71] -we have also analyzed the association between T ko and wing size from our previous experiment where both traits were recorded [30,34] and found no relationship whatsoever (results not shown). Most likely, 22°C was a suboptimal and potentially stressful temperature for our flies, making them to be weaker and less resistant to the heat shock. Note, however, that this conclusion might not be extrapolated to wild flies that harbour higher levels of genetic variability than our chromosomal lines.
To interpret the interplay between thermal preference and heat stress resistance, an understanding of the environmental temperatures experienced by D. subobscura along climatic gradients is required. As far as we are aware, the only data available on T b for active flies along a latitudinal gradient (spanning 12°) come from recent work by Huey and Pascual [12] in western North America. They found that mean T b varies by as much as 21°C (from 8°C to 29°C), and that the temporal activity of flies during the day did not match predictions from optimal temperature regulation or desiccation avoidance. Temperatures of maximum activity in summer ( Figure 2 in [12]) -when wild flies are smaller probably due to their higher developmental temperatures and/or crowding conditions [72]; and crowding is known to affect adult thermal stress resistance in Drosophila [73] -are dangerously close to the T ko obtained here for the outbred flies raised at 22°C. This suggests that active D. subobscura flies can experience extreme conditions in the wild, and one would expect flies' activity to be correlated with heat resistance under these conditions if behaviour and physiology were coadapted. Some evidence indicates that diurnal activity patterns in summer can vary according to inversion polymorphism, and chromosome arrangements on the O chromosome seem to behave as expected from our data: O st is more frequent towards the evening while chromosomes carrying gene arrangement O 3+4 are most frequently sampled at midday [74]. This behavioural thermoregulation, however, would not confer less susceptibility to high temperatures because the genetic basis of both traits does not seem to allow for the building up of "coadaptation". It is well known from basic population genetics theory that genetic covariance between traits can arise when alleles at different loci are associated (linkage disequilibrium), and this critically depends on relatively low recombination rates [75]. The lack of association between T p and T ko in D. subobscura is fully consistent with their genetic basis as independently segregating chromosomes are involved [34]. Yet, a correlation between these traits can be expected at the interpopulational level due to patterns of correlated selection (rather than genetic correlations) across a latitudinal gradient because of the congruent latitudinal clinal variation for chromosome arrangements on the E (which influences T ko [34]), and on the A and O chromosomes (which influence T p [ [34], this work]).
We now speculate that the mismatch between T p and T ko could apparently generate an interesting dynamics in the population frequencies of different chromosome arrangements on chromosome O. Suppose the daily activity of flies in the warmest months of the year follows the previously described pattern; that is, flies carrying gene arrangement O 3+4 are more active at midday and, therefore, have a higher risk of a heat shock than O st and are selected against. On the other hand, assuming T p corresponds closely with temperatures that maximize fitness O 3+4 flies likely enjoy a fitness advantage in summer. The net effect would be a compromise between "behaviour unresponsiveness" and general performance, which means that chromosome arrangements on chromosome O may or may not cycle seasonally according to average environmental temperature (i.e., O 3+4 could be expected to increase in frequency in summer and decrease in winter if general performance is what matters). Interestingly, both patterns have been detected: consistent seasonal cycling at a north-western population in Spain [8] and apparently no seasonal variation at a north-eastern population also in Spain [76]. The point here is that parallel seasonal changes should also be detected for chromosome A since it also affects T p [34]. In accordance with this prediction, no seasonal cycling was detected for chromosome A in the northeastern population, but unfortunately no information is available for the other population because chromosome O was the only chromosome scored. It would be very interesting to see what happens for chromosome A in the cycling population.

Conclusions
For ectotherms facing spatiotemporal variation in environmental temperature theory predicts that a coevolution between thermal preference and physiological performance can occur [1]. In the widespread species D. subobscura behavioural thermoregulation and heat tolerance are "coadapted" in the sense that flies carrying cold-climate (warm-climate) chromosome arrangements tend to choose colder (warmer) temperatures and have lower (higher) heat stress tolerance [34]. We have analyzed the genetic basis of these thermal traits using isochromosomal lines for the O chromosome. This chromosome was known to affect thermal preference [34], and also harbours several genes involved in the heat shock response (Hsp68 and Hsp70) [42,43]. These genes are located inside of, or close to, the chromosome regions covered by inversions that show conspicuous northwest-southwest latitudinal clines in Palaearctic populations, as well as seasonal fluctuations that are in agreement with the latitudinal patterns [22]. Our results corroborate that arrangements on chromosome O affect adult thermal preference: flies inheriting the cold-climate O st chromosome are predicted to choose a temperature around 0.31°C -0.45°C below the average temperature chosen by the population and, conversely, flies inheriting the warm-climate O 3+4 and O 3+4+8 chromosomes are expected to choose a temperature ranging from around 0.03°C -0.52°C above the average. However, these chromosome arrangements did not have any differential effect on adult heat tolerance. We conclude that thermal preference and heat tolerance in D. subobscura appear to be genetically independent and, therefore, any latitudinal correlation between both traits would likely reflect a pattern of correlated selection across populations rather than within-population genetic correlations.

Methods
Origin of flies and experimental procedures D. subobscura wild flies were collected near Barcelona (41°43'N, 2°13'E) in October 2007. More than 200 isofemale lines were derived and used to obtain isochromosomal lines for the O chromosome in an otherwise homogeneous genetic background following standard protocols [77]. Briefly, one offspring male from each isofemale line was crossed to three or four virgin females from the ch-cu marker strain, which is homozygous for the morphological recessive markers on the O chromosome cherry eyes (ch) and curled wings (cu) and the chromosomal arrangement O 3+4 . A single wild-type male from each cross was repeatedly backcrossed to three or four ch-cu females for at least five generations in order to homogenize the genetic background, and the chromosomal arrangement carried by the wild chromosome was identified after the second backcross. To derive the isochromosomal lines, wild-type males from each line were crossed with the Va/Ba balancer stock [78], which has the same genetic background as the ch-cu strain. Once obtained, the isochromosomal lines were genotyped for 13 microsatellite loci located on the O chromosome to check that no recombination events occurred during the different crosses. The 18 independent isochromosomal lines used in this study (see Experimental settings) were found to be homozygous for all the loci. The lines were kept at 18°C (12:12 light/ dark cycle) in 130-mL bottles with low adult density (around 20 pairs/bottle) to standardize the rearing conditions before egg collections.
To obtain the experimental flies, all 54 crosses (inbred and outbred) were performed at 18°C by mating 4 daysold virgin males and females from the corresponding isochromosomal lines. After six days the males were discarded and the females (an equal number from each reciprocal cross in the outbred combinations) were transferred to egg-laying chambers containing fresh food and charcoal colouring. Eggs were placed in vials (45 eggs/vial containing 6 mL of food) at two rearing temperatures: 18°C and 22°C. Non-anaesthetized emerging flies were stored in bottles at low adult density and used to evaluate laboratory thermal preference (T p ) and knock out temperature (T ko ) for each cross (see below). All fly handling was done at room temperature using CO 2 anaesthesia only to sort virgin flies and to place females in the egg-laying chambers.
Thermal preference behaviour in a laboratory gradient and heat resistance Laboratory T p was measured as previously described [34]. Briefly, adult flies (about 7 days old) were individually placed in separate lanes on an aluminium base plate where a thermal gradient with temperatures ranging from 11°C to 29°C was generated. Adults were given approximately 1 h to adjust, and afterwards their positions were recorded four times every 10 min. We used the median of the four measurements to estimate T p of each fly. Measurements were performed in a room with constant temperature (22°C -23°C), and the flies were assayed under white light illumination. This protocol renders a repeatable assessment of flies' thermal preferences [34]. After the thermal preference assay, each fly was gently removed from the lane and individually placed in a vial with fresh food for the subsequent assay of heat stress tolerance.
One day after measurements of thermal preference flies were assayed for heat resistance also as previously described [34]. Adults were individually placed in sealed empty vials and immersed in water-baths at T min = 24°C . Every 10 min individuals were scored for mobility (fly active or knocked out) and the temperature of the water was increased by ΔT = +1°C. The procedure was repeated until the water-baths reached T max , defined as the temperature when the last active fly was knocked out (T max = 38°C was the upper limit in the assays; median T max = 33°C). For each fly T ko was estimated as the temperature taken to knock it out (defined as the onset of muscle spasms; [79]).

Statistical methods
The experimental setup was devised to assay one male and one female from each cross and temperature per day (five blocks) for both T p and T ko , amounting to 1,080 flies in total. Some mishaps (e.g. individuals flew away or just died during the assays) were, however, unavoidable and the final data set contains a few more than or a few less than 10 flies in several crosses (the harmonic means of flies per cross and temperature were: T p assay, 5.04 females and 4.80 males; T ko assay, 4.89 females and 4.37 males). Statistical analysis with and without block design qualitatively yielded the same results. Therefore, to simplify matters blocks were not considered in the linear models below.

a) Consanguinity and temperature effects
Inbreeding and temperature effects were simultaneously analyzed by contrasting isogenic vs. outbred homokaryotypic flies reared at both developmental temperatures. The linear model used was: where μ is the overall grand mean, i is the fixed effect of the karyotype (i =1, 2, 3), C j(i) is the random effect of the jth cross (j = 1, 2, ..., 6) within karyotype i, τ k is the fixed effect of the developmental temperature (18°C or 22°C), ι l is the fixed effect of inbreeding (isogenic or outbred homokaryotypic flies), ς m is the fixed effect of sex, and ε ijklmn is the residual error associated with the thermal preference (T p ) of the nth fly from the mth sex with the ith karyotype from the jth cross that was derived from the ιth group of crosses and assayed at the kth temperature. The covariate plate-hour was also introduced in the model to control for differences in circadian activity since several trials were conducted during each day. A similar linear model was used for knock out temperature, also introducing water-bath as a covariate since T ko was assessed in different water-baths.
Notice that for the main effect "karyotype" the linear model (1) can be conveniently reduced to the following two-level nested ANOVA model: where the sum of squares for the error term e ijk is simply the sum of the sum of squares for the remainder terms in (1). The usefulness of this model reduction is to efficiently perform randomization tests to test the null hypothesis about karyotype effects in a randomized (i.e., random assignment) experiment [80]. Permutation tests are far less sensitive to the presence of outliers than parametric tests. The null hypothesis of no karyotype effect was tested here after performing random permutations among replicate and selection temperature for the among selection temperature F-statistics. Each test used 10,000 random permutations.

b) Karyotype variation
To asses the effect of O chromosome karyotypes on T p and T ko we have focused in the outbred crosses, including both structural homo-and heterokaryotypes. The linear model used was similar to (1) including the fixed effect of karyotype ( i ; i = 1, 2, ..., 6), the random effect of cross within karyotypes (C j(i) ; j = 1, 2, ..., 6), the fixed effect of developmental temperature, and the fixed effect of sex. The covariate plate-hour was also introduced in the model. As above, a similar linear model was used for knock out temperature, also introducing water-bath as a covariate.
In the original Palaearctic populations chromosome arrangements O 3+4 and O 3+4+8 have a higher frequency at lower latitudes than arrangement O st , and the converse is true a higher latitudes [35,36]. For this reason, the variation explained by the six karyotypes was further decomposed after pooling the first two arrangements