A computational lens for sexual-stage transmission, reproduction, fitness and kinetics in Plasmodium falciparum
© The Author(s) 2016
Received: 12 August 2016
Accepted: 14 September 2016
Published: 21 September 2016
The burden of falciparum malaria remains unacceptably high in much of sub-Saharan Africa and massive efforts are underway to eliminate the parasite. While symptoms of malaria are caused by asexual reproduction of the parasite, transmission to new human hosts relies entirely on male and female sexual-stage parasites, known as gametocytes. Successful transmission can be observed at very low gametocyte densities, which raises the question of whether transmission-enhancing mechanisms exist in the human host, the mosquito, or both.
A new computational model was developed to investigate the probability of fertilization over a range of overdispersion parameters and male gamete exploration rates. Simulations were used to fit a likelihood surface for data on rates of mosquito infection across a wide range of host gametocyte densities.
The best fit simultaneously requires very strong overdispersion and faster gamete exploration than is possible with random swimming in order to explain typical prevalence levels in mosquitoes. Gametocyte overdispersion or clustering in the human host and faster gamete exploration of the mosquito blood meal are highly probably given these results.
Density-dependent gametocyte clustering in the human host, and non-random searching (e.g., chemotaxis) in the mosquito are probable. Future work should aim to discover these mechanisms, as disrupting parasite development in the mosquito will play a critical role in eliminating malaria.
In order to investigate overdispersion, gamete exploration of the blood meal, and other mechanisms that could influence parasite fertilization rates in the blood meal, a new computational model was developed as described in the “Methods” section. The resulting model generates a distribution for the number of zygotes formed per feed under given conditions and parameters, as well as the fraction of feeds with at least one zygote. Model outputs to experimental data for oocyst counts and prevalence (i.e., fraction of mosquitoes infected) were compared over a wide range of gametocyte densities [10, 11]. It is challenging to fit these data and others [12–14] because models based on random Poisson draws tend to underpredict success at low gametocyte densities and either rise too rapidly towards complete success or plateau artificially at a lower success rate. The new model structure allows for the inclusion and exploration of transmission-enhancing mechanisms such as aggregation and bloodmeal exploration.
The key model parameters are gametocyte density g in the human host, female sex ratio f, viable male gametes per male gametocyte m, blood meal size b, overdispersion parameter k, and blood meal coverage per male gamete c.
For Poisson gametocyte draws, a Poisson draw for female gametocytes X is made with mean bfg and a Poisson draw for male gametocytes Y with mean b(1 − f)g. For negative binomial draws, a single draw for number of gametocytes is made with mean bg and overdispersion parameter k. Each gametocyte then has probability f of being female and a series of random draws results in the number of female and male gametocytes in that blood meal, represented by X and Y, respectively.
Within the blood meal, the number of male gametocytes is multiplied by m (range 0–8) to get the number of male gametes in the blood meal. The female gametocytes are assumed to be uniformly distributed throughout the volume, although if female and male gametocytes cluster in the human host and then retain proximity in the blood meal, that would facilitate a higher effective blood meal coverage per male gamete.
The blood meal coverage c can be estimated through various means. The swimming speed of exflagellated male gametes in ookinete culture medium has been measured as 5 µm/s, which is only sufficient to map out 0.0001 of a 2 µL blood meal in 30 min . This is due to the motion-orthogonal oscillation exploring a 5-µm diameter, which carves out an explored cylinder at 5 µm/s. Thirty minutes of independent exploration at this rate corresponds to approximately 0.0001 of a 2 µL blood meal. Wilson et al. also observed that the wave speed down the male gamete is a factor of ten faster (50 µm/s) with a wavelength similar in dimension to a red blood cell (RBC), so if the male gamete could move at wave speed in the RBC-dense environment of a blood meal, it could potentially explore up to 0.001 of the blood meal in 30 min. Coverage values above 0.001 would tend to mean that effective coverage rates were occurring faster than wave swimming speed, which would invoke alternative mechanisms such as chemotaxis or conserved spatial proximity.
For the varying-k model, k = 0.5 for gametocyte densities below 1/µL, and then k = 0.5 + 0.5 log 10 (g) above 1/µL.
The probability of mosquito infection at low gametocytaemias has been fit previously using experimental infections , and the equations describing the fits were used here, with a standard deviation of 0.01. The probability of successful infection and number of oocysts at higher gametocytaemias using experimental infections  were manually digitized here from Fig. 1 in the manuscript, with standard deviations taken from the error bars for probability of infection. These previously published data show a steady increase in infection success with gametocyte density, rising above 80 % by 80 gametocytes/µL, with an infection rate of almost 8 % down at the microscopic threshold of 5 gametocytes/µL.
For each set of parameters, 10,000 mosquito feeds were simulated for each of various values of gametocyte densities—0.1, 0.2, 0.5, 1, 2, 5, 10, 20, 50, and 100 gametocytes/µL. The fit quality was Σ −(simulation mean-data mean) 2 /(data std dev) 2 for each density, with Schneider  data used for gametocyte densities of 4 or less and Da  data used for gametocyte densities of 4 or more.
Likelihood of zygote formation under different overdispersion scenarios and male gamete swimming speeds
Figure 2 (black) shows the results for a Poisson draw of male and female gametocytes, with male gametes exploring the blood meal at measured free solution swimming speeds, as estimated in . Infection success is negligible below 10 gametocytes/µL, after which fertilization increases rapidly to over 80 % success by 100 gametocytes/µL. This baseline model fails to capture the non-trivial infection rates achieved by sub-microscopic gametocyte densities . Note that since the model outputs zygote counts, these should necessarily be higher than experimental oocyst counts. Previous work counting gametocytes inside mosquito guts estimated overdispersion in blood meals ranging from k < 1 to k = 3, with k = 3 the best overall value  and from k = 0.5 at low gametocyte densities to approximately k = 3 at high gametocyte densities . Such overdispersion values are far from Poisson-distributed. When overdispersion k = 1 is used in the model (Fig. 2, blue), success relative to baseline increases at low gametocyte densities and decreases above 50 gametocytes/µL, but the increase at low densities remains below observations . This strong overdispersion of gametocyte draws, such as would come from within-host clustering of gametocytes, reshapes the infection success curve to be much more like experimental data, but computed zygote counts remain well below experimental oocyst counts, instead of being substantially higher. This is problematic given that 100+ fold  to 1000+ fold  reduction in parasite numbers are estimated in the transition from gametocyte to oocyst.
If male gametes travel at their oscillatory wave speed as described in “Methods” section, blood meal coverage per gamete is multiplied by ten. For a Poisson draw of male and female gametocytes and this faster exploration, fertilization rates rise uniformly over baseline (Fig. 2, green). However, infection success remains far too low at sub-microscopic gametocyte densities, infection rates rise too rapidly above 50 gametocytes/µL, and below 50 gametocytes/µL zygote densities are below oocyst data.
The formation of zygotes requires more than one explanatory factor
The likelihood surface in Fig. 3 shows that both fast exploration (coverage > 0.001) and strong overdispersion (k < 1) are jointly favoured, with best values approximately k = 0.6 and coverage = 0.004. Reduced overdispersion, regardless of blood meal coverage, results in lower likelihoods, as do all physiologically possible random swimming speeds independent of clustering. This strongly suggests that at least two separate mechanisms are responsible for observed transmission success at low gametocyte densities.
Finally, letting k vary as described in “Methods” section and previously observed  results in the best likelihood value observed. The maximum likelihood for this function of k corresponds to a gamete coverage of 0.003, and the resulting zygote distributions and success rates are seen in Fig. 4 (orange). The fit is improved, and mean zygote counts remain about an order of magnitude above the experimental oocyst data, but the distribution at high gametocyte densities does not rise to such extreme outlier values as the constant-k best fit. This makes sense as an optimal strategy for the parasite, as strong overdispersion at low gametocyte densities increases the fraction of mosquitoes infected, while less overdispersion at high densities increases transmission success while reducing the parasite load in the most infected mosquitoes.
Given that gametocyte sex ratio varies and male gametocytes can produce 0-8 viable gametes, sensitivity was explored, while holding other parameters at best-fit values from Fig. 3. The results are seen in Additional file 1: Figure S4: rather than a single likelihood peak, there is a ridge of likely parameter values from moderate female sex ratios of 0.5–0.75 and two male gametes, to higher female sex ratios with higher numbers of male gametes. Given that there is known female bias in gametocyte counts and expected numbers of male gametes tend to be two to three, this likelihood ridge matches observations. Finally, male gametes versus coverage for the varying-k optimal fit were varied in Additional file 1: Figure S5. There is a likelihood ridge conserving blood meal exploration: two male gametes with coverage = 0.003 is approximately as likely as six to seven male gametes and coverage = 0.001.
Presence of a likelihood peak in both dimensions of Fig. 3 was robust to sex ratio, male gametes per gametocyte, and the likelihood function, although the optimal values of k could vary from 0.6 to 1.2 and the optimal values of c could vary from 0.001 to 0.006 depending on the assumptions. In each set of conditions, the conclusion of both mechanisms being jointly favoured is conserved.
Aggregation has often been cited as a likely phenomenon to explain infection rates of mosquitoes [4, 18] yet it still has not been adequately explained from a mechanistic standpoint. Beyond the need to discover how gametocytes carry out aggregation, the model presented here shows that aggregation on its own still remains insufficient to explain realistic probabilities of zygote formation. Additionally, the space that a male gamete can explore under the fastest realistic swimming speeds  is also shown to be insufficient. Results presented here show that in order to explain natural rates of infection at low gametocytaemia, gametocytes must use both a within-host clustering mechanism such as cytoadherence-style aggregation and at least one mechanism that increases fertilisation rates in the bloodmeal such as chemotaxis. Cytoadherence is common among different stages of malaria parasites: RBCs infected with asexual stage parasites can adhere to uninfected RBCs (rosetting) and to endothelial cells (sequestration) , egressed male gametes stick to RBCs in species-specific manners , and immature gametocytes also sequester through unknown mechanisms in the bone marrow . Although circulating mature gametocytes are more deformable than sequestered immature gametocytes , it is possible that mature gametocytes also have the capacity to adhere, either to capillary surfaces or to each other as they pass through the capillaries. The shape and/or buoyancy of gametocytes might also influence their movement through capillaries and thus their uptake by mosquitoes. Once gametes egress from the RBC, the results presented here suggest that additional mechanisms are required to reach infection levels observed in nature. Plasmodium berghei male gametes have recently been shown to have some chemotaxis towards female gametes . In P. falciparum, chemotaxis has not yet been observed but gametes in the mosquito midgut form sticky nanotubes that could facilitate sex either through cytoadherence or chemical gradients . It is also possible that spatial proximity of gametocytes in the host circulation could be maintained after entering the mosquito, or that the process of feeding and diuresis creates dynamics in blood movement inside the mosquito’s gut that brings gametocytes closer to each other inside the mosquito’s gut. Flow dynamics have been observed for broadcast spawners and can greatly enhance fertilization rates in seemingly improbable situations .
Gametocytes are responsible for transmission of P. falciparum from humans to mosquitoes, and they can successfully infect a mosquito even at sub-microscopic densities in the human host . Existing models have difficulty capturing the functional form of the rate of mosquito infection versus gametocyte density, and mechanistic models with non-clustered draws of gametocytes and random searching of the blood meal cannot reach the infection levels observed in existing data. In order to explore the existence of possible transmission-enhancing mechanisms, a new mechanistic model of sexual-stage transmission was developed and explored here. Both high degrees of clustering in the human host and faster than the fastest possible random exploration in the mosquito blood meal are independently required in order to recreate existing data on mosquito infection. These results strongly indicate that transmission-enhancing mechanisms are likely to exist, while providing estimates of their nature and magnitude.
The present work examines the probability of successful gamete fertilization in the mosquito midgut for varying gametocyte densities at different fixed values of sex ratio, male gametes per gametocyte, and bloodmeal volume. The average infectiousness over the course of a single infection would be influenced by time-varying parasite densities  as well as possible variations in parasite sex ratio [26, 27]. The infectiousness over a local human population would depend on local transmission intensity and seasonality , with individual infectiousness influenced by age and prior exposure effects on immunity . Future work could investigate the net effect of these variations over the course of an infection in different transmission settings, as well as incorporating the effect of transmission-blocking immunity  to expand the modelled range of gametocyte densities.
Exploring the behaviours of these fascinating cells likely has practical relevance. Gametocytes are extremely difficult parasite stages to access due to their low numbers, yet they are absolutely critical for malaria elimination campaigns. Discovering how the sexual stages of these deadly parasites find each other, both in the human host and in the mosquito, will improve human ability to disrupt this process and thus break the cycle of transmission.
MKNL and PAE contributed to concept and design. PAE wrote the model and carried out the simulations. PAE and MKNL analysed the results and wrote the paper. Both authors read and approved the final manuscript.
We thank Arthur Talman for discussions on this topic and Edward Wenger for comments on the manuscript.
The authors declare that they have no competing interests.
Availability of data and materials
All data used in these analyses are available at the cited sources.
MKNL is supported by an MRC Career Development Award (G1100339) and by the Wellcome Trust (098051). PAE is supported by the Global Good Fund.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.
- Anderson RM, May RM. Regulation and stability of host-parasite population interactions. J Anim Ecol. 1978;47:219–47.View ArticleGoogle Scholar
- Vazquez-Prokopec GM, Perkins TA, Waller LA, Lloyd AL, Reiner RC Jr, Scott TW, et al. Coupled heterogeneities and their impact on parasite transmission and control. Trends Parasitol. 2016;32:356–67.View ArticlePubMedGoogle Scholar
- Pichon G, Awono-Ambene HP, Robert V. High heterogeneity in the number of Plasmodium falciparum gametocytes in the bloodmeal of mosquitoes fed on the same host. Parasitology. 2000;121(Pt 2):115–20.View ArticlePubMedGoogle Scholar
- Gaillard FO, Boudin C, Chau NP, Robert V, Pichon G. Togetherness among Plasmodium falciparum gametocytes: interpretation through simulation and consequences for malaria transmission. Parasitology. 2003;127:427–35.View ArticlePubMedGoogle Scholar
- Bousema T, Dinglasan RR, Morlais I, Gouagna LC, van Warmerdam T, Awono-Ambene PH, et al. Mosquito feeding assays to determine the infectiousness of naturally infected Plasmodium falciparum gametocyte carriers. PLoS ONE. 2012;7:e42821.View ArticlePubMedPubMed CentralGoogle Scholar
- Bousema T, Churcher TS, Morlais I, Dinglasan RR. Can field-based mosquito feeding assays be used for evaluating transmission-blocking interventions? Trends Parasitol. 2013;29:53–9.View ArticlePubMedGoogle Scholar
- Van Den Berghe L, Chardome M, Peel E. Superiority of preparations from skin scarification over preparations of peripheral blood for the diagnosis of malaria. An Inst Med Trop. 1952;9:553–62.Google Scholar
- Chardome M, Janssen PJ. Inquiry on malarial incidence by the dermal method in the region of Lubilash, Belgian Congo. Ann Soc Belg Med Trop. 1952;32:209–11 (in French).Google Scholar
- Lin Ouédraogo A, Gonçalves BP, Gnémé A, Wenger EA, Guelbeogo MW, Ouédraogo A, et al. Dynamics of the human infectious reservoir for malaria determined by mosquito feeding assays and ultrasensitive malaria diagnosis in Burkina Faso. J Infect Dis. 2015;213:90–9.View ArticleGoogle Scholar
- Schneider P, Bousema JT, Gouagna LC, Otieno S, van de Vegte-Bolmer M, Omar SA, et al. Submicroscopic Plasmodium falciparum gametocyte densities frequently result in mosquito infection. Am J Trop Med Hyg. 2007;76:470–4.PubMedGoogle Scholar
- Da DF, Churcher TS, Yerbanga RS, Yaméogo B, Sangaré I, Ouedraogo JB, et al. Experimental study of the relationship between Plasmodium gametocyte density and infection success in mosquitoes; implications for the evaluation of malaria transmission-reducing interventions. Exp Parasitol. 2015;149:74–83.View ArticlePubMedGoogle Scholar
- 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–9.PubMedGoogle Scholar
- Eckhoff P. P. falciparum Infection durations and infectiousness are shaped by antigenic variation and innate and adaptive host immunity in a mathematical model. PLoS ONE. 2012;7:e44950.View ArticlePubMedPubMed CentralGoogle Scholar
- Churcher TS, Bousema T, Walker M, Drakeley C, Schneider P, Ouédraogo AL, et al. Predicting mosquito infection from Plasmodium falciparum gametocyte density and estimating the reservoir of infection. Elife. 2013;2:e00626.View ArticlePubMedPubMed CentralGoogle Scholar
- Wilson LG, Carter LM, Reece SE. High-speed holographic microscopy of malaria parasites reveals ambidextrous flagellar waveforms. Proc Natl Acad Sci USA. 2013;110:18769–74.View ArticlePubMedPubMed CentralGoogle Scholar
- Sinden RE, Dawes EJ, Alavi Y, Waldock J, Finney O, Mendoza J, et al. Progression of Plasmodium berghei through Anopheles stephensi is density-dependent. PLoS Pathog. 2007;3:e195.View ArticlePubMedPubMed CentralGoogle Scholar
- Vaughan JA, Noden BH, Beier JC. Population dynamics of Plasmodium falciparum sporogony in laboratory-infected Anopheles gambiae. J Parasitol. 1992;78:716–24.View ArticlePubMedGoogle Scholar
- Paul REL, Bonnet S, Boudin C, Tchuinkam T, Robert V. Aggregation in malaria parasites places limits on mosquito infection rates. Infect Genet Evol. 2007;7:577–86.View ArticlePubMedGoogle Scholar
- Rowe JA, Claessens A, Corrigan RA, Arman M. Adhesion of Plasmodium falciparum-infected erythrocytes to human cells: molecular mechanisms and therapeutic implications. Expert Rev Mol Med. 2009;11:e16.View ArticlePubMedPubMed CentralGoogle Scholar
- Templeton TJ, Keister DB, Muratova O, Procter JL, Kaslow DC. Adherence of erythrocytes during exflagellation of Plasmodium falciparum microgametes is dependent on erythrocyte surface sialic acid and glycophorins. J Exp Med. 1998;187:1599–609.View ArticlePubMedPubMed CentralGoogle Scholar
- Joice R, Nilsson SK, Montgomery J, Dankwa S, Egan E, Morahan B, et al. Plasmodium falciparum transmission stages accumulate in the human bone marrow. Sci Transl Med. 2014;6:2445.View ArticleGoogle Scholar
- Tibúrcio M, Niang M, Deplaine G, Perrot S, Bischoff E, Ndour PA, et al. A switch in infected erythrocyte deformability at the maturation and blood circulation of Plasmodium falciparum transmission stages. Blood. 2012;119:e172–80.View ArticlePubMedPubMed CentralGoogle Scholar
- Carter LM, Pollitt LC, Wilson LG, Reece SE. Ecological influences on the behaviour and fertility of malaria parasites. Malar J. 2016;15:220.View ArticlePubMedPubMed CentralGoogle Scholar
- Rupp I, Sologub L, Williamson KC, Scheuermayer M, Reininger L, Doerig C, et al. Malaria parasites form filamentous cell-to-cell connections during reproduction in the mosquito midgut. Cell Res. 2011;21:683–96.View ArticlePubMedGoogle Scholar
- Crimaldi JP. The role of structured stirring and mixing on gamete dispersal and aggregation in broadcast spawning. J Exp Biol. 2012;215:1031–9.View ArticlePubMedGoogle Scholar
- Reece SE, Drew DR, Gardner A. Sex ratio adjustment and kin discrimination in malaria parasites. Nature. 2008;453:609–14.View ArticlePubMedGoogle Scholar
- 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–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Laneri K, Paul RE, Tall A, Faye J, Diene-Sarr F, Sokhna C, et al. Dynamical malaria models reveal how immunity buffers effect of climate variability. Proc Natl Acad Sci USA. 2015;112:8786–91.View ArticlePubMedPubMed CentralGoogle Scholar