- Research article
- Open Access
Genetic constraints for thermal coadaptation in Drosophila subobscura
© Dolgova et al; licensee BioMed Central Ltd. 2010
- Received: 12 July 2010
- Accepted: 25 November 2010
- Published: 25 November 2010
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.
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.
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.
- Genetic Correlation
- Heat Tolerance
- Heat Shock Response
- Isofemale Line
- Thermal Preference
Ectotherms exhibit a suite of behavioural and physiological strategies to cope with spatiotemporal variation in ambient temperature . For instance, behavioural adjustments (e.g. modifying daily activity patterns and selecting favourable microclimates; ) can buffer the impact of sub-optimal temperatures, and are the main means of thermoregulation in small insects [3–5]. Although such adjustments can enable ectotherms to maintain relatively constant body temperatures (Tb) 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–11] suggest that behavioural thermoregulation may be insufficient to fully compensate shifts in environmental temperature .
If behavioural thermoregulation is not fully compensatory and climate variation influences the actual Tb 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 Tb 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" , 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–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 . Contrasting outcomes were also observed for other clinally varying traits - wing cell size and cell number , 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  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 . Besides, different patterns of linkage disequilibrium could result from variability in migration rates between genetically differentiated populations in the various latitudinal clines . 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 (; but see ), it is essential to analyze the genetic basis of thermal preference and temperature-related 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 (Tp: the preferred body temperature in a laboratory thermal gradient, which we expect to correlate with the thermal optimum for performance; ) and knock out temperatures (Tko: the temperature required to knock out a fly in a water-bath) of their carriers . 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 Tp and had a lower Tko, in accordance with the natural patterns; (ii) different chromosomes were responsible for the bulk of the genetic variation in Tp (chromosomes A and O) and Tko (chromosome E); and (iii) Tp and Tko 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 intra- and 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 Tp and Tko 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 .
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 . 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.
In south-western European populations, the most frequent chromosome arrangements for chromosome O are Ost, O3+4, O3+4+7, and O3+4+8 . The first two arrangements show a clear contrasting clinal pattern in original Palaearctic populations, with Ost increasing and O3+4 decreasing in frequency with increasing latitude [35, 56] (Figure 1a). Arrangement O3+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 . 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., ; j = st, 3 + 4, 3 + 4 + 8) were used in the present experiments. Extensive genetic differentiation of up to 4 Mb (i.e., about 15% of the euchromatic portion) has been detected among these arrangements . 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.  the experimental flies were obtained from 54 crosses, which will be referred to as inbred (isogenic: with 18 crosses in total), or as outbred including both structural homokaryotypes ( with 18 cyclically permuted reciprocal crosses in total) and heterokaryotypes (; 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 Tp (pooled average 16.6°C; ) in D. subobscura flies raised at 18°C, and that obtained by Huey and Pascual (23.7°C; ) where flies were raised at 22°C. Even though the flies assayed came from different sources - south-western Europe in Rego et al. , and North America in Huey and Pascual  -, which could account for the observed difference because thermal responses can vary between populations , 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 Tp and Tko 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: inbred crosses , t = 1.21, df = 347, P = 0.226; outbred crosses , t = 0.79, df = 701, P = 0.429. Furthermore, as expected from the low values of the phenotypic correlation, the genetic (karyotypic) correlation for the outbred flies was also close to zero (r k = -0.068, P = 0.914). The conclusion is that both traits are nearly orthogonal to each other (pooled , 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 temperature effects on thermal preference.
Source of variation
Covariate (plate hour)
Cross ⊂ κ
κ × ι
κ × τ
κ × ς
ι × τ
ι × ς
τ × ς
κ × ι × τ
κ × ι × ς
κ × τ × ς
ι × τ × ς
κ × ι × τ ×ς
Scheffé post hoc tests using the mean square of the nested "cross" effect as the error term showed that the thermal preference of Ost/Ost flies was significantly lower when compared to those of O3+4/O3+4 and O3+4+8/O3+4+8 homokaryotypes, which did not differ between 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" (Tset) ranges (central 50% of preferred body temperatures; ) are bounded by 15.1°C - 20.5°C for Ost/Ost karyotypes, and 16.6°C - 22.2°C for the other two karyotypes.
The karyotype × inbreeding interaction arises from the somewhat different behaviour between Ost/Ost and O3+4/O3+4 karyotypes on one side, and O3+4+8/O3+4+8 on the other: for the first two karyotypes Tp was slightly higher in inbred crosses when compared to their outbred counterparts, whereas the opposite was true for the O3+4+8/O3+4+8 karyotype. Average Tp was, however, almost identical for inbred (18.9°C ± 3.6°C) and outbred (18.5°C ± 3.6°C) flies. On the other hand, Ost/Ost flies raised at 22°C had a higher Tp than those raised at 18°C, but no clear trend was observed for O3+4/O3+4 and O3+4+8/O3+4+8 karyotypes.
b) Knock out temperature
Inbreeding and temperature effects on knockout temperature.
Source of variation
Covariate (water bath)
Cross ⊂ κ
κ × ι
κ × τ
κ × ς
ι × τ
ι × ς
τ × ς
κ × ι ×τ
κ × ι × ς
κ × τ ×ς
ι ×τ × ς
κ × ι × τ × ς
Gene arrangement effects in the outbred lines
a) Thermal preference
Karyotype and temperature effects on thermal preference.
Source of variation
Covariate (plate hour)
Cross ⊂ κ
κ × τ
κ × ς
τ × ς
κ × τ × ς
b) Knock out temperature
Karyotype and temperature effects on knockout temperature.
Source of variation
Covariate (water bath)
Cross ⊂ κ
κ × τ
κ × ς
τ × ς
κ × τ × ς
The genetic correlation between Tp and Tko after pooling O3+4 and O3+4+8was 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 Tp and Tko, together with their statistical significance obtained from the appropriate contrasts (Table 3, 4).
c) Average effects on thermal preference
Average effect of chromosome O on thermal preferences (°C).
36°43'N - 4°25'W
Punta Umbría (Spain)
37°10'N - 6°57'W
39°33'N - 2°29'E
39°33'N - 0°34'W
42°13'N - 2°10'E
43°05'N - 2°37'E
43°36'N - 3°53'E
45°26'N - 0°44'E
46°19'N - 7°39'E
48°13'N - 16°22'E
48°32'N - 9°04'E
50°43'N - 4°37'E
Groningen (The Netherlands)
53°13'N - 6°35'E
Combined with our previous results with chromosome A (which is the sex chromosome and additive values can be estimated using males' Tp; ), where gene arrangement Ast exhibits a similar latitudinal pattern than Ost and flies carrying Ast also display a laboratory thermal preference towards colder temperature, the conclusion is that flies inheriting simultaneously Ast and Ost will choose temperatures ranging from approximately 0.5°C - 1.0°C below the average (these estimates assume perfect additivity).
The present results with isogenic lines and their crosses corroborate and extend our previous work with wild flies from south-western Europe . They confirm that arrangements on chromosome O have a biometrical effect on thermal preference in a laboratory temperature gradient, with cold-climate Ost carriers displaying a lower Tp than their warm-climate O3+4 and O3+4+8 counterparts. In addition, Tp and Tko 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 Tko (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 Tp and Tko; (ii) a lack of phenotypic plasticity for Tp according to the temperature at which the flies were raised (18°C and 22°C); and (iii) a substantial effect of developmental temperature on Tko. The absence of inbreeding depression for Tp agrees with the genetic analysis from outbred flies, where a dominance effect after pooling chromosome arrangements O3+4 and O3+4+8 into a single class () was absent (Table 3; note that the differences detected among the three karyotypes, and attributed to some under-dominance, could not be appreciated in the inbreeding analysis because it only included inbred and outbred homokaryotypes). On the other hand, the lack of inbreeding depression for Tko is expected and does not mean anything here, simply because no "gene" effects linked to chromosomal arrangements on chromosome O were detected. At first sight this might be surprising because a well-characterized cellular defence mechanism once environmental temperature approaches the upper 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 . 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 Tp and Tko. 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 . 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 , 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 . 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 . 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 Tp here was (mean ± SD) 18.4°C ± 3.6°C (Tset: 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, Tset:12.4°C - 20.4°C; ) 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 Tp; ). 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, Tset: 21.2°C - 25.9°C; ), 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 Tko makes us confidently conclude that our estimates are indeed closer to the actual Tp 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  - we have also analyzed the association between Tko 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 Tb for active flies along a latitudinal gradient (spanning 12°) come from recent work by Huey and Pascual  in western North America. They found that mean Tb 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 ) - when wild flies are smaller probably due to their higher developmental temperatures and/or crowding conditions ; and crowding is known to affect adult thermal stress resistance in Drosophila  - are dangerously close to the Tko 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: Ost is more frequent towards the evening while chromosomes carrying gene arrangement O3+4 are most frequently sampled at midday . 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 . The lack of association between Tp and Tko in D. subobscura is fully consistent with their genetic basis as independently segregating chromosomes are involved . 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 Tko ), and on the A and O chromosomes (which influence Tp [, this work]).
We now speculate that the mismatch between Tp and Tko 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 O3+4 are more active at midday and, therefore, have a higher risk of a heat shock than Ost and are selected against. On the other hand, assuming Tp corresponds closely with temperatures that maximize fitness O3+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., O3+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  and apparently no seasonal variation at a north-eastern population also in Spain . The point here is that parallel seasonal changes should also be detected for chromosome A since it also affects Tp . In accordance with this prediction, no seasonal cycling was detected for chromosome A in the north-eastern 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.
For ectotherms facing spatiotemporal variation in environmental temperature theory predicts that a coevolution between thermal preference and physiological performance can occur . 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 . 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 , 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 . Our results corroborate that arrangements on chromosome O affect adult thermal preference: flies inheriting the cold-climate Ost 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 O3+4 and O3+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.
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 . 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 O3+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 , 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 days-old 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 (Tp) and knock out temperature (Tko) for each cross (see below). All fly handling was done at room temperature using CO2 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 Tp was measured as previously described . 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 Tp 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 . 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 . Adults were individually placed in sealed empty vials and immersed in water-baths at Tmin = 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 Tmax, defined as the temperature when the last active fly was knocked out (Tmax = 38°C was the upper limit in the assays; median Tmax = 33°C). For each fly Tko was estimated as the temperature taken to knock it out (defined as the onset of muscle spasms; ).
The experimental setup was devised to assay one male and one female from each cross and temperature per day (five blocks) for both Tp and Tko, 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: Tp assay, 5.04 females and 4.80 males; Tko 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
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 (Tp) 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 Tko was assessed in different water-baths.
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 . 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 Tp and Tko 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 O3+4 and O3+4+8 have a higher frequency at lower latitudes than arrangement Ost, 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 into a single class () as follows: between the two heterokaryotypes; among the three karyotypes; and among Ost/Ost, , . The karyotypic values for Tp and Tko were also estimated in the additive-dominance scale [81, 82] after pooling the two chromosome arrangements that share O3+4 (each comparison or contrast between two means has one degree of freedom).
The square of this correlation can be interpreted as that fraction of the total variation among karyotypes that is explained by Ost/Ost, , .
c) Computer software for statistical analysis
The computer programs used for statistical data analyses were MATLAB algebra program environment (ver. 7.0.4 ) together with the collection of tools supplied by the Statistics Toolbox. The statistical software packages STATISTICA version 9  and SPSS version 15  were also used.
We thank three anonymous reviewers for comments on the manuscript. OD is supported by a pre-doctoral fellowship (BES-2007-17438) from the Ministerio de Educación y Ciencia (Spain). CR is supported by a post-doctoral fellowship (SFRH/BPD/39998/2007) from Fundação para a Ciência e a Tecnologia (Portugal). GC is supported by a pre-doctoral fellowship (BES-2007-15096) from the Ministerio de Educación y Ciencia (Spain). ELR is supported by a Ramón y Cajal contract from the Ministerio de Ciencia e Innovación (Spain). This research was supported by grants CGL2006-13423-C01 and CGL2006-13423-C02-02 from the Ministerio de Ciencia y Tecnología (Spain); CGL2009-12912-C03-01, CGL2010-15395 and BFU2009-07564 from the Ministerio de Ciencia e Innovación; and 2009SGR 636 from Generalitat de Catalunya to the Grup de Biologia Evolutiva.
- Angilletta MJ: Thermal Adaptation. A Theoretical and Empirical Synthesis. 2009, Oxford: Oxford Univ PressView ArticleGoogle Scholar
- Stevenson RD: The relative importance of behavioral and physiological adjustments controlling body temperature in terrestrial ectotherms. Am Nat. 1985, 126: 362-386. 10.1086/284423.View ArticleGoogle Scholar
- Casey TM, Hegel JR: Caterpillar setae: insulation for an ectotherm. Science. 1981, 214: 1131-1133. 10.1126/science.214.4525.1131.View ArticlePubMedGoogle Scholar
- Kingsolver JG, Watt WB: Thermoregulatory strategies in Colias butterflies: thermal stress and limits to adaptation in temporally varying environments. Am Nat. 1983, 121: 32-35. 10.1086/284038.View ArticleGoogle Scholar
- Willmer PG, Unwin S: Field analyses of insect heat budget: reflectance, size and heating rates. Oecologia. 1981, 50: 250-255. 10.1007/BF00348047.View ArticleGoogle Scholar
- Huey RB, Hertz PE, Sinervo B: Behavioral drive versus behavioral inertia in evolution: a null model approach. Am Nat. 2003, 161: 357-366. 10.1086/346135.View ArticlePubMedGoogle Scholar
- Dobzhansky Th: Genetics of the Evolutionary Process. 1970, New York: Columbia Univ PressGoogle Scholar
- Rodríguez-Trelles F, Alvarez G, Zapata C: Time-series analysis of seasonal changes of the O inversion polymorphism of Drosophila subobscura. Genetics. 1996, 142: 179-187.PubMed CentralPubMedGoogle Scholar
- Coyne JA, Bundgaard J, Prout T: Geographic variation of tolerance to environmental stress in Drosophila pseudoobscura. Am Nat. 1983, 122: 474-488. 10.1086/284150.View ArticleGoogle Scholar
- Hoffmann AA, Sørensen JG, Loeschcke V: Adaptation of Drosophila to temperature extremes: bringing together quantitative and molecular approaches. J Therm Biol. 2003, 28: 175-216. 10.1016/S0306-4565(02)00057-8.View ArticleGoogle Scholar
- Hoffmann AA: Physiological climatic limits in Drosophila: patterns and implications. J Exp Biol. 2010, 213: 870-880. 10.1242/jeb.037630.View ArticlePubMedGoogle Scholar
- Huey RB, Pascual M: Partial thermoregulatory compensation along a latitudinal cline in a rapidly evolving invasive species. Drosophila subobscura in North America. Ecology. 2009, 90: 1715-1720. 10.1890/09-0097.1.View ArticlePubMedGoogle Scholar
- Sinclair BJ, Vernon P, Klok CJ, Chown SL: Insects at low temperatures: an ecological perspective. Trends Ecol Evol. 2003, 18: 257-262. 10.1016/S0169-5347(03)00014-4.View ArticleGoogle Scholar
- Hofmann GE, Todgham AE: Living in the now: physiological mechanisms to tolerate a rapidly changing environment. Annu Rev Physiol. 2010, 72: 127-145. 10.1146/annurev-physiol-021909-135900.View ArticlePubMedGoogle Scholar
- Huey RB, Bennett AF: Phylogenetic studies of coadaptation: preferred temperatures versus optimal performance temperatures of lizards. Evolution. 1987, 41: 1098-1115. 10.2307/2409194.View ArticleGoogle Scholar
- Angilletta MJ, Bennett AF, Guderley H, Navas CA, Seebacher F, Wilson RS: Coadaptation: a unified principle in evolutionary thermal biology. Physiol Biochem Zool. 2006, 79: 282-294. 10.1086/499990.View ArticlePubMedGoogle Scholar
- Mayr E: Animal Species and Evolution. 1963, Cambridge, Harvard Univ PressView ArticleGoogle Scholar
- Ayala FJ, Serra L, Prevosti A: A grand experiment in evolution Drosophila subobscura colonization of the Americas. Genome. 1989, 31: 246-255.View ArticleGoogle Scholar
- Hoffmann AA, Weeks AR: Climatic selection on genes and traits after 100 year-old invasion: a critical look at the temperate-tropical clines in Drosophila melanogaster from eastern Australia. Genetica. 2007, 129: 133-147. 10.1007/s10709-006-9010-z.View ArticlePubMedGoogle Scholar
- Huey RB, Gilchrist GW, Carlson ML, Berrigan D, Serra L: Rapid evolution of a geographical cline in size in an introduced fly. Science. 2000, 287: 308-309. 10.1126/science.287.5451.308.View ArticlePubMedGoogle Scholar
- Balanyà J, Serra L, Gilchrist GW, Huey RB, Pascual M, Mestres F, Solé E: Evolutionary pace of chromosomal polymorphism in colonizing populations of Drosophila subobscura: an evolutionary time series. Evolution. 2003, 57: 1837-1845.View ArticlePubMedGoogle Scholar
- Rezende EL, Balanyà J, Rodríguez-Trelles F, Rego C, Fragata I, Matos M, Serra L, Santos M: Climate change and chromosomal inversions in Drosophila subobscura. Clim Res. 2010, 43: 103-114. 10.3354/cr00869.View ArticleGoogle Scholar
- Gilchrist G, Jeffers L, West B, Folk D, Suess J, Huey RB: Clinal patterns of desiccation and starvation resistance in ancestral and invading populations of Drosophila subobscura. Evol Appl. 2008, 1: 513-23. 10.1111/j.1752-4571.2008.00040.x.PubMed CentralView ArticlePubMedGoogle Scholar
- Calboli FCF, Gilchrist GW, Partridge L: Different cell size and cell number contribution in two newly established and one ancient body size cline of Drosophila subobscura. Evolution. 2003, 57: 566-573.View ArticlePubMedGoogle Scholar
- Gilchrist GW, Huey RB, Serra L: Rapid evolution of wing size clines in Drosophila subobscura. Genetica. 2001, 112-113: 273-286. 10.1023/A:1013358931816.View ArticlePubMedGoogle Scholar
- Betrán E, Santos M, Ruiz A: Antagonistic pleiotropic effect of second-chromosome inversions on body size and early life-history traits in Drosophila buzzatii. Evolution. 1998, 52: 144-154. 10.2307/2410929.View ArticleGoogle Scholar
- Etterson JR, Shaw RG: Constraint to adaptive evolution in response to global warming. Science. 2001, 294: 151-154. 10.1126/science.1063656.View ArticlePubMedGoogle Scholar
- Fontdevila A, Zapata C, Alvarez G, Sanchez L, Méndez J, Enriquez I: Genetic coadaptation in the chromosomal polymorphism of Drosophila subobscura. I. Seasonal changes of gametic disequilibrium in a natural population. Genetics. 1983, 105: 935-955.PubMed CentralPubMedGoogle Scholar
- Rodríguez-Trelles F: Seasonal cycles of allozyme-by-chromosomal-inversion gametic disequilibrium in Drosophila subobscura. Evolution. 2003, 57: 839-848.View ArticlePubMedGoogle Scholar
- Fragata I, Balanyà J, Rego C, Matos M, Rezende EL, Santos M: Contrasting patterns of phenotypic variation linked to chromosomal inversions in native and colonizing populations of Drosophila subobscura. J Evol Biol. 2010, 23: 112-123. 10.1111/j.1420-9101.2009.01873.x.View ArticlePubMedGoogle Scholar
- Nosil P, Crespi BJ, Sandoval CP, Kirkpatrick M: Migration and the genetic covariance between habitat preference and performance. Am Nat. 2006, 167: E66-E78. 10.1086/499383.View ArticlePubMedGoogle Scholar
- Duckworth RE: The role of behavior in evolution: a search for mechanism. Evol Ecol. 2009, 23: 513-531. 10.1007/s10682-008-9252-6.View ArticleGoogle Scholar
- Dillon ME, Wang G, Garrity PA, Huey RB: Thermal preference in Drosophila. J Therm Biol. 2009, 34: 109-119. 10.1016/j.jtherbio.2008.11.007.PubMed CentralView ArticlePubMedGoogle Scholar
- Rego C, Balanyà J, Fragata I, Matos M, Rezende EL, Santos M: Clinal patterns of chromosomal inversion polymorphisms in Drosophila subobscura are partly associated with thermal preferences and heat stress resistance. Evolution. 2010, 64: 385-397. 10.1111/j.1558-5646.2009.00835.x.View ArticlePubMedGoogle Scholar
- Menozzi P, Krimbas CB: The inversion polymorphism of Drosophila subobscura revisited: synthetic maps of gene arrangement frequencies and their interpretation. J Evol Biol. 1992, 5: 625-641. 10.1046/j.1420-9101.1992.5040625.x.View ArticleGoogle Scholar
- Krimbas CB: The inversion polymorphism of Drosophila subobscura. Drosophila Inversion Polymorphism. Edited by: Krimbas CB, Powell JR. 1992, Boca Raton: CRC Press, 127-220.Google Scholar
- Sperlich D, Feuerbach-Mravlag H, Lange P, Michaelidis A, Pentzos-Daponte A: Genetic load and viability distribution in central and marginal populations of Drosophila subobscura. Genetics. 1977, 86: 835-848.PubMed CentralPubMedGoogle Scholar
- Pegueroles C, Araúz PA, Pascual M, Mestres F: A recombination survey using microsatellites: the O chromosome of Drosophila subobscura. Genetica. 2010, 138: 795-804. 10.1007/s10709-010-9461-0.View ArticlePubMedGoogle Scholar
- Powell JR: Progress and Prospects in Evolutionary Biology. The Drosophila Model. 1997, New York: Oxford Univ PressGoogle Scholar
- Santos J, Serra L, Solé E, Pascual M: FISH mapping of microsatellite loci from Drosophila subobscura and its comparison to related species. Chromosome Res. 2010, 18: 213-226. 10.1007/s10577-010-9112-4.View ArticlePubMedGoogle Scholar
- Laayouni H, García-Franco F, Chávez-Sandoval BE, Trotta V, Beltran S, Corominas M, Santos M: Thermal evolution of gene expression profiles in Drosophila subobscura. BMC Evol Biol. 2007, 7: 42-10.1186/1471-2148-7-42.PubMed CentralView ArticlePubMedGoogle Scholar
- Moltó MD, Pascual L, Martínez-Sebastián MJ, de Frutos R: Genetic analysis of the heat shock response in three Drosophila species of the obscura group. Genome. 1992, 35: 870-880.View ArticlePubMedGoogle Scholar
- Cuenca JB, Galindo MI, Saura AO, Sorsa V, de Frutos R: Ultrastructure of regions containing homologous loci in polytene chromosomes of Drosophila melanogaster and Drosophila subobscura. Chromosoma. 1998, 107: 113-126. 10.1007/s004120050287.View ArticlePubMedGoogle Scholar
- Kunze-Mühl E, Müller E: Weitere Untersuchungen über die chromosomale Struktur und natürlichen Strukturtypen von D. subobscura. Chromosoma. 1958, 9: 559-570. 10.1007/BF02568093.View ArticlePubMedGoogle Scholar
- Parsell DA, Taulien J, Lindquist S: The role of heat-shock proteins in thermotolerance. Phil Trans Roy Soc London Series B. 1993, 339: 279-286. 10.1098/rstb.1993.0026.View ArticleGoogle Scholar
- Krebs RA: A comparison of Hsp70 expression and thermotolerance in adults and larvae of three Drosophila species. Cell Stress Chaperones. 1999, 4: 243-249. 10.1379/1466-1268(1999)004<0243:ACOHEA>2.3.CO;2.PubMed CentralView ArticlePubMedGoogle Scholar
- Bettencourt BR, Kim I, Hoffmann AA, Feder ME: Response to natural and laboratory selection at the Drosophila hsp70 genes. Evolution. 2002, 56: 1796-1801.View ArticlePubMedGoogle Scholar
- McColl G, Hoffmann AA, McKechnie SW: Response of two heat shock genes to selection for knockdown heat resistance in Drosophila melanogaster. Genetics. 1996, 143: 1615-1627.PubMed CentralPubMedGoogle Scholar
- Santos M, Iriarte PF, Céspedes W, Balanyà J, Fontdevila A, Serra L: Swift laboratory thermal evolution of wing shape (but not size) in Drosophila subobscura and its relationship with chromosomal inversion polymorphism. J Evol Biol. 2004, 17: 841-855. 10.1111/j.1420-9101.2004.00721.x.View ArticlePubMedGoogle Scholar
- Santos M, Céspedes W, Balanyà J, Trotta V, Calboli FCF, Fontdevila A, Serra L: Temperature-related genetic changes in laboratory populations of Drosophila subobscura: evidence against simple climatic-based explanations for latitudinal clines. Am Nat. 2005, 165: 258-273. 10.1086/427093.View ArticlePubMedGoogle Scholar
- Santos M: Evolution of total net fitness in thermal lines: Drosophila subobscura likes it "warm". J Evol Biol. 2007, 20: 2361-2370. 10.1111/j.1420-9101.2007.01408.x.View ArticlePubMedGoogle Scholar
- Solé E, Balanyà J, Sperlich D, Serra L: Long-term changes in the chromosomal inversion polymorphism of Drosophila subobscura. I. Mediterranean populations from southwestern Europe. Evolution. 2002, 56: 830-835.View ArticlePubMedGoogle Scholar
- Prevosti A: Chromosomal inversion polymorphism in the southwestern range of Drosophila subobscura distribution area. Genetica. 1974, 45: 111-124. 10.1007/BF01508935.View ArticleGoogle Scholar
- Balanyà J, Solé E, Oller JM, Sperlich D, Serra L: Long-term changes in the chromosomal inversion polymorphism of Drosophila subobscura. II. European populations. J Zool Syst Evol Res. 2004, 42: 191-201. 10.1111/j.1439-0469.2004.00274.x.View ArticleGoogle Scholar
- Munté A, Rozas J, Aguadé M, Segarra C: Chromosomal inversion polymorphism leads to extensive genetic structure: a multilocus survey in Drosophila subobscura. Genetics. 2005, 69: 1573-1581.Google Scholar
- Andolfatto P, Depaulis F, Navarro A: Inversion polymorphism and nucleotide variability in Drosophila. Genet Res. 2001, 77: 1-8. 10.1017/S0016672301004955.View ArticlePubMedGoogle Scholar
- Laayouni H, Hasson E, Santos M, Fontdevila A: Inversion polymorphism and nucleotide variability in different regions of the second-chromosome in Drosophila buzzatii. Mol Biol Evol. 2003, 20: 931-944. 10.1093/molbev/msg099.View ArticlePubMedGoogle Scholar
- Santos M, Iriarte PF, Cespedes W: Genetics and geometry of canalization and developmental stability in Drosophila subobscura. BMC Evol Biol. 2005, 5: 7-10.1186/1471-2148-5-7.PubMed CentralView ArticlePubMedGoogle Scholar
- Yamamoto AH: Temperature preference of Drosophila immigrans and D. virilis: Intra- and inter-population genetic variation. Jpn J Genet. 1994, 69: 67-76. 10.1266/jjg.69.67.View ArticlePubMedGoogle Scholar
- Hertz PE, Huey RB, Stevenson RD: Evaluating temperature regulation by field-active ectotherms: the fallacy of the inappropriate question. Am Nat. 1993, 142: 796-818. 10.1086/285573.View ArticlePubMedGoogle Scholar
- Falconer DS, Mackay TF: Introduction to Quantitative Genetics. 1996, Essex, England: Longman, 4Google Scholar
- Balanyà J, Oller JM, Huey RB, Gilchrist GW, Serra L: Global genetic change tracks global climate warming in Drosophila subobscura. Science. 2006, 313: 1773-1775. 10.1126/science.1131002.View ArticlePubMedGoogle Scholar
- Sørensen JG, Kristensen TN, Loeschcke V: The evolutionary and ecological role of heat shock proteins. Ecol Lett. 2003, 6: 1025-1037. 10.1046/j.1461-0248.2003.00528.x.View ArticleGoogle Scholar
- Jensen LT, Cockerell FE, Kristensen TN, Rako L, Loeschcke V, McKechnie SW, Hoffmann AA: Adult heat tolerance variation in Drosophila melanogaster is not related with Hsp70 expression. J Exp Zool. 2010, 313A: 35-44. 10.1002/jez.573.View ArticleGoogle Scholar
- Feder ME, Blair N, Figueras H: Oviposition site selection: unresponsiveness of Drosophila to cues of potential thermal stress. Anim Behav. 1997, 53: 585-588. 10.1006/anbe.1996.0333.View ArticleGoogle Scholar
- Chown SL, Jumbam KR, Sørensen JG, Terblanche JS: Phenotypic variance, plasticity and heritability estimates of critical thermal limits depend on methodological context. Funct Ecol. 2009, 23: 133-140. 10.1111/j.1365-2435.2008.01481.x.View ArticleGoogle Scholar
- Rezende EL, Tejedo M, Santos M: Estimating the adaptive potential of critical thermal limits: methodological problems and evolutionary implications. Funct Ecol. 2010Google Scholar
- Mitchell KA, Hoffmann AA: Thermal ramping rate influences evolutionary potential and species differences for upper thermal limits in Drosophila. Funct Ecol. 2010, 24: 694-700. 10.1111/j.1365-2435.2009.01666.x.View ArticleGoogle Scholar
- Moreteau B, Morin J-P, Gibert P, Petavy G, Pla E, David JR: Evolutionary changes of nonlinear reaction norms according to thermal adaptation: a comparison of two Drosophila species. C R Acad Sci Paris. 1997, 320: 833-841.View ArticlePubMedGoogle Scholar
- Santos M, Brites D, Laayouni L: Thermal evolution of pre-adult life history traits, geometric size and shape, and developmental stability in Drosophila subobscura. J Evol Biol. 2006, 19: 2006-2021. 10.1111/j.1420-9101.2006.01139.x.View ArticlePubMedGoogle Scholar
- Bubli OL, Imasheva AG, Loeschcke V: Selection for knockdown resistance to heat in Drosophila melanogaster at high and low larval densities. Evolution. 1998, 52: 619-625. 10.2307/2411097.View ArticleGoogle Scholar
- Kari JS, Huey RB: Size and seasonal temperature in free-ranging Drosophila subobscura. J Therm Biol. 2000, 25: 267-272. 10.1016/S0306-4565(99)00096-0.View ArticlePubMedGoogle Scholar
- Sørensen JG, Loeschcke V: Larval crowding in Drosophila melanogaster induces Hsp70 expression, and leads to increased adult longevity and adult thermal stress resistance. J Insect Physiol. 2001, 47: 1301-1307. 10.1016/S0022-1910(01)00119-6.View ArticlePubMedGoogle Scholar
- Savković V, Stamenkovic-Radak M, Andjelkovic M: Diurnal variability of gene arrangement frequencies in Drosophila subobscura populations from two habitats. J Zool Syst Evol Res. 2004, 42: 208-214. 10.1111/j.1439-0469.2004.00273.x.View ArticleGoogle Scholar
- Hedick PW: Genetics of Populations. 2000, Boston: Jones and Bartlett Publ, 2Google Scholar
- de Frutos R, Prevosti A: Temporal changes of chromosomal polymorphism in natural populations of Drosophila subobscura. Genetica. 1984, 63: 181-187. 10.1007/BF00128411.View ArticleGoogle Scholar
- Iriarte FP, Céspedes W, Santos M: Quantitative-genetic analysis of wing form and bilateral asymmetry in isochromosomal lines of Drosophila subobscura using the Procrustes methods. J Genet. 2003, 82: 95-113. 10.1007/BF02715812.View ArticlePubMedGoogle Scholar
- Mestres F, Pegueroles G, Prevosti A, Serra L: Colonization of America by Drosophila subobscura: lethal genes and the problem of the O5 inversion. Evolution. 1990, 44: 1823-1836. 10.2307/2409511.View ArticleGoogle Scholar
- Lutterschmidt WI, Hutchinson VH: The critical thermal maximum: data to support the onset of spasms as the definitive end point. Can J Zool. 1997, 75: 1553-1560. 10.1139/z97-782.View ArticleGoogle Scholar
- Edgington ES: Randomization Tests. 1995, New York: Marcel Dekker, 3Google Scholar
- Mather K, Jinks JL: Introduction to Biometrical Genetics. 1977, London: Chapman and HallView ArticleGoogle Scholar
- Mather K, Jinks JL: Biometrical Genetics: The Study of Continuous Variation. 1982, London: Chapman and Hall, 3View ArticleGoogle Scholar
- Sokal RR, Rohlf FJ: Biometry. 1995, New York: Freeman, 3Google Scholar
- MathWorks Inc: Matlab, V. 184.108.40.2062 (R2007b). The Language of Technical Computing. 2007, [http://www.mathworks.com]Google Scholar
- StatSoft Inc: STATISTICA (data analysis software system), version 9. 2009, [http://www.statsoft.com]Google Scholar
- SPSS Inc: SPSS for Windows. 2006, [http://www.spss.com]Google Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<url>http://creativecommons.org/licenses/by/2.0</url>), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.