Predicting mosquito infection from Plasmodium falciparum gametocyte density and estimating the reservoir of infection
 Cited 101
 Views 2,185
 Annotations
Abstract
Transmission reduction is a key component of global efforts to control and eliminate malaria; yet, it is unclear how the density of transmission stages (gametocytes) influences infection (proportion of mosquitoes infected). Human to mosquito transmission was assessed using 171 direct mosquito feeding assays conducted in Burkina Faso and Kenya. Plasmodium falciparum infects Anopheles gambiae efficiently at low densities (4% mosquitoes at 1/µl blood), although substantially more (>200/µl) are required to increase infection further. In a site in Burkina Faso, children harbour more gametocytes than adults though the nonlinear relationship between gametocyte density and mosquito infection means that (per person) they only contribute slightly more to transmission. This method can be used to determine the reservoir of infection in different endemic settings. Interventions reducing gametocyte density need to be highly effective in order to halt human–mosquito transmission, although their use can be optimised by targeting those contributing the most to transmission.
https://doi.org/10.7554/eLife.00626.001eLife digest
Malaria is one of the world’s most deadly infectious diseases. The most severe form is caused by the parasite Plasmodium falciparum, which can reside within red blood cells and thus evade the human immune system.
Plasmodium is transmitted between humans by mosquitoes. When a mosquito takes a blood meal from an individual infected with the parasite, the insect ingests Plasmodium gametocytes (i.e., eggs and sperm), and these go on to reproduce in the gut of the mosquito. These parasites then move to the mosquito’s salivary glands, to be injected into the next person whom the mosquito bites.
Although malaria is both preventable and curable, the mortality rates in many African countries remain high, especially among children. Reducing the transmission of malaria to mosquitoes is one of the primary goals in the global effort to control and eliminate the disease. While a range of drugs and vaccines that specifically try to reduce transmission are in development, nonmedical interventions such as mosquito nets and insecticide spraying can quickly and effectively reduce infection rates.
Here, Churcher et al. examine the dynamics of human to mosquito transmission of P. falciparum, and report that the ease with which mosquitoes become infected is not directly proportional to the density of parasite gametocytes in human blood. They found that the transmission occurs readily at very low gametocyte densities. Moreover, the transmission rate remains relatively stable as the density increases, before increasing significantly when the density reaches around 200 cells per microlitre.
Churcher et al. also challenge the assumption that children are mostly responsible for transmitting the malaria parasite by suggesting that, in certain locations, there is a more significant role for adults than previously assumed. By identifying the groups that contribute most to transmission, and targeting resources to reduce gametocyte density in those individuals, it could be possible to greatly reduce the number of infected mosquitoes and, therefore, the number of infected humans.
https://doi.org/10.7554/eLife.00626.002Introduction
The malaria parasite is transmitted among humans by anopheline mosquitoes. Male and female transmission stages (gametocytes) are ingested by the mosquito and reproduce sexually in its stomach before developing into oocysts. Once oocysts establish, it is assumed that mosquitoes will form infectious sporozoites, so the proportion of mosquitoes developing oocysts is used as a measure of mosquito infectivity. In the human malaria parasite Plasmodium falciparum, mosquito infection is thought to increase with the number of gametocytes ingested by the mosquito (Jeffery and Eyles, 1955; Graves et al., 1988; Bousema and Drakeley, 2011). However, several studies have failed to find an association (Boudin et al., 1993; Haji et al., 1996), and the precise shape of the relationship has never been rigorously quantified, particularly at very low gametocyte densities. Difficulties arise because estimates of gametocyte density have relied on microscopy, which may miss up to 80% of parasites (Dowling and Shute, 1966). More sensitive molecular methods such as Pfs25mRNA quantitative nucleic acid sequence–based amplification (QTNASBA) have been developed (Schneider et al., 2004). Unlike conventional microscopy, this technique enables gametocyte densities to be quantified over the entire epidemiologically relevant range.
Transmission reduction is now a key component of global efforts to control and eliminate malaria (Alonso et al., 2011). A wide range of novel transmissionreducing drugs and vaccines are currently under development, which aim to reduce malaria incidence by restricting human to mosquito transmission. The transmissionreducing ability of frontline therapeutics (such as artemisininbased combination therapies [ACTs], or potential combinations of ACTs with gametocytocidal drugs) is also gaining increased attention (WHO, 2012; White, 2013). It is likely that these interventions will be partially effective, but it remains unclear how much they need to reduce transmission from humans to mosquitoes in a particular location, and to which age groups they should be delivered in order to halt malaria transmission.
Mathematical models of malaria transmission tend not to include explicitly the relationship between gametocytaemia and mosquito infectivity, opting, for simplicity, to assume densityindependent transmission probabilities (Smith et al., 2007; Griffin et al., 2010), or fit functions to the relationship between the number of asexual parasites and the probability of a bloodfeeding mosquito becoming infected (Ross et al., 2006). Greater complexity, however, will be required to capture fully malaria population dynamics following the introduction of drugs and vaccines that specifically target transmission stages. For example, the World Health Organization (WHO) has recently recommended that in preelimination or elimination malaria programmes, singledose primaquine (0·25 mg base per kg) with an ACT should be given to all patients with falciparum malaria except pregnant women and infants <1 year old (WHO, 2012). More detailed mathematical models will allow the full impact of this change in strategy to be assessed in a range of different endemic settings.
The optimal use of transmissionreducing interventions will also require a better understanding of the human reservoir of infection, which is likely to vary between different endemic settings. Various studies have attempted to estimate the relative contribution of different age groups to mosquito infection using skin or membrane feeding assays (MuirheadThomson, 1957; Graves et al., 1988; Boudin et al., 1991; Githeko et al., 1992; Drakeley et al., 2000; Bonnet et al., 2003). These studies are logistically complicated and expensive, so more efficient methods of determining the reservoir of infection are required in order to target interventions optimally. A thorough understanding of the factors determining mosquito infection would allow mathematical models to predict the relative contribution of different groups to transmission from crosssectional surveys, both before and after the introduction of transmissionreducing interventions.
This article uses data from mosquito feeding assays conducted in Burkina Faso (Ouédraogo et al., 2009, Dryad: Ouédraogo et al., 2013) and Kenya (Schneider et al., 2007, Dryad: Schneider et al., 2013) to estimate the shape of the relationship between gametocyte density and mosquito infection (processed data available at Dryad, Churcher et al., 2013). Other covariates that may influence infection such as asexual parasite density (measured by microscopy) and host age were also included to improve predictions of human to mosquito transmission. These results are combined with data from a crosssectional survey conducted in a high transmission setting in Burkina Faso to predict the relative contribution of different age groups to overall malaria transmission.
Results
The relationship between the number of gametocytes in the blood and mosquito infection was found to be highly nonlinear (Figure 1A). Plasmodium falciparum infects mosquitoes at the very low gametocyte densities that predominate in natural infections and may not be detected using standard microscopy (Bousema and Drakeley, 2011). Infection rises rapidly with increasing gametocyte density, and by 1 gametocyte per microlitre, ∼4% (95% Bayesian credible interval [CI], 3–5%) of all mosquitoes develop oocysts. The subsequent increase in infection with increasing gametocyte density is best described by the Gompertz model (deviance information criterion [DIC] = 1034), which gave a significantly better fit than the linear (DIC = 1073), power (DIC = 1059), or hyperbolic (DIC = 1062) functions (Figure 1A). The bestfit model predicts that increasing density from 1 to 200 gametocytes per microlitre does not appreciably increase infection. Beyond 200 gametocytes per microlitre, infection rises again to finally plateau at ∼18% infected mosquitoes. Gametocyte density on an arithmetic scale was a better predictor of mosquito infection than gametocyte density on a logarithmic scale. Including information on the host’s asexual parasite density significantly improved model fit. Children with asexual parasite densities between 1 and 1000 parasites per microlitre were on average 27% (CI, 19–58%) less infectious to mosquitoes than those with no detectable asexual parasites, whereas those with asexual densities >1000 µl^{−1} were 77% (CI, 14–140%) more infectious. Including age improves the fit of the model suggesting that age is an important confounder, although the Bayesian credible intervals include 0. Children >6 years old were 15% (CI, 24–64%) more likely to infect mosquitoes than those of younger ages, whereas no difference in the relationship between gametocyte density and mosquito infection was detected between Burkina Faso and Kenya.
The complex shape of the relationship between gametocyte density and mosquito infection will influence the success of transmissionreducing interventions and may explain why their ability to reduce the proportion of infected mosquitoes has been shown to depend on the gametocyte density of the host (Churcher et al., 2012). For example, reducing gametocyte density by 99% in a host with 200 gametocytes per microlitre may not have much effect on their immediate contribution to transmission (Figure 1B), although it will probably reduce the duration of infectivity. The same intervention efficacy at reducing gametocyte density in a host with 300 gametocytes per microlitre would cause an appreciable reduction in human to mosquito transmission from that individual. Figure 1B can be used to estimate how a reduction in the number of gametocytes will equate to a reduction in the proportion of mosquitoes becoming infected, and hence mosquito to human transmission. How this relates to the incidence of malaria and subsequent disease will depend, in part, on the degree of immunity in the human population.
In the crosssectional survey in Burkina Faso, children harbour more gametocytes than adults, with 10 year olds having five times as many gametocytes compared to the 50 year olds (Figure 2A). However, the relationship between gametocyte density and mosquito infection established here means that adults are only 37% less likely to infect mosquitoes (infecting on average 3.5% of them, Figure 2B). Results indicate that, at the time of the survey conducted in Burkina Faso, the peak in the reservoir of infection occurs at an earlier age than the peak in gametocyte density, which will improve the (cost) effectiveness of schoolbased programmes.
Discussion
The complex shape of the relationship between gametocyte density and mosquito infection elucidated here will influence our understanding of the population dynamics of falciparum malaria and will determine the success of transmissionreducing interventions. The efficiency with which P. falciparum gametocytes can infect mosquitoes means that transmissionreducing interventions, which reduce gametocyte density, will need to be highly effective in order to reduce human–mosquito transmission. Gametocytes can infect mosquitoes at very low densities, despite the mosquito needing to ingest both male and female parasites in the same blood meal. Substantial transmission was seen from hosts with <1 gametocyte per microlitre of blood. After this initial increase in infection, considerably more gametocytes are required in order to increase further the proportion of mosquitoes developing oocysts. The causes of this are unclear, but experimental systems may be informative. Data from the rodent malaria model Plasmodium berghei (Sinden et al., 2007) suggest that the phenomenon could be associated with parasite mortality during penetration of the mosquito gut wall or due to the mosquito’s immune response. Beyond 200 gametocytes per microlitre mosquito infection rises again (although QTNASBA measurement error makes the exact gametocyte density for this transition relatively uncertain). If an intervention can keep gametocyte density beneath 200 gametocytes per microlitre in a low transmission area, it may be sufficient to push the basic reproduction number of malaria beneath one and eliminate the disease. However, the highly overdispersed distribution of gametocytes between hosts (adequately described by the negative binomial distribution) means that in the population from Burkina Faso, 30% of transmission comes from hosts with densities beneath the microscopy detection threshold of 16 gametocytes per microlitre. This means that evaluating interventions solely on their ability to reduce microscopically detectable gametocytes, as has been suggested (Graves et al., 2012), may give misleading results.
Including information on asexual parasite density (as measured by microscopy) significantly improved the fit of the model, with hosts harbouring intermediate densities having the highest infectivity. It is unclear whether asexual parasites directly influence parasite infectivity or are associated with other (here unmeasured) variables such as the multiplicity of infection (Nsango et al., 2012) or factors influencing the blood environment such as immunity (Bousema and Drakeley, 2011). Plasmodium falciparum has asynchronous waves of asexual and sexual parasites. Circulating gametocytes take 2–3 days to mature before they can infect mosquitoes (Lensen et al., 1999), which may explain the complicated relationship between infectivity and asexual density (i.e., intermediate asexual parasite densities might occur at the start of gametocyte production and have a relatively low proportion of infectious mature gametocytes). On average, older children were more infectious, supporting previous indications that transmissionreducing immunity may predominate in young children (Drakeley et al., 2006).
Figure 1A indicates that individual QTNASBA gametocyte density estimates are relatively uncertain. Despite this, they are considerably more accurate than conventional microscopy. Molecular techniques may be highly precise in a stateoftheart laboratory, but the practicalities of collecting and processing samples in nonideal field settings increases the risks of incurring measurement error. Even though uncertainty is generally accepted by those conducting the experiments, measures of precision such as confidence intervals or standard errors around such density estimates are rarely reported in the scientific literature. Methods such as those presented in this article can allow robust quantitative insight to be gained from uncertain molecular methods. Care should be taken when interpreting gametocyte density estimates as it has been suggested that the marker used in this analysis (Pfs25 mRNA) might be parasite female specific (Schneider, 2006). Even if the current QTNASBA method does only detect female gametocytes, it will neither change the qualitative conclusions of this study nor their application if other studies use the same technique. However, direct comparison of our results with those of studies using different methods for quantifying gametocyte density should be aware that the shape of the relationship may differ.
Gametocyte density measurement error is likely to cause the majority of the uncertainty seen in model outputs, although the complexity of the membrane feeding assay is also likely to contribute. The proportion of mosquitoes infected by Plasmodiuminfected blood is known to vary substantially within and between studies. This is in part due to methodological issues with the membrane feeding assay (Bousema et al., 2013) but also related to biological differences in the parasitevector combination and the blood environment (Bousema and Drakeley, 2011). Further standardising the membrane feeding assay, improving the accuracy of QTNASBA technique and the inclusion of additional covariates (such as estimates of gametocyte maturity, which have been investigated in Plasmodium vivax; Chansamut et al., 2012), would improve the accuracy of the relationship between gametocyte density and mosquito infection. It would also permit the patterns described here to be checked for consistency across time and space, and enable a wider range of functional forms to be tested. Care should also be taken when interpreting the results of feeding assays as the mosquitoes used, and their biting behaviour is likely to be different from that seen in wild mosquitoes in field situations (Bousema et al., 2013).
Transmission reduction will become increasingly important as areas approach local elimination. Identifying host age groups that contribute most to the reservoir of infection will allow the optimal targeting of malaria control. In the high transmission site in Burkina Faso, young children harbour the majority of gametocytes but are only slightly more infective to mosquitoes than adults. The per person contribution of adults estimated here is considerably greater than that predicted by other mathematical models (Ross et al., 2006) and may increase further once agedependent biting rates are taken into consideration (Carnevale et al., 1978; Ross et al., 2006). The contribution of different age groups to overall transmission will depend on local demography, agedependent protection by malaria interventions (such as the use of bed nets or time spent in the house and therefore personally protected by indoor residual spraying), and human to mosquito contact patterns. The last two of these are difficult to measure and poorly understood, reducing our ability to accurately predict the sources of infection. The reservoir of infection is likely to vary between endemicity settings and over time so multiple crosssectional surveys may be required. The relationship between gametocyte density and mosquito infection will allow estimates of the reservoir of infection from crosssectional gametocytaemia surveys without the need for logistically complicated mosquito feeding assays. The results of this article should be included within mathematical models that capture changes in gametocyte density with time (Lawpoolsri et al., 2009) to evaluate the full impact that different transmissionreducing interventions will have on transmission and prioritise the most appropriate candidates according to transmission setting.
This article provides insights into the relationship between gametocyte density and mosquito infection that are needed to predict the outcome of transmissionreducing interventions. It shows that, given the ability of very low P. falciparum gametocyte densities to establish infection in A. gambiae, transmissionreducing interventions will need to be highly efficacious at reducing gametocyte density in order to halt human to mosquito transmission. The complex and nonlinear shape of the relationship between gametocyte density and mosquito infection may explain why the ability of an interventions to reduce the proportion of infected mosquitoes has been shown to depend on the parasite load of the host (Churcher et al., 2012). Gametocyte density also changes with host age so the effectiveness of transmissionreducing interventions that target gametocytes will also vary with age. This should be considered in clinical trials of transmissionreducing candidate drugs and vaccines. The highly overdispersed distribution of gametocytes in the host population and the nonlinear relationship between gametocytes/asexual parasite density and mosquito infection means that the impact of different transmissionreducing interventions on overall transmission at the population level will be far from intuitive.
Material and methods
Estimating gametocyte density and its associated uncertainty
Request a detailed protocolQuantitative nucleic acid sequencebased amplification (QTNASBA) is routinely used to estimate pathogen density. Like all diagnostic methods, it is prone to measurement error. To understand fully the associated uncertainty, it is important to appreciate the causes of the variability and how the quantification process might magnify uncertainty of point density estimates. Here, the uncertainty is quantified by repeatedly testing samples with known gametocyte density and fitting a hierarchical mathematical model.
Nucleic acids are extracted from 50 µl of blood and then amplified in the presence of a fluorescence probe. The assay measures time to positivity (TTP), which is the time it takes for the number of target amplicons detected to exceed a defined threshold (Schneider et al., 2004). The relationship between TTP and gametocyte density is estimated by fitting a linear regression to TTP estimates generated using a sample with known gametocyte density (a 10fold dilution series of in vitro cultured gametocytes ranging from 10^{6} to 10^{1} gametocytes per millilitre). Let the observed TTP be denoted by Y then,
where β_{0} and β_{1} are regression coefficients estimates, x is the (known) parasite density from the dilution series and ε represents a normally distributed random error with mean equal to 0 and constant variance, that is ε ∼ N(0,σ^{2}). A full list of the parameters is given in Table 1. The socalled statistical calibration or inverse regression problem (Osborne, 1991) concerns the issue of making statistical inference on the value of an unknown (logtransformed) gametocyte density), denoted as ln x′, from a new TTP observation, Y′ The classical approach (Eisenhart, 1939) involves rearranging the regression model (equation 1) so that
and substituting the regression parameters with their estimates ${\widehat{\beta}}_{0}$ and ${\widehat{\beta}}_{1}$ to yield an estimate $\mathrm{ln}{\widehat{x}}^{\prime}$. Here, we adapt this approach for use in a Bayesian hierarchical model. A number of different methods have been used to do this (see Hoadley, 1970 and Hunter and Lamboy, 1981 and for discussion see Osborne, 1991), each of which has a different method for dealing with the problem of high ${\widehat{\beta}}_{1}$ values (i.e., a gentle gradient, which when used in equation 2 can generate infinitely large confidence interval estimates). Here, we use the most parsimonious approach that does not require assignment of previous distributions to (the unknown) gametocyte densities. Rather, uncertainty in the estimated regression coefficients of equation 1 is propagated numerically via equation 2 to yield uncertainty in the estimated gametocyte densities. Since all the calibration line data have a relatively steep gradient, the difference between the different methods will be relatively minor, and all will generate sensibly tight confidence interval estimates.
Calibration lines can vary between runs and batches of reagents. Here, we define an experiment as being a single plate run with the same batch of reagents each of which will have its own dilution series and samples with unknown density. The accuracy with which the individual TTP estimates fit the loglinear calibration line can be used to estimate assay measurement error (the intraassay variability, σ) by fitting a hierarchical mixedeffects model to multiple dilution series (allowing estimates of β_{0} and β_{1} to vary among experiments to determine whether this significantly improves the fit of the model). The accuracy of gametocyte density estimates can also be improved by running multiple assays on the same sample and then taking the mean of all the ln x′ estimates. The application of equation 2 yields an estimate of gametocyte density, ${\widehat{x}}^{\prime},$ on the logarithmic scale. Estimates on the arithmetic scale are generated by taking the exponent, that is $\widehat{{x}^{\prime}}=\mathrm{exp}\left(\mathrm{ln}{\widehat{x}}^{\prime}\right)$. These estimates are used in the functions below to determine the relationship between gametocyte density and the proportion of mosquitoes developing oocysts.
Mosquito infection
Request a detailed protocolPlasmodium falciparum–infected blood was collected from children and fed through a membrane to A. gambiae sensu stricto mosquitoes that were dissected 7–9 days later to determine infection (oocyst carriage). Data from 171 mosquito feeds on patients’ blood conducted in Burkina Faso (Ouédraogo et al., 2009) and Kenya (Schneider et al., 2007) were combined and used to fit a quantitative relationship between gametocyte density, asexual parasite density, host age, and mosquito infection (proportion of mosquitoes developing oocysts). A full description of these data is given in Figure 1—source data 1. Gametocyte density was quantified using QTNASBA, and the uncertainty in the density estimates (the intraassay variability) was quantified by fitting a hierarchical model to the 16 independent dilution series. In the Burkina Faso dataset, multiple assays were carried out on the same blood sample (an average of 2.63 assays per unknown blood sample). The mean gametocyte density from these multiple assays was taken to increase the accuracy of the estimates. The precise shape of the relationship was determined by fitting a range of different functional forms (a modified constant, linear, power, hyperbolic, sigmoid, and Gompertz functions) to data presented in Figure 1—source data 1 and statistically determining which gave the best fit. The proportion of mosquitoes developing oocysts, denoted by g, can be described by equation 3
Function ${f}_{i}\left({\widehat{x}}^{\prime}\right)$ describes how infection changes with increasing estimated gametocyte density with subscript i indicating the different functions tested (whose equations are given in Table 1); μ is a vector of regression coefficients and z_{1}, z_{2} and z_{3} are dummy variables denoting asexual parasite density, age of the blood donor for the membrane feeing assays, and location of provenance, respectively. Asexual parasite density (estimated by microscopy) was categorized as being either undetectable (z_{1} = 0), low (<1000 parasites per microlitre of blood, z_{1} = 1) or high (≥1000 parasites per microlitre of blood, z_{1} = 2). Hosts were classified into two different age categories (<6 years [children, z_{2} = 0] or ≥ 6 yr [older children, z_{2} = 1). Mosquito infection was allowed to vary between Burkina Faso and Kenya (z_{3} = 0 for Burkina Faso, z_{3} = 1 for Kenya). The function $\varphi \left({\widehat{x}}^{\prime},\eta \right)$ determines the shape of the relationship at very low gametocyte densities and is motivated by the observation that the malaria parasite can adjust its sex ratio to optimise transmission (Reece et al., 2008). Evidence indicates that transmission is possible even at very low gametocyte densities (Schneider et al., 2007; Ouédraogo et al., 2009), so $\varphi \left({\widehat{x}}^{\prime},\eta \right)$ allows infection to rise very quickly with increasing gametocyte density at a rate determined by parameter η,
Equation 4 was originally derived to describe the probability that a host would contain both male and female parasites according to the mean number of parasites and the aggregation (overdispersion) parameter of the negative binomial distribution (May, 1977). Reproduction in malaria is more complex than in the (helminth) system for which equation 4 was devised; hence, parameter η is not biologically interpretable.
To investigate whether mosquito infection is best predicted by gametocyte density on the arithmetic or logarithmic scale, both hypotheses were tested in the model, substituting an estimate of gametocyte density on the logarithmic scale ($\mathrm{ln}\widehat{{x}^{\prime}}$) for ${\widehat{x}}^{\prime}$ in equations 3 and 4 and comparing model fits. The model quantifying the uncertainty in gametocyte density estimates was fitted at the same time as the regression determining the relationship between gametocyte density (and other covariates) and the proportion of mosquitoes infected using Bayesian Markov Chain Monte Carlo methods. Fitting the models simultaneously enables the uncertainty in the gametocyte density estimates to be reflected in the uncertainty of the shape of the relationship.
Age profile of gametocyte density
Request a detailed protocolBlood samples were randomly collected from 412 hosts from a single village in Burkina Faso (Ouédraogo et al., 2010). QTNASBA assays were run on the samples, and the methods described above were used to convert assay results into estimates of gametocyte density. A description of these data is given in Figure 2—source data 1. To facilitate visual inspection of the age profile of gametocyte density, and to generate summary statistics, the function h(A) was fitted to these data to describe how the mean number of gametocytes per microlitre of blood changes with host age (A)
Parameters τ, ψ, and ω were estimated assuming a negative binomial error structure to account for the high degree of overdispersion (aggregation) observed in gametocyte density estimates (using an overdispersion parameter that did not change with host age or gametocyte density). As above, a hierarchical model was used to estimate the gametocyte density and its associated uncertainty at the same time as the age profile allowing the full uncertainty of the shape of the function to be expressed.
Age profile of the infectious reservoir
Request a detailed protocolThe probability that a mosquito biting a host of a particular age will become infected can be estimated using the above model together with estimates of the host gametocyte and asexual parasite density. This was done using the crosssectional data from Burkina Faso to generate a proxy for the human reservoir of infection at this specific time in this location (although it does not take into account how vector biting rate may vary with host age). It is important to use individual host estimates of gametocyte density instead of mean estimates as the high degree of gametocyte overdispersion among hosts may accentuate the nonlinear relationship between gametocyte density and mosquito infection. For each of the 412 hosts, point estimates were obtained of the percentage of feeding mosquitoes (taking a single full blood meal) that would develop oocysts according to the age and the gametocyte/asexual parasite density of the host (equation 3). Equation 5 was used to fit a relationship between host age and their contribution towards human to mosquito transmission. This bestfit line was used to predict the percentage of overall transmission that originated from hosts who had gametocyte density estimates below the detection threshold of microscopy. Here, we define the detection threshold as the minimum (positive) density that can be estimated by sampling a certain volume of blood. It is assumed that gametocytes are counted against 500 white blood cells (WBC) and that there are on average 8000 WBC per µl of blood. This gives an average detection threshold of 16 gametocytes per microlitre of blood. This is a conservative estimate of the threshold density since densities above this value may still give false negative results both due reading errors or random sampling of gametocytes on the microscope slide.
Care should be taken when interpreting the results as mosquito membrane feeding experiments were not performed with blood from adult hosts. Evidence indicates that antibodies associated with transmissionblocking activity may be lower in older age groups (Drakeley et al., 2006), which could mean that adult hosts could have a greater potential to contribute to overall transmission.
Fitting procedure
Request a detailed protocolBayesian Markov Chain Monte Carlo techniques were used to fit the models in OpenBUGS (Lunn et al., 2009). This approach allows the uncertainties in the assessments of gametocyte density and mosquito infection rates (generated by different numbers of mosquitoes being dissected) to be taken into consideration, propagated into parameter posterior distributions. All parameters were assigned uninformative priors and were run until convergence was reached using standard methodology (Gelman and Rubin, 1996). The most parsimonious models were selected by comparing deviance information criteria (DIC) values (the lowest value giving the most parsimonious yet adequate fit) (Spiegelhalter et al., 2002). Uncertainly around the bestfit line is indicated by showing the 95% range of 10,000 runs randomly sampled from the posterior distribution. The distribution of gametocytes between hosts was highly overdispersed, with an aggregation parameter of the negative binomial distribution of 0.185 (95% Bayesian credible interval [CI], 0.16–0.21).
References
 1
 2

3
Epidemiology of Plasmodium falciparum in a rice field and a savanna area in Burkina Faso: seasonal fluctuations of gametocytaemia and malarial infectivityAnn Trop Med Parasitol 85:377–385.

4
High human malarial infectivity to laboratorybred Anopheles gambiae in a village in Burkina FasoAm J Trop Med Hyg 48:700–706.
 5
 6

7
The aggressiveness of Anopheles gambiae A in relation to the age and sex of the human subjectsBull World Health Organ 56:147–154.
 8
 9

10
Data from: Predicting mosquito infection from Plasmodium falciparum gametocyte density and estimating the reservoir of infectionDryad Digital Repository, 10.5061/dryad.0k402.

11
A comparative study of thick and thin blood films in the diagnosis of scanty malaria parasitaemiaBull World Health Organ 34:249–267.
 12
 13
 14

15
Markov chain Monte Carlo methods in biostatisticsStat Methods Med Res 5:339–355.https://doi.org/10.1177/096228029600500402

16
The reservoir of Plasmodium falciparum malaria in a holoendemic area of western KenyaTrans R Soc Trop Med Hyg 86:355–358.https://doi.org/10.1016/00359203(92)90216Y
 17

18
Primaquine for reducing Plasmodium falciparum transmissionCochrane Database Syst Rev 9:CD008152.https://doi.org/10.1002/14651858.CD008152.pub2
 19
 20
 21

22
A Bayesiananalysis of the linear calibrationproblem—responseTechnometrics 23:344–350.https://doi.org/10.2307/1268224

23
Infectivity to mosquitoes of Plasmodium falciparum as related to gametocyte density and duration of infectionAm J Trop Med Hyg 4:781–789.
 24
 25
 26
 27

28
The malarial infectivity of an African village population to mosquitoes (Anopheles gambiae): a random xenodiagnostic surveyAm J Trop Med Hyg 6:971–979.
 29
 30
 31
 32

33
Data from: Substantial contribution of submicroscopical Plasmodium falciparum gametocyte carriage to the infectious reservoir in an area of seasonal transmissionDryad Digital Repository, 10.5061/dryad.hv01f.
 34

35
Relationships between host infectivity to mosquitoes and asexual parasite density in Plasmodium falciparumAm J Trop Med Hyg 75:32–37.

36
Submicroscopic Plasmodium falciparum gametocytaemia and the contribution to malaria transmissionPhD Thesis, Radboud University Nijmegen Medical Centre, Nijmegen, The Netherlands.

37
Submicroscopic Plasmodium falciparum gametocyte densities frequently result in mosquito infectionAm J Trop Med Hyg 76:470–474.
 38

39
Data from: Submicroscopic Plasmodium falciparum gametocyte densities frequently result in mosquito infectionDryad Digital Repository, 10.5061/dryad.589ft.
 40
 41

42
Bayesian measures of model complexity and fitJ Roy Stat Soc B Stat Meth 64:583–616.https://doi.org/10.1111/14679868.00353

43
Primaquine to prevent transmission of falciparum malariaLancet Infect Dis 13:175–181.https://doi.org/10.1016/S14733099(12)701986

44
Updated WHO policy recommendation (October 2012): single dose primaquine as a gametocytocide in Plasmodium falciparum malariahttp://www.who.int/malaria/pq_updated_policy_recommendation_en_102012.pdf, (accessed Dec 14, 2012).
Decision letter

Prabhat JhaReviewing Editor; University of Toronto, Canada
eLife posts the editorial decision letter and author response on a selection of the published articles (subject to the approval of the authors). An edited version of the letter sent to the authors after peer review is shown, indicating the substantive concerns or comments; minor concerns are not usually shown. Reviewers have the opportunity to discuss the decision before the letter is sent (see review process). Similarly, the author response typically shows only responses to the major concerns raised by the reviewers.
Thank you for sending your work entitled “Targeting Malaria Transmission Control: Predicting Mosquito Infectivity from Plasmodium falciparum Gametocyte Density” for consideration at eLife. Your article has been favorably evaluated by Prabhat Jha as Senior editor and 3 reviewers, one of whom was Vasee Moorthy.
The reviewers discussed their comments before we reached this decision and the Senior editor has assembled the following major substantive comments to help you prepare a revised submission.
1) Estimating Gametocyte Density From QTNASBA Assays
I do not know the inverse regression problem well and can appreciate that there are alternative approaches. The supplementary file indicates some uncertainty in the choice of this method “for discussion see Osbourne, 1991” but this is not explained. It would be useful to give brief reasons why the particular method was chosen.
The authors say that QTNASBA may detect only female gametocytes and elsewhere that the sex ratio of P falciparum can change to optimise transmission. This may alter the shape of the relationship between QTNASBA results and the estimated gametocyte densities.
It would be useful to validate the model for estimating gametocyte densities from QTNASBA assays for data other than that to what was fitted.
Figure 1A shows a double hump in the relationship between gametocytes in the blood sample and mosquito infectivity. This seems unintuitive for a biological process. I wondered to what extent this might be forced by the functions chosen. Does adding extra flexibility result in the same double hump?
Figure 1: An elegant presentation of this data. In the inset to part A there is a large clump of points on the x axis and it is not clear that the fit is good enough to be useful. For me this graph highlights the importance of understanding how different the line is in this model to other lines, and how much better this model fits, than the models that are the next best fit.
Please provide some discussion of this in the text (in addition to supplementary information). A listing of alternative models with their goodness of fit on a quantitative basis should be included in the manuscript itself.
The uncertainty range is huge, and ranges from very little increase in infectivity with increasing gametocyte density, to very steep increases in infectivity with increasing gametocyte density. What steps do the authors propose are taken to reduce the uncertainty range?
2) Mosquito Infectivity
Equation A4 includes a parameter called kappa. This is confusing since kappa is conventionally used to represent the infectious reservoir.
Results: The results for age are presented as if there was evidence of an effect, but the confidence intervals suggest otherwise.
The rationale for adjusting for asexual parasite density was not clear. It “improved model fit” but the aim here is not the best fitting model so much as to describe the relationship between estimated gametocyte density and infectivity. Including asexual parasite density as a main effect modifies the interpretation of the parameters.
3) The Paper Estimates the Relationship Between Gametocyte Densities and Age
The data comes from a crosssectional survey of 412 people in a village in Burkina Faso. There is likely to be seasonal variation in parasite densities and the relative densities in adults and children may vary seasonally due to different levels of acquired immunity. It should be stated when the survey was carried out and also that a single survey is a potential limitation.
4) Efficacy Definition
A major issue is to be quite clear what is meant by efficacy throughout the paper. This lack of clarity permeates this entire technical area (i.e., not only this manuscript) and causes a lot of confusion. I believe most or all instances of use of the word efficacy in this manuscript refers to a % reduction in the gametocyte density. How does this relate to % reduction in proportion of infected mosquitoes and how will this translate to reduced incidence of human infection? This is a critical discussion that is omitted. Discussion of efficacy without some idea what this translates to in terms of reduction in human infection is rather misleading.
If the authors mean reduced transmission from humans to mosquito, I suggest they use a narrative phrase, such as this and do not use the word efficacy at all.
5) Data/Model Description
The data the model was fitted to was not well described (“repeatedly testing samples with known gametocyte density”). I also had trouble working out whether the fitting was simultaneous with the model for mosquito infectivity or not. The text mentions that the “estimates were used in the functions below” but the table presents all parameters together. My concern is that that the error in the estimated gametocyte densities, from the QTNASBA results and also the chance of gametocytes being contained in the blood sample, is incorporated.
The model and parameters are not explained. The methods section merely mentions “a mathematical model” and the supplementary file gives an equation involving three parameters, but it is not clear how the equation was derived or what the parameters represent. A smooth curve is produced, but the relationship does not necessarily have to be smooth: for example, pregnant women may have higher parasitaemia and affect the agecurves.
Figure 2A suggests that the relationship does not fit well for two of the age groups. This should be acknowledged and possible reasons mentioned.
A negative binomial error structure was used. The Figure 2—source data 1 table suggests that there were some zeros in the data and that the densities are skewed. Did the residuals suggest a good fit? If there are many zeros, then a zeroinflated negative binomial model may be more appropriate.
6) Acknowledge Limitations
On a couple of occasions, the authors make statements that are not backed up by evidence from this paper.
Discussion. It is implied that this work gives direct evidence of the effectiveness of parasite strategies at low gametocyte densities such altering the sex ratio, but no direct information on strategies is obtained.
Although the paper does include an elegant sampling approach to quantify uncertainties, Although the paper does include an elegant sampling approach to quantify uncertainties, there should be a beefed up discussion of the uncertainties underlying this work. The authors should provide more discussion on what is known about which parameters in the input data are driving uncertainty, how much confidence they have that the model they chose is superior to alternative models, and on what basis they made this decision. There should also be a discussion on the limitations related to the differences between mosquitoes used in feeding assays, and wildtype mosquitoes.
Although uncertainty is presented in one of the figure panels, and this is very helpful, uncertainty ranges/confidence intervals should be quoted in the text where figures are provided.
We suggest that the authors consider how the paper can act as an aid to conceptualising the components underlying persontoperson transmission in field settings, and to gaining a better understanding of the uncertainties related to the component of transmission that their results speak to, and also to gain a better understanding of what the existing data gaps are that must be filled in order that the research community can determine how to test new transmissionreducing interventions, and quantify their effect on persontoperson transmission.
https://doi.org/10.7554/eLife.00626.010Author response
1) Estimating Gametocyte Density From QTNASBA Assays
I do not know the inverse regression problem well and can appreciate that there are alternative approaches. The supplementary file indicates some uncertainty in the choice of this method “for discussion see Osbourne, 1991” but this is not explained. It would be useful to give brief reasons why the particular method was chosen.
Inverse regression has been widely discussed in the statistical literature over the last 30 years. The majority of the controversy stems from when the gradient of the calibration line is very small (i.e., a small difference in TTP spans the whole range of gametocyte densities). This is not the case for our QTNASBA results where the gradient is sufficiently large and no appreciable differences in the results are observed irrespective of the method used. To explain this fully to readers, we have included the following section in Supplementary file 1:
“Here we adapt this approach for use in a Bayesian hierarchical model. A number of different methods have been used to do this (see Hoadley (1970) and Hunter & Lamboy (1981) and for discussion see Osborne (1991)), each of which has a different approach for dealing with the problem of low beta1 values (i.e., a gentle gradient, which when used in equation A2 can generate infinitely large confidence interval estimates). Here we use the most parsimonious approach that does not require assignment of prior distributions to (the unknown) gametocyte densities. Rather, uncertainty in the estimated regression coefficients of equation A1 is propagated numerically via equation A2 to yield uncertainty in the estimated gametocyte densities. Since each of the calibration lines have a relatively steep gradient the difference between the different methods will be relatively minor and all will generate sensibly tight confidence interval estimates.”
The authors say that QTNASBA may detect only female gametocytes and elsewhere that the sex ratio of P falciparum can change to optimise transmission. This may alter the shape of the relationship between QTNASBA results and the estimated gametocyte densities.
The reviewers are correct in saying that if QTNASBA does only detect female gametocytes then the relationship between gametocyte density (as estimated using another method which does not use Pfs25 mRNA such as microscopy) and mosquito infectivity might be different. However, we would suggest that in such a situation the number of female gametocytes would be a better prediction of infectivity than total gametocytes as each male gametocyte can fertilise up to 7 female gametocytes and the sex ratio will likely change to maximise female fertilization. Therefore, female gametocyte density represents the maximum potential infectivity.
We have changed the text accordingly:
“Even if the current QTNASBA method does only detect female gametocytes, it will not change the qualitative conclusions of this study, nor their application if other studies use the same technique. However, direct comparison of our results with those of studies using different methods for quantifying gametocyte density should be aware that the shape of the relationship may differ.”
It would be useful to validate the model for estimating gametocyte densities from QTNASBA assays for data other than that to what was fitted.
We are unsure of the exact purpose of this validation procedure. If the reviewers meant to validate the QTNASBA against another highly accurate method of determining gametocyte density, then we agree that this is a necessary process for us to have full confidence in the technique and our method of correcting for its uncertainty. However, we feel that our dilution procedure (whereby a gametocyte sample is diluted 5 times and each generates a TTP) provides a strong validation procedure within the model. In this example the fit of the calibration line to these dilutions is relatively poor, which is the reason behind the uncertainty round the individual gametocyte density estimates being so high. We are in the process of refining the nucleic acid extraction and amplification technique and we are generating increasingly precise QTNASBA estimates, as validated by the dilution procedure. Since this has been done on samples not related to the current study we feel that including these here would confuse the message. If the reviewers are suggesting that we repeat the experiment with a different dataset then we do agree, and we have highlighted this in the text. However, to our knowledge no dataset currently exists.
Figure 1A shows a double hump in the relationship between gametocytes in the blood sample and mosquito infectivity. This seems unintuitive for a biological process. I wondered to what extent this might be forced by the functions chosen. Does adding extra flexibility result in the same double hump?
We would argue that the double hump is consistent with our current knowledge of the biological process and what has been observed in laboratory models. The first hump is a necessity for the line to go through point 0,0 (which it has to do as zero gametocytes cannot generate any oocysts) since ∼4% of mosquitoes are infected at very low gametocyte densities. The second hump is less certain since there are relatively few data points at very high gametocyte densities. There is evidence of a second plateau in these data using differently shaped functions, which is why the power and hyperbolic functions give better fits than the linear function. This is also consistent with what has been observed in animal models (see Sinden et al. 2009). That said, we agree that the functions chosen might influence the shape and that it could be more complicated than shown here.
More complex functions were tested (such as asymmetrical logistic functions), though they failed to converge adequately due to the additional parameters and the high uncertainty in gametocyte density. We have added a comment to stress that a wider range of functions should be tested when the accuracy of the data allows:
Discussion: “Further standardising the membrane feeding assay, improving the accuracy of QTNASBA technique and the inclusion of additional covariates (such as estimates of gametocyte maturity which have been investigated in P. vivax (Chansamut et al., 2012)) would improve the accuracy of the relationship between gametocyte density and mosquito infectivity. It would permit to see if the patterns described here are consistent across time and space, and enable a wider range of functional forms to be tested.”
Figure 1: An elegant presentation of this data. In the inset to part A there is a large clump of points on the x axis and it is not clear that the fit is good enough to be useful. For me this graph highlights the importance of understanding how different the line is in this model to other lines, and how much better this model fits, than the models that are the next best fit.
Please provide some discussion of this in the text (in addition to supplementary information). A listing of alternative models with their goodness of fit on a quantitative basis should be included in the manuscript itself.
We agree that this information should be presented in the main text and we have added the following sentences to the Results section:
“Infectivity rises rapidly with increasing gametocyte density and by 1 gametocyte µl1 ∼4% of all mosquitoes develop oocysts. The subsequent increase in infectivity with increasing gametocyte density is best described by the Gompertz model (deviance information criterion, DIC = 1034), which gave a significantly better fit than the linear (DIC=1073), power (DIC=1059) or hyperbolic (DIC=1062) functions (Figure 1A). The best fit model predicts that increasing density from 1 to 200 gametocytes µl1 does not appreciably increase infectivity.”
The uncertainty range is huge, and ranges from very little increase in infectivity with increasing gametocyte density, to very steep increases in infectivity with increasing gametocyte density. What steps do the authors propose are taken to reduce the uncertainty range?
We agree that there is considerable uncertainty in the shape of the fitted curves reflecting the uncertainty in gametocyte density and the variability seen in mosquito feeding assays. We have included a few sentences outlining how the uncertainty can be minimised:
“The infectivity of mosquitoes to Plasmodiuminfected blood is known to vary substantially within and between studies. This is in part due to methodological issues with the membrane feeding assay (Bousema, Churcher, Morlais, & Dinglasan, 2013) but also related to biological differences in the parasitevector combination and the blood environment (Bousema & Drakeley, 2011). Further standardising the membrane feeding assay, improving the accuracy of QTNASBA technique and the inclusion of additional covariates (such as estimates of gametocyte maturity which have been investigated in P. vivax (Chansamut et al., 2012)) would improve the accuracy of relationship between gametocyte density and mosquito infectivity.”
2) Mosquito Infectivity
Equation A4 includes a parameter called kappa. This is confusing since kappa is conventionally used to represent the infectious reservoir.
We have changed this parameter to eta to avoid confusion.
Results: The results for age are presented as if there was evidence of an effect, but the confidence intervals suggest otherwise.
Including age significantly improves the fit of the model though we appreciate that the 95% Bayesian credible intervals span 0. We have augmented the text in order to reflect this: “Including age improves the fit of the model suggesting age is an important confounder though the Bayesian credible intervals include zero.”
The rationale for adjusting for asexual parasite density was not clear. It “improved model fit” but the aim here is not the best fitting model so much as to describe the relationship between estimated gametocyte density and infectivity. Including asexual parasite density as a main effect modifies the interpretation of the parameters.
We appreciate that adjusting for asexual parasite density (as with age) may alter the shape of the relationship between gametocytes and mosquito infectivity, as asexual parasite density changes with increasing gametocyte density. However, the relationship between the density of gametocytes and asexual parasites is unlikely to remain constant (certainly not after certain drug treatments), so we feel that it is important to understand the impact of increasing gametocyte density without the influence of confounding variables. It also allows more accurate predictions of mosquito infectivity from crosssectional surveys where these covariates are collected (such as in Burkina Faso). Since asexual parasite density and age are routinely collected we felt that the extra complexity of the model was warranted. We thank the reviewers for bringing this to our attention and we have changed the end of the Introduction accordingly to stress the rationale for doing so: “Other covariates that may influence infectivity such as asexual parasite density (measured by microscopy) and host age were also included to increase the ability of the model to predict human to mosquito transmission.”
3) The Paper Estimates the Relationship Between Gametocyte Densities and Age
The data comes from a crosssectional survey of 412 people in a village in Burkina Faso. There is likely to be seasonal variation in parasite densities and the relative densities in adults and children may vary seasonally due to different levels of acquired immunity. It should be stated when the survey was carried out and also that a single survey is a potential limitation.
Thank you for highlighting this oversight. We have included a sentence to highlight the dangers of crosssectional surveys in the main text. In the Results: “In the crosssectional survey in Burkina Faso children harbour more gametocytes than adults, with 10year olds having 5 times as many gametocytes compared to 50year olds (Figure 2A).” In the Discussion: “The reservoir of infection is likely to vary between endemicity settings and according to local malaria control practices. It is also likely to vary over time so multiple crosssectional surveys may be required.”
4) Efficacy Definition
A major issue is to be quite clear what is meant by efficacy throughout the paper. This lack of clarity permeates this entire technical area (i.e., not only this manuscript) and causes a lot of confusion. I believe most or all instances of use of the word efficacy in this manuscript refers to a % reduction in the gametocyte density. How does this relate to % reduction in proportion of infected mosquitoes and how will this translate to reduced incidence of human infection? This is a critical discussion that is omitted. Discussion of efficacy without some idea what this translates to in terms of reduction in human infection is rather misleading.
If the authors mean reduced transmission from humans to mosquito, I suggest they use a narrative phrase, such as this and do not use the word efficacy at all.
We appreciate the ambiguity of the term “efficacy” in the literature and we thank the reviewers for highlighting it here. As suggested we have changed the terminology so that each time the word efficacy is used we clarify which lifestage is targeted using a narrative phrase. We have also included the following sentences in the Discussion, stating:
“Figure 1B can be used to estimate how a reduction in the number of gametocytes will equate to a reduction in the proportion of mosquitoes becoming infected, and hence mosquito to human transmission. How this relates to the incidence of malaria and subsequent disease will depend, in part, on the degree of immunity in the human population.”
5) Data/Model Description
The data the model was fitted to was not well described (“repeatedly testing samples with known gametocyte density”). I also had trouble working out whether the fitting was simultaneous with the model for mosquito infectivity or not. The text mentions that the “estimates were used in the functions below” but the table presents all parameters together. My concern is that that the error in the estimated gametocyte densities, from the QTNASBA results and also the chance of gametocytes being contained in the blood sample, is incorporated.
The uncertainty in the QTNASBA technique includes the variability in the number of gametocytes in the blood sample. We have added a number of sentences stressing this and that the model for quantifying gametocyte density was fit at the same time as the model determining mosquito infectivity. This allows the uncertainty in gametocyte density to be reflected in the uncertainty of the bestfit line (and other covariates).
We have changed the “Mosquito infectivity” section of the methods:
“Gametocyte density was quantified using QTNASBA and the uncertainty in the density estimates (the intraassay variability) was quantified by fitting a hierarchical model to the 16 independent dilution series. In the Burkina dataset multiple assays were carried out on the same blood sample (an average of 2.63 assays per unknown blood sample). The mean gametocyte density from these multiple assays was taken to increase the accuracy of the estimates. The precise shape of the relationship was determined by fitting a range of different functional forms (a modified constant, linear, power, hyperbolic, sigmoid and Gompertz functions) and statistically determining which gave the best fit (a full description of the model is given in the Supplementary files). The model quantifying the uncertainty in gametocyte density estimates was fitted at the same time as the regression determining the relationship between gametocyte density (and other covariates) and the proportion of mosquitoes infected using Bayesian Markov Chain Monte Carlo methods. Fitting the models simultaneously enables the uncertainty in the gametocyte density estimates to be reflected in the uncertainty of the best fit model.”
The model and parameters are not explained. The methods section merely mentions “a mathematical model” and the supplementary file gives an equation involving three parameters, but it is not clear how the equation was derived or what the parameters represent. A smooth curve is produced, but the relationship does not necessarily have to be smooth: for example, pregnant women may have higher parasitaemia and affect the agecurves.
We have explained the age profile model in greater detail in the main text. We appreciate that certain groups, such as pregnant women, may have higher parasite densities, though we do not have more detailed information from the crosssectional survey to include additional covariates within the model. Nevertheless, we feel that looking at things at a population level using smoothed curves is important for policy makers. The age profile takes into account how these heterogeneities, such as pregnancy status, changes with age and will result in elevated parasite levels in particular age groups.
The following section was included in the “Age profile of gametocyte density” section of the methods:
“A three parameter functional form was used which can describe an age profile which increases with age or peaks at an intermediate age. As above, a hierarchical model was used to estimate the gametocyte density and its associated uncertainty at the same time as the age profile allowing the full uncertainty of the shape of the function to be expressed (a full description of the data and model is given in the Supplementary files).”
Figure 2A suggests that the relationship does not fit well for two of the age groups. This should be acknowledged and possible reasons mentioned.
We agree that these two data points could be closer to the line, though the confidence intervals of the point estimate and the model do overlap. The functional form is flexible enough to pass much closer to these points though doing so generates a worse fit to these data. We feel that the main reason for the apparent discrepancy is the sampling from the highly overdispersed data, where missing the rare very high parasite densities by chance would result in lower mean estimates. Neither of these data points have particularly large sample sizes, which increases the probability of this happening. This overdispersion is taken into account in the fitting process and though other methods of central tendency could be presented in the figure (such as the geometric mean or median), since the model is showing the arithmetic mean we felt this was the best comparison. To highlight this we have added the following sentence to the legend of Figure 2: “In Figure 1A the [1015) and [1520) groups appear lower than the best fit line due to sampling artefacts generated by the highly overdispersed data.”
A negative binomial error structure was used. The Figure 2—source data 1 table suggests that there were some zeros in the data and that the densities are skewed. Did the residuals suggest a good fit? If there are many zeros, then a zeroinflated negative binomial model may be more appropriate.
Though there are a lot of zero counts, this is adequately described by the simple negative binomial distribution. A zeroinflated negative binomial distribution was fitted to the estimates of gametocyte density (used to generate the gametocyte age profile) using maximum likelihood, though the improved fit didn’t warrant the extra complexity (likelihood ratio test p value = 0.7). Diagnostics on the residuals indicate no major issues other than the high level of uncertainty. We have changed the following section of the Discussion to reflect this: “However, the highly overdispersed distribution of gametocytes between hosts (adequately described by the negative binomial distribution) means that in the population from Burkina Faso…”
6) Acknowledge Limitations
On a couple of occasions, the authors make statements that are not backed up by evidence from this paper.
We have removed the misleading statements such as those highlighted below.
Discussion. It is implied that this work gives direct evidence of the effectiveness of parasite strategies at low gametocyte densities such altering the sex ratio, but no direct information on strategies is obtained.
We were trying to suggest that transmission at very low gametocyte densities is evidence that the strategies of the parasite to maximise transmission were working. However we agree that this sentence could be misleading and so we have deleted it from the Discussion.
Although the paper does include an elegant sampling approach to quantify uncertainties, there should be a beefed up discussion of the uncertainties underlying this work. The authors should provide more discussion on what is known about which parameters in the input data are driving uncertainty, how much confidence they have that the model they chose is superior to alternative models, and on what basis they made this decision. There should also be a discussion on the limitations related to the differences between mosquitoes used in feeding assays, and wildtype mosquitoes.
We agree that the uncertainties and the difference between laboratory and wild mosquitoes should be stressed further. We have added the following section to the Discussion:
“Gametocyte density measurement error is likely to cause the majority of the uncertainty seen in model outputs though the complexity of the membrane feeding assay is also likely to contribute. The infectivity of mosquitoes to Plasmodium infected blood is known to vary substantially within and between studies. This is in part due to methodological issues with the membrane feeding assay (Bousema, Churcher, Morlais, & Dinglasan, 2013) but also related to biological differences in the parasitevector combination and the blood environment (Bousema & Drakeley, 2011). Further standardising the membrane feeding assay, improving the accuracy of QTNASBA technique and the inclusion of additional covariates (such as estimates of gametocyte maturity which have been investigated in P. vivax (Chansamut et al., 2012)) would improve the accuracy of the relationship between gametocyte density and mosquito infectivity. It would also permit to see if the patterns described here are consistent across time and space, and enable a wider range of functional forms to be tested. Care should also be taken when interpreting the results of feeding assays as the mosquitoes used and their biting behaviour is likely to be different from that seen in wild mosquitoes in field situations (Bousema et al., 2013).”
Although uncertainty is presented in one of the figure panels, and this is very helpful, uncertainty ranges/confidence intervals should be quoted in the text where figures are provided.
We agree and have added 95% Bayesian Credible Intervals where appropriate.
We suggest that the authors consider how the paper can act as an aid to conceptualising the components underlying persontoperson transmission in field settings, and to gaining a better understanding of the uncertainties related to the component of transmission that their results speak to, and also to gain a better understanding of what the existing data gaps are that must be filled in order that the research community can determine how to test new transmissionreducing interventions, and quantify their effect on persontoperson transmission.
We very much agree with the reviewers that this is an important topic. Persontoperson transmission is highly complex due to the complexities of the human immune response. Since this is not the topic of this manuscript, we concentrate our answers on the reservoir of infection (i.e., transmission, not necessarily incidence of malaria). We have augmented the following section of the Discussion to reflect our opinions on the subject:
“The contribution of different age groups to overall transmission will depend on local demography, agedependent protection by malaria interventions (such as the use of bednets or time spent in the house and therefore personally protected by indoor residual spraying), and human to mosquito contact patterns. The last two of these are difficult to measure and poorly understood, reducing our ability to predict accurately the sources of infection. The reservoir of infection is likely to vary between endemicity settings and over time so multiple crosssectional surveys may be required.”
https://doi.org/10.7554/eLife.00626.011Article and author information
Author details
Funding
European Commission FP7 Collaborative Project (HEALTHF32008223736)
 Thomas S Churcher
 MaríaGloria Basáñez
AFIRM grant from the Bill & Melinda Gates Foundation (OPP1034789)
 Teun Bousema
 Chris Drakeley
 André Lin Ouédraogo
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Reviewing Editor
 Prabhat Jha, University of Toronto, Canada
Publication history
 Received: February 11, 2013
 Accepted: April 9, 2013
 Version of Record published: May 21, 2013 (version 1)
Copyright
© 2013, Churcher et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics

 2,185
 Page views

 399
 Downloads

 101
 Citations
Article citation count generated by polling the highest count across the following sources: Crossref, Scopus, PubMed Central.
Download links
Downloads (link to download the article as PDF)
Download citations (links to download the citations from this article in formats compatible with various reference manager tools)
Open citations (links to open the citations from this article in various online reference manager services)
Further reading

 Epidemiology and Global Health

 Epidemiology and Global Health
 Microbiology and Infectious Disease