Previous Article | Next Article ![]()
Antimicrobial Agents and Chemotherapy, March 2007, p. 1071-1077, Vol. 51, No. 3
0066-4804/07/$08.00+0 doi:10.1128/AAC.01186-06
Copyright © 2007, American Society for Microbiology. All Rights Reserved.

Department of International Affairs and Tropical Medicine, Tokyo Women's Medical University, Tokyo, Japan,1 Laboratory of Malariology, Research Institute for Microbial Diseases, Osaka University, Osaka, Japan,2 The National Center for Parasitology, Entomology and Malaria Control, Phnom Penh, Cambodia,3 Department of Parasitology, National Institute of Infectious Diseases, Tokyo, Japan,4 Department of Biomedical Chemistry, Graduate School of Medicine, University of Tokyo, Tokyo, Japan,5 The Asian Centre of International Parasite Control, Faculty of Tropical Medicine, Mahidol University, Thailand,6 Department of Medicine, Malaria Research Laboratory, Karolinska Institutet, Stockholm, Sweden7
Received 22 September 2006/ Returned for modification 1 December 2006/ Accepted 18 December 2006
|
|
|---|
|
|
|---|
Pyrimethamine and sulfadoxine inhibit two separate enzymes in the folate synthesis pathway of Plasmodium falciparum: dihydrofolate reductase (DHFR) and dihydropteroate synthase (DHPS), respectively. Point mutations at amino acid positions 16, 50, 51, 59, 108, and 164 in the DHFR gene (dhfr) are the major causes of resistance to pyrimethamine (3, 17, 18, 24). The mutation at position 108 (Ser
Asn) appears to be an initial prerequisite for a significant (10-fold) increase in in vitro resistance (24). Additional mutations at other amino acid positions within the gene are associated with stepwise increases in resistance. Isolates harboring four mutations at positions 51, 59, 108, and 164 (CIRNL at positions 50, 51, 59, 108, and 164 [mutations are indicated by underlining]) show the highest pyrimethamine resistance so far described.
Various dhfr alleles have been observed in regions of endemicity (30). A dhfr triple mutant (CIRNI) represents the most common type in Africa and Southeast Asia, while the dhfr quadruple mutant (CIRNL) is observed predominantly in Thailand and some other regions in Southeast Asia where SP resistance is very high (1, 12, 30). Two distinct triple dhfr mutant genotypes (RICNI and CICNL) are prevalent in South America (2, 18). A five-amino-acid insertion after position 30, termed the Bolivia repeat, is also exclusive to South America, suggesting two unique and different evolutionary origins of pyrimethamine resistance in South America (2).
The migration of drug-resistant alleles can be traced by the analysis of microsatellite markers closely linked to the gene conferring resistance. Microsatellite analysis flanking pfcrt has revealed that chloroquine resistance evolved independently in at least four different regions: Southeast Asia, two regions in South America, and New Guinea (31). Meanwhile, all dhfr triple (CIRNI) and quadruple (CIRNL) mutants from Southeast countries displayed the same or nearly identical microsatellite haplotypes flanking dhfr (12). Strikingly, pyrimethamine-resistant isolates in Africa also harbored microsatellite haplotypes identical to those found in Southeast Asia (21), suggesting a single origin of pyrimethamine resistance in Southeast Asia, which subsequently spread to Africa. However, whether the Melanesian dhfr mutants originated in Southeast Asia or arose independently remains unclear.
In the present study, we determined dhfr and microsatellite haplotypes flanking the gene in P. falciparum isolates from Melanesia (Papua New Guinea, Vanuatu, and the Solomon Islands) and Southeast Asia (Thailand and Cambodia). Our results show two major lineages of pyrimethamine resistance in Melanesia. One has apparently originated in Melanesia, while the other originated in Southeast Asia and spread to Africa and Melanesia. This is the first unambiguous demonstration of the unique evolution of P. falciparum pyrimethamine resistance in Melanesia.
|
|
|---|
DNA preparation. Finger-prick blood (75 µl) was spotted onto chromatography filter paper ET31CHR (Whatman Limited, Kent, United Kingdom). Venous blood was transferred into heparin-containing test tubes. Parasite DNA was purified using a QIAamp DNA blood mini kit (QIAGEN, Germany) according to the manufacturer's instructions with some modifications (22).
dhfr genotyping. dhfr was amplified by PCR, and amplified products were directly sequenced with a BigDye Terminator 1.1 cycle sequencing kit in the ABI 377 DNA sequencer (GE Healthcare UK Ltd., Buckinghamshire, England) as previously reported (10, 19).
Microsatellite haplotyping. In order to determine the evolutionary history of pyrimethamine-resistant alleles of dhfr, we measured variation in the number of TA repeats at six microsatellite loci closely linked to the gene. These were located on chromosome 4, 0.1, 3.87, and 4.49 kb upstream and 0.52, 1.48, and 5.87 kb downstream of dhfr. In some cases, in order to estimate the limit of genetic hitchhiking, which is defined as a valley of reduced variation around dhfr, an additional six-microsatellite markers were analyzed at 7.55, 29.61, 57.68, and 363.33 kb upstream and 30.31 and 299.72 kb downstream of dhfr. Polymorphisms in these microsatellite markers were determined as previously described (12). Briefly, seminested PCR was performed using fluorescent end-labeled primers. Size variations in the amplified products were determined by electrophoresis on an ABI 377 sequencer and analyzed with GeneScan software (GE Healthcare UK Ltd.). Samples with two or more peaks at the same locus in the electropherogram were considered to be mixed infections and were excluded from further analysis.
Polymorphism between microsatellite markers is measured as variation in nucleotide length derived from various numbers of TA repeats. Microsatellite haplotypes harboring an association of bp 200-194-176-106-203-111 at microsatellite positions 4.49, 3.87, and 0.1 kb upstream and 0.52, 1.48, and 5.87 kb downstream of dhfr were designated "SEA" haplotypes, and those harboring an association of bp 220-202-156-100-205-111 were designated "Melanesia" haplotypes. Microsatellite haplotypes showing slight differences at one or two microsatellite markers from the SEA haplotype, e.g., at bp 204- 200-176-106-203-111 (underlining indicates the differences), were considered SEA variation haplotypes. Similarly, haplotypes displaying minor variation at one or two markers from the Melanesia haplotype, e.g., bp 220-202-156-100-205-113, were considered Melanesia variation haplotypes. Isolates showing mixed dhfr genotypes and/or microsatellite haplotypes were excluded from analysis.
Statistical analysis.
We calculated the expected heterozygosity (h) at each microsatellite locus as h = [n/(n 1)][1
pi2], where n is the number of infections sampled and pi2 is the frequency of the ith allele. The sampling variance of h was calculated according to the following formula (23), a slight modification of the standard diploid variance (13), [2/n(n 1)]{2(n 2)[
pi3 (
pi2)2] +
pi2 (
pi2)2}. A P value of <0.05 was considered statistically significant.
Nucleotide sequence accession number. The complete sequence of the allele identified has been submitted to the DDBJ and assigned accession number AB271908.
|
|
|---|
![]() View larger version (26K): [in a new window] |
FIG. 1. Frequency of dhfr genotypes in P. falciparum isolates from Papua New Guinea, the Solomon Islands, Vanuatu, Cambodia, and Thailand.
|
|
View this table: [in a new window] |
TABLE 1. Microsatellite polymorphisms in 15 P. falciparum isolates with wild-type dhfr or single-mutant dhfr
|
![]() View larger version (17K): [in a new window] |
FIG. 2. Repeat length variations of six microsatellite markers flanking dhfr in P. falciparum isolates with dhfr double (CNRNI) and triple (CIRNI)/quadruple (CIRNL) mutants. x axes, size (bp) of microsatellite markers. y axes, frequency of microsatellite alleles. a, location of microsatellite marker (distance from dhfr).
|
|
View this table: [in a new window] |
TABLE 2. Expected heterozygosity of microsatellite markers in P. falciparum isolates
|
![]() View larger version (26K): [in a new window] |
FIG. 3. Frequency of microsatellite haplotype flanking dhfr in P. falciparum isolates with dhfr double, triple, or quadruple mutants from Papua New Guinea, the Solomon Islands, Vanuatu, Cambodia, and Thailand. a, isolate (n = 1) carrying a hybrid of the SEA and Melanesia haplotypes (it was of SEA haplotype upstream and Melanesia haplotype downstream of dhfr, bp 200-194-176-100-205-111).
|
Genetic hitchhiking in dhfr double-mutant parasites from Papua New Guinea. These results suggest that the dhfr double mutants present today in Melanesia emerged independently in Southeast Asia and Melanesia. To determine the history of these two lineages, we measured the extent of genetic hitchhiking, which is determined by the distance of reduced microsatellite variation around dhfr. For this purpose, the variance of h at 12 microsatellite markers spanning 363.33 kb upstream to 299.72 kb downstream of dhfr was measured for Papua New Guinean isolates with a dhfr double mutant displaying a SEA/SEA variation haplotype (n = 17) and those displaying the Melanesia/Melanesia variation haplotype (n = 64) (Fig. 4). The patterns of genetic hitchhiking in both haplotypes were similar within a distance of 58 kb upstream and 30 kb downstream of dhfr. These results suggest that these two lineages, both carrying the same point mutations (CNRNI), appeared coincidentally in Papua New Guinea.
![]() View larger version (11K): [in a new window] |
FIG. 4. Reduced microsatellite polymorphism near dhfr in P. falciparum isolates with CNRNI genotype from Papua New Guinea.
|
|
|
|---|
A recent study has reported multiple origins of dhfr resistance within Kenya (6). However, care must be taken when basing conclusions about the origins of drug resistance on microsatellite variation from areas of high endemicity, such as Kenya. Two factors are likely to affect microsatellite polymorphism in areas of intense transmission. First, new microsatellite haplotypes are easily generated by meiotic recombination because of a high recombination rate and high prevalence of mixed-haplotype infections. Second, interallelic recombination within a microsatellite may generate new microsatellite alleles by an unequal crossing-over mechanism. Indeed, in a study by Roper et al. (21), nonuniform microsatellite haplotypes were noticeable around dhfr in pyrimethamine-resistant African isolates. These factors may be less important in areas of low transmission, such as Southeast Asia, and so do not compromise the conclusions of the present study.
Genetic hitchhiking reduces the expected heterozygosity of microsatellite markers around a selected gene, resulting in a valley of reduced variation. However, this association is easily broken down by recombination, resulting in a narrowing of the selection valley as the number of generations increases. In this study, we compared the selection valleys around the dhfr gene in two dhfr double mutants carrying SEA and Melanesia microsatellite haplotypes from Papua New Guinea. In both haplotypes, the microsatellite patterns within the valley were very similar from a distance of 58 kb upstream to 30 kb downstream of dhfr. The size of a selection valley is determined by several different parameters: the strength of the selection pressure on the mutant allele, the frequency of recombination, the transmission intensity, and the number of parasite generations since the emergence of the selected allele (16). In this analysis, these four parameters could be considered equal because all isolates were sampled from the same area. Thus, these results indicate that two ancestors of the dhfr double mutant in Papua New Guinea emerged coincidentally: one came from Southeast Asia, and the other arose independently within Melanesia. Although the appearance of the two resistant lineages emerged nearly simultaneously, we consider that the Melanesian-resistant type might have appeared slightly earlier than the influx of the SEA-resistant type. This is because if the SEA-resistant type migrated to Papua New Guinea earlier, it would have swept away microsatellite polymorphisms linked to the wild dhfr. Therefore, the possibility that a novel dhfr-resistant type having distinct microsatellite haplotypes appeared soon after the sweep in Papua New Guinea seems very unlikely.
The way drugs are used within regions of endemicity affects the generation and selection of resistant alleles. SP was widely used in Thailand and Cambodia during the 1970s and 1980s as the first-line treatment for uncomplicated malaria. In Melanesian countries, SP was introduced as a first-line treatment during the mid-1990s. Up until this time, pyrimethamine monotherapy was infrequently used for the treatment of malaria and other infections. It is difficult, therefore, to explain the nearly simultaneous emergence of pyrimethamine-resistant parasites in this area. It is, however, possible to attribute the emergence and spread of the triple and quadruple dhfr mutants to the widespread use of SP in Southeast Asia during the 1970s and 1980s.
We consider that weak and persistent pyrimethamine pressure by medicated salt projects, a form of mass drug administration, may explain the first selection of pyrimethamine-resistant parasites (15, 25). In the late 1950s and early 1960s, medicated salt projects were carried out in four different endemic regions: Indonesian Papua, the Thailand-Cambodia border, Ghana, and Iran (15, 25). In Southeast Asia and Melanesia, this project was carried out from 1959 to 1962 and from 1960 to 1961, respectively. Pyrimethamine resistance in P. falciparum was reported within 3 months of the start of the project (1960) in Indonesian Papua (7). At the Thailand-Cambodia border, resistance also developed quickly and pyrimethamine was replaced by chloroquine beginning in 1961. A large number of people treated with subcurative doses of antimalarial drugs present ideal conditions for the emergence of drug resistance (4, 27, 28). This, combined with the long half-life of pyrimethamine (116 h) (29), would have facilitated the emergence of drug resistance within the areas covered by the medicated salt project: Indonesian Papua and the Thailand-Cambodia border.
The dhfr double mutant (CNRNI), which confers moderate resistance to pyrimethamine, is widely distributed in Africa, West Asia, and Melanesia. However, whether this mutant is regularly selected de novo or whether it has spread from a limited number of foci of emergence is not known. In the present study, only two distinct microsatellite haplotypes (SEA and Melanesia haplotypes) were observed in a total of 184 P. falciparum isolates harboring the dhfr double mutant from Papua New Guinea, Vanuatu, and the Solomon Islands, suggesting that the generation of two mutations at positions 59 and 108 in dhfr is not frequent. In laboratory isolates, key point mutations in dhfr have occurred at frequencies as high as 2.5 x 109 per parasite replication, which predicts the generation of one mutant parasite in every malaria patient, assuming the number of parasites to be 1010 to 1012 in every infection (14). Consistently, the expected heterozygosities at microsatellite markers around dhfr were comparable between the wild-type and single dhfr mutant parasites. Thus, the initial mutation at position 108 in dhfr may occur relatively frequently (12, 20), but the generation and selection of an additional mutation at position 59 appear to be considerably less frequent. Mutations that render pathogens resistant to drug treatment are often associated with a loss of fitness (8, 11, 26). Resistant mutants may themselves develop compensatory mutations, which could then allow them to grow and survive in competition with wild-type forms (5, 26). The discrepancy between the frequent generation of the mutation at position 108 in dhfr and the rare occurrence of the dhfr double mutant as observed in this study may thus be reconciled by the requirement of complex compensatory mutations in a locus other than dhfr for restoring parasite fitness in natural populations.
In the present study, our samples of Papua New Guinea and Thailand were from individuals with clinical malaria, while samples from other sites were from cross-sectional studies of asymptomatic individuals. Symptomatic patients usually have higher parasite densities than do asymptomatic individuals. Thus, we cannot exclude the possibility that prevalences of microsatellite haplotypes may differ between isolates that cause disease and those that do not produce symptoms. However, we do not consider it very likely because there was no significant difference in the frequency distribution of genotypes of an antigen gene (msp1) between clinical patients and asymptomatic individuals in Melanesia (23).
In conclusion, this study provides strong evidence for the unique and independent origin of pyrimethamine resistance in Melanesia. The dhfr mutant, perhaps emerging from West Papua, has the same double mutations found in other geographic areas but distinct microsatellite haplotypes flanking the gene. Our results also show that the generation of double mutants with mutations at positions 59 and 108 of dhfr is a rare event, and this double mutation may be a first rate-determining step for the stable persistence of pyrimethamine resistance in P. falciparum.
This work was supported by a grant-in-aid of the Ministry of Health, Labor and Welfare (H17-Sinkou-ippan-019); grants-in-aid for scientific research from the Ministry of Education, Culture, Sports, Science, and Technology of Japan (17590380, 18390131, 18073013, and 18GS031); the Japan International Cooperation Agency (JICA) under the Integrated Cooperative Research for Malaria Control (ICRMC) project between the University of Papua New Guinea and Tokyo Women's Medical University; the Core University System Exchange Programme under the Japan Society for the Promotion of Science, coordinated by the University of Tokyo and Mahidol University; an institutional grant from the Swedish Foundation for International Cooperation in Research and Higher Education (STINT); and a project grant of the Swedish Research Council (VR, 2005-6836).
Published ahead of print on 8 January 2007. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»