Previous Article | Next Article ![]()
Antimicrobial Agents and Chemotherapy, December 2008, p. 4233-4240, Vol. 52, No. 12
0066-4804/08/$08.00+0 doi:10.1128/AAC.00578-08
Copyright © 2008, American Society for Microbiology. All Rights Reserved.

EA4170, Parasitology and Tropical Medicine, Lyon 1 University, Lyon,1 Plate-forme Génomique, Institut Pasteur de Paris, Paris, France,4 Epidemiology Unit,2 Malaria Research Unit, Institut Pasteur de Madagascar, Antananarivo, Madagascar3
Received 2 May 2008/ Returned for modification 30 August 2008/ Accepted 14 September 2008
|
|
|---|
|
|
|---|
P. falciparum developed resistance to chloroquine (CQ) in various disease foci in the 1950s (47), but the first case of P. vivax resistance to CQ was not reported until 1989, in Papua New Guinea (37). Further sporadic cases were subsequently observed in the Western Pacific (Indonesia, Philippines) (1, 2, 12, 28, 35, 42, 44), Southeast Asia (Myanmar, India, Vietnam) (11, 15, 22, 24, 29, 31, 40), South America (Colombia, Guyana) (32, 41), and the Middle East (Turkey) (21). Despite these reports, it remains difficult to estimate the worldwide prevalence of P. vivax resistance to CQ. Genotyping by using molecular markers is required for in vivo studies of the clinical efficacy of antimalarial drugs, so that the recurrence of parasitemia due to reinfection can be distinguished from recrudescence. However, P. vivax genotyping in the context of antimalarial drug treatment is confounded by the occurrence of relapses due to hypnozoite stages (which can be the same isolate as that at baseline or a different one) or recrudescence (8, 9, 19). In vitro assays, which should provide drug susceptibility data free from the effects of confounding factors, such as host immunity, are still difficult to conduct (45). The molecular mechanisms underlying CQ resistance in P. vivax malaria remain unclear and may involve multigenic loci, but two genes orthologous to the pfmdr1 and pfcrt genes, which encode putative transporters (pvmdr1 [5, 43] and pvcrt-o [39]), have been identified as possible genetic markers of CQ resistance.
Malaria is endemic over three-quarters of Madagascar, with almost 1 million clinical cases reported per year, making this disease a major public health problem (20, 46). Four of the five malaria parasites known to infect human beings are present. Plasmodium vivax, the second most prevalent species after P. falciparum (causing 6.3% of all malaria cases), has a heterogeneous geographical distribution, ranging from 0% in the east coast area to 17.5% to the west of the foothills of the Central Highlands (4). CQ therapeutic efficacy has been assessed in Madagascar for P. falciparum (26) and Plasmodium malariae (3) but not for P. vivax malaria. CQ was replaced by artemisinin-based combination therapy (artesunate plus amodiaquine) as the recommended first-line treatment for uncomplicated malaria in 2005, but CQ is nonetheless still widely used. It is still recommended by the National Malaria Control Programme for the home management of presumed malaria in children (36).
The aims of this study were (i) to assess the therapeutic efficacy of CQ in P. vivax malaria and (ii) to explore the prevalences of mutations in the pvcrt-o and pvmdr1 genes in samples taken before treatment and the association between the presence of a mutant parasite and the clinical response of the patient to CQ treatment.
|
|
|---|
All patients with fever or a history of fever in the 48 h before their arrival at the health center were screened with a rapid diagnostic test, based on the detection of Plasmodium-specific lactate dehydrogenase (OptiMAL-IT; DiaMed AG, Cressier sur Morat, Switzerland). Giemsa-stained thin and thick blood films were prepared for each patient who was positive by the rapid diagnostic test. Plasmodium was identified to species level, and parasite density was assessed by a skilled microscopist. Blood from a finger prick was collected on filter paper once informed consent had been obtained from each adult subject and from at least one of the parents of each minor.
Clinical efficacy of CQ in P. vivax malaria. (i) Patient enrollment, treatment, and follow-up procedures. The clinical trial was carried out in 2006, at previously identified sentinel sites for assessing the therapeutic efficacy of CQ in P. vivax malaria. Patients with microscopically confirmed P. vivax infections were enrolled according to the inclusion criteria of the 2001 WHO guidelines (48). All patients included were weighed and their medical histories (including previous antimalarial medication) recorded. Clinical examination, including axillary temperature recording and Giemsa staining of serial thick and thin blood films, was carried out on days 0, 1, 2, 3, 7, 14, 21, and 28. Thick and thin films were read by a skilled microscopist. Patients were treated with the standard CQ regimen (25 mg/kg of body weight per day for 3 days) and monitored for 28 days. Patients were directly observed for 30 min after treatment, and the dose was readministered if vomiting occurred. Patients who repeatedly vomited their first dose of study medication were excluded from the study. Blood was blotted onto filter paper during follow-up and was stored at 4°C for DNA analysis. The hemoglobin concentration was determined on days 0 and 28 (HemoCue, Ängelholm, Sweden).
(ii) Outcome measures.
Treatment outcomes were assessed according to WHO 2001 guidelines (48) as "treatment failure" (TF; clinical deterioration due to P. vivax illness requiring hospitalization, with parasitemia and an axillary temperature of
37.5°C any time between days 3 and 28, or parasitemia on any day between days 7 and 28, regardless of clinical condition) or "adequate clinical and parasitological response" (ACPR; absence of parasitemia on day 28 without the criteria for TF having been met previously).
Patients with TF were treated with artesunate (4 mg/kg on days 0, 1, and 2) plus amodiaquine (10 mg/kg on days 0, 1, and 2), but their response to rescue therapy was not assessed.
Patients were excluded after enrollment if any of the following occurred: (i) use of antimalarial drugs outside of the study protocol, (ii) detection during follow-up of mixed malarial infections, except at day 28, (iii) parasitemia in the presence of a concomitant febrile illness potentially interfering with the classification of treatment outcome, (iv) withdrawal of consent, (v) loss to follow-up, (vi) protocol violation, or (vii) death due to a nonmalaria illness.
(iii) Laboratory procedures. Blood smears were stained by incubation with 10% Giemsa stain for 10 min. Asexual- and sexual-stage parasite counts were determined per 200 white blood cells. Parasitemia was calculated based on a white blood cell count of 8,000/µl. Thick-film examinations were considered negative if no parasite was found in 100 high-power fields.
DNA was extracted from blood spots on filter paper with Instagene Matrix resin (Bio-Rad, Marnes la Coquette, France), according to the manufacturer's instructions. The identification of the parasite species was confirmed by real-time PCR using species-specific primers, as described by de Monbrison et al. (10), with a protocol adapted for the RotorGene 3000 thermocycler (Corbett Life Science, Sydney, Australia).
Molecular genotyping techniques were used to distinguish recrudescence from new infection or relapse of liver stages for all patients with TF occurring after day 7. Briefly, blood samples collected on the day of enrollment, on day 1, and on the day of TF were analyzed for polymorphisms in the genes encoding circumsporozoite protein (pvcsp) (18) and merozoite surface protein-3 (pvmsp3) (6) and in six different microsatellite markers, as previously described (4). The PCR products of the genes encoding surface proteins were sequenced, whereas microsatellite PCR products were genotyped on the basis of size, using a GeneScan 500 LIZ size standard on an ABI Prism 3730 XL DNA analyzer. Genotype patterns on the day of TF were compared with those at treatment initiation and on day 1. An outcome was defined as recrudescence if all alleles present at the time of failure were present at the time of treatment initiation or on day 1; otherwise, it was defined as a reinfection or relapse.
Sequencing of pvcrt-o and pvmdr1. The pvcrt-o (exons 1 to 6) and pvmdr1 genes were amplified by nested PCR using gene-specific primers (Table 1). For the outer pvcrt-o PCR, amplification took place in the following reaction mixture: 2 µl of 10x buffer, 2.5 mM MgCl2, 0.4 mM each deoxynucleoside triphosphate, 0.25 µM each primer, 1 U of FirePol Taq polymerase (Solis Biodyne), and 2 µl of genomic DNA. PCR was performed under the following conditions: heating at 94°C for 5 min, followed by 40 cycles of heating at 94°C for 30 s, 56°C for 30 s, and 72°C for 2 min, and a final extension period at 72°C for 10 min. For the inner pvcrt-o PCR, the following reaction mixture was used: 5.5 µl of 10x buffer, 2.5 mM MgCl2, 0.4 mM each deoxynucleoside triphosphate, 0.27 µM each primer, 2.5 U of FirePol Taq polymerase (Solis Biodyne), and 3 µl of the amplicon. PCR was performed under the following conditions: 94°C for 5 min, followed by 45 cycles of 94°C for 30 s, 60°C for 30 s, and 72°C for 2 min, and a final extension period at 72°C for 10 min.
|
View this table: [in a new window] |
TABLE 1. Sequences of primers used to amplify the pvcrt-o (exons 1 to 6) and pvmdr1 genes
|
PCR products were purified by 96-well plate filtration (Millipore, St. Quentin en Yvelines, France) using a polyacrylamide gel (Bio-Gel P-100; Bio-Rad, Marnes-la-Coquette, France). Sequencing reactions were carried out with the ABI Prism BigDye Terminator cycle sequencing ready reaction kit run on a 3730 XL genetic analyzer (Applied Biosystems, Courtaboeuf, France). Electropherograms were visualized and analyzed with Seqscape software (version 2.0; Applied Biosystems, Courtaboeuf, France). Amino acid sequences were compared with wild-type sequences (GenBank accession no. AF314649 for pvcrt-o and GenBank accession no. AY571984.1 for pvmdr1) using BioEdit Sequence Alignment Editor software (16). Parasites with mixed alleles (both wild-type and mutant alleles present) were considered mutant.
Statistical analysis.
Data were entered and verified using EpiInfo software (version 6.04; Centers for Disease Control and Prevention, Atlanta, GA) and were analyzed with MedCalc software (version 9.1.0.1; MedCalc Software, Mariakerke, Belgium). The frequency of TF at the sentinel sites was unknown. We therefore included all patients with microscopically confirmed P. vivax malaria during the duration of the in vivo study, between February and June 2006 (26). Efficacy data were assessed by a per-protocol analysis including all patients who completed the study. An age-stratified analysis for patients <5 years old, 5 to 14 years old, and
15 years old was planned. Categorical variables were compared using chi-square or Fisher's exact tests, and continuous variables were compared using an independent-samples ANOVA (analysis of variance) or Mann-Whitney test.
The primary end point was the therapeutic response to CQ, based on parasitological and clinical cure by day 28 and adjusted by genotyping. The secondary end points included fever clearance, parasite clearance, and a change in hemoglobin levels between day zero and the last day of follow-up. Hypotheses were tested by determining the difference in risk, exact 95% confidence intervals (95% CI), and P values. A P value (two-tailed) of <0.05 was considered statistically significant.
Ethical approval. The study protocol was reviewed and approved by the Ethics Committee of the Ministry of Health of Madagascar (007/SANPF/2007). Informed written consent was provided by all patients or their parents/guardians before inclusion in the study.
Nucleotide sequence accession numbers. The exact sequence of each mutant allele has been submitted to GenBank (accession no. EU683813 to EU683819).
|
|
|---|
|
View this table: [in a new window] |
TABLE 2. Baseline characteristics of the enrolled patients by site, Madagascar, 2006
|
|
View this table: [in a new window] |
TABLE 3. Baseline characteristics of the enrolled patients by age group, Madagascar, 2006
|
Results for response to CQ treatment are presented by study site in Table 4 and by age group in Table 5. Eighty-eight patients (89.8%) successfully cleared P. vivax parasitemia after CQ treatment, whereas 10 patients (10.2%; 95% CI, 5.0% to 18.0%) experienced recurrences: 1 (10%) on day 7, 3 (30%) on day 21 (2 with mixed P. falciparum-P. vivax infections), and 6 (60%) on day 28 (1 with fever). After PCR correction, five cases of recurrence were confirmed as TF, giving a failure rate of 5.1% (95% CI, 1.7% to 11.5%). The frequency of TF was significantly higher in the Central Highlands (14.8%; 95% CI, 4.2% to 33.7%) than in the South (0%; 95% CI, 0% to 45.9%) and West (1.5%; 95% CI, 0% to 8.3%) (P = 0.02), regardless of the age group considered.
|
View this table: [in a new window] |
TABLE 4. Response to CQ treatment by site, Madagascar, 2006
|
|
View this table: [in a new window] |
TABLE 5. Response to CQ treatment by age group, Madagascar, 2006
|
![]() View larger version (8K): [in a new window] |
FIG. 1. Cumulative incidence of treatment failure over the 28-day follow-up period, Madagascar, 2006.
|
Analysis of the pvcrt-o gene. The pvcrt-o gene was successfully sequenced in 95 of the 105 P. vivax isolates collected (90.5%). All isolates were wild type, including those from patients displaying recrudescence. Only synonymous mutations or insertions (in introns) were found.
Analysis of the pvmdr1 gene. The pvmdr1 gene was successfully sequenced in 80 of the 105 isolates collected (76.2%), and at least four nonsynonymous mutations were detected in each case.
Ten nonsynonymous mutations were identified: five mutations described in previous studies (G698S, M908L, T958M, Y976F, and F1076L) and five new mutations (F194Y, S510T, S513R, I636T, and A829V) (Table 6). Six mutations were highly frequent and were found in almost all isolates tested (from 96.3% to 100%): the five previously identified mutations and the S513R mutation. The new mutations (other than S513R) were found only in isolates from the West and had a low prevalence (from 1.3% to 7.5%). Seven different haplotypes were identified, with four to seven amino acid substitutions (Table 7). The most common haplotype was the 513R 698S 908L 958M 976F 1076L sextuple-mutant type (85.0%). This haplotype was the only haplotype found at the South and Central Highlands sites. All seven haplotypes were found at sites in western Madagascar.
|
View this table: [in a new window] |
TABLE 6. Frequency distribution of mutations in the pvmdr1 gene by sampling site, Madagascar, 2006a
|
|
View this table: [in a new window] |
TABLE 7. Frequency distribution of pvmdr1 haplotypes by sampling site and treatment outcome
|
|
View this table: [in a new window] |
TABLE 8. Demographic, clinical, and parasitological characteristics of baseline isolates from patients with and without recrudescence
|
|
|
|---|
This study was the first to investigate the susceptibility of P. vivax to CQ in Madagascar. We assessed the efficacy of CQ at the national level by carrying out our clinical trial at six sentinel sites, selected as representative of the various epidemiological strata for malaria. Treatment outcome was assessed after adjustment based on molecular genotyping techniques to distinguish between recurrences of parasitemia due to new infection and those due to recrudescence or relapse. Since primaquine is not available in Madagascar, this drug was not used to treat hypnozoite stages in this study. The occurrence of relapse can be due to the same isolate as that observed at baseline or a different one. It can lead to misclassification of the recurrent parasitemia, as recrudescence (if the antimalarial drug level is under the MIC) in the first case or as a new infection in the second case. The interpretation of genotyping in this context is uncertain, and it is recommended to measure the concentration of the antimalarial drug in the blood (25). The absence of measurements of blood CQ concentrations in patients displaying recrudescence is the main limitation of this study. Overall, our data suggest that CQ-resistant P. vivax isolates are present in Madagascar, particularly in the foothills of the Central Highlands. We have previously found a high frequency of sulfadoxine-pyrimethamine TF in patients with P. vivax malaria in the same area (4).
In vivo studies remain the "gold standard" method for assessing antimalarial drug efficacy. However, due to the difficulties involved in determining TF rates, particularly for P. vivax malaria, molecular markers seem to be useful for monitoring the drug resistance of malaria parasites. The K76T mutation in the P. falciparum crt gene is known to be involved in CQ resistance, but very few data are available concerning the possible relationship between the P. vivax crt-o and mdr1 genes and CQ resistance. Only a few studies have been carried out, and associations between TF and nonsynonymous mutations in isolates obtained before treatment have not yet been clearly established (5, 30, 38). However, the Y976F substitution in the pvmdr1 gene is thought to be involved in CQ resistance in P. vivax (23, 43). Suwanarusk et al. found that the geometric mean 50% inhibitory concentration of CQ was significantly higher in P. vivax isolates carrying the Y976F mutation than in isolates with the wild-type allele (43).
We used a sequencing approach to screen for the presence of mutations in the target genes in order to confirm the involvement of these two genetic markers, pvcrt-o and pvmdr1, in CQ resistance in P. vivax. The sequences analyzed in pvcrt-o were those of exons 1 to 6, in which mutations have already been described (30). All the Malagasy P. vivax isolates sequenced were wild type, with only single-nucleotide polymorphisms or insertions in intron sequences. However, isolates were more polymorphic for the pvmdr1 gene, with 11 synonymous mutations and 10 nonsynonymous mutations, 5 of which had not been identified previously. Six of these mutations, including Y976F, were highly frequent (from 96.3% to 100%). These findings may reflect the CQ resistance situation of P. falciparum, in which the pfmdr1 mutation N86Y is highly prevalent even in the absence of the K76T pfcrt mutation (34). Comparison with the SalI reference strain showed all samples to be of the mutant type, with haplotype patterns comprising four to seven mutations. We also found that the 510R and 513T amino acid substitutions were mutually exclusive, suggesting that two different stepwise mutant selection processes had occurred. For the 510R amino acid substitution, the 510T 908L 958M 1076L haplotype must have been selected first and the 510T 698S 908L 958M 976F 1076L haplotype second, whereas for the 513T amino acid substitution, the 513R 908L 958M 976F 1076L haplotype must have been selected first, followed by the 513R 698S 908L 958M 976F 1076L haplotype and, finally, the 513R 698S 829V 908L 958M 976F 1076L or 513R 636T 698S 908L 958M 976F 1076L haplotype.
Unlike the investigators in another recent study (43), we found no significant correlation between the pvmdr1 Y976F mutation and CQ treatment outcome, essentially because this mutation was already fixed in the Malagasy P. vivax population. There are two possible reasons for this. First, the Y976F mutation may have been selected by CQ drug pressure over the last 50 years in Madagascar. This possibility implies a role of this mutation in CQ resistance. Second, the Malagasy P. vivax population may have been introduced into Madagascar at the beginning of the Christian era from Indonesia (13), where the prevalence of the Y976F mutation is high (43). If that is the case, this mutation is unlikely to play any role in CQ resistance.
The monitoring of the resistance of malaria parasites to CQ through molecular markers remains complex. The usefulness of marker sets may differ for different areas and needs to be evaluated for a given area before use. These findings call into question the validity of in vitro tests. Major efforts have been made to develop a continuous system for the maintenance of P. vivax isolates in culture (45), but this remains difficult. Furthermore, the determination of antimalarial drug activity based on in vitro tests and the schizont maturation method is more complicated for P. vivax than for P. falciparum, due to the asynchrony of P. vivax blood parasites. Furthermore, the parasite response to antimalarial dugs may depend on the parasite development stage (43).
Our study provides new data concerning the susceptibility of P. vivax to CQ in Madagascar and suggests that CQ-resistant isolates have been selected in Madagascar, particularly in the foothills of the Central Highlands. However, levels of resistance to CQ remain much lower in P. vivax than in P. falciparum (5.1% versus 44%) (26). The 976Y pvmdr1 mutation was found not to be a useful tool for monitoring CQ resistance in Madagascar, and further efforts are required to develop surveillance tools for the temporal and geographical mapping of emerging drug resistance in P. vivax malaria.
This study was supported by grants from Natixis/Impact Malaria through the Observatoire de la Résistance aux Antipaludiques Project and the Genomics Platform, Pasteur Génopôle, Pasteur Institute, France. Sample collection was funded by the Global Fund project for Madagascar, round 3 (Community Action to Roll Back Malaria grant MDG-304-G05-M). C. Barnadas is a Ph.D. student supported by the Fondation Jeunesse Internationale (Fondation de France), BioMérieux (Prix BioMérieux Infectiologie 2006), Association des Internes et Anciens Internes en Pharmacie des Hôpitaux de Lyon (Prix R. Rizard), and Hospices Civils de Lyon.
Published ahead of print on 22 September 2008. ![]()
|
|
|---|
locus of Plasmodium vivax: global and local diversity. Am. J. Trop. Med. Hyg. 61:518-525.[Abstract]
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»