Parasit VectorsParasit VectorsParasites & Vectors1756-3305BioMed Central2433072038672191756-3305-6-35110.1186/1756-3305-6-351ResearchModelling adult Aedes aegypti and Aedes albopictus survival at different temperatures in laboratory and field settingsBradyOliver J1oliver.brady@zoo.ox.ac.ukJohanssonMichael A2mjohansson@cdc.govGuerraCarlos A1carlos.guerraloaiza@gmail.comBhattSamir1samir.bhatt@zoo.ox.ac.ukGoldingNick1nick.golding@zoo.ox.ac.ukPigottDavid M1david.pigott@zoo.ox.ac.ukDelatteHélène3Helene.delatte@cirad.frGrechMarta G4grechmarta@gmail.comLeisnhamPaul T5leisnham@umd.eduMaciel-de-FreitasRafael6freitas@ioc.fiocruz.brStyerLinda M7lms24@health.state.ny.usSmithDavid L810dlsmith@jhsph.eduScottThomas W910twscott@ucdavis.eduGethingPeter W1peter.gething@zoo.ox.ac.ukHaySimon I110simon.hay@zoo.ox.ac.ukSpatial Ecology and Epidemiology Group, Department of Zoology, University of Oxford, Tinbergen Building, South Parks Road, Oxford, UKDengue Branch, Division of Vector-Borne Diseases, Centers for Disease Control and Prevention, San Juan, Puerto Rico, USACIRAD, UMR PVBMT, 97410, Saint Piarre, la Réunion, FranceLaboratorio de Investigaciones en Ecología y Sistemática Animal (LIESA), Universidad Nacional de la Patagonia San Juan Bosco, FCN-Sede, Esquel, Chubut, ArgentinaDepartment of Environmental Science & Technology, University of Maryland, College Park, MD 20742, USALaboratório de Transmissores de Hematozoários, Instituto Oswaldo Cruz, Fiocruz, Rio de Janeiro, RJ, BrazilWadsworth Center, New York State Department of Health, Albany, NY 12208, USADepartment of Epidemiology, Johns Hopkins Bloomberg School of Public Health, Baltimore MD, USADepartment of Entomology, University of California Davis, Davis, CA, USAFogarty International Center, National Institutes of Health, Bethesda, MD 20892, USA201312122013635135118920136122013Copyright © 2013 Brady et al.; licensee BioMed Central Ltd.2013Brady et al.; licensee BioMed Central Ltd.This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.Background

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.

Methods

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.

Results

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.

Conclusions

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.

AedesSurvivalTemperatureMortalityLongevityModellingMark-release-recaptureDengueTransmissionGeneralised additive models
Background

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 [17]. 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 [23]. 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 [24] and sensitivities of spatial mosquito sampling limit the level of survival detail that can be elucidated through the majority of MRR experiments [25]. 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 [26]. 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 [31]. 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 [32]. 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 [21]. 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 [33].

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.

MethodsData collection

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 [44]. 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 [45]. 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.

Summary of laboratory and field data

  Ae. albopictus laboratory data Ae. albopictus MRR data Ae. aegypti laboratory data Ae. aegypti MRR data
Number of experiments
210
9
141
50
Mean temperature (oC) (±SD)
25.9 (22.3-29.5)
20.3 (17.0-23.7)
25.5 (21.9-29.1)
25.5 (22.9-28.3)
Minimum/Maximum temperature (oC)
15/35
16.7/26
10/35
20.7/30.1
Median number of mosquitoes observed/released (IQR)
29 (15–40)
552 (249–1007)
70 (25–382)
602 (493–798)
Dietary regime (%)
Blood: 7.1
-
Blood: 3.0
-
Sugar: 1.0
Sugar: 38.2
Blood + Sugar: 91.9
Blood + Sugar: 58.8
Age at release (Days) (±SD)-3.2 (0.6-5.8)-3.6 (1.3-5.9)

SD = Standard deviation, IQR = Interquartile range, MRR = Mark-release-recapture.

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 [56] 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) [57] of each model to the model with the lowest AIC:

RLj=eAICminAICj2

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 [58]. 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:

Sij~BinomialNij,pij

logitpij=fDi,Ti+ϵj+ϵd

ϵj~N0,θj2

ϵd~N0,θd2

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 [59]. 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 [60].

For the model to fit biologically appropriate responses, additional data defining the limits of prediction were required. Observations from Christophers [18], 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 [60]. 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.

Schematic overview of the methods. Green parallelograms indicate input data, orange rectangles show processing or modelling steps, blue diamonds show decision steps and red boxes show output analysis and models (dotted shows intermediate, unbroken line shows final outputs). MRR = Mark-release-recapture.

ResultsParametric 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).

Relative likelihood of four different parametric models for Aedes adult female survival data over a range of constant temperatures. The models included are i) a two parameter Log-logistic model (shape and scale), ii) a two parameter Gompertz model (shape and rate), iii) a one parameter Exponential model (rate) and iv) a two parameter Weibull model (shape and scale).

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.

Evaluation of parametric and non-parametric model fit to laboratory data

ModelMedian increase in AIC relative to GAM
Optimal model choice (%)
( Ae. aegypti / Ae. albopictus )( Ae. aegypti / Ae. albopictus )
Log-Logistic
320.40/22.620
4.4/1.5
Gompertz
304.83/0.005
2.2/21.9
Exponential
302.89/2.969
1.5/13.8
Weibull
302.90/0.023
0.7/13.3
GAM-91.1/49.5

AIC = Akaike information criterion, GAM = Generalised additive model.

Examples of parametric and non-parametric model fit. Open circles show Ae. aegypti survival data under controlled laboratory conditions from Joy et al.[64](A) and Yang et al.[21](B), two experiments that show contrasting survival curve shape. Parametric models are shown as dashed lines and the non-parametric GAM is shown as a solid orange line.

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 distribution of adult female Aedes aegypti and Aedes albopictus survival across a range of temperatures under laboratory conditions (A and B) and field conditions (C and D). Colours from red to yellow show survival from 100% - 1% of the population remaining. Grey indicates <1% of the population remaining. Dotted blue lines show the limits for 50% and 95% of the original population remaining.

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.

The distribution of uncertainty of the laboratory model prediction. Colours from blue to beige show the interquartile range (IQR) in predictions from 200 bootstrap runs of the laboratory model (A and B). This uncertainty is then combined with the field data uncertainty quantified by 200 bootstrap runs of MRR data to give the IQR predictions for the field survival model (C and D). Red dotted lines of the 50% and 95% of the population remaining are added for reference.

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].

Discussion

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 [70] or measuring mosquito infection rates [71].

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 [73], 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 [78]. 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 [12] with an estimated 390 million new infections per year [10], 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 [83].

Conclusions

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.

Abbreviations

DENV: Dengue virus; MRR: Mark-release-recapture; AIC: Akaike’s information criterion; RL: Relative likelihood; GAM: Generalised additive model.

Competing interests

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.

Authors’ contributions

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.

Supplementary MaterialAdditional file 1

Data bibliography and summary statistics.

Click here for file

Acknowledgements

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.

Garrett-JonesCShidrawiGMalaria vectorial capacity of a population of Anopheles gambiae: an exercise in epidemiological entomologyBull World Health Organ19694045315306719ReisenWKMahmoodFParveenTAnopheles culicifacies Giles: a release-recapture experiment with cohorts of known age with implications for malaria epidemiology and genetical control in PakistanTrans R Soc Trop Med Hyg198074330731710.1016/0035-9203(80)90089-97434425MacdonaldGEpidemiological basis of malaria controlBull World Health Organ1956153–561313404439ParhamPEMichaelEModeling the effects of weather and climate change on malaria transmissionEnviron Health Perspect2010118562020435552BarbazanPGuiserixMBoonyuanWTuntaprasartWPontierDGonzalezJPModelling the effect of temperature on transmission of dengueMed Vet Entomol2010241667310.1111/j.1365-2915.2009.00848.x20377733SmithDLBattleKEHaySIBarkerCMScottTWMcKenzieFERoss, Macdonald, and a theory for the dynamics and control of mosquito-transmitted pathogensPLoS Pathog201284e100258810.1371/journal.ppat.100258822496640LundeTMBayohMNLindtjørnBHow malaria models relate temperature to malaria transmissionParasite Vector201362010.1186/1756-3305-6-20GageKLBurkotTREisenRJHayesEBClimate and vectorborne diseasesAm J Prev Med200835543645010.1016/j.amepre.2008.08.03018929970OteroMSolariHGSchweigmannNA stochastic population dynamics model for Aedes aegypti: formulation and application to a city with temperate climateBull Math Biol20066881945197410.1007/s11538-006-9067-y16832731BhattSGethingPWBradyOJMessinaJPFarlowAWMoyesCLDrakeJMBrownsteinJSHoenAGSankohOThe global distribution and burden of dengueNature2013496744650410.1038/nature1206023563266JentesESPoumerolGGershmanMDHillDRLemarchandJLewisRFStaplesJETomoriOWilder-SmithAMonathTPThe revised global yellow fever risk map and recommendations for vaccination, 2010: consensus of the Informal WHO Working Group on Geographic Risk for Yellow FeverLancet Infect Dis201111862263210.1016/S1473-3099(11)70147-521798462BradyOJGethingPWBhattSMessinaJPBrownsteinJSHoenAGMoyesCLFarlowAWScottTWHaySIRefining the global spatial limits of dengue virus transmission by evidence-based consensusPLoS Negl Trop Dis201268e176010.1371/journal.pntd.000176022880140GratzNCritical review of the vector status of Aedes albopictusMed Vet Entomol200418321522710.1111/j.0269-283X.2004.00513.x15347388LambrechtsLScottTWGublerDJConsequences of the expanding global distribution of Aedes albopictus for dengue virus transmissionPLoS Negl Trop Dis201045e64610.1371/journal.pntd.000064620520794ShroyerDAedes albopictus and arboviruses: a concise review of the literatureJ Am Mosq Control Assoc1986244243333791TurellMHorizontal and Vertical Transmission of Viruses by Insect and Tick Vectors, Volume 11988Boca Raton, FL: CRC PressPerrauJCatteauCMichaultAParainAFavierFFin 2006, 300 000 personnes avaient été atteintes par le chikungunyaEcon Reun20071291617ChristophersSAedes aegypti (L.) the Yellow Fever Mosquito: Its Life History, Bionomics and Structure19601Cambridge: Cambridge University PressKunoGReview of the factors modulating dengue transmissionEpidemiol Rev19951723213358654514AltoBWBettinardiDTemperature and dengue virus infection in mosquitoes: independent effects on the immature and adult stagesAm J Trop Med Hyg201388349750510.4269/ajtmh.12-042123382163YangHMacorisMGalvaniKAndrighettiMWanderleyDAssessing the effects of temperature on the population of Aedes aegypti, the vector of dengueEpidemiol Infect20091378118810.1017/S095026880900204019192322NawrockiSHawleyWEstimation of the northern limits of distribution of Aedes albopictus in North AmericaJ Am Mosq Control Assoc1987323143173504917HarringtonLCBuonaccorsiJPEdmanJDCosteroAKittayapongPClarkGGScottTWAnalysis of survival of young and old Aedes aegypti (Diptera: Culicidae) from Puerto Rico and ThailandJ Med Entomol200138453754710.1603/0022-2585-38.4.53711476334HarringtonLCScottTWLerdthusneeKColemanRCCosteroAClarkGGJonesJJKitthaweeSKittayapongPSithiprasasnaRDispersal of the dengue vector Aedes aegypti within and between rural communitiesAm J Trop Med Hyg200572220922015741559ServiceMWMosquito Ecology: Field Sampling Methods19932London: Elsevier Applied ScienceMacdonaldGThe analysis of the sporozoite rateTropical Diseases Bulletin195249656914958825StyerLMCareyJRWangJ-LScottTWMosquitoes do senesce: departure from the paradigm of constant mortalityTrop Dis Med2007761111LambrechtsLPaaijmansKPFansiriTCarringtonLBKramerLDThomasMBScottTWImpact of daily temperature fluctuations on dengue virus transmission by Aedes aegyptiProc Natl Acad Sci USA2011108187460746510.1073/pnas.110137710821502510CarringtonLBArmijosMVLambrechtsLScottTWFluctuations at a low mean temperature accelerate dengue virus transmission by Aedes aegyptiPLoS Negl Trop Dis201374e219010.1371/journal.pntd.000219023638208ClementsANPatersonGDThe analysis of mortality and survival rates in wild populations of mosquitoesJ Appl Ecol198118237339910.2307/2402401DegallierNServainJLucioPSHannartADurandBde SouzaRNRibeiroZMThe influence of local environment on the aging and mortality of Aedes aegypti (L.): case study in Fortaleza-CE, BrazilJ Vector Ecol201237242810.1111/j.1948-7134.2012.00247.x23181868StrickmanDLongevity of Aedes aegypti (Diptera: Culicidae) compared in cages and field under ambient conditions in rural ThailandSoutheast Asian J Trop Med Public Health200637345646217120964WaldockJChandraNLLelieveldJProestosYMichaelEChristophidesGParhamPEThe role of environmental variables on Aedes albopictus biology and chikungunya epidemiologyPathog Glob Health2013107522424110.1179/2047773213Y.000000010023916332FocksDADanielsEHaileDGKeeslingJEA simulation model of the epidemiology of urban dengue fever: literature analysis, model development, preliminary validation, and samples of simulation resultsAm J Trop Med Hyg19955354895067485707KearneyMPorterWPWilliamsCRitchieSHoffmannAAIntegrating biophysical models and evolutionary theory to predict climatic impacts on species’ ranges: the dengue mosquito Aedes aegypti in AustraliaFunct Ecol200923352853810.1111/j.1365-2435.2008.01538.xOteroMSolariHStochastic eco-epidemiological model of dengue disease transmission by Aedes aegypti mosquitoMath Biosci20102231324610.1016/j.mbs.2009.10.00519861133EdmanJDStrickmanDKittayapongPScottTWFemale Aedes aegypti (Diptera: Culicidae) in Thailand rarely feed on sugarJ Med Entomol1992296103510381460619HarringtonLCEdmanJDScottTWWhy do female Aedes aegypti (Diptera: Culicidae) feed preferentially and frequently on human blood?J Med Entomol200138341142210.1603/0022-2585-38.3.41111372967DayJFEdmanJDScottTWReproductive fitness and survivorship of Aedes aegypti (Diptera: Culicidae) maintained on blood, with field observations from ThailandJ Med Entomol19943146116177932609NaksathitATScottTWEffect of female size on fecundity and survivorship of Aedes aegypti fed only human blood versus human blood plus sugarJ Am Mosq Control Assoc19981421481529673914StyerLMMinnickSLSunAKScottTWMortality and reproductive dynamics of Aedes aegypti (Diptera: Culicidae) fed human bloodVector Borne Zoonot Dis200771869810.1089/vbz.2007.0216CosteroAEdmanJDClarkGGScottTWLife table study of Aedes aegypti (Diptera: Culicidae) in Puerto Rico fed only human blood versus blood plus sugarJ Med Entomol19983558098139775613ScottTWNaksathitADayJFKittayapongPEdmanJDA fitness advantage for Aedes aegypti and the viruses it transmits when females feed only on human bloodAm J Trop Med Hyg19975722352399288822GetData Graph Digitizergetdata-graph-digitizer.comWorldClim global climate data, data for current conditionshttp://www.worldclim.org/currentYamarBADialloDKebeCMDiaIDialloMAspects of bioecology of two Rift Valley Fever Virus vectors in Senegal (West Africa): Aedes vexans and Culex poicilipes (Diptera: Culicidae)J Med Entomol200542573975010.1603/0022-2585(2005)042[0739:AOBOTR]2.0.CO;216363157DavidMRLourenco-de-OliveiraRFreitasRMContainer 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 biologyMem Inst Oswaldo Cruz2009104692793210.1590/S0074-0276200900060001919876569LacroixRDelatteHHueTDehecqJSReiterPAdaptation of the BG-Sentinel trap to capture male and female Aedes albopictus mosquitoesMed Vet Entomol200923216016210.1111/j.1365-2915.2009.00806.x19493196Maciel-de-FreitasRCodecoCTLourenco-de-OliveiraRDaily survival rates and dispersal of Aedes aegypti females in Rio de JaneiroBrazil. Am J Trop Med Hyg2007764659665Maciel-de-FreitasRMarquesWAPeresRCCunhaSPde OliveiraRLVariation in Aedes aegypti (Diptera: Culicidae) container productivity in a slum and a suburban district of Rio de Janeiro during dry and wet seasonsMem Inst Oswaldo Cruz2007102448949610.1590/S0074-0276200700500005617612770MuirLEKayBHAedes aegypti survival and dispersal estimated by mark-release-recapture in northern AustraliaAm J Trop Med Hyg19985832772829546403TsudaYTakagiMWangSWangZTangLMovement of Aedes aegypti (Diptera: Culicidae) released in a small isolated village on Hainan Island, ChinaJ Med Entomol2001381939810.1603/0022-2585-38.1.9311268697DelatteHGimonneauGTriboireAFontenilleDInfluence of temperature on immature development, survival, longevity, fecundity, and gonotrophic cycles of Aedes albopictus, vector of chikungunya and dengue in the Indian OceanJ Med Entomol2009461334110.1603/033.046.010519198515WilsonDLThe analysis of survival (mortality) data: fitting Gompertz, Weibull, and logistic functionsMech Ageing Dev1994741–215337934205LeeETGoOTSurvival analysis in public health researchAnnu Rev Public Health19971810513410.1146/annurev.publhealth.18.1.1059143714R Development Core TeamR: A language and environment for statistical computing20122142Vienna, Austria: R Foundation for Statistical ComputingBurnhamKPAndersonDRModel Selection and Multi-Model Inference: A Practical Information-Theoretic Approach20022New York: SpringerHastieTTibshiraniRGeneralized Additive ModelsStatist Sci198612297WoodSNGeneralized Additive Models: An Introduction with R, vol. 662006Boca Raton, FL: Chapman & HallBuonaccorsiJPHarringtonLCEdmanJDEstimation and comparison of mosquito survival rates with release-recapture-removal dataJ Med Entomol200340161710.1603/0022-2585-40.1.612597647GaoJZhenZXueJHuangPZhaoJCaoNStudies on the longevity of adult Aedes (S.) albopictus (Skuse): the longevity of caged females under laboratory conditionsActa Entomol Sin198427182188SmithGEliasonDMooreCIhenachoEUse of elevated temperatures to kill Aedes albopictus and Ae. aegyptiJ Am Mosq Control Assoc1988445575583225579PartonWJLoganJAA model for diurnal variation in soil and air temperatureAgr Met198123205216JoyTKArikAJCorby-HarrisVJohnsonAARiehleMAThe impact of larval and adult dietary restriction on lifespan, reproduction and growth in the mosquito Aedes aegyptiExp Gerontol201045968569010.1016/j.exger.2010.04.00920451597Model of adult Aedes albopictus survival/mortality at different temperatureshttp://dx.doi.org/10.6084/m9.figshare.865035Model of adult Aedes aegypti survival/mortality at different temperatureshttp://dx.doi.org/10.6084/m9.figshare.865034GrechMGLuduenaAFranciscoAAlmironWRBionomics of Aedes aegypti subpopulations (Diptera: Culicidae) from ArgentinaJ Vector Ecol201035227710.1111/j.1948-7134.2010.00083.x21175932LansdowneCHackerCSThe effect of fluctuating temperature and humidity on the adult life table characteristics of five strains of Aedes aegyptiJ Med Entomol19751167237331123832LeisnhamPSalaLJulianoSGeographic variation in adult survival and reproductive tactics of the mosquito Aedes albopictusJ Med Entomol200845221010.1603/0022-2585(2008)45[210:GVIASA]2.0.CO;218402136JoyTKJeffrey GutierrezEHErnstKWalkerKRCarriereYTorabiMRiehleMAAging field collected Aedes aegypti to determine their capacity for dengue transmission in the southwestern United StatesPLoS One2012710e4694610.1371/journal.pone.004694623077536YoonIKGetisAAldstadtJRothmanALTannitisupawongDKoenraadtCJFansiriTJonesJWMorrisonACJarmanRGFine scale spatiotemporal clustering of dengue virus transmission in children and Aedes aegypti in rural Thai villagesPLoS Negl Trop Dis201267e173010.1371/journal.pntd.000173022816001RowleyWAGrahamCLEffect of temperature and relative humidity on flight ability of female Aedes aegyptiLaboratory Studies of Mosquito Flight1967Fort Detrick, Frederick, MD: Department of the ArmyCosteroAEdmanJDClarkGGKittayapongPScottTWSurvival of starved Aedes aegypti (Diptera: Culicidae) in Puerto Rico and ThailandJ Med Entomol199936327227610337096PumpuniCKneplerJCraigGJrInfluence of temperature and larval nutrition on the diapause inducing photoperiod of Aedes albopictusJ Am Mosq Control Assoc1992832232271402857HansonSMCraigGBJrRelationship between cold hardiness and supercooling point in Aedes albopictus eggsJ Am Mosq Control Assoc199511135387616187FischerSAlemISDe MajoMSCamposRESchweigmannNCold season mortality and hatching behavior of Aedes aegypti L. (Diptera: Culicidae) eggs in Buenos Aires City, ArgentinaJ Vector Ecol2011361949910.1111/j.1948-7134.2011.00145.x21635646ThomasSMObermayrUFischerDKreylingJBeierkuhnleinCLow-temperature threshold for egg survival of a post-diapause and non-diapause European aedine strain, Aedes albopictus (Diptera: Culicidae)Parasit Vectors2012511710.1186/1756-3305-5-122212459ChanMJohanssonMAThe incubation periods of dengue virusesPLoS One2012711e5097210.1371/journal.pone.005097223226436BradyOJMessinaJPScottTWHaySIGubler DJ, Ooi EE, Vasudevan SG, Farrar JMapping the global epidemiology of dengueDengue and Dengue Hemorrhagic Fever20142New York: CABI InternationalIn pressWHODengue: Guidelines For Diagnosis, Treatment, Prevention and Control2009Geneva: World Health OrganizationReinerRCPerkinsTABarkerCMNiuTChavesLFEllisAMGeorgeDBLe MenachAPulliamJRBisanzioDA systematic review of mathematical models of mosquito-borne pathogen transmission: 1970–2010J Royal Soc Interface20131081GethingPWVan BoeckelTPSmithDLGuerraCAPatilAPSnowRWHaySIModelling the global constraints of temperature on transmission of Plasmodium falciparum and P. vivaxParasit Vectors201149210.1186/1756-3305-4-9221615906BadurdeenSValladaresDBFarrarJGozzerEKroegerAKuswaraNRanzingerSRTinhHTLeitePMahendradhataYSharing experiences: towards an evidence based model of dengue surveillance and outbreak response in Latin America and AsiaBMC Public Health201313160710.1186/1471-2458-13-60723800243