{'fraction_non_alphanumeric': 0.02985510065156083, 'fraction_numerical': 0.012933968686181076, 'mean_word_length': 4.899885255306942, 'pattern_counts': {'":': 0, '<': 4, '<?xml version=': 0, '>': 2, 'https://': 0, 'lorem ipsum': 0, 'www.': 6, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Chill coma onset and recovery fail to reveal true variation in thermal performance among populations of Drosophila melanogaster

Species from colder climates tend to be more chill tolerant regardless of the chill tolerance trait measured, but for Drosophila melanogaster, population-level differences in chill tolerance among populations are not always found when a single trait is measured in the laboratory. We measured chill coma onset temperature, chill coma recovery time, and survival after chronic cold exposure in replicate lines derived from multiple paired African and European D. melanogaster populations. The populations in our study were previously found to differ in chronic cold survival ability, which is believed to have evolved independently in each population pair; however, they did not differ in chill coma onset temperature and chill coma recovery time in a manner that reflected their geographic origins, even though these traits are known to vary with origin latitude among Drosophila species and are among the most common metrics of thermal tolerance in insects. While it is common practice to measure only one chill tolerance trait when comparing chill tolerance among insect populations, our results emphasise the importance of measuring more than one thermal tolerance trait to minimize the risk of missing real adaptive variation in insect thermal tolerance.

more chill tolerant overall (i.e. we often see correlations among chill tolerance traits as well as between individual traits and geographic distribution) than species or populations from warmer regions 1,25 . Among Drosophila species, CCO correlates strongly with two measures of survival (lethal temperature and lethal time at low temperature), weakly with CCRT, and these traits also correlate with geographical distribution (albeit only weakly, in the case of CCRT) 2 . Within a single species, Drosophila melanogaster, populations from colder regions are sometimes, but not always, more chill tolerant based on one or more of these traits. For example, higher latitude Australian populations are more chill tolerant based on CCO, CCRT and cold shock survival [26][27][28] . Similarly, in a large study of D. melanogaster populations from around the world, temperate populations were more chill tolerant based on CCRT (the only measure used) than tropical populations 29 ; however, in a separate study of African and European D. melanogaster populations, critical thermal minimum was not related to latitude (although it was in a related species, D. subobscura) 30 , and in Japan, there is minimal variation in cold shock survival and CCO among northerly/southerly populations within Drosophila species (although species with more northerly distributions are more chill tolerant based on those same measures) 3,31 . Even focusing on a single population and a single trait can lead to conflicting evidence of variation in thermal performance. For example, survival following cold stress at a given temperature among 40 lines of the Drosophila Genetic Reference Panel (DGRP) does not always predict survival at another low temperature 32 .
Patterns of variation of chill tolerance within Drosophila species can be further complicated by not knowing how previously-tested population sets are related: CCO, at least, has a very strong phylogenetic signal 33 , which can cloud interpretation if not accounted for. One system for which phylogenetic relationships are known is the set of African and European populations of D. melanogaster used by Pool et al. 34 to study the parallel evolution of chill tolerance traits. The system consists of paired, closely-related populations and a single outgroup (Zambia; D. melanogaster originally evolved in the warm, tropical forests of Zambia and Zimbabwe 35,36 ). In each pair, one population is derived from flies caught in a relatively cold region, whereas the other is derived from flies caught in a relatively warm region (Fig. 1). Pool et al. (2017) used a chronic cold exposure assay (4 days at 4 °C) to demonstrate that populations from cold regions are more chill tolerant, and concluded that this tolerance likely evolved independently in each coldregion population. This system provides an opportunity to also examine whether and how other chill tolerance traits have evolved in cold vs. warm climates, without the concern that observed chill tolerance in separate regions might be inherited from a common ancestor rather than independently evolved. We thus set out to further characterise these populations using two of the most common chill tolerance measures-chill coma onset temperature (CCO) and chill coma recovery time (CCRT). Because chill tolerance traits often do correlate both with each other and with climate of origin in natural systems, we expected that cold-climate populations would be more chill tolerant than their warm-climate counterparts based on one or both of these traits. To our surprise, the populations differed very little based on either the CCO or CCRT assays-within or among pairs. To ensure that differences in chill tolerance had not simply been lost via lab adaptation, we further characterised survival www.nature.com/scientificreports/ during a 4 day exposure to 4 °C in French and Egyptian lines (the pair previously found to exhibit the greatest difference in recovery after chronic cold exposure 34 ) and confirmed clear differences in at least one measure of chill tolerance persist in these populations.

Results
Chill coma onset temperature (CCO) and chill coma recovery time (CCRT). Population was a significant predictor of CCO (p < 0.001), but this was driven entirely by the Zambian population (our outgroup), which had a significantly lower CCO than all other populations (p < 0.001; Fig. 2). There were no significant differences among the other populations (p > 0.05 in each pairwise comparison). Similarly, population was a significant predictor of CCRT (p = 0.003), once fluctuations in room temperature were taken into account (ambient temperature has a significant impact on chill coma recovery time (p < 0.001)). However, in the pairwise comparisons, CCRTs differed significantly only between Zambia and both the French www.nature.com/scientificreports/ (p = 0.01) and Egyptian (p = 0.01) populations; there were no significant differences within population pairs (Fig. 3). There was also no significant pattern in the number of flies that did not stand during the observation period (p > 0.05; Fig. 3). Likewise, CCRT did not correlate with CCO (R 2 = 0.00, p = 0.40), even when CCRTs were adjusted based on room temperature (R 2 = 0.02, p = 0.30) (Fig. 4).

Survival.
French lines had, on the whole, higher survival scores over time at 4 °C than the Egyptian lines. The separation between populations was statistically significant on days one through four (p < 0.001) and was largest after two days at 4 °C (Fig. 5).

Discussion
In this study, we confirmed that chill tolerance traits can vary independently in Drosophila melanogaster. Similar results have emerged from laboratory selection experiments (where selection pressure is controlled) and using the DRGP system to focus on genetic variation within a single D. melanogaster population 20,21 , but never, to our knowledge, in multiple wild-derived populations like those used here. We were surprised to find that chill tolerance traits vary independently even among populations that are known to differ significantly in one chill tolerance trait 34 -in a way that matched the climatic conditions where their ancestors were collected. In some www.nature.com/scientificreports/ ways, this finding is not surprising: it is increasingly clear that the underlying mechanisms of several commonly measured chill tolerance traits differ 6,9 , and it is also possible to select on such traits separately 21 . However, it remains unclear why some wild-derived populations from colder regions display greater chill tolerance across a range of traits [26][27][28] , whereas others (e.g. the populations studied here) can only be differentiated if we measure the "right" trait. In other words, in order to draw meaningful conclusions about differences in stress tolerance (or the lack thereof), it is important to measure more than one trait 20 . Another surprise was the relatively low chill coma onset temperature (but unremarkable chill coma recovery time) of the Zambian population, derived from flies caught in the very warm ancestral range of this species. We do not have an explanation for this finding, but it may warrant further study. There are several possibilities for why some cold climate populations might not evolve improvements in specific chill tolerance traits. First, it is possible that after evolving better chill tolerance in one trait, or other adaptations to a cold climate that we were unable to detect, these populations did not face the necessary selection pressure to drive improvements in other traits. Notably, however, at an interspecific level this does not appear to be the case: Drosophila species from cold climates do not evolve one chill tolerance trait and then stop 2 . Second,  www.nature.com/scientificreports/ there may not have been enough evolutionary time for the populations to fully differentiate: D. melanogaster only left Sub-Saharan Africa some ten thousand years ago 37 , and it is possible that some chill tolerance traits may respond more quickly to the selection pressures of a cold climate than others. However, Australian populations have been shown to differ across multiple chill tolerance traits 26-28 -and the species only reached that continent within the past few hundred years 37 . Finally, populations that do not appear to differ based on multiple chill tolerance traits may differ in their plastic responses to cold. Adult and developmental cold acclimation improve multiple chill tolerance traits (e.g. reducing chill coma onset temperatures, increasing survival after cold shock, etc.) 23 . Plastic responses can be very strong and may be more ecologically important than differences in basal (i.e. non-plastic) chill tolerance 27 . We note that the one trait that distinguished our populations, survival after 4 days at 4 °C, was one that may have permitted plasticity (specifically acclimation) during the assay. When held at temperatures slightly below the chill coma onset temperature, some Drosophila species, including D. melanogaster, show evidence of acclimation by recovering the ability to stand shortly after falling into a chill coma 38 or displaying improved cold shock resistance after exposure to 4 °C for several days 39 . Differences in chill tolerance plasticity have been studied on an intraspecific level before: Australian D. melanogaster populations, which have been shown to differ in multiple chill tolerance traits along a latitudinal gradient, do not differ in their response to either adult or developmental acclimation treatments 40 . However, a much larger study of D. melanogaster populations from around the world, including Africa and Europe, found that flies from temperate regions do indeed respond more strongly to a developmental acclimation treatment (i.e. their chill tolerance improved more) than flies from tropical regions 29 . That study measured only one trait, CCRT, and also found that flies from colder regions tended to have lower CCRTs than flies from warmer regions-a pattern that we did not observe in our study populations-but it demonstrates that intraspecific variation in the capacity for acclimation is possible.
It remains to be shown whether plasticity is behind the discrepancies that we observed. We saw essentially no differentiation in two commonly-and rapidly-measured traits (CCO and CCRT), and clear differentiation in one trait (survival after 4 days at 4 °C) measured under conditions that would normally induce some form of acclimation. Nevertheless, we feel that this system, which contains populations that adapted separately to cold environments in Africa and Europe, would lend itself well to future studies of how different forms of plasticity affect different thermal tolerance traits. Within a single population, genetic correlations of CT min can depend on rearing temperature 41 , so adaptive variation in chill tolerance traits, or correlations among traits, within or among species could appear depending on rearing conditions. Such an approach may also offer an opportunity to explore potential trade-offs between basal and plastic cold tolerance, which have previously been explored in drosophilids 42,43 (but see also van Heerwaarden and Kellermann 44 , who call previous work into question and discuss the challenges of accurately identifying basal vs. plastic thermal tolerance trade-offs).
An important caveat here is that we worked with laboratory-reared flies. However, past research has shown that chill tolerance traits are remarkably stable in lab reared flies over time 29,45,46 , and since we were able to confirm that the large difference in chill tolerance previously observed by Pool et al. (2017) between the French and Egyptian populations still exists. Therefore, we are confident that differences in chill tolerance seen in these populations reflect differences in their wild progenitors and are worthy of further study.

Conclusion
A cold environment can exert selection pressure on an insect and impact diverse chill tolerance traits in unexpected ways. To our surprise, traits that would allow higher activity in the cold (i.e. a low chill coma onset temperature) or faster recovery from chill coma have apparently not been selected for in the cold-climate African and European D. melanogaster populations tested here-but the ability to survive a chronic, mild cold stress has. These findings support the notion that chill tolerance traits are mediated by different physiological mechanisms and highlight the importance of measuring more than one chill tolerance trait in wild-derived populations. A benefit of such paired population study systems is that they can allow us to see how chill tolerance and its underlying mechanisms can evolve under similar low temperature selection pressures.

Methods
Study system. The study system used here (supplied and previously studied by Pool et al. (2017)) consists of seven Drosophila melanogaster populations: three pairs of closely-related populations plus one outgroup (Zambia) ( Table 1). Within each population pair, one population is derived from flies collected in a relatively warm location (e.g. Ethiopian lowlands), while the other is derived from flies collected in a relatively cold location (e.g. Ethiopian highlands).
To compare populations' cold tolerance, Pool et al. (2017) exposed flies to 4 °C for 96 h (4 days), returned them to room temperature, and then recorded the number that regained the ability to move or stand within 30 min. Within each population pair in that study, the population from the colder climate was significantly more cold tolerant than its warm-climate counterpart. This was especially the case for the France/Egypt pair 34 .
Fly rearing. Flies were reared at 25 °C on Bloomington medium with a 12 h:12 h light:dark cycle. Larval density can impact cold tolerance 47 , so when flies were to be used in experiments, egg density was standardised to approximately 50 per vial. Sexually mature females (presumed to have mated) were isolated under CO 2 on the day after emergence. CO 2 exposure was brief (less than 10 min) and flies were then given a minimum of 48 h to recover 48 . All flies used in experiments were females between 5 and 7 days post-ecdysis.

CCO.
A fly's chill coma onset temperature (CCO) is the temperature at which it becomes completely immobilized 6,49 . To measure this, we used a temperature ramping assay 50,51 . Briefly, the flies (eight to ten females www.nature.com/scientificreports/ per line) were placed individually in 3.7 mL glass screw-top bottles, which were clipped to a frame in a custom aquarium filled with a mix of ethylene glycol and water (1:1). The aquarium was connected to a 28 L programmable cooling bath (Model AP28R-30, VWR International, Radnor, USA) filled with the same ethylene glycol/ water mixture. For the first 15 min of the experiment, the flies were held at their rearing temperature (25 °C), and then the temperature was ramped down at 0.1 °C/min. The temperature inside the aquarium was monitored with three type K thermocouples connected to a TC-08 recorder (Pico Technology, St Neots, UK), and the average of the three thermocouples was used to estimate fly body temperature. To test whether flies were still capable of movement at temperatures approaching the chill coma onset temperature, the vials and metal frame were periodically tapped with a metal rod, and the CCO was recorded as the temperature when this failed to stimulate any response.
CCRT . Chill coma recovery time (CCRT) was measured by holding flies individually in 3.7 mL glass screwtop vials at 0 °C (submerged in an ice-water slurry) for six hours, returning them to room temperature, and then measuring the time that it took for each one to stand on all six legs 38,51 . Room temperature was recorded and accounted for in subsequent analysis. We measured eight to ten flies in most (26) lines, but in six, we measured between eleven and fifteen. In order to ensure a similar number of replicates in all lines, we randomly sampled ten replicates from each of those six lines (using the sample function in R) for the statistical analysis.
Survival. Because of the large number of flies needed to accurately measure survival and the number of lines in our study system, and because this trait has already been characterised across a large number of lines in this system 34 , only two populations (one pair) were used in the survival experiment: France (5 lines total) and Egypt (4 lines). Flies were held in groups of 10 in vials (30 mL; containing c. 7 mL of medium) that were kept horizontal www.nature.com/scientificreports/ to reduce the risk of flies getting stuck in the food during exposure or recovery. The groups of flies were held at 4 °C for zero to four days (4 to 9 replicates per line per day, where each replicate is a vial with 10 flies), given 24 h to recover at 25 °C, and then evaluated for survival using a 5 point scale (0: unmoving; 1: twitching but not standing; 3: standing but not walking; 4: walking but not climbing/flying; 5: climbing and/or flying) modified from MacMillan et al. 52 .
Statistics. All analyses were performed in R version 3.6.3 53 . To examine the effect of population (e.g. the population derived from flies collected in the South African highlands) on the chill coma onset temperature (CCO), we fit a general linear model with CCO as the response variable and line nested within population as predictor variables. In all analyses, normality of the residuals was confirmed through visual inspection. We then used Tukey's HSD (Honestly Significant Difference) test for pairwise comparisons among populations. Likewise, to examine the effect of population on chill coma recovery time (CCRT), we fit a general linear model with the log-transformed CCRT as the response variable (log transformation was necessary to ensure homogeneity of variance) with room temperature and line nested within population as predictor variables. Flies that did not recover during the 90 min observation time were given a CCRT of 90 min (17 out of 311 flies). Excluding these animals from the dataset entirely would artificially lower the mean, and we confirmed that doing so has no effect on our conclusions. Tukey's HSD test was then used for pairwise comparisons among populations.
To test for a correlation between CCO and CCRT, we calculated the mean CCO and CCRT of each line and fit a general linear model with mean CCO as the response variable, and population and mean CCRT as predictor variables. Because room temperature may affect CCRT, we also ran this analysis after adjusting CCRTs for the effect of room temperature extracted from the general linear model above. This adjustment, however, had no effect on our conclusion, so we opted to show the uncorrected data.
Mann-Whitney tests were used to compare population survival scores (which were not normally distributed) on each day of the chronic cold survival assay.

Data availability
All data is provided as a supplementary file for review and the same file will be included as supplementary material should the manuscript be accepted for publication.