|
|
||||||||
| ABSTRACT |
|
|
|---|
3-fold. Methods presented are applicable to both observational studies and clinical trials assessing the effect of therapies on transmission potential. | INTRODUCTION |
|
|
|---|
Development of gametocytemia is evaluated by serial microscopic examination of thick blood smears after drug administration. Gametocytemia is usually measured on the day of treatment and 7, 14, 21, and 28 days after treatment, during which it usually peaks and then starts to decrease. Studies frequently evaluate the impact of drugs on gametocytemia after treatment by measuring the gametocyte carriage rate on these days.26 Typically, two descriptors of the gametocyte distribution have been used: 1) the proportion of individuals with detectable gametocytes at a particular day after treatment (percentage > 0 or discrete component), and 2) mean gametocytemias, excluding those with no gametocytes (mean for those that are > 0 or continuous component).
We propose the use of an alternative outcome for transmission potential that can incorporate magnitude and duration of gametocytemia: the area under the curve (AUC) of the gametocyte levels after treatment. The main advantage of our approach is that the outcome incorporates the longitudinal data of the gametocytes levels collected at several days. We examined this analytic approach using data from a study designed to evaluate determinants of treatment response to sulfadoxine-pyrimethamine (SP) among P. falciparum-infected patients in the Pacific coast region of Colombia.
In this region, resistance to SP has remained low in contrast to high rates of SP failure elsewhere in South America. However, we have reported that known parasite mutations that confer low-level resistance to SP, even though insufficient to cause treatment failure, may contribute to the potential for the transmission of P. falciparum and the spread of resistance.5 Specifically, point mutations in parasite dihydrofolate reductase (DHFR) that confer in vitro resistance to pyrimethamine were associated with longer parasite clearance time (PCT) and the presence of gametocytes after SP treatment.5 This finding suggested that even before clinical SP resistance is apparent drug treatment may be responsible, by increased gametocytemia, for selecting resistant parasites and contributing to the spread of resistance.
This report evaluates the use of the AUC of the gametocyte distribution after treatment as an alternative way to measure transmission potential. We describe methods for characterizing the joint effects of PCT and DHFR mutations at codons 108 and 51 on two outcomes: presence of gametocytes and magnitude of gametocytemia among treated individuals in whom gametocytes are present. We characterized the role of the explanatory variables (i.e., PCT and mutations) on the two types of outcomes, and this characterization can be helpful in assessing the ability of antimalarial treatment regimens to prevent the occurrence of gametocytes and/or to reduce the magnitude of gametocytemia.
| MATERIALS AND METHODS |
|
|
|---|
Study population. One hundred twenty subjects with acute P. falciparum malaria were enrolled in the study, received standard SP treatment, and were followed-up for 21 days or until treatment failure. The study protocol was reviewed and approved by the Universidad del Valle Ethics Committee in compliance with national regulations governing the protection of human subjects. All subjects provided written informed consent. Subjects were 170 years old (median = 21.2 years), 63% were male, and 88% lived in the urban area of the study. Four individuals had treatment failure, received rescue treatment, and were not included in the analysis of the distribution of gametocytes. In addition, parasite DNA from blood samples of 20 individuals could not be amplified by polymerase chain reaction, which prevented the assessment of mutations for DHFR and dihydropteroate synthase (DHPS). The 96 individuals with available data on presence of DHFR mutations and with parasites clearance time on days 1, 2, and 3 comprised the study population.
Outcome variable: AUC of gametocytemia after treatment. Gametocyte density was determined on days 0, 3, 7, 14, and 21 after treatment from thick blood films by counting the number of sexual parasites per 200 white blood cells (WBCs) (or per 500 if the count was less than 10 parasites/200 WBCs) and calculating parasites per microliter assuming a WBC count of 8,000/µL.
The AUC is a summary calculation used when serial measurements on each subject under study are carried out. Alternative outcomes include mean of all measurements, height of peak, time to reach peak and time to return to baseline level.8 The AUC is frequently used in clinical pharmacology, where estimates of the area made from drug serum levels over time can be interpreted as the total uptake of the drug. In the case of assessing the response to antimalarial treatment, we propose using AUC of the gametocyte distribution after treatment because it incorporates both the magnitude and the duration of transmission potential.
The AUC from days 3 to 21 was calculated as AUC = [(7 3) x (g3 + g7)/2 + (14 7) x (g7 + g14)/2 + (21 14) x (g14 + g21)/2]/(21 3) where gd represents gametocyte density on day d. We scaled this area by 18 (21 3) so that it represents AUC per day and transformed this by log10 to make its distribution to be amenable to methods based on the Gaussian bell shaped distribution. Specifically, the outcomes used in all analyses that follow were 1) a binary outcome indicating presence of gametocytes, and 2) the log10 (AUC) if AUC > 0.
The four values of gametocytes at days 3, 7, 14, and 21 were required to have the AUC as a complete outcome measure. For instances when a value was missing between two available values (e.g., missing on day 14 with values available on days 7 and 21), simple interpolation (i.e., mean) of the available data was used to complete the missing data. The other type of missing data was that of not having data on any subsequent days. Fifteen (16%) of the 96 individuals had missing data on day 21 and only 3 (3%) had missing data on days 14 and 21. We completed the data for each of these 18 individuals by using the observed data in individuals with similar values of gametocytemia at all the visits if the 18 individuals have data available. Thus, we used the standard and simple "hot deck" method for completing missing data.
To estimate the effect of these imputations we repeated the analysis using only observed data with heterogeneous total length of follow-up. In particular, for the 18 individuals with missing data the AUC corresponds to the previous formula until days 7 or 14.
Determinants of AUC. Analysis of determinants of AUC was based on two main exposures: PCT and presence of mutations in vitro associated with SP resistance at DHFR and DHPS codons. Polymerase chain reaction methods to assess parasite mutations were applied according to protocols described elsewhere9 and available from http://medschool.umaryland.edu/CVD/plowe.html.
Since we restricted our study population to the individuals who cleared parasites on days 1, 2, or 3, this exposure preceded the AUC from days 3 to 21, which was used as the outcome variable. The PCT (1, 2, or 3 days) and mutations (none, one [108], or 2 [108 and 51]) were treated as indicator variables in all regression models using no mutation and PCT = 1 as the reference categories.
Statistical analysis.
Given that the number of subjects who maintain 0 gametocytes over the post-treatment period (i.e., AUC = 0) was not uncommon, the true distribution of AUC had a mass at 0 with probability 1
and we modeled the
values greater than 0 as a normal distribution (in the log scale). To characterize differences in AUC according to exposures, we performed analysis in two stages. In the first stage, we used the full population and implemented logistic regression methods for the odds of gametocytes being present at any time to assess the predictors of AUC being greater than 0 (i.e., presence of gametocytes). In the second stage, we restricted to those with AUC > 0 and implemented standard linear regression methods for the log10 (AUC) to describe the heterogeneity of the magnitude of gametocytemia. Both models were run as univariate by either mutations or clearance time and also as multivariate with both mutations and clearance time in the models. For both regressions we used standard methods of maximum likelihood estimation and we estimated the 95% confidence intervals for their corresponding coefficients.
| RESULTS |
|
|
|---|
|
|
Table 2
shows the univariate and bivariate analyses of the effects of PCT and mutations on presence of gametocytes. In the bivariate model, although the effect of mutations suggest a dose response, it was not statistically significant. In contrast, PCT showed that the presence of gametocytes was significantly increased by PCT > 1 day, but was not affected by PCT increasing from two to three days.
|
|
| DISCUSSION |
|
|
|---|
To overcome potential limitations due to missing data, we applied simple methods of imputation of missing data based on the observed data to achieve complete data for each individual. Comparative analysis with available data showed that these imputations had not influence our results.
We have shown how AUC can be used to study the transmission potential of Plasmodium after treatment in an observational setting. Clinical trials that compare the efficacy of antimalarial drugs used alone or in combination could incorporate the estimation of AUC and help to define which treatments would potentially be better able to deter the rapid appearance of resistance.
In this study, we showed that a prolonged PCT is the primary explanatory variable of the presence of gametocytes. These results are consistent with a previous report by Price and others2 that showed a significant trend between longer time to clear parasitemia and gametocyte carriage rate. Our results suggest that a PCT of two or three days is associated with a similar substantial and highly significant odds ratio for presence of gametocytes. This suggests that rapidly acting treatments that clear parasites by the first day after treatment may effectively reduce transmission of resistance. In addition, for magnitude of gametocytemia we showed that both PCT and mutations are important, whereby a PCT prolonged to three days is significantly associated to a high magnitude, and that the presence of any mutations is equally deleterious towards an increase in the magnitude of gametocytemia.
Antimalarial drugs that inhibit gametocytogenesis or kill mature gametocytes have the potential to interrupt transmission of malaria. Given in combination, these therapies have been proposed as strategies to secure a longer life for available drugs.16 However, a comprehensive evaluation of the efficacy of the different drug combinations requires valid and comprehensive methods of analysis, as those we propose in this report.
A limitation of our study is that we did not assess the actual transmission of gametocytes detected in blood smears but just presence and magnitude. The transmission capacity of gametocytes should be confirmed when possible through direct-feeding assays or membrane-feeding assays. In a study in The Gambia, gametocytes emerging after successful treatment of P. falciparum-infected children with chloroquine were significantly more likely to be those with resistance genotypes and were infective to mosquitoes.10
Based on the shape of the dose response curve between gametocytemia levels and transmission to mosquitoes, it is expected that at similar conditions of malaria endemicity, higher levels of gametocytemias among treated individuals should increase the source of new infections. In particular, infectivity to mosquitoes was shown to present a log-sigmoid relationship with gametocyte density in original malaria therapy studies,15 but this may vary considerably according to several factors, including the level of transmission-blocking activity among the population, which is affected by previous malaria episodes17 and variable infectivity of circulating gametocytes.18
The data analysis methods presented provided means to jointly characterize the effect of PCT and mutations on the presence and magnitude of gametocytes over the period of 18 days after the third day after treatment with SP. These methods are applicable to both observational studies and clinical trials assessing the effect of therapies on transmission potential.
Received January 25, 2006. Accepted for publication May 31, 2006.
Financial support: This research was supported in part by a grant AI055687-02 from the National Institutes of Health, and the Universidad del Valle, Cali, Colombia.
* Address correspondence to Fabián Méndez, Escuela de Salud Pública, Universidad del Valle, San Fernando, Calle 4B No. 36-140, Edificio 118, Cali, Colombia. E-mail: famendez{at}univalle.edu.co ![]()
Authors addresses: Fabián Méndez, Escuela de Salud Pública, Universidad del Valle, San Fernando, Calle 4B No. 36-140, Edificio 118, Cali, Colombia, Fax: 57-2-557-0425, E-mail: famendez{at}univalle.edu.co. Álvaro Muñoz, Department of Epidemiology, The Johns Hopkins Bloomberg School of Public Health, 615 North Wolfe Street, E-7648, Baltimore, MD 21205, E-mail: jvaldez{at}jhsph.edu. Christopher V. Plowe, Center for Vaccine Development, University of Maryland School of Medicine, 685 West Baltimore Street, HSF 480, Baltimore, MD 21201, E-mail: cplowe{at}medicine.umaryland.edu.
| REFERENCES |
|
|
|---|
This article has been cited by other articles:
![]() |
J. T. Bousema, C. J. Drakeley, P. F. Mens, T. Arens, R. Houben, S. A. Omar, L. C. Gouagna, H. Schallig, and R. W. Sauerwein Increased Plasmodium falciparum Gametocyte Production in Mixed Infections with P. malariae Am J Trop Med Hyg, March 1, 2008; 78(3): 442 - 448. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |