- Research
- Open Access
Detectability of Plasmodium falciparum clones
- Michael T Bretscher^{1, 2}Email author,
- Francesca Valsangiacomo^{1, 2},
- Seth Owusu-Agyei^{3},
- Melissa A Penny^{1, 2},
- Ingrid Felger^{1, 2} and
- Tom Smith^{1, 2}
https://doi.org/10.1186/1475-2875-9-234
© Bretscher et al; licensee BioMed Central Ltd. 2010
- Received: 24 May 2010
- Accepted: 18 August 2010
- Published: 18 August 2010
Abstract
Background
In areas of high transmission people often harbour multiple clones of Plasmodium falciparum, but even PCR-based diagnostic methods can only detect a fraction (the detectability, q) of all clones present in a host. Accurate measurements of detectability are desirable since it affects estimates of multiplicity of infection, prevalence, and frequency of breakthrough infections in clinical drug trials. Detectability can be estimated by typing repeated samples from the same host but it has been unclear what should be the time interval between the samples and how the data should be analysed.
Methods
A longitudinal molecular study was conducted in the Kassena-Nankana district in northern Ghana. From each of the 80 participants, four finger prick samples were collected over a period of 8 days, and tested for presence of different Merozoite Surface Protein (msp) 2 genotypes. Implications for estimating q were derived from these data by comparing the fit of statistical models of serial dependence and over-dispersion.
Results
The distribution of the frequencies of detection for msp2 genotypes was close to binomial if the time span between consecutive blood samples was at least 7 days. For shorter intervals the probabilities of detection were positively correlated, i.e. the shorter the interval between two blood collections, the more likely the diagnostic results matched for a particular genotype. Estimates of q were rather insensitive to the statistical model fitted.
Conclusions
A simple algorithm based on analysing blood samples collected 7 days apart is justified for generating robust estimates of detectability. The finding of positive correlation of detection probabilities for short time intervals argues against imperfect detection being directly linked to the 48-hour periodicity of P. falciparum. The results suggest that the detectability of a given parasite clone changes over time, at an unknown rate, but fast enough to regard blood samples taken one week apart as statistically independent.
Keywords
- Parasite Density
- Deviance Information Criterion
- Merozoite Surface Protein
- Survey Round
- Imperfect Detection
Background
In areas of high endemicity of Plasmodium falciparum, human hosts are often superinfected with multiple clones of the parasite [1]. Identification of these concurrent infections is important for understanding patterns of drug resistance [2] and of the transmission of the parasite. PCR-based methods for detecting parasites not only have lower detection limits than blood smear microscopy, but also make it possible to distinguish genetically distinct clones, and hence to compute multiplicity of infection. But at least two diagnostic problems remain: i) the same host might be infected with more than one parasite clone of the same genotype, which can introduce bias into estimates of multiplicity of infection [3]. ii) PCR detection can be negative because the sample taken does not contain any parasites. This may happen due to effects of acquired immunity or synchronization of the parasite population. Failure to account for imperfect detection biases several standard epidemiological measures, such as prevalence and multiplicity of infection. Most critically, analysis of drug failure rates using molecular typing may overlook breakthrough parasite clones or conversely misclassify them as new infections after treatment. Repeated blood samples from the same host can be analysed to estimate the probability that a clone is detected in any given sample (the detectability, q). For microscopy data, where individual parasite clones cannot be distinguished, the statistical methods of [4, 5] are applicable. Both assume infections are neither acquired nor cleared during the study. For molecular data, several pieces of work aiming at estimating infection duration and force of infection also yielded measurements of detectability and its dependence on age of the host [1, 6–8]. These methods make use of data collected over longer time periods (several months up to a year), with surveys every 1 to 2 months, the kind of which may not be easily available in practice. Moreover, the obtained estimates of detectability depend on simultaneous estimates of infection and recovery rates as well as on assumptions concerning these processes. A simple method is therefore presented, to estimate the detectability of infecting clones from molecular data with short inter-survey intervals. It makes use of pairs of surveys sufficiently close in time, such that reinfection with the same parasite genotype can be safely excluded. The method is similar to the one presented in [4], but adapted for the context of molecular diagnostic methods. This implies that the maximal number of "infections" is not limited by the number of hosts in the study, but rather represents individual parasite clones. The methods of [1, 6–8] as well as the one presented here assume that the detections of an infecting clone at different time points are independent from each other. While it seems reasonable to make such an assumption, provided intervals between surveys are long enough, it is not clear how long these intervals need to be. Numerous publications report complicated periodic behaviour of fevers or parasitaemia [9], or detection events [10], which creates a need to establish the circumstances under which the methods mentioned above can be applied.
In order to evaluate the effect of possible "nonrandom" behaviour of clonal infections on estimates of detectability, a longitudinal study comprising 80 individuals was conducted in northern Ghana. From each participant, four blood samples were collected over a period of 8 days. Using these data, various statistical models are compared with respect to their goodness of fit, and a series of hypothesis tests is performed. The resulting statistical description of the within-host dynamics of P. falciparum clones, as observed by molecular typing methods, allows us to justify a simple algorithm for obtaining reasonably robust estimates of q and specify the circumstances under which this method is applicable.
Methods
Study site and sample collection
Genotyping
DNA was eluted from "ISOCode™Stix" filter paper and screened for presence of P. falciparum by polymerase chain reaction (PCR). Sample processing and PCR conditions have been described in detail [14]. In brief, all samples were subjected to PCR using primers specific for the merozoite surface protein (msp) 2 locus. Genotypes were distinguished on the basis of length polymorphism and PCR fragments were precisely sized by automated capillary electrophoresis and GeneMapper^{®} software. An inhouse generated software identified all genotypes per sample and transformed the data into different formats suitable for data management and statistical analysis. Given the high number of msp2 genotypes in the population, re-infection with the same genotype was assumed to be a rare event. As a consequence of this, for any given host, msp2 genotype is assumed to be synonymous with "infecting clone" in all analyses.
Data analysis
Only data of those participants who were present at all four survey rounds, and where at least one genotype was found, were included in the analysis. This reduced the number of individuals in the data set to 69. Patterns of appearance and disappearance of specific parasite genotypes depend on rates of infection and clearance as well as on detectability. However, for the purpose of the present analysis, acquisition and loss of infections were neglected. It was assumed that there are no false positive results and that an infecting clone is present throughout all four surveys if detected at least once. This is justified by the comparatively short time interval between the first and the last survey, and by previously published estimates of infection and clearance rates from the dataset of the main study [8]:
Data coding
R1 | R2 | R3 | R4 | sequence no. | count |
---|---|---|---|---|---|
0 | 0 | 0 | 0 | 0 | - |
0 | 0 | 0 | 1 | 1 | 43 |
0 | 0 | 1 | 0 | 2 | 42 |
0 | 0 | 1 | 1 | 3 | 26 |
0 | 1 | 0 | 0 | 4 | 54 |
0 | 1 | 0 | 1 | 5 | 10 |
0 | 1 | 1 | 0 | 6 | 13 |
0 | 1 | 1 | 1 | 7 | 19 |
1 | 0 | 0 | 0 | 8 | 28 |
1 | 0 | 0 | 1 | 9 | 21 |
1 | 0 | 1 | 0 | 10 | 41 |
1 | 0 | 1 | 1 | 11 | 34 |
1 | 1 | 0 | 0 | 12 | 64 |
1 | 1 | 0 | 1 | 13 | 22 |
1 | 1 | 1 | 0 | 14 | 41 |
1 | 1 | 1 | 1 | 15 | 61 |
A series of χ^{2} tests and Spearman's rank correlation analysis yielded qualitative information on thetemporal behavior of detectability. Further, a series of models for the dynamics of detectability were fit to the data using Bayesian MCMC. These are de-scribed in detail below. The models and their estimates of detectability are compared using Deviance Information Criterion (DIC) as measure of goodness of fit [15]. That only sequences with at least one positive result were included in the data, and therefore the data are biased, was accounted for in all analyses. The software Winbugs [16] was used for all Bayesian model fitting, whereas for all other analyses the software package R was used [17].
Models of detection
In addition, a simple method of directly measuring detectability from pairs of surveys (M0) is used. Applying this method to all available survey pairs in the data set and comparing the estimates of q with the model results allows us to develop criteria for the circumstances under which the method may be used.
M0: Direct estimation of detectability
The binomial likelihood model implies statistical independence of detections at different time points. This assumption might be violated if detectability exhibits temporally structured behavior. The detectability model underlying this method is identical to M1 below, but this method only uses two observations, and model fitting in order to estimate q is done analytically.
This simple method is compared to models M1 to M3 (below) in order to justify it's use, and to establish the conditions where it can be applied. For a formal derivation of (1) and confidence intervals for $\stackrel{\wedge}{q}$ please refer to the appendix (additional file 1).
M1: Binomial model
This implies independence of detecting a clone at time t from whether it was detected at other time points, and homogeneity of the infection population with respect to detectability.
M2: Beta-binomial model
where a and b are the shape parameters of the beta distribution.
M3: First order Markov Chain
An important feature of this simple model is that it not only represents a random walk in "detection space" (i.e. switching between being detected and not being detected), but that it can also be interpreted as a random walk in "detectability space" (switching between the two detectabilites). This can be illustrated as follows: The probability that a clone changes its internal state from q_{0} to q_{1} is equal to the probability that it is detected while being in state q_{0}, which is equal to q_{0}. Likewise, the probability of a transition from q_{1} to q_{0} is 1 - q_{1}. The resulting transition matrix for such a process is identical to T. However, since this model is fitted to a population of sequences, as opposed to just fitting it to one single time-series, one has to be careful in interpreting a possible best fit of this model: simple heterogeneity in detectability among clonal infections would also result in different estimates for q_{1} and q_{0}, even if there is no random walk in detectability within a single clone. Therefore, only in combination with the results of M2, which is also able to capture such heterogeneity, is one able to interpret t_{10} and t_{01} as transition probabilities between q_{0} and q_{1}.
Bias correction of detectability estimates
with N_{ obs } = 519.
Results
In the complete study population (80 individuals), the average prevalence across all four survey rounds was 46% by microscopy and 0.69% by PCR. The dataset used for statistical analyses comprised 69 parasite-positive individuals between 6 months and 20 years of age, with a median age of 5.2 years (interquartile range 3.5-9.7). The median multiplicity of infection (MOI) among these was 10 (inter-quartile range 7-13), when pooling all four observations from each individual. This differs from standard practice when reporting MOI, but was justified given the very short interval between the surveys. The obtained value is expected to be a better estimate of the true MOI. When only considering single survey rounds, the median MOI among PCR-positives was 4. The 519 detected clones belonged to 77 different msp2 genotypes, with the most common allele reaching a frequency of 9.2%.
Tests of proportion and correlation
A series of hypothesis tests was performed in order to gain insight into the statistical properties of the data-generating process. These do not relate to the models M1-M3 directly, but rather aim to look at similar questions using a completely different methodology. Any conclusions would need to be consistent with both approaches. In this analysis, the detection bias is accounted for by adding a total of 26 all-zero sequences to the data set, such that the total number of sequences equals 545. This is the "true" number of clonal infections, as estimated robustly by models M1 to M3. In the following list, H_{1}-H_{5} indicate the hypotheses tested, and the corresponding p-values obtained using χ^{2} tests (with the exception of the Spearman's Rank Correlation analysis) are given:
H_{1}: All 4 surveys have an equal proportion of positive results, i.e. $\sum {d}_{i,1}=\sum {d}_{i,2}=\sum {d}_{i,3}=\sum {d}_{i,4}$ (Data: 312, 284, 277, 236). This hypothesis of stationarity is rejected by a χ^{2} test with 3 degrees of freedom: P-value < 0.0001.
H_{2}: The frequencies s_{ i } of the 16 binary sequences (including the added all-zero sequence) are multinomially distributed with expectations ${s}_{i}=545{q}^{{o}_{i}}{(1-q)}^{4-{0}_{i}}$, where o_{ i } is the number of positive testing results in sequence i). This hypothesis, effectively proposing that a Bernoulli-process is able to perfectly describe the data, is rejected by a χ^{2} test with 14 degrees of freedom: P-Value < 0.0001.
H_{3}: The number of sequences with i = 0, 1, .., 4 detections are multinomially distributed with expectations ${s}_{i}=545{(}_{{0}_{i}}^{4}){q}^{{o}_{i}}{(1-q)}^{4-{0}_{i}}.$ This is a slightly relaxed version of H_{2}, such that the time order of detections is neglected, and sequences with a certain number of detections are pooled. However, this hypothesis is rejected by a χ^{2} test with 3 degrees of freedom: P-Value < 0.0001.
H_{4}: The frequencies of all four possible results of a survey pair (i.e. "00","01","10" and "11") are multinomially distributed with expectations ${s}_{i}=545{q}^{{o}_{i}}{(1-q)}^{2-{0}_{i}}.$ This is a special case of H_{2}, only applied to a pair of surveys, instead of the whole dataset. Except for survey pair 2-3 (p-value 0.06) this hypothesis is rejected on all survey pairs by χ^{2} tests with 2 degrees of freedom: P-values < 0.0001.
Direct estimation of q on all survey pairs, using M0
survey pair | interval (days) | n _{1} | n _{2} | $\stackrel{\wedge}{q}$ | 95% CI | p-value |
---|---|---|---|---|---|---|
1-2 | 1 | 220 | 188 | 0.63 | 0.59-0.68 | 9.4e-12 |
3-4 | 1 | 233 | 140 | 0.55 | 0.49-0.60 | 9.1e-7 |
2-3 | 4 | 293 | 134 | 0.48 | 0.43-0.53 | 1.1e-2 |
1-3 | 5 | 235 | 177 | 0.60 | 0.55-0.65 | 2.4e-7 |
2-4 | 5 | 276 | 112 | 0.45 | 0.39-0.50 | 1.7e-3 |
1-4 | 7 | 272 | 138 | 0.50 | 0.45-0.55 | 0.83 |
A Spearman's Rank Correlation analysis showed significant positive correlation between all pairs of surveys which are 24 h apart, and no correlation for all other pairs, with the exception of the pair formed of surveys 1 and 3.
Model comparison
Comparison of Models M1-M3
model | age groups | $\overline{q}$ | DIC | ||
---|---|---|---|---|---|
M1 (binomial) | 1 | 0.50 | 2870.2 | ||
M1 (binomial) | 5 | 0.55, 0.52, 0.49, 0.53, 0.41 | 2872.3 | ||
model | Age groups | $\overline{q}$ | var(q) | DIC | |
M2 (beta-bin.) | 1 | 0.51 | 0.003 | 2871.5 | |
M2 (beta-bin.) | 5 | 0.51, 0.46, 0.46, 0.47, 0.46 | 0.004,0.005,0.005,0.005,0.005 | 2871.3 | |
Model | age groups | $\overline{q}$ | q _{0} | q _{1} | DIC |
M3 (Markov) | 1 | 0.50 | 0.47 | 0.59 | 2856.4 |
M3 (Markov) | 5 | 0.55, 0.53,0.49, 0.53, 0.41 | 0.45, 0.44, 0.47, 0.56, 0.46 | 0.65, 0.64, 0.56, 0.55, 0.45 | 2860.0 |
Estimates of q
The estimates of detectability (Table 3) were found to be similar for models M1 through M3, especially when common parameters for all age groups were estimated. All values were approximately 0.5, and showed little prior sensitivity. Although measuring separate detectabilities for every age group did not improve model fit, a decreasing trend of detectability with age was observed. Estimates of q for the youngest group are between 0.51 and 0.55, and decrease to values between 0.41 an 0.46 for the oldest age group. This is consistent with the findings of other authors [8]. Estimates of q obtained using M0, however, show some variation, with values ranging from 0.45 to 0.63. Table 2 shows the corresponding estimates obtained from all available pairs of surveys. Only for the measurement using survey pair 1-4 were the criteria for using the method fulfilled, as the corresponding p-value of 0.83 indicates, and the value of q estimated from this pair matches the estimates from the models very well.
Discussion
The short term dynamics of asymptomatic P. falciparum clonal infections in vivo were characterized in order to find a simple way of measuring detectability in the field. A series of statistical tests as well as a progression through three simple models provided insight into some statistical properties of within-host dynamics monitored by molecular typing. Classical PCR ignores absolute parasite densities, but length polymorphic amplicons make it possible to distinguish between co-infecting parasite clones. Detectability, however, can be used as a proxy for parasite densities, as the two must be correlated. Since key epidemiological measures such as prevalence and multiplicity of infection (MOI) depend on the numerical value of detectability, planning and monitoring of malaria interventions rely on accurate measurements of detectability. It is important to note that detectability may not only depend on host or parasite factors, but also on the methods for generation of genotyping data. These include the method for collecting blood samples, the actual volumes of blood collected, storage conditions, PCR conditions, and competition within the PCR assay limiting the detection of minority clones. The present analysis did not consider any of these factors, but rather assumed that their impact is more or less identical for all samples and clones.
Within-host dynamics
Measurement of detectability from longitudinal data will for practical reasons rely on binomial models of detection. This creates a need to establish under what conditions such binomial models are applicable. Different hypotheses tested on the typing data showed complicated dynamics of clonal infections for short timescales. These dynamics could not be described by a binomial model (rejection of hypotheses H_{1} to H_{3}). A hypothesis which could not always be rejected was H_{5}. This hypothesis stated that the number of successful detections in sample pairs were binomially distributed. It was not rejected for the sample pair collected at the most distant dates, i.e. from surveys 1 and 4 with interval of 7 days (Table 2). This finding could indicate that the processes governing detectability on short time scales are prone to stochastic variation such that the effect of the initial state of a clone vanishes after some time, and the two observations become independent. In other words, the corresponding test could be interpreted as test of an ergodic hypothesis, which implies that the system under investigation "forgets" it's initial state after enough time has passed. That the frequencies of "01" and "10" sequences are not equal, and therefore H_{4} (stricter than H_{5}) is rejected on all survey pairs except pair 2-3, questions this interpretation, and can not be explained it in a satisfactory way. Since there are consistently more "10" pairs than '01", one could presume that the detectability of clones simply decreases with time, which would also explain the decreasing trend in the number of detections per survey (H_{1}). However, this is mere speculation and can hardly be shown from this dataset. Nevertheless, the method of directly measuring detectability is presumably little affected by this phenomenon, as it's estimates of q only depend on the sum of single positive pairs and the obtained numerical values of q agree very well with the results of the other models.
Intuitively, one would expect a certain amount of variation in detectability among clonal infections, especially since these were pooled across individuals. It is therefore surprising that M2, which would allow for such variation, measured zero variance of q and effectively reduced to the binomial model M1. The best fitting model M3 offers a possible interpretation, as it is capable of capturing change in detectability over time. M3 models the time evolution of detectability as a Markov chain, which is equivalent to assume that a clone has detectability q_{0} if it was not detected on the preceding survey, and detectability q_{1} if it was. The obtained estimates of q_{0} and q_{1} as 0.47 and 0.59, respectively, could either indicate variation in the dataset with respect to detectability or that the detectability of a clone performs a random walk in detectability space, alternating between the two states q_{0} and q_{1}. Since M2 reduced to M1, and estimated practically no variation in detectability, we have to assume the latter.
One might expect parasite densities to fluctuate with a period of approximately 48 hours, as observed in malariatherapy-data [9], and in good agreement with in-vitro measurements of a 48 hour erythrocytic cycle. In fact, such periodic behaviour of asymptomatic infections has been reported [10, 19]. The present analysis does not find a 48-hour periodicity, rather the opposite: both the best fitting model as well as the results of the Spearman's rank correlation analysis indicate positive autocorrelation between time points which are 24 hours apart. A process with a periodicity of 48 hours, on the contrary, should show negative correlation. A possible explanation for the difference between malariatherapy data and the data presented here could be that malariatherapy patients were not immune and therefore had fever more often. The question of periodicity in symptomatic malaria should be considered separately, and it's causes are thought to be well explained [20, 21]: high temperature (fever) differentially affects the intra-erythrocytic stages of parasite development, and nearly stops development in some of these. This leads to "queuing" of the parasite population, and when the fever goes down, all parasites continue their development in a synchronized way. Fever can by definition not be operating in asymptomatic individuals, but at least in simian and avian malaria an effect of normal diurnal changes in body temperature on synchronization has been demonstrated, alongside with the observation that sometimes the parasite population is split into "two broods [..], coming to schizogony on alternate days" [22]. Two broods, synchronized within themselves, appearing in the peripheral blood with a 48 hour periodicity, yet with a 24 hour phase-shift, would appear in the data as having a 24 hour periodicity. This would be consistent with the finding that detection results one day apart are positively correlated. As the data does not contain smaller time intervals, however, any such periodicity cannot be distinguished from a simple gradual change in detectability.
Analysis of periodicity of clonal infections would ideally make use of long series of parasitological observations of untreated infections with short intervals, but few studies have collected such data, partly for ethical reasons. Exceptions include the malariatherapy datasets [23], the studies of Farnert et al[10, 24] and Magesa et al[25] in Tanzania, and Bruce et al[19, 26] from Papua New Guinea. Bruce et al aggregated data for paired observations with identical interval length and calculated the probability of detecting an infection at the second occasion, conditional on it being detected at the first occasion. This analysis suggests values of detectability similar to the estimates in the present study, with a six day periodicity. This periodicity was interpreted as signal of a 48-hour underlying cycle because the sampling interval was three days, which meant that six-day and two-day periodicity could not be distinguished. Similar analyses of the other available datasets would be of value.
Measurement of detectability
This approach was found to give very similar results as the more sophisticated methods, provided the underlying assumption is met: the number of successful detections must follow a binomial distribution. The statistical properties of the data, as assessed by a series of tests, suggest that if there is an interval of at least 7 days between consecutive surveys, it is safe to make these assumptions (see H_{5}). Alternative methods [1, 6–8] rely on the same two assumptions, yet are further incorporating models for the processes of acquisition and loss of infections. While those may themselves be of interest, the associated measurements of q may be affected by assumptions about acquisition and loss of infections. Direct estimation using M0 is therefore recommended as a simple and practical alternative, if only detectability is of interest, and if the interval between two surveys is short enough so acquisition and loss of infection clones can safely be excluded.
Epidemiological significance of detectability
This implies that the measurement error for prevalence, when neglecting detectability, should be highest at the lowest multiplicities of infection - a situation to be expected when approaching local elimination. It is therefore desirable to routinely report q together with other epidemiological measures, if possible.
In drug efficacy trials, the phenomenon of imperfect detection complicates the task of distinguishing new from breakthrough infections, and therefore must have an influence on drug efficacy estimates. In addition, residual drug levels may keep parasite densities at undetectable levels for some time, which is usually taken into account when designing drug efficacy trials. No satisfactory statistical methodology for analysis of such trials appears to exist, taking into account both imperfect detection and residual drug levels. It is suspected that many recrudescent infections, i.e. infection clones which survive treatment and are detected several days or weeks later, might in fact often be detected earlier if multiple testing took place. This is strongly supported by the findings of [27], who note that consecutive-day blood sampling changes the results of a drug efficacy trial compared to single-day blood sampling.
Conclusions
The presented work demonstrates the importance of paying attention to the phenomenon of imperfect detection not only in the sense of assessing the sensitivity of diagnostic tests, but also looking at it as a property of infections or individuals. Various epidemiological measures, such as prevalence or MOI, are affected by imperfect detection. Failure to account for it may severely distort the outcome of measurements and may even lead to wrong conclusions. As an example, the decrease in prevalence with age, as it is frequently observed in malaria, might actually mean that detectability decreases with age, while prevalence remains constant or even increases. Indeed, some publications suggest that this might be the case [8]. In addition, it is likely that underestimation of prevalence may be substantial in situations where the multiplicity of infection is low. This is because the chance of missing every single clone in a host is highest when there are only few.
A simple method of estimating detectability from molecular data, using pairs of surveys, was presented. It is a modification of existing methods, which can deal with data on multiple infections within one host. The numerical estimates of detectability obtained using said formula appeared remarkably robust. Through comparison of the detectability estimates with estimates from different models, and through a series of statistical tests, the conditions under which the underlying assumptions of the method are fulfilled could be established. Its use is recommended when the time interval between the two surveys is one week or more, but discouraged on data with shorter time intervals, if possible. Both the method itself as well as the way of addressing it's applicability are not restricted to malaria, but may in a similar way be used for other infectious diseases where molecular data on individual clones is available.
The restrictions on applicability stem from the complicated dynamics of detectability on short time scales. These were investigated and it was found that treating individual detections as statistically independent is only an acceptable approximation for time intervals longer than one week. Contrary to expectation, however, no changes in detectability indicative of a 48 hour cycle were found, as is reported from malariatherapy data. This suggests that not the 48 hour erythrocytic cycle of P. falciparum is dominating detectability in vivo, but that other factors, such as e.g. the dynamics of the immune system, may be important. As the participants of the study must be considered partly immune, it is presumed that the within-host dynamics of infections differ between immune and non-immune individuals. This questions the use of malariatherapy data for fitting of within-host models for the immune host, and encourages further collection of relevant data as well as development of analysis methods in order to gain better insight into the within-host dynamics of P. falciparum in immune individuals.
Declarations
Acknowledgements
This work was supported by Bill & Melinda Gates Grant no. 39777 and Swiss National Science Foundation, project no. 320030 - 125316. Sincere thanks are extended to the study participants and the staff of the Navrongo Health Research Centre involved in the field work.
Authors’ Affiliations
References
- Smith T, Felger I, Fraser-Hurt N, Beck H: Effect of insecticide-treated bed nets on the dynamics of multiple Plasmodium falciparum infections. Trans R Soc Trop Med Hyg. 1999, 93 (Suppl 1): 53-57. 10.1016/S0035-9203(99)90328-0.View ArticlePubMedGoogle Scholar
- Hastings IM, Smith TA: MalHaploFreq: a computer programme for estimating malaria haplotype frequencies from blood samples. Malar J. 2008, 7: 130-10.1186/1475-2875-7-130.PubMed CentralView ArticlePubMedGoogle Scholar
- Hill W, Babiker H: Estimation of numbers of malaria clones in blood samples. Proc Biol Sci. 1995, 262: 249-257. 10.1098/rspb.1995.0203.View ArticlePubMedGoogle Scholar
- Mullen K, Prost A: Decreased microfilarial load and its effect on the calculation of prevalence and the rate of false negatives in the detection of on- chocerciasis. Int J Epidemiol. 1983, 12: 102-104. 10.1093/ije/12.1.102.View ArticlePubMedGoogle Scholar
- Aron J: Malaria epidemiology and detectability. Trans R Soc Trop Med Hyg. 1982, 76: 595-601. 10.1016/0035-9203(82)90219-X.View ArticlePubMedGoogle Scholar
- Smith T, Vounatsou P: Estimation of infection and recovery rates for highly polymorphic parasites when detectability is imperfect, using hidden Markov models. Stat Med. 2003, 22: 1709-1724. 10.1002/sim.1274.View ArticlePubMedGoogle Scholar
- Sama W, Owusu-Agyei S, Felger I, Vounatsou P, Smith T: An immigration-death model to estimate the duration of malaria infection when detectability of the parasite is imperfect. Stat Med. 2005, 24: 3269-3288. 10.1002/sim.2189.View ArticlePubMedGoogle Scholar
- Sama W, Owusu-Agyei S, Felger I, Dietz K, Smith T: Age and seasonal variation in the transition rates and detectability of Plasmodium falciparum malaria. Parasitology. 2006, 132: 13-21. 10.1017/S0031182005008607.View ArticlePubMedGoogle Scholar
- Molineaux L, Diebner H, Eichner M, Collins W, Jeffery G, Dietz K: Plasmodium falciparum parasitaemia described by a new mathematical model. Parasitology. 2001, 122: 379-391. 10.1017/S0031182001007533.View ArticlePubMedGoogle Scholar
- Farnert A, Snounou G, Rooth I, Bjorkman A: Daily dynamics of Plasmodium falciparum subpopulations in asymptomatic children in a holoendemic area. Am J Trop Med Hyg. 1997, 56: 538-547.PubMedGoogle Scholar
- Owusu-Agyei S, Smith T, Beck H, Amenga-Etego L, Felger I: Molecular epidemiology of Plasmodium falciparum infections among asymptomatic inhabitants of a holoendemic malarious area in northern Ghana. Trop Med Int Health. 2002, 7: 421-428. 10.1046/j.1365-3156.2002.00881.x.View ArticlePubMedGoogle Scholar
- Falk N, Maire N, Sama W, Owusu-Agyei S, Smith T, Beck H, Felger I: Comparison of PCR-RFLP and Genescan-based genotyping for analyzing infection dynamics of Plasmodium falciparum. Am J Trop Med Hyg. 2006, 74: 944-950.PubMedGoogle Scholar
- Binka F, Morris S, Ross D, Arthur P, Aryeetey M: Patterns of malaria morbidity and mortality in children in northern Ghana. Trans R Soc Trop Med Hyg. 1994, 88: 381-385. 10.1016/0035-9203(94)90391-3.View ArticlePubMedGoogle Scholar
- Felger I, Irion A, Steiger S, Beck H: Genotypes of merozoite surface protein 2 of Plasmodium falciparum in Tanzania. Trans R Soc Trop Med Hyg. 1999, 93: 3-9. 10.1016/S0035-9203(99)90320-6.View ArticlePubMedGoogle Scholar
- Spiegelhalter DJ, Best NG, Carlin BP, van der Linde A: Bayesian measures of model complexity and fit. J R Stat Soc Ser B. 2002, 64: 583-639. 10.1111/1467-9868.00353.View ArticleGoogle Scholar
- Lunn DJ, Thomas A, Best N, Spiegelhalter D: WinBUGS A Bayesian modelling framework: Concepts, structure, and extensibility. Stat Comput. 2000, 10: 325-337. 10.1023/A:1008929526011.View ArticleGoogle Scholar
- Team RDC: R: A Language and Environment for Statistical Computing. 2008, Vienna, Austria, ISBN 3-900051-07-0, [http://www.R-project.org]Google Scholar
- Smith T, Beck H, Kitua A, Mwankusye S, Felger I, Fraser- Hurt N, Irion A, Alonso P, Teuscher T, Tanner M: Age dependence of the multiplicity of Plasmodium falciparum infections and of other malariological indices in an area of high endemicity. Trans R Soc Trop Med Hyg. 1999, 93: 15-20. 10.1016/S0035-9203(99)90322-X.View ArticlePubMedGoogle Scholar
- Bruce M, Donnelly C, Packer M, Lagog M, Gibson N, Narara A, Walliker D, Alpers M, Day K: Age- and species-specific duration of infection in asymptomatic malaria infections in Papua New Guinea. Parasitology. 2000, 121: 247-256. 10.1017/S0031182099006344.View ArticlePubMedGoogle Scholar
- Kwiatkowski D: Febrile temperatures can synchronize the growth of Plasmodium falciparum in vitro. J Exp Med. 1989, 169: 357-361. 10.1084/jem.169.1.357.View ArticlePubMedGoogle Scholar
- Kwiatkowski D, Nowak M: Periodic and chaotic host- parasite interactions in human malaria. Proc Natl Acad Sci USA. 1991, 88: 5111-5113. 10.1073/pnas.88.12.5111.PubMed CentralView ArticlePubMedGoogle Scholar
- Hawking F, Worms M, Gammage K: 24- and 48-hour cycles of malaria parasites in the blood; their purpose, production and control. Trans R Soc Trop Med Hyg. 1968, 62: 731-765. 10.1016/0035-9203(68)90001-1.View ArticlePubMedGoogle Scholar
- Collins WE, Jeffery GM: A retrospective examination of secondary sporozoite- and trophozoite- induced infections with Plasmodium falciparum: development of parasitologic and clinical immunity following secondary infection. Am J Trop Med Hyg. 1999, 61: 20-35.View ArticlePubMedGoogle Scholar
- Farnert A, Lebbad M, Faraja L, Rooth I: Extensive dynamics of Plasmodium falciparum densities, stages and genotyping profiles. Malar J. 2008, 7: 241-10.1186/1475-2875-7-241.PubMed CentralView ArticlePubMedGoogle Scholar
- Magesa SM, Mdira KY, Babiker HA, Alifrangis M, Färnert A, Simonsen PE, Bygbjerg IC, Walliker D, Jakobsen PH: Diversity of Plasmodium falciparum clones infecting children living in a holoendemic area in north-eastern Tanzania. Acta Tropica. 2002, 84: 83-92. 10.1016/S0001-706X(02)00179-1.View ArticlePubMedGoogle Scholar
- Bruce MC, Galinski MR, Barnwell JW, Donnelly CA, Walmsley M, Alpers MP, Walliker D, Day KP: Genetic diversity and dynamics of Plasmodium falciparum and P. vivax populations in multiply infected children with asymptomatic malaria infections in Papua New Guinea. Parasitology. 2000, 121: 257-272. 10.1017/S0031182099006356.View ArticlePubMedGoogle Scholar
- Martensson A, Ngasala B, Ursing J, Veiga MI, Wiklund L, Membi C, Montgomery SM, Premji Z, Farnert A, Bjorkman A: Influence of consecutive-day blood sampling on polymerase chain reaction-adjusted parasitological cure rates in an antimalarial-drug trial conducted in Tanzania. J Infect Dis. 2007, 195: 597-601. 10.1086/510910.View ArticlePubMedGoogle Scholar
Copyright
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.
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.