|
|
||||||||
| ABSTRACT |
|
|
|---|
| INTRODUCTION |
|
|
|---|
Artemisinin combination therapy (ACT) has been promoted as an effective strategy to combat the emergence and spread of resistance.3 Currently, 38 countries have adopted ACT as first-line regimens for the treatment of malaria.4 However, a recent in vivo efficacy study on Thailands border with Cambodia reported an artesunate-mefloquine 28-day cure rate of 78.6%.5 In addition, an in vivo study on the Cambodian side of the border documented a 28-day cure rate of 85.7% in 2002.6 These studies suggest declining susceptibility to the artesunate-mefloquine combination in this subregion. However, it is not yet known whether these reported ACT failures are caused by true parasite resistance.
Because antimalarial resistance can emerge quickly, surveillance is a key element to a successful malaria control program.7,8 While the gold standard for monitoring drug efficacy is the in vivo trial, molecular markers could provide an inexpensive complementary tool. Before they can be used in a particular area, the relationship between specific molecular markers and in vivo resistance must be established.
One possible molecular marker is Plasmodium falciparum multidrug resistance gene-1 (pfmdr1). Single nucleotide polymorphisms (SNPs) and gene copy number have previously been associated with in vitro resistance to a wide array of antimalarials, including chloroquine, lumefantrine, quinine, halofantrine, mefloquine, and artesunate.916 However, in vitro and in vivo resistance are not always correlated.13 The purpose of this study was to determine whether pfmdr1 genetic changes are associated with in vivo evidence of resistance to artesunate-mefloquine in Cambodia.
| MATERIALS AND METHODS |
|
|
|---|
100,000 asexual parasites/µL, having a measured axillary temperature
37.5°C, providing informed consent (by parent or guardian, when appropriate), and willing to return for follow-up. In total, 25 children and 56 adults were enrolled. At enrollment, a brief clinical exam was performed, and a questionnaire was administered. Approximately 2 mL of venous blood was collected and frozen. Subjects were treated according to the national policy: they were given 25 mg/kg mefloquine and 12 mg/kg artesunate over 3 days (Days 02). The mefloquine was divided into two equal doses and given in the morning and evening on the day of enrollment (Day 0). The daily 4-mg/kg dose of artesunate was divided into two doses on Day 0 and given in one dose on the next 2 days. Subjects stayed at the clinic for 34 days (depending on when the parasites cleared) and had follow-up visits at 7, 14, 21, 28, 35, and 42 days after the date of enrollment. People who did not return on their own were actively sought out. At these follow-up visits, temperature was measured, and thick blood smears were used to screen for peripheral parasitemia. In instances of reoccurrence of parasitemia, an additional blood sample was acquired and stored on 3M Whatman filter paper (Whatman, Middlesex, UK).
Genotyping. DNA was extracted using the Qiamp DNA mini kit (Qiagen, Hilden, Germany). Samples taken at enrollment (Day 0), Day 1, and recurrence of parasitemia were used to distinguish between recrudescence and reinfection. The number of variants in three polymorphic genes (msp1, msp2, and glurp) was measured using size fractionization and ethidium bromide staining.18,19 The number of unique bands was determined visually and confirmed by measuring band size using Quantity One software on Gel Doc 2000 (Bio-Rad, Hercules, CA). If the recurrence specimen contained one or more of the variants seen in the enrollment or Day 1 specimen, the infection was classified as a recrudescence. If the recurrence specimen contained no variants present in the enrollment or Day 1 sample, it was classified as a reinfection.
Enrollment and recurrence samples were genotyped at pfmdr1 codons 86, 184, 1,034, and 1,042 using real-time polymerase chain reaction (PCR). A previously published protocol was followed, except that the concentration of probes was changed to 200 nmol/L.20
Pfmdr1 copy number was assessed for all samples using a modification of two previous assays.11,21 The primers and a FAM-TAMRA (6-carboxyfluorescein 6-carboxy-tetra-methylrhodamine) probe specific to a conserved region of pfmdr1 and the primers and a VIC-TAMRA (chemical structure not released by Applied Biosystems) probe specific to ß-tubulin were multiplexed so that both genes could be assayed in the same well. The primers were synthesized by MWG Biotech (High Point, NC), and the probes were synthesized by Applied Biosystems (ABI, Foster City, CA).
PCR reactions were performed on either the ABI Prism 7000 or ABI Prism 7300. Each reaction consisted of pfmdr1 probe (150 nmol/L), pfmdr1 primers (300 nmol/L), ß-tubulin probe (100 nmol/L), ß-tubulin primers (100 nmol/L), Abgene QPCR Mastermix (1x), DNA (2 µL), and water up to 25 µL. If amplification of a sample was unsuccessful (no amplification, replicates being > 50% different, or a copy number < 0.6), it was repeated with 4 µL DNA. The reaction conditions consisted of 95°C for 15 minutes and then 50 cycles of 95°C for 15 seconds and 60°C for 1 minute. The cycle threshold (CT) was calculated with SDS software (ABI). DNA from strains 3D7 and Dd2 were included on each plate. The 3D7 DNA was extracted from a laboratory grown culture and Dd2 DNA was acquired from MR4 (MRA-387; ATCC, Manassas, VA).
Pfmdr1 copy number was calculated according to the following formula: copy number = (Eßtubulin)CT(ßtubulin)/(Epfmdr1)CT(pfmdr1). The efficiency (E) of ß-tubulin, which was higher than that of pfmdr1, was assumed to be 2. Pfmdr1s efficiency, relative to that of ß-tubulin, was calculated for each plate by assuming the 3D7 control has one pfmdr1 copy. Dd2, an additional control, was previously determined to have approximately four pfmdr1 copies.22,23 In our assay, Dd2 had a mean copy number of 4.27 and a coefficient of variation of 5.77%.
Statistical analysis. The clinical data were entered into Excel (Microsoft, Redman, WA) and analyzed with Stata 8.2 (StataCorp., College Station, TX). The clinical outcome of each patient was classified according to the WHO protocol.17 Because there were no early treatment failures, the main outcome for this analysis was the recurrence of P. falciparum parasites in the peripheral blood. Recurrences of parasitemia were further categorized as reinfections or recrudescence by the msp1, msp2, and glurp genotyping.
The copy number of samples taken before treatment and those taken at recurrence were compared using the Wilcoxon signed rank test.
The relationship between molecular changes in pfmdr1 and treatment failure was estimated using survival analysis. Because the subjects were assessed at distinct time-points, the failure time is interval-censored. Discrete Cox proportional hazards model was used to take this censoring in account. This model was estimated using regression with a complementary log-log link.24 Indicator variables representing four time periods (021, 2228, 2935, and 3642 days) were included, whereas a constant term was not estimated. All variables were assessed for the proportional hazard assumption. The main outcome was time to recrudescence. Subjects who were lost to follow-up were censored on their last visit. Subjects who developed a P. vivax infection or a P. falciparum reinfection were censored on the day of that second diagnosis. The precision of all hazard ratio estimates were evaluated by calculating the confidence limits ratio (CLR).25
The genotype at codons 86, 184, 1,034, and 1,042 in pfmdr1 were combined into haplotypes, as previously described.11 The comparison between haplotypes I and III was made because these haplotypes had the largest sample size and because a previous in vitro study in Cambodia found haplotype III was associated with mefloquine resistance and haplotype I with mefloquine sensitivity.9 Samples that were a mixture of haplotypes I and III were coded as haplotype III. The relationship between copy number and time to recrudescence was not linear on the log hazard scale. Therefore, copy number was coded as a binary variable of less than three copies and at least three copies to reflect the relationship between copy number and time to recrudescence and to meet the proportional hazards assumption of the model.
The following covariates were evaluated as potential effect measure modifiers and confounders: previous use of antimalarials, initial parasite density, age, sex, and hematocrit. Initial parasite density was coded as a binary variable, with the cutoff point of 80,000 parasites/µL. Effect measure modification was assessed with the Wald test for the interaction term in the model (P
0.1 was considered significant). All non-effect measure modifiers were evaluated in a directed acyclic diagram to select a sufficient set of covariates to control for confounding.26
The population attributable rate fraction (AFp) and 95% confidence intervals (CIs) were calculated.27 The specific formula for AFp was [(HRadj 1) x p)/HRadj, where HRadj = adjusted hazard ratio, and p = proportion of cases that were exposed.
The relationship between molecular change in pfmdr1 and parasite clearance time was assessed using discrete Cox proportional hazards model in the same manner as described above.
Ethics. Informed consent was obtained from all adult subjects and from the parent or guardian of all minors. The pfmdr1 genotyping and data analysis were approved by the Institutional Review Boards at the United States Naval Medical Research Unit No. 2 and at UNC Chapel Hill School of Public Health. Both the genotyping and the efficacy study was approved by the National Ethics Committee for Health Research, Ministry of Health, Cambodia.
| RESULTS |
|
|
|---|
|
|
|
Comparison of enrollment and recurrence pfmdr1 genotypes. pfmdr1 copy number was determined for all recrudescent (N = 12) and all reinfection (N = 12) samples; 75.0% (9/12) of the reinfection samples and 66.7% (8/12) of the recrudescent samples had a higher copy number than the corresponding enrollment sample. These differences were not statistically significant (Wilcoxon signed rank test, for recrudescent: z = 1.26, P = 0.209; for reinfection: z = 1.57, P = 0.117). However, when comparing all enrollment and recurrent samples, the difference was statistically significant (Wilcoxon signed rank test, z = 2.06, P = 0.040). Therefore, artesunate-mefloquine treatment seems to select for increased pfmdr1 copy number.
Four recrudescent and 10 reinfection samples were genotyped for pfmdr1 polymorphisms. When comparing the haplotypes at enrollment and recurrence for the reinfection samples, five were identical, one sample changed from haplotype IV to III, two samples changed from being mixed (haplotypes I and III) to haplotype III, and two samples changed from having haplotype III to being mixed (haplotypes I and III). The four recrudescent samples had identical haplotypes to the enrollment sample (three of these paired samples had haplotype III and one had haplotype I). Thus, there seems to be no selection of a particular pfmdr1 haplotype because of artesunate-mefloquine treatment.
pfmdr1 and time to recrudescence.
Subjects with increased pfmdr1 copy number recrudesced sooner and more often than subjects with low copy number (Figure 2
; Table 3
). The crude hazard ratio for increased copy number (at least three copies compared with less than three copies) on time to recrudescence was 8.30 (95% CI: 2.6026.43). Adjusting for confounders resulted in a hazard ratio of 7.91 (95% CI: 2.3826.29). These results were not dependent on the statistical model, because the hazard ratio from a continuous Cox proportional model was similar in size and precision. Thus, pfmdr1 copy number is strongly related to time to recrudescence.
|
|
To investigate the impact of increased pfmdr1 copies on recrudescence in the entire cohort, the population attributable fraction was calculated. Using the prevalence of greater than three pfmdr1 copies in the cases (58.3%) and the adjusted hazard ratio (7.91), the population attributable fraction of recrudescence for increased pfmdr1 copy number was 50.9% (95% CI: 4.174.9%). pfmdr1 copy number seems to contribute to the cause of one half the cases of recrudescence observed in this study.
When the outcome of all recurrences of parasites (recrudescence + reinfection) was used, the estimate of effect for copy number was lower (unadjusted HR = 4.51, 95% CI: 1.9910.23; adjusted HR = 4.02, 95% CI: 1.739.34). In contrast, the hazard ratio for time to recurrence for haplotype III was higher but very imprecise (unadjusted HR = 4.78, 95% CI: 0.6435.64; adjusted HR = 4.40, 95% CI: 0.5933.02). Therefore, the estimated effect of pfmdr1 haplotype but not copy number changes depending on whether the treatment failures were PCR corrected. This suggests that if some misclassification of recrudescence and reinfection had occurred, the interpretation of the effects would not be substantially different.
The pfmdr1 copy number is a potential screening tool for drug resistance. Having parasites with greater than three pfmdr1 copies before treatment identified cases of recrudescence with a sensitivity and specificity of 58.3% and 88.2%, respectively, and a positive predictive value of 66.7%. Including information on either pfmdr1 haplotype or initial parasite density increased the specificity but decreased the sensitivity of predicting recrudescence (at least three copies + haplotype III: sensitivity = 50.0%, specificity = 91.2%; at least three copies + parasitemia
80,000: sensitivity = 50.0%, specificity = 91.2%).
pfmdr1 and parasite clearance time. Of the 80 subjects, 4 (5.0%) cleared the parasitemia by Day 1, 41 (51.3%) cleared by Day 2, 27 (33.8%) cleared by Day 3, and 8 (10.0%) cleared by Day 4. The parasite clearance times observed in this study are longer than reported previously for artesunate-mefloquine.28
Increased copy number was not associated with delayed parasite clearance when controlling for parasite density and hematocrit (Table 4
). There was a slight association between pfmdr1 haplotype and delayed parasite clearance time, although the estimate was imprecise.
|
| DISCUSSION |
|
|
|---|
Our observed association between pfmdr1 and failure is consistent with previous in vivo studies that looked at mefloquine monotherapy and artesunate-mefloquine combination therapy.13,21 In both studies, associations between pfmdr1 copy number and drug failure were found, but no association was seen between pfmdr1 haplotype and drug failure. In contrast, numerous in vitro studies have found associations between both copy number and haplotype and in vitro IC50.13,21,2931 There are four possible reasons for this. First, the effect of these mutations might be too small to be seen in vivo. Second, determinations of SNPs in complex polyclonal infections could result in misclassification. Third, in Pailin, the pfmdr1 genotype had minimal variation, with some sensitive genotypes being either rare or absent, so the effect might have been missed. Last, the factors responsible for resistance might differ by geographic region. Therefore, the inconsistency of results could be caused by the fact that several of the in vitro studies were executed outside of Southeast Asia.
It has been suggested that increased pfmdr1 copy number might cause resistance to both artesunate and mefloquine.1113,16 The observed drug failures are probably caused by mefloquine resistance, because a study in Pailin in 20012002 found the clinical isolates had high mefloquine IC50s and low artesunate IC50s.32 Also, parasite clearance is predominantly caused by artesunate, and no effect of pfmdr1 copy number on parasite clearance time was observed in this study. Therefore, the effect pfmdr1 copy number on recrudescence is most likely caused by the effect of pfmdr1 on mefloquine sensitivity and not artesunate sensitivity.21
We observed an increase in pfmdr1 copy number in parasites isolated after artesunate-mefloquine treatment. Because mefloquine has a much longer half-life than artesunate (1441 days versus 1 hour),33,34 resistant parasites were probably selected by exposure to subtherapeutic concentrations of the single drug mefloquine.35 Selection for mefloquine resistance in Pailin may have also been facilitated by the pre-existence of resistance to this drug.36,37 This observation suggests that ACT selects for resistance to the non-artemisinin partner drug.
Strengths of this study include the fact that genotyping was done directly on peripheral blood samples, which negates the potential bias caused by culturing.20 Another strength is its 42-day follow-up with PCR correction.38 The main limitation of this study is the small sample size, which affected the power and the precision of the effect estimates, especially for pfmdr1 haplotypes. Another limitation is that not all predictors of recrudescence were measured, such as pharmacokinetics, host genetics, immunity, and nutritional status.
Molecular markers of drug-resistant malaria, such as pfmdr1, are promising tools for the surveillance of drug resistance. pfmdr1 copy number is not only strongly associated with recrudescence to artesunate-mefloquine but also with Coartem failure (artemether-lumefantrine).39 Surveillance for increased pfmdr1 copy number could aid malaria control efforts by pinpointing areas where these drugs may be failing.
Received September 12, 2006. Accepted for publication January 1, 2007.
Acknowledgments: We thank the study subjects for their participation. We appreciate the assistance and support provided by Dr. Socheat Duong, Augustina Ika Susanti, Jesse Kwiek, William Miller, Anne Purfield, Paul Wilson, Nareth Kong, Wini Kania, and Dr. J. Kevin Baird. We thank MR4 for providing us with malaria DNA.
Financial support: This work was funded by the United States Navy Medical (NHRC/NMRC) ILIR Program.
Disclaimer: The views and opinions are those of the authors and do not purport to represent those of the US Navy or Department of Defense.
* Address correspondence to Steven R. Meshnick, Department of Epidemiology, UNC Chapel Hill School of Public Health, CB # 7435, Chapel Hill, NC 27599. E-mail: meshnick{at}email.unc.edu ![]()
Authors addresses: Alisa P. Alker, Naman K. Shah, and Steven R. Meshnick, Department of Epidemiology, UNC Chapel Hill School of Public Health, CB # 7435, Chapel Hill, NC 27599, E-mails: alker{at}unc.edu, naman{at}email.unc.edu, and meshnick{at}email.unc.edu. Pharath Lim, Rithy Sem, and Frederic Ariey, Institut Pasteur du Cambodge, 5 Monivong Boulevard, BP 983, Phnom Penh, Cambodia, E-mails: pharath{at}pasteur-kh.org, srithy{at}pasteur-kh.org, and fariey{at}pasteur-kh.org. Poravuth Yi, National Center for Parasitology, Entomology and Malaria Control (CNM), 372 Boulevard Monivong/Rue 322, Phnom Penh, Cambodia, E-mail: yiporavuth{at}yahoo.com. Denis Mey Bouth, World Health Organization, PO Box 1217, Phnom Penh, Cambodia, E-mail: denism{at}cam.wpro.who.int. Reiko Tsuyuoka, WHO Office in the Lao Peoples Democratic Republic, PO Box 343, Vientiane, Laos, E-mail: Tsuyuokar{at}lao.wpro.who.int. Jason D. Maguire, United States Naval Medical Research Unit No.2, U.S. Embassy, Jakarta, Indonesia, E-mail: JDMaguire{at}mar.med.navy.mil. Thierry Fandeur, Unité dImmunologie Moléculaire des Parasites CNRS URA 2581, Institut Pasteur, 28 rue du Dr ROUX, 75724 Paris Cedex 15, France, E-mail: tfandeur{at}pasteur.fr. Chansuda Wongsrichanalai, National Institute of Public Health/Naval Medical Research Unit No. 2 Laboratory, Phnom Penh, Cambodia, E-mail: chansuda{at}namru2.med.navy.mil.
Reprint requests: Steven R. Meshnick, Department of Epidemiology, UNC Chapel Hill School of Public Health, CB # 7435, Chapel Hill, NC 27599. E-mail: meshnick{at}email.unc.edu.
| REFERENCES |
|
|
|---|
This article has been cited by other articles:
![]() |
K. Nagamune, S. N. J. Moreno, and L. D. Sibley Artemisinin-Resistant Mutants of Toxoplasma gondii Have Altered Calcium Homeostasis Antimicrob. Agents Chemother., November 1, 2007; 51(11): 3816 - 3823. [Abstract] [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |