Thermal limits of wild and laboratory strains of two African malaria vector species, Anopheles arabiensis and Anopheles funestus
Malaria Journalvolume 11, Article number: 226 (2012)
Malaria affects large parts of the developing world and is responsible for almost 800,000 deaths annually. As climates change, concerns have arisen as to how this vector-borne disease will be impacted by changing rainfall patterns and warming temperatures. Despite the importance and controversy surrounding the impact of climate change on the potential spread of this disease, little information exists on the tolerances of several of the vector species themselves.
Using a ramping protocol (to assess critical thermal limits - CT) and plunge protocol (to assess lethal temperature limits - LT) information on the thermal tolerance of two of Africa’s important malaria vectors, Anopheles arabiensis and Anopheles funestus was collected. The effects of age, thermal acclimation treatment, sex and strain (laboratory versus wild adults) were investigated for CT determinations for each species. The effects of age and sex for adults and life stage (larvae, pupae, adults) were investigated for LT determinations.
In both species, females are more tolerant to low and high temperatures than males; larvae and pupae have higher upper lethal limits than do adults. Thermal acclimation of adults has large effects in some instances but small effects in others. Younger adults tend to be more tolerant of low or high temperatures than older age groups. Long-standing laboratory colonies are sufficiently similar in thermal tolerance to field-collected animals to provide reasonable surrogates when making inferences about wild population responses. Differences between these two vectors in their thermal tolerances, especially in larvae and pupae, are plausibly a consequence of different habitat utilization.
Limited plasticity is characteristic of the adults of these vector species relative to others examined to date, suggesting limited scope for within-generation change in thermal tolerance. These findings and the greater tolerance of females to thermal extremes may have significant implications for future malaria transmission, especially in areas of current seasonal transmission and in areas on the boundaries of current vector distribution.
Malaria affects large parts of Africa and Asia and is responsible for nearly 800,000 deaths annually. Despite interventions resulting in a reduction in global malaria mortality in the last 10 years, much concern still exists that in regions where malaria is either endemic, seasonal or has been present in the recent past, climate change might affect its presence and/or prevalence. Forecasts of the effects of climate change on the disease are controversial. Some sources indicate a possible spread of malaria at its current distribution margins, whilst others suggest that climate change will decrease the disease burden in many parts of its current range. In southern Africa, malaria already presents a significant health risk, and how climate change will influence malaria incidence in this region[5–7], depends on several factors which remain poorly understood. These include the form of the change in climate, the environmental responses of the vectors[9, 10], parasite-host interactions[11, 12], and how interventions might interact with these changes[13, 14].
In southern Africa, Plasmodium falciparum, the causative agent of cerebral malaria and the most common of the malarias in Africa, is transmitted by three primary vector species – Anopheles gambiae, Anopheles arabiensis and Anopheles funestus[15, 16]. Current climate change forecasts for the parts of the region where malaria is endemic suggest an increase in both temperature and rainfall, both of which could increase the numbers of mosquitoes and hence the number of cases of the disease. However, many factors remain to be clarified, including how the vectors will respond to such changed climatic conditions. Information on the response of vectors in southern Africa to a variety of conditions is necessary to forecast any change in malaria burden due to changing climates.
Understanding the likely future abundance and distribution of free-living organisms (including malaria vectors) usually involves some form of species distribution modelling, either using environmental niche modelling or a more mechanistic approach. Both approaches have been used to estimate the impacts of climate change on mosquito vectors[18–20], and it has been suggested that a combination of the two can provide the most insight because both the fundamental and realized niches can be estimated (or a sound assessment made of all the factors influencing abundance and distribution). For mechanistic models, typically a range of basic physiological information is required, such as thermal tolerance limits, desiccation resistance and development rate[9, 21].
Because many insect species show phenotypic plasticity, because the sexes often differ in their thermal response, and because tolerances may change with age, and age is an important feature of structured population models[24, 25], these aspects should ideally be investigated too. A further complication is the fact that for many vectors, populations long-established in the laboratory are used for assessments, but laboratory adaptation might affect the outcome of the assays[26–28]. In consequence, these factors must be considered when providing information that can be used for mechanistic niche modelling.
For Anopheles mosquitoes, information on physiological tolerances required for such species distribution modelling is largely lacking. Of the three primary southern African malaria vectors, An. gambiae has been the most widely studied from this perspective, followed by An. arabiensis, but information on the physiological responses of An. funestus is largely absent (a review of the published information is available from the authors on request). Furthermore, the immature forms of An. arabiensis and An. funestus have rarely been considered.
Here, comprehensive assessments of the thermal tolerances of these species, their phenotypic responses to short-term changes in the thermal environment, and an estimate of the extent of laboratory adaptation of these thermal tolerance traits are provided. Additionally, information on the upper and lower lethal temperature limits of the larval and pupal stages of both species is provided. Finally, how climate change might affect vector populations, and hence malaria transmission is briefly considered.
Two long-established laboratory colonies held at the Vector Control Reference Unit in Johannesburg were used for all investigations of thermal tolerance. Anopheles arabiensis was taken from the KGB colony, originally established in 1975 from Kanyemba in the Zambezi Valley, Zimbabwe (R.H. Hunt, pers. comm.) and An. funestus from the FUMOZ colony established in 2000 from southern Mozambique. These colonies are maintained under an insectary temperature of 25 °C (± 1 °C) and 80% relative humidity (verified using repeated measures with a Masons Hygrometer, Brannan, UK) with 12:12 light/dark cycle and 45 min dusk/dawn simulation. Larvae are fed a mixture of ground-up dog biscuits and yeast extracts and females are offered a blood meal three times weekly and allowed to lay eggs two to three times weekly. All adults are provided with a 10% sugar water solution ad libitum.
Anopheles arabiensis females were collected from Malahapanga in the Kruger National Park, South Africa (22° 53.23 S, 31° 02.22 E) in October 2010. Wild An. funestus females were collected from villages surrounding the Maragra Sugar Estate in southern Mozambique (25° 27.41 S, 32° 46.59 E) in April 2011. Adult anophelines were collected using active-search techniques from inside huts or houses or from indoor animal dwellings using a flashlight and 30 cm glass aspirator. Females were transported back to the laboratory within three days for egg-laying in polystyrene cups with rough surfaces at a density of 20 females per 250 ml and were provided with a ball of cotton wool moistened with 10% sugar water solution. Egg batches from these females were kept separate until positive species identifications of the wild adults were made using standard PCR methods[31, 32]. The progeny of at least 80 individual females was used to establish a laboratory colony of the wild strains, with the fifth to seventh generations being used in experiments on An. arabiensis, and the first generation used in experiments on An. funestus. Different generations were used as a result of the inherent difficulties associated with establishing An. funestus colonies compared with An. arabiensis colonies (R.H. Hunt, pers. comm.). These colonies were kept under the same conditions as the laboratory strains.
Critical thermal limits (CTL)
Three age groups for each of the laboratory strains were used. An. arabiensis adults were 10-, 15- and 20-day olds, while An. funestus adults were 10-, 20- and 30-day olds. These ages were chosen because of the different lengths of the gonotrophic cycle and different adult longevities of the two species[30, 33]. Only two adult age comparisons for the wild An. arabiensis strain (10- and 15-day olds) and wild An. funestus strain (10- and 20-day olds) were possible due to low colony numbers and the requirement to make assessments before 10 generations in the laboratory.
Between 20 and 40 individual males and females from all age groups were exposed to each of three acclimation treatments prior to CT determinations. Adult mosquitoes were acclimated for a period of five to seven days at 20 °C, 25 °C or 30 °C and a RH > 80% at either insectary conditions (25 °C) or using PTC-1 Peltier portable temperature control cabinets (Sable Systems, Las Vegas, Nevada, USA, 20 ± 1 °C and 30 ± 1 °C). Humidity in the insectary was checked using a Masons hygrometer (Brannan, UK). At 20 °C and 30 °C, humidity was maintained through the use of distilled water (checked using a Hygrochron i-button, DS 1923-F5, Maxim/Dallas Semiconductor, Sunnyvale, CA, USA). Each acclimation treatment was maintained on a 12L:12D cycle for the five or seven day period. Most insect species show acclimation responses in less than seven days. Following these acclimation treatments, 10 individuals (comprising five individuals of each sex) per individual trial were placed into a double-jacketed insulated chamber connected to a programmable water bath (Grant LTC-12 Series, Grant Instruments, Ltd., Cambridge, UK). For each age group and acclimation treatment of each species, a total of four replicate trials were completed. CTmin experiments started at 20 °C while CTmax started at 25 °C, decreasing or increasing at a rate of 0.25 °C/min, respectively after an equilibration period of 10 min. While it has been shown that rate of temperature change can significantly alter the upper thermal tolerances of various insect species, the current rate was chosen as one comparable with many other studies. The CTmin was regarded as the point where individuals displayed reduced motor function (ie, onset of spasms) and could not cling to the tip of a paint brush, while CTmax was regarded as the point where individuals displayed reduced motor function following a period of rapid flight. At no point were individuals removed from the trial chambers for assessments of motor function (ie, individuals were continuously subjected to the thermal assay).
Lethal temperature limits
Lethal temperature (LT) determinations of larvae and pupae, most appropriate for less mobile stages, for both species were carried out on six groups of 10 individuals each, per life stage (n = 60 per exposure temperature). The plunging technique was used instead of a ramping protocol[37, 38]. Each replicate (ie, group of 10 individuals) was exposed for a period of two hours to temperatures ranging from −12 °C to 8 °C for LLT (lower lethal temperature) and from 34 °C to 44 °C for ULT (upper lethal temperature) in 2 °C increments to ensure that 0% and 100% survival of test individuals was recorded. A water temperature of 24 ± 0.5 °C was used as a control and survival at this temperature was 100%. Temperatures were maintained through the use of programmable water baths (Grant LTD-20 and GR150 R4 Series, Grant Instruments, Ltd., Cambridge, UK). Following the two-hour exposure, experimental groups were returned to water at 24 °C (± 1 °C) and survival was scored every 24 hours until either eclosion to adulthood or complete mortality occurred. Percentage survival was then scored as the percentage of the 10 individuals that eclosed. Larvae were fed daily on the same larval food as the colony strains.
Adult lethal temperature experiments were performed on five groups of 10 individual males and females each (n = 50 individuals per sex, per temperature), acclimated at only one temperature (25 °C, RH 80%). The three age groups (An. arabiensis: 10-, 15- and 20-days; An. funestus 10-, 20- and 30-days), were used in the upper lethal temperature and the LLT experiments, with the exception of the LLT determinations for An. arabiensis adults where only two age groups (10- and 15-day olds) were used due to unexpected mortality in the colony. Each replicate (ie, group of 10 individuals) was exposed to a given temperature in the range −6 °C to 16 °C for LLT determinations and 24 °C to 38 °C for ULT determinations, for a period of four hours to ensure that 0% and 100% survival temperature was measured for both LLT and ULT. This four-hour temperature exposure was chosen as an estimate of the length of time of the hottest period in the day, to which mosquitoes would be exposed, based on generalized daily temperature profile data which show that for many regions, including those of tropical Africa, high daytime temperatures are maintained for approximately four hours[39, 40]. Experiments were conducted in a SANYO incubator (MIR-154, SANYO Electric Co. Ltd., Osaka, Japan). A temperature of 25 ± 1 °C was chosen as a control and survival at this temperature was close to 100%. Adults were immediately removed from the exposure temperature following the four-hour period, given sugar water and left to recover at 25 °C (± 1 °C) and relative humidity of 80%. Survival was scored as the percentage of the 10 adults still living, 24 hours after the experiment concluded.
Normality and homogeneity of variances were examined using Shapiro-Wilk’s tests and Levene’s tests, respectively (Statistica v. 11, StatSoft, Tulsa, Oklahoma, USA). Some deviations from normality were observed, but the model assumptions were generally met (supplementary materials, Additional file1) and the sample sizes sufficiently large to allow for the use of parametric general linear models (GLM), as implemented in R (v. 2.13.1) (R Foundation for Statistical Computing, Vienna, Austria). The first model examined the effects of age, acclimation, sex and strain on the variables CTmin and CTmax for each species. Because significant effects of strain or an interaction with strain were found for CTmin/max, models were then run separately for each strain incorporating age, sex and acclimation as predictor variables. As an estimate of effect size, the mean percent deviation in CTmin/max from the grand mean per group was calculated by subtracting from each factor mean, the grand mean, and dividing this by the grand mean, multiplied by 100 to obtain a percentage (Additional file2). The sign of this % deviation from the mean provides an indication of whether or not each factor had on average a lower (negative) or higher (positive) CTmin/max than the grand mean.
The mean (± S.E.) lethal temperature at which 50% of the sample population died (LT50) for each species and life stage, in relation to age and sex (for adults) was determined through the use of logistic regression with binomial distributions (logit link) in R (v.2.13.1). Using Hochberg’s GT-2 method as described in, lower and upper 95% confidence limits for each group were calculated using the means and standard errors obtained from logistic regression analyses. Mean LT50 (± 95% C.I.) for each group was plotted. Overlapping confidence intervals indicate no significant difference between groups.
Critical thermal limits of Anopheles funestus
No significant differences in CTmax were found between the wild and laboratory strains of An. funestus mosquitoes, and the interactions involving strain were generally not significant, except in a single case (Table 1, Figure 1). Only acclimation affected CTmax values significantly in both strains (Table 1), although the effect size was typically ≤2 °C, with the significant 3-way interaction between strain, age and acclimation not being clearly interpretable (Figure 1). However, it is clear that the overall acclimation response is less in the laboratory than wild strain in both males and females, explaining the significant two-way interaction between strain and acclimation. Clearly, some difference in the effects of acclimation, sex and age exists among strains and therefore the models were run separately for each strain (Table 2). Acclimation and age have much greater effects on CTmax in the wild than in the laboratory strains, with some differences in the interactions too. However, the total variation in CTmax was c. 3 °C (Figure 1, Additional file2). Generally, higher acclimation treatments resulted in higher CTmax, and younger adults and females tend to have higher CTmax values (Figures 1 and2). CTmin values differed between the An. funestus strains, which also showed significant differences in response to acclimation (Table 1). The strongest acclimation response was found in the colony strain and specifically in 10-day old males and females, whereas by comparison differences among other ages and among genders and acclimation treatments at other ages were much reduced (Figure 1). Maximum effect size of c. 4 °C was found following different acclimation treatments in 10-day old colony females (Figure 1, Additional file2). When the models were run independently for the two strains it became clear that in each strain, sex, age and acclimation temperature had significant effects, but in somewhat different ways among strains, with the effects tending to be most pronounced in the laboratory strains (Table 3). Across the full set of treatments, the maximum difference in CTmin was c. 6 °C (Figure 1).
Critical thermal limits of Anopheles arabiensis
In An. arabiensis, with the exception of two-way interactions between sex and age, and sex and acclimation treatment, as well as several three-way interactions between strain, sex, age and acclimation treatment, no other effects on CTmax were significant, and especially not the main effects in the model (Table 4). It does appear that females have higher CTmax values than males (Figure 2), but these effects were not readily distinguished in the full model. When the models were implemented separately for each strain, the sex effect was significant (Table 5, Figure 2), as was the effect of acclimation for the wild strain, largely reflecting the large effect of the 30 °C acclimation treatment on 15-day old males (Figure 3). By contrast, age and various interactions did not have significant effects on the laboratory strain. The overall range of CTmax values was c. 3 °C (Figure 3).
CTmin responded strongly to acclimation treatments, and age, sex and strain were also all significant in An. arabiensis (Table 4). The wild strain tended to have lower CTmin values than the laboratory strain, while 10-day old females in the laboratory colony showed the strongest response to acclimation (Table 6, Figure 3), just as was the case in An. funestus. In the wild strain, females tended to have a lower CTmin than males (Figure 2), and acclimation had a strong, generally linear effect on CTmin (Table 6, Figure 3). However, in the laboratory strain, although all of the main effects and interactions were significant (Table 6), the responses were non-linear among acclimation treatments, and the variation among age groups at a given acclimation ≤ 1.5 °C (Figure 3). Overall, among strains, ages, sexes and acclimation treatments the variation in CTmin was c. 5 °C (Figure 3, Additional file2).
Lethal temperature limits
Lower lethal temperature (LLT) in An. funestus was approximately −1 °C to −2 °C for all stages and age groups examined, with the exception of the larvae (mean ± 95% C.I., 1.94 °C ± 0.62 °C), and 30-day old adult males (mean ± 95% C.I., 0.68 °C ± 0.83 °C), which were less tolerant of low temperature (Figure 4). In An. arabiensis, the situation was similar, with larvae likewise showing the least tolerance of low temperatures (mean ± 95% C.I., 1.59 °C ± 0.71 °C), and adult males being the least resistant of all groups (10-day old males mean ± 95% C.I., 3.66 °C ± 0.98 °C; 15-day old males mean ± 95% C.I., 3.48 °C ± 0.83 °C). Lower lethal limits in the adults were generally 8-11 °C less than the CTmin. The full range of LLT values spanned c. 6 °C (Figure 4).
Upper lethal temperatures (ULT) across the full range of stages, ages and species varied by c. 11 °C. In both species, larvae and pupae had the highest ULT, with An. arabiensis having more tolerant immature stages than An. funestus (Figure 5). Females of both species tended to have higher ULT than males, with the most heat sensitive group being the males of An. arabiensis. The lethal temperature estimates were typically 8-10 °C lower than the CTmax estimates, indicating a much reduced scope for long-term tolerance of high temperature in the adults.
Laboratory colonies are used for a wide range of investigations of insect responses to changing environmental conditions. These include investigations of the responses of mosquitoes to various thermal conditions (eg,), and to pathogens and insecticides. However, as has now been demonstrated in a range of arthropod taxa, laboratory adaptation and acclimation can take place rapidly, affecting some traits, but not others and affecting sexes differentially[26–28]. In consequence, extrapolations to the field situation, such as is required for mechanistic niche modelling or assessments of the outcomes of control interventions, may be compromised, making estimations of the extent of differences among laboratory and field strains essential.
The current results indicate that differences in mean CTmin or CTmax among the wild and laboratory strains of An. arabiensis and An. funestus typically did not exceed 2 °C. In most instances differences between strains were approximately 1 °C. The 2 °C difference among wild and laboratory strains was observed for CTmin in the longest-lived colony (35 years) of An. arabiensis and might indicate a loss of thermal tolerance after extensive exposure to constant laboratory colony conditions. Differences in the acclimation responses between the wild and laboratory strains were also evident. However, the range of acclimation responses over all treatments was similar for both strains except in younger females of the laboratory colonies. Thus, results suggest that although caution is required when extrapolating laboratory thermal tolerance data to the field, as is recommended for other aspects of malaria biology (eg,), at least for the species examined here, using thermal tolerance data from laboratory colonies will provide a reasonable approximation of expected responses in the field (but see also).
Other biologically significant sources of variation in thermal tolerance limits, especially in the context of understanding and forecasting responses to environmental change, are those associated with age, sex and short-term responses to change (phenotypic plasticity)[22, 25, 46]. Several recent studies have shown that upper lethal limits or limits to activity in insects and other ectotherms are typically much less variable, both among populations and species, and over time (through plasticity or responses to selection) than are lower limits[47–49]. The same pattern was found here for the adults and in addition the extent of variation amongst the age groups in CTmax and ULT tended to be fairly narrow. Thus, whilst increasing temperatures may benefit the species in cooler areas (contributing perhaps to rising malaria incidence as is the case in East Africa, eg,), where they are close to their thermal limits rising temperatures may act to suppress populations. Indeed, constrained upper thermal limits may be the mechanistic basis, together with the thermal sensitivity of immature development (see, unpublished data), for the forecast range declines of An. gambiae and An. arabiensis in northern and west Africa and increases in south-eastern Africa. In this respect, males might be more sensitive than females given the 1-2 °C difference in CTmax, which may well be associated with the blood-feeding habits of females.
Transmission of malaria is dependent on the effects of ambient temperature on the Plasmodium parasite, and on the effect of ambient temperature on the vector species. Lower limits to Plasmodium development are c. 16 °C. However, although parasite development rate increases with increasing temperatures, temperatures above c. 30 °C are detrimental to parasite development and could therefore, have consequences for transmission. Transmission of malaria is also dependent on the ability of adults to withstand high temperatures and the greater sensitivity of older mosquitoes to high temperatures, as found here in most cases, may cause female death before the parasite migrates to the salivary glands and can be transmitted. The complexity of this interaction between the sensitivity to temperature of malaria parasites and their vectors has been noted previously. In consequence, rising temperatures may not only reduce mosquito population densities, but also the extent to which the malaria parasite is transmitted. The dynamics of this interaction are likely to be complicated by the habits of the vector (indoor or outdoor species), temperature variability, and the nature of the host-parasite interaction[12, 50]. Nonetheless, the finding that temperature sensitivity, at least for critical limits, increases with age is in keeping with other studies of thermal responses in insects. For mosquitoes, mortality is highly age- and infection-dependent. Even in the absence of Plasmodium ookinetes, mortality of females increases with an increase in age, suggesting the potential for female anophelines to senesce. As a confounding factor to malaria transmission, the presence of large numbers of parasite ookinetes in the mosquito midgut greatly increases the mortality experienced within a population and reduces overall mosquito longevity adding to the potential for reduced overall mosquito populations, and hence, the potential for reduced malaria transmission, with increasing environmental temperatures as a consequence of climate change.
Variation found among the lethal and critical thermal limits for the two anopheline species is typical of that found in a range of other taxa[46, 55]. Activity tends to cease well before the lower lethal limit in adults, whilst the upper lethal limits tend to be somewhat lower than the short-term tolerances represented by CTmax. The latter may in part be explained by the differences between the two techniques used to measure these variables and the rate at which temperature was changed during the ramping method used for CTmax estimation. Slower rates often, though not always, result in lower CTmax values[35, 56]. Nonetheless, these thermal traits might also be under different genetic control. Irrespective, it is clear that the most pronounced differences in ULT were found among the stages, with the immatures having ULTs 2-10 °C higher than those of the adults. Such among-stage variation is common in other insects and usually reflects their exposure to different conditions. For An. arabiensis and An. funestus, as with many other species where the adults are more mobile than the immatures, greater tolerance to high temperatures can be expected in the immature stages. Behavioural regulation is more straightforward for a highly mobile individual living in air than for a much less active individual living in a thermally conductive medium such as water (see also). In particular, the adults of both species are highly anthropophilic, with An. arabiensis displaying more exophilic behaviour than An. funestus[15, 59]. This behaviour of the adults, combined with their mobility, means that they are able to escape unfavourable temperatures and make use of indoor-resting behaviour during the hottest or coldest parts of the day. However, behavioural avoidance of temperature extremes is likely limited for larvae of An. funestus and is probably largely absent for larval An. arabiensis, owing to their breeding habits (see Table 7).
In absolute terms, the larvae of both An. arabiensis (ULT50c. 41 °C) and An. funestus (ULT50c.38 °C) were able to survive higher temperatures (to eclosion) than are those of An. gambiae s.s. (ULT50c.32 °C), although generally the lethal limits were within the range found for anophelines[62–65]. Anopheles arabiensis breeds in shallow, temporary pools or puddles, while An. funestus prefers to breed in semi-permanent to permanent water bodies (Table 7). The smaller water bodies are likely to show much greater thermal variation than the latter simply on the grounds of volume alone, and are also likely to offer less opportunity for microhabitat selection. Thus, the high upper thermal tolerances of An. arabiensis in the immature stages are not unexpected. Nonetheless, how the lethal limits determined here relate to thermal limits to development over the entire immature stage, given that the latter are typically narrower than the former needs to be explored, especially in determining the environmental limits to distribution both at the upper and lower temperature extremes. Interactions between changing climates and lower development limits may account for forecasts of expansion of An. arabiensis into cooler areas as climates warm, given that such interactions can reasonably account for current coarse-scale distributions of An. gambiae s.s.. Furthermore, interactions between climate and upper development limits may change the seasonality of occurrence or lead to range limitation, depending on interactions with rainfall (see eg). Investigations of the relationship between lethal and development limits for both An. arabiensis and An. funestus are currently underway (unpublished data), and should provide insights into changing distribution patterns and the extent to which they match those forecast on the basis of environmental niche modelling alone (see).
This study has shown that with the necessary caution, laboratory colonies provide an initial basis for investigating physiological tolerances of An. arabiensis and An. funestus to both high and low temperatures. In addition, it suggests that limited variation in upper thermal limits may well account for forecasts of declining distributions in already warm areas as temperatures rise, whilst sensitivity of development may be more significant a limiting factor in cool areas, given low, lower lethal limits. Finally, this study has demonstrated substantial physiological differences in tolerance between two of the main malaria vectors in southern Africa, which will have to be taken into account when forecasting responses to environmental change of all kinds, including the ways in which water bodies are manipulated to account for expected changes in rainfall regimes.
World Health Organization: World Malaria Report. 2010, Geneva, Switzerland,http://www.who.int/malaria,
Martens WJM, Jetten TH, Focks DA: Sensitivity of malaria, schistosomiasis and dengue to global warming. Clim Chang. 1997, 35: 145-156.
Small J, Goetz SJ, Hay SI: Climatic suitability for malaria transmission in Africa, 1911–1995. Proc Natl Acad Sci USA. 2003, 100: 15341-15345.
Da Silva J, Garanganga B, Teveredzi V, Marx SM, Mason SJ, Connor SJ: Improving epidemic malaria planning, preparedness and response in Southern Africa. Malar J. 2004, 3: 37-
Rogers DJ, Randolph SE: The global spread of malaria in a future, warmer world. Science. 2000, 289: 1763-1766.
Tanser FC, Sharp B, le Sueur D: Potential effect of climate change on malaria transmission in Africa. Lancet. 2003, 362: 1792-1798.
Van Lieshout M, Kovats RS, Livermore MTJ, Martens P: Climate change and malaria: analysis of the SRES climate and socio-economic scenarios. Global Environ Chang. 2004, 14: 87-99.
New M, Hewitson B, Stephenson DB, Tsiga A, Kruger A, Manhique A, Gomez B, Coelho CAS, Masisi DN, Kululanga E, Mbambalala E, Adesina F, Saleh H, Kanyanga J, Adosi J, Bulane L, Fortunata L, Mdoka ML, Lajoie R: Evidence of trends in daily climate extremes over southern and west Africa. J Geophys Res. 2006, 111: D14102-10.1029/2005JD006289.
Kearney M, Porter WP, Williams C, Ritchie S, Hoffmann AA: Integrating biophysical models and evolutionary theory to predict climatic impacts on species’ ranges: the dengue mosquito Aedes aegypti in Australia. Funct Ecol. 2009, 23: 528-538.
Kearney M, Porter W: Mechanistic niche modelling: combining physiological and spatial data to predict species’ ranges. Ecol Lett. 2009, 12: 334-350.
Paaijmans KP, Read AF, Thomas MB: Understanding the link between malaria risk and climate. Proc Natl Acad Sci USA. 2009, 106: 13844-13849.
Pascual M, Dobson AP, Bouma MJ: Underestimating malaria risk under variable temperatures. Proc Natl Acad Sci USA. 2009, 106: 13645-13646.
Rogers DJ: Models for vectors and vector-borne diseases. Adv Parasit. 2006, 62: 1-35.
Rogers DJ, Randolph SE: Climate change and vector-borne disease. Adv Parasit. 2006, 62: 345-381.
Gillies MT, Coetzee M: A supplement to the Anophelinae of Africa south of the Sahara (Afrotropical Region). 1987, Johannesburg: Publications of the South African Institute of Medical Research
Collins FH, Besansky NJ: Vector biology and the control of malaria in Africa. Science. 1994, 264: 1874-1875.
Boko M, Niang I, Nyong A, Vogel C, Githeko A, Medany M, Osman-Elasha B, Tabo R, Yanda P: Climate Change 2007: Impacts, Adaptation and Vulnerability. Contribution of Working Group II to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change. Edited by: Parry ML, Canziani OF, Palutikof JP, Linden PJ, Hanson CE. 2007, Cambridge: Cambridge University Press, 433-467.
Martin PH, Lefebvre MG: Malaria and climate: sensitivity of malaria potential transmission to climate. Ambio. 1995, 24: 200-207.
Hay SI, Rogers DJ, Randolph SE, Stern DI, Cox J, Shanks GD, Snow RW: Hot topic or hot air? Climate change and malaria resurgence in East African highlands. Trends Parasitol. 2002, 18: 530-534.
Lindsay SW, Bayoh MN: Mapping members of the Anopheles gambiae complex using climate data. Physiol Entomol. 2004, 29: 204-209.
Crozier L: Warmer winters drive butterfly range expansion by increasing survivorship. Ecology. 2004, 85: 231-241.
Chown SL, Terblanche JS: Physiological diversity in insects: ecological and evolutionary contexts. Adv Insect Physiol. 2007, 33: 50-152.
Krebs RA, Loeschcke V: Acclimation and selection for increased resistance to thermal stress in Drosophila buzzatii. Genetics. 1996, 142: 471-479.
Kareiva P, Mullen A, Southwood R: Population dynamics in spatially complex environments: Theory and data [and discussion]. Phil Trans R Soc B. 1990, 330: 175-190.
Bowler K, Terblanche JS: Insect thermal tolerance: what is the role of ontogeny, ageing and senescence?. Biol Rev. 2008, 83: 339-355.
Harshman LG, Hoffmann AA: Laboratory selection experiments using Drosophila: what do they really tell us?. Trends Ecol Evol. 2000, 15: 32-36.
Sgrò CM, Partridge L: Evolutionary responses of the life history of wild-caught Drosophila melanogaster to two standard methods of laboratory culture. Am Nat. 2000, 156: 341-353.
Huho BJ, Ng’habi KR, Killeen GF, Nkwengulila G, Knols BGJ, Ferguson HM: Nature beats nurture: a case study of the physiological fitness of free-living and laboratory-reared male Anopheles gambiae s.l. J Exp Biol. 2007, 210: 2939-2947.
Kirby MJ, Lindsay SW: Effect of temperature and inter-specific competition on the development and survival of Anopheles gambiae sensu stricto and An. arabiensis larvae. Acta Trop. 2009, 109: 118-123.
Hunt RH, Brooke BD, Pillay C, Koekemoer LL, Coetzee M: Laboratory selection for and characteristics of pyrethroid resistance in the malaria vector Anopheles funestus. Med Vet Entomol. 2005, 19: 271-275.
Scott JA, Brogdon WG, Collins FH: Identification of single specimens of the Anopheles gambiae complex by the polymerase chain reaction. Am J Trop Med Hyg. 1993, 49: 520-529.
Koekemoer LL, Kamau L, Hunt RH, Coetzee M: A cocktail polymerase chain reaction assay to identify members of the Anopheles funestus (Diptera: Culicidae) group. Am J Trop Med Hyg. 2002, 6: 804-811.
Munhenga G, Brooke BD, Chirwa TF, Hunt RH, Coetzee M, Govender D, Koekemoer LL: Evaluating the potential of the sterile insect technique for malaria control: relative fitness and mating compatibility between laboratory colonized and wild population of Anopheles arabiensis from the Kruger National Park. South Africa. Parasit Vectors. 2011, 4: 208-
Weldon CW, Terblanche JS, Chown SL: Time-course for attainment and reversal of acclimation to constant temperature in two Ceratitis species. J Thermal Biol. 2011, 36: 479-485.
Chown SL, Jumbam KR, Sørenson JG, Terblanche JS: Phenotypic variance, plasticity and heritability estimates of critical thermal limits depend on methodological context. Funct Ecol. 2009, 23: 133-140.
Lutterschmidt WI, Hutchison VH: The critical thermal maximum: data to support the onset of spasms as the definitive endpoint. Can J Zool. 1997, 75: 1553-1560.
Chown SL, Nicolson SW: Insect Physiological Ecology. Mechanisms and Patterns. 2004, Oxford: Oxford University Press
Mitchell KA, Hoffmann AA: Thermal ramping rate influences evolutionary potential and species differences for upper thermal limits in Drosophila. Funct Ecol. 2010, 24: 694-700.
Kingsolver JG: Thermal and hydric aspects of environmental heterogeneity in the pitcher plant mosquito. Ecol Monogr. 1979, 49: 357-376.
Bonan GB: Ecological Climatology: Concepts and Applications. 2002, Cambridge: Cambridge University Press
Quinn GP, Keough MJ: Experimental Design and Data Analysis for Biologists. 2002, Cambridge: Cambridge University Press
Sokal RR, Rohlf FJ: Biometry: The Principles and Practice of Statistics in Biological Research. 1995, New York: W.H. Freeman and Company
Bayoh MN, Lindsay SW: Effect of temperature on the development of the aquatic stages of Anopheles gambiae sensu stricto (Diptera: Culicidae). Bull Ent Res. 2003, 93: 375-381.
Blanford S, Read AF, Thomas MB: Thermal behaviour of Anopheles stephensi in response to infection with malaria and fungal entomopathogens. Malar J. 2009, 8: 72-
Aguilar R, Dong Y, Warr E, Dimopoulos G: Anopheles infection responses; laboratory models versus field malaria transmission systems. Acta Trop. 2005, 95: 285-291.
Hoffmann AA, Sørensen JG, Loeschcke V: Adaptation of Drosophila to temperature extremes: bringing together quantitative and molecular approaches. J Thermal Biol. 2003, 28: 175-216.
Chown SL: Physiological variation in insects: hierarchical levels and implications. J Insect Physiol. 2001, 47: 649-660.
Griffiths JA, Schiffer M, Hoffmann AA: Clinal variation and laboratory adaptation in the rainforest species Drosophila birchii for stress resistance, wing size, wing shape and development time. J Evol Biol. 2005, 18: 213-222.
Sunday JM, Bates AE, Dulvy NK: Global analysis of thermal tolerance and latitude in ectotherms. Proc R Soc B. 2011, 278: 1823-1830.
Alonso D, Bouma MJ, Pascual M: Epidemic malaria and warmer temperatures in recent decades in an East African highland. Proc R Soc B. 2011, 278: 1661-1669.
Townsend Peterson A: Shifting suitability for malaria vectors across Africa with warming climates. BMC Infect Dis. 2009, 9: 1-6.
Benoit JB, Lopez-Martinez G, Patrick KR, Phillips ZP, Krause TB, Denlinger DL: Drinking a hot blood meal elicits a protective heat shock response in mosquitoes. Proc Natl Acad Sci USA. 2011, 108: 8026-8029.
Craig MH, Snow RW, le Sueur D: A climate-based distribution model of malaria transmission in sub-Saharan Africa. Parasitol Today. 1999, 15: 105-111.
Dawes EJ, Churcher TS, Zhuang S, Sinden RE, Basáñez M-G: Anopheles mortality is both age- and Plasmodium-density dependent: implications for malaria transmission. Malar J. 2009, 8: 228-
Terblanche JS, Hoffmann AA, Mitchell KA, Rako L, le Roux PC, Chown SL: Ecologically relevant measures of tolerance to potentially lethal temperatures. J Exp Biol. 2011, 214: 3713-3725.
Allen JL, Clusella-Trullas S, Chown SL: The effects of acclimation and rates of temperature change on critical thermal limits in Tenebrio molitor (Tenebrionidae) and Cyrtobagous salviniae (Curculionidae). J Insect Physiol. 2012, 58: 669-678.
Hoffmann AA: Physiological climatic limits in Drosophila: patterns and implications. J Exp Biol. 2010, 213: 870-880.
Huey RB: Physiological consequences of habitat selection. Am Nat. 1991, 137: S91-S115.
Coetzee M, Craig M, le Sueur D: Distribution of African malaria mosquitoes belonging to the Anopheles gambiae complex. Parasitol Today. 2000, 16: 74-77.
Paaijmans KP, Thomas MB: The influence of mosquito resting behaviour and associated microclimate for malaria risk. Malar J. 2011, 10: 183-
Bayoh MN, Lindsay SW: Temperature-related duration of aquatic stages of the Afrotropical malaria vector Anopheles gambiae in the laboratory. Med Vet Entomol. 2004, 18: 174-179.
Muirhead Thomson RC: The reactions of mosquitoes to temperature and humidity. Bull Ent Res. 1938, 29: 125-140.
Love GJ, Whelchel JG: Lethal effects of high temperatures on the immature stages of Anopheles quadrimaculatus. Ecology. 1957, 38: 570-576.
Benedict MQ, Cockburn AF, Seawright JA: Heat-shock mortality and induced thermotolerance in larvae of the mosquito Anopheles albimanus. J Am Mosquito Contr. 1991, 7: 547-550.
Raghavendra K, Barik TK, Adak T: Development of larval thermotolerance and its impact on adult susceptibility to malathion insecticide and Plasmodium vivax infection in Anopheles stephensi. Parasitol Res. 2010, 107: 1291-1297.
Paaijmans KP, Wandago MO, Githeko AK, Takken W: Unexpected high losses of Anopheles gambiae larvae due to rainfall. PLoS One. 2007, 2: e1146-10.1371/journal.pone.0001146.
This study was supported by the National Research Foundation through the DST-NRF Centre of Excellence for Invasion Biology, the DST-NRF Research Chair award to MC, and the Stellenbosch University Hope Project. Students and staff of the Vector Control Reference Unit (VCRU) in Johannesburg are thanked for assistance with colony maintenance and establishment, field work and laboratory processing. Field assistance in Mozambique was provided by E Nhamahango from the Maragra Sugar Estate. Colony maintenance was undertaken by the Vector Control Reference Unit and all permits for animal maintenance issued to them were valid for these trials (NHLS Animal Ethics Clearance Certificate #1993047). Two anonymous reviewers are thanked for their helpful comments.
The authors declare that they have no competing interests.
CLL, MC, JST and SLC designed the research. CLL collected the data. CLL and SLC analysed the data. CLL, MC, JST and SLC wrote the manuscript. All authors read and approved the final manuscript.