Skip to main content

Population biology of malaria within the mosquito: density-dependent processes and potential implications for transmission-blocking interventions

Abstract

Background

The combined effects of multiple density-dependent, regulatory processes may have an important impact on the growth and stability of a population. In a malaria model system, it has been shown that the progression of Plasmodium berghei through Anopheles stephensi and the survival of the mosquito both depend non-linearly on parasite density. These processes regulating the development of the malaria parasite within the mosquito may influence the success of transmission-blocking interventions (TBIs) currently under development.

Methods

An individual-based stochastic mathematical model is used to investigate the combined impact of these multiple regulatory processes and examine how TBIs, which target different parasite life-stages within the mosquito, may influence overall parasite transmission.

Results

The best parasite molecular targets will vary between different epidemiological settings. Interventions that reduce ookinete density beneath a threshold level are likely to have auxiliary benefits, as transmission would be further reduced by density-dependent processes that restrict sporogonic development at low parasite densities. TBIs which reduce parasite density but fail to clear the parasite could cause a modest increase in transmission by increasing the number of infectious bites made by a mosquito during its lifetime whilst failing to sufficiently reduce its infectivity. Interventions with a higher variance in efficacy will therefore tend to cause a greater reduction in overall transmission than a TBI with a more uniform effectiveness. Care should be taken when interpreting these results as parasite intensity values in natural parasite-vector combinations of human malaria are likely to be significantly lower than those in this model system.

Conclusions

A greater understanding of the development of the malaria parasite within the mosquito is required to fully evaluate the impact of TBIs. If parasite-induced vector mortality influenced the population dynamics of Plasmodium species infecting humans in malaria endemic regions, it would be important to quantify the variability and duration of TBI efficacy to ensure that community benefits of control measures are not overestimated.

Background

Density-dependent processes that regulate population growth are common in host-parasite systems and can influence the resilience of an infection to control interventions [1]. In a model malaria-mosquito system, the progression of Plasmodium berghei through Anopheles stephensi depends non-linearly on parasite density [2]. Two types of density dependence operate during sporogony in this system. Firstly, the transitions from the female (macro-)gametocyte to ookinete, ookinete to oocyst, and oocyst to sporozoites are all restricted at high parasite densities. These negative density-dependent processes limit sporogony at high parasite densities (when the per parasite rate of transition to the next stage tends to zero), but are relaxed as density decreases. Secondly, an additional, positive density dependence impedes the transformation from ookinete to oocyst in mosquitoes with a low number of ookinetes [2]. This mechanism (the Allee effect) initially facilitates transmission as ookinete density increases, but will make it unstable at low parasite densities.

In addition to this density-dependent sporogonic development, the mosquito's (and thus parasite's) survival could be influenced by parasite density. Indeed, there is growing evidence that the mosquitoes survival depends on its age [3, 4] and that Plasmodium may influence the life-expectancy of the mosquito [5] (although a meta-analysis of laboratory experiments investigating the effect of the malaria parasite on mosquito survival yielded inconclusive results [6]). In these experiments, the mortality rate of An. stephensi was dependent on its age and on the presence and density of P. berghei [7]. These density-dependent processes, identified in the P. berghei-An. stephensi model system, may operate in other parasite-vector combinations [8, 9], including those relevant to human malaria [1014]. Therefore, the interactions between these different positive and negative density-dependent processes may have important implications for the control of human malaria.

The intensity of malaria transmission in endemic areas is typically measured by the entomological inoculation rate (EIR); the annual number of infectious bites received by a person living in such areas. This metric makes the assumption that the infectiousness of a mosquito is the same irrespective of how many sporozoites are within the salivary glands. However, recent evidence indicates that there is a correlation between the number of P. berghei sporozoites at the bite site and the probability that the mouse host will go on to develop blood infection [15]. The importance of parasite density as a determinant of the potential for malaria transmission from the vector to the human population would depend on whether there is a correlation between the number of sporozoites within the salivary glands, the number injected into the human host, and ultimately the probability of acquiring the infection. Using quantitative PCR, Medica & Sinnis [16] recently showed that there was a statistically significant and positive correlation in the Plasmodium yoelii-An. stephensi system, despite substantial variability in the number of sporozoites injected between different mosquitoes and over time from the same mosquito. Notwithstanding such a correlation, the presence of parasite-induced, density-dependent vector mortality would make parasite density important so both sporozoite presence and density shall be presented in this paper.

Density-dependent parasite growth may impact upon various malaria control strategies. A number of promising control measures are being developed with the goal of reducing the incidence of human malaria by blocking transmission to and from the mosquito vectors. These include transmission-blocking vaccines (TBVs), which target antigens expressed on the malaria parasite within the mosquito [17], biological agents to prime the mosquitoes' immune system to reduce its infectivity [18], and the use of refractory, genetically modified mosquitoes to reduce infection of the mosquito [19]. However, most of these transmission-blocking interventions (TBIs) are, at present, only partially effective, so it is important to understand how their interaction with the non-linear processes taking place in the vector would influence overall transmission. Potential vaccines under investigation aim to target different Plasmodium life-stages within the vector [20]. Understanding how the population dynamics of the parasite within the mosquito may enhance or hinder the impact of a TBI could guide decisions as to which stage(s) to prioritize as target(s).

In this paper, mathematical models are used to investigate the potential effect on malaria transmission of the multiple density-dependent processes that have been identified in the P. berghei-An. stephensi model system. Since deterministic models can underestimate the cumulative effect of multiple non-linear functions, an individual-based stochastic model is used to capture the changes in parasite density over the different stages of sporogonic development. The highly overdispersed distribution of parasites recorded among mosquitoes [2, 16, 2123] is also modelled explicitly, as parasite aggregation will, on average, increase the influence of density-dependent regulatory processes [24]. The results of the model are used to re-fit the mortality data from mosquito experimental infections reported by Dawes et al [7] to investigate how different parasite life-stages may influence vector mortality (in [7] it was assumed that the ookinete stage was responsible for excess mosquito mortality). The full model is then used to investigate how interventions targeting different within-vector parasite life-stages would influence the transmission dynamics of Plasmodium.

Methods

A stochastic individual-based model was used to calculate the number of macrogametocytes, ookinetes, oocysts and salivary gland sporozoites at different times post-feeding in every mosquito of the hypothetical population. To facilitate comparison with published data the same post-feeding times are used as in [2], as these represent typical time-points at which mosquitoes are dissected to assess sporogonic development. Let L i j indicate the number of parasites of life-stage j within mosquito i, be it macrogametocytes ingested ( L i 1 ) at cessation of feeding, ookinetes at 15 hours ( L i 2 ), oocysts at day 10 ( L i 3 ), or sporozoites in the salivary glands at day 21 ( L i 4 ). The life-expectancy of each mosquito was then estimated from data in Dawes et al [7] and used to quantify the contribution of each mosquito to overall transmission. A graphical representation of the model and the different density-dependent functions describing sporogonic development are given in Figure 1. A full description of the model and parameter values can be found in Additional file 1.

Figure 1
figure 1

Graphical description of the mathematical model describing sporogonic development within the mosquito prior to the introduction of an intervention. The model charts the number of macrogametocytes, L i 1 , ookinetes, L i 2 , oocysts, L i 3 , and salivary gland sporozoites, L i 4 , within mosquito i. For a full list of notation see Additional file 1. The graphs under each arrow indicate the density-dependent processes describing the number of parasites which develop into the subsequent life-stage. Shaded areas indicate the 95% confidence intervals for the best-fit model fitted in [2].

Mosquito life-expectancy

It is unknown which stage(s) of Plasmodium development cause(s) the excess mortality in the mosquito. Ookinetes could increase susceptibility to bacterial infection as they perforate the mosquito midgut [25], whilst all three stages, ookinetes, oocysts, and sporozoites, could cause physiological disruption [26], nutrient depletion [26, 27] or costly immune responses [28]. Assessing the parasite density in dead mosquitoes can be difficult (unless a mosquito is dissected shortly after its death, parasite counts are unreliable), so studies investigating parasite-induced vector mortality typically explore the association between mosquito mortality and the estimated number of parasites ingested. If later sporogonic stages cause mosquito mortality, any statistical investigation must take into consideration the highly non-linear relationship between the number of parasites ingested and the number developing into the stage causing the mortality. Otherwise parasite-induced vector mortality would be harder to detect. The relative success of different TBIs will be influenced by which sporogonic stage(s) cause(s) vector mortality. Therefore, scenarios are investigated in which vector mortality is associated with ookinete load; oocyst burden at day 10, or number of salivary gland sporozoites at day 21 post-feeding by re-fitting the empirical (parabolic) hazard function of [7] to the mosquito mortality rate obtained by these authors varying the life-stage that is responsible for the excess mortality. The Akaike information criterion (AIC) was used to distinguish between model fits (the best supported model being the one with the lowest AIC) [29].

Variability in intervention efficacy

If vector mortality were determined by parasite load and not just parasite presence or absence, any variability in intervention efficacy would influence its overall effectiveness. Vaccinated populations typically generate a wide range of antibody responses [30] and the efficacy of other interventions may also vary. A partially efficacious intervention, which reduces parasite density but fails to clear infection, may increase the life-expectancy of a mosquito without substantially reducing its infectivity. This could lead to the intervention actually increasing rather than decreasing overall transmission.

By way of illustration, the variability in efficacy of all the interventions investigated in this paper was modelled on the range of responses expected from a single immunogenic TBV. Though the distribution of antibody titres generated by any TBV is not well known, it is expected to be similar to that in other human vaccines, which tend to show a log-normal distribution of antibody concentrations [30, 31]. The relationship between antibody concentration and TBV efficacy is typically highly non-linear (hyperbolic) and can be described using the Hill equation [31] (see Protocol S1). This generates a skewed distribution of intervention efficacies whose mean is always closer to 50% than the median (i.e. a negative skew at high mean efficacies and a positive skew at low mean efficacies). For comparability it was assumed that all interventions reduce the production of the specific life-stage at which they are targeted, and not additionally the production of a subsequent stage. This is analogous to TBV antibodies attacking surface proteins of a particular life-stage as soon as they appear.

Lifetime mosquito contribution to transmission

To aid interpretation of the complex concepts investigated in the paper, a measurement of parasite transmission is devised based on the contribution of an individual mosquito to onwards transmission over its lifetime (assuming that infection is acquired during its first bloodmeal). This method of estimating the lifetime transmission potential of a mosquito has been used in other mathematical approaches, though the exact definitions vary [32, 33]. In this paper, the definition will depend on whether it is assumed that the potential for malaria transmission from vectors to humans is determined by the presence or the density of salivary gland sporozoites. For the former, onwards transmission is defined as "the mean number of infectious bites made per mosquito during its lifetime". For the latter, onwards transmission is "the mean number of salivary gland sporozoites available to be injected per mosquito during its lifetime". These metrics allow transmission to be compared between mosquitoes that ingest different gametocyte densities, irrespective of other contributing factors such as the mosquito to human ratio, the human blood index, and the endemicity of malaria.

Results

Cumulative impact of density-dependent sporogonic development

The cumulative impact of the multiple density-dependent processes acting on the progression of P. berghei through the mosquito makes the relationship between the number of macrogametocytes ingested and the presence or number of salivary gland sporozoites highly non-linear (Figures 2A, D). In highly infected vertebrate hosts, reducing the number of macrogametocytes ingested by 50%, from 300,000 to 150,000, reduces sporozoite prevalence by only ~6% and sporozoite density by ~24%. The number of vertebrate hosts with very high gametocytaemia is likely to be relatively small as the majority of infected individuals tend to have intermediate or low gametocyte densities (see [34], though care should be taken when comparing gametocyte densities between species as P. falciparum typically has higher mosquito infectivity [35]). To illustrate the epidemiological importance of vertebrate hosts with relatively low gametocytaemia, all additional graphs have a logged x-axis. Figures 2B and 2E show how TBIs targeting different life-stages will influence the final presence and density of salivary gland sporozoites within the mosquito. Precisely which parasite life-stage is best to target to achieve the greatest reduction in transmission will depend on intervention efficacy, number of gametocytes ingested, and on whether malaria transmission intensity is assumed to depend on the presence or the density of salivary gland sporozoites (Figures 2C and 2F). Targeting the ookinete life-stage is the most successful strategy when the intervention is highly efficacious or the mosquito ingests a low/intermediate number of parasites (Figure 2C). This is because the Allee-effect (positive density dependence) acts on the production of oocysts at low ookinete densities (Figure 1). Reducing ookinete density per mosquito beneath a certain threshold (212 in this model system) would enhance parasite control by making the development of ookinetes to oocysts less efficient, giving additional benefits to TBIs. An intervention reducing ookinete density by 60% (blue line) in mosquitoes ingesting less than 100,000 macrogametocytes has a greater than 60% reduction in the density of salivary gland sporozoites (Figure 2E).

Figure 2
figure 2

The impact of transmission-blocking interventions which target different parasite life-stages on the prevalence and density of salivary gland sporozoites. The Figure shows the relationship between the number of macrogametocytes ingested and the mean prevalence of infectious mosquitoes (dashed lines, A and B) or the mean number of salivary gland sporozoites per mosquito (solid lines, D and E). The model (described in Protocol S1) was run with either no intervention (black line) or representing an intervention which reduced the production of a Plasmodium life-stage by 60%: macrogametocytes (yellow line); ookinetes (blue line); oocysts (green line); salivary gland sporozoites (red line). The grey dotted-dashed line in panel B and panel E (where it lies on top of the red line) indicates an overall reduction in sporozoite density/prevalence of 60% as a benchmark for comparison. The shaded areas of panels A and D depict the 95% confidence intervals for the best-fit model. Panel D illustrates the importance of using an individual-based model (with which to account for parasite aggregation), as simply combining the three density-dependent functions within a mean-based, deterministic model (thin dotted-dashed line) underestimates the severity of the non-linear relationship. Panels C and F show which is the best life-stage to target to reduce transmission for a range of gametocyte densities and intervention efficacies: macrogametocytes (yellow surface); ookinetes (blue surface); sporozoites (red surface).

Parasite life-stage causing vector mortality

The results of the above model were used to re-fit the relationship between parasite density and mosquito mortality rate using the same data and methodology presented in [7] (see Additional file 1). The model that used the number of salivary gland sporozoites as the explanatory variable gave the best fit to the experimental data (AIC = -3384), followed by oocysts (AIC = -3376) and finally ookinetes (AIC = -3370). However, this result should be interpreted with caution. Though the AIC values were able to clearly distinguish between the models that used sporozoites or ookinetes as the explanatory variable, the difference between the sporozoite and oocyst models was not so great as to exclude either model. The best-fit models were used to show the relationship between the number of gametocytes ingested and the number of infectious bites made by a mosquito during its lifetime (see Additional file 2).

The complex relationship between the aggregated parasite distribution and the range of efficacies generated by a potential TBV makes the impact of parasite-induced vector mortality far from intuitive. Ultimately, the impact of parasite-induced vector mortality will depend on the time it takes for an infected mosquito to become infectious, the shape of the relationship between parasite density and mosquito mortality and the distribution of parasites within the population (both before and after an intervention). If vector mortality is determined by the number of oocysts, a mosquito that ingests 100,000 macrogametocytes will bite, on average, 21% fewer times from 16 days after blood-feeding onwards (i.e., once it is infectious, see Additional file 1) than an uninfected mosquito (a drop from 3.8 to 3.0). This figure drops to 13% (from 3.8 to 3.3) if vector mortality is caused by the number of sporozoites within the salivary glands.

Overall transmission prior to intervention

The cumulative impact of density-dependent sporogonic development and vector mortality on overall transmission is shown in Figure 3. If the intensity of malaria transmission were determined by parasite density, then the reduction in the number of infectious bites at high parasite densities would be offset by a net increase in the number of sporozoites reaching the salivary glands. Therefore, in this scenario the highest contribution to transmission would be made by mosquitoes that ingest the greatest number of gametocytes (Figures 3E and 3G). However, if malaria transmission were determined solely by the number of infectious bites, then increases in parasite density above a certain threshold would have only a minimal impact on mosquito infectivity. Mosquitoes that feed on hosts with intermediate gametocytaemia densities could therefore make the highest contribution to overall transmission. This is seen in the absence of an intervention in Figure 3A, in which transmission peaks in mosquitoes that ingest ~100,000 macrogametocytes per bloodmeal, if vector mortality is determined by oocyst density. However, if vector mortality were determined by sporozoite density (Figure 3C), then the small increase in the prevalence of infectious mosquitoes as the number of macrogametocytes rises above 200,000 per bloodmeal is enough to counteract the reduction in the mean number of infectious bites made by the mosquito population during their lifetime. This is in part because the highly aggregated distribution of parasites within the vector population ensures that a high percentage of the parasites are within heavily infected mosquitoes, which will not survive long enough to contribute to transmission. Re-running the model assuming that parasites are Poisson-distributed within the mosquito population increases the modal relationship between gametocytaemia and transmission (see Additional file 3).

Figure 3
figure 3

The impact on transmission of interventions which target different parasite life-stages. The Figure shows the relationship between the number of macrogametocytes ingested and transmission according to whether the latter is dependent on the presence (dashed lines, panels A, B, C and D) or density (solid lines, panels E, F, G and H) of salivary gland sporozoites. The mean number of infectious bites per mosquito and the lifetime mean number of salivary gland sporozoites available to be injected per mosquito are both estimated assuming the mosquito is infected at its first bloodmeal. Parasite-induced vector mortality is dependent on the number of oocysts at day 10 (A, B, E and F) or salivary gland sporozoites on day 21 (C, D, G and H). In panels A, C, E and G the model was run with either no intervention (black line) or representing an intervention which reduced the production of a life-stage by 60%: gametocytes (yellow line); ookinetes (blue line); oocysts (green line); salivary gland sporozoites (red line). Note that the black and red lines overrun each other in panel A. The thin grey dotted-dashed line indicates an overall efficacy of 60% for comparison. The contour plots (panels B, D, F and H) show the relationship between gametocytaemia, the efficacy of an intervention which targeting oocysts and the overall reduction in transmission that this intervention causes.

Overall transmission after an intervention

A TBI with low efficacy could elevate transmission intensity above that seen prior to the intervention. For example, if parasite-induced vector mortality were determined by oocyst density in the mosquito, then according to the model, the maximal average number of infectious bites during the lifetime of a mosquito prior to an intervention would be 2.3. This could increase to 2.9 if a TBI reducing oocyst density by approximately 60% was introduced (Figure 3A). This is because mosquitoes that ingest a high number of gametocytes have a greater chance of becoming infected, but a lower life expectancy. If an intervention reduced the number of parasites of the life-stage that causes vector mortality, the life-expectancy of the mosquito would increase, boosting transmission.

Again, the impact of TBIs that target different sporogonic stages on the overall level of transmission will depend on the efficacy of the intervention, the number of gametocytes ingested by the mosquito and the biology of the infection (Figure 3). The likelihood of an intervention increasing transmission will depend on which parasite life-stage the intervention is targeting and which life-stage vector mortality is dependent on. Interventions are less likely to have a negative effect in mosquitoes which ingest a low number of gametocytes. However, if vector mortality is caused by later sporogonic stages then interventions with low efficacies which target oocysts or sporozoites may have a perverse outcome irrespective of how many gametocytes the mosquito ingests, though the size of the effect is likely to be modest (Figures 3B and 3D).

Parasite-induced vector mortality will make the effectiveness of an intervention depend on its ability to reduce parasite prevalence as well as parasite density. A partially effective intervention with a high variance in efficacy will tend to cause a greater reduction in transmission than one with a lower variance, as it will have a the greatest chance of clearing infection in some mosquitoes whilst minimising the number of bites made by those still contributing to transmission (Figure 4).

Figure 4
figure 4

The relationship between gametocytaemia and transmission for three different theoretical transmission blocking vaccines which have the same mean efficacy but differ in their variance. The two panels show the impact of the interventions which reduce the production of ookinetes by an average of 60% if malaria transmission is dependent on the presence (dashed line, panel A) or density (solid line, pane B) of salivary gland sporozoites. Predictions are shown for no intervention (black line) or for a TBV which generates a range of antibody responses in the vaccinated population corresponding to a ratio of the 97.5% percentile to the 2.5% percentile of 9 fold (low variance, light green), 68 fold (medium variance, dark green) or 6664 fold (high variance, purple) [31].

Discussion

Density-dependent processes regulating malaria parasite development within the mosquito are likely to influence the success of partially effective TBIs and should be considered when making decisions on which stage(s) to prioritize as target(s). The complex relationships between different positive and negative density-dependent processes mean that the life-stage that is best to target may also vary between settings. However, some general conclusions can be drawn from this analysis. Interventions that reduce ookinete density beneath a threshold are likely to have auxiliary benefits as they utilize positive density-dependent processes naturally restricting sporogonic development at low densities. TBIs which fail to reduce ookinete density beneath this threshold will have a lower overall impact on transmission than the efficacy achieved at the targeted life-stage. Unless the intervention successfully clears the infection, any reductions in parasite density could potentially increase the life-expectancy of the mosquito, enhancing transmission. Therefore partially effective interventions with a high variance in effectiveness will be more effective than those that have a more uniform efficacy.

Model systems provide a powerful tool for identifying testable hypotheses that may be applicable to the transmission of human malaria. The processes investigated in this paper have been identified and quantified in the P. berghei-An. stephensi laboratory model and this work needs to be repeated in naturally found parasite-vector combinations to draw definitive conclusions about human malaria. Similar density-dependent processes influencing sporogonic development within the mosquito have been identified in other parasite-vector combinations [8, 36]. Indeed, a number of studies of P. falciparum indicate a sigmoidal relationship between gametocyte density and mosquito infectivity [1014], indicating that positive and negative regulatory processes within the mosquito may influence the population dynamics of human malaria.

This paper has examined the effect of different potential TBIs on mosquitoes fed on blood with known gametocyte concentrations. More accurate estimation of an intervention's impact at the population level will require information on the distribution of gametocytes within the bitten human population, and this is likely to vary temporally [37], geographically [14], and with treatment [38]. As a result, the community impact of different TBIs will likely vary from setting to setting. Gametocyte densities of Plasmodium falciparum are lower than those in the experimental system investigated here, with a significant proportion of human hosts having densities too low to detect through standard microscopy [34]. This may improve the efficacy of TBIs if ookinete densities were low enough to restrict transmission as a consequence of the positive density-dependent processes that may operate on oocyst development.

The chance that a TBI actually increases the average number of infectious bites a human host population is exposed to will depend on the severity of parasite-induced vector mortality, the distribution of parasites within the human host and mosquito population before and after the intervention, and the shape of the relationship between gametocyte density and mosquito infectivity. The probability of infecting a feeding mosquito plateaus at moderate gametocyte densities, though the exact shape of this gametocytaemia-infectivity relationship may vary [3941]. Therefore, partially decreasing gametocyte density in human hosts with high gametocytaemia may not substantially decrease the host's ability to infect a mosquito. In such a scenario, even weak density-dependent parasite-induced vector mortality would reduce the number of infectious bites made by vectors ingesting a high number of gametocytes. The consequences of this for malaria transmission will depend on the distribution of gametocytes within the human host population. An intervention may reduce the contribution to transmission of hosts with low gametocytaemia whilst at the same time increasing the number of infectious bites made by mosquitoes (during their lifetime) that feed on highly gametocytaemic hosts. The overall change in the potential for malaria transmission will therefore depend on whether the cumulative reduction in transmission from human hosts with low number of gametocytes outweighs the increase in transmission caused by heavily infected hosts.

The model has been parameterized using mosquito mortality data from laboratory experiments where the insects will live much longer than they would do under natural conditions. In the wild, mosquitoes are likely to make relatively few bites more than 10 days after becoming infected, so malaria transmission will be highly sensitive to small changes in mosquito mortality. In addition, older mosquitoes are much more susceptible to changes in parasite density [7]. Since these are the mosquitoes contributing most to transmission (given the long extrinsic incubation period of Plasmodium within Anopheles) even a small increase in mortality of this age group may have a significant impact on transmission.

The results of this analysis indicate that interventions with an intermediate or low efficacy are more likely to cause an increase in transmission. Although it is unlikely that any TBI with a low efficacy would be deployed in a public health programme, the effectiveness of an intervention may wane over time so it is important to consider how these low efficacies will influence transmission dynamics. For example, TBV efficacy may decline relatively rapidly as the concentration of antibodies within the bloodstream falls. This means that unless vaccine efficacy is regularly boosted, antibody concentration may drop, causing the level of transmission to increase above pre-intervention levels. This is in additions to any enhancement of mosquito infectivity that may occur as the concentration of serum antibodies falls [4244].

It is interesting to note that if parasite-induced vector mortality operates in natural combinations of medical importance, then Plasmodium may benefit from eliciting a mild transmission-blocking immune response. If malaria transmission is determined by sporozoite presence (rather than density), then reducing parasite density once the mosquito is infected will tend to increase overall transmission by reducing the mortality of the insect vector. Pre-fertilization antigens are known to be a target of the natural immune response and it may be more plausible to ascribe this type of immunity to an evolutionary response by the parasite and not the human host.

Further work is required to identify whether parasite-induced vector mortality acts in malaria endemic regions so that any perverse outcome of TBIs can be thoroughly understood and contained. Taking into account the possible non-linear relationships between different sporogonic life-stages will increase the chance of identifying any changes in mosquito mortality caused by the parasite. Re-fitting vector mortality data indicates that later sporogonic stages are the most likely cause of parasite-induced vector mortality. However, it is likely that more than one parasite life-stage will play a role, as oocysts or sporozoites cannot have caused the elevated mortality seen in mosquitoes immediately after blood feeding.

Research into possible TBIs should consider how the development of the parasite within the mosquito may influence study results. The standard method of assessing the efficacy of a TBV is to compare the number of oocysts within mosquitoes fed on immunised/unimmunised blood [45, 46]. These efficacy estimates therefore already include the net effects of all density-dependent mechanisms acting upon the production of ookinetes and oocysts. Interventions targeting gametocytes and ookinetes may have had a different efficacy than that measured at the oocyst stage, depending on the number of gametocytes ingested. Understanding this may improve the accuracy of molecular methods for assessing TBV efficacy [46]. Researchers should also be aware that reductions in oocyst intensity may overestimate reductions in malaria transmission due to the actions of density-dependent sporozoite development and parasite-induced vector mortality.

In light of the possible impact of parasite-induced vector mortality, lab-based studies measuring intervention efficacy should compare the survival of mosquitoes between treatment groups over their whole lifetime. This is because vectors that ingest a lower number of parasites may have had a higher probability of surviving long enough to be dissected. Failing to control for this may cause the intervention efficacy to be overestimated. It is also important to investigate the impact of an intervention on the survival of older mosquitoes as they make the greatest contribution to overall transmission. Studies should quantify the distribution of efficacies caused by a TBI (both within the human host and the mosquito) to allow the full impact of a control measure to be estimated. It is also important to understand how the efficacy of an intervention may change over time. Phase II clinical trials of possible TBV candidates should quantify the rate and variance of antibody decay in a range of environmental settings (i.e. with different degrees of natural immunological boosting from ongoing infection). In addition, ideally the duration of both Phase II and Phase III trials should be lengthened until target antibody concentrations have returned to pre-intervention levels, allowing the full impact of starting an intervention but failing to sustain it to be quantified.

There is evidence that other regulatory processes acting upon the malaria parasite within the mosquito may also be density-dependent in natural populations. For example, Plasmodium may influence mosquito biting frequency [47], feeding persistence [48] and physiology [49]. Other factors such as parasite sex ratio may also affect the likelihood of onwards transmission [50], and therefore should be considered when quantifying the overall impact of an intervention. Similar methods should be used to understand the influence of multiple density-dependent processes on other vector-borne diseases as parasite development and vector survival have been shown to be influenced by the parasite density in both human onchocerciasis [51, 52] and lymphatic filariasis [5355].

Conclusions

This paper has highlighted the importance of understanding the population dynamics of the malaria parasite within the mosquito as processes regulating the development of the parasite may enhance or impede the effectiveness of a control strategy. Mathematical models should be used to investigate how an intervention influences both parasite development and vector mortality in order to fully evaluate the effectiveness of the new transmission blocking interventions currently under development.

Abbreviations

TBI:

(Transmission-blocking interventions)

TBV:

(Transmission-blocking vaccine).

References

  1. Dietz K: Density-dependence in parasite transmission dynamics. Parasitol Today. 1988, 4: 91-97. 10.1016/0169-4758(88)90034-8.

    Article  CAS  PubMed  Google Scholar 

  2. Sinden RE, Dawes EJ, Alavi Y, Waldock J, Finney O, Mendoza J, Butcher GA, Andrews L, Hill AV, Gilbert SC, Basáñez MG: Progression of Plasmodium berghei through Anopheles stephensi is density-dependent. PLoS Pathog. 2007, 3: e195-10.1371/journal.ppat.0030195.

    Article  PubMed Central  PubMed  Google Scholar 

  3. Harrington LC, Vermeylen F, Jones JJ, Kitthawee S, Sithiprasasna R, Edman JD, Scott TW: Age-dependent survival of the dengue vector Aedes aegypti (Diptera: Culicidae) demonstrated by simultaneous release-recapture of different age cohorts. J Med Entomol. 2008, 45: 307-313. 10.1603/0022-2585(2008)45[307:ASOTDV]2.0.CO;2.

    Article  PubMed  Google Scholar 

  4. Styer LM, Carey JR, Wang JL, Scott TW: Mosquitoes do senesce: departure from the paradigm of constant mortality. Am J Trop Med Hyg. 2007, 76: 111-117.

    PubMed Central  PubMed  Google Scholar 

  5. Anderson RA, Knols BG, Koella JC: Plasmodium falciparum sporozoites increase feeding-associated mortality of their mosquito hosts Anopheles gambiae s.l. Parasitology. 2000, 120: 329-333. 10.1017/S0031182099005570.

    Article  PubMed  Google Scholar 

  6. Ferguson HM, Read AF: Why is the effect of malaria parasites on mosquito survival still unresolved?. Trends Parasitol. 2002, 18: 256-261. 10.1016/S1471-4922(02)02281-X.

    Article  PubMed  Google Scholar 

  7. Dawes EJ, Churcher TS, Zhuang S, Sinden RE, Basáñez MG: Anopheles mortality is both age- and Plasmodium-density dependent: implications for malaria transmission. Malar J. 2009, 8: 228-10.1186/1475-2875-8-228.

    Article  PubMed Central  PubMed  Google Scholar 

  8. Poudel SS, Newman RA, Vaughan JA: Rodent Plasmodium: population dynamics of early sporogony within Anopheles stephensi mosquitoes. J Parasitol. 2008, 94: 999-1008. 10.1645/GE-1407.1.

    Article  PubMed  Google Scholar 

  9. Vaughan JA: Population dynamics of Plasmodium sporogony. Trends Parasitol. 2007, 23: 63-70. 10.1016/j.pt.2006.12.009.

    Article  PubMed  Google Scholar 

  10. Barnes KI, White NJ: Population biology and antimalarial resistance: The transmission of antimalarial drug resistance in Plasmodium falciparum. Acta Trop. 2005, 94: 230-240.

    Article  CAS  PubMed  Google Scholar 

  11. Drakeley CJ, Secka I, Correa S, Greenwood BM, Targett GA: Host haematological factors influencing the transmission of Plasmodium falciparum gametocytes to Anopheles gambiae s.s. mosquitoes. Trop Med Int Health. 1999, 4: 131-138. 10.1046/j.1365-3156.1999.00361.x.

    Article  CAS  PubMed  Google Scholar 

  12. Paul RE, Bonnet S, Boudin C, Tchuinkam T, Robert V: Aggregation in malaria parasites places limits on mosquito infection rates. Infect Genet Evol. 2007, 7: 577-586. 10.1016/j.meegid.2007.04.004.

    Article  CAS  PubMed  Google Scholar 

  13. Pichon G, Robert V, Tchuinkam T, Mulder B, Verhave JP: A quantitative analysis of the distribution of Plasmodium falciparum oocysts in Anopheles gambiae. Parasite. 1996, 3: 161-167.

    Article  Google Scholar 

  14. Stepniewska K, Price RN, Sutherland CJ, Drakeley CJ, von Seidlein L, Nosten F, White NJ: Plasmodium falciparum gametocyte dynamics in areas of different malaria endemicity. Malar J. 2008, 7: 249-10.1186/1475-2875-7-249.

    Article  PubMed Central  PubMed  Google Scholar 

  15. Kebaier C, Voza T, Vanderberg J: Kinetics of mosquito-injected Plasmodium sporozoites in mice: fewer sporozoites are injected into sporozoite-immunized mice. PLoS Pathog. 2009, 5: e1000399-10.1371/journal.ppat.1000399.

    Article  PubMed Central  PubMed  Google Scholar 

  16. Medica DL, Sinnis P: Quantitative dynamics of Plasmodium yoelii sporozoite transmission by infected anopheline mosquitoes. Infect Immun. 2005, 73: 4363-4369. 10.1128/IAI.73.7.4363-4369.2005.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  17. Saul A: Mosquito stage, transmission blocking vaccines for malaria. Curr Opin Infect Dis. 2007, 20: 476-481. 10.1097/QCO.0b013e3282a95e12.

    Article  PubMed  Google Scholar 

  18. Koella JC, Lorenz L, Bargielowski I: Microsporidians as evolution-proof agents of malaria control?. Adv Parasitol. 2009, 68: 315-327. full_text.

    Article  PubMed  Google Scholar 

  19. Christophides GK: Transgenic mosquitoes and malaria transmission. Cell Microbiol. 2005, 7: 325-333. 10.1111/j.1462-5822.2005.00495.x.

    Article  CAS  PubMed  Google Scholar 

  20. Carter R, Mendis KN, Miller LH, Molineaux L, Saul A: Malaria transmission-blocking vaccines--how can their development be supported?. Nat Med. 2000, 6: 241-244. 10.1038/73062.

    Article  CAS  PubMed  Google Scholar 

  21. Beier JC, Onyango FK, Ramadhan M, Koros JK, Asiago CM, Wirtz RA, Koech DK, Roberts CR: Quantitation of malaria sporozoites in the salivary glands of wild Afrotropical Anopheles. Med Vet Entomol. 1991, 5: 63-70. 10.1111/j.1365-2915.1991.tb00522.x.

    Article  CAS  PubMed  Google Scholar 

  22. Medley GF, Sinden RE, Fleck S, Billingsley PF, Tirawanchai N, Rodriguez MH: Heterogeneity in patterns of malarial oocyst infections in the mosquito vector. Parasitology. 1993, 106: 441-449. 10.1017/S0031182000076721.

    Article  PubMed  Google Scholar 

  23. Vaughan JA, Noden BH, Beier JC: Sporogonic development of cultured Plasmodium falciparum in six species of laboratory-reared Anopheles mosquitoes. Am J Trop Med Hyg. 1994, 51: 233-243.

    CAS  PubMed  Google Scholar 

  24. Churcher TS, Ferguson NM, Basáñez MG: Density dependence and overdispersion in the transmission of helminth parasites. Parasitology. 2005, 131: 121-132. 10.1017/S0031182005007341.

    Article  CAS  PubMed  Google Scholar 

  25. Jahan N, Docherty PT, Billingsley PF, Hurd H: Blood digestion in the mosquito, Anopheles stephensi: the effects of Plasmodium yoelii nigeriensis on midgut enzyme activities. Parasitology. 1999, 119: 535-541. 10.1017/S0031182099005090.

    Article  CAS  PubMed  Google Scholar 

  26. Beier JC: Malaria parasite development in mosquitoes. Annu Rev Entomol. 1998, 43: 519-543. 10.1146/annurev.ento.43.1.519.

    Article  CAS  PubMed  Google Scholar 

  27. Hurd H, Hogg JC, Renshaw M: Interactions between bloodfeeding, fecundity and infection in mosquitos. Parasitology Today. 1995, 11: 411-416. 10.1016/0169-4758(95)80021-2.

    Article  Google Scholar 

  28. Dimopoulos G, Christophides GK, Meister S, Schultz J, White KP, Barillas-Mury C, Kafatos FC: Genome expression analysis of Anopheles gambiae: responses to injury, bacterial challenge, and malaria infection. Proc Natl Acad Sci USA. 2002, 99: 8814-8819. 10.1073/pnas.092274999.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  29. Burnham K, Anderson D: Model Selection and Multi-Model Inference: A Practical Information-Theoretic Approach. 2004, New York: Springer, 2

    Chapter  Google Scholar 

  30. Saul A, Fay MP: Human immunity and the design of multi-component, single target vaccines. PLoS One. 2007, 2: e850-10.1371/journal.pone.0000850.

    Article  PubMed Central  PubMed  Google Scholar 

  31. Saul A: Efficacy model for mosquito stage transmission blocking vaccines for malaria. Parasitology. 2008, 135: 1497-1506. 10.1017/S0031182008000280.

    Article  CAS  PubMed  Google Scholar 

  32. Saul AJ, Graves PM, Kay BH: A cyclical feeding model for pathogen transmission and its application to determine vectorial capacity from vector infection-rates. J Appl Ecol. 1990, 27: 123-133. 10.2307/2403572.

    Article  Google Scholar 

  33. Smith DL, McKenzie FE: Statics and dynamics of malaria infection in Anopheles mosquitoes. Malar J. 2004, 3: 13-10.1186/1475-2875-3-13.

    Article  PubMed Central  PubMed  Google Scholar 

  34. Okell LC, Ghani AC, Lyons E, Drakeley CJ: Submicroscopic infection in Plasmodium falciparum-endemic populations: a systematic review and meta-analysis. J Infect Dis. 2009, 200: 1509-1517. 10.1086/644781.

    Article  PubMed  Google Scholar 

  35. Schneider P, Bousema JT, Gouagna LC, Otieno S, van de Vegte-Bolmer M, Omar SA, Sauerwein RW: Submicroscopic Plasmodium falciparum gametocyte densities frequently result in mosquito infection. Am J Trop Med Hyg. 2007, 76: 470-474.

    PubMed  Google Scholar 

  36. Rosenberg R, Koontz LC, Carter R: Infection of Aedes aegypti with zygotes of Plasmodium gallinaceum fertilized in vitro. J Parasitol. 1982, 68: 653-656. 10.2307/3280924.

    Article  CAS  PubMed  Google Scholar 

  37. Drakeley CJ, Carneiro I, Reyburn H, Malima R, Lusingu JP, Cox J, Theander TG, Nkya WM, Lemnge MM, Riley EM: Altitude-dependent and -independent variations in Plasmodium falciparum prevalence in northeastern Tanzania. J Infect Dis. 2005, 191: 1589-1598. 10.1086/429669.

    Article  PubMed  Google Scholar 

  38. Adjuik M, Babiker A, Garner P, Olliaro P, Taylor W, White N: Artesunate combinations for treatment of malaria: meta-analysis. Lancet. 2004, 363: 9-17. 10.1016/S0140-6736(03)15162-8.

    Article  CAS  PubMed  Google Scholar 

  39. Draper CC: Observations on the infectiousness of gametocytes in hyperendemic malaria. Trans R Soc Trop Med Hyg. 1953, 47: 160-165. 10.1016/0035-9203(53)90072-8.

    Article  CAS  PubMed  Google Scholar 

  40. Jeffery GM, Eyles DE: Infectivity to mosquitoes of Plasmodium falciparum as related to gametocyte density and duration of infection. Am J Trop Med Hyg. 1955, 4: 781-789.

    CAS  PubMed  Google Scholar 

  41. Tchuinkam T, Mulder B, Dechering K, Stoffels H, Verhave JP, Cot M, Carnevale P, Meuwissen JH, Robert V: Experimental infections of Anopheles gambiae with Plasmodium falciparum of naturally infected gametocyte carriers in Cameroon: factors influencing the infectivity to mosquitoes. Trop Med Parasitol. 1993, 44: 271-276.

    CAS  PubMed  Google Scholar 

  42. Naotunne TD, Rathnayake KD, Jayasinghe A, Carter R, Mendis KN: Plasmodium cynomolgi: serum-mediated blocking and enhancement of infectivity to mosquitoes during infections in the natural host, Macaca sinica. Exp Parasitol. 1990, 71: 305-313. 10.1016/0014-4894(90)90035-B.

    Article  CAS  PubMed  Google Scholar 

  43. Peiris JS, Premawansa S, Ranawaka MB, Udagama PV, Munasinghe YD, Nanayakkara MV, Gamage CP, Carter R, David PH, Mendis KN: Monoclonal and polyclonal antibodies both block and enhance transmission of human Plasmodium vivax malaria. Am J Trop Med Hyg. 1988, 39: 26-32.

    CAS  PubMed  Google Scholar 

  44. Tirawanchai N, Winger LA, Nicholas J, Sinden RE: Analysis of immunity induced by the affinity-purified 21-kilodalton zygote-ookinete surface antigen of Plasmodium berghei. Infect Immun. 1991, 59: 36-44.

    PubMed Central  CAS  PubMed  Google Scholar 

  45. Chowdhury DR, Angov E, Kariuki T, Kumar N: A potent malaria transmission blocking vaccine based on codon harmonized full length Pfs48/45 expressed in Escherichia coli. PLoS One. 2009, 4: e6352-10.1371/journal.pone.0006352.

    Article  PubMed Central  PubMed  Google Scholar 

  46. Miura K, Keister DB, Muratova OV, Sattabongkot J, Long CA, Saul A: Transmission-blocking activity induced by malaria vaccine candidates Pfs25/Pvs25 is a direct and predictable function of antibody titer. Malar J. 2007, 6: 107-10.1186/1475-2875-6-107.

    Article  PubMed Central  PubMed  Google Scholar 

  47. Koella JC, Sorensen FL, Anderson RA: The malaria parasite, Plasmodium falciparum, increases the frequency of multiple feeding of its mosquito vector, Anopheles gambiae. Proc Biol Sci. 1998, 265: 763-768. 10.1098/rspb.1998.0358.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  48. Anderson RA, Koella JC, Hurd H: The effect of Plasmodium yoelii nigeriensis infection on the feeding persistence of Anopheles stephensi Liston throughout the sporogonic cycle. Proc Biol Sci. 1999, 266: 1729-1733. 10.1098/rspb.1999.0839.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  49. Hogg JC, Hurd H: Malaria-induced reduction of fecundity during the first gonotrophic cycle of Anopheles stephensi mosquitoes. Med Vet Entomol. 1995, 9: 176-180. 10.1111/j.1365-2915.1995.tb00175.x.

    Article  CAS  PubMed  Google Scholar 

  50. Mitri C, Thiery I, Bourgouin C, Paul RE: Density-dependent impact of the human malaria parasite Plasmodium falciparum gametocyte sex ratio on mosquito infection rates. Proc Biol Sci. 2009, 276: 3721-3726. 10.1098/rspb.2009.0962.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  51. Basáñez MG, Remme JHF, Alley ES, Bain O, Shelley AJ, Medley GF, Anderson RM: Density-dependent processes in the transmission of human onchocerciasis: Relationship between the numbers of microfilariae ingested and successful larval development in the simuliid vector. Parasitology. 1995, 110: 409-427. 10.1017/S0031182000064751.

    Article  PubMed  Google Scholar 

  52. Basáñez MG, Townson H, Williams JR, Frontado H, Villamizar NJ, Anderson RM: Density-dependent processes in the transmission of human onchocerciasis: Relationship between microfilarial intake and mortality of the simuliid vector. Parasitology. 1996, 113: 331-355. 10.1017/S003118200006649X.

    Article  PubMed  Google Scholar 

  53. Krishnamoorthy K, Subramanian S, Van Oortmarssen GJ, Habbema JD, Das PK: Vector survival and parasite infection: the effect of Wuchereria bancrofti on its vector Culex quinquefasciatus. Parasitology. 2004, 129: 43-50. 10.1017/S0031182004005153.

    Article  CAS  PubMed  Google Scholar 

  54. Snow LC, Bockarie MJ, Michael E: Transmission dynamics of lymphatic filariasis: vector-specific density dependence in the development of Wuchereria bancrofti infective larvae in mosquitoes. Med Vet Entomol. 2006, 20: 261-272. 10.1111/j.1365-2915.2006.00629.x.

    Article  CAS  PubMed  Google Scholar 

  55. Stolk WA, Van Oortmarssen GJ, Subramanian S, Das PK, Borsboom GJ, Habbema JD, de Vlas SJ: Assessing density dependence in the transmission of lymphatic filariasis: uptake and development of Wuchereria bancrofti microfilariae in the vector mosquitoes. Med Vet Entomol. 2004, 18: 57-60. 10.1111/j.0269-283X.2004.0470.x.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

This work has been supported by the European Commission FP7 Collaborative project TransMalariaBloc (HEALTH-F3-2008-223736). EJD also thanks the Medical Research Council UK, for a Doctoral Training Account within the Faculty of Medicine, Imperial College London.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Thomas S Churcher.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

TSC carried out the mathematical modelling and drafted the manuscript. TSC, EJD and MGB conceived the study. EJD, RES, GKC, JCK and MGB participated in the design of the study and helped to draft the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

Additional file 1: Description of the mathematical model and table of parameter values. (DOC 340 KB)

12936_2010_1393_MOESM2_ESM.DOC

Additional file 2: The impact of varying which parasite life-stage causes mosquito mortality on the mean number of potentially infectious bites made by a mosquito in its lifetime. Parasite mortality is dependent on the mean number of ingested ookinetes (blue line), oocysts at day 10 (green line) or salivary gland sporozoites on day 21 (red line). All bites made 10 days after the first bloodmeal are defined as being potentially infectious. The relationship between ookinete density and the later life-stages was estimated using the full model and used to re-fit the mosquito mortality data from [7]. (DOC 78 KB)

12936_2010_1393_MOESM3_ESM.DOC

Additional file 3: The relationship between gametocytaemia and transmission for a Poisson-distributed parasite population. The potential for malaria transmission is dependent on either salivary gland sporozoite density (solid lines, panel A) or prevalence (dashed lines, panel B). The mean number of salivary gland sporozoites available to be injected per mosquito and the mean number of infectious bites per mosquito are both estimated assuming the mosquito is infected at its first bloodmeal and correspond to the total during the mosquito's lifetime. The model was run with either no intervention (black line) or representing a TBI which reduced the production of a life-stage by 60%: gametocytes (yellow line); ookinetes (blue line); oocysts (green line); salivary gland sporozoites (red line). The thin dotted-dashed grey line (which overlaps with the red line in A) indicates an overall reduction in sporozoite density/prevalence efficacy of 60% as a benchmark for comparison. The number of parasites of each life-stage among the mosquitoes follows a Poisson distribution (i.e. k i j = 10). (DOC 61 KB)

Authors’ original submitted files for images

Rights and permissions

This article is published under license to 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.

Reprints and permissions

About this article

Cite this article

Churcher, T.S., Dawes, E.J., Sinden, R.E. et al. Population biology of malaria within the mosquito: density-dependent processes and potential implications for transmission-blocking interventions. Malar J 9, 311 (2010). https://doi.org/10.1186/1475-2875-9-311

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1475-2875-9-311

Keywords