 Methodology
 Open Access
 Published:
Measures of clinical malaria in field trials of interventions against Plasmodium falciparum
Malaria Journalvolume 6, Article number: 53 (2007)
Abstract
Background
Standard methods for defining clinical malaria in intervention trials in endemic areas do not guarantee that efficacy estimates will be unbiased, and do not indicate whether the intervention has its effect by modifying the force of infection, the parasite density, or the risk of pathology at given parasite density.
Methods
Three different sets, each of 500 Phase IIb or III malaria vaccine trials were simulated corresponding to each of a preerythrocytic, blood stage, and antidisease vaccine, each in a population with 80% prevalence of patent malaria infection. Simulations considered only the primary effects of vaccination in a homogeneous trial population. The relationships between morbidity and parasite density and the performance of different case definitions for clinical malaria were analysed using conventional likelihood ratio tests to compare incidence of episodes defined using parasite density cutoffs. Bayesian latent class models were used to compare the overall frequencies of clinical malaria episodes in analyses that did not use diagnostic cutoffs.
Results
The different simulated interventions led to different relationships between clinical symptoms and parasite densities. Consequently, the operating characteristics of parasitaemia cutoffs in general differ between vaccine and placebo arms of the simulated trials, leading to different patterns of bias in efficacy estimates depending on the type of intervention effect. Efficacy was underestimated when low parasitaemia cutoffs were used but the efficacy of an asexual blood stage vaccine was overestimated when a high parasitaemia cutoff was used. The power of a trial may be maximal using case definitions that are associated with substantial bias in efficacy.
Conclusion
Secondary analyses of the data of malaria intervention trials should consider the relationship between clinical symptoms and parasite density, and attempt to estimate overall numbers of clinical episodes and the degree of bias of the primary efficacy measure. Such analyses would help to clarify whether the effect of an intervention corresponds to that anticipated on the basis of the parasite stage that is targeted, and would highlight whether the primary measure of efficacy results from unexpected behaviour in the parasitological and clinical data used to estimate it.
Background
In endemic settings malaria usually presents with rather nonspecific symptoms, such as fever, and not all sick individuals with malaria parasites are really suffering from clinical malaria. This is because most of the population may be infected with Plasmodium falciparum parasites, without these causing any acute illness. It follows that the presence of parasites in a sick person does not necessarily mean that malaria is the cause of the illness. In field trials of novel interventions, estimates of efficacy need to be made using case definitions with high specificity. Otherwise efficacy will be underestimated.
The greater the parasite density in the blood the more reasonable it is to assume that an illness is caused by malaria, so a case definition for clinical malaria for use in a trial can be obtained by defining a parasite density cutoff specific for the surveillance mechanism of choice (local health centre, hospital, active case detection). The sensitivity and specificity of different parasite density cutoffs can be obtained by modeling the excess risk of fever as a function of parasite density [1, 2], where the comparator is the risk in aparasitaemic individuals. This has been used in a number of trials to decide upon an appropriate cutoff [3–5]. Ideally this analysis is carried out in the same population (and age groups) as the vaccine trial and using the same morbidity surveillance system, since relationships between morbidity and infection depend on age and immune status [6–8].
This approach has been endorsed by the WHO Study Group on Measures of Malaria Vaccine Efficacy for obtaining case definitions for use in pivotal trials of malaria vaccine [9]. Efficacy estimates based on this algorithm can be easily obtained using standard software and are appropriate for defining the primary outcome for trials aiming to achieve registration.
Applying such an algorithm though does not guarantee that efficacy estimates will be unbiased, and does not provide an interpretation of how a vaccine is acting. As secondary objectives of malaria vaccine trials, investigators should be interested also in drawing inferences about whether the vaccine acts in accordance with its design and how it interacts with natural immunity. This paper uses simulations of trials to consider the theoretical performance of this method for different kinds of vaccines, and suggests a range of additional exploratory analyses that can be carried out in order to better understand vaccine action. The simulations consider the likely effects of different kinds of vaccines but the same approach is applicable to the analysis of the clinical impact of any effective intervention against malaria.
Methods
Simulations
Phase IIb or III malaria vaccine trials were simulated assuming the study population to comprise children with an 80% prevalence of P. falciparum malaria, a distribution of parasite densities in the population, θ_{ c }(x), (in the absence of vaccine or in a placebo group) as shown in the thick line Figure 1a. The simulated frequency of disease in the placebo group at different parasite densities (relative to the risk in aparasitaemic individuals) is shown in the thick line in Figure 1b, with 45.5% of the clinical episodes attributable to malaria (Table 1). Three different hypothetical vaccine effects were simulated (Table 1):
Vaccine A: (Preerythrocytic vaccine). This vaccine reduces the force of infection by 50%. For simplicity, this is assumed to be reflected in a 50% reduction in the proportion of individuals who are infected across the whole range of parasite densities, and a random sample of 50% of the clinical malaria episodes are assumed to be averted.
Vaccine B: (Asexual blood stage vaccine). This vaccine reduces parasite densities by 50%, but does not affect the number of individuals who are infected, or the parasite densities at which they become ill.
Vaccine C: (Antidisease vaccine). This vaccine is assumed to have no effect on parasite densities but to lead to an increase in the parasite density at which clinical malaria occurs.
For each of these models of vaccination 500 simulated trials were run. Each simulated trial had a standard design with (i) an equal number of placebo and vaccine recipients; (ii) a total of 100 parasitological slides in each arm used to estimate the effect on the parasite density distribution in the population, and (iii) a standardised clinical surveillance to detect fever cases, with a total of 100 cases expected in unvaccinated individuals.
Datasets were constructed from these simulations using 19 categories of parasite density. The community parasitaemia data were generated by defining multinomial probabilities for each of the 19 categories based on the data of Figure 1a and randomly drawing a sample of 100 parasite densities from this distribution. The relative frequencies in Figure 1b were used to define the expected number of clinical cases for each of these 19 categories, on the assumption that the expected total number of clinical cases in the placebo arm was 100. The simulated number of cases in each category was then drawn from a Poisson distribution.
Analyses of effects of simulated vaccines
The effects of the different hypothetical vaccines were summarized both by considering the theoretical effects on θ_{ c }(x), the parasite density distribution in the community (Figure 1a), and on θ_{ s }(x) the density in sick individuals (Figure 1b). The latter arises as a mixture:
θ_{ s }(x) = λθ_{ m }(x) + (1λ)θ_{ c }(x)
where θ_{ m }(x) is the parasite density distribution in patients with clinical malaria and λ is the malaria attributable fraction. θ_{ s }(0), θ_{ c }(0) are then the frequencies of the uninfected (aparasitaemic) classes, and the relative risk of given parasite density x in sick individuals relative to controls (Figure 1c) is given by $\frac{{\theta}_{s}(x){\theta}_{c}(0)}{{\theta}_{c}(x){\theta}_{s}(0)}$. At each parasite density x, λ(x) is the malaria attributable proportion of clinical cases at x, equal to:
Figure 1d gives values for this quantity for each of the simulated populations and Figure 1e gives the distributions of θ_{ m }(x). The sensitivity of a cutoff X is then given by:
(Figure 1f) and the specificity (Figure 1g) by:
To compare the efficacy in averting clinical episodes that might be estimated in the simulated trials two distinct analyses of all 500 sets of simulated trials were carried out:

1.
In the first analysis, efficacy was estimated as the median of the sample of 1I_{V}(X)/I_{P}(X) in the 500 simulated trials, where where I_{V}(X) is the number of simulated cases in vaccinees, with parasite density>X; and I_{P}(X) is the number of simulated cases in the placebo arm (Figure 1h). Each of the category boundaries (X) used to sample the data was used in turn as a cutoff, and the average of the resulting efficacy estimates plotted against the cutoff (Figure 1h).
The power of these analyses was then estimated by determining the proportion of the sample of 500 simulated trials that gave significant efficacy as assessed using using binomial likelihood ratio tests (two sided) of the null hypothesis I_{V}(X) = I_{P}(X) (significance level α = 0.05).

2.
The second analysis sought to estimate the number of clinical malaria cases in each arm of each of the 500 simulated trials, by assigning a probability to each fever case, as a function of parasite density, rather than by using a dichotomous classification. These probabilities were estimated using a latent class model [10, 11] as previously described. This algorithm, implemented separately for vaccine and placebo, involved comparing the parasite density distributions in the simulated cases with those in the community samples drawn from θ_{ c }(x).
Results
The parasite density distributions in the community (Figure 1a) vary as straightforward consequences of the primary effects of vaccination. The preerythrocytic vaccine (A) halves the number of slide positive individuals in each category (a simplification of what we expect in a field study, where superinfection may occur); the effect of the simulated asexual blood stage vaccine (B) is more complicated, as it disproportionately reduces the frequency of high parasite densities, and slightly increases the frequency of very low parasite densities by shifting each individual to a lower density (Figure 1a). One consequence of this is that the highest density class is not represented among individuals who receive vaccine B, since any individual who would have been in this density class is now in the second highest class. The simulated antidisease vaccine (C) has no effect on the parasite density distribution in the community.
The distributions of parasite densities in the clinical cases differ among the three vaccines (Figure 1b). For all the vaccines there is a background incidence of nonmalaria disease, which is assumed to occur independently of the parasite density (left hand side of figure 1b), corresponding to nonmalaria illness and is the same in all groups. For all three vaccines fewer cases are expected at each positive value of the density distribution than occur in the placebo group.
These differences in parasite density distributions lead in turn to different relationships between incidence of disease and the community parasite density distribution, depending on the action of the vaccine (Figure 1c), and hence to different curves for the relationship between the attributable fraction, the frequencies of clinical cases with different densities, and the operating characteristics of case definitions (Figure 1defg, Table 2). The relative risk of a given parasite density among cases, relative to the risk in controls, is the same for the antiblood stage vaccine B as for the placebo arm (because the risk of disease, conditional on the parasite density, is the same in both arms, and the number of disease cases with no parasites is unchanged by the vaccine). For vaccine A, the proportion of cases at any given positive density is lower than in the corresponding proportion of cases in placebo recipients, because more of the cases with nonmalaria etiology are now aparasitaemic, so the relative risk of a given parasite density among cases relative to controls is lower than in placebo (Figure 1c). For vaccine C the curve in Figure 1c is also lower than that for placebo, but this is because there are fewer parasitaemic cases there is no change in the number of aparasitaemic ones.
Similarly, and as a direct consequence of the curves shown in Figure 1c, for these vaccines A and C, but not for vaccine B, the proportion of disease cases attributable to malaria at any given density is less in the active than in the placebo arm (Figure 1d).
The relative frequencies of malaria cases at different parasite densities (Figure 1e) show similar patterns to those of the relative frequencies of all disease cases (Figure 1d), but instead of intersecting at a nonzero point on the vertical axis, the plots pass through the origin, since clinical malaria cannot occur in the absence of parasites.
Integration of the curves in Figure 1e (Equation 3) then gives the curves for the sensitivity of parasite density cutoffs. There are clear differences between the vaccines. At any parasite density the sensitivity for vaccine A is equivalent to that for placebo, but for vaccines B and C it falls below that of placebo. The specificity, in contrast, is the same as placebo for vaccines A and C, but is higher than placebo in vaccine B (Figure 1g).
These differences in sensitivity and specificity have effects on the estimation of efficacy. At high values of x, corresponding to high specificity, the mean efficacy estimate (of the 500 simulated trials) for vaccines A and C approaches the true efficacy, while for vaccine B (where specificity in the vaccine arm is higher than in the placebo arm) the efficacy is overestimated. The proportion of trials giving statistically significant results (Figure 1i) (assuming them to have been analysed using a fixed cutoff) gives the power of the study. The power of the trials of vaccines A and C showed maxima at relatively low cutoffs, indicating that different cutoffs must be used if the aim is to avoid bias in the estimate of efficacy, from those used to optimise power.
Discussion
Field trials of interventions against malaria need to have easily interpretable primary outcome measures in order to make an impact on regulatory and policy decisions. At the same time, field trials represent the main opportunity for experimental study of immunoepidemiology of malaria and need to be fully exploited to further understanding of the mechanisms of action of the interventions. The analyses demonstrated in this paper are intended to contribute to plans for such secondary analyses.
The three hypothetical vaccines simulated in this study represent limiting cases of the effects of different interventions on clinical malaria. They do not correspond on a onetoone basis to real vaccines, but rather to possible intervention effects. Any real intervention might have secondary effects on the other measures in addition to a primary effect on force of infection, asexual parasite growth, or on pyrogenic thresholds. Analyses of trial datasets should aim to identify contributions of an intervention to each of these dimensions of protective efficacy.
The primary outcome of most trials is likely to use a single parasite density cutoff that is chosen to give a high specificity in order to reduce underestimation of efficacy since decisions to develop a vaccine depend on the magnitude of protection. However a highly specific case definition does not necessarily result in optimization of study power (Figure 1i) and in early stages of vaccine development it might be most important to to test whether there is any effect at all so a lower cutoff would be more appropriate. There is no reason why a threshold chosen to reduce bias in efficacy should be particularly appropriate for any other purpose and in particular a diagnostic threshold optimized for use in a trial is not necessarily appropriate as a tool in clinical management [12].
Analyses of trial data using such parasitaemia cutoffs have generally not quantified the bias that remains. The true efficacy is defined as E = 1  I_{ V }/I_{ P }where I_{ V }is the case incidence in vaccinees and I_{ P }is the case incidence in the placebo arm, and the usual estimate of efficacy is $\tilde{E}$ = 1  I_{ V }(x)/I_{ P }(x) where I_{ V }(x) is the incidence of cases at or above cutoff (x) in the vaccine arm, and I_{ P }(x) the corresponding incidence in the placebo arm. Assuming the specificity (ψ) of the diagnostic cutoff to be the same in both arms then an estimate of E adjusted for the effects of the imperfect case definition is:
where λ is the attributable fraction in the placebo arm. A potential improvement in efficacy estimates is to thus to estimate $\widehat{E}$ from $\tilde{E}$, λ, and ψ and to use $\widehat{E}$ as an estimate of E. If ψ is sufficiently close to unity, then the difference between these two estimates is small.
Exploratory analyses of the behaviour of $\widehat{E}$ suggest that it can be sensitive to x (Aponte, pers. comm), though it should not be so if the assumptions underlying its estimation are correct. The nonlinear logistic regression model most widely used for defining the parasitaemia cutoff[1] assumes a specific parametric form for the relationship between relative risk and parasite density. This can lead to severely biased estimates of the specificity of the cutoff if the relationship happens not to conform to this pattern [10]. This assumption is avoided in the latent class models that we have used in this paper which fit nonparametric curves for this relationship.
It is evident from Figure 1h though that the bias in efficacy does not only arise from lack of specificity in cutoff, and need not always be in the direction of underestimating efficacy. Bias also arises because of the specificity of cutoffs can differ between vaccine and placebo. Our model indicates that this is particularly a problem for asexual blood stage vaccines (vaccine B) (Figure 1g). This leads to the idea that perhaps different cutoffs should be used for vaccine and for placebo groups[7]. To justify this in practice though, it would be necessary to demonstrate a statistically significant difference between trial arms in the specificity vs cutoff relationship. This would be a difficult statistical exercise, (because the specificity is estimated only indirectly), and would lead to considerable difficulties in describing the results convincingly especially if the efficacy proved highly sensitive to the choices of cutoffs in the different groups. Because sample size is determined in order to give adequate power to measure the primary outcome (effect on case incidence), most trials are too small to conclusively demonstrate whether the specificity vs cutoff relationship varies between arms. The decision of the WHO Study Group on Measures of Malaria Vaccine Efficacy not to recommend trialarm specific cutoffs [9] is therefore probably wellfounded.
The most satisfactory alternative to using a single cutoff would probably be to estimate the total number of clinical malaria cases in each arm of the trial by assigning a probability to each fever case, rather than classifying each case dichotomously as above, or below, cutoff. This approach has not so far been used in analyses of clinical trials though it has been proposed as an alternative to the arbitrary choice of a cutoff [12]. It has been used in observational epidemiological studies [13, 14]. The preferred estimation method is to use a Bayesian latent class model to estimate the probabilities [11] carrying out this analysis separately for both placebo and vaccine arms. The simulations of this approach presented in Table 2 suggest that it has comparable power to that of the cutoff method. Moreover, interval estimates for all the quantities involved are readily available using software written in Winbugs [15] available from the author.
Such secondary analyses using latent class models, or considering the whole range of possible parasite density cutoffs will also help to identify possible biases in efficacy estimates made using single case definitions, at the same time as analysing the kind of protection. Where multicentre trials give heterogeneous efficacy estimates, it will be important to examine whether this can be accounted for by differential bias in the primary outcome measurements.
In a real trial the reduction in proportion of individuals infected varies over the trial period, depending on the time course of incidence, patterns of treatment with antimalarial drugs, and on the variation between individuals in exposure to vectors and responses to vaccination. These factors significantly complicate the analysis of relationships between infection and morbidity because, strictly speaking, the comparison should always be between contemporaneous data. This problem is particularly acute if parasites are cleared at the start of the trial, leading to complicated dynamics of infection and disease during the trial followup period. The present simulations do not address the implications of subpatent parasitaemia. This especially complicates analysis of effects of asexual blood stage vaccines because reduction of parasite densities differentially inflates the proportion of falsenegative blood slides in vaccinated individuals.
It follows that the analyses illustrated in this paper represent considerable simplifications of those that might be carried out in a real trial, where these complicating factors need to be taken into account. Nevertheless, when feasible, it would be logical to carry out secondary analyses corresponding to the different panels in Figure 1. Such analyses would help to clarify whether the effect of an intervention corresponds to that anticipated on the basis of the parasite stage that is targeted, and would highlight whether the primary measure of efficacy results from any unexpected behaviour in the parasitological and clinical data used to estimate it.
References
 1.
Smith T, Schellenberg JA, Hayes R: Attributable fraction estimates and case definitions for malaria in endemic areas. Stat Med. 1994, 13: 23452358. 10.1002/sim.4780132206.
 2.
Schellenberg JR, Smith T, Alonso PL, Hayes R: What is clinical malaria? Finding case definitions for field research in highly endemic areas. Parasitol Today. 1994, 10: 439442. 10.1016/01694758(94)901791.
 3.
Alonzo PL, Smith T, Schellenberg JR, Masanja H, Mwankusye S, Urassa H, Bastos DA, Chongela J, Kobero S, Menendez C, Teuscher T, Tanner M: Randomised trial of SPf66 vaccine against Plasmodium falciparum malaria in children in southern Tanzania. Med Trop (Mars ). 1995, 55: 4146.
 4.
Genton B, Betuela I, Felger I, AlYaman F, Anders R, Saul A, Rare L, Baisor M, Lorry K, Brown G, Pye D, Irving DO, Smith TA, Beck HP, Alpers MP: A recombinant bloodstage malaria vaccine reduces Plasmodium falciparum density and exerts selective pressure on parasite populations in a phase 12b trial in Papua New Guinea. J Infect Dis. 2002, 185: 820827. 10.1086/339342.
 5.
Alonso PL, Sacarlal J, Aponte J, Leach A, Macete E, Milman J, Mandomando I, Spiessens B, Guinovart C, Espasa M, Bassat Q, Aide P, OforiAnyinam O, Navia MM, Corachan S, Ceuppens M, Dubois MC, Demoitie MA, Dubovsky F, Menendez C, Tornieporth N, Ballou WR, Thompson R, Cohen J: Efficacy of the RTS,S/AS02A vaccine against Plasmodium falciparum infection and disease in young African children: randomised controlled trial. Lancet. 2004, 364: 14111420. 10.1016/S01406736(04)172231.
 6.
Smith T, Genton B, Baea K, Gibson N, Taime J, Narara A, AlYaman F, Beck HP, Hii J, Alpers MP: Relationships between Plasmodium falciparum infection and morbidity in a highly endemic area. Parasitology. 1994, 109: 539549.
 7.
Mwangi TW, Ross A, Snow RW, Marsh K: Case definitions of clinical malaria under different transmission conditions in Kilifi district, Kenya. J Infect Dis. 2005, 191: 19321939. 10.1086/430006.
 8.
Whitworth J, Morgan D, Quigley M, Smith A, Mayanja B, Eotu H, Omoding N, Okongo M, Malamba S, Ojwiya A: Effect of HIV1 and increasing immunosuppression on malaria parasitaemia and clinical episodes in adults in rural Uganda: a cohort study. Lancet. 2000, 356: 10511056. 10.1016/S01406736(00)027276.
 9.
Moorthy VS, Reed Z, Smith PG: Measurement of malaria vaccine efficacy in Phase III trials, report of a WHO consultation. Vaccine. 2007, in press:
 10.
Smith T, Vounatsou P: Logistic regression and latent class models for estimating positivities in diagnostic assays with poor resolution. Communications In StatisticsTheory And Methods. 1997, 26: 16771700.
 11.
Vounatsou P, Smith T, Smith AFM: Bayesian analysis of twocomponent mixture distributions applied to estimating malaria attributable fractions. J Roy Stat Soc C Applied Statistics. 1998, 575587. 10.1111/14679876.00129.
 12.
Rogers WO, Atuguba F, Oduro AR, Hodgson A, Koram KA: Clinical case definitions and malaria vaccine efficacy. J Infect Dis. 2006, 193: 467473. 10.1086/499314.
 13.
Smith T, Hurt N, Teuscher T, Tanner M: Is fever a good sign for clinical malaria in surveys of endemic communities?. Am J Trop Med Hyg. 1995, 52: 306310.
 14.
Vounatsou P, Smith T, Kitua A, Alonso PL, Tanner M: Apparent tolerance of Plasmodium falciparum in infants in a highly endemic area. Parasitology. 2000, 120: 19. 10.1017/S0031182099005211.
 15.
Spiegelhalter DJ, Thomas A, Best N, Lunn D: Winbugs Version 1.4. 2003, Cambridge, England, MRCBSU
Acknowledgements
The author acknowledges helpful discussions with Amanda Ross, and financial support from the Bill & Melinda Gates Foundation (project #39777). The author has no conflicts of interest concerning the work reported in this paper.
Author information
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
About this article
Received
Accepted
Published
DOI
Keywords
 Malaria
 Parasite Density
 Clinical Malaria
 Latent Class Model
 Malaria Vaccine
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate. Please note that comments may be removed without notice if they are flagged by another user or do not comply with our community guidelines.