|Home | About | Journals | Submit | Contact Us | Français|
Treatment of Plasmodium falciparum is complicated by the emergence and spread of parasite resistance to many of the first-line drugs used to treat malaria. Antimalarial drug resistance has been associated with specific point mutations in several genes, suggesting that these single nucleotide polymorphisms can be useful in tracking the emergence of drug resistance. In India, P. falciparum infection can manifest itself as asymptomatic, mild, or severe malaria, with or without cerebral involvement. We tested whether chloroquine- and antifolate drug-resistant genotypes would be more commonly associated with cases of cerebral malaria than with cases of mild malaria in the province of Jabalpur, India, by genotyping the dhps, dhfr, pfmdr-1, and pfcrt genes using pyrosequencing, direct sequencing, and real-time PCR. Further, we used microsatellites surrounding the genes to determine the origins and spread of the drug-resistant genotypes in this area. Resistance to chloroquine was essentially fixed, with 95% of the isolates harboring the pfcrt K76T mutation. Resistant genotypes of dhfr, dhps, and pfmdr-1 were found in 94%, 17%, and 77% of the isolates, respectively. Drug-resistant genotypes were equally likely to be associated with cerebral malaria as with mild malaria. We found evidence of a selective sweep in pfcrt and, to a lesser degree, in dhfr, indicating high levels of resistance to chloroquine and evolving resistance to pyrimethamine. Microsatellites surrounding pfcrt indicate that the resistant genotypes (SVMNT) were most similar to those found in Papua New Guinea.
Malaria is arguably the most important vector-borne disease in the world, with annual morbidity and mortality estimates surpassing 300 and 1 million, respectively (20, 61). Over 90% of the total malaria incidence is reported from sub-Saharan and tropical Africa; however, each year Southeast Asia, including the Indian subcontinent, reports approximately 2.5 million malaria cases, 75% of which are from India (20). Further, Plasmodium falciparum incidence in India has increased dramatically over the past few years, including the spread of drug-resistant strains (1, 2, 23, 46, 47).
Efforts to control malaria have been hindered by the rapid rise and spread of drug-resistant P. falciparum strains. Chloroquine (CQ)-resistant strains of P. falciparum first appeared in the late 1950s, almost simultaneously in Southeast Asia and South America (51, 58, 64), and subsequently spread through most regions where P. falciparum is endemic. Sulfadoxine-pyrimethamine (SP) was next used as the drug of choice against CQ-resistant malaria; however, resistance quickly emerged on the Thai-Cambodian border around 1980 and is now found throughout most of Southeast Asia, the Amazonian basin of South America, and Africa (1, 2, 7, 17, 40, 53). In India, CQ and SP resistance was first documented in 1973 (42) and 1979 (10), respectively, in the northeast region of the country. Now, studies using molecular markers suggest that CQ and SP resistance is widespread across India (1, 23, 55). However, in India CQ still remains the first line of treatment for Plasmodium vivax malaria and for P. falciparum in low-risk and CQ-sensitive areas. In light of reports of CQ treatment failures, artesunate plus SP (artesunate combination therapy [ACT]) has been introduced in states with high burdens of P. falciparum malaria (6, 45) and is being implemented for other districts with high prevalences of P. falciparum.
Resistance to chloroquine has been associated with point mutations in the P. falciparum chloroquine resistance transporter (pfcrt) gene (16), while resistance to SP has been linked to the dihydrofolate reductase (dhfr) and dihydropteroate synthase (dhps) genes (32, 52). Point mutations in P. falciparum multidrug resistance gene 1 (pfmdr-1) have been reported to modulate resistance to different antimalarial drugs, and variations in copy number appear to be associated with mefloquine resistance (12, 36).
In order to combat drug resistance in Plasmodium falciparum, it is important to understand the genetic basis and the evolutionary forces affecting loci governing resistance. The discovery of new drug targets and the development of effective drugs and vaccines require careful study of the population genetics of P. falciparum. Investigations regarding point mutations in genes conferring drug resistance and the microsatellite loci that surround these genes can provide information on selection pressures, rates of recombination, and the potential origin of the resistant alleles or mutations.
P. falciparum infection can manifest as asymptomatic, mild (uncomplicated) malaria or severe malaria, with or without cerebral involvement; little is known about the factors involved in these clinical manifestations. Cerebral malaria (CM) is one of the most common complications of P. falciparum infection in India, besides severe malaria anemia and multiorgan failure (26, 27, 60, 61). It is not known whether drug-resistant parasites also contribute to the increased risk for CM, especially because patients may receive inadequate treatment with drugs of reduced efficacy. In this context, we were interested in determining whether parasites with resistant genotypes were more often associated with patients diagnosed with CM than with patients with mild malaria (MM). We hypothesized that individuals harboring resistant parasites may be more likely to progress to severe disease due to treatment failure than those with wild-type parasites. Additionally, we wanted to determine whether drug-resistant genotypes in India have evolved locally or have been influenced by gene flow from other regions. To this end, we genotyped four genes associated with drug resistance (pfcrt, dhfr, dhps, and pfmdr-1) and assessed the genetic diversity of microsatellites surrounding pfcrt, dhfr, and dhps from P. falciparum-positive blood samples taken from patients enrolled in a hospital-based study to assess neurological disorders associated with cerebral malaria in central India.
This study utilized samples collected from a hospital study to assess the neurological deficits associated with CM patients in India. The study details have been reported previously (18, 19). Here we provide brief details about the study site and study populations. The study was conducted in Jabalpur province of Madhya Pradesh state, located in central India. The study was carried out at two sites: the Nethaji Subhash Chandra Bose (NSCB) Hospital (a regional referral hospital) in Jabalpur and the Civil Hospital (a primary care hospital) in Maihar, Satna district. Both P. vivax and P. falciparum are prevalent in this area, and P. falciparum transmission occurs primarily during the monsoon and postmonsoon seasons (July to January). Previous studies revealed that malaria is present in all age groups in this region, with the highest prevalence occurring in children between 8 and 14 years of age (48).
The study samples used in this investigation were obtained from all patients enrolled between October 2004 and December 2006, as reported by Jain et al. (19). This included 52 cerebral malaria patients and 52 mild malaria patients. All subjects were enrolled after providing written informed consent. Patients were enrolled as they presented to one of the two hospitals with either mild malaria or severe malaria that progressed to cerebral malaria. The human experimentation guidelines of the National Institutes of Health (NIH) and of the ethical committees of the Morehouse School of Medicine, the National Institute of Malaria Research (India), and the Centers for Disease Control and Prevention (CDC) were followed.
To be considered to have a case of CM, a patient had to fulfill the World Health Organization's definition of CM (59), have a Glasgow coma score of ≤10, have P. falciparum parasitemia, and have no other clinically evident cause of impaired consciousness.
Fifty-two patients who had fever with P. falciparum parasitemia of <25,000 parasites/μl of blood (detected microscopically from blood smears); no evidence of impaired consciousness or seizures; and no history of mental illness, meningitis, or accidental head injury were included.
Venous blood samples from children (2 to 5 ml) and adults (10 ml) were collected soon after enrollment from patients before the start of antimalarial treatment or transfusions. Patients were required to complete a questionnaire detailing any drug treatment history. The blood was separated after centrifugation in Becton Dickinson cell preparation tubes (catalogue no. 362753; BD Pharmingen, Franklin Lakes, NJ); red cell pellets and plasma were separated, aliquoted, and frozen at −80°C for long-term storage.
DNA was isolated from red cell pellets by using the QIAamp DNA minikit (Qiagen, Valencia, CA) according to the manufacturer's recommendations. To confirm the diagnosis and to rule out any mixed infection with other species of malaria parasites, we tested all the samples using a nested diagnostic PCR as previously described (50) with Promega (Madison, WI) Taq PCR master mix according to the manufacturer's instructions. All PCRs were conducted in a Bio-Rad iCycler thermocycler. All the samples were confirmed to be P. falciparum, and no other malaria parasite species was detected.
Mutations in dhfr and dhps were determined through pyrosequencing as previously described (65). Briefly, primary and nested PCR amplicons were generated in 50-μl reaction volumes containing 1 μl of DNA, 0.5 μM forward and reverse primers, and 25 μl Promega Taq PCR master mix; the cycling conditions are given in reference 65. Single-stranded biotinylated PCR products were prepared for pyrosequencing as described by Zhou et al. in 2006 (65), and reactions were performed according to the manufacturer's instructions using the PSQ 96 single nucleotide polymorphism (SNP) reagent kit (Biotage AB). Genotypes were determined using SNP software (Biotage AB).
Mutations in pfcrt were determined through direct sequencing of codons 72 to 76, 220, and 271 using primers and PCR parameters described previously (16). PCR products were purified using Montage PCR centrifugal filter devices (Millipore Corporation, Billerica, MA) according to the manufacturer's instructions. Purified PCR products were sequenced using the BigDye Terminator Cycle Sequencing Ready Reaction kit (Applied Biosystems, Foster City, CA) according to the manufacturer's recommended protocol. Unincorporated fluorescence-labeled deoxynucleoside triphosphates were removed with Centri-Sep columns (Princeton Separations, Adelphia, NJ) according to the manufacturer's recommendations. Both strands were sequenced using an ABI 3130xl genetic analyzer (Applied Biosystems, Foster City, CA). Nucleotide sequences were edited and assembled with the Pregap and Gap programs of the STADEN sequence analysis package (51a). Multiple sequence alignments were made with the PILEUP program of the GCG package (Wisconsin package, version 10.3; Accelrys [GCG], San Diego, CA).
Mutations in codons 86 and 184 of the pfmdr-1 gene, located on chromosome 5, were detected using a Stratagene MX3005P real-time PCR system as previously described (37). Briefly, 2 μl of DNA was used as a template for real-time PCR using 1× Platinum quantitative PCR (qPCR) master mix (Invitrogen, Carlsbad, CA), 300 nM forward and reverse primers, and 100 nM each probe. TaqMan MGB probes (Applied Biosystems, Foster City, CA) were used for both codons. Wild-type probes were labeled with a FAM (5-carboxyfluorescein) reporter dye and MGB-NFQ (minor groove binder-nonfluorescent quencher), while mutant probes were labeled with a VIC reporter dye (Applied Biosystems, Foster City, CA) and MGB-NFQ. The following cycling conditions were used for PCR of codons 86 and 184: 96°C for 10 min, followed by 50 cycles of 96°C for 15 s and either 59°C for 1 min (codon 86) or 61°C for 1 min (codon 184). Data were collected at the end of each annealing step.
Samples were assayed for nine microsatellites that span 25 kb around dhfr on chromosome 4, seven microsatellites that span approximately 9 kb around dhps on chromosome 8, and seven microsatellites that span 31 kb around pfcrt on chromosome 7. Microsatellite PCR primer sequences for dhfr, dhps, and pfcrt are provided as supplemental material and are adopted from the work of Nair et al. (24), Roper et al. (41), Pearce et al. (31), and McCollum et al. (21). The distance from each gene was recalculated using the current version of PlasmoDB (version 5.5). Single-round PCR thermal cycling conditions are detailed by Nair et al. (24), and nested PCR thermal cycling conditions are given by Roper et al. (41). PCR products were separated on an ABI 3130xl genetic analyzer and were scored using GeneMapper software, version 3.7 (Applied Biosystems, Foster City, CA). Samples that were found by genotyping to be multiply infected were not considered in the analysis (n = 2), and no samples generated multiple alleles for any of the loci used for haplotype determination. Missing data (no amplification) were reported but not considered in defining haplotypes. Microsatellites surrounding pfcrt were calibrated to the sizes reported by Wootton et al. (62) by running the microsatellite panel on both Dd2 and 7G8 and adjusting the sizes of the loci by the difference between our sequencer's results and those previously reported for those two isolates.
An exact test of linkage disequilibrium (38) was used to test for significant association between microsatellite markers around pfcrt, dhfr, and dhps using Arlequin, version 3.11 (14) with 10,000 Monte Carlo steps. The null hypothesis is no linkage disequilibrium between loci located on different chromosomes in a panmictic population with no genetic structure.
The genetic diversity for each microsatellite locus was measured by calculating the expected heterozygosity (He) as [n/(n − 1)][1 − Σpi2], where n is the number of isolates sampled and pi is the frequency of the ith allele. The sampling variance for He was calculated as 2(n − 1)/n3[2(n − 2)][Σpi3 − (Σpi2)2] (25). The Excel microsatellite tool kit was used to compute allele frequencies and the number of alleles per locus (A) and to format the data for Arlequin (30). Chi-square tests were used to test for differences between the numbers of mutant and wild-type CM and MM cases at each locus (JMP version 8.0, 2008; SAS Institute, Inc., Cary, NC), without correction for multiple comparisons.
The genetic relationship between the SVMNT genotypes found in India and the SVMNT, CVIET, CVMNT, and CVMET genotypes from other parts of the world (62) was investigated by a median-joining network constructed using 10-locus haplotypes in NETWORK, version 18.104.22.168 (http://www.fluxus-engineering.com/sharenet.htm). Median-joining networks are used for reconstructing the phylogeny of regions with reticulate evolution (5, 33). We constructed the network using all Indian haplotypes and excluding the singletons. For simplicity, only the analysis including Indian haplotypes that occurred more than once is given.
The distributions of drug-resistant genotypes in pfcrt, dhfr, dhps, and pfmdr-1 in CM and MM patients are reported in Fig. Fig.1.1. We found no significant differences in the frequencies of dhfr (P = 0.36), pfcrt (P = 0.65), or pfmdr-1 codon 86 (P = 0.37) or codon 184 (P = 0.12) genotypes between the CM and MM groups. Mutations in dhps were more likely to be found in patients with mild malaria than in those with cerebral malaria (P = 0.04), though this may be due to chance. Only one sample, from a CM patient, was found to be wild type (susceptible) at pfcrt, dhfr, dhps, and pfmdr-1. Most of the samples had mutations at two or more loci associated with drug resistance. In fact, 80 samples (77%) were found to have mutations in either three or all four of the resistance loci (Fig. (Fig.2).2). Ninety-nine (95%) of the samples contained the K76T mutation for chloroquine resistance in combination with at least one other mutation in dhfr, dhps, or pfmdr-1.
A large proportion of P. falciparum parasites in this Indian population harbor the mutation for chloroquine resistance. The CQ resistant genotype SVMNT was found in 95% (n = 104) of the samples (Fig. (Fig.1A;1A; Table Table1).1). The remaining five were of the chloroquine-sensitive ancestral genotype CVMNK. Three of the five chloroquine-sensitive genotypes were from CM patients, and two of the sensitive genotypes had the double dhfr mutation. Additionally, we found 100 parasites with the A220S mutation; three CVMNK parasites and one SVMNT parasite retained the ancestral state (Table (Table1).1). None of the samples harbored the Q271E mutation (Table (Table11).
Mutations encoding pyrimethamine resistance were found in 94.2% of the samples (Fig. (Fig.1B).1B). Twenty-two samples (21%) were mutant at one amino acid position, 59R or 108N, while the remaining 76 patients (73%) had the 59R 108N double mutant. Importantly none of the samples contained the highly resistant triple or quadruple mutant genotype in this region. Six patients had parasite isolates with the ancestral pyrimethamine-sensitive dhfr genotype: one had drug-sensitive genotypes at all four genes tested; three had mutations at the pfcrt and pfmdr-1 loci; one had mutations at the dhps, pfcrt, and pfmdr-1 loci; and one had mutations at the dhps and pfcrt loci (Fig. (Fig.22).
Parasites with mutations associated with sulfadoxine resistance in dhps were found in 18 (17.3%) of the samples, 14 of which were mutant at only one amino acid position (Fig. (Fig.1C).1C). A single sample was mutant at three amino acid positions (437G, 540E, and 581G), while three others were mutated at two (437G and 581G, 437G and 540E, or 540E and 581G).
Mutations in pfmdr-1 were found in 14 samples (16%) at codon 86 and in 74 samples (71%) at codon 184 (Fig. (Fig.1D).1D). We found the presence of both wild-type and mutant alleles in two samples, one from a CM patient with both alleles at codon 86 and one from an MM patient with both alleles at codon 184. Forty-six percent of the samples in the real-time PCR assay for codon 184 failed to generate data, requiring us to sequence the samples. We found an additional silent mutation at codon 182, located within the probe binding site and explaining the failure of the assay to produce results for these samples (data not shown).
The molecular signatures of new, strongly selected advantageous mutations, such as genes that confer drug resistance on the malaria parasite, are distinct. In a process known as genetic hitchhiking, flanking nucleotides are linked by proximity to the allele (mutation) under selection (49) and these nucleotides increase in frequency as the advantageous allele increases in frequency. A selective sweep occurs when a mutation becomes fixed in the population, resulting in a decrease in the level of variation in surrounding nucleotides. Increases in the levels of linkage disequilibrium will also be found in a transient phase of a selective sweep (28).
Seven microsatellites spanning 35 kb around pfcrt were characterized in order to assess the relationship among pfcrt alleles reported from other areas of the world as well as to assess the impact of selection on this locus. The specific microsatellite profiles surrounding the chloroquine resistance transporter were compared to those found from other parts of the world as reported by Wootton et al. (62) (Table (Table1).1). Patterns of allele sizes were most similar to those found in Papua New Guinea (PNG). Indeed, 18% of the resistant alleles had allele sizes identical to those reported from PNG, and 21% of the resistant alleles differed at only one locus. Four of the five parasite isolates with the CVMNK genotype had a more diverse pattern of microsatellite allele sizes, differing from the PNG type at five or more loci (Table (Table1).1). The median-joining network indicated a close relationship between the SVMNT genotypes from India and PNG, with one hypothetical node between the most common Indian genotype and that from PNG (Fig. (Fig.3).3). Although in Fig. Fig.33 we have included data from haplotypes with more than one isolate, the median-joining network that included all Indian haplotypes was similar, and for the purpose of simplicity, we have not shown this figure. However, it is important that the relationships between the Indian/PNG cluster, the South American cluster, and the African-Asian cluster were retained (data not shown). Heterozygosity estimates were lowest around the pfcrt gene, indicating strong selection pressure (Fig. (Fig.4A).4A). We found a much greater reduction in heterozygosity around the SVMNT genotype than around the ancestral CVMNK genotype; however, there was a limited number of CVMNK genotypes for comparison (n = 5) (Fig. (Fig.4A).4A). Significant linkage disequilibrium was found between 7 microsatellite markers surrounding pfcrt (Fig. (Fig.55).
We also found a reduction in heterozygosity immediately surrounding dhfr along chromosome 4 (Fig. (Fig.4B).4B). Samples with resistant genotypes were less diverse, on average, than those with wild-type genotypes. Heterozygosity estimates ranged from 0.28 to 0.81 (mean He, 0.638), and it is noteworthy that triple mutant dhfr genotypes have not become established in this region, suggesting the possible beginning stages of a sweep for pyrimethamine resistance. Significant linkage disequilibrium was found between 11 pairwise combinations of loci within dhfr and between 7 pairwise combinations of loci between dhfr and dhps (Fig. (Fig.55).
We again found reduced variation around dhps; however, there was no obvious difference in microsatellite variation between resistant and susceptible genotypes (Fig. (Fig.4C).4C). Heterozygosity estimates ranged from 0.186 to 0.993 (mean He, 0.797). Linkage disequilibrium within dhps was limited to six combinations of loci (Fig. (Fig.55).
One of the objectives of this study was to determine whether malaria parasites in central India with mutations conferring drug resistance evolved locally or whether gene flow from other regions of the world where malaria is endemic has shaped the diversity of these parasites. Mutations conferring drug resistance were found in 95% (pfcrt), 94% (dhfr), 17% (dhps), and 77% (pfmdr-1) of the parasites. Allele sizes for microsatellites surrounding pfcrt were most closely related to those previously reported from PNG, and a median-joining network indicated a close relationship between PNG SVMNT parasites and those from central India. There is strong selection on pfcrt based on the heterozygosity estimates of the surrounding microsatellites. Similarly, we found reduced variation surrounding dhfr and dhps, albeit to a much lesser extent. Additionally, we investigated whether drug-resistant parasites were more likely to be associated with patients diagnosed with CM than with patients diagnosed with mild malaria; however, no association was found.
We used direct sequencing for pfcrt, pyrosequencing for dhfr and dhps, and real-time PCR for pfmdr-1 genotyping. While three different genotyping methods were used in this study, each at different time points, this difference is unlikely to have affected the results of this study. Pyrosequencing has been used to genotype P. falciparum parasites previously and was shown to be more cost-effective, less time-consuming, and more efficient at detecting mixed parasite infections than sequencing (65). We have previously shown that pyrosequencing does not bias the data compared to direct sequencing (65). Similarly, real-time PCR has the advantage of detecting mixed infections and, in our case, led to the discovery of a novel, silent mutation not previously reported. Direct sequencing is the most cost-effective method for genotyping pfcrt, since the critical mutations occur in a 15-bp region. Pyrosequencing and real-time PCR were the most cost-effective and time-effective for genotyping dhfr, dhps, and pfmdr-1, since the mutations do not occur adjacent to one another.
CM patients were no more likely to have mutations in pfcrt, dhps, dhfr, or pfmdr-1 than patients exhibiting mild malaria. This suggests that, at least in this hospital study, the occurrence of CM is not associated with the overrepresentation of parasites with drug-resistant mutations in P. falciparum. One limitation of this study is the fact that the CQ-resistant pfcrt genotype is fixed in this population; therefore, it may be difficult to discern to what extent the resistant genotypes could have influenced the CM outcome in this genetic background. However, resistant genotypes of dhfr and dhps were not fixed, and no differences were found in the probability of infection with dhfr or pfmdr-1 resistant parasites. Thus, individuals with resistant parasites were no more likely to progress to severe malaria than individuals without resistant parasites.
All but five of the individuals genotyped for pfcrt had the K76T mutation, which is critical for resistance to chloroquine. Resistance to CQ was first detected in the early 1970s and is currently found throughout the country, though the highly resistant CVIET genotype is found only in the northeast and southeast regions of India (23). Interestingly, quite high frequencies of the CQ-resistant genotype were found in samples from India dating back to 1996, and these reports indicated the presence of SVMNT and CVIET genotypes in some states; however, it appears that the CVIET genotypes have not spread to the central Indian region (54). Of note, one resistant isolate (SVMNT) did not harbor the A220S mutation commonly associated with CQ-resistant alleles. While rare, the ancestral state at codon 220 has been found in other isolates with the K76T mutation in India, China, and the Philippines (8, 23, 63). We also found two isolates that were wild type at codon 76 but harbored the A220S mutation, which could potentially be due to a recombination event between wild-type and CQ-resistant parasites. The microsatellite loci surrounding pfcrt (mean He, 0.492) have much less variation than those found around dhfr (mean He, 0.638) and dhps (mean He, 0.797), and the selection valley of reduced variation around pfcrt is also wider than we find around dhfr and dhps resistance genotypes, suggestive of longer and stronger selective pressure on pfcrt, as reported previously (54). The specific microsatellite profiles surrounding the chloroquine resistance transporter were compared to those found in other parts of the world as reported by Wootton et al. (62) (Table (Table1).1). The microsatellite allele sizes are most similar to those found in PNG over 21 kb surrounding pfcrt. Most of the differences between Indian and PNG parasites occurred at kb −12.3, 5.97, and 18.8, and the majority of Indian parasites differed only by a few base pairs from the PNG parasites for loci within 5 kb of pfcrt. It is possible that the loci more than 5 kb from pfcrt are behaving neutrally, given that the heterozygosity has increased over 6-fold at these loci (Fig. (Fig.4A).4A). The median-joining diagram also suggested a close relationship between the SVMNT haplotypes in India and PNG, suggesting a common origin for these haplotypes. Similarities in microsatellite haplotypes from India and PNG have been observed previously (9). Although one can speculate, based on this study, that this genotype may have been introduced from PNG, it is difficult to confirm such a hypothesis without studying the parasite isolates from other parts of India and the neighboring region, and sampling in a population-based manner. In this context, it should be noted that, unlike the pattern in Africa, where the CVIET genotype with origins in Southeast Asia is widely prevalent, the CQ-resistant genotype found in central India resembles a South American or PNG pattern, where the SVMNT genotype is more commonly found. However, the SVMNT genotype in South America has evolved independently and is distinct from the PNG and Indian parasite types (62).
The pfmdr-1 gene is speculated to confer resistance to chloroquine, quinine, and mefloquine. Mutations at codon 86 have been associated with CQ resistance (4, 11, 34, 35, 56), and codons 184, 1034, 1042, and 1246 have been implicated to various degrees in resistance to mefloquine and artesunate (13, 15, 39). It has been suggested that mutations at codon 86 result in an increase in sensitivity to mefloquine, while codons 1034, 1042, and 1246 may confer resistance (13, 35, 39). A study from Cambodia found an increased association between artesunate-mefloquine failure and a mutation at codon 184 (43). In our study, a higher proportion of individuals had mutations at codon 184; however, it is difficult to associate this with mefloquine resistance, since this drug is not commonly used in India. A correlation between CQ resistance and the pfmdr-1 N86Y mutation has been found in Mali (11) and PNG (22), and it was suggested that the pfmdr-1 mutation in conjunction with the pfcrt K76T mutation yields enhanced levels of resistance to CQ (57). We found few parasites with pfmdr-1 mutations at codon 86. Similarly, Vathsala et al. (54) failed to find coselection of the pfcrt K76T mutation with the pfmdr-1 N86Y mutation in Indian isolates, and they suggested that the hypothesis that pfmdr-1 N86Y mutation confers a compensatory advantage for coping with CQ pressure may not be valid in all geographic regions.
In contrast to the near-fixation of CQ-resistant genotypes in central India, it appears that there has been recent strong selection for mutations in dhfr associated with pyrimethamine resistance, as indicated by the narrow selective valley seen in Fig. Fig.4B.4B. Additionally only narrow chromosomal regions (within 1 kb) are affected by hitchhiking around dhfr. Only 11 individuals lacked the critical mutation at codon 108 in dhfr that confers resistance to pyrimethamine, suggestive of strong selection for the 108N mutant allele. Both wild-type and single mutant dhps genotypes exist in the population, with no double or triple mutants. No large differences in heterozygosity were seen between wild-type and resistant genotypes in dhps, while in dhfr there is a noticeable reduction, as would be expected after strong selective pressure. These data, combined with the low frequency of mutations in dhps, suggest that selection is operating only weakly on dhps, or that there has been insufficient time for selection to leave a molecular signature. Evidence of selection is generally found in dhfr earlier than in dhps (29, 44); thus, it is likely that we are seeing the beginning stages of the evolution of SP resistance in the state of Madhya Pradesh. Slightly higher levels of heterozygosity are found around the microsatellite markers surrounding dhps than around dhfr (mean He, 0.797 and 0.638, respectively). The increased reduction in heterozygosity in markers surrounding dhfr supports our theory of the beginning stages of a selective sweep for pyrimethamine resistance in India. Furthermore, the majority of the resistant dhfr genotypes are 59R 108N double mutants, rather than the highly resistant triple mutant found in Africa and Southeast Asia. It has been reported that parasites from other parts of India, particularly from the northeast state of Assam, carry a triple mutant dhfr genotype (59R 108N 164L) (2) and that isolates from the Andaman and Nicobar Islands carry a quadruple mutant dhfr genotype (51I 59R 108N 164L) (2). Other studies have also found predominantly double (59R 108N) dhfr mutants on the mainland of India (3); thus, it appears that highly resistant genotypes have not spread extensively to central India. We speculate that the valley of reduced heterozygosity surrounding these two genes will likely become wider as time progresses if SP selection pressure is maintained or increases.
India has adopted combination therapy using artesunate plus SP to treat P. falciparum malaria in high-burden states where resistance to CQ is confirmed. Our data indicate that resistance to SP may be in the early stages of evolution, and it remains to be determined whether the introduction of artesunate plus SP will slow down the evolution of SP-resistant genotypes. Therefore, continued molecular surveillance to monitor changes in dhfr and dhps mutant genotypes will provide warning signals for potential changes in the efficacy of SP. Because CQ is still used to treat P. falciparum as well as P. vivax infections, a decline in the frequency of CQ-resistant alleles in India may not happen soon.
In summary, our study has shown that there is no increased presence of drug-resistant genotypes in CM patients compared to MM patients. The CQ-resistant SVMNT genotype found in central India is most closely related to the SVMNT genotypes found in PNG, and this genotype is under strong selective pressure. Selection has only recently begun to alter the frequencies of the SP-resistant dhfr and dhps genotypes in this area of India. A higher proportion of isolates was found with mutations at codon 184 of pfmdr-1 rather than codon 86, and the importance of these mutations in this population has yet to be determined, since mefloquine is not commonly used in India. Continued molecular surveillance in India will provide useful information about evolving drug-resistant genotypes.
T.M.-H. was supported by the American Society for Microbiology, Coordinating Center for Infectious Diseases Postdoctoral Fellowship Program in Infectious Disease and Public Health Microbiology. Funding for this work was provided in part by the CDC Antimalarial Resistance Working Group, the Atlanta Research and Education Foundation, the VA Medical Center, and NIH/FIC grant R21TW006804-01, awarded to J.K.S.
We thank Naomi Lucchi, Sumiti Vinayak, and Md Tauqeer Alam for discussion and critical review of the manuscript and the two anonymous reviewers for valuable comments. We thank the study staff in India and the volunteers who participated in this study.
The findings and conclusions in this report are those of the authors and do not necessarily represent the official position of the Centers for Disease Control and Prevention.
Published ahead of print on 28 December 2009.
†Supplemental material for this article may be found at http://aac.asm.org/.