Am. J. Trop. Med. Hyg., 74(5), 2006, pp. 786-794
Copyright © 2006 by The American Society of Tropical Medicine and Hygiene
PHENOTYPIC PLASTICITY AND GEOGRAPHIC VARIATION IN THERMAL TOLERANCE AND WATER LOSS OF THE TSETSE GLOSSINA PALLIDIPES (DIPTERA: GLOSSINIDAE): IMPLICATIONS FOR DISTRIBUTION MODELLING
JOHN S. TERBLANCHE*,
C. JACO KLOK,
ELLIOT S. KRAFSUR, AND
STEVEN L. CHOWN
Centre for Invasion Biology, Department of Botany and Zoology, Stellenbosch University, Stellenbosch, South Africa; Department of Entomology, Iowa State University, Ames, Iowa
 |
ABSTRACT
|
Using the tsetse, Glossina pallidipes, we show that physiologic plasticity (resulting from temperature acclimation) accounts for among-population variation in thermal tolerance and water loss rates. Critical thermal minimum (CTMin) was highly variable among populations, seasons, and acclimation treatments, and the full range of variation was 9.3°C (maximum value = 3.1 x minimum). Water loss rate showed similar variation (max = 3.7 x min). In contrast, critical thermal maxima (CTMax) varied least among populations, seasons, and acclimation treatments, and the full range of variation was only approximately 1°C. Most of the variation among the four field populations could be accounted for by phenotypic plasticity, which in the case of CTMin, develops within 5 days of temperature exposure and is lost rapidly on return to the original conditions. Limited variation in CTMax supports bioclimatic models that suggest tsetse are likely to show range contraction with warming from climate change.
 |
INTRODUCTION
|
Climate change and its interaction with human land use patterns are among the most significant threats facing biodiversity and society.1 Therefore, forecasts of the interactive effects of climate change and land use are of considerable importance for predicting the future distribution and abundance of species and for enabling society to mitigate and cope with the effects of these changes.1 Typically, predictions of climate change effects on species ranges are based on bioclimatic models,2,3 although there is substantial variety in the way in which this modeling is done. These models are often able to explain much of the variation in species distributions and abundances.47 They are therefore considered one of the most effective ways to obtain a rapid estimate of likely alterations in species distributions.2 Nonetheless, they have been criticized on a variety of grounds.8
From a physiologic perspective, the bioclimatic modeling approach raises concerns from three principal perspectives. Spatial variation in population responses to the environment is often not considered9; the rapid alterations to phenotypes that might take place through phenotypic plasticity in the form of developmental plasticity, acclimation, and hardening are typically ignored10; and the likely outcome of covariation among abiotic variables is often not adequately assessed.11 In consequence, it has been proposed that physiologic studies and biophysical modeling should be used in concert with large scale bioclimatic studies of species responses to understand what the future might hold for various taxa under a changing climate.10
These concerns are equally, if not more, pertinent to understanding changes in the distributions of arthropod-borne diseases.12 Nonetheless, it has been shown that climate change is likely to alter, or already has altered, the geographic distribution of at least some of these diseases,12,13 although in other cases, the effects are less straightforward to discern.11,14 Given that the distributions of vector species are important in influencing changing patterns of disease transmission,12,13 considerable attention has been given to the likely effects of environmental change on vectors. For Africa, much of the attention has focused on malaria and its mosquito vectors.14 However, tsetse (Glossina spp. Diptera: Glossinidae), which are long-lived, hematophagous flies, with low reproductive rates, have also been the subject of consideration because they carry trypanosomes that have serious implications for human and animal health and are considered a major limitation to the socio-economic development of the continent.7,15
Several studies have predicted that the ranges of African trypanosomiases and their tsetse vectors are likely to change under forecast climate scenarios.1618 These predictions arise from the fact that abiotic variables that will be affected by climate change, such as land surface temperature and saturation deficit, and other integrated environmental variables (e.g., the normalized difference vegetation index [NDVI]), are excellent correlates of tsetse distributions and abundances,47 explaining much of their variation. In some cases, the mechanisms underlying these relationships are likely to be indirect.7,4,19 However, the direct role of abiotic environmental variables in affecting tsetse life history parameters is well established.20 Temperature has major effects on birth rates, development, and mortality, and especially strong relationships exist between adult mortality and ambient temperature, which in some species are markedly non-linear.20,21 Water availability is also known to affect mortality rates; puparia and teneral adults are most sensitive to dry conditions, whereas mature adults are most sensitive to temperature.20
Despite the foregoing information on how responses to climate might vary among different tsetse species and subspecies, based both on remote sensing work5,6 and laboratory-based physiologic studies,22 inter-population physiologic variability and the rate of development and extent of phenotypic plasticity in tsetse are poorly understood (although some behavioral variation has been described23). Nonetheless, it is clear that if the responses of these species to climate change are to be fully comprehended, such information, combined with estimates of gene flow24,25 and dispersal rates,20,25,26 is essential. We therefore undertook such a study by investigating between-population variability in seasonal acclimatization and the extent and time-course of phenotypic plasticity of thermal tolerance and desiccation resistance in G. pallidipes Austen.
 |
MATERIALS AND METHODS
|
Field sampling and study populations.
Adult flies were sampled from four locations representing different thermal habitats in Kenya, East Africa (Table 1
). To determine seasonal variation in physiology, the Nguruman District (mid-altitude) and Lambwe District (high-altitude) populations were sampled twice within the same year in mid-July to mid-August 2003 (the end of the long rainy season) and in mid-November to mid-December 2003 (the end of the short rainy season). The Narok (highest altitude) and Kwale Districts (low-altitude) were sampled in November-December 2003. For each of the field experiments, flies were collected from odor-baited, biconical traps15 set out in the respective localities (key attractive components: 4-methyl-phenol, 3-n-propynol, and acetone). After removal from the traps, the flies were transported in an insulated container, lined with moist paper towel, to field laboratories. These were located a maximum of 2 hours away from each of the sampling sites of a given district. At the laboratories, the flies were separated into groups for each of the experimental procedures, after which experiments started immediately. Depending on the time of the experiment, flies were collected either in the morning or in the late afternoon, usually from 0800 to 1000 and from 1700 to 1800 hours. Thermal limit assessments were performed during the day using flies sampled in the morning, whereas desiccation experiments took place during the night using afternoon-collected flies.
View this table:
[in this window]
[in a new window]
|
TABLE 1 Field sampling sites (coordinates in decimal degreees) and time of sampling of G. pallidipes for this study
|
|
Critical thermal limits.
An insulated system of 11 double-jacketed isolation chambers was connected to a programmable water bath (LTD 20 with PZ1 programmer; Grant Instruments, Cambridge, UK) that regulated water flow around the chambers. A single fly was placed into each chamber, and a 40 SWG type-T thermocouple was inserted into the control chamber to measure chamber temperature. The flies were allowed to equilibrate for 10 minutes at either 16°C or 35°C before either minimum or maximum critical thermal limit assessments started, respectively.27 After equilibration, chamber temperature was lowered or raised at 0.25°C/min. Preliminary experiments confirmed that body temperature did not lag behind or exceed chamber temperature at this rate of temperature change. Critical thermal minimum (CTMin) was defined as the loss of coordinated muscle function, and critical thermal maximum (CTMax) was defined as the onset of muscle spasms. The temperature at which either of these reactions occurred was recorded for each individual. These endpoints are readily identifiable for any species once an observer is practiced,27 variance about the endpoints is typically low, and here a single observer (C.J.K.) undertook all of this work. Moreover, this experimental procedure has been verified using thermolimit respirometry,28 is widely used to assess thermal limits,29 and the observer typically was not informed which acclimation treatment was being assessed. The procedure was repeated three times for both CTMin and CTMax to give a total N = 30 per population and trait. Preliminary, replicated experiments found no effect of sex, age, or feeding status on critical thermal limits (C. J. Klok and J. S. Terblanche, unpublished data). Critical thermal maxima were lethal, whereas minima were not. However, even in the case of CTMin, all experimental animals were discarded after a trial.
Water loss rates.
Male flies (N = 16), held individually in 5-mL cuvettes, were subjected to desiccation in flowing air (< 2.5% relative humidity [RH]) for 10 hours at 24.0 ± 1.0°C (PTC-1 Peltier-controlled temperature cabinet; Sable Systems, Las Vegas, NV).30 Air was pumped through a scrubbing column containing silica gel and Drierite to remove residual water and into a mass flow controller (MFC). The MFC outlet was connected to a Sable Systems MF8 airflow manifold. Each outflow channel of the manifold was further split into two streams so that two cuvettes, each containing a fly, were attached to a single manifold channel. The air flow rate through each cuvette, tested with a second MFC, was regulated to produce a rate of 100 mL/min. Experiments took place during the night (21000700 hours), a period of minimal activity in tsetse, and were repeated twice per location (N = 16 per experiment), each time on individuals collected in the field on that particular day, to give a total N = 32. Each fly was weighed before and after a trial (Avery Berkel FA 304T, Fairmont, MN, electronic microbalance, 0.1-mg resolution), and water loss rate was expressed as milligrams of water per hour. Excretion usually takes place after a blood meal, and here, experiments took place well after such events. In addition, even in field-collected flies, where the time since feeding could not be accurately determined, no excretion in the cuvettes was observed. Therefore, mass loss was considered representative of water loss even though a small amount of mass change took place as a consequence of metabolism.22 After the trials, flies were dried to constant mass (~5060°C for ~72 hours) and were re-weighed to determine water content. Solvent-based lipid extractions were used to determine lipid mass. Flies used in desiccation trials were dismembered and soaked in 2 mL chloroform:methanol (2:1) solution for approximately 48 hours to dissolve internal (body) lipids. Lipid content was estimated by determining dry body mass before and after lipid extraction.31
Activity levels.
Several authors have suggested that differences in desiccation resistance among taxa might simply reflect divergence in activity levels.22,32 Because little is known about inter-population variation in tsetse activity levels, particularly under stressful conditions, the proportion of time individuals were active was estimated over a 30-min period for each population. Time active at 20°C, 24°C, 28°C, and 32°C was determined using a single electronic activity detector (AD-1; Sable Systems) that exported voltage data to a computer through Sable Systems DATACAN V software. These experiments were performed during the same field sampling periods reported above. Voltage data were sliced from the raw Datacan V (Sable Systems) data files, copied to MS-Excel, converted to absolute values, and sorted from lowest to highest. Based on preliminary observations, we considered activity as any signal more than 0.01 V measured by the AD-1 (electronic spikes of 5 V were ignored and can easily be distinguished from true activity readings). The total time spent above this level of activity was expressed as a proportion of the total recording time. While it is well known that activity in tsetse is sensitive to feeding state,33 the objective here was to determine whether differences in activity level were a likely cause of any differences in water loss rates between populations and not whether the flies activity overall differs between populations.
Within-generation physiologic plasticity.
Logistic constraints precluded the use of field populations for investigations of phenotypic plasticity. Therefore, G. pallidipes from a laboratory colony maintained at the Entomology Unit, FAO/IAEA Agriculture and Biotechnology Laboratory, Seibersdorf, International Atomic Energy Agency, Vienna, Austria, was used. The IAEA laboratory colony used in this study was established in 1975 from a natural population in Tororo, Uganda, close to the Kenyan border (and < 200 km from our Lambwe field site). Gene diversities over two mitochondrial loci were within the range of 18 field populations from Ethiopia and East and Southern Africa.34 Therefore, the laboratory-reared flies may be considered genetically representative of flies from field populations. Puparia were immediately placed inside plastic containers and transferred to a climate chamber set to 24°C (24.6 ± 2.5°C, 12:12-hour L:D) and 76% RH as described previously.35 Feeding of the newly emerged adults took place using a membrane-tray system36 every alternative day, and at this time, container locations were randomized within the climate chamber. Care was taken to ensure that all treatment groups were handled for the same duration in transferal from the climate chamber to the feeding area and spent a similar amount of time outside the climate chamber while feeding (~25 minutes on each feeding day). Flies were held at 19°C, 24°C, and 29°C for 10 days in climate chambers with a synchronized photoperiod and constant 76% RH levels. The 24°C and 29°C groups were stored in temperature-regulated climate chambers (Labcon, Pretoria, South Africa; mean ± SD: 24.0 ± 2.3 [N = 1029]; 28.9 ± 1.7°C [N = 937]) and the 19°C group was stored in a Peltier-controlled temperature cabinet (mean ± SD: 19.4 ± 0.4°C [N = 905]). Acclimation temperatures were chosen to represent the mean annual temperatures of the respective geographic locations sampled (reported mean annual temperaturewarm-temperate midlands [< 1,850 m]: 1820°C; warm lowlands [9001,200 m]: 2224°C; hot lowlands [0900 m]: 2430°C).37 All experimental and rearing facilities were contained within a quarantine approved, air-conditioned laboratory. The same system and methods as those used in the field experiments for determining geographic variation were used in the laboratory determination of critical thermal limits and desiccation rate.
A separate batch of colony-reared flies was used to determine whether the temperature acclimation effect observed in the CTMin of G. pallidipes is time- and/or temperature-dependent and whether these changes in CTMin are reversible. Flies were allowed to emerge in our laboratory as described above and were fed twice within 4 days after eclosion. However, a small number of flies seemed only to feed during the second feeding round. Thus, typically after two blood meals, flies were randomly assigned to four temperature acclimation groups (24°C, 21°C, 19°C, and 16°C; N = 180 ~ 220 per group, five to six cages per group). Humidity and photo-period were regulated as described above, and randomization of cage locations within chambers followed the protocol adopted previously. Therefore, handling times were kept constant across all groups during feeding sessions and experimental treatments. The CTMin of 10 randomly selected flies from each treatment group was determined on acclimation days 1, 3, 5, 7, and 9. After 9 days of acclimation, the remaining flies were all returned to 24°C for an additional 9-day period, after which another sample of flies was assessed on the tenth day. In this experiment, once a fly had been used for determination of CTMin, it was discarded by bottling in 90% ethanol (i.e., no animals were re-used). Virgin females, separated from males into cages immediately after eclosion, were used to study effects of temperature acclimation on body mass. Feeding took place using the same methods described above only on days between experimental studies (i.e., flies were always in a fasted state during experiments).
Statistical analyses.
Data were assessed for normality using a Shapiro-Wilk test before analyses. Because critical thermal limit data were normally distributed, analyses of variance (ANOVAs) using the raw data were used to compare critical thermal limits for both field and laboratory studies. Analyses of covariance (ANCOVAs) were used to test for differences in rates of water loss and lipid content between populations and seasons (i.e., mass-independent water loss rates), using mean experimental or dry mass as the covariate. Pearson product-moment correlations were used to assess the relationships between time spent active and temperature. Data are presented as means ± SE unless otherwise stated. Significance was set at P = 0.05. In no instances for any trait were experimental assessments performed on the same individual on more than one occasion (i.e., animals were not re-used).
 |
RESULTS
|
Seasonal variation.
In the only populations that were sampled on two occasions, season had a significant effect on the majority of the variables examined (Table 2
). In the Lambwe and Nguruman populations, CTMin increased and water loss rate declined after the short rains. Although lipid mass also declined in the Lambwe population, the change was not significant in flies from Nguruman. Seasonal changes in body mass and CTMax were opposite in sign in the two populations, although the seasonal difference in CTMax was only 0.4°C in both populations. Body mass increased in the Lambwe population but declined at Nguruman, leading to an overall difference in mass between the two populations after the short rains (Table 2
; F1,91 = 20.35; P < 0.0001) but not after the long rains (Table 2
; F1,126 = 0.414; P = 0.52).
View this table:
[in this window]
[in a new window]
|
TABLE 2 Summary statistics for seasonal variation of pysiological traits in G. pallidipes from the Lambwe and Nguruman populations
|
|
Geographic variation.
Using data collected within the same season after the short rains, CTMax and CTMin differed among the four populations (Table 3
). The mid-altitude and low-altitude populations differed only marginally, whereas both CTMin and CTMax declined with a further increase in elevation. The slope of the two altitude-critical limit relationships also differed significantly (ANCOVA, F1,234 = 59.42; P < 0.0001), with CTMin showing a steeper slope and stronger relationship with altitude than CTMax. Overall, CTMin varied by 3.3° C among the four populations, whereas CTMax differed by only 0.4°C.
View this table:
[in this window]
[in a new window]
|
TABLE 3 Geographic variation in mean population critical thermal maxima (CTMax) and minima (CTMin) for G. pallidipes
|
|
Water loss rate also varied significantly among the populations (Figure 1
). It seems unlikely that differences in activity levels among the populations could account for the variation in water loss rate. Total activity time was positively, although in one case weakly, correlated with temperature across the range of 2032°C in all four populations (Narok: r = 0.64, P < 0.01; Lambwe: r = 0.47, P = 0.068; Nguruman: r = 0.64, P < 0.01; Kwale: r = 0.50, P < 0.05). Because there were no sex-related differences in activity time in any of the populations, sexes were pooled for the inter-population comparisons. There were no among-population differences for time spent active (ANCOVA [covariate: temperature]: F3,59 = 1.38; P = 0.26). Therefore, differences in activity probably did not account for differences in water loss rates.

View larger version (11K):
[in this window]
[in a new window]
|
FIGURE 1. Geographic variation of water loss rate in G. pallidipes. Note that desiccation rate is highest in the low altitude population from Kwale (388 m). All flies were sampled in the same season (November to December 2003; Narok: 1,691 m; Lambwe: 1,353 m; Nguruman: 670 m).
|
|
Temperature acclimation of laboratory-bred flies.
Ten days of laboratory acclimation at three different temperatures resulted in significant among-treatment variation in CTMax, CTMin, desiccation rate, and body water content (Table 4
). increased significantly with treatment temperature, CTMax but the change was less than 0.5°C. In contrast, at 19°C, declined by approximately 3°C relative to the other CTMin treatments. Water loss rate showed similar threshold effects, declining significantly at the highest acclimation temperature, whereas body water content showed little change with acclimation temperature (Table 4
).
View this table:
[in this window]
[in a new window]
|
TABLE 4 Upper (CTMax) and lower (CTMin) critical thermal limits, water loss rate at 24°C and initial body water content after 10 days of temperature acclimation in laboratory-reared G. pallidipes adults
|
|
In the time-course experiments, changes in CTMin were induced in the 16°C and 19°C treatments only, and the magnitude of the induced change was quite similar (Figure 2
). In contrast, the 21°C and 24°C treatments had no effect. In the former treatments, CTMin had reached its nadir by the fifth day of treatment and typically remained at a similar level for the remainder of the acclimation period (day 9). Ten days at 24°C were sufficient to completely reverse the effects of the acclimation treatments (i.e., at day 19, no groups differed from the 24°C group). Thus, significant day x treatment interactions were found for CTMin (F15,216 = 3.99; P < 0.0001). The low treatment temperatures clearly came at some cost to normal adult development. Whereas flies gained mass in the 21°C and 24°C treatments and showed low and typically consistent among-individual variation, the flies at the lower temperatures did not show any steady progression in mass, and variation among individuals was considerable (Figure 3
). The differential effect of treatment on change in mass over the 19 days is evidenced by a significant day x treatment interaction in a two-way ANOVA examining the effects of day and treatment on mass (F15,120 = 2.7, P = 0.0013).

View larger version (21K):
[in this window]
[in a new window]
|
FIGURE 2. The time-course of acclimation for CTMin in G. pallidipes adults. After 9 days of acclimation to 16°C, 19°C, 21°C, and 24°C, all acclimation groups were returned to 24°C for 10 days of recovery (tested on day 19). (Recorded climate chamber temperatures [mean ± SD]: 16: 15.7 ± 0.37°C; 19: 18.6 ± 0.22°C; 21: 20.6 ± 0.50°C; 24: 24.4 ± 1.38°C). N = 10 flies for all days and all acclimation temperature groups. There was no effect of acclimation on 21°C and 24°C acclimation groups, but both 16°C and 19°C had significantly reduced CTMin by day 5 (ANOVA: 16°C: F5,54 = 8.70; P < 0.0001; 19°C: F5,54 = 9.31; P < 0.0001; 21°C: F5,54 = 1.99; P = 0.095; 24°C: F5,54 = 1.38; P = 0.25). By day 19, all acclimation groups were similar (F3,36 = 2.64; P = 0.064).
|
|

View larger version (16K):
[in this window]
[in a new window]
|
FIGURE 3. Body mass (grams; mean ± SE) of virgin female G. pallidipes during the course of acclimation to (A) 24°C, (B) 21°C, (C) 19°C, and (D) 16°C. All flies were returned to 25°C after day 9 of acclimation. All mass measurements were obtained on non-feeding days (i.e., 24 h after feeding). Although no differences were found in body mass on completion of the acclimation treatments, there was a significant day x treatment interaction (see text for statistics).
|
|
Colony versus field flies.
Upper and lower critical thermal limits differed among the laboratory-acclimated (colony) flies and field populations (Figure 4
, A and B). However, the extent of this difference was only approximately 1°C for CTMax but from 4°C to 10°C for CTMin. Water loss rates also varied among the field and laboratory populations, and here, the maximum difference among groups (Kwale and the 29°C treatment group) was substantial, although none of the other groups differed from each other (Figure 4C
). In contrast, dry mass varied little among the populations and groups (Figure 4D
).

View larger version (24K):
[in this window]
[in a new window]
|
FIGURE 4. Traits means (±SE) in field and laboratory populations of G. pallidipes of (A) critical thermal maxima (ANOVA, F6,161 = 22.38; P < 0.0001); (B) critical thermal minima (ANOVA, F6,195 = 453.20; P < 0.0001); (C) water loss rate (ANCOVA, covariate: dry body mass; F6,126 = 12.378; P < 0.0001; least squares adjusted means ± 95% confidence limits); and (D) body size as dry body mass from all desiccation experiment flies (ANOVA; F6,239 = 8.43; P < 0.0001; means ± 95% confidence limits). The numerals in Lab19, Lab24, and Lab29 refer to the temperature of acclimation in laboratory-reared flies for periods of 910 days.
|
|
 |
DISCUSSION
|
Physiologic implications.
The temperature and water balance traits differed substantially from each other in the extent of their variability. Of the three traits, varied least CTMax among populations (0.6°C), seasons (0.4°C), and acclimation treatments (0.5°C), and the full range of variation among all of these samples was only 1.1°C (maximum value, 1.02 times larger than the minimum). In contrast, CTMin was more variable, differing by 3.3°C, 12°C, and 3.3°C among the populations, seasons, and acclimation treatments, respectively, and the full range of variation was 9.3°C (maximum value = 3.1 x minimum). Water loss rate showed similar, substantial variation, with the overall maximum value, across laboratory and field populations, being 3.7 times the minimum. These results are consistent with work across a range of insect taxa. In insects, upper critical temperatures show much less variation than lower critical temperatures among populations and acclimation treatments,38 between selection regimens,39 and among species at a variety of spatial scales.40 Substantial differences in water loss rates among populations, among seasons within a population, and between humidity acclimation treatments are also common.38,41,42 However, only a few studies have examined temperature acclimation or acclimatization effects on desiccation resistance, with some finding notable effects of temperature,38,43 and others finding no influence.30 Thus, upper thermal critical limits that are not especially variable, and water loss rates and lower thermal limits that are more flexible are typical of insects in general. In consequence, their occurrence in G. pallidipes is not entirely unexpected.44
A further significant finding is that the magnitude of the among-population differences was identical to or smaller than the magnitude of the acclimation effects in the laboratory-reared flies. Phenotypic plasticity might therefore account for all of the among-population variation in G. pallidipes. However, because of logistic constraints, each population was not examined to determine the extent of its phenotypic plasticity. Therefore, the extent to which among-population variation really is accounted for by phenotypic plasticity could not be fully determined. Nonetheless, seasonal acclimatization of CTMin, CTMax, and water loss rates in the Lambwe and Nguruman populations encompassed 50100% of the among-population variation in these traits. Together with the results from the laboratory acclimation trials, these findings provide substantive evidence for the idea that most of the variation among populations could be accounted for by phenotypic rather than genetic differences. Similar results have been obtained for two independent assessments of Drosophila melanogaster.38,45 Thus, it seems likely that much of the among-population variation in these traits examined in G. pallidipes is a consequence of phenotypic plasticity.
In the case of CTMin these plastic changes take place rapidly. Within 5 days of exposure to a different temperature, flies held at the lower temperatures had substantially altered their CTMin, and by 9 days after their return to a common, higher temperature environment, CTMin had returned to its original value. Whether the rapid phenotypic change found in G. pallidipes is typical of CTMin or similar traits in other insect species is not clear, largely because the time-course of thermal acclimation (rather than developmental plasticity) has not been thoroughly studied in the group. Nonetheless, several studies have documented similar rapid changes in physiologic traits.4648 Therefore, reasonably rapid acclimation effects are probably not unusual in insects.
Finally, whereas CTMax, water loss rate, and dry mass were quite similar across all of the populations and acclimation treatments, CTMin differed substantially between the field-collected and laboratory-reared flies (Figure 4
). The latter had CTMins that were 39°C lower than those of the field-collected flies. In consequence, it seems likely that differences in the phenotypic plasticity of the upper and lower critical thermal limits extend to the ease with which these traits can evolve. Higher acclimation temperatures in laboratory flies, or lower acclimation temperatures in field flies, might have further altered CTMin. However, this seems unlikely given that more extreme temperatures are debilitating to flies over the longer term.7,20 Therefore, the differences between the field and laboratory flies might have evolved as a consequence of conditions in the laboratory. The IAEA laboratory colony used in this study was established in 1975. Assuming linear evolutionary change, approximately 175 generations in the laboratory, and CTMin values of the original Ugandan field population similar to those found at Lambwe, this would equate to an approximately 1°C change every 29 generations (i.e., 6°C/175 generations). While the reasons for laboratory evolution are not always clear, substantial changes in animals held in colonies have been documented in several studies.49,50 Whatever the reason for the evolved difference, the most noteworthy feature thereof is that this evolution has been restricted to CTMin and has not taken place in the other traits. In consequence, it seems likely that most of the variation in the other traits is a consequence of phenotypic plasticity, whereas the variation found in CTMin is a mixture of plasticity and evolved differences between laboratory and field-collected G. pallidipes. In addition, it raises questions about the extent to which the laboratory flies are genetically representative of field-collected flies, at least in the case of this trait. Future work exploring the underlying genetic basis of CTMin would be required to resolve this question.51
Implications for distributional modeling.
Although the critical thermal maxima of adult G. pallidipes documented here are higher than the temperatures that induce mortality in the field,20 these findings nonetheless substantiate correlative work on the determinants of the distribution of this species and provide several insights into the likely responses of tsetse to climate change. In the first instance, univariate and multivariate analysis of climate and productivity (NDVI) variables showed that temperature variables are much better predictors of G. pallidipes presence and absence than NDVI variables.5,6 These results are in keeping with population level analyses that show a strong relationship between increasing temperature and elevated mortality in G. pallidipes adults.20 Our results suggest a causal basis for the determination of G. pallidipes distributions by temperature. It seems that, in this species, tolerance to high temperatures cannot change much either by phenotypic plasticity or by evolution. Such change in tolerance would result in a weaker relationship between temperature and population dynamics, and in doing so, would reduce the strength of the relationships between temperature and the presence/absence of the fly species.52 The limited variation in CTMax also provides a causal explanation for why, in statistical models of tsetse distribution, the mean temperature difference between areas of fly presence and absence may be less than 1°C.19
Limited scope for change in upper critical temperatures might also explain the narrow variation in this trait among populations despite the fact that gene flow between G. pallidipes populations is typically low.25,34 From a theoretical perspective, limited gene flow between populations should promote diversification as long as gene flow is not so low that novel variants that have successfully survived similar conditions elsewhere are seldom introduced into a population.53,54 However, if the trait itself cannot evolve because it is globally constrained,55 changes in gene flow are unlikely to affect trait diversification.
From these data, it might seem difficult to discern why flies are not able to colonize lower temperature areas. However, lower temperature thresholds are likely set by development of the larvae rather than by adult capabilities.56 Moreover, in our acclimation experiments, we showed that, although flies are able to survive longer-term maintenance at low temperatures, their muscle growth and development (which accounts for the largest change in mass in adult tsetse57) is compromised, given that mass declines substantially by days 79, and variance in mass between individuals was larger in the lower than in the high temperature treatments (Figure 3
). Therefore, low temperatures are likely to compromise adults also.
To date, most climate change predictions suggest that tsetse species would show a reduction in the range of the flies owing to increases in temperature.17 Our data support these findings because they show limited plasticity in upper critical temperatures. If upper thermal tolerances are inflexible, increases in temperature are likely to lead to a reduction in range size as suggested by modeling approaches.7 What altered rainfall patterns are likely to mean is more difficult to discern because water loss rates in adult flies respond to altered ambient conditions,44 water balance effects are likely to be most significant for puparia,20,56 which were not investigated here, and fine-scale predictions of rainfall change are typically not available.1 Therefore, it is difficult to be able to predict what the outcome of the interaction between fly desiccation resistance and changing water availability might be. Nonetheless, our results provide further evidence in support of the idea that, owing to changes in land use and climate, at least some tsetse may be brought under autonomous control in parts of their natural ranges.58
Received October 18, 2005.
Accepted for publication January 18, 2006.
Acknowledgments: The authors thank Andrew Parker and the IAEA/FAO Laboratories for contributing laboratory-rearing equipment, puparia, and source colony information. The fieldwork in Kenya was supported by KETRI who aided us with staff and logistics. Russ Jurenka performed the lipid extractions at ISU. In Stellenbosch, Henry Davids and Larissa Heyns provided laboratory support. Sue Jackson, Allen Gibbs, Saskia Goldberg, and Brent Sinclair are thanked for comments and discussion. Two anonymous referees are thanked for comments.
Financial support: This work was funded by NIH Grant AI-52456 to E. S. Krafsur.
* Address correspondence to John S. Terblanche, Centre for Invasion Biology, Department of Botany and Zoology, Stellenbosch University, Private Bag X1, Matieland 7602, South Africa. E-mail: jst{at}sun.ac.za 
Authors addresses: John S. Terblanche and Steven L. Chown, Centre for Invasion Biology, Department of Botany and Zoology, Stellenbosch University, Private Bag X1, Matieland 7602, South Africa, E-mail: jst{at}sun.ac.za and slchown{at}sun.ac.za. C. Jaco Klok,School of Life Sciences, Box 874501, Arizona State University, Tempe, AZ 85287-4501, E-mail: cjklok{at}asu.edu. Elliot S. Krafsur, Department of Entomology, Iowa State University, Ames, IA 50011, E-mail: ekrafsur{at}iastate.edu.
 |
REFERENCES
|
- Watson RT, 2002. Climate Change 2001: Synthesis Report. Cambridge: Cambridge University Press.
- Thomas CD, Cameron A, Green RE, Bakkenes M, Beaumont LJ, Collingham YC, Erasmus BFN, de Siqueira MF, Grainger A, Hannah L, Hughes L, Huntley B, van Jaarsveld AS, Midgley GF, Miles L, Ortega-Huerta MA, Peterson AT, Phillips OL, Williams SE, 2004. Extinction risk from climate change. Nature 427: 145148.[Medline]
- Pearson RG, Dawson TP, 2003. Predicting the impacts of climate change on the distribution of species: are bioclimate envelope models useful? Global Ecol Biogeog 12: 361371.
- Rogers DJ, Randolph SE, 1991. Mortality rates and population density of tsetse flies correlated with satellite imagery. Nature 351: 739741.[Medline]
- Robinson T, Rogers D, Williams B, 1997. Univariate analysis of tsetse habitat in the common fly belt of Southern Africa using climate and remotely sensed vegetation data. Med Vet Entomol 11: 223234.[Medline]
- Robinson T, Rogers D, Williams B, 1997. Mapping tsetse habitat suitability in the common fly belt of Southern Africa using multivariate analysis of climate and remotely sensed vegetation data. Med Vet Entomol 11: 235245.[Web of Science][Medline]
- Rogers DJ, 2000. Satellites, space, time, and the African trypanosomiases. Adv Parasitol 47: 129171.[Web of Science][Medline]
- Davis AJ, Jenkinson LS, Lawton JH, Shorrocks B, Wood S, 1998. Making mistakes when predicting shifts in species range in response to global warming. Nature 391: 783786.
- Davis MB, Shaw RG, 2001. Range shifts and adaptive responses to Quaternary climate change. Science 292: 673679.[Web of Science]
- Helmuth B, Kingsolver JG, Carrington E, 2005. Biophysics, physiological ecology, and climate change: does mechanism matter? Ann Rev Physiol 67: 177201.[Web of Science][Medline]
- Rogers DJ, Randolph SE, 2000. The global spread of malaria in a future, warmer world. Science 289: 17631766.[Abstract/Free Full Text]
- Kovats RS, Campbell-Lendrum DH, McMichael AJ, Woodward A, Cox J, St H, 2001. Early effects of climate change: do they include changes in vector-borne disease? Phil Trans R Soc Lond B 356: 10571068.[Abstract/Free Full Text]
- Burse BV, Mellor PS, Rogers DJ, Samuel AR, Mertens PPC, Baylis M, 2005. Climate change and the recent emergence of bluetongue in Europe. Nature Rev Microbiol 3: 171181.
- Rogers DJ, Randolph SE, Snow RW, Hay SI, 2002. Satellite imagery in the study and forecast of malaria. Nature 415: 710715.[Medline]
- Leak SGA, 1999. Tsetse Biology and Ecology. Wallingford: CABI Publishing.
- Watson RT, Zinyowera MC, Moss RH, Dokken DJ, 1998. The Regional Impacts of Climate Change. An Assessment of Vulnerability. Cambridge: Cambridge University Press.
- Hulme M, 1996. Climate change in Southern Africa: An Exploration of Some Potential Impacts and Implications in the SADC Region. Norwich: Climatic Research Unit, University of East Anglia.
- World Health Organization, 2004. Using Climate to Predict Disease Outbreaks. Geneva, Switzerland: World Health Organization.
- Rogers DJ, Robinson TP, 2004. Tsetse distribution. In: Maudlin I, Holmes PH, Miles MA, eds. The Trypanosomiases. Wallingford: CABI Publishing; 139179.
- Hargrove JW, 2004. Tsetse population dynamics. In: Maudlin I, Holmes PH, Miles MA, eds. The Trypanosomiases. Wallingford: CABI Publishing; 113138.
- Hargrove JW, 2001. The effect of temperature and saturation deficit on mortality in populations of male Glossina m. morsitans (Diptera: Glossinidae) in Zimbabwe and Tanzania. Bull Ent Res 91: 7986.[Web of Science][Medline]
- Bursell E, 1959. The water balance of tsetse flies. Trans R Ent Soc Lond 111: 205235.
- Langley PA, Maudlin I, Leedham MP, 1984. Genetic and behavioural differences between Glossina pallidipes from Uganda and Zimbabwe. Entomol Exp Appl 35: 5560.
- Krafsur ES, 2003. Tsetse fly population genetics: an indirect approach to dispersal. Trends Parasitol 19: 162166.[Medline]
- Ouma JO, Marquez JG, Krafsur ES, 2005. Macrogeographic population structure of the tsetse fly, Glossina pallidipes (Diptera: Glossinidae). Bull Ent Res 95: 437447.[Medline]
- Williams B, Dransfield R, Brightwell R, 1992. The control of tsetse flies in relation to fly movement and trapping efficiency. J Appl Ecol 29: 163179.
- Lutterschmidt WI, Hutchison VH, 1997. The critical thermal maximum: history and critique. Can J Zool 75: 15611574.
- Klok CJ, Sinclair BJ, Chown SL, 2004. Upper thermal tolerance and oxygen limitation in terrestrial arthropods. J Exp Biol 207: 23612370.[Abstract/Free Full Text]
- Chown SL, Nicolson SW, 2004. Insect Physiological Ecology. Mechanisms and Patterns. Oxford: Oxford University Press.
- Terblanche JS, Sinclair BJ, Klok CJ, McFarlane ML, Chown SL, 2005. The effects of acclimation on thermal tolerance, desiccation resistance and metabolic rate in Chirodica chalcoptera (Coleoptera: Chrysomelidae). J Ins Physiol 51: 10131023.[Web of Science][Medline]
- Hoffmann AA, Hallas R, Anderson AR, Telonis-Scott M, 2005. Evidence for a robust sex-specific trade-off between cold resistance and starvation resistance in Drosophila melanogaster. J Evol Biol 18: 804810.[Medline]
- Gibbs AG, Chippindale AK, Rose MR, 1997. Physiological mechanisms of evolved desiccation resistance in Drosophila melanogaster. J Exp Biol 200: 18211832.[Abstract]
- Brady J, 1988. Circadian ontogeny in the tsetse fly: a permanent major phase change after the first feed. J Ins Physiol 34: 743749.[Web of Science]
- Krafsur ES, Wohlford DL, 1999. Breeding structure of Glossina pallidipes populations evaluated by mitochondrial variation. J Hered 90: 635642.[Abstract/Free Full Text]
- Terblanche JS, Klok CJ, Chown SL, 2004. Metabolic rate variation in Glossina pallidipes (Diptera: Glossinidae): gender, ageing and repeatability. J Ins Physiol 50: 419428.[Medline]
- Gooding RH, Feldmann U, Robinson AS, 1997. Care and maintenance of tsetse colonies. In: Crampton JM, Beard CB, Louis C, eds. The Molecular Biology of Insect Disease Vectors. London: Chapman & Hall; 4155.
- Sombroek WG, Braun HMH, van der Pouw BJA, 1980. Exploratory Soil Map and Agro-Climatic Zone Map of Kenya. Nairobi, Kenya: Ministry of Agriculture, National Agricultural Laboratories.
- Hoffmann AA, Shirriffs J, Scott M, 2005. Relative importance of plastic vs genetic factors in adaptive differentiation: geographical variation for stress resistance in Drosophila melanogaster from eastern Australia. Funct Ecol 19: 222227.
- Gilchrist GW, Huey RB, Partridge L, 1997. Thermal sensitivity of Drosophila melanogaster: evolutionary responses of adults and eggs to laboratory natural selection at different temperatures. Physiol Zool 70: 403414.[Medline]
- Chown SL, 2001. Physiological variation in insects: hierarchical levels and implications. J Ins Physiol 47: 649660.[Web of Science][Medline]
- Addo-Bediako A, Chown SL, Gaston KJ, 2001. Revisiting water loss in insects: a large scale view. J Ins Physiol 47: 13771388.[Web of Science][Medline]
- Gibbs AG, Matzkin LM, 2001. Evolution of water balance in the genus Drosophila. J Exp Biol 204: 23312338.[Abstract/Free Full Text]
- Parkash R, Tyagi P, Sharma I, Rajpurohit S, 2005. Adaptations to environmental stress in altitudinal populations of two Drosophila species. Physiol Entomol 30: 353361.
- Rogers DJ, 1990. A general model for tsetse populations. Ins Sci Appl 11: 331346.
- Ayrinhac A, Debat V, Gibert P, Kister AG, Legout H, Moreteau B, Vergilino R, David JR, 2004. Cold adaptation in geographical populations of Drosophila melanogaster: phenotypic plasticity is more important than genetic variability. Funct Ecol 18: 700706.
- Zeilstra I, Fischer K, 2005. Cold tolerance in relation to developmental and adult temperature in a butterfly. Physiol Entomol 30: 9295.
- Fischer K, Eenhoorn E, Bot ANM, Brakefield PM, Zwaan BJ, 2003. Cooler butterflies lay larger eggs: developmental plasticity vs acclimation. Proc R Soc Lond Biol Sci 270: 20512056.[Medline]
- Gibert P, Huey RB, Gilchrist GW, 2001. Locomotor performance of Drosophila melanogaster: interactions among developmental and adult temperatures, age, and geography. Evolution 55: 205209.[Web of Science][Medline]
- Sgrò CM, Partridge L, 2000. Evolutionary responses of the life history of wild-caught Drosophila melanogaster to two standard methods of laboratory culture. Am Nat 156: 341353.[Web of Science]
- Terblanche JS, Klok CJ, Marais E, Chown SL, 2004b. Metabolic rate in the whip-spider, Damon annulatipes (Arachnida: Amblypygi). J Ins Physiol 50: 637645.[Medline]
- Loeschcke V, Sorensen JG, Kristensen TN, 2004. Ecologically relevant stress resistance: from microarrays and quantitative trait loci to candidate genes - A research plan and preliminary results using Drosophila as a model organism and climatic and genetic stress as model stresses. J Biosci 29: 503511.[Medline]
- Hoffmann AA, Blows MW, 1994. Species borders: ecological and evolutionary perspectives. Trends Ecol Evol 9: 223227.
- Case TJ, Taper ML, 2000. Interspecific competition, environmental gradients, gene flow, and the coevolution of species borders. Am Nat 155: 583605.
- Forde SE, Thompson JN, Bohannan BJM, 2004. Adaptation varies through space and time in a coevolving host-parasitoid interaction. Nature 431: 841844.[Medline]
- Hoffmann AA, Hallas RJ, Dean JA, Schiffer M, 2003. Low potential for climatic stress adaptation in a rainforest Drosophila species. Science 301: 100102.[Abstract/Free Full Text]
- Langley PA, 1977. Physiology of tsetse flies (Glossina ssp.) (Diptera: Glossinidae): a review. Bull Ent Res 67: 523574.
- Bursell E, 1961. Post-teneral development of the thoracic musculature in tsetse flies. Proc R Ent Soc Lond A 36: 6974.
- Rogers DJ, Randolph SE, 2002. A response to the aim of eradicating tsetse from Africa. Trends Parasitol 18: 534536.[Medline]
This article has been cited by other articles:

|
 |

|
 |
 
E. Kleynhans and J. S Terblanche
The evolution of water balance in Glossina (Diptera: Glossinidae): correlations with climate
Biol Lett,
February 23, 2009;
5(1):
93 - 96.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
S. L Chown and K. J Gaston
Macrophysiology for a changing world
Proc R Soc B,
July 7, 2008;
275(1642):
1469 - 1478.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
P. Calosi, D. T Bilton, and J. I Spicer
Thermal tolerance, acclimatory capacity and vulnerability to global climate change
Biol Lett,
February 23, 2008;
4(1):
99 - 102.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
J. S Terblanche, J. A Deere, S. Clusella-Trullas, C. Janion, and S. L Chown
Critical thermal limits depend on methodological context
Proc R Soc B,
December 7, 2007;
274(1628):
2935 - 2943.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
S. L Chown, S. Slabber, M. A McGeoch, C. Janion, and H. P. Leinaas
Phenotypic plasticity mediates climate change responses among invasive and indigenous arthropods
Proc R Soc B,
October 22, 2007;
274(1625):
2531 - 2537.
[Abstract]
[Full Text]
[PDF]
|
 |
|