Developmental thermal plasticity among Drosophila melanogaster populations

Many biotic and abiotic variables influence the dispersal and distribution of organisms. Temperature has a major role in determining these patterns because it changes daily, seasonally and spatially, and these fluctuations have a significant impact on an organism's behaviour and fitness. Most ecologically relevant phenotypes that are adaptive are also complex and thus they are influenced by many underlying loci that interact with the environment. In this study, we quantified the degree of thermal phenotypic plasticity within and among populations by measuring chill‐coma recovery times of lines reared from egg to adult at two different environmental temperatures. We used sixty genotypes from six natural populations of Drosophila melanogaster sampled along a latitudinal gradient in South America. We found significant variation in thermal plasticity both within and among populations. All populations exhibit a cold acclimation response, with flies reared at lower temperatures having increased resistance to cold. We tested a series of environmental parameters against the variation in population mean thermal plasticity and discovered the mean thermal plasticity was significantly correlated with altitude of origin of the population. Pairing our data with previous experiments on viability fitness assays in the same populations in fixed and variable environments suggests an adaptive role of this thermal plasticity in variable laboratory environments. Altogether, these data demonstrate abundant variation in adaptive thermal plasticity within and among populations.     

Dissecting chill coma recovery as a measure of cold resistance: evidence for a biphasic response in Drosophila melanogaster

Cold resistance in insects has traditionally been measured in terms of survival following a stress, but alternative methods are increasingly being used because of their relevance to the ecology of organisms and their utility in characterizing variation among species, populations and individuals. One such method capable of discriminating among Drosophila species and conspecific Drosophila populations from different environments is adult chill coma recovery time, the time taken for adults to become active again after being knocked down by a cold stress. Here we characterized the chill coma response of D.melanogaster in detail. Adults were exposed to a range of temperatures and stressful periods prior to measuring recovery. Recovery from chill coma in D.melanogaster was biphasic; as flies were stressed under cooler temperatures, recovery times leveled off and then decreased before sharply increasing again as mortality starts to occur. This biphasic response has previously been observed in D.subobscura where it has a somewhat different shape. A second mechanism therefore acts at relatively lower temperatures to ameliorate the effects of the cold stress. When D.melanogaster were reared at 19 and 25 °C for two generations, the shape of the curve relating temperature to recovery time was similar, but flies from the warmer temperature had longer recovery times and showed responses that leveled off and then decreased at relatively higher temperatures. As exposure time to cold stress was increased, recovery times also increased except at mild stress levels. Chill coma recovery in D.melanogaster is a complex trait and likely to reflect multiple underlying components.     

Intraspecific variation in thermal acclimation and tolerance between populations of the winter ant,Prenolepis imparis

Thermal phenotypic plasticity, otherwise known as acclimation, plays an essential role in how organisms respond to short‐term temperature changes. Plasticity buffers the impact of harmful temperature changes; therefore, understanding variation in plasticity in natural populations is crucial for understanding how species will respond to the changing climate. However, very few studies have examined patterns of phenotypic plasticity among populations, especially among ant populations. Considering that this intraspecies variation can provide insight into adaptive variation in populations, the goal of this study was to quantify the short‐term acclimation ability and thermal tolerance of several populations of the winter ant, Prenolepis imparis. We tested for correlations between thermal plasticity and thermal tolerance, elevation, and body size. We characterized the thermal environment both above and below ground for several populations distributed across different elevations within California, USA. In addition, we measured the short‐term acclimation ability and thermal tolerance of those populations. To measure thermal tolerance, we used chill‐coma recovery time (CCRT) and knockdown time as indicators of cold and heat tolerance, respectively. Short‐term phenotypic plasticity was assessed by calculating acclimation capacity using CCRT and knockdown time after exposure to both high and low temperatures. We found that several populations displayed different chill‐coma recovery times and a few displayed different heat knockdown times, and that the acclimation capacities of cold and heat tolerance differed among most populations. The high‐elevation populations displayed increased tolerance to the cold (faster CCRT) and greater plasticity. For high‐temperature tolerance, we found heat tolerance was not associated with altitude; instead, greater tolerance to the heat was correlated with increased plasticity at higher temperatures. These current findings provide insight into thermal adaptation and factors that contribute to phenotypic diversity by revealing physiological variance among populations.     

Effects of acclimation and latitude on the activity thresholds of the aphid Myzus persicae in Europe

Activity thresholds were measured in nine anholocyclic clones of the peach‐potato aphid Myzus persicae collected along a latitudinal cline of its European distribution from Sweden to Spain. The effects of collection origin and intra‐ and intergenerational acclimation on these thresholds were investigated. Low‐temperature (10°C) acclimation for one generation depressed the movement threshold and chill coma temperatures, with the largest reduction in movement threshold recorded for clone UK 1 (8.8–2.5°C) and in chill coma for UK 2 (4.8–2.0°C). High‐temperature (25°C) acclimation for one generation increased the heat movement threshold and heat coma temperature with the largest increase in the movement threshold (40.1–41.1°C) and heat coma (41.4–42.3°C) recorded for clone Swed 1. There was no further intergenerational acclimation over three generations. High‐temperature activity thresholds were less plastic than low‐temperature thresholds, and, consequently, thermal activity ranges were expanded following low‐temperature acclimation. No constant affect of acclimation was observed on chill coma recovery, although clonal differences were observed with Swed 1 and 3 requiring some of the longest complete recovery times. There was no relationship between latitude and activity thresholds with the exception of heat coma data where Scandinavian clones Swed 2 and 3 consistently displayed some of the lowest heat coma temperatures (e.g. 41.3°C for both clones at 20°C) and Mediterranean clones Span 1, 2 and 3 displayed some of the highest (e.g. 42.1, 41.9 and 42.5°C, respectively, at 20°C). These data suggest that clonal mixing could occur over a large scale across Europe, limiting local adaptation to areas where conditions enable long‐term persistence of populations, e.g. adaptation to higher temperatures in the Mediterranean region. It is suggested that aphid thermal tolerance could be governed more by clonal type than the latitudinal origin.     

Rapid growth reduces cold resistance: evidence from latitudinal variation in growth rate, cold resistance and stress proteins

Background

Physiological costs of rapid growth may contribute to the observation that organisms typically grow at submaximal rates. Although, it has been hypothesized that faster growing individuals would do worse in dealing with suboptimal temperatures, this type of cost has never been explored empirically. Furthermore, the mechanistic basis of the physiological costs of rapid growth is largely unexplored.

Methodology/Principal Finding

Larvae of the damselfly Ischnura elegans from two univoltine northern and two multivoltine southern populations were reared at three temperatures and after emergence given a cold shock. Cold resistance, measured by chill coma recovery times in the adult stage, was lower in the southern populations. The faster larval growth rates in the southern populations contributed to this latitudinal pattern in cold resistance. In accordance with their assumed role in cold resistance, Hsp70 levels were lower in the southern populations, and faster growing larvae had lower Hsp70 levels. Yet, individual variation in Hsp70 levels did not explain variation in cold resistance.

Conclusions/Significance

We provide evidence for a novel cost of rapid growth: reduced cold resistance. Our results indicate that the reduced cold resistance in southern populations of animals that change voltinism along the latitudinal gradient may not entirely be explained by thermal selection per se but also by the costs of time constraint-induced higher growth rates. This also illustrates that stressors imposed in the larval stage may carry over and shape fitness in the adult stage and highlights the importance of physiological costs in the evolution of life-histories at macro-scales.     

Cold adaptation shapes the robustness of metabolic networks in Drosophila melanogaster

When ectotherms are exposed to low temperatures, they enter a cold‐induced coma (chill coma) that prevents resource acquisition, mating, oviposition, and escape from predation. There is substantial variation in time taken to recover from chill coma both within and among species, and this variation is correlated with habitat temperatures such that insects from cold environments recover more quickly. This suggests an adaptive response, but the mechanisms underlying variation in recovery times are unknown, making it difficult to decisively test adaptive hypotheses. We use replicated lines of Drosophila melanogaster selected in the laboratory for fast (hardy) or slow (susceptible) chill‐coma recovery times to investigate modifications to metabolic profiles associated with cold adaptation. We measured metabolite concentrations of flies before, during, and after cold exposure using nuclear magnetic resonance (NMR) spectroscopy to test the hypotheses that hardy flies maintain metabolic homeostasis better during cold exposure and recovery, and that their metabolic networks are more robust to cold‐induced perturbations. The metabolites of cold‐hardy flies were less cold responsive and their metabolic networks during cold exposure were more robust, supporting our hypotheses. Metabolites involved in membrane lipid synthesis, tryptophan metabolism, oxidative stress, energy balance, and proline metabolism were altered by selection on cold tolerance. We discuss the potential significance of these alterations.     

Evolution of thermotolerance in seasonal environments: the effects of annual temperature variation and life-history timing in Wyeomyia smithii

In organisms with complex life cycles, the adaptive value of thermotolerance depends on life-history timing and seasonal temperature profiles. We illustrate this concept by examining variation in annual thermal environments and thermal acclimation among four geographic populations of the pitcher plant mosquito. Only diapausing larvae experience winter, whereas both postdiapause and nondiapause adults occur only during the growing season. Thus, adults experience transient cold stress primarily during the spring. We show that adult cold tolerance (chill coma recovery) is enhanced in spring-like conditions via thermal acclimation but is unaffected by diapause state. Moreover, adult mosquitoes from northern populations were more cold tolerant than those from southern populations largely because acclimation responses were steeper in the north. In contrast to cold tolerance, there was no significant acclimation of heat tolerance (heat knockdown), and no significant differences in heat tolerance between northern and southern populations. Field temperature data show that because of evolved differences in diapause timing, adult exposure to cold stress is remarkably consistent across geography. This suggests that geographic variation in cold tolerance may not be the result of direct selection on adults. Our results illustrate the importance of the interplay between phenological and thermal adaptation for understanding variation along climatic gradients.     

Thermal tolerance in a south-east African population of the tsetse fly Glossina pallidipes (Diptera, Glossinidae): implications for forecasting climate change impacts

For tsetse (Glossina spp.), the vectors of human and animal trypanosomiases, the physiological mechanisms linking variation in population dynamics with changing weather conditions have not been well established. Here, we investigate high- and low-temperature tolerance in terms of activity limits and survival in a natural population of adult Glossina pallidipes from eastern Zambia. Due to increased interest in chilling flies for handling and aerial dispersal in sterile insect technique control and eradication programmes, we also provide further detailed investigation of low-temperature responses. In wild-caught G. pallidipes, the probability of survival for 50% of the population at low-temperatures was at 3.7, 8.9 and 9.6 degrees C (95% CIs: +/-1.5 degrees C) for 1, 2 and 3 h treatments, respectively. At high temperatures, it was estimated that treatments at 37.9, 36.2 and 35.6 degrees C (95% CIs: +/-0.5 degrees C) would yield 50% population survival for 1, 2 and 3 h, respectively. Significant effects of time and temperature were detected at both temperature extremes (GLZ, p<0.05 in all cases) although a time-temperature interaction was only detected at high temperatures (p<0.0001). We synthesized data from four other Kenyan populations and found that upper critical thermal limits showed little variation among populations and laboratory treatments (range: 43.9-45.0 degrees C; 0.25 degrees C/min heating rate), although reduction to more ecologically relevant heating rates (0.06 degrees C/min) reduce these values significantly from approximately 44.4 to 40.6 degrees C, thereby providing a causal explanation for why tsetse distribution may be high-temperature limited. By contrast, low-temperature limits showed substantial variation among populations and acclimation treatments (range: 4.5-13.8 degrees C; 0.25 degrees C/min), indicating high levels of inter-population variability. Ecologically relevant cooling rates (0.06 degrees C/min) suggest tsetses are likely to experience chill coma temperatures under natural conditions (approximately 20-21 degrees C). The results from acute hardening experiments in the Zambian population demonstrate limited ability to improve low-temperature tolerance over short (hourly) timescales after non-lethal pre-treatments. In flies which survived chilling, recovery times were non-linear with plateaus between 2-6 and 8-12 degrees C. Survival times ranged between 4 and 36 h and did not vary between flies which had undergone chill coma by comparison with flies which had not, even after factoring body condition into the analyses (p>0.5 in all cases). However, flies with low chill coma values had the highest body water and fat content, indicating that when energy reserves are depleted, low-temperature tolerance may be compromised. Overall, these results suggest that physiological mechanisms may provide insight into tsetse population dynamics, hence distribution and abundance, and support a general prediction for reduced geographic distribution under future climate warming scenarios.     

Thermal tolerance and recovery behaviour of Thorectes lusitanicus (Coleoptera,Geotrupidae)

1. Population differences in physiological responses are examined in Thorectes lusitanicus, an endemic Iberian dung beetle species, by submitting individuals of different populations to the same experimental and acclimation conditions. 2. An infrared thermography protocol was used, consisting of three assays: start of activity, cold response, and heat response. Individuals of 12 populations were studied and the comparative explanatory capacities of several environmental factors in relation to the observed inter‐population differences were examined. 3. The heating rate from chill coma to the beginning of activity was the variable with the highest discrimination power among the studied populations, accounting for 94% of the observed variance. Regarding the heat response, only six of the 16 thermal variables reached significance (inter‐population differences accounted for 52–74% in these six thermal parameters). 4. From the three considered environmental factors (Mediterranean climate, land cover, and trophic characteristics) only land cover characteristics remain statistically significant, affecting the cold response of individuals. 5. Thorectes lusitanicus is a species characterised by a high diversity of thermotolerance and recovery traits across populations with a low degree of association with broad environmental factors. Finally, it is suggested that the apterous character of this species could be a determinant factor explaining the high diversity of ecophysiological traits related to thermal stress tolerance and the recovery time.     

Complexity of the cold acclimation response in Drosophila melanogaster

Insects can increase their resistance to cold stress when they are exposed to non-lethal conditions prior to the stress; these plastic responses are normally described only in terms of immediate effects on mortality. Here we examine in Drosophila melanogaster the short- and longer-term effects of different conditions on several measures of cold resistance, but particularly chill coma recovery. Short-term exposure to sublethal temperature (cold hardening) did not decrease chill coma recovery times even though it decreased mortality. Exposure to 12 degrees C for 2 days (acclimation) decreased chill coma recovery times for a range of stressful temperatures when flies were cultured at 25 degrees C, but did not usually affect recovery times when flies were cultured at 19 degrees C. In contrast, 2-day exposure to 12 degrees C decreased mortality regardless of rearing temperature. Rearing at 19 degrees C decreased mortality and chill coma recovery time relative to rearing at 25 degrees C. Acclimation increased the eclosion rate of eggs from stressed females, but did not affect development time or size of the offspring. These results indicate that plastic responses to cold in D. melanogaster are complex when resistance is scored in different ways, and that effects can extend across generations.     

Intraspecific differences in thermal tolerance of the diamondback watersnake (Nerodia rhombifer): effects of ontogeny, latitude, and sex

Ontogenetic shifts in microhabitat use are widespread among taxa and can result in drastic shifts in thermal habitat among age classes. Likewise, geographic variation in climate along latitudinal gradients can cause differences in thermal environments among populations of a species. Using a common garden design, we examined four populations of a single species of semi-aquatic snake, Nerodia rhombifer, to determine whether ontogenetic shifts in habitat use (and/or body size) and latitudinal differences in ambient temperature have resulted in evolutionary changes in thermal tolerance. We found ontogenetic differences in thermal tolerance for all populations, with neonates tolerating temperatures 2 degrees C higher than adults, a pattern that is consistent with ontogenetic shifts in body size and microhabitat use in this species. There were differences in thermal tolerance among latitudes in neonates, suggesting genetic differences among populations, but adults showed no latitudinal differences. In combination, the increased thermal tolerance of neonates and the age-specific response to latitude suggest individuals may be most sensitive to selection on thermal tolerance as neonates. Although latitudinal differences exist in neonates, their tolerances were not ranked according to latitude, suggesting the effects of some other local factor (e.g., microclimate) may be important. Lastly, among neonates, females tolerate higher temperatures than males.     

Accounting for Genetic Architecture Improves Sequence Based Genomic Prediction for a Drosophila Fitness Trait

The ability to predict quantitative trait phenotypes from molecular polymorphism data will revolutionize evolutionary biology, medicine and human biology, and animal and plant breeding. Efforts to map quantitative trait loci have yielded novel insights into the biology of quantitative traits, but the combination of individually significant quantitative trait loci typically has low predictive ability. Utilizing all segregating variants can give good predictive ability in plant and animal breeding populations, but gives little insight into trait biology. Here, we used the Drosophila Genetic Reference Panel to perform both a genome wide association analysis and genomic prediction for the fitness-related trait chill coma recovery time. We found substantial total genetic variation for chill coma recovery time, with a genetic architecture that differs between males and females, a small number of molecular variants with large main effects, and evidence for epistasis. Although the top additive variants explained 36% (17%) of the genetic variance among lines in females (males), the predictive ability using genomic best linear unbiased prediction and a relationship matrix using all common segregating variants was very low for females and zero for males. We hypothesized that the low predictive ability was due to the mismatch between the infinitesimal genetic architecture assumed by the genomic best linear unbiased prediction model and the true genetic architecture of chill coma recovery time. Indeed, we found that the predictive ability of the genomic best linear unbiased prediction model is markedly improved when we combine quantitative trait locus mapping with genomic prediction by only including the top variants associated with main and epistatic effects in the relationship matrix. This trait-associated prediction approach has the advantage that it yields biologically interpretable prediction models.     

Variation in Thermal Sensitivity and Thermal Tolerances in an Invasive Species across a Climatic Gradient: Lessons from the Land Snail Cornu aspersum

The ability of organisms to perform at different temperatures could be described by a continuous nonlinear reaction norm (i.e., thermal performance curve, TPC), in which the phenotypic trait value varies as a function of temperature. Almost any shift in the parameters of this performance curve could highlight the direct effect of temperature on organism fitness, providing a powerful framework for testing thermal adaptation hypotheses. Inter-and intraspecific differences in this performance curve are also reflected in thermal tolerances limits (e.g., critical and lethal limits), influencing the biogeographic patterns of species’ distribution. Within this context, here we investigated the intraspecific variation in thermal sensitivities and thermal tolerances in three populations of the invasive snail Cornu aspersum across a geographical gradient, characterized by different climatic conditions. Thus, we examined population differentiation in the TPCs, thermal-coma recovery times, expression of heat-shock proteins and standard metabolic rate (i.e., energetic costs of physiological differentiation). We tested two competing hypotheses regarding thermal adaptation (the “hotter is better” and the generalist-specialist trade-offs). Our results show that the differences in thermal sensitivity among populations of C. aspersum follow a latitudinal pattern, which is likely the result of a combination of thermodynamic constraints (“hotter is better”) and thermal adaptations to their local environments (generalist-specialist trade-offs). This finding is also consistent with some thermal tolerance indices such as the Heat-Shock Protein Response and the recovery time from chill-coma. However, mixed responses in the evaluated traits suggest that thermal adaptation in this species is not complete, as we were not able to detect any differences in neither energetic costs of physiological differentiation among populations, nor in the heat-coma recovery.     

The effect of selection for desiccation resistance on cold tolerance of Drosophila melanogaster

Abstract.  Low temperature and desiccation stress are thought to be mechanistically similar in insects, and several studies indicate that there is a degree of cross-tolerance between them, such that increased cold tolerance results in greater desiccation tolerance and vice versa . This assertion is tested at an evolutionary scale by examining basal cold tolerance, rapid cold-hardening (RCH) and chill coma recovery in replicate populations of Drosophila melanogaster selected for desiccation resistance (with controls for both selection and concomitant starvation) for over 50 generations. All of the populations display a RCH response, and there is no effect of selection regime on RCH or basal cold tolerance, although there are differences in basal cold tolerance between sampling dates, apparently related to inter-individual variation in development time. Flies selected for desiccation tolerance recover from chill coma slightly, but significantly, faster than control and starvation-control flies. These findings provide little support for cross-tolerance between survival of near-lethal cold and desiccation stress in D. melanogaster .     

Impact of different acclimation temperatures and duration on the chill coma temperature and oxygen consumption in the tenebrionid beetle Alphitobius diaperinus

When the ambient temperature is lowered to an insect's lower thermal limit, the insect enters into chill coma. Chill coma temperature and chill coma recovery can vary within species as a result of thermal acclimation, although the physiological basis of the onset of chill coma remains poorly understood. The present study investigates how the temperature of acclimation (0, 5, 10, 15 and 20 °C for 2 or 7 days) affects chill coma temperature and oxygen consumption in adult Alphitobius diaperinus Panzer (Coleoptera: Tenebrionidae). It is hypothesized that the threshold decline in metabolic rate corresponds to the entry into chill coma. Oxygen consumption (as a proxy of metabolism) is measured across the chill coma temperature threshold, and a strong decline in oxygen consumption is expected at entry into chill coma. The acclimation decreases the chill coma temperature significantly from 6.6 ± 1.1 °C in control insects to 3.1 ± 0.7 °C in those acclimated to 10 °C. The change in metabolic rate (Q₁₀) after acclimation to temperatures ranging from 10 to 20 °C is 3.7. Despite acclimation, the metabolic rate of A. diaperinus conforms to Arrhenius kinetics, suggesting that the response of this beetle does not show metabolic compensation. The data suggest the existence of a threshold decline in metabolic rate during cooling that coincides with the temperature at which an insect goes into chill coma.     

Latitudinal and voltinism compensation shape thermal reaction norms for growth rate

Latitudinal variation in thermal reaction norms of key fitness traits may inform about the response of populations to climate warming, yet their adaptive nature and evolutionary potential are poorly known. We assessed the contribution of quantitative genetic, neutral genetic and environmental effects to thermal reaction norms of growth rate for populations of the damselfly Ischnura elegans. Among populations, reaction norms differed primarily in elevation, suggesting that time constraints associated with shorter growth seasons in univoltine, high-latitude as well as multivoltine, low-latitude populations selected for faster growth rates. Phenotypic divergence among populations is consistent with selection rather than drift as Q(ST) was greater than F(ST) in all cases. Q(ST) estimates increased with experimental temperature and were influenced by genotype by environment interactions. Substantial additive genetic variation for growth rate in all populations suggests that evolution of trait means in different environments is not constrained. Heritability of growth rates was higher at high temperature, driven by increased genetic rather than environmental variance. While environment-specific nonadditive effects also may contribute to heritability differences among temperatures, maternal effects did not play a significant role (where these could be accounted for). Genotype by environment interactions strongly influenced the adaptive potential of populations, and our results suggest the potential for microevolution of thermal reaction norms in each of the studied populations. In summary, the observed latitudinal pattern in growth rates is adaptive and results from a combination of latitudinal and voltinism compensation. Combined with the evolutionary potential of thermal reaction norms, this may affect populations' ability to respond to future climate warming.     

Influence of developmental temperature on cold shock and chill coma recovery in Drosophila ananassae: Acclimation and latitudinal variations among Indian populations

Exposure of various Drosophila species to mild increase or decrease in temperature has consistently been shown to result in increased resistance to subsequent temperature extremes. We investigated cold tolerance in 45 Indian natural populations of Drosophila ananassae collected from all over India by monitoring the time taken by adults to recover from chill-coma after a treatment for 16 h at 4 °C. Significant latitudinal and altitudinal differentiation was observed for chill coma recovery in D. ananassae. Chill-coma recovery was closely associated with local climatic factors like average annual temperature and relative humidity of origin of populations.     

Clinal variation in Drosophila serrata for stress resistance and body size

Clines for size and stress resistance traits have been described for several Drosophila species and replicable clines across different species may indicate climatic selection. Here we consider clines in stress resistance traits in an Australian endemic species, D. serrata, by comparing levels of variation within and among isofemale lines initiated with flies collected from the eastern coast of Australia. We also consider clinical variation in chill coma recovery, a trait that has recently been shown to exhibit high levels of variation among Drosophila species. Patterns were compared with those in the cosmopolitan species D. melanogaster from the same area. Both desiccation and starvation resistance showed no clinical pattern despite heritable variation among isofemale lines. In contrast chill coma resistance exhibited a linear cline in the anticipated direction, resistance increasing with latitude. Body size was measured as wing length and body weight. Both traits showed geographic variation and strong non-linear clines with a sharp reduction in size in the tropics. These results are discussed in the context of climatic selection and evolutionary processes limiting species borders.     

Contrasting dynamics of cold resistance traits in field-fresh Myrmica ants during the active season

As a result of acclimation populations of long-lived ectotherms should display lowered ability to counter cold stress in warmer periods of active season, and increased resistance in colder ones. We tested this proposition by investigating dynamics of cold resistance in Myrmica ants during most of the active season in two types of habitats. Resistance of ants to knock-down by cold and their rate of recovery after chill coma were expected to be lower in summer.Cooled at a rate of 0.17 °C min⁻¹, the ants showed lower capability to resist knock-down in summer, and a significant lowering in knock-down temperature in response to colder weather both in spring and autumn as confirmed by linear regression against air temperatures. In a more eurytopic species M. rubra the responses were significantly faster in meadow than in forest habitats. However, times of recovery of the ants after 10 min at −3 °C did not change in parallel to air temperatures. Whereas M. rubra from forest habitats took less time to recover in early summer and early autumn, in their conspecifics from meadow habitats the contrary was the case. Regardless of habitat, recoveries tended to be faster in other investigated species, of which M. ruginodis (a forest stenotopic) recovered faster in early summer than later.According to the knock-down data, in warmer months the ants are indeed less resistant to cold stress, whilst the recovery data do not always support the proposition. The contrasting seasonal dynamics of the two measures of low-temperature resistance in field-fresh Myrmica suggest that knock-down (chill coma onset) is a better index of thermal acclimation, whilst the rate of recovery from chill coma is more indicative of interspecific differences and, possibly, behavioural thermoregulation.