Measuring the efficacy of anti-malarial drugs in vivo: quantitative PCR measurement of parasite clearance
© Beshir et al; licensee BioMed Central Ltd. 2010
Received: 15 July 2010
Accepted: 5 November 2010
Published: 5 November 2010
Artemisinin-based combination therapy, currently considered the therapy of choice for uncomplicated Plasmodium falciparum malaria in endemic countries, may be under threat from newly emerging parasite resistance to the artemisinin family of drugs. Studies in Southeast Asia suggest some patients exhibit an extended parasite clearance time in the three days immediately following treatment with artesunate monotherapy. This phenotype is likely to become a more important trial endpoint in studies of anti-malarial drug efficacy, but currently requires frequent, closely spaced blood sampling in hospitalized study participants, followed by quantitation of parasite density by microscopy.
A simple duplex quantitative PCR method was developed in which distinct fluorescent signals are generated from the human and parasite DNA components in each blood sample. The human amplification target in this assay is the β tubulin gene, and the parasite target is the unique methionine tRNA gene (pgmet), which exhibits perfect sequence identity in all six Plasmodium species that naturally infect humans. In a small series of malaria cases treated as hospital in-patients, the abundance of pgmet DNA was estimated relative to the human DNA target in daily peripheral blood samples, and parasite clearance times calculated.
The qPCR assay was reproducibly able to replicate parasite density estimates derived from microscopy, but provided additional data by quantification of parasite density 24 hours after the last positive blood film. Robust estimates of parasite clearance times were produced for a series of patients with clinical malaria.
Large studies, particularly in Africa where children represent a major proportion of treated cases, will require a simpler blood sample collection regime, and a method capable of high throughput. The duplex qPCR method tested may fulfil these criteria, and should now be evaluated in such field studies.
Artemisinin-based combination therapy (ACT) is currently considered the best treatment for uncomplicated Plasmodium falciparum malaria . The artemisinin component of these combinations produces a characteristic rapid reduction in parasite biomass immediately after treatment, but these compounds are metabolized in hours, and thus combination with a partner drug is required to provide complete parasite clearance with short treatment regimens, and to minimize the opportunity for evolution of parasites resistant to either component drug [2, 3]. Recently, however, there have been reports of a decline of ACT efficacy, manifested as slower parasite clearance than expected in a significant subset of patients receiving artemisinin monotherapy, and ACT, in Western Cambodia [4–6].
In areas where ACT is deployed as the major anti-malarial treatment option, early warning of any decrease in the rate of parasite clearance in treated patients would enable adequate countermeasures to be taken to prevent onward transmission of these potentially resistant parasites, and to facilitate the search for genetic markers associated with this newly recognized parasite phenotype. In an individual, the persistence of parasites to day 3 following chemotherapy is a reliable predictor of subsequent treatment failure , and thus a marker for this would assist in the rational use of alternative treatment regimens to prevent prolonged parasite carriage, and enhanced transmission of novel parasite genotypes .
Parasite Clearance Time (PCT) is estimated by measuring change in peripheral parasitaemia in a sequence of samples taken after treatment, typically for 72 hours, and currently this requires microscopic examination of blood films. However, accurate microscopic estimation of parasite density requires considerable expertise; false negatives can easily be read especially at low-level parasitaemia; it is labour intensive and suffers from low throughput. Further, as the majority of ACT-treated patients are free of parasitaemia by microscopy within 48 hours of treatment, many closely-spaced venous blood samples in catheterized in-patients have been taken in previous studies to obtain a useful estimate of PCT by this method [5, 6]. This approach is not suitable for large clinical trials, particularly in African settings, where the majority of participants are likely to be outpatients under the age of 10 years; such excessive sampling would be considered unethical in most such studies. A simple, rapid method for estimating PCT that did not require multiple venous blood samples, but could be achieved with relatively few finger-prick blood samples, spotted and dried on filter paper for storage, would thus be a useful new tool for large African studies of anti-malarial efficacy. However, absolute quantitation of parasite density from dried filter-paper blood spots is difficult, as the volume of blood in each spot is usually not known, and the efficiency of DNA extraction and target amplification differs between paper types and even among individual patients .
This paper reports the development and use of a new duplex quantitative real-time PCR (qPCR) assay to measure parasite clearance after anti-malarial treatment. To overcome problems with absolute parasite quantitation in the intended future application to DNA extracted from dried filter paper blood spots, parasite DNA abundance is normalized to a human white blood cell DNA target in the same sample, permitting calculation of relative parasite density in each sequential patient sample by the delta-delta CT (ΔΔCT) method . The time required to clear 95% of the parasites (PCT95), and the ratio in parasite density between admission and 48 hours post-treatment, or parasite reduction ratio (PRR48), are then easily derived. These parameters were used to compare treatment response among patients treated with different anti-malarial regimens at the Hospital for Tropical Diseases (HTD), London.
Samples and patients
Parasitized peripheral blood was obtained from malaria patients admitted to HTD, as part of routine in-patient care . Samples were collected on admission and approximately every 24 hours until cleared from the peripheral blood after treatment. Written informed consent was given by patients 1, 2 and 3 for additional biological studies of malaria parasites isolated from their blood samples to be carried out, according to a protocol approved by the University College London Hospital (UCLH) Research Ethics Committee (ref. 07/Q0505/60). Parasite isolates from patient 4 were investigated in detail at the specific request of the attending physician, as a suspected case of quinine treatment failure.
For each sample, DNA was purified from an aliquot of 200 μL of peripheral blood collected in EDTA, using the Qiagen QIAamp blood extraction kit and was eluted in a volume of approximately 100 μl . In addition, for patients 1, 2 and 3 between 2.5 mL and 4.0 mL of pre-treatment (day 0) leuco-depleted blood was prepared and DNA extracted using the Qiagen maxi extraction kit. P. falciparum genome sequence data was generated from these samples using the Solexa Illumina platform at the Wellcome Trust Sanger Institute (WTSI), Hinxton, Cambridgeshire. These data, for patients 1, 2 and 3, are freely available via the MapSeq software platform at WTSI, designated OX007, OX005 and OX006, respectively .
PCT using microscopy
Microscopic examinations of thick and thin smears, stained by the rapid Field's method for diagnosis and estimate of parasite density, and Giemsa-stained for discrimination of Plasmodium species, were performed according to standard methods by at least two experienced microscopists. This was performed in the Department of Clinical Parasitology at HTD, a fully accredited diagnostic laboratory that undergoes external quality assurance as a participant in the UK National External Quality Assessment Scheme for Parasitology. Parasite density is expressed as the percentage of erythrocytes which carry asexual stages of P. falciparum. The density of Plasmodium ovale curtisi parasites in patient 1 was not specifically measured by the microscopists, but was estimated to be less than 1% of total parasitaemia, P. falciparum being dominant at all time points.
Designing primers and probes
Specific primer pairs and double-labelled fluorescent hydrolysis probes were designed using Primer Express software (AP Applied Biosystems) to amplify part of the Plasmodium tRNA methionine (PgMET) gene, which exhibits 100% sequence identity between Plasmodium vivax and P. falciparum (Genbank accession numbers: AL844509, nt 1575551-1575622; AZ570363, nt 127-56). Primers based on these sequences successfully amplified a fragment of pgmet from P. ovale curtisi, Plasmodium ovale wallikeri, Plasmodium malariae and Plasmodium knowlesi DNA. The identity of the PCR target was confirmed by direct sequencing for each species. Thus it was concluded that this amplicon was suitable for cross-species parasite quantitation, under the assumption that each species had a single copy of the pgmet gene.
Amplification primers and double-labelled hydrolysis probes for PgMET and HumTuBB.
A 20 μl PCR master mix was prepared as follows: 5.5 mM MgCl2 buffer (Bioline), 0.3 mM dNTPs (Bioline), 0.3 μM forward and reverse primers, 0.2 μM of each probes (MWG-Biotech), 10× NH4 buffer (Bioline), 1.0 unit of BIOTAQ™ DNA polymerase (Bioline) and molecular grade water. 5 μl DNA was added and a total of 25 μl reaction was run in a Rotorgene 3000 thermocycler (Corbett, Australia). A cut off of 40 cycles was used to define positive samples. International Standard (INT) for P. falciparum DNA , which contains white blood cell DNA from the original donor specimen, was used as a positive control while parasite negative blood and water were used as negative controls. All reactions were performed in triplicate. The amplification and florescence detection was performed with the following two-step conditions: 95°C for 6 min; 40 cycles of 95°C for 15 s and 68°C for 1 min.
Patient 1 (OX007) presented with symptomatic malaria at HTD after returning from attending school in Lagos, Nigeria, and was found to harbour a mixed infection of P. falciparum (at 3.2% parasitaemia) and both gametocytes and trophozoites of P. ovale curtisi (at less than 0.03% parasitaemia), confirmed by PCR. The patient was treated with intravenous (i.v.) quinine (QN) for 1 day, oral QN for 5 days and given a full dose of sulphadoxine-pyrimethamine on discharge.
Patient 2 (OX005) developed symptomatic malaria after visiting Ghana, presenting to a local clinic in outer London one week after symptoms began. Unfortunately, treatment was further delayed due to a misdiagnosis of pandemic influenza (H1N1). The patient was seen again after four days and was immediately hospitalized, with the following signs of severe falciparum malaria: renal function impairment, hepatic function impairment, acidosis, antigen test for P. falciparum positive. The patient was started on i.v. QN, and blood films were taken and sent by courier to HTD. A diagnosis of falciparum malaria was confirmed by the on-call parasitology biomedical scientist, with a parasitaemia of 3%. The patient's condition had further deteriorated and transfer to the Intensive Care Unit (ICU) at HTD was arranged. Therapy with i.v. QN was continued for five days, and metabolic derangements gradually improved under ICU management. By the fourth day of treatment, only gametocytes of P. falciparum were observed on blood films. Five sequential daily blood samples were available.
Patient 3 (OX006) spent two weeks in Kenya based in Nairobi, and visiting Lake Naivasha, an area considered low risk for malaria. On advice, no prophylaxis was taken. The patient presented at night to a hospital outside London, and was antigen test positive for P. falciparum. Blood-films were taken and immediately couriered to the on-call parasitology biomedical scientist who identified extreme hyperparasitaemia of P. falciparum (approximately 2 × 106 parasites μl-1; 39% parasitaemia). To ensure early commencement of treatment, prior to transfer of the patient, the same motorcycle courier was dispatched in haste back to the referring hospital with i.v. artesunate (AS), which was administered just prior to departure for HTD. The patient was transferred to ICU, and treatment with i.v. AS continued. The patient exhibited the following signs of severe malaria: severe acidosis, renal impairment, deranged hepatic function, respiratory distress and haematuria. Remarkably, the patient remained conscious and lucid throughout his period of care in ICU. Haemofiltration was required for renal failure, but the initially high parasitaemia fell rapidly under treatment. AS was given up to day 5, followed by a full oral course of doxycycline. The patient was transferred from ICU to a standard medical ward on day 5, and discharged on day 7. Six sequential samples were available for analysis.
Patient 4, a British oil industry worker, presented with uncomplicated falciparum malaria in Port Harcourt, Nigeria, while working as a contractor in Nigeria and Angola. Quinine therapy administered in Port Harcourt failed to resolve a reported high parasitaemia, and the patient was transferred to HTD, London, by his employer, as a precaution. The patient reported previous bouts of malaria, but dates were not specified. P. falciparum infection was confirmed, with an initial parasite density of 0.5% on arrival in London. Oral QN therapy was continued at HTD for 48 hours (five doses in total) but poor response was again observed, necessitating a switch to oral atovaquone-proguanil (AP) (four tablets daily, for three days). Five sequential daily blood samples were available for analysis. The patient was clinically well by day 4 and was discharged before complete parasite clearance, with the remaining AP dose to be taken upon arrival in Nigeria.
All patients responded to treatment and were discharged after clinical improvement and reduction of parasite density to at or near zero.
Microscopic estimates of parasite density
Parasite density by microscopy of daily blood samples from each patient
Estimating PCT by qPCR
where n is time after treatment and time 0 is defined as the time at which the pre-treatment diagnostic sample was taken, or if that was unavailable, as for patients 3 and 4, the time at which the first sample in our series was taken.
The parasite density in each sample relative to the calibrator sample (day 0) was then estimated using the formula 2-ΔΔCT. The relative parasite abundance was plotted over time (Figure 1b). The qPCR plots are similar to those based on microscopy data (Figure 1a) but provide additional quantitative data in the later samples (days 3, 4 post-treatment) due to the higher sensitivity of qPCR.
Comparison of measured and estimated parameters of parasite clearance
A. PCT (hr):
B. PRR 48 h (ratio, no units):
qPCR * observed
C. PCT 95 (hr):
qPCR * observed
The current study presents the design and validation of a rapid, highly sensitive and potentially high-throughput assay to detect and quantify malaria parasites. Compared to expert double-read microscopy, the assay successfully detected and accurately measured malaria parasite densities in four individual patients treated with different anti-malarial regimens. The sensitivity, accuracy and reproducibility of the PCT assay have been further validated against the WHO International Standard for P. falciparum DNA, which is essentially lyophilized whole blood from a malaria patient, and thus provides both human and parasite targets for our duplex qPCR . The assay can detect Plasmodium densities as low as five parasites per μl of peripheral blood; parasite DNA was detected in two of the post-treatment samples that were reported negative by microscopy, and in a third sample reported to harbour scanty gametocytes of P. falciparum. This appears to be the first real-time PCR assay specifically designed to measure parasite clearance in follow-up of treated malaria patients, and has distinct advantages over previous methods used for this purpose, such as sybrG detection of 18 S ribosomal DNA sequences . The pgmet target has the advantage of being equally suitable for amplifying all Plasmodium species, and quantifies total parasite load in multi-species infections. It is a relatively simple matter to adapt the assay to any specific single parasite species of choice.
The assay was designed specifically for use with dried blood-spot samples collected simply in field situations, without measuring blood sample volumes, which would then be transferred to a laboratory with qPCR capacity. As a first step to validating this approach, the sensitivity, specificity and quantitative accuracy of the assay has been demonstrated here in a UK reference laboratory for parasitology, against high quality microscopic diagnosis. The efficacy of anti-malarial treatment is typically assessed by monitoring the clearance of parasites from the peripheral blood, and the speed at which symptoms and signs resolve. The utility of the assay to compare anti-malarial treatment response in different patients in this small pilot study has been demonstrated, and clear differences shown in PCT, PRR48 h and PCT95 among the four patients (Table 3, Figure 3). The ability to discern details of parasite clearance dynamics in individuals demonstrates the likely utility of this assay for detection of P. falciparum isolates with reduced sensitivity to ACT treatment. The qPCR method is now being successfully used to estimate PCT and PRR48 h using DNA extracted from dried filter paper blood spots. These have been collected during follow-up of a randomised efficacy study of two ACT regimens in 300 children with malaria in Mbita, western Kenya in 2009, similar to earlier studies in this area .
Dondorp and colleagues  used closely-spaced microscopic examination of blood films in small groups of malaria patients to demonstrate unusually slow parasite clearance time in some P. falciparum infected individuals from western Cambodia treated with AS. This parasite clearance phenotype could be explained by heritable changes in parasite growth rate, but may also indicate evolution of new mechanisms of resistance to artemisinin compounds [16, 17]. The duplex qPCR assay, capable of relatively high throughput, and thus incorporation in larger studies, should assist in performing association studies to identify candidate genetic markers in P. falciparum of the artemisinin slow clearing parasite phenotype. Such markers, once validated, would be of enormous benefit in monitoring and containing resistance to ACT therapy.
Study protocols are now being developed that incorporate the qPCR PCT assay as an endpoint in large scale clinical trials of ACT efficacy in Africa. Initial testing of the assay on filter-paper DNA has shown promising results, and thus there is optimism that the method can be adopted to detect the phenotypic characteristics of parasites in clinical trials during the three days immediately after treatment. As the proportion of positive microscopy readings for ACT-treated patients on day 3 is generally less than 3% , a particular advantage of the qPCR PCT assay is its potential to detect and quantify parasites in many more patients on that day (Figure 3).
It is perhaps not surprising that the performance of the qPCR assay was markedly different in patient 3, who presented with extreme hyperparasitaemia, but recovered under AS monotherapy. There appeared to be a persistence of parasitized erythrocytes as detected by PCR for longer, and at greater estimated densities, than was expected from the microscopy results. It is unlikely the presence of peripheral gametocytes were responsible for this qPCR result, as the microscopists did not indicate the presence of sexual stages of P. falciparum, an observation routinely recorded. The higher parasite density estimated by qPCR could be due to the well-described circulation of more mature parasite forms, with significantly higher DNA content per cell, in severe malaria. Dead parasites may also contain detectable parasite DNA, but are not expected to circulate in a patient with an intact spleen, although pigment-positive phagocytes may also contain some detectable parasite DNA in parasites with very high parasitaemia. Dead parasites appear to be of no consequence in comparisons of qPCR-derived parasite densities with microscopy in uncomplicated P. falciparum infections in which parasite density is in the normal range . Thus, it is expected that the qPCR assay will be most useful in studies of treatment response in uncomplicated malaria patients in the absence of hyperparasitaemia.
This study provides proof of principle that a duplex fluorescent qPCR assay can be deployed to robustly estimate the rate of parasite clearance in treated malaria patients. The particular assay used here has the advantage of deploying a parasite target sequence that is identical in all Plasmodium species investigated and thus provides an estimate of total parasite burden. The method could be easily adapted to use specific amplification targets for each Plasmodium species.
List of abbreviations
artemisinin combination therapy
Hospital for Tropical Diseases, London
human b-tubulin gene sequence
intensive care unit
parasite clearance time
Plasmodium genus methionine tRNA gene sequence
parasite reduction ratio
quantitative (real-time) PCR
This work was supported by the EDCTP project IP.2007.31060.002, "An integrated approach to clinical trials, capacity building and networking in West Africa", co-ordinated by Abdoulaye Djimde, University of Bamako. CJS is supported by the UK Health Protection Agency. SDP is supported by the Foundation for Innovative New Diagnostics. PLC is supported by the UCL Hospitals Comprehensive Biomedical Research Centre Infection Theme. These funding sources played no role in the collection, analysis or interpretation of the data; nor in the writing of the manuscript, or in the decision to submit the manuscript for publication. We thank Maggie Armstrong for assistance in accessing patients' notes.
- WHO: Guidelines for the treatment of malaria. 2010, World Health Organization, Geneva, ISBN: 9789241547925, SecondGoogle Scholar
- White NJ: Assessment of the pharmacodynamic properties of antimalarial drugs in vivo. Antimicrob Agents Chemother. 1997, 41: 1413-22.PubMed CentralPubMedGoogle Scholar
- White NJ, Nosten F, Looareesuwan S, Watkins WM, Marsh K, Snow RW, Kokwaro G, Ouma J, Hien TT, Molyneux ME, Taylor TE, Newbold CI, Ruebush TK, Danis M, Greenwood BM, Anderson RM, Olliaro P: Averting a malaria disaster. Lancet. 1999, 353: 1965-1967. 10.1016/S0140-6736(98)07367-X.View ArticlePubMedGoogle Scholar
- Alker AP, Lim P, Sem R, Shah NK, Yi P, Bouth DM, Tsuyuoka R, Maguire JD, Fandeur T, Ariey F, Wongsrichanalai C, Meshnick SR: Pfmdr1 and in vivo resistance to artesunate-mefloquine in falciparum malaria on the Cambodian-Thai border. Am J Trop Med Hyg. 2007, 76: 641-647.PubMedGoogle Scholar
- Noedl H, Se Y, Schaecher K, Smith BL, Socheat D, Fukuda MM, Artemisinin Resistance in Cambodia 1 (ARC1) Study Consortium: Evidence of artemisinin-resistant malaria in western Cambodia. N Engl J Med. 2008, 359: 2619-2620. 10.1056/NEJMc0805011.View ArticlePubMedGoogle Scholar
- Dondorp AM, Nosten F, Yi P, Das D, Phyo AP, Tarning J, Lwin KM, Ariey F, Hanpithakpong W, Lee SJ, Ringwald P, Silamut K, Imwong M, Chotivanich K, Lim P, Herdman T, An SS, Yeung S, Singhasivanon P, Day NP, Lindegardh N, Socheat D, White NJ: Artemisinin resistance in Plasmodium falciparum malaria. N Engl J Med. 2009, 361: 455-467. 10.1056/NEJMoa0808859.PubMed CentralView ArticlePubMedGoogle Scholar
- Stepniewska K, Ashley E, Lee SJ, Anstey N, Barnes KI, Binh TQ, D'Alessandro U, Day NP, de Vries PJ, Dorsey G, Guthmann JP, Mayxay M, Newton PN, Olliaro P, Osorio L, Price RN, Rowland M, Smithuis F, Taylor WR, Nosten F, White NJ: In vivo parasitological measures of artemisinin susceptibility. J Infect Dis. 2010, 201: 570-579. 10.1086/650301.PubMed CentralView ArticlePubMedGoogle Scholar
- Hallett RL, Dunyo S, Ord R, Jawara M, Pinder M, Randall A, Alloueche A, Walraven G, Targett GAT, Alexander N, Sutherland CJ: Treatment of malaria in Gambian children with chloroquine plus sulphadoxine-pyrimethamine favours survival and transmission to mosquitoes of multi-drug-resistant Plasmodium falciparum. PLoS Clin Trials. 2006, 1: e15-10.1371/journal.pctr.0010015.PubMed CentralView ArticlePubMedGoogle Scholar
- Gadalla NB, Elzaki SE, Mukhtar E, Warhurst DC, El-Sayed B, Sutherland CJ: Dynamics of pfcrt alleles CVMNK and CVIET in chloroquine-treated Sudanese patients infected with Plasmodium falciparum. Malaria J. 2010, 9: 74-10.1186/1475-2875-9-74.View ArticleGoogle Scholar
- Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2-ΔΔCT Method. Method. 2001, 24: 402-408. 10.1006/meth.2001.1262.View ArticleGoogle Scholar
- Sutherland CJ, Fifer H, Pearce RJ, bin Reza F, Nicholas M, Haustein T, Njimgye-Tekumafor NE, Doherty JF, Gothard P, Polley SD, Chiodini PL: Novel pfdhps haplotypes among imported cases of Plasmodium falciparum malaria in the UK. Antimicr Agents Chemother. 2009, 53: 3405-3410. 10.1128/AAC.00024-09.View ArticleGoogle Scholar
- MapSeq/pf. (last accessed 03/11/2010), [http://www.sanger.ac.uk/MapSeq/]
- Sutherland CJ, Tanomsing N, Nolder D, Oguike M, Jennison C, Pukrittayakamee S, Dolecek C, Hien TT, do Rosário VE, Arez AP, Pinto J, Michon P, Escalante AA, Nosten F, Burke M, Lee R, Blaze M, Otto TD, Barnwell JW, Pain A, Williams J, White NJ, Day NP, Snounou G, Lockhart PJ, Chiodini PL, Imwong M, Polley SD: Two nonrecombining sympatric forms of the human malaria parasite Plasmodium ovale occur globally. J Infect Dis. 2010, 201: 1544-1550. 10.1086/652240.View ArticlePubMedGoogle Scholar
- Padley DJ, Heath AB, Sutherland C, Chiodini PL, Baylis SA, the Collaborative Study Group: Establishment of the 1st World Health Organization International Standard for Plasmodium falciparum DNA for nucleic acid amplification technique (NAT)-based assays. Malar J. 2008, 7: 139-10.1186/1475-2875-7-139.PubMed CentralView ArticlePubMedGoogle Scholar
- Bousema JT, Schneider P, Gouagna LC, Drakeley CJ, Tostmann A, Houben R, Githure JI, Ord R, Sutherland CJ, Omar SA, Sauerwein RW: Moderate effect of artemisinin-based combination therapy on transmission of Plasmodium falciparum. J Infect Dis. 2006, 193: 1151-1159. 10.1086/503051.View ArticlePubMedGoogle Scholar
- White NJ: Qinghaosu (artemisinin): the price of success. Science. 2008, 320: 330-334. 10.1126/science.1155165.View ArticlePubMedGoogle Scholar
- Anderson TJ, Nair S, Nkhoma S, Williams JT, Imwong M, Yi P, Socheat D, Das D, Chotivanich K, Day NP, White NJ, Dondorp AM: High heritability of malaria parasite clearance rate indicates a genetic basis for artemisinin resistance in western Cambodia. J Infect Dis. 2010, 201: 1326-1330. 10.1086/651562.PubMed CentralView ArticlePubMedGoogle Scholar
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.