- Open Access
Plasmodium falciparum resistance to anti-malarial drugs in Papua New Guinea: evaluation of a community-based approach for the molecular monitoring of resistance
Malaria Journal volume 9, Article number: 8 (2010)
Molecular monitoring of parasite resistance has become an important complementary tool in establishing rational anti-malarial drug policies. Community surveys provide a representative sample of the parasite population and can be carried out more rapidly than accrual of samples from clinical cases, but it is not known whether the frequencies of genetic resistance markers in clinical cases differ from those in the overall population, or whether such community surveys can provide good predictions of treatment failure rates.
Between 2003 and 2005, in vivo drug efficacy of amodiaquine or chloroquine plus sulphadoxine-pyrimethamine was determined at three sites in Papua New Guinea. The genetic drug resistance profile (i.e., 33 single nucleotide polymorphisms in Plasmodium falciparum crt, mdr1, dhfr, dhps, and ATPase6) was concurrently assessed in 639 community samples collected in the catchment areas of the respective health facilities by using a DNA microarray-based method. Mutant allele and haplotype frequencies were determined and their relationship with treatment failure rates at each site in each year was investigated.
PCR-corrected in vivo treatment failure rates were between 12% and 28% and varied by site and year with variable longitudinal trends. In the community samples, the frequencies of mutations in pfcrt and pfmdr1 were high and did not show significant changes over time. Mutant allele frequencies in pfdhfr were moderate and those in pfdhps were low. No mutations were detected in pfATPase6. There was much more variation between sites than temporal, within-site, variation in allele and haplotype frequencies. This variation did not correlate well with treatment failure rates. Allele and haplotype frequencies were very similar in clinical and community samples from the same site.
The relationship between parasite genetics and in vivo treatment failure rate is not straightforward. The frequencies of genetic anti-malarial resistance markers appear to be very similar in community and clinical samples, but cannot be used to make precise predictions of clinical outcome. Thus, indicators based on molecular data have to be considered with caution and interpreted in the local context, especially with regard to prior drug usage and level of pre-existing immunity. Testing community samples for molecular drug resistance markers is a complementary tool that should help decision-making for the best treatment options and appropriate potential alternatives.
The development and spread of Plasmodium falciparum resistance to the most commonly used anti-malarial drugs is a major challenge in the control of malaria . Quantifying the level of drug resistance through regular monitoring provides essential information to the health authorities responsible for ensuring ready access to effective drugs. Methods for assessing drug resistance in malaria include in vivo drug efficacy studies (the gold standard), in vitro sensitivity testing of patient isolates, and surveys designed to detect molecular markers of drug resistance. Each approach has advantages and disadvantages (discussed in [2–4]) and herein field data evaluating the usefulness of the molecular marker approach are presented.
Molecular screening to detect drug resistance was made possible by studies identifying the genetic basis of drug resistance in P. falciparum and this approach has subsequently become an integral part for the evaluation of resistance to treatment [5, 6]. Chloroquine resistance (CQR) is attributable to single nucleotide polymorphisms (SNPs) in pfcrt and pfmdr1 (reviewed in [7–9]). Resistance to sulphadoxine-pyrimethamine (SP) is associated with a stepwise accumulation of mutations in pfdhfr and pfdhps (reviewed in ). Artemisinin derivates have been shown to inhibit pfATPase6 and SNPs in pfATPase6 may become associated with resistance in vitro[12, 13]. It has also become apparent that these mutations do not act in isolation, but often act synergistically to encode or enhance resistance. Resistance may arise through sequential accumulation of mutations in a single gene, such as anti-folate resistance in the pfdhfr gene [4, 14]. Mutations in different genes may also act synergistically. Linkage disequilibrium is often observed between pfcrt K76T and pfmdr1 N86Y in vivo[15–17] and the commonly accepted hypothesis is that CQR is mediated by multigenic processes , primary resistance being encoded by pfcrt K76T with SNPs in pfmdr1 playing an important role in modulating levels of CQR. The requirement to consider intra- and inter-genic genotypes impose limitations on the use of molecular methods such as PCR-RFLP which can identify the presence of mutations at individual markers, but not individual malaria genotypes defined at multiple loci in large sample sets.
The link between molecular markers and resistance can be clearly demonstrated in the laboratory, but the usefulness of these markers in the field has been controversial, primarily because the strength of their associations with in vivo treatment outcome is not consistent in different epidemiological settings. There are two plausible reasons for this . Firstly, drug failure may occur for reasons other than parasite genetics. Patients may take incomplete drug courses, may have poor drug absorption or metabolism and so on. When marker frequencies are low, these human sources of drug failure may obscure the relationship between marker and resistance. Secondly, most studies have primarily focused on single genes and markers , rather than multiple markers which may act synergistically. A further limitation is that most analyses of molecular markers have been conducted as part of clinical studies and have, therefore, been limited to sentinel sites with good access to health care facilities and mostly examined clinical cases within a restricted age group . Moreover, most such studies used parasite genetic measures based on the proportions of patients carrying a given genetic marker, or composite measures formed from ad hoc functions of such proportions . This ignores the fact that patients often carry multiple clones of parasites and so may carry a mutation even if most of their parasites are wild-type. In principle, the evolution of drug resistance should be assessed by genetic profiling of representative samples of the circulating parasite population, which includes the reservoir in asymptomatic carriers. This should be summarized using population-genetic measures of allele frequencies (for single alleles) and haplotype frequencies (for multiple linked alleles). However, genetic data from community samples have rarely been compared with in vivo treatment responses [21, 22]. In addition, it is not clear whether the molecular profile of parasites in symptomatic patients seen at health facilities matches that of the large circulating parasite reservoir.
The current study has investigated the role and applicability of the molecular drug resistance profiles in community samples for the monitoring of drug resistant malaria in Papua New Guinea (PNG). For this purpose, in vivo drug efficacy studies with the first-line regimen of amodiaquine or chloroquine plus sulphadoxine-pyrimethamine (AQ+SP or CQ+SP) were conducted at three different sites. A DNA microarray-based technology was used to compare the molecular drug resistance profile in clinical samples and blood samples collected in the communities from the catchment areas of the corresponding health facilities. Finally, the relationships of the allele and haplotype frequencies of the molecular markers with the rates and time-trends in clinical failure at the different sites were investigated.
Study areas and design
The three study sites in Papua New Guinea (PNG) were 1) the Sigimaru health centre (HC) in the Karimui area (Simbu Province), 2) the Kunjingini HC in the South Wosera area (East Sepik Province), and 3) the Mugil HC in the North Coast area of Madang (Madang Province). In Karimui, a rural region in the highland fringe area of PNG, the studies were run between October and April in three consecutive years (2003, 2004, and 2005). In the Wosera, located in the floodplain of the Sepik river in the north-eastern part of the country bordering Indonesian Papua, the study period was between December and June in 2003 and 2004. The study in the rainforest area at the North Coast of Madang was conducted between April 2004 and February 2005.
Malaria transmission is perennial with limited variations between wet (October to April) and dry (May to September) season at all sites. Transmission intensity in PNG decreases significantly with increasing altitude  and is higher in the lowland regions of the Wosera and the North Coast than in the Karimui area, an elevated plateau situated at an altitude of 700 to 1,200 m. There is little socioeconomic stratification between and within sites, although there are differences in health care provision and drug use patterns [24–28]. Baseline characteristics of the study sites are summarized in Table 1.
Assessment of in vivo drug efficacy
Drug efficacy studies were conducted according to the standardized WHO protocol for low to moderate transmission areas  and are described in detail elsewhere . Patients were classified according to their clinical and parasitological responses into early treatment failure (ETF), late clinical failure (LCF), late parasitological failure (LPF), or adequate clinical and parasitological response (ACPR) and crude failure rates were PCR-corrected by comparing the msp2 (merozoite surface protein 2) genotyping patterns of paired Day 0 and treatment failure samples . Annual treatment failure rates (TFR = ETF + LCF + LPF) of AQ+SP and CQ+SP were pooled for the analysis in the current study.
Community-based cross-sectional surveys
Cross-sectional surveys were conducted in the catchment areas of each of the health centres using a randomized household approach. To obtain a representative random sample of the parasite population circulating in the corresponding communities (i.e., approximately 100 PCR-positive P. falciparum samples from each location), between 300 and 350 blood specimens were collected per community . Overall, a total of 2,013 individuals from randomly selected villages and households were recruited for the community surveys between 2003 and 2005.
Apart from demographic characteristics, collected information included history of sickness (onset, type and duration of symptoms), health facility attendance, purchase or consumption of drugs outside health facilities, and anti-malarial treatment courses received in the preceding year (extracted from health books when available). Axillary temperature was measured with an electronic thermometer and spleen size was assessed in the recumbent position using Hackett's grading system . Blood samples for parasitological examination by microscopy, haemoglobin (Hb) level determination (HemoCue®, Ångelholm, Sweden), and molecular identification of parasite genotypes were collected by venepuncture using 2 ml EDTA-Vacutainer™ tubes (BD Biosciences, Allschwil, Switzerland). Plasma was separated by centrifugation and red blood cell pellets were stored frozen until further processing.
DNA was extracted using QIAamp® DNA Blood Kit (Qiagen, Hombrechtikon, Switzerland) according to the manufacturer's instructions. Assessment of SNPs for drug resistant malaria was done for pfmdr1 mutations N86Y, Y184F, S1034C, N1042D and D1246Y, pfcrt mutations K76T, H97Q, T152A, S163R, A220S, Q271E, N326D/S, I356L/T and R371I, pfdhfr mutations A16V, N51I, C59R, S108N/T and I164L, pfdhps mutations S436A, A437G, K540E, A581G and A613T/S, and pfATPase6 mutations S538R, Q574P, A623E, N683K and S769N. The method is based on parallel PCR amplification of the target sequences followed by primer extension mediated mini-sequencing using fluorochrome-labelled ddNTPs. Subsequent base calling occurs on a microarray upon sequence specific hybridization [33, 34]. Multiplicity of infection (MOI) for each sample was assessed by determining the number of msp2 genotypes .
Data were double entered in EpiData software (version 3.02, Odense, Denmark) and analysis was performed using STATA software (version 10.1; Stata Corp., College Station, Texas). Allele frequencies were calculated separately for each site for the community samples collected in each year, and separately for the health centre samples from both, Karimui and South Wosera collected in 2003. The estimates were very similar to those obtained using a corresponding Bayesian algorithm . Confidence intervals were calculated using bootstrap sampling, in each case based on 10,000 bootstrap samples of the data.
As a measure of genetic differentiation between populations and sample sets , the estimate of the GST statistic , obtained using the maximum likelihood estimates of the allele frequencies, was used. Overall tests of the null hypothesis that allele frequencies at all four loci result from a sample from a single population were carried out using randomization tests of , in each case using 10,000 random re-assortments of the data to sub-populations.
The frequencies of multi-SNP haplotypes were calculated using the MalHaploFreq programme as described in detail elsewhere . Haplotype frequencies were estimated for each of the main combinations of SNPs possible within each locus. In addition, two-locus haplotype frequencies were estimated for SNPs at the loci involved in SP resistance (i.e., pfdhfr and pfdhps). Locus-specific comparisons of haplotype frequencies in the different sample sets were carried out using likelihood ratio tests. The likelihood ratio statistic was computed as twice the difference between the log-likelihood of the full model (estimating distinct haplotype frequencies for each sample set) and that of a model estimating common haplotype frequencies for all sample sets included in the comparisons.
Scientific approval and ethical clearance for the study was obtained from the Ethikkommission beider Basel (EKBB) and the Medical Research and Advisory Committee (MRAC) of the Ministry of Health in PNG. Informed consent was first requested from all the communities involved and prior to recruitment, individual consent was obtained from each study participant and parents or legal guardians.
Characteristics of community survey population
Key characteristics of the sample sets collected from each site between 2003 and 2005 are given in Table 1. Age distribution, sex ratio, and mean axillary temperature were similar between sites, and mean haemoglobin levels were lower in the high transmission areas of the Wosera and the North Coast. Microscopic P. falciparum prevalence rates ranged between 13% and 27%, depending on the site and the respective year, and from a total of 2013 community blood samples, 695 (34.5%) were found to be positive for P. falciparum msp 2 by PCR and prevalence rates determined by a positive msp2 PCR ranged from 20% to 41%. MOI determined by msp2 PCR was not significantly different between and within sites over time (range: 1.5 to 1.8), but the proportion of multiple versus single infections varied between high and moderate transmission areas. Multiple infections in either msp 2 or in one or more resistance markers were detected in 436 out of 643 community samples tested, implying the presence of at least 1,247 parasite clones overall.
In vivo drug efficacy
Clinical and parasitological monitoring up to Day 28 between 2003 and 2005 at all three sites was accomplished for a study group of 649 children (Table 1, ). In Karimui, the pooled (AQ+SP and CQ+SP) PCR-corrected TFR tended to decrease over the three-year period from 28% to 18% and 16%, respectively, but this was not statistically significant ( = 4.8, p = 0.09). In the Wosera, the TFR tended to increase, changing from 16% in 2003 to 22% in 2004, but again this was not statistically significant ( = 1.2, p = 0.3). In 2004, the TFR at the North Coast area of Madang was 11.5%.
Allele and haplotype frequencies in community surveys
Analysis of drug resistance mutations was carried out for 639 (92%) of the community samples. Polymorphisms were found in pfmdr1 SNPs N86Y, Y184F and N1042D, pfcrt SNPs K76T, A220S, N326D, I356L and S163R, pfdhfr SNPs C59R and S108N, and pfdhps SNPs A437G and K540E. None of the other SNPs was detected as mutated allele (Table S1, Additional file 1). The proportions of samples found to contain mutant alleles depends on both, the proportion of parasite clones with the mutation (the mutant allele frequency) and the MOIs of the samples [35, 37]. The estimated mutant allele frequencies were calculated using the typing results (Table S1, Additional file 1) and the MOIs determined by msp 2 typing and are given in table S2 (Additional file 2).
High mutant allele frequencies for the CQ relevant markers in pfcrt were observed for K76T, N326D and I356L at all three sites (Table S2, Additional file 2). Frequencies of A220S were highest at the North Coast and lower in Karimui and the Wosera. Mutated pfcrt S163R alleles were rarely detected and only in Karimui and the Wosera. Frequencies for pfmdr1 N86Y reached quasi fixed levels in Karimui and at the North Coast, but were substantially lower in the Wosera. Pfmdr1 polymorphisms Y184F and N1042D were never detected in samples from Karimui and frequencies were low at the other two sites. Mutant allele frequencies for the SP relevant markers pfdhfr S108N and C59R were high at all three sites, but were moderate for pfdhps A437G and very low for pfdhps K540E.
Haplotype frequencies were estimated for the three SNP combinations pfmdr1 N86Y+Y184F+N1042D and pfcrt K76T+A220S+I356L, and for each of the possible haplotypes of pfdhfr plus pfdhps arising from the polymorphic SNPs (Table S3, Additional file 3).
Comparison between sites
The analysis of between-site genetic variation, based on comparing the allele frequencies in all community samples available from each of the three sites, gave a value of 0.082, which was greater than the value obtained in any of the 10,000 randomizations of the data. This thus corresponds to a P-value < 0.01 suggesting rejection of the null hypothesis that the community samples come from a single parasite population. This value was considerably higher than the values obtained by comparing between years within sites (Table 2).
At the pfcrt locus, haplotype frequencies of the wild-type, single K76T or double K76T+I356L mutants were observed with frequencies between 0.02 to 0.36 in Karimui and the Wosera, but were significantly lower in samples from the North Coast, where the haplotype with the triple mutation K76T+A220S+I356L was close to fixation (Table S3, Additional file 3 and Figure 1). At pfmdr1, the haplotype comprising the single mutation N86Y predominated at all three sites and was close to fixation in Karimui and the North Coast, but less frequent (0.68) in the Wosera. Fully wild-type pfdhfr/pfdhps haplotypes were infrequent at all three sites (0.01 to 0.15) and the predominant haplotype at pfdhfr was the double mutant S108N+C59R. The pfdhps A437G mutation (which always occurred alongside the pfdhfr double mutation S108N+C59R) was found in a substantial proportion of samples from the Karimui area, but was much less frequent at the other two sites.
Comparison of clinical and community samples
Comparison of the drug resistance profile in community and clinical samples was only possible in the Karimui and South Wosera areas in 2003. In Karimui, the allele frequencies in health centre samples were close to those from community samples for the same year (Table S4, Additional file 4 and Figure 2). In the Wosera, there were small effects of both, year and source of the samples on the allele and haplotype frequencies (both values 0.011). Despite the rather similar allele frequencies in community and health centre samples, the differences in haplotype frequencies between health centre and community samples were statistically significant for each locus except pfmdr1 (Table 2).
Comparison between years
In Karimui, there was a significant difference between years in the vector of allele frequencies in the community samples, and the value of 0.035 for measuring inter-annual variation was relatively high for each locus except pfmdr1, where a single haplotype appeared to be at fixation. The most obvious changes were related to the allele frequencies of pfcrt N326D, I356L and A220S, which all decreased over time, as did pfdhps A437G, while the frequencies of the pfdhfr mutants increased (Table S2, Additional file 2). In terms of haplotypes, pfdhfr and pfdhps appeared to be moving towards fixation of the same haplotypes as were already fixed at the North Coast, while there was no indication that any pfcrt haplotype was moving towards fixation (Table S3, Additional file 3 and Figure 1).
In the Wosera, the value of measuring inter-annual variation was lower (0.011), but the variation in haplotype frequencies was statistically significant for all polymorphic loci. While mutations in pfmdr1 appeared to be relatively stable, those in the other loci appeared to be moving towards fixation of the same haplotypes as at the North Coast. Triple mutations in pfcrt increased in frequency between 2003 and 2004, as did the frequencies of the pfdhfr S108N+C59R double and pfdhps A437G single mutants (Table S3, Additional file 3 and Figure 1).
Relationship of allele and haplotype frequencies with treatment failure rates
These differences in allele and single-locus haplotype frequencies between sites were not reflected in treatment failure rates (TFR) (Figure 3). In particular, the North Coast site had the lowest TFR, but the highest frequencies of triple mutants in pfcrt, of double mutants in pfdhfr and of the main polymorphic mutation in pfmdr1. Despite the substantial differences in allele and haplotype frequencies between the Wosera and Karimui, the average TFRs at these two sites were similar.
The allelic and single-locus haplotype frequencies in the Wosera seemed to be converging over time with those of the North Coast parasite population (Figure 1), while the TFR in the Wosera was increasing (i.e., diverging from that in the North Coast). At the same time, the general trend in single-locus pfdhfr and pfdhps mutant haplotype frequencies in Karimui towards the North Coast values accompanied a decrease in TFRs. Thus, the allele frequencies in the Wosera and Karimui could not account for the patterns in TFR. However, a strikingly different pattern was seen with the two-locus haplotypes. While the clinical TFR correlated negatively with the pfdhfr double mutant haplotype when this was not linked to the pfdhps A437G mutation, the frequency of the pfdhfr C59R+S108N/pfdhps A437G mutant haplotype increased with the TFR (Figure 3).
The present study was specifically designed to analyse the relationship between parasite genetics in the community and in vivo failure rates, but the complexity of this relationship was not initially appreciated. There is substantial evidence linking specific genetic markers to drug resistance in vitro and the presence of genetic resistance mutations often predicts treatment failure at the level of the individual patient, but it was observed that these individual level relationships are not reflected in straightforward relationships at the population level between molecular typing results and in vivo tests of resistance. A number of possible reasons for this discrepancy were investigated.
One limitation of many studies is that they have included only small numbers of SNPs. The chosen SNPs might prove to be the wrong mutations in the given population or secondary mutations with indirect relationships to resistance might be important. This was addressed by using a DNA microarray-based technology that allowed parallel assessment of many SNPs in several genes . The analysis of a total of 12 polymorphic SNPs, as well as a further 21 SNPs that proved to be uniformly wild-type, enabled to consider all candidate SNPs so far identified. The study populations had highly mutated genes involved in CQR resistance (pfcrt and mdr1) combined with pfdhfr mutations consistent with a moderately pyrimethamine resistant phenotype and the emergence of key mutations in pfdhps. The key CQR markers pfcrt K76T, N326D and N326D were almost fixed at all three sites, in keeping with other recent findings [38–40]. Similarly, the high allele frequencies for pfdhfr S108N and C59R and the resulting high frequency of the double-mutant haplotype are consistent with a recent report from East Sepik Province . Despite this substantial amount of information about each sample, the allele and single-locus haplotype frequencies could still not easily be related to in vivo failure rates.
The choice of strategy for sampling the parasite population does not seem to be an important determinant of whether a relationship is seen with treatment failure. The differences between corresponding genetic profiles in health facility and community samples were small in comparison with the differences between sites. In Karimui, the profiles were near identical, and the small differences between the two profiles in the Wosera, which may reflect heterogeneities in access and hence drug pressure for different parts of the catchment area, are too small to influence the choice of sampling strategy. The decision as to whether to sample the community or health centres should therefore be made on logistic grounds. Community sampling can be carried out during a single short visit to a site, while health facility sampling is likely to take longer, but may be more convenient especially if blood sampling for drug resistance can be combined with other procedure(s).
A more important factor in being able to predict treatment failure from parasite genetic data may be the appropriate analysis of the data. Only by analysing two-locus pfhfr/pfdhps haplotypes coherent relationships between the genetic data and TFRs could be seen, corroborating the previous finding that the pfdhfr C59R+S108N/pfdhps A437G triple mutant predicts treatment failure best on individual level in PNG . While this haplotype was absent in the North Coast and so cannot account for the treatment failures observed there, increases in haplotype frequency in the other sample sets matched remarkably well the increased TFRs. This is in agreement with studies from Uganda  and Laos , that also report pfdhps mutations to be important markers for unsuccessful treatment response to combination therapy with CQ plus SP, although the pfdhps genotype was not indicative for treatment failure with SP monotherapy in Tanzania [44, 45] or Ghana .
It has been proposed [47–49] that genetic resistance can be measured by the ratio between the prevalence of mutated genotypes to the TFR (the genetic failure index, GFI). A low, reasonably stable GFI would indicate a good marker of resistance. The theoretical properties of GFI are discussed in detail elsewhere  where it was pointed out that the GFI for any given SNP could vary widely depending on genotypes at other SNPs modulating the parasites ability to survive therapy. This study quantified this effect, showing that the relationship between SNPs and clinical outcome needs to be based on haplotype frequencies rather than crude prevalence of a single SNP. The estimation of further multi-locus haplotype frequencies may be needed before the full potential of such analysis becomes clear.
Currently, there are technical limitations in the ability to estimate multi-locus haplotype frequencies. Where there are no polyclonal infections, haplotype frequencies can be estimated by simple gene counting. However, 68% of the PCR-positive community samples had multiple msp2 types, implying that only 17% of the parasite clones could validly be used to estimate allele frequencies by simple gene counting, making this a very inefficient and imprecise procedure. The estimation of haplotype frequencies from multi-clone samples using the MalHaploFreq programme makes use of substantially more information, but with the current implementation it is not feasible to identify the 3- and 4-locus haplotype frequencies.
Fortunately, there are two ways of avoiding this limitation. Firstly, to date, most resistance appears to have arisen through a few single key SNPs such as pfcrt K76T, pfcrt A220S and pfmdr1 N86Y, or by a series of mutations such as the pfdhfr 108-51-59-164 sequence. Thus, most key multi-locus haplotypes can be characterized at a few numbers of SNPs (e.g. pfcrt K76T+pfmdr1 N86Y, pfcrt K76T+dhfr I164L, and so on). Secondly, in areas of moderate to high transmission there is likely to be large numbers of mixed infections and hence large amounts of sexual recombination which tends to bring separate genes into linkage disequilibrium (LD). If the genes are in LD then the frequency of multi-gene haplotypes can be obtained simply by multiplying their individual frequencies. Alternately, if mixed infections are rare then haplotypes can be directly observed in single-clone infections and MalHaploFreq is not required.
The inter-site differences in TFR need to be understood in the context of the long-term histories of malaria transmission and drug use at each site, possibly including locality-specific factors, such as clinical practice, genetic background of human population and degree of isolation. After a long history of 4-aminoquinoline use, treatment policy against uncomplicated malaria in PNG changed to combination therapy with AQ or CQ plus SP only in 2000. Both, the TFR and allele/haplotype frequencies, including high frequencies of important markers in pfdhfr and pfdhps, are therefore the outcomes after only two years of implementation of combination therapy and are probably highly dependent on the baseline frequency of resistant mutations at the time of the policy change.
The sites differed in the frequency of pfcrt A220S corresponding to the highly CQ resistant pfcrt quadruple mutant haplotype (K76T+N326D+N326D+A220S), which was close to fixation only at the North Coast. Similarly, the pfmdr1 N86Y mutant was close to fixation in Karimui and at the North Coast (>92%), but was less frequent in the Wosera (71%), and the pfmdr1 polymorphisms Y184F and N1042D were absent in Karimui. These genetic differences might simply reflect stochastic differences in source parasite populations, but they could also reflect differences in historical drug pressure. In particular, while pfcrt mutations presumably reflect a history of CQ use, SNPs in pfmdr1 conferring both, resistance to CQ and altered sensitivity to other drug classes, including quinine (which was used for treatment of severe malaria in PNG), arylaminoalcohols and artemisinin derivates [37, 50], are presumably selected by exposure to these. The fixation of a highly CQ resistant haplotype at the North Coast may reflect a history of very high CQ pressure due to good health care provision as well as easy (and unregulated) access to drugs in the nearby town of Madang. In contrast, former 4-aminoquinoline use is assumed to be lower in the remote regions in Karimui and the Wosera, so it is possible that the North Coast pattern is one towards which the others are progressing as a result of drug pressure, yet the North Coast had the lowest TFR. However, former drug pressure with anti-folates in the course of previous mass drug administration campaigns was higher in the two more remote regions [51, 52], so a reservoir of anti-folate resistance may have been present there for a long time. Parasite CQR at all three sites may have reached a plateau, while resistance to SP evolved at a different speed depending on both, historical and current anti-folate drug pressure. The frequency of the pfcrt A220S mutation increased in the Wosera between 2003 and 2004, most likely as a consequence of inappropriate drug use exerted on a large and heterogeneous parasite population. This arose because of 1) health systems management problems in the area, which lead to an undersupply with SP (local health extension officer at the Kunjingini health centre, personal communication) and 2) drug consumption outside health facilities which was relatively frequent by local standards (10% and 12% of the survey populations in 2003 and 2004, respectively) with AQ or CQ monotherapy distributed by inadequately trained community health workers. These factors may have further enhanced CQR leading to a reduced capacity of CQ to protect SP. Within the same time period, decreasing trends of TFRs were accompanied by decreasing frequencies of pfcrt A220S and pfdhps A437G in Karimui. Health care provision generally improved due to better medical supply to this remote highland area and the study increased awareness of malaria. The result appears to have been effective delivery of combination therapy to a higher proportion of malaria cases, acting on a relatively small parasite population, therefore preventing the further development and spread of parasite resistance in this moderate transmission area [53, 54].
The in vivo and molecular data suggest that in the current first-line regimen in PNG, SP, and more specifically sulphadoxine is the effective component. Hence, molecular monitoring of resistance to this component is important under constant treatment policy. However, significant levels of in vivo failure are a strong argument for the introduction of an artemisinin-based combination regimen. Once such a regimen is introduced, rapid, easy-to-use and affordable surveillance systems will be needed to monitor emergence and spread of resistance to the components of the new combination. At the same time, surveillance of resistance to withdrawn drugs should also be continued, with the prospect that decreases in resistance may make it possible to reuse one or more of these safe and cheap drugs as partner compounds in combination regimens [55, 56].
In summary, despite there being considerable interest in using molecular data to predict drug failure rates, and hence inform and guide policy choice [5, 57], the current analyses showed that the relationship was not as simple as might be hoped. Molecular drug resistance profiles were similar in parasites from community and clinical samples, but these did not closely reflect the longitudinal and geographical variation in in vivo drug efficacy observed at health centres. The data confirm that the genetic drug resistance background of the parasite is only one of many factors determining clinical outcome and one which has evolved differently according to epidemiological characteristics as well as history and patterns of drug use in a given area. In particular it was shown that (i) using genetic markers to predict failure rates is much better if based on haplotype frequencies rather than crude prevalence of individual SNPs, and (ii) SNP marker sets will have to be constantly updated and adapted to future monitoring purposes, such as the inclusion of markers for newly introduced or withdrawn drug classes. Finally, this community-based molecular monitoring approach will have to be further evaluated in geographical areas at both extremes of transmission intensity and different drug use patterns (e.g. rural versus urban areas) in order to test its validity as a complementary resistance monitoring tool.
- pfmdr1 Plasmodium falciparum :
multidrug resistance gene 1
- pfcrt Plasmodium falciparum :
chloroquine resistance transporter
- pfdhfr Plasmodium falciparum :
- pfdhps Plasmodium falciparum :
- pfATPase6 :
sarcoendoplasmic reticulum Ca2+-ATPase (SERCA) of P. falciparum
polymerase chain reaction-restriction fragment length polymorphism.
Olliaro P: Drug resistance hampers our capacity to roll back malaria. Clin Infect Dis. 2005, 41 (Suppl 4): S247-S257. 10.1086/430785.
Hastings IM, Korenromp EL, Bloland PB: The anatomy of a malaria disaster: drug policy choice and mortality in African children. Lancet Infect Dis. 2007, 7: 739-748. 10.1016/S1473-3099(07)70214-1.
Hastings IM: Molecular markers as indicators of antimalarial drug failure rates. Trop Med Int Health. 2007, 12: 1298-1301.
Sibley CH, Hyde JE, Sims PF, Plowe CV, Kublin JG, Mberu EK, Cowman AF, Winstanley PA, Watkins WM, Nzila AM: Pyrimethamine-sulfadoxine resistance in Plasmodium falciparum: what next?. Trends Parasitol. 2001, 17: 582-588. 10.1016/S1471-4922(01)02085-2.
Plowe CV, Roper C, Barnwell JW, Happi CT, Joshi HH, Mbacham W, Meshnick SR, Mugittu K, Naidoo I, Price RN, Shafer RW, Sibley CH, Sutherland CJ, Zimmerman PA, Rosenthal PJ: World Antimalarial Resistance Network (WARN) III: molecular markers for drug resistant malaria. Malar J. 2007, 6: 121-10.1186/1475-2875-6-121.
World Health Organization (WHO): Susceptibility of Plasmodium falciparum to antimalarial drugs. Report on global monitoring 1996-2004. WHO/HTM/MAL/2005.1103. 2005
Cooper RA, Hartwig CL, Ferdig MT: pfcrt is more than the Plasmodium falciparum chloroquine resistance gene: a functional and evolutionary perspective. Acta Trop. 2005, 94: 170-180.
Duraisingh MT, Cowman AF: Contribution of the pfmdr1 gene to antimalarial drug-resistance. Acta Trop. 2005, 94: 181-190.
Picot S, Olliaro P, de Monbrison F, Bienvenu AL, Price RN, Ringwald P: A systematic review and meta-analysis of evidence for correlation between molecular markers of parasite resistance and treatment outcome in falciparum malaria. Malar J. 2009, 8: 89-10.1186/1475-2875-8-89.
Gregson A, Plowe CV: Mechanisms of resistance of malaria parasites to antifolates. Pharmacol Rev. 2005, 57: 117-145. 10.1124/pr.57.1.4.
Eckstein-Ludwig U, Webb RJ, Van Goethem ID, East JM, Lee AG, Kimura M, O'Neill PM, Bray PG, Ward SA, Krishna S: Artemisinins target the SERCA of Plasmodium falciparum. Nature. 2003, 424: 957-961. 10.1038/nature01813.
Jambou R, Legrand E, Niang M, Khim N, Lim P, Volney B, Ekala MT, Bouchier C, Esterre P, Fandeur T, Mercereau-Puijalon O: Resistance of Plasmodium falciparum field isolates to in-vitro artemether and point mutations of the SERCA-type PfATPase6. Lancet. 2005, 366: 1960-1963. 10.1016/S0140-6736(05)67787-2.
Uhlemann AC, Cameron A, Eckstein-Ludwig U, Fischbarg J, Iserovich P, Zuniga FA, East M, Lee A, Brady L, Haynes RK, Krishna S: A single amino acid residue can determine the sensitivity of SERCAs to artemisinins. Nat Struct Mol Biol. 2005, 12: 628-629. 10.1038/nsmb947.
Lozovsky ER, Chookajorn T, Brown KM, Imwong M, Shaw PJ, Kamchonwongpaisan S, Neafsey DE, Weinreich DM, Hartl DL: Stepwise acquisition of pyrimethamine resistance in the malaria parasite. Proc Natl Acad Sci USA. 2009, 106: 12025-12030. 10.1073/pnas.0905922106.
Adagut IS, Warhurst DC: Plasmodium falciparum: linkage disequilibrium between loci in chromosomes 7 and 5 and chloroquine selective pressure in Northern Nigeria. Parasitology. 2001, 123: 219-224.
Happi CT, Gbotosho GO, Folarin OA, Sowunmi A, Bolaji OM, Fateye BA, Kyle DE, Milhous W, Wirth DF, Oduola AM: Linkage disequilibrium between two distinct loci in chromosomes 5 and 7 of Plasmodium falciparum and in vivo chloroquine resistance in Southwest Nigeria. Parasitol Res. 2006, 100 (1): 141-8. 10.1007/s00436-006-0246-4.
Mita T, Kaneko A, Hombhanje F, Hwaihwanje I, Takahashi N, Osawa H, Tsukahara T, Masta A, Lum JK, Kobayakawa T, Ishizaki T, Björkman A: Role of pfmdr1 mutations on chloroquine resistance in Plasmodium falciparum isolates with pfcrt K76T from Papua New Guinea. Acta Trop. 2006, 98: 137-144. 10.1016/j.actatropica.2006.03.002.
Mu J, Ferdig MT, Feng X, Joy DA, Duan J, Furuya T, Subramanian G, Aravind L, Cooper RA, Wootton JC, Xiong M, Su XZ: Multiple transporters associated with malaria parasite responses to chloroquine and quinine. Mol Microbiol. 2003, 49: 977-989. 10.1046/j.1365-2958.2003.03627.x.
Djimdé A, Doumbo OK, Cortese JF, Kayentao K, Doumbo S, Diourte Y, Dicko A, Su XZ, Nomura T, Fidock DA, Wellems TE, Plowe CV, Coulibaly D: A molecular marker for chloroquine-resistant falciparum malaria. N Engl J Med. 2001, 344: 257-263. 10.1056/NEJM200101253440403.
Fevre EM, Barnish G: Malaria-treatment policies: when and how should they be changed?. Ann Trop Med Parasitol. 1999, 93: 549-560. 10.1080/00034989958050.
Talisuna AO, Kyosiimire-Lugemwa J, Langi P, Mutabingwa TK, Watkins W, Van Marck E, Egwang T, D'Alessandro U: Role of the pfcrt codon 76 mutation as a molecular marker for population-based surveillance of chloroquine (CQ)-resistant Plasmodium falciparum malaria in Ugandan sentinel sites with high CQ resistance. Trans R Soc Trop Med Hyg. 2002, 96: 551-556. 10.1016/S0035-9203(02)90439-6.
Talisuna AO, Langi P, Mutabingwa TK, Watkins W, Van Marck E, Egwang TG, D'Alessandro U: Population-based validation of dihydrofolate reductase gene mutations for the prediction of sulfadoxine-pyrimethamine resistance in Uganda. Trans R Soc Trop Med Hyg. 2003, 97: 338-342. 10.1016/S0035-9203(03)90163-5.
Muller I, Bockarie M, Alpers M, Smith T: The epidemiology of malaria in Papua New Guinea. Trends Parasitol. 2003, 19: 253-259. 10.1016/S1471-4922(03)00091-6.
Benet A, Mai A, Bockarie F, Lagog M, Zimmerman P, Alpers MP, Reeder JC, Bockarie MJ: Polymerase chain reaction diagnosis and the changing pattern of vector ecology and malaria transmission dynamics in Papua New Guinea. Am J Trop Med Hyg. 2004, 71: 277-284.
Genton B, Al Yaman F, Beck HP, Hii J, Mellor S, Narara A, Gibson N, Smith T, Alpers MP: The epidemiology of malaria in the Wosera area, East Sepik Province, Papua New Guinea, in preparation for vaccine trials. I. Malariometric indices and immunity. Ann Trop Med Parasitol. 1995, 89: 359-376.
Hii JL, Smith T, Vounatsou P, Alexander N, Mai A, Ibam E, Alpers MP: Area effects of bednet use in a malaria-endemic area in Papua New Guinea. Trans R Soc Trop Med Hyg. 2001, 95: 7-13. 10.1016/S0035-9203(01)90315-3.
Mehlotra RK, Kasehagen LJ, Baisor M, Lorry K, Kazura JW, Bockarie MJ, Zimmerman PA: Malaria infections are randomly distributed in diverse holoendemic areas of Papua New Guinea. Am J Trop Med Hyg. 2002, 67: 555-562.
Mueller I, Kundi J, Bjorge S, Namuigi P, Saleu G, Riley ID, Reeder JC: The epidemiology of malaria in the Papua New Guinea highlands: 3. Simbu Province. PNG Med J. 2004, 47: 159-173.
World Health Organization (WHO): Assessment and monitoring of antimalarial drug efficacy for the treatment of uncomplicated falciparum malaria. WHO/HTM/RBM/2003.50. 2003
Marfurt J, Mueller I, Sie A, Maku P, Goroti M, Reeder JC, Beck HP, Genton B: Low efficacy of amodiaquine or chloroquine plus sulfadoxine-pyrimethamine against Plasmodium falciparum and P. vivax malaria in Papua New Guinea. Am J Trop Med Hyg. 2007, 77: 947-954.
Felger I, Beck HP: Genotyping of Plasmodium falciparum. PCR-RFLP analysis. Methods Mol Med. 2002, 72: 117-129.
Gilles HM, Warrell DA: Bruce-Chwatt's essential malariology. 1993, London, UK: Arnold (Hodder Headline Group)
Crameri A, Marfurt J, Mugittu K, Maire N, Regos A, Coppee JY, Sismeiro O, Burki R, Huber E, Laubscher D, Puijalon O, Genton B, Felger I, Beck HP: Rapid microarray-based method for monitoring of all currently known single-nucleotide polymorphisms associated with parasite resistance to antimalaria drugs. J Clin Microbiol. 2007, 45: 3685-3691. 10.1128/JCM.01178-07.
Marfurt J, Muller I, Sie A, Oa O, Reeder JC, Smith TA, Beck HP, Genton B: The usefulness of twenty-four molecular markers in predicting treatment outcome with combination therapy of amodiaquine plus sulphadoxine-pyrimethamine against falciparum malaria in Papua New Guinea. Malar J. 2008, 7: 61-10.1186/1475-2875-7-61.
Schneider AG, Premji Z, Felger I, Smith T, Abdulla S, Beck HP, Mshinda H: A point mutation in codon 76 of pfcrt of P. falciparum is positively selected for by Chloroquine treatment in Tanzania. Infect Genet Evol. 2002, 1: 183-189. 10.1016/S1567-1348(01)00021-1.
Nei M, Kumar S: Molecular Evolution and Phylogenetics. 2002, Oxford University Press, New York
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.
Chen N, Russell B, Staley J, Kotecka B, Nasveld P, Cheng Q: Sequence polymorphisms in pfcrt are strongly associated with chloroquine resistance in Plasmodium falciparum. J Infect Dis. 2001, 183: 1543-1545. 10.1086/320206.
Mehlotra RK, Mattera G, Bhatia K, Reeder JC, Stoneking M, Zimmerman PA: Insight into the early spread of chloroquine-resistant Plasmodium falciparum infections in Papua New Guinea. J Infect Dis. 2005, 192: 2174-2179. 10.1086/497694.
Nagesha HS, Casey GJ, Rieckmann KH, Fryauff DJ, Laksana BS, Reeder JC, Maguire JD, Baird JK: New haplotypes of the Plasmodium falciparum chloroquine resistance transporter (pfcrt) gene among chloroquine-resistant parasite isolates. Am J Trop Med Hyg. 2003, 68: 398-402.
Mita T, Kaneko A, Hwaihwanje I, Tsukahara T, Takahashi N, Osawa H, Tanabe K, Kobayakawa T, Bjorkman A: Rapid selection of dhfr mutant allele in Plasmodium falciparum isolates after the introduction of sulfadoxine/pyrimethamine in combination with 4-aminoquinolines in Papua New Guinea. Infect Genet Evol. 2006, 6 (6): 447-52. 10.1016/j.meegid.2006.02.004. Epub 2006 Apr 5
Dorsey G, Dokomajilar C, Kiggundu M, Staedke SG, Kamya MR, Rosenthal PJ: Principal role of dihydropteroate synthase mutations in mediating resistance to sulfadoxine-pyrimethamine in single-drug and combination therapy of uncomplicated malaria in Uganda. Am J Trop Med Hyg. 2004, 71: 758-763.
Berens N, Schwoebel B, Jordan S, Vanisaveth V, Phetsouvanh R, Christophel EM, Phompida S, Jelinek T: Plasmodium falciparum: correlation of in vivo resistance to chloroquine and antifolates with genetic polymorphisms in isolates from the south of Lao PDR. Trop Med Int Health. 2003, 8: 775-782. 10.1046/j.1365-3156.2003.01099.x.
Mugittu K, Ndejembi M, Malisa A, Lemnge M, Premji Z, Mwita A, Nkya W, Kataraihya J, Abdulla S, Beck HP, Mshinda H: Therapeutic efficacy of sulfadoxine-pyrimethamine and prevalence of resistance markers in Tanzania prior to revision of malaria treatment policy: Plasmodium falciparum dihydrofolate reductase and dihydropteroate synthase mutations in monitoring in vivo resistance. Am J Trop Med Hyg. 2004, 71: 696-702.
Mutabingwa T, Nzila A, Mberu E, Nduati E, Winstanley P, Hills E, Watkins W: Chlorproguanil-dapsone for treatment of drug-resistant falciparum malaria in Tanzania. Lancet. 2001, 358: 1218-1223. 10.1016/S0140-6736(01)06344-9.
Mockenhaupt FP, Teun BJ, Eggelte TA, Schreiber J, Ehrhardt S, Wassilew N, Otchwemah RN, Sauerwein RW, Bienzle U: Plasmodium falciparum dhfr but not dhps mutations associated with sulphadoxine-pyrimethamine treatment failure and gametocyte carriage in northern Ghana. Trop Med Int Health. 2005, 10: 901-908. 10.1111/j.1365-3156.2005.01471.x.
Djimde A, Doumbo OK, Steketee RW, Plowe CV: Application of a molecular marker for surveillance of chloroquine-resistant falciparum malaria. Lancet. 2001, 358: 890-891. 10.1016/S0140-6736(01)06040-8.
Kublin JG, Dzinjalamala FK, Kamwendo DD, Malkin EM, Cortese JF, Martino LM, Mukadam RA, Rogerson SJ, Lescano AG, Molyneux ME, Winstanley PA, Chimpeni P, Taylor TE, Plowe CV: Molecular markers for failure of sulfadoxine-pyrimethamine and chlorproguanil-dapsone treatment of Plasmodium falciparum malaria. J Infect Dis. 2002, 185: 380-388. 10.1086/338566.
Kyabayinze D, Cattamanchi A, Kamya MR, Rosenthal PJ, Dorsey G: Validation of a simplified method for using molecular markers to predict sulfadoxine-pyrimethamine treatment failure in African children with falciparum malaria. Am J Trop Med Hyg. 2003, 69: 247-252.
Sidhu AB, Valderramos SG, Fidock DA: pfmdr1 mutations contribute to quinine resistance and enhance mefloquine and artemisinin sensitivity in Plasmodium falciparum. Mol Microbiol. 2005, 57: 913-926. 10.1111/j.1365-2958.2005.04729.x.
McMahon JE: Malaria endemicity amongst the semi-nomadic people of the Karimui area in Papua New Guinea. P N G Med J. 1974, 17: 99-107.
Spencer M: The history of malaria control in the southwest Pacific region, with particular reference to Papua New Guinea and the Solomon Islands. P N G Med J. 1992, 35: 33-66.
Talisuna AO, Langi P, Bakyaita N, Egwang T, Mutabingwa TK, Watkins W, Van Marck E, D'Alessandro U: Intensity of malaria transmission, antimalarial-drug use and resistance in Uganda: what is the relationship between these three factors?. Trans R Soc Trop Med Hyg. 2002, 96: 310-317. 10.1016/S0035-9203(02)90108-2.
Talisuna AO, Erhart A, Samarasinghe S, Van Overmeir C, Speybroeck N, D'Alessandro U: Malaria transmission intensity and the rate of spread of chloroquine resistant Plasmodium falciparum: Why have theoretical models generated conflicting results?. Infect Genet Evol. 2006, 6: 241-248. 10.1016/j.meegid.2005.06.003.
Laufer MK, Plowe CV: Withdrawing antimalarial drugs: impact on parasite resistance and implications for malaria treatment policies. Drug Resist Updat. 2004, 7: 279-288. 10.1016/j.drup.2004.08.003.
World Health Organization (WHO): Facts on ACTs (artemisinin-based combination therapies). 2006, [http://rbm.who.int/cmc_upload/0/000/015/364/RBMInfosheet_9.htm]
Sibley CH, Barnes KI, Watkins WM, Plowe CV: A network to monitor antimalarial drug resistance: a plan for moving forward. Trends Parasitol. 2008, 24: 43-48. 10.1016/j.pt.2007.09.008.
We thank all the study participants and their parents or legal guardians. We also thank the field and laboratory staff at the Papua New Guinea Institute of Medical Research branches in Goroka, Maprik and Madang, for their assistance with patient recruitment and microscopical analysis. Jean-Yves Coppee and Odile Sismeiro of the Plate-forme puces à ADN, Institut Pasteur, Paris, France provided technical support for DNA microarray analysis.
Financial Support: This study was supported by the Swiss National Science Foundation (Grant no: 3100-AO-103968).
The authors declare that they have no competing interests.
JM participated in the coordination of the field and laboratory studies, performed data acquisition and molecular and statistical analyses and drafted the manuscript. TAS, IMH, HPB and BG participated in the design of the study, the statistical analysis and the drafting of the manuscript. IM, AS, OO, MB and JCR participated in the coordination of the field studies. All authors read and approved the final manuscript.