Age- and bite-structured models for vector-borne diseases

Incorporating adult age into mosquito population models: Implications for predicting abundances in changing climates

Mosquito-borne diseases (MBDs) pose increasing threats under future climate change scenarios and an understanding of mosquito population dynamics is pivotal to predicting future risk of MBDs. Most models that describe mosquito population dynamics often assume that adult life-history is independent of adult age and yet mosquito senescence is known to affect mosquito mortality, fecundity and other key biological traits. Despite this, little is known about the effects of adult age at the level of the mosquito population, especially under varying temperature scenarios. We develop a stage-structured delayed differential equations (DDEs) model incorporating the effects of the abiotic environment and adult age to shed light on the complex interactions between age, temperature, and mosquito population dynamics. Taking Culex pipiens, a major vector of West Nile Virus, as our study species our results show that failing to consider mosquito senescence can lead to underestimates of future mosquito abundances predicted under climate change scenarios. We also find that the age-dependent mechanisms combined with the effects of density-dependent mortality on the immature stages can result in mosquito abundances decreasing at extreme temperatures. With our work, we underscore the need for more studies to consider the effects of mosquito age. Not accounting for senescence can compromise the accuracy of abundance estimates and has implications for predicting the risk of future MBD outbreaks.

Aedes albopictus is not an arbovirus aficionado when feeding on cynomolgus macaques or squirrel monkeys

Viruses transmitted by Aedes mosquitoes (e.g., dengue [DENV], Zika [ZIKV]) have demonstrated high potential to spill over from their ancestral, sylvatic cycles in non-human primates to establish transmission in humans. Epidemiological models require accurate knowledge of the contact structure between hosts and vectors, which is highly sensitive to any impacts of virus infection in mosquitoes or hosts on mosquito feeding behavior. Current evidence for whether these viruses affect vector behavior is mixed. Here we leveraged a study on sylvatic DENV-2 and ZIKV transmission between two species of monkey and Aedes albopictus to determine whether virus infection of either host or vector alters vector feeding behavior. Engorgement rates varied from 0% to 100%, but this was not driven by vector nor host infection, but rather by the individual host, host species, and host body temperature. This study highlights the importance of incorporating individual-level heterogeneity of vector biting in arbovirus transmission models.

On the necessity of accounting for age structure in human malaria transmission modeling

Malaria is one of the most common mosquito-borne diseases widespread in tropical and subtropical regions, causing thousands of deaths every year in the world. In a previous paper, we formulated an age-structured model containing three structural variables: (i) the chronological age of human and mosquito populations, (ii) the time since they are infected, and (iii) humans waning immunity (i.e. the progressive loss of protective antibodies after recovery). In the present paper, we expand the analysis of this age-structured model and focus on the derivation of entomological and epidemiological results commonly used in the literature, following the works of Smith and McKenzie. We generalize their results to the age-structured case. In order to quantify the impact of neglecting structuring variables such as chronological age, we assigned values from the literature to our model parameters. While some parameters values are readily accessible from the literature, at least those about the human population, the parameters concerning mosquitoes are less commonly documented and the values of a number of them (e.g. mosquito survival in the presence or in absence of infection) can be discussed extensively. Our analysis, informed by parameter values from the literature, demonstrates that overlooking those structural variables of human and mosquito populations may result in inaccurate epidemiological predictions and suboptimal control strategies. We highlight the epidemiological implications of these findings and emphasize the necessity of considering age structure in future malaria control programs.

Seasonal Patterns of Resource Use Within Natural Populations of Burying Beetles

For organisms in temperate environments, seasonal variation in resource availability and weather conditions exert fluctuating selection pressures on survival and fitness, resulting in diverse adaptive responses. By manipulating resource availability on a local spatial scale, we studied seasonal patterns of resource use within natural populations of burying beetles Nicrophorus vespilloides in a Norfolk woodland. Burying beetles are necrophagous insects that breed on vertebrate carcasses. They are active in Europe between April and October, after which they burrow into the soil and overwinter. Using breeding and chemical analyses, we compared the fecundity and physiological state of beetles that differed in their seasonal resource use. We found seasonal variation in carrion use by wild burying beetles and correlated differences in their reproductive success and cuticular hydrocarbon profiles. Our results provide novel insight into the seasonal correlates of behaviour, physiology and life history in burying beetles.

Aedes albopictus is not an arbovirus aficionado - Impacts of sylvatic flavivirus infection in vectors and hosts on mosquito engorgement on non-human primates

The contact structure between vertebrate hosts and arthropod vectors plays a key role in the spread of arthropod-borne viruses (arboviruses); thus, it is important to determine whether arbovirus infection of either host or vector alters vector feeding behavior. Here we leveraged a study of the replication dynamics of two arboviruses isolated from their ancestral cycles in paleotropical forests, sylvatic dengue-2 (DENV-2) and Zika (ZIKV), in one non-human primate (NHP) species from the paleotropics (cynomolgus macaques, Macaca fascicularis) and one from the neotropics (squirrel monkeys, Saimiri boliviensis) to test the effect of both vector and host infection with each virus on completion of blood feeding (engorgement) of the mosquito Aedes albopictus. Although mosquitoes were starved and given no choice of hosts, engorgement rates varied dramatically, from 0% to 100%. While neither vector nor host infection systematically affected engorgement, NHP species and body temperature at the time of feeding did. We also interrogated the effect of repeated mosquito bites on cytokine expression and found that epidermal growth factor (EGF) and macrophage migration inhibitory factor (MIF) concentrations were dynamically associated with exposure to mosquito bites. This study highlights the importance of incorporating individual-level heterogeneity of vector biting in arbovirus transmission models.

An analytically tractable, age-structured model of the impact of vector control on mosquito-transmitted infections

Vector control is a vital tool utilised by malaria control and elimination programmes worldwide, and as such it is important that we can accurately quantify the expected public health impact of these methods. There are very few previous models that consider vector-control-induced changes in the age-structure of the vector population and the resulting impact on transmission. We analytically derive the steady-state solution of a novel age-structured deterministic compartmental model describing the mosquito feeding cycle, with mosquito age represented discretely by parity-the number of cycles (or successful bloodmeals) completed. Our key model output comprises an explicit, analytically tractable solution that can be used to directly quantify key transmission statistics, such as the effective reproductive ratio under control, Rc, and investigate the age-structured impact of vector control. Application of this model reinforces current knowledge that adult-acting interventions, such as indoor residual spraying of insecticides (IRS) or long-lasting insecticidal nets (LLINs), can be highly effective at reducing transmission, due to the dual effects of repelling and killing mosquitoes. We also demonstrate how larval measures can be implemented in addition to adult-acting measures to reduce Rc and mitigate the impact of waning insecticidal efficacy, as well as how mid-ranges of LLIN coverage are likely to experience the largest effect of reduced net integrity on transmission. We conclude that whilst well-maintained adult-acting vector control measures are substantially more effective than larval-based interventions, incorporating larval control in existing LLIN or IRS programmes could substantially reduce transmission and help mitigate any waning effects of adult-acting measures.

Mosquito ageing modulates the development, virulence and transmission potential of pathogens

Host age variation is a striking source of heterogeneity that can shape the evolution and transmission dynamic of pathogens. Compared with vertebrate systems, our understanding of the impact of host age on invertebrate-pathogen interactions remains limited. We examined the influence of mosquito age on key life-history traits driving human malaria transmission. Females of Anopheles coluzzii, a major malaria vector, belonging to three age classes (4-, 8- and 12-day-old), were experimentally infected with Plasmodium falciparum field isolates. Our findings revealed reduced competence in 12-day-old mosquitoes, characterized by lower oocyst/sporozoite rates and intensities compared with younger mosquitoes. Despite shorter median longevities in older age classes, infected 12-day-old mosquitoes exhibited improved survival, suggesting that the infection might act as a fountain of youth for older mosquitoes specifically. The timing of sporozoite appearance in the salivary glands remained consistent across mosquito age classes, with an extrinsic incubation period of approximately 13 days. Integrating these results into an epidemiological model revealed a lower vectorial capacity for older mosquitoes compared with younger ones, albeit still substantial owing to extended longevity in the presence of infection. Considering age heterogeneity provides valuable insights for ecological and epidemiological studies, informing targeted control strategies to mitigate pathogen transmission.

Estimating malaria transmission risk through surveillance of human-vector interactions in northern Ghana

BACKGROUND

Vector bionomics are important aspects of vector-borne disease control programs. Mosquito-biting risks are affected by environmental, mosquito behavior and human factors, which are important for assessing exposure risk and intervention impacts. This study estimated malaria transmission risk based on vector-human interactions in northern Ghana, where indoor residual spraying (IRS) and insecticide-treated nets (ITNs) have been deployed.

METHODS

Indoor and outdoor human biting rates (HBRs) were measured using monthly human landing catches (HLCs) from June 2017 to April 2019. Mosquitoes collected were identified to species level, and Anopheles gambiae sensu lato (An. gambiae s.l.) samples were examined for parity and infectivity. The HBRs were adjusted using mosquito parity and human behavioral observations.

RESULTS

Anopheles gambiae was the main vector species in the IRS (81%) and control (83%) communities. Indoor and outdoor HBRs were similar in both the IRS intervention (10.6 vs. 11.3 bites per person per night [b/p/n]; z = -0.33, P = 0.745) and control communities (18.8 vs. 16.4 b/p/n; z = 1.57, P = 0.115). The mean proportion of parous An. gambiae s.l. was lower in IRS communities (44.6%) than in control communities (71.7%). After adjusting for human behavior observations and parity, the combined effect of IRS and ITN utilization (IRS: 37.8%; control: 57.3%) on reducing malaria transmission risk was 58% in IRS + ITN communities and 27% in control communities with ITNs alone (z = -4.07, P < 0.001). However, this also revealed that about 41% and 31% of outdoor adjusted bites in IRS and control communities respectively, occurred before bed time (10:00 pm). The mean directly measured annual entomologic inoculation rates (EIRs) during the study were 6.1 infective bites per person per year (ib/p/yr) for IRS communities and 16.3 ib/p/yr for control communities. After considering vector survival and observed human behavior, the estimated EIR for IRS communities was 1.8 ib/p/yr, which represents about a 70% overestimation of risk compared to the directly measured EIR; for control communities, it was 13.6 ib/p/yr (16% overestimation).

CONCLUSION

Indoor residual spraying significantly impacted entomological indicators of malaria transmission. The results of this study indicate that vector bionomics alone do not provide an accurate assessment of malaria transmission exposure risk. By accounting for human behavior parameters, we found that high coverage of ITNs alone had less impact on malaria transmission indices than combining ITNs with IRS, likely due to observed low net use. Reinforcing effective communication for behavioral change in net use and IRS could further reduce malaria transmission.

How do i bite thee? let me count the ways: Exploring the implications of individual biting habits of Aedes aegypti for dengue transmission

In models of mosquito-borne transmission, the mosquito biting rate is an influential parameter, and understanding the heterogeneity of the process of biting is important, as biting is usually assumed to be relatively homogeneous across individuals, with time-between-bites described by an exponentially distributed process. However, these assumptions have not been addressed through laboratory experimentation. We experimentally investigated the daily biting habits of Ae. aegypti at three temperatures (24°C, 28°C, and 32°C) and determined that there was individual heterogeneity in biting habits (number of bites, timing of bites, etc.). We further explored the consequences of biting heterogeneity using an individual-based model designed to examine whether a particular biting profile determines whether a mosquito is more or less likely to 1) become exposed given a single index case of dengue (DENV) and 2) transmit to a susceptible human individual. Our experimental results indicate that there is heterogeneity among individuals and among temperature treatments. We further show that this results in altered probabilities of transmission of DENV to and from individual mosquitoes based on biting profiles. While current model representation of biting may work under some conditions, it might not uniformly be the best fit for this process. Our data also confirm that biting is a non-monotonic process with temperatures around 28°C being optimum.

Omitting age-dependent mosquito mortality in malaria models underestimates the effectiveness of insecticide-treated nets

Mathematical models of vector-borne infections, including malaria, often assume age-independent mortality rates of vectors, despite evidence that many insects senesce. In this study we present survival data on insecticide-resistant Anopheles gambiae s.l. from experiments in Côte d’Ivoire. We fit a constant mortality function and two age-dependent functions (logistic and Gompertz) to the data from mosquitoes exposed (treated) and not exposed (control) to insecticide-treated nets (ITNs), to establish biologically realistic survival functions. This enables us to explore the effects of insecticide exposure on mosquito mortality rates, and the extent to which insecticide resistance might impact the effectiveness of ITNs. We investigate this by calculating the expected number of infectious bites a mosquito will take in its lifetime, and by extension the vectorial capacity. Our results show that the predicted vectorial capacity is substantially lower in mosquitoes exposed to ITNs, despite the mosquitoes in the experiment being highly insecticide-resistant. The more realistic age-dependent functions provide a better fit to the experimental data compared to a constant mortality function and, hence, influence the predicted impact of ITNs on malaria transmission potential. In models with age-independent mortality, there is a great reduction for the vectorial capacity under exposure compared to no exposure. However, the two age-dependent functions predicted an even larger reduction due to exposure, highlighting the impact of incorporating age in the mortality rates. These results further show that multiple exposures to ITNs had a considerable effect on the vectorial capacity. Overall, the study highlights the importance of including age dependency in mathematical models of vector-borne disease transmission and in fully understanding the impact of interventions.

Interventions against malaria are most commonly targeted on the adult mosquitoes, which transmit the infection from person to person. One of the most important interventions are bed-nets, treated with insecticides. Unfortunately, extensive exposure of mosquitoes to insecticide has led to widespread evolution of insecticide resistance, which might threaten control strategies. Piecing together the overall impact of resistance on the efficacy of insecticide-treated nets is complex, but can be informed by the use of mathematical models. However, there are some assumptions that the models frequently use which are not realistic in terms of the mosquito biology. In this paper, we formulate a model that includes age-dependent mortality rates, an important parameter in vector control since control strategies most commonly aim to reduce the lifespan of the mosquitoes. By using novel data collected using field-derived insecticide-resistant mosquitoes, we explore the effects that the presence of insecticides on nets have on the mortality rates, as well as the difference incorporating age dependency in the model has on the results. We find that including age-dependent mortality greatly alters the anticipated effects of insecticide-treated nets on mosquito transmission potential, and that ignoring this realism potentially overestimates the negative impact of insecticide resistance.

Vector-Borne Disease Models with Active and Inactive Vectors: A Simple Way to Consider Biting Behavior

Vector-borne diseases are a serious public health problem, mosquitoes being one of the most important vectors. To analyze the dynamics of this type of disease, Ross-Macdonald models are commonly used. In its simplest formulation and the most common in scientific literature, it is assumed that all mosquitoes are biting at a given rate. To improve this general assumption, we developed a vector-borne disease model with active and inactive vectors as a simple way to incorporate the more general characteristics of mosquito feeding behavior into disease dynamics. Our objective is to obtain an estimate of the Ross-Macdonald biting rate from the feeding parameters that reproduce the same dynamics as the model with active and inactive vectors. Two different cases were analyzed: a SIS-SI model and a SIR-SI model with a single epidemic. Different methods to estimate the biting rate in the Ross-Macdonald model were proposed and analyzed. To compare the results of the models, different epidemiological indicators were considered. When the biting rate is estimated considering that both models have the same basic reproduction number, very similar disease dynamics are obtained. This method is a simple way to incorporate the mosquito feeding behavior into the standard Ross-Macdonald model.

Controlling Wolbachia Transmission and Invasion Dynamics among Aedes Aegypti Population via Impulsive Control Strategy

This work is devoted to analyzing an impulsive control synthesis to maintain the self-sustainability of Wolbachia among Aedes Aegypti mosquitoes. The present paper provides a fractional order Wolbachia invasive model. Through fixed point theory, this work derives the existence and uniqueness results for the proposed model. Also, we performed a global Mittag-Leffler stability analysis via Linear Matrix Inequality theory and Lyapunov theory. As a result of this controller synthesis, the sustainability of Wolbachia is preserved and non-Wolbachia mosquitoes are eradicated. Finally, a numerical simulation is established for the published data to analyze the nature of the proposed Wolbachia invasive model.

AN IMMUNO-EPIDEMIOLOGICAL VECTOR–HOST MODEL WITH WITHIN-VECTOR VIRAL KINETICS

A current challenge for disease modeling and public health is understanding pathogen dynamics across scales since their ecology and evolution ultimately operate on several coupled scales. This is particularly true for vector-borne diseases, where within-vector, within-host, and between vector–host populations all play crucial roles in diversity and distribution of the pathogen. Despite recent modeling efforts to determine the effect of within-host virus-immune response dynamics on between-host transmission, the role of within-vector viral dynamics on disease spread is overlooked. Here, we formulate an age-since-infection-structured epidemic model coupled to nonlinear ordinary differential equations describing within-host immune-virus dynamics and within-vector viral kinetics, with feedbacks across these scales. We first define the within-host viral-immune response and within-vector viral kinetics-dependent basic reproduction number [Formula: see text] Then we prove that whenever [Formula: see text] the disease-free equilibrium is locally asymptotically stable, and under certain biologically interpretable conditions, globally asymptotically stable. Otherwise, if [Formula: see text] it is unstable and the system has a unique positive endemic equilibrium. In the special case of constant vector to host inoculum size, we show the positive equilibrium is locally asymptotically stable and the disease is weakly uniformly persistent. Furthermore, numerical results suggest that within-vector-viral kinetics and dynamic inoculum size may play a substantial role in epidemics. Finally, we address how the model can be utilized to better predict the success of control strategies such as vaccination and drug treatment.

Age-structured vectorial capacity reveals timing, not magnitude of within-mosquito dynamics is critical for arbovirus fitness assessment

Background

Transmission dynamics of arboviruses like Zika virus are often evaluated by vector competence (the proportion of infectious vectors given exposure) and the extrinsic incubation period (EIP, the time it takes for a vector to become infectious), but vector age is another critical driver of transmission dynamics. Vectorial capacity (VC) is a measure of transmission potential of a vector-pathogen system, but how these three components, EIP, vector competence and vector age, affect VC in concert still needs study.

Methods

The interaction of vector competence, EIP, and mosquito age at the time of infection acquisition (Age_acquisition) was experimentally measured in an Aedes aegypti-ZIKV model system, as well as the age-dependence of probability of survival and the willingness to bite. An age-structured vectorial capacity framework (VC_age) was then developed using both EIP_Min and EIP_Max, defined as the time to first observed minimum proportion of transmitting mosquitoes and the time to observed maximum proportion of transmitting mosquitoes.

Results

The within-mosquito dynamics of vector competence/EIP were not significant among treatments where mosquitoes were exposed at different ages. However, VC_age revealed: (i) age-dependence in vector-virus interactions is important for transmission success; (ii) lower vector competence but at shorter EIPs was sufficient for transmission perpetuation; and (iii) R₀ may be overestimated by using non-age-structured VC.

Conclusions

The results indicate that ultimately the temporal component of the virus-vector dynamics is most critical, especially when exposure occurred at advanced mosquito age. While our study is limited to a single virus-vector system, and a multitude of other factors affect both vector competence and mosquito mortality, our methods can be extrapolated to these other scenarios. Results indicate that how ‘highly’ or ‘negligibly’ competent vectors are categorized may need adjustment.

Robustness of the reproductive number estimates in vector-borne disease systems

BACKGROUND

The required efforts, feasibility and predicted success of an intervention strategy against an infectious disease are partially determined by its basic reproduction number, R0. In its simplest form R0 can be understood as the product of the infectious period, the number of infectious contacts and the per-contact transmission probability, which in the case of vector-transmitted diseases necessarily extend to the vector stages. As vectors do not usually recover from infection, they remain infectious for life, which places high significance on the vector's life expectancy. Current methods for estimating the R0 for a vector-borne disease are mostly derived from compartmental modelling frameworks assuming constant vector mortality rates. We hypothesised that some of the assumptions underlying these models can lead to unrealistic high vector life expectancies with important repercussions for R0 estimates.

METHODOLOGY AND PRINCIPAL FINDINGS

Here we used a stochastic, individual-based model which allowed us to directly measure the number of secondary infections arising from one index case under different assumptions about vector mortality. Our results confirm that formulas based on age-independent mortality rates can overestimate R0 by nearly 100% compared to our own estimate derived from first principles. We further provide a correction factor that can be used with a standard R0 formula and adjusts for the discrepancies due to erroneous vector age distributions.

CONCLUSION

Vector mortality rates play a crucial role for the success and general epidemiology of vector-transmitted diseases. Many modelling efforts intrinsically assume these to be age-independent, which, as clearly demonstrated here, can lead to severe over-estimation of the disease's reproduction number. Our results thus re-emphasise the importance of obtaining field-relevant and species-dependent vector mortality rates, which in turn would facilitate more realistic intervention impact predictions.

The Interaction between Vector Life History and Short Vector Life in Vector-Borne Disease Transmission and Control

Epidemiological modelling has a vital role to play in policy planning and prediction for the control of vectors, and hence the subsequent control of vector-borne diseases. To decide between competing policies requires models that can generate accurate predictions, which in turn requires accurate knowledge of vector natural histories. Here we highlight the importance of the distribution of times between life-history events, using short-lived midge species as an example. In particular we focus on the distribution of the extrinsic incubation period (EIP) which determines the time between infection and becoming infectious, and the distribution of the length of the gonotrophic cycle which determines the time between successful bites. We show how different assumptions for these periods can radically change the basic reproductive ratio (R₀) of an infection and additionally the impact of vector control on the infection. These findings highlight the need for detailed entomological data, based on laboratory experiments and field data, to correctly construct the next-generation of policy-informing models.

The basic reproductive ratio (R₀) is a crucial measure of transmission intensity, lying at the interface between mathematical modelling and policy decision making. If control measures can induce a situation where R₀ ≤ 1 for a sustained period of time then the pathogen must be eradicated. For diseases spread by short-lived insect vectors a modeller can not calculate R₀ without addressing questions of chance such as, “What percentage of vectors will survive their extrinsic incubation period (EIP) to become infectious?”. Classical Ross-Macdonald theory provides answers for the modeller by making certain concrete assumptions, such as a fixed length EIP and exponentially distributed times between vector blood-meals. Using bluetongue virus spread by biting midges as an exemplar we demonstrate that biologically plausible alterations to the classical assumptions can significantly change the modeller’s prediction of R₀ with both serious over-estimation and under-estimation being possible. The important modelling/control question, “How does R₀ respond to increased vector per-capita mortality?”, is also found to depend strongly on details of the vector life-cycle expressed in the language of probability distributions.

Dengue and chikungunya: modelling the expansion of mosquito-borne viruses into naïve populations

With the recent global spread of a number of mosquito-borne viruses, there is an urgent need to understand the factors that contribute to the ability of viruses to expand into naïve populations. Using dengue and chikungunya viruses as case studies, we detail the necessary components of the expansion process: presence of the mosquito vector; introduction of the virus; and suitable conditions for local transmission. For each component we review the existing modelling approaches that have been used to understand recent emergence events or to assess the risk of future expansions. We identify gaps in our knowledge that are related to each of the distinct aspects of the human-mosquito transmission cycle: mosquito ecology; human-mosquito contact; mosquito-virus interactions; and human-virus interactions. Bridging these gaps poses challenges to both modellers and empiricists, but only through further integration of models and data will we improve our ability to better understand, and ultimately control, several infectious diseases that exert a significant burden on human health.

Mosquito age and avian malaria infection

Background

The immune system of many insects wanes dramatically with age, leading to the general prediction that older insects should be more susceptible to infection than their younger counterparts. This prediction is however challenged by numerous studies showing that older insects are more resistant to a range of pathogens. The effect of age on susceptibility to infections is particularly relevant for mosquitoes given their role as vectors of malaria and other diseases. Despite this, the effect of mosquito age on Plasmodium susceptibility has been rarely explored, either experimentally or theoretically.

Methods

Experiments were carried out using the avian malaria parasite Plasmodium relictum and its natural vector in the field, the mosquito Culex pipiens. Both innate immune responses (number and type of circulating haemocytes) and Plasmodium susceptibility (prevalence and burden) were quantified in seven- and 17-day old females. Whether immunity or Plasmodium susceptibility are modulated by the previous blood feeding history of the mosquito was also investigated. To ensure repeatability, two different experimental blocks were carried out several weeks apart.

Results

Haemocyte numbers decrease drastically as the mosquitoes age. Despite this, older mosquitoes are significantly more resistant to a Plasmodium infection than their younger counterparts. Crucially, however, the age effect is entirely reversed when old mosquitoes have taken one previous non-infected blood meal.

Conclusions

The results agree with previous studies showing that older insects are often more resistant to infections than younger ones. These results suggest that structural and functional alterations in mosquito physiology with age may be more important than immunity in determining the probability of a Plasmodium infection in old mosquitoes. Possible explanations for why the effect is reversed in blood-fed mosquitoes are discussed. The reversal of the age effect in blood fed mosquitoes implies that age is unlikely to have a significant impact on mosquito susceptibility in the field.

Electronic supplementary material

The online version of this article (doi:10.1186/s12936-015-0912-z) contains supplementary material, which is available to authorized users.