- Open Access
Modelling adult Aedes aegypti and Aedes albopictus survival at different temperatures in laboratory and field settings
Parasites & Vectorsvolume 6, Article number: 351 (2013)
The survival of adult female Aedes mosquitoes is a critical component of their ability to transmit pathogens such as dengue viruses. One of the principal determinants of Aedes survival is temperature, which has been associated with seasonal changes in Aedes populations and limits their geographical distribution. The effects of temperature and other sources of mortality have been studied in the field, often via mark-release-recapture experiments, and under controlled conditions in the laboratory. Survival results differ and reconciling predictions between the two settings has been hindered by variable measurements from different experimental protocols, lack of precision in measuring survival of free-ranging mosquitoes, and uncertainty about the role of age-dependent mortality in the field.
Here we apply generalised additive models to data from 351 published adult Ae. aegypti and Ae. albopictus survival experiments in the laboratory to create survival models for each species across their range of viable temperatures. These models are then adjusted to estimate survival at different temperatures in the field using data from 59 Ae. aegypti and Ae. albopictus field survivorship experiments. The uncertainty at each stage of the modelling process is propagated through to provide confidence intervals around our predictions.
Our results indicate that adult Ae. albopictus has higher survival than Ae. aegypti in the laboratory and field, however, Ae. aegypti can tolerate a wider range of temperatures. A full breakdown of survival by age and temperature is given for both species. The differences between laboratory and field models also give insight into the relative contributions to mortality from temperature, other environmental factors, and senescence and over what ranges these factors can be important.
Our results support the importance of producing site-specific mosquito survival estimates. By including fluctuating temperature regimes, our models provide insight into seasonal patterns of Ae. aegypti and Ae. albopictus population dynamics that may be relevant to seasonal changes in dengue virus transmission. Our models can be integrated with Aedes and dengue modelling efforts to guide and evaluate vector control, better map the distribution of disease and produce early warning systems for dengue epidemics.
Survival of arthropod vectors is one of the most important components of transmission of a vector-borne pathogen [1–3]. Increased survival allows the vector to produce more offspring, to increase the chances of them becoming infected, to disperse over greater distances, to survive long enough to become infectious, and then to deliver more infective bites during the remainder of its lifetime. As a result, small changes in survival rate cause large changes in the rate of pathogen transmission [3–7]. Furthermore, seasonality and the geographic distribution of vector borne diseases are often constrained by differences in vector survival [8, 9].
Adult female Aedes aegypti (L.) are the principal vectors for dengue and urban yellow fever viruses, two globally important human arboviruses [10–12]. Ae. albopictus (Skuse) can also act as a secondary vector for dengue viruses (DENV) [13, 14], along with 22 other arboviruses [13, 15, 16], some of which have an increasing public health burden . There has, therefore, been considerable interest in investigating and quantifying factors affecting survival of adult females of both species and how this contributes to virus transmission [18, 19]. Because Aedes mosquitoes are small-bodied poikilotherms it is logical that temperature is consistently observed as a principal factor affecting survival [20, 21] and the current global range of both species align broadly with separate critical limits imposed by winter isotherms [9, 22]. While other factors such as humidity and photoperiod are important, the effects of temperature have been most rigorously quantified and most frequently identified as limiting factors for survival [10, 18].
Reconciling the expected difference between laboratory and field survival estimates has been complicated by the lack of precision of available field techniques and the limited temperature ranges over which field experiments have been undertaken. In mosquito cages in the laboratory, conditions can be controlled and the effects of temperature quantified accurately, but additional causes of mortality experienced in the field, such as predation or disease are absent. The most commonly used method for observing Aedes survival in the field is mark-release-recapture (MRR) where marked individuals of a known age are released and the number of recaptures observed over subsequent days . While MRR experiments expose mosquitoes to all causes of mortality in the field, experiments with Ae. aegypti and Ae. albopictus typically have small sample sizes  and sensitivities of spatial mosquito sampling limit the level of survival detail that can be elucidated through the majority of MRR experiments . Furthermore, MRR experiments are limited in number and cover only a narrow range of environmental temperatures (see Additional file 1). There is, therefore, considerable uncertainty over variation in Aedes survival in the field.
Using laboratory estimates of survival to predict field survival is complicated by the different mortality risks in the two settings, unknown importance of senescence in the field, and the transition from constant laboratory to fluctuating natural temperature regimes. It has long been assumed that the high rate of mortality in the field, due to external causes, ensures few mosquitoes live long enough to experience age-dependent mortality . Observations from laboratory and field experiments on Ae. aegypti, however, have begun to challenge this assumption [23, 27]. As a result, the conditions under which age-dependence or age-independence dominates mortality remain unresolved, yet are potentially of considerable epidemiological importance because older mosquitoes are more likely to have survived beyond the virus’s extrinsic incubation period. Furthermore, controlled laboratory conditions are not expected to have the same effect as what is observed in the field. This is particularly apparent with observations at constant temperature in the laboratory in contrast to the fluctuating temperature experienced in the wild. At lower mean temperatures, fluctuation may increase survival because at least part of the day could permit significant increases in survival. The converse is true for optimal mean temperatures, where any positive or negative departure from the optimum could lead to decreased survival [28, 29]. Both of these issues have widened the gap between field and laboratory observed survival in terms of the average and distribution of survival times.
Observations from temperature-controlled laboratory experiments have produced a range of candidate parametric functions that are suitable for modelling age-dependent mortality, such as the Gompertz and Logistic functions [27, 30]. Early experiments were extended by Degallier et al. who observed age-dependent mortality in caged mosquitoes exposed to uncontrolled field temperature regimes . Degallier et al. also extended the parametric analysis by using non-parametric Cox proportional hazard models, which allowed them to analyse the significance of both age and environmental conditions on adult mortality. The importance of seasonal changes in temperature was emphasised by Strickman et al. who conducted year-round field experiments on caged mosquitoes in Thailand . Yang et al. attempted to explicitly model the relationship between temperature and survival, however, their model assumed age-independent mortality and relies on relatively limited data . The paucity of data and the challenges of fitting robust statistical relationships to Ae. albopictus adult longevity data are discussed by Waldock et al. who highlight the disparity relative to the data available for development of the mosquito’s immature stages at a range of different temperatures .
In this study we compiled published observations from a variety of studies and modelled the effect of temperature on survival of adult female Ae. aegypti and Ae. albopictus under laboratory and field conditions. Our models incorporated the effects of age-dependent mortality and predict survival under fluctuating temperature regimes. They provide a more detailed understanding of Aedes adult female mortality and present a statistical solution to investigating age-dependent survival. The resulting models may enhance current DENV transmission models [5, 33–36], which could improve predictions that guide vector control, identify areas suitable for DENV transmission and contribute to outbreak early warning systems.
Relevant publications were collected by searching the databases of PubMed, Google Scholar and the Armed Forces Pest Management Board Literature Retrieval System using the search terms Aedes, aegypti, albopictus, survival, mortality, longevity, stegomyia and albopicta. Resulting abstracts were examined for their likelihood to contain laboratory-based adult female survival experiments (closed cage with controlled or semi-controlled environmental and feeding regimes). A more detailed inspection of the manuscript allowed us to apply further inclusion criteria for the final database of laboratory survival experiments. Included experiments monitored mosquito survival of at least five individuals at least every two days and gave all mosquitoes access to nutrition (sugar or any form of blood feeding). In the field it appears that female Ae. aegypti feed almost exclusively on human blood [37, 38], however, laboratory experiments use a range of diets, such as blood plus sugar, sugar only or blood only, artificial blood meals (blood drawn from a variety of different vertebrates presented to mosquitoes in different ways), and direct feeding on a variety of different vertebrates that can confer significant differences in survival [38–43]. As a result this information was collected and incorporated as a random effect in our analysis and all predictions were made for mosquitoes fed human blood. Due to factors other than mosquito diet, laboratory survival observations still showed a significant amount of variation (see Additional file 1), as a result we decided to include a random effect at the experiment level to control for the wide variety of experimental set ups. All experiments analysed used mosquitoes colonised in the laboratory for less than 2 years. Experiments where mosquitoes were chemically or genetically treated were excluded. Finally, Ae. aegypti used in the survival experiments were geographically distributed across five continents incorporating both tropical and temperate latitudes. There was a slight bias towards North American samples in the Ae. albopictus survival data, however, strains from tropical areas, such as Réunion and Malaysia, were also included.
Literature searches for MRR field experiments were completed independently under a parallel project to investigate mosquito movement (Guerra et al. unpublished observations). Briefly, existing literature was searched using a similar key word strategy, adult female mosquito MRR and ancillary geographic data were extracted from suitable articles. Using these references and ancillary data, daily recapture numbers were extracted from relevant articles. Of the Ae. aegypti MRR studies, 14% released field caught mosquitoes, 36% released laboratory reared mosquitoes from a recently caught strain and 50% released a colonised mosquito strain. For Ae. albopictus, the corresponding breakdown was 4%, 20% and 76%.
Where available, the initial number of mosquitoes (for laboratory experiments) and number of deaths/recaptures on each subsequent day of observation (for MRR experiments) were extracted. Where only graphical summaries were available, these observations were recreated using GetData Graph Digitizer . For any article where observations could not be recorded or recreated the authors were contacted and primary data was requested. For laboratory data the mean temperature was recorded and for field data the maximum and minimum temperatures were recorded. Where unavailable, study site temperatures were estimated using the average climate for the particular location and time of year using WorldClim global climate data . The resulting number of articles and data points included are summarised in Table 1. The full reference list and summary boxplots are available in Additional file 1.
Comparability of data and fitting parametric models to laboratory survival data
A range of parametric equations have been used to model the survival function of laboratory and field adult mosquitoes. For field data, the exponential function, adjusted for recapture number, was the most frequently used model [23, 32, 46–52]. For laboratory data, a range of functional forms were applied with the aim of detecting age-dependent mortality under conditions where overall mortality is lower. These have included the Gompertz function which can allow a non-symmetric mortality rate to increase or decrease exponentially with age [31, 53, 54], the Weibull function where the mortality rate is a power function of age [31, 48, 53, 54] and the log-logistic function where the mortality rate is not restricted to a monotonic function [31, 54, 55]. Here we aimed to compare the fit of these different parametric models to laboratory survival data at each temperature regime across all of the studies to see if any one function offers a significantly improved fit and whether this changes with temperature. Using the laboratory data we re-examined the fit of the log-logistic, Gompertz, exponential and Weibull survival functions. Function fitting was performed in R version 2.14.2  using the “flexsurv” library which fits a variety of parametric models using maximum likelihood. The Akaike Information Criterion (AIC) was calculated for each model and used to calculate the relative likelihood (RL)  of each model to the model with the lowest AIC:
The relative likelihood values for each model were averaged for experiments of equivalent temperature to give a balanced estimate of the each model’s suitability at different temperatures.
Fitting non-parametric models to laboratory survival data
As the distribution of survival times may undergo complex changes across a range of temperatures we also chose to construct a non-parametric model that would not share the same restrictions as the aforementioned parametric models. For this we chose regression spline generalised additive models (GAMs) which apply a smoothing variable to the explanatory variables in order to model the response variable . This method has the advantage of being able to model unknown and non-linear effects of covariates and thus elucidate the potentially complex effect of temperature on adult mosquito survival.
To evaluate the improvement of using GAMs over the parametric alternatives, we fitted parametric and GAM models to each laboratory experiment and calculated the difference in AIC between parametric and non-parametric models across all experiments.
A second GAM was then formulated to use the data from all experiments in one model to recreate the relationship between survival, time and temperature. The GAM was formulated as follows:
S ij = number of mosquitoes surviving at observation i in experiment j
N ij = number of mosquitoes at start of time step at i, j
P ij = survival probability for a mosquito at i, j
f() = smooth term
D i = day of observation i
T i = temperature of observation i
ϵ j = random error term for experiment j
ϵ d = random error term for mosquito diet d
θ j 2 = variance across experiments
θ d 2 = variance across mosquito diets
Smoothing parameters were selected by restricted maximum likelihood with a data-driven basis dimension choice of k D = 8 and 5 and k T = 5 and 5 for Ae. aegypti and Ae. albopictus respectively . Confidence intervals for the interquartile range of predictions were obtained by bootstrapping with 200 repeats, each the size of the original dataset. This model was fit using D ≥ 1 to be consistent with the experimental observations that record mortality. Extrapolated model predictions for 0 ≤ D < 1 were scaled proportionally to ensure 100% survival at D = 0. All GAMs were implemented using the “mgcv” package in R .
For the model to fit biologically appropriate responses, additional data defining the limits of prediction were required. Observations from Christophers , suggest 4°C and 42-43°C as suitable minimum and maximum critical temperatures at which survival of Ae. aegypti is minimal (<24 hours). Similar observations for Ae. albopictus suggest values between -5°C and 40–40.6°C [61, 62]. To constrain mortality in the model, all non right censored experimental observations were extended to 120 days at 0% survival and a maximum lifetime of 120 days was imposed at all temperatures, the maximum longevity observed in our dataset. Furthermore, to produce meaningful estimates of longevity, survival of less than 0.1% of the initial mosquito population was considered sufficient to indicate complete mortality.
Reconciling laboratory survival predictions with survival measurements in the field
Because MRR data was too limited to fit models over a range of temperatures as we did for the laboratory mosquitoes, we used the laboratory-based models above to estimate temperature and age-matched expected mortality rates for mosquitoes in the field and compared those estimates to observed field mortality estimates from MRR data. To simulate the realistic effects of temperature in the field, in contrast to the constant temperatures of the laboratory, we recreated 59 separate daily fluctuating temperature regimes using the maximum and minimum temperatures of each MRR experiment and assumed sinusoidal progression in the day with a decreasing exponential curve at night [28, 63]. These matching fluctuating temperatures were then applied to the laboratory survival model to estimate hourly survival for each of the MRR experiments. For each MRR experiment, we directly estimated the hourly mortality rate using the Buonaccorsi nonlinear regression method applied using the “stats” package in R . This assumes that the mortality rate in the field is constant for each experiment and, therefore, our estimation of external mortality, due to causes such as predation or disease, also has this assumption. We then matched the laboratory-based predicted mortality with the directly estimated mortality and calculated the average difference, which we attribute to mortality due to external factors. We added this estimated hourly external mortality to the temperature-dependent laboratory model to produce the field survival model, which estimates mortality at all temperatures. Confidence intervals for the field survival model were derived by propagating the uncertainty of the laboratory model and combining this with 200 bootstrap samples of the MRR experiments to test the sensitivity of the calculated external mortality value. The overall schematic of the methods appears in Figure 1.
Parametric models of laboratory data
Of the four parametric models tested, no one model was consistently most suitable across the range of temperatures tested (Figure 2). The log-logistic model was most suitable for Ae. aegypti with a relative likelihood of 0.441, whereas the exponential model was most suitable for Ae. albopictus with a relative likelihood of 0.397. The Weibull model was the least suitable for both species with relative likelihoods of 0.07 and 0.05 for Ae. aegypti and Ae. albopictus, respectively. While some models may be more suitable at specific temperatures, such as the exponential at 15°C for Ae. aegypti, they are not consistently more suitable at other temperatures, nor is there a clear trend suggesting some models are more suitable at higher or lower temperatures (Figure 2).
Non-parametric models of laboratory data
When non-parametric GAMs were used to fit the same data, both the overall model likelihood and the number of experiments for which it was optimal increased (Table 2). For both Ae. aegypti and Ae. albopictus, respectively, median difference in AIC between parametric and non-parametric models was 307.6 and 4.38. While still significant, this improvement was less pronounced in Ae. albopictus datasets due to a lower average number of mosquitoes under observation (Table 1), which limited the ability of some experiments to detect finer scale changes in survival. The number of experiments for which GAMs were the most suitable, however, (49.5%, Table 2) showed its improved fit over the entire dataset. The fits of both parametric and non-parametric models to two of the Ae. aegypti experiments are shown in Figure 3 as an example of the fit of each of the five models. Due to its more accurate quantification of mortality rate and its improved fit across a range of temperatures we chose a non-parametric GAM to construct our laboratory temperature survival model.
The non-parametric GAM model allowed us to quantify the effects of age and temperature on mortality whilst still taking into account the random effects at the experiment and mosquito diet level. Figures 4A and 4B show the model fit for each species using the laboratory data. Overall, Ae. albopictus had a greater longevity than Ae. aegypti, however, Ae. aegypti tolerated a wider range of temperatures (Figures 4A and 4B). The survival patterns also differ at the limits of survival, Ae. aegypti had optimal survival over only a narrow window around 21°C, while Ae. albopictus had optimal survival over a much wider range with only minor differences observed between 20-30°C. The mean and 95% survival limits for both species show that mortality was widely distributed across different ages with the majority of mosquitoes dying well before reaching maximum longevity. For Ae. aegypti the predicted median longevity was only 38 days even at optimal temperatures despite some individuals living over 100 days. This same pattern was observable in Ae. albopictus data, but to a lesser extent: median longevity at optimal temperatures was over 60 days despite some individuals living up to 120 days. Random effects for mosquito diet and study level were highly significant for both Ae. aegypti and Ae. albopictus models (Chi square, p < 0.001).
The greatest uncertainty around these predictions occurs around the 50% survival time (Figure 5A and 5B). For the majority of survival curves this is the point at which survival has the greatest rate of change; i.e. highest mortality rate, and it is, therefore, not surprising that uncertainty is higher around these points. Peaks in uncertainty can also be identified at low temperatures for Ae. albopictus (Figure 5B). Overall the uncertainty was higher for Ae. albopictus reflecting the limited data, particularly between 0 and 15°C.
Reconciliation with field data
To estimate longevity in the field, the daily mortality in the laboratory model was compared to field observations of mortality from MRR experiments where other, external factors may contribute to mortality. For Ae. aegypti, the laboratory model-predicted daily mortality rate varied between 0.05 and 0.082 (average 0.064) compared to an observed MRR mortality rate of between 0.033 and 0.595 (average 0.288). The average difference between these predictions for each experiment, and thus the calculated additional daily mortality, was 0.179 (standard deviation 0.012 – 0.346). For Ae. albopictus the external daily mortality was lower. Laboratory model predictions of daily mortality rate varied from 0.0024 to 0.289 (average 0.040) compared to observed MRR mortality rate of between 0.031 and 0.231 (average 0.121) giving an average difference of 0.0621 (standard deviation −0.030 – 0.154). Due to the techniques used to measure external mortality, we had to assume it acted independently of mosquito age or temperature and so predicted survival in the field was calculated as the sum of laboratory model mortality and external mortality and is shown in Figure 4C and 4D.
Clear differences in survival can be observed between laboratory models (Figure 4A and 4B) and field models (Figure 4C and 4D). Survival is decreased for both species, but less so for Ae. albopictus due to the lower estimated external mortality. The effect of incorporating external mortality also reduces the relative influence of temperature and as a result less variation in survival is observed at mid-range temperatures (approximately 15-35°C). The most significant difference between laboratory and field models is the change in median longevity, which is reduced considerably in both species.
When field data uncertainty was incorporated with the existing laboratory model uncertainty, overall uncertainty was reduced due to the higher value of external mortality and the lower uncertainty in estimating this parameter (Figure 5). Uncertainty was higher at extreme temperatures, particularly for Ae. albopictus at lower temperatures, however this uncertainty was confined to a limited range of temperatures and days, ensuring that overall prediction had low uncertainty. The predicted hourly survivorship from the models for Ae. aegypti and Ae. albopictus field survival and their associated uncertainty profiles are available for free download on the online data repository figshare (http://www.figshare.com) [65, 66].
In this study we used data from 410 Ae. aegypti and Ae. albopictus adult female survival experiments to construct a model of temperature and age-dependent survival. The best-fitting model was an age and temperature dependent GAM, which we combined with field-based mortality to estimate the distribution of longevity of mosquitoes in the field. These models should be interpreted in light of the limitations to the data used and the modelling methods. Of the data compiled, there were relatively few experiments conducted at more extreme temperatures, resulting in higher uncertainty in these temperature ranges. Considering the implications these extreme temperatures have for the geographic range of the mosquito, more data is necessary for better defining that range. Mosquito survival may also vary geographically due to variation in the genetic background of local mosquito populations. Our analysis incorporates experimental data from a wide variety of locations and conditions but, given the paucity of MRR data at climatic extremes, it is possible that estimates of field-based mortality may not be representative, particularly in locations where mosquitoes may be more tolerant of extreme climates. However, observations suggest that adult survival is one of the few bionomics of Ae. aegypti that remains consistent among geographically disparate populations [67, 68]. North American samples of Ae. albopictus were over-represented in our laboratory data (91%) which may overestimate thermal tolerance at the lower temperature limits due to evidence of Ae. albopictus adaptation in temperate climates [53, 69]. Perhaps the most important data limitation is that there are very few field mortality experiments upon which to estimate mortality outside of the laboratory, where arbovirus transmission actually occurs. More data would allow better estimation of field mortality and allow relaxation of the assumption that external contributors to mortality are not age-related. An extension should include alternatives to MRR approaches that might be more sensitive for estimating survival of older potentially infectious mosquitoes, for example, using transcript levels of an age-associated gene  or measuring mosquito infection rates .
It is also worth considering the scope of the model we developed. While it may suggest that adult survival is possible even at low temperatures, factors limiting the development of the immature stages may preclude the establishment of an adult population. For adult Aedes, there may also be additional temperature-based limits to survival. Below 14-15°C, Ae. aegypti experiences reduced mobility and struggles to imbibe blood [18, 21, 72]. As observations (at tropical temperatures) suggest Ae. aegypti cannot survive longer than 2-3 days without a blood meal , extended periods of time below 14-15°C are likely to result in complete mortality. We extended our predictions below these temperatures because limited time spent at these colder temperatures, such as a few hours in the early morning, may not confer a lasting negative effect on survival, however, these effective mortality limits must be considered when integrating this into a wider Aedes/DENV transmission model.
Here we have shown that GAMs captured more of the variation in survival between different experiments than conventional, parametric models of mosquito mortality, including the Gompertz and logistic models which previously provided the best fit of the options examined [27, 30, 31]. In our analysis we did not consider frailty-based mortality models, however non-parametric approaches, such as GAM, offer more flexibility than parametric models, suitable predictive power and straight-forward implementation in a variety of statistical software packages. We have also emphasised the importance of including random effects in analyses of this type as it has allowed us to comment on the effect of specific experimental treatments, such as feeding regime, and thus guide laboratory and field experimental design. Furthermore adding a random effect at the study level allowed us to comment on the generalisability of any one experiment to a wider context.
The relative fit of the GAM indicated the presence of age-dependent survival. This age-related effect may be related to senescence or other factors, but appears to be less important in the wild. On average we found that mortality due to external causes was significantly greater (4.5 times for Ae. aegypti and 3.0 times for Ae. albopictus) than mortality due to temperature and age at mid-ranging temperatures (20-30°C) indicating that few mosquitoes will be alive at older ages by the time senescence measurably impacts mortality. Results from our analysis indicate that across different temperature regimes the assumption of age-independent mortality for Ae. aegypti and Ae. albopictus in the field is likely to be appropriate. It may be possible that senescence acts at different ages in the laboratory than it does in the field, however, quantifying this would require unfeasibly large mosquito cohorts, or innovative new techniques to measure mosquito survival in the field. As a result, the statistical model presented here provides a novel approach for quantifying the effects of age-dependent survival among wild mosquitoes.
The greater tolerance of lower temperatures observed for Ae. aegypti compared to Ae. albopictus appears at odds with their observed geographic distribution. This finding can be explained, however, if we consider the adaptations of the egg stages of each mosquito. Aedes albopictus eggs are able to undergo diapause allowing the species to persist during cold winter temperatures that are unfavourable to adult survival [74, 75]. Aedes aegypti shows only limited adaptation to egg stage survival in unfavourable periods [76, 77] and there may, therefore, be greater selection pressure for thermal tolerance at the adult stage to resist diurnal and inter-seasonal variations in temperature.
The seasonal variation and geographic limits of Aedes transmitted viruses are intrinsically linked to the seasonality and geographic limits of mosquito populations required for their transmission. It is, therefore, important that the survival of adult Aedes is well estimated using appropriate data and considering the variety of conditions they may encounter in nature. Quantifying how changes in survival translate to dengue transmission potential is also important and requires integration with models of the DENV extrinsic incubation period, which is also sensitive to temperature . Characterising the interaction between mortality and extrinsic incubation could be used further to explore entomological components of seasonal forcing among Aedes-borne viruses, which may help inform strategies for disease prevention and control.
Given that dengue is a disease that affects over half the world’s population  with an estimated 390 million new infections per year , improvements in existing approaches and new strategies to control the disease are needed. Depending on the situation, both the type of vector control used and how it is implemented must be optimised to achieve the biggest reduction in dengue burden. Modelling and mapping are two important tools for addressing both of these issues, identifying risks and opportunities for control [34, 79, 80]. Considerable efforts have already been undertaken to re-examine many of the fundamental assumptions of early Aedes and dengue transmission models [6, 78, 81, 82]. The next generation of models is likely to incorporate significant advances, including new frameworks for modelling survival of adult mosquitoes. The adult Aedes survival models developed here can be integrated into these new modelling approaches with the aim of refining the global distribution of Aedes vectors, guiding seasonal vector control efforts, modelling dengue transmission and developing early warning systems to prevent dengue epidemics .
We created an explicit model of Ae. aegypti and Ae. albopictus adult survival across a range of temperatures under both laboratory and field conditions. Commonly used parametric models did not capture the variation in survival to the same extent as our approach across a variety of temperatures and experiments. The use of non-parametric GAMs allowed the effects of senescence and temperature to be captured explicitly, revealing the importance of temperature and age-dependent mortality under laboratory conditions and fluctuating temperatures and external mortality in the field. The models developed here can be used for investigating seasonal variation in dengue transmission and will strengthen components to dengue transmission models that investigate various options for control.
Akaike’s information criterion
Generalised additive model.
Garrett-Jones C, Shidrawi G: Malaria vectorial capacity of a population of Anopheles gambiae: an exercise in epidemiological entomology. Bull World Health Organ. 1969, 40 (4): 531-
Reisen WK, Mahmood F, Parveen T: Anopheles culicifacies Giles: a release-recapture experiment with cohorts of known age with implications for malaria epidemiology and genetical control in Pakistan. Trans R Soc Trop Med Hyg. 1980, 74 (3): 307-317. 10.1016/0035-9203(80)90089-9.
Macdonald G: Epidemiological basis of malaria control. Bull World Health Organ. 1956, 15 (3–5): 613-
Parham PE, Michael E: Modeling the effects of weather and climate change on malaria transmission. Environ Health Perspect. 2010, 118 (5): 620-
Barbazan P, Guiserix M, Boonyuan W, Tuntaprasart W, Pontier D, Gonzalez JP: Modelling the effect of temperature on transmission of dengue. Med Vet Entomol. 2010, 24 (1): 66-73. 10.1111/j.1365-2915.2009.00848.x.
Smith DL, Battle KE, Hay SI, Barker CM, Scott TW, McKenzie FE: Ross, Macdonald, and a theory for the dynamics and control of mosquito-transmitted pathogens. PLoS Pathog. 2012, 8 (4): e1002588-10.1371/journal.ppat.1002588.
Lunde TM, Bayoh MN, Lindtjørn B: How malaria models relate temperature to malaria transmission. Parasite Vector. 2013, 6: 20-10.1186/1756-3305-6-20.
Gage KL, Burkot TR, Eisen RJ, Hayes EB: Climate and vectorborne diseases. Am J Prev Med. 2008, 35 (5): 436-450. 10.1016/j.amepre.2008.08.030.
Otero M, Solari HG, Schweigmann N: A stochastic population dynamics model for Aedes aegypti: formulation and application to a city with temperate climate. Bull Math Biol. 2006, 68 (8): 1945-1974. 10.1007/s11538-006-9067-y.
Bhatt S, Gething PW, Brady OJ, Messina JP, Farlow AW, Moyes CL, Drake JM, Brownstein JS, Hoen AG, Sankoh O: The global distribution and burden of dengue. Nature. 2013, 496 (7446): 504-10.1038/nature12060.
Jentes ES, Poumerol G, Gershman MD, Hill DR, Lemarchand J, Lewis RF, Staples JE, Tomori O, Wilder-Smith A, Monath TP: The revised global yellow fever risk map and recommendations for vaccination, 2010: consensus of the Informal WHO Working Group on Geographic Risk for Yellow Fever. Lancet Infect Dis. 2011, 11 (8): 622-632. 10.1016/S1473-3099(11)70147-5.
Brady OJ, Gething PW, Bhatt S, Messina JP, Brownstein JS, Hoen AG, Moyes CL, Farlow AW, Scott TW, Hay SI: Refining the global spatial limits of dengue virus transmission by evidence-based consensus. PLoS Negl Trop Dis. 2012, 6 (8): e1760-10.1371/journal.pntd.0001760.
Gratz N: Critical review of the vector status of Aedes albopictus. Med Vet Entomol. 2004, 18 (3): 215-227. 10.1111/j.0269-283X.2004.00513.x.
Lambrechts L, Scott TW, Gubler DJ: Consequences of the expanding global distribution of Aedes albopictus for dengue virus transmission. PLoS Negl Trop Dis. 2010, 4 (5): e646-10.1371/journal.pntd.0000646.
Shroyer D: Aedes albopictus and arboviruses: a concise review of the literature. J Am Mosq Control Assoc. 1986, 2 (4): 424-
Turell M: Horizontal and Vertical Transmission of Viruses by Insect and Tick Vectors, Volume 1. 1988, Boca Raton, FL: CRC Press
Perrau J, Catteau C, Michault A, Parain A, Favier F: Fin 2006, 300 000 personnes avaient été atteintes par le chikungunya. Econ Reun. 2007, 129: 16-17.
Christophers S: Aedes aegypti (L.) the Yellow Fever Mosquito: Its Life History, Bionomics and Structure. 1960, Cambridge: Cambridge University Press, 1
Kuno G: Review of the factors modulating dengue transmission. Epidemiol Rev. 1995, 17 (2): 321-335.
Alto BW, Bettinardi D: Temperature and dengue virus infection in mosquitoes: independent effects on the immature and adult stages. Am J Trop Med Hyg. 2013, 88 (3): 497-505. 10.4269/ajtmh.12-0421.
Yang H, Macoris M, Galvani K, Andrighetti M, Wanderley D: Assessing the effects of temperature on the population of Aedes aegypti, the vector of dengue. Epidemiol Infect. 2009, 137 (8): 1188-10.1017/S0950268809002040.
Nawrocki S, Hawley W: Estimation of the northern limits of distribution of Aedes albopictus in North America. J Am Mosq Control Assoc. 1987, 3 (2): 314-317.
Harrington LC, Buonaccorsi JP, Edman JD, Costero A, Kittayapong P, Clark GG, Scott TW: Analysis of survival of young and old Aedes aegypti (Diptera: Culicidae) from Puerto Rico and Thailand. J Med Entomol. 2001, 38 (4): 537-547. 10.1603/0022-2585-38.4.537.
Harrington LC, Scott TW, Lerdthusnee K, Coleman RC, Costero A, Clark GG, Jones JJ, Kitthawee S, Kittayapong P, Sithiprasasna R: Dispersal of the dengue vector Aedes aegypti within and between rural communities. Am J Trop Med Hyg. 2005, 72 (2): 209-220.
Service MW: Mosquito Ecology: Field Sampling Methods. 1993, London: Elsevier Applied Science, 2
Macdonald G: The analysis of the sporozoite rate. Tropical Diseases Bulletin. 1952, 49 (6): 569-
Styer LM, Carey JR, Wang J-L, Scott TW: Mosquitoes do senesce: departure from the paradigm of constant mortality. Trop Dis Med. 2007, 76 (1): 111-
Lambrechts L, Paaijmans KP, Fansiri T, Carrington LB, Kramer LD, Thomas MB, Scott TW: Impact of daily temperature fluctuations on dengue virus transmission by Aedes aegypti. Proc Natl Acad Sci USA. 2011, 108 (18): 7460-7465. 10.1073/pnas.1101377108.
Carrington LB, Armijos MV, Lambrechts L, Scott TW: Fluctuations at a low mean temperature accelerate dengue virus transmission by Aedes aegypti. PLoS Negl Trop Dis. 2013, 7 (4): e2190-10.1371/journal.pntd.0002190.
Clements AN, Paterson GD: The analysis of mortality and survival rates in wild populations of mosquitoes. J Appl Ecol. 1981, 18 (2): 373-399. 10.2307/2402401.
Degallier N, Servain J, Lucio PS, Hannart A, Durand B, de Souza RN, Ribeiro ZM: The influence of local environment on the aging and mortality of Aedes aegypti (L.): case study in Fortaleza-CE, Brazil. J Vector Ecol. 2012, 37 (2): 428-10.1111/j.1948-7134.2012.00247.x.
Strickman D: Longevity of Aedes aegypti (Diptera: Culicidae) compared in cages and field under ambient conditions in rural Thailand. Southeast Asian J Trop Med Public Health. 2006, 37 (3): 456-462.
Waldock J, Chandra NL, Lelieveld J, Proestos Y, Michael E, Christophides G, Parham PE: The role of environmental variables on Aedes albopictus biology and chikungunya epidemiology. Pathog Glob Health. 2013, 107 (5): 224-241. 10.1179/2047773213Y.0000000100.
Focks DA, Daniels E, Haile DG, Keesling JE: A simulation model of the epidemiology of urban dengue fever: literature analysis, model development, preliminary validation, and samples of simulation results. Am J Trop Med Hyg. 1995, 53 (5): 489-506.
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 (3): 528-538. 10.1111/j.1365-2435.2008.01538.x.
Otero M, Solari H: Stochastic eco-epidemiological model of dengue disease transmission by Aedes aegypti mosquito. Math Biosci. 2010, 223 (1): 32-46. 10.1016/j.mbs.2009.10.005.
Edman JD, Strickman D, Kittayapong P, Scott TW: Female Aedes aegypti (Diptera: Culicidae) in Thailand rarely feed on sugar. J Med Entomol. 1992, 29 (6): 1035-1038.
Harrington LC, Edman JD, Scott TW: Why do female Aedes aegypti (Diptera: Culicidae) feed preferentially and frequently on human blood?. J Med Entomol. 2001, 38 (3): 411-422. 10.1603/0022-2585-38.3.411.
Day JF, Edman JD, Scott TW: Reproductive fitness and survivorship of Aedes aegypti (Diptera: Culicidae) maintained on blood, with field observations from Thailand. J Med Entomol. 1994, 31 (4): 611-617.
Naksathit AT, Scott TW: Effect of female size on fecundity and survivorship of Aedes aegypti fed only human blood versus human blood plus sugar. J Am Mosq Control Assoc. 1998, 14 (2): 148-152.
Styer LM, Minnick SL, Sun AK, Scott TW: Mortality and reproductive dynamics of Aedes aegypti (Diptera: Culicidae) fed human blood. Vector Borne Zoonot Dis. 2007, 7 (1): 86-98. 10.1089/vbz.2007.0216.
Costero A, Edman JD, Clark GG, Scott TW: Life table study of Aedes aegypti (Diptera: Culicidae) in Puerto Rico fed only human blood versus blood plus sugar. J Med Entomol. 1998, 35 (5): 809-813.
Scott TW, Naksathit A, Day JF, Kittayapong P, Edman JD: A fitness advantage for Aedes aegypti and the viruses it transmits when females feed only on human blood. Am J Trop Med Hyg. 1997, 57 (2): 235-239.
GetData Graph Digitizer. getdata-graph-digitizer.com
WorldClim global climate data, data for current conditions.http://www.worldclim.org/current,
Yamar BA, Diallo D, Kebe CM, Dia I, Diallo M: Aspects of bioecology of two Rift Valley Fever Virus vectors in Senegal (West Africa): Aedes vexans and Culex poicilipes (Diptera: Culicidae). J Med Entomol. 2005, 42 (5): 739-750. 10.1603/0022-2585(2005)042[0739:AOBOTR]2.0.CO;2.
David MR, Lourenco-de-Oliveira R, Freitas RM: Container productivity, daily survival rates and dispersal of Aedes aegypti mosquitoes in a high income dengue epidemic neighbourhood of Rio de Janeiro: presumed influence of differential urban structure on mosquito biology. Mem Inst Oswaldo Cruz. 2009, 104 (6): 927-932. 10.1590/S0074-02762009000600019.
Lacroix R, Delatte H, Hue T, Dehecq JS, Reiter P: Adaptation of the BG-Sentinel trap to capture male and female Aedes albopictus mosquitoes. Med Vet Entomol. 2009, 23 (2): 160-162. 10.1111/j.1365-2915.2009.00806.x.
Maciel-de-Freitas R, Codeco CT, Lourenco-de-Oliveira R: Daily survival rates and dispersal of Aedes aegypti females in Rio de Janeiro. Brazil. Am J Trop Med Hyg. 2007, 76 (4): 659-665.
Maciel-de-Freitas R, Marques WA, Peres RC, Cunha SP, de Oliveira RL: Variation in Aedes aegypti (Diptera: Culicidae) container productivity in a slum and a suburban district of Rio de Janeiro during dry and wet seasons. Mem Inst Oswaldo Cruz. 2007, 102 (4): 489-496. 10.1590/S0074-02762007005000056.
Muir LE, Kay BH: Aedes aegypti survival and dispersal estimated by mark-release-recapture in northern Australia. Am J Trop Med Hyg. 1998, 58 (3): 277-282.
Tsuda Y, Takagi M, Wang S, Wang Z, Tang L: Movement of Aedes aegypti (Diptera: Culicidae) released in a small isolated village on Hainan Island, China. J Med Entomol. 2001, 38 (1): 93-98. 10.1603/0022-2585-38.1.93.
Delatte H, Gimonneau G, Triboire A, Fontenille D: Influence of temperature on immature development, survival, longevity, fecundity, and gonotrophic cycles of Aedes albopictus, vector of chikungunya and dengue in the Indian Ocean. J Med Entomol. 2009, 46 (1): 33-41. 10.1603/033.046.0105.
Wilson DL: The analysis of survival (mortality) data: fitting Gompertz, Weibull, and logistic functions. Mech Ageing Dev. 1994, 74 (1–2): 15-33.
Lee ET, Go OT: Survival analysis in public health research. Annu Rev Public Health. 1997, 18: 105-134. 10.1146/annurev.publhealth.18.1.105.
R Development Core Team: R: A language and environment for statistical computing. 2012, Vienna, Austria: R Foundation for Statistical Computing, 2142
Burnham KP, Anderson DR: Model Selection and Multi-Model Inference: A Practical Information-Theoretic Approach. 2002, New York: Springer, 2
Hastie T, Tibshirani R: Generalized Additive Models. Statist Sci. 1986, 1 (2): 297-
Wood SN: Generalized Additive Models: An Introduction with R, vol. 66. 2006, Boca Raton, FL: Chapman & Hall
Buonaccorsi JP, Harrington LC, Edman JD: Estimation and comparison of mosquito survival rates with release-recapture-removal data. J Med Entomol. 2003, 40 (1): 6-17. 10.1603/0022-2585-40.1.6.
Gao J, Zhen Z, Xue J, Huang P, Zhao J, Cao N: Studies on the longevity of adult Aedes (S.) albopictus (Skuse): the longevity of caged females under laboratory conditions. Acta Entomol Sin. 1984, 27: 182-188.
Smith G, Eliason D, Moore C, Ihenacho E: Use of elevated temperatures to kill Aedes albopictus and Ae. aegypti. J Am Mosq Control Assoc. 1988, 4 (4): 557-558.
Parton WJ, Logan JA: A model for diurnal variation in soil and air temperature. Agr Met. 1981, 23: 205-216.
Joy TK, Arik AJ, Corby-Harris V, Johnson AA, Riehle MA: The impact of larval and adult dietary restriction on lifespan, reproduction and growth in the mosquito Aedes aegypti. Exp Gerontol. 2010, 45 (9): 685-690. 10.1016/j.exger.2010.04.009.
Model of adult Aedes albopictus survival/mortality at different temperatures.http://dx.doi.org/10.6084/m9.figshare.865035,
Model of adult Aedes aegypti survival/mortality at different temperatures.http://dx.doi.org/10.6084/m9.figshare.865034,
Grech MG, Luduena A, Francisco A, Almiron WR: Bionomics of Aedes aegypti subpopulations (Diptera: Culicidae) from Argentina. J Vector Ecol. 2010, 35 (2): 277-10.1111/j.1948-7134.2010.00083.x.
Lansdowne C, Hacker CS: The effect of fluctuating temperature and humidity on the adult life table characteristics of five strains of Aedes aegypti. J Med Entomol. 1975, 11 (6): 723-733.
Leisnham P, Sala L, Juliano S: Geographic variation in adult survival and reproductive tactics of the mosquito Aedes albopictus. J Med Entomol. 2008, 45 (2): 210-10.1603/0022-2585(2008)45[210:GVIASA]2.0.CO;2.
Joy TK, Jeffrey Gutierrez EH, Ernst K, Walker KR, Carriere Y, Torabi M, Riehle MA: Aging field collected Aedes aegypti to determine their capacity for dengue transmission in the southwestern United States. PLoS One. 2012, 7 (10): e46946-10.1371/journal.pone.0046946.
Yoon IK, Getis A, Aldstadt J, Rothman AL, Tannitisupawong D, Koenraadt CJ, Fansiri T, Jones JW, Morrison AC, Jarman RG: Fine scale spatiotemporal clustering of dengue virus transmission in children and Aedes aegypti in rural Thai villages. PLoS Negl Trop Dis. 2012, 6 (7): e1730-10.1371/journal.pntd.0001730.
Rowley WA, Graham CL: Effect of temperature and relative humidity on flight ability of female Aedes aegypti. Laboratory Studies of Mosquito Flight. 1967, Fort Detrick, Frederick, MD: Department of the Army
Costero A, Edman JD, Clark GG, Kittayapong P, Scott TW: Survival of starved Aedes aegypti (Diptera: Culicidae) in Puerto Rico and Thailand. J Med Entomol. 1999, 36 (3): 272-276.
Pumpuni C, Knepler J, Craig G: Influence of temperature and larval nutrition on the diapause inducing photoperiod of Aedes albopictus. J Am Mosq Control Assoc. 1992, 8 (3): 223-227.
Hanson SM, Craig GB: Relationship between cold hardiness and supercooling point in Aedes albopictus eggs. J Am Mosq Control Assoc. 1995, 11 (1): 35-38.
Fischer S, Alem IS, De Majo MS, Campos RE, Schweigmann N: Cold season mortality and hatching behavior of Aedes aegypti L. (Diptera: Culicidae) eggs in Buenos Aires City, Argentina. J Vector Ecol. 2011, 36 (1): 94-99. 10.1111/j.1948-7134.2011.00145.x.
Thomas SM, Obermayr U, Fischer D, Kreyling J, Beierkuhnlein C: Low-temperature threshold for egg survival of a post-diapause and non-diapause European aedine strain, Aedes albopictus (Diptera: Culicidae). Parasit Vectors. 2012, 5 (1): 1-7. 10.1186/1756-3305-5-1.
Chan M, Johansson MA: The incubation periods of dengue viruses. PLoS One. 2012, 7 (11): e50972-10.1371/journal.pone.0050972.
Brady OJ, Messina JP, Scott TW, Hay SI: Mapping the global epidemiology of dengue. Dengue and Dengue Hemorrhagic Fever. Edited by: Gubler DJ, Ooi EE, Vasudevan SG, Farrar J. 2014, New York: CABI International, In press, 2
WHO: Dengue: Guidelines For Diagnosis, Treatment, Prevention and Control. 2009, Geneva: World Health Organization
Reiner RC, Perkins TA, Barker CM, Niu T, Chaves LF, Ellis AM, George DB, Le Menach A, Pulliam JR, Bisanzio D: A systematic review of mathematical models of mosquito-borne pathogen transmission: 1970–2010. J Royal Soc Interface. 2013, 10 (81):
Gething PW, Van Boeckel TP, Smith DL, Guerra CA, Patil AP, Snow RW, Hay SI: Modelling the global constraints of temperature on transmission of Plasmodium falciparum and P. vivax. Parasit Vectors. 2011, 4: 92-10.1186/1756-3305-4-92.
Badurdeen S, Valladares DB, Farrar J, Gozzer E, Kroeger A, Kuswara N, Ranzinger SR, Tinh HT, Leite P, Mahendradhata Y: Sharing experiences: towards an evidence based model of dengue surveillance and outbreak response in Latin America and Asia. BMC Public Health. 2013, 13 (1): 607-10.1186/1471-2458-13-607.
We would like to thank Lauren Carrington, Louis Lambrechts, Dalva Wanderley, Maria Macoris, Ephantus J. Muturi and Barry Alto for kindly sharing primary survival data.
O.J.B. is funded by a BBSRC studentship. DMP is funded by a Sir Richard Southwood Graduate Scholarship from the Department of Zoology at the University of Oxford. P.W.G is a Medical Research Council (UK) Career Development Fellow (#K00669X) and receives support from the Bill and Melinda Gates Foundation (#OPP1068048) which also supports S.B. N.G. is funded by a grant from the Bill & Melinda Gates Foundation (OPP1053338). S.I.H. is funded by a Senior Research Fellowship from the Wellcome Trust (095066). This study was partially funded by EU grant 21803 IDAMS (http://www.idams.eu). The contents of this publication are the sole responsibility of the authors and don't necessarily reflect the views of the European Commission. S.I.H. C.A.G. and T.W.S. also acknowledge funding support from the RAPIDD program of the Science & Technology Directorate, Department of Homeland Security, and the Fogarty International Center, National Institutes of Health.
Disclaimer: The findings and conclusions are those of the authors and do not necessarily represent the official position of the Centers for Disease Control and Prevention.
OJB designed the experiment, wrote the manuscript and collected and analysed the data. MAJ, DLS, TWS, PWG and SIH also helped conceive and design the experiments. MAJ, SB, NG, DMP, DLS and TWS helped with data analysis and interpretation. CAG, HD, MGG, PL, RM and LMS collected data. All authors were involved in drafting and revising the manuscript and all authors approved the final version.