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 ). 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. ). 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. . 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.