ANTIMICROBIAL AGENTS AND CHEMOTHERAPY, Aug. 2002, p. 2482–2489 0066-4804/02/$04.00⫹0 DOI: 10.1128/AAC.46.8.2482–2489.2002 Copyright © 2002, American Society for Microbiology. All Rights Reserved.
Vol. 46, No. 8
Sulfadoxine-Pyrimethamine Resistance in the Rodent Malaria Parasite Plasmodium chabaudi Karen Hayton,1† Lisa C. Ranford-Cartwright,2 and David Walliker1* Institute of Cell, Animal and Population Biology, University of Edinburgh, Edinburgh EH9 3JT,1 and Division of Infection and Immunity, Institute of Biomedical and Life Sciences, University of Glasgow, Glasgow G12 8QQ,2 United Kingdom Received 1 April 2002/Returned for modification 25 April 2002/Accepted 9 May 2002
We have studied resistance to sulfadoxine-pyrimethamine (S/P) in the rodent malaria parasite Plasmodium chabaudi. A stable S/P-resistant mutant, AS(50S/P), was selected by drug treatment of a clone, AS(PYR), already resistant to pyrimethamine. The sequences of the P. chabaudi dhfr and dhps genes were obtained and found to be identical in AS(50S/P) and AS(PYR), showing that resistance to S/P in AS(50S/P) was not due to additional mutations in either gene. AS(50S/P) was crossed with a drug-sensitive clone, AJ, and 16 independent recombinant progeny were obtained. These clones were phenotyped for their susceptibility to S/P and to sulfadoxine and pyrimethamine separately. Pyrimethamine resistance was invariably associated with S/P resistance, but no correlation was found between resistance to S/P and resistance to sulfadoxine. Quantitative trait locus analysis of the progeny with 31 chromosome-specific markers showed that mutant P. chabaudi dhfr, or one or more genes closely linked to it, was a major determinant of S/P resistance. In addition, the inheritance of genes on chromosomes 5 and 13 from the sensitive parent appeared to contribute to the level of resistance observed. These results demonstrate that the S/P resistance of the AS(50S/P) mutant of P. chabaudi does not involve mutation in dhps and is not due simply to a combination of two genes determining resistance to pyrimethamine and sulfadoxine separately. (2, 22, 43). A major problem in studying S/P resistance in P. falciparum is the lack of a reliable in vitro test for phenotyping parasite responses to the S/P combination, due to variations in levels of the antagonists PABA and folic acid in serum and erythrocytes used for routine culture (21, 41, 42). Therefore, the identification of alleles associated with S/P resistance in natural populations has been based only on molecular assays to detect genotypes associated with therapeutic failure (15, 29, 39). However, predictions of the outcome of S/P treatment based on mutations in P. falciparum dhfr and dhps have had only limited success (15, 22, 23, 28). We report here a genetic study on S/P resistance using the rodent malaria species Plasmodium chabaudi, which is probably the best rodent model for studies on drug resistance in P. falciparum (7). Pyrimethamine resistance in P. chabaudi is known to be associated with the presence of Asn-106 in its DHFR, which is equivalent to Asn-108 in P. falciparum DHFR (9). We first describe the selection of a mutant exhibiting stable S/P resistance, starting from a P. chabaudi clone already resistant to pyrimethamine. We then report the isolation and sequence of the P. chabaudi pppk-dhps gene, which encodes the bifunctional protein hydroxymethylpterin pyrophosphokinase (PPPK)-DHPS, as well as those of the P. chabaudi dhfr-ts gene, encoding DHFR-thymidylate synthase, from the S/P-sensitive and -resistant clones. No additional mutations are seen in either gene following the development of S/P resistance. Genetic crossing work shows that S/P resistance is linked to pyrimethamine resistance but is independent of resistance to sulfadoxine alone. Quantitative trait locus (QTL) analysis of the cross progeny shows that P. chabaudi dhfr on chromosome 7 and genes on chromosomes 5 and 13 are strongly linked to S/P resistance and that another gene determining resistance may be located on chromosome 4.
The spread of chloroquine resistance in Plasmodium falciparum has favored the combination of sulfadoxine and pyrimethamine (S/P) (Fansidar) as one of the most important and widely used treatments for malaria (3). However, S/P resistance is now widespread in Southeast Asia and South America and is spreading through Africa (44). While the genetic mechanisms involved in resistance to pyrimethamine are well understood, only circumstantial evidence on the genetic basis of S/P resistance is available. Pyrimethamine is an antifolate drug which inhibits dihydrofolate reductase (DHFR), an essential enzyme in the parasite’s folic acid pathway. Sulfadoxine, a sulfa drug, is an analogue of p-aminobenzoic acid (PABA) and competitively inhibits dihydropteroate synthase (DHPS), which also is required for folate biosynthesis. In combination, these drugs exhibit a high level of synergy (10). In P. falciparum, pyrimethamine resistance is clearly due to amino acid changes in DHFR, in particular a change at position 108 of serine (Ser) or threonine (Thr) to asparagine (Asn) (27, 30, 45). Similarly, genetic and transfection studies have shown that sulfadoxine resistance in vitro is associated with mutations in P. falciparum dhps (4, 33–35, 40). No experimental work has been carried out on the genetics of S/P resistance. In general, it is assumed that parasites exhibiting S/P resistance contain mutations in both P. falciparum dhfr and dhps, although some studies suggest that S/P failure may be determined simply by several mutations in dhfr alone
* Corresponding author. Mailing address: Institute of Cell, Animal and Population Biology, University of Edinburgh, Edinburgh EH9 3JT, United Kingdom. Phone: (44) 131 650 5548. Fax: (44) 131 650 6564. E-mail:
[email protected]. † Present address: LMVR/NIAID, National Institutes of Health, Bethesda, MD 20892. 2482
VOL. 46, 2002
SULFADOXINE-PYRIMETHAMINE RESISTANCE IN MALARIA MATERIALS AND METHODS
Parasites and hosts. The parasites used in this work were cloned lines of P. chabaudi isolates AS and AJ originating from independent wild-caught Thamnomys rutilans thicket rats from the Central African Republic (8). Clone AS(PYR) is pyrimethamine resistant and was selected from a sensitive AS clone designated AS(SENS) (38). Clone AJ is a wild-type parasite which is sensitive to pyrimethamine and innately resistant to sulfadoxine, as seen in all other isolates of this species so far examined. Parasites were maintained essentially as described previously (37), using adult CBA/Ca mice (University of Edinburgh, Edinburgh, United Kingdom, or Bantim and Kingman, Hull, United Kingdom) fed on SDS Formula Number 1 diet (Special Diets Services Ltd., Essex, United Kingdom), with the drinking water supplemented with 0.05% PABA to enhance parasite growth (14). Animal procedures were conducted under license following the United Kingdom Animals Scientific Procedures. Drugs. Pyrimethamine base was obtained from Wellcome Research Laboratories (Kent, United Kingdom). Sulfadoxine and Fansidar were donated by Roche Pharmaceuticals. Drugs were made up as solutions in dimethyl sulfoxide (Sigma, Poole, United Kingdom) and administered orally to mice. Selection for S/P resistance. One hundred outbred mice (Bantim and Kingman) were each inoculated intraperitoneally (i.p.) with 106 AS(PYR)-infected red blood cells on day 0. A high curative dose of S/P (125/5.25 mg kg⫺1) was administered to each mouse daily from day 3 for a total of 4 days. Blood smears were then examined every 2 days from day 10. Blood was taken from any mice exhibiting parasitemia and stored as stabilates in liquid nitrogen. Lines chosen for further studies were cloned by limiting dilution (38). One of these clones, designated AS(50S/P), was used in the present work. Tests for drug resistance. A total of 106 parasitized red blood cells were inoculated i.p. into groups of five 6- to 9-week-old CBA/Ca inbred mice (Bantim and Kingman). Dilutions of pyrimethamine, sulfadoxine, or S/P were administered to three of the mice in each group at 3 h after inoculation and at the same time of day for the next 3 days. The remaining two mice in each group were undrugged controls. Drug doses were expressed as milligram per kilogram of mouse body weight. Blood smears were taken from day 6 onwards, and parasitemias were counted microscopically. Each parasite clone was tested for responses to each drug at least twice. Stability of resistance in the absence of drug. Parasites were examined for the stability of their drug response by weekly passage through mice without drug pressure. Each week, the parasites were also subpassaged into further groups of mice and tested for their responses to pyrimethamine, sulfadoxine, and S/P. Their resistance was also tested following transmission through mosquitoes (Anopheles stephensi) (38). Sequencing of the P. chabaudi pppk-dhps and dhfr-ts genes. To sequence the P. chabaudi pppk-dhps gene, known sequences of dhps (4, 12, 17, 24, 31, 33, 36) were aligned using the University of Wisconsin Genetics Computer Group DNA analysis software to identify two regions of sequence conservation. Degenerate oligonucleotides designated DHPS1 [TGATTG(G/T)TT(A/T)(C/T)AGATGG] and DHPS2 [G(A/T)GG(A/C)TTC(A/T)CC(A/T)CC(A/T)ATATC] were designed for use in PCR. Genomic DNA was isolated from AS(SENS) as described previously (32). A PCR product of approximately 100 bp was obtained and purified using Concert (Gibco BRL) according to the manufacturer’s instructions and sequenced using BigDye (Perkin-Elmer) chemistry on an ABI 377 Prism Sequencer. Sequences were extended to include the complete P. chabaudi dhps and pppk genes by PCR walking using P. chabaudi Vectorette II (Genosys) libraries generated from genomic AS(SENS) DNA as described previously (1). Gene-specific primers were designed from the initial PCR product and used to generate additional sequence from the libraries, according to the manufacturer’s instructions. The resulting sequence was used to design oligonucleotides (DHPS3, GAAAAAACAAATATAGTTGGC; DHPS4, AATTCTATACCTAGGTATG CC; DHPS5, TTTCTTCGACAATTGCATG; DHPS6, TATGAGGATGTTGT ATATG; PPPK1, GAAAAGCATAAAAATAAATACTACAC; PPPK2, GAAA CTTTTAAAAACTACGGC; and PPPK3, CTCAGATAGTCTATTAAATGTA AG) to obtain the sequence of the P. chabaudi pppk-dhps gene from AS(PYR), AS(50S/P), AJ, and each of the lines that survived S/P selection in order to detect any mutation events. To sequence P. chabaudi dhfr-ts, fragments of the gene were amplified from parasite DNA by using primers designed from available sequences (9, 11). The primers used to amplify the gene from genomic DNAs of AS(PYR) and AS(50S/P) were DHFR1 (GAAGATATCTCTGAAATATTCGATATATAT G), DHFR2 (TTCATTGATATCTGGAAATAAAACATCAC), DHFR4 (CAA GTATAAGCATTATTTATTCTTG), DHFR5 (GATTTTGTAATGGCTTAA ATTTTGAG), TS1 (AGAGAAGTTAATGACTTAGGTCC), and TS2 (TAAG
2483
CTGCCATATCCATACTG). Amplification products were purified and sequenced as described for P. chabaudi pppk-dhps above. Crosses and recovery of progeny. In order to determine the genetic basis of resistance to S/P, AS(50S/P) was crossed with the unrelated drug-sensitive clone AJ (5). Two crosses were made using the basic procedure of Walliker et al. (38), by allowing mosquitoes (A. stephensi) to take up blood containing gametocytes of a mixed infection of the two clones, thus allowing cross-fertilization of gametes and recombination to occur. The resulting sporozoites were used to infect mice. The blood forms which developed in these animals, designated the progeny of the cross, were cloned by injecting inocula containing an average of 0.75 parasitized cell i.p. into mice (38). Inheritance of chromosome-specific markers and QTL analysis. Chromosome markers developed by Carlton et al. (5) were used in linkage analysis of the cross progeny. Details of these markers can be found at http://www.ncbi.nlm.nih.gov /Malaria/Rodent/ctabgenmarker.html. Details of marker 100 (pccg10) were provided by P. Hunt (personal communication). The inheritance pattern of each marker was compared with the inheritance of S/P, sulfadoxine, and pyrimethamine responses of each progeny clone. The data were analyzed by QTL methodology (16, 18) to determine whether inheritance patterns of the three types of resistance were consistent with a single or multigenic basis for each character. Since malaria parasites are haploid, this was simply carried out by calculating the mean difference in log parasitemias on day 6 following treatment with either 25/1.25 mg of S/P kg⫺1, 10 mg of pyrimethamine kg⫺1, or 150 mg of sulfadoxine kg⫺1 between groups of progeny clones inheriting either the AS(50S/P) or the AJ allele of each marker. The statistical significance was determined using Student’s t test. The Holm procedure was used to correct genome-wide thresholds to account for the analysis of multiple chromosomal markers and drug resistance traits (13). Nucleotide sequence accession number. A 2,489-bp sequence encoding a portion of the P. chabaudi pppk-dhps gene has been deposited in the DDBJ/EMBL/ GenBank databases under accession number AJ302077.
RESULTS Selection for S/P resistance. An S/P-resistant mutant of P. chabaudi was selected by single-step exposure of a clone already resistant to pyrimethamine, AS(PYR), to a high dose of the drug combination (125/5.25 mg of S/P kg⫺1) given over 4 days. Fifteen of 100 mice infected and treated with S/P developed patent parasitemias on day 18 after inoculation, 15 days after the last dose of S/P. The growth patterns of each of the S/P-selected lines were compared under S/P pressure to that of the parent clone AS(PYR). Seven had drug responses comparable to that of AS(PYR), while eight grew more quickly following treatment with 50/2.5 mg of S/P kg⫺1 (data not shown). These parasite lines were classified as S/P resistant. One S/P-resistant line was cloned by limiting dilution and designated AS(50S/P). This clone was tested for its response to four daily doses of either (i) 25/1.25 mg of S/P kg⫺1, (ii) 10 mg of pyrimethamine alone kg⫺1, or (iii) 150 mg of sulfadoxine alone kg⫺1 (Fig. 1). Under S/P treatment, AS(50S/P) grew almost as well as parasites in undrugged mice; AS(PYR) grew much more slowly, with recrudescence being observed only on day 10; and AS(SENS) failed to grow during the 10-day period (Fig. 1A). Under pyrimethamine treatment, AS(50S/P) grew slightly faster than AS(PYR), while AS(SENS) was completely sensitive (Fig. 1B). Under sulfadoxine treatment, AS(50S/P) was slightly resistant, recrudescing on day 8 following drug treatment (Fig. 1C). In contrast, AS(PYR) was completely eliminated by sulfadoxine, a result which was not unexpected given that selection for pyrimethamine resistance in rodent malaria parasites invariably results in an increased PABA requirement of the parasite and a consequent increased sensitivity to sulfa drugs (14). AS(SENS) was unaffected by sulfadoxine, and its growth was comparable to that observed in
2484
HAYTON ET AL.
ANTIMICROB. AGENTS CHEMOTHER.
FIG. 1. Assessment of the drug sensitivities of AS(SENS) 䊐, AS(PYR) E, and AS(50S/P) F. CBA/Ca mice were infected on day 0 and treated for 4 days with either 25/1.25 mg of S/P kg⫺1 (A), 10 mg of pyrimethamine kg⫺1 (B), or 150 mg of sulfadoxine kg⫺1 (C). The data shown are the means ⫾ standard errors from representative experiments.
undrugged mice. The responses of AJ to all three drug regimens were comparable to those of AS(SENS). Stability of resistance in the absence of drug. AS(50S/P) was passaged weekly for 20 weeks through uninfected mice without drug pressure, and each week parasites were tested for their drug responses. AS(50S/P) retained its original level of resistance to S/P, pyrimethamine alone, and sulfadoxine alone. It was also transmitted through A. stephensi into mice and was found to have retained identical responses. Sequences of the P. chabaudi pppk-dhps and dhfr-ts genes. Using degenerate oligonucleotides and Vectorette PCR (1), a 2,489-bp sequence which appeared to encode a portion of the P. chabaudi pppk-dhps gene was identified (data not shown). The sequence lacked approximately 10 amino acids at the amino terminus of the P. chabaudi PPPK, based on comparison with the P. falciparum gene. The putative P. chabaudi dhps gene was analyzed in all lines that survived drug selection, including AS(50S/P) (data not shown). Sequences from all of these lines were identical to that from AS(PYR). We also sequenced the pppk portions of the P. chabaudi pppk-dhps genes from AS(PYR), AS(50S/P), and the drug-sensitive clone AJ. Again, all three sequences were identical. These results suggest strongly that resistance to the S/P combination is not due to mutations in P. chabaudi pppk-dhps. Although we were unable to obtain sequences for approximately 10 codons at the amino terminus of P. chabaudi PPPKDHPS, no polymorphisms have been reported in this region among alleles of P. falciparum pppk-dhps. P. chabaudi dhfr-ts was sequenced from AS(SENS), AS (PYR), AS(50S/P), and AJ and compared to sequences previously published for P. chabaudi (data not shown) (9, 11). In addition, P. chabaudi dhfr was sequenced from all the S/Pselected lines. The AS(SENS) dhfr and the AJ dhfr contained a Ser (AGT) at position 106. In contrast, as expected, the AS(PYR) dhfr and AS(50S/P) dhfr were identical and possessed Asn (AAT) at this position (9). Cross results. Two crosses were carried out between AS(50S/P) and AJ. Thirty-three parasite lines were obtained in mice by limiting dilution of the progeny of the crosses, 14 from cross 1 and 19 from cross 2. Examination of these lines for 31 markers distinguishing each parent revealed that six contained two alleles of several parent clone markers and were thus mixtures of more than one clone. The remaining 27 lines were classified as pure clones by their possession of single alleles of all markers. Two of these were pure AJ parent type, one was pure AS parent type, and the remaining 24 clones were recom-
binants, exhibiting nonparental combinations of parent alleles. Eight of these clones had genotypes identical to that of at least one other recombinant clone and were therefore considered to be derived from duplicate cloning events. In total, therefore, 16 independent recombinants were identified, 8 from cross 1 and 8 from cross 2, and these were used in the subsequent linkage analysis work (Fig. 2). Drug-responses of progeny clones. The growth rates of each progeny clone treated with S/P, with pyrimethamine alone, and with sulfadoxine alone were assessed over a 14-day period and were compared with those of each of the two parents and AS(PYR) (Table 1; Fig. 2). Each progeny clone was tested in this way twice. For almost all clones, the drug response phenotypes proved to be similar in each test. Exceptions occurred when, for example, parasites in one of the three drug-treated mice in a test underwent a course of infection different from that in the other two animals. In such instances, further tests were made until at least two successive tests gave comparable results in all mice. Drug responses were assessed over a 14-day period after drug treatment. The parasitemias on day 6 proved to be a clear measure of the different phenotypes, and for simplicity only these data are shown (Table 1). With regard to S/P, three phenotypes were detected, a sensitive phenotype (similar to that of parent clone AJ), a highly resistant phenotype [similar to that of AS(50S/P)], and a nonparental low-resistance phenotype [like that seen in AS (PYR)]. For pyrimethamine, two phenotypes were observed, sensitive (similar to AJ) or resistant [similar to AS(50S/P) and AS(PYR)]. For sulfadoxine, three phenotypes were identified, a high-resistance phenotype (like that of AJ), a low-resistance phenotype [like that of AS(50S/P)], and a nonparental sensitive phenotype [like that of AS(PYR)]. Four progeny (1107/14, 1109/11, 1384/10, and 1538/12) exhibited the sulfadoxine-sensitive phenotype, and all of these possessed a high resistance to S/P as well as resistance to pyrimethamine (Table 1). Inheritance of chromosome-specific markers and QTL analysis. The inheritance patterns of the 31 markers and the drug responses of the 16 progeny clones are shown in Fig. 2. The marker distribution and genetic map data will be presented elsewhere. The S/P and pyrimethamine responses of the progeny were compared. Figure 3B shows that progeny inheriting the AJtype dhfr and progeny inheriting the AS(50S/P)-type dhfr segregate within distinct groups of pyrimethamine-sensitive and -resistant phenotypes. While all progeny that inherited the AJ-type dhfr were S/P sensitive, progeny that inherited the AS
VOL. 46, 2002
SULFADOXINE-PYRIMETHAMINE RESISTANCE IN MALARIA
2485
FIG. 2. Inheritance of 31 chromosomal markers among AJ, AS(50S/P), and 16 progeny clones from two AJ ⫻ AS(50S/P) crosses. Progeny clones are listed vertically. Roman numerals refer to chromosome number, marker numbers are shown along the bottom, and markers are ordered such that the number of crossovers per chromosome was minimal (5). An asterisk marks the inheritance profile for P. chabaudi dhfr-ts (marker 6). Black boxes indicate inheritance of AS(50S/P)-type alleles, open boxes indicate inheritance of AJ-type alleles, and grey boxes indicate data not available. Drug phenotypes are shown on the right, where black boxes indicate the inheritance of a resistant phenotype, white boxes indicate the inheritance of a sensitive phenotype, and hatched boxes identify progeny with an intermediate phenotype.
(50S/P)-type dhfr are distributed across a range of S/P values, suggesting that at least two genes are required for S/P resistance. Therefore, the data were analyzed by QTL methodology to identify loci linked to the S/P resistance phenotype. The results of the QTL analysis are shown in Table 2. For each marker, the difference between the mean log parasitemia of grouped progeny clones inheriting the AS(50S/P)-type allele and that of clones with the AJ-type allele was first tested for statistical significance by using Student’s t test. For example, for marker 4 (msp1) the mean log parasitemia of progeny inheriting the AS(50S/P)-type allele was 0.73 on day 6 following treatment with 25/1.25 mg of S/P kg⫺1, while the mean log
parasitemia of all progeny inheriting the AJ-type allele was 0.82 (Table 2 and Fig. 3A). The difference between the two means (0.09) was insignificant (P ⫽ 0.46) at n ⫽ 14. It could be concluded, therefore, that marker 4 is unlikely to be linked to a locus conferring resistance to 25/1.25 mg of S/P kg⫺1. Following application of the Holm procedure to correct for false positives due to multiple testing, two loci showed evidence of strong linkage to drug responses. The AS(50S/P) allele of P. chabaudi dhfr (marker 6) on chromosome 7 was strongly linked not only to resistance to pyrimethamine (as expected) (P ⬍ 0.001) but also to resistance to S/P (P ⬍ 0.01) (Table 2 and Fig. 3B). In addition, mutant P. chabaudi dhfr was
TABLE 1. Drug resistance phenotypes of AJ, AS(PYR), AS(50S/P), and the AJ ⫻ AS(50S/P) progeny clones Clone
AJ AS(PYR) AS(50S/P) 1107/10 1107/14 1108/8 1109/11 1110/4 1110/8 1110/10 1110/13 1351/2 1384/10 1385/2 1396/1 1398/14 1538/12 1538/14 1539/1 a b
% Parasitemia (mean ⫾ SE)a
Phenotypeb
Untreated control
S/P (25/1.25 mg kg⫺1)
Pyrimethamine (10 mg kg⫺1)
Sulfadoxine (150 mg kg⫺1)
S/P
Pyrimethamine
Sulfadoxine
40.2 ⫾ 3.3 29.6 ⫾ 3.3 31.4 ⫾ 3.3 44.0 ⫾ 9.8 27.8 ⫾ 5.4 38.3 ⫾ 8.5 46.2 ⫾ 3.7 58.6 ⫾ 1.3 52.6 ⫾ 0.5 36.7 ⫾ 15.0 35.4 ⫾ 11.6 18.5 ⫾ 7.8 17.7 ⫾ 9.3 37.1 ⫾ 11.1 43.3 ⫾ 9.0 39.5 ⫾ 7.6 28.2 ⫾ 2.4 26.2 ⫾ 7.0 40.9 ⫾ 7.4
0.0 ⫾ 0.0 0.3 ⫾ 0.2 14.3 ⫾ 3.4 1.9 ⫾ 1.9 9.6 ⫾ 3.2 0.0 ⫾ 0.0 23.0 ⫾ 22.3 40.9 ⫾ 2.1 0.0 ⫾ 0.0 0.0 ⫾ 0.0 0.0 ⫾ 0.0 0.0 ⫾ 0.0 10.0 ⫾ 6.6 0.2 ⫾ 0.2 19.9 ⫾ 1.4 14.6 ⫾ 10.0 17.7 ⫾ 1.1 24.5 ⫾ 1.5 39.2 ⫾ 0.9
0.0 ⫾ 0.0 3.8 ⫾ 2.7 15.0 ⫾ 3.1 10.3 ⫾ 1.65 18.4 ⫾ 15.9 0.0 ⫾ 0.0 24.0 ⫾ 21.4 13.3 ⫾ 4.6 0.0 ⫾ 0.0 0.0 ⫾ 0.0 0.0 ⫾ 0.0 0.0 ⫾ 0.0 5.3 ⫾ 0.5 29.8 ⫾ 23.6 56.6 ⫾ 3.2 32.3 ⫾ 1.7 11.1 ⫾ 6.7 13.3 ⫾ 2.1 32.8 ⫾ 1.0
21.0 ⫾ 2.0 0.0 ⫾ 0.0 0.7 ⫾ 0.5 0.0 ⫾ 0.0 0.4 ⫾ 0.4 33.7 ⫾ 1.6 0.0 ⫾ 0.0 0.4 ⫾ 0.2 20.0 ⫾ 4.5 11.8 ⫾ 2.3 33.3 ⫾ 18.8 26.6 ⫾ 4.6 0.0 ⫾ 0.0 0.1 ⫾ 0.1 0.3 ⫾ 0.3 5.8 ⫾ 0.8 0.0 ⫾ 0.0 14.0 ⫾ 3.4 8.0 ⫾ 2.6
S LR R R R S R R S S S S R LR R R R R R
S R R R R S R R S S S S R R R R R R R
R S LR S LR R S LR R R R R S LR LR R S R R
On day 6 of the indicated treatment. Drug phenotypes were assessed over 14 days after drug treatment. R, resistant; S, sensitive; LR, low-level resistance.
2486
HAYTON ET AL.
ANTIMICROB. AGENTS CHEMOTHER.
FIG. 3. Scatter plots of S/P versus pyrimethamine (PYR) responses for the two marker classes at two loci, msp1 (marker 4) (A) and dhfr (marker 6) (B). Progeny that inherited the AS 50S/P-type allele are shown as open circles, while progeny that inherited the AJ-type allele are shown as filled circles. Note that the single filled circle in panel B represents six completely sensitive progeny. treatm, treatment.
linked strongly to sensitivity to sulfadoxine (P ⬍ 0.001). The differences between the mean log parasitemias of the two groups of progeny clones inheriting either the AJ or AS(50S/P) dhfr allele were highly significant in each analysis. In contrast, the AS(50S/P) allele of marker 19 (g6pd, chromosome 13) was linked to sensitivity to S/P (although not after Holm correction), to sensitivity to pyrimethamine (P ⬍ 0.01), and to resistance to sulfadoxine (P ⬍ 0.01). Marker 23 (pfcrk3) on this chromosome showed similar strong linkage to pyrimethamine sensitivity and sulfadoxine resistance before correction, although there was no association with S/P response. Other loci which appeared to show association with drug responses before Holm correction were as follows. The AS (50S/P) allele of marker 88 (anonymous marker P.1) on chromosome 4 showed linkage to S/P resistance. Conversely, the AS(50S/P) allele of marker 57 (Ag3035) on chromosome 5 was associated with sulfadoxine resistance but also with sensitivity to pyrimethamine and S/P. AJ-type alleles on chromosomes 5 and 13 showed a significant linkage with both S/P and pyrimethamine resistance and sulfadoxine sensitivity. DISCUSSION We have demonstrated here that stable resistance to S/P can be produced in the malaria parasite P. chabaudi by single-step treatment with this drug combination of a clone already resistant to pyrimethamine. The resistance was caused by a gene or genes distinct from dhfr or dhps, since neither of these genes contained novel mutations in the S/P-resistant mutant AS (50S/P) compared to the clone from which it was selected. Our results also demonstrate clearly that S/P resistance in this parasite does not necessarily depend on the parasite being resistant to pyrimethamine and sulfadoxine separately. This is shown by the appearance of S/P-resistant clones that were sensitive to sulfadoxine among the progeny of the cross between AS(50S/P) and AJ. The genetic basis of resistance to S/P has not been investigated experimentally in P. falciparum, mainly because of the difficulties of testing the response to this drug combination in
vitro. Sulfadoxine and pyrimethamine have a synergistic action against malaria parasites (10). However, it should be noted that sulfadoxine by itself has limited efficacy; when it was used alone for clinical treatment of P. falciparum, failure to clear parasites commonly occurred (26). Indeed, this is seen in our present work by the lack of sensitivity to sulfadoxine of the wild-type P. chabaudi clones AS(SENS) and AJ. Also, several previous studies using rodent malaria models have demonstrated that mutants resistant to pyrimethamine exhibit an increased sensitivity to sulfa drugs, due to an increased dependence on PABA (14, 26). This was evident in our present work in the increased sensitivity of pyrimethamine-resistant AS (PYR) to sulfadoxine in comparison to sensitive wild-type clones. AS(50S/P) was derived from AS(SENS) in two selection steps, first with pyrimethamine to produce AS(PYR) (38), which in the present work was treated with S/P to produce AS (50S/P). The simplest explanation for the development of AS (50S/P)’s high S/P resistance is that successive mutations occurred at each step. The initial selection of AS(PYR) by pyrimethamine treatment resulted in a mutation in P. chabaudi dhfr in which Ser was replaced by Asn at position 106 (9), causing resistance to pyrimethamine. In the present study, we have shown that AS(PYR) also possessed a slight resistance to S/P and a marked decrease in resistance to sulfadoxine compared to AS(SENS). The selection of AS(50S/P) from AS (PYR) produced high S/P resistance and an increase in sulfadoxine resistance, although not to the level seen in wild-type parasites. The precise mutation event or events causing the S/P resistance of AS(50S/P) are unknown. We have shown that there are no extra mutations in either dhfr or dhps in this mutant. However, our genetic analysis has shown that its resistance is probably determined by several genes. This can be deduced from the finding that progeny clones exhibiting various parental and nonparental combinations of resistance to S/P, sulfadoxine, and pyrimethamine were obtained in the cross. Certain patterns of resistance were apparent. Sensitivity to S/P was invariably associated with sensitivity to pyrimethamine and re-
VOL. 46, 2002
SULFADOXINE-PYRIMETHAMINE RESISTANCE IN MALARIA
TABLE 2. Analysis for linkage between inheritance of markers and drug resistance Chromosome
Marker
1
49
2
Mean difference in log parasitemiaa with: S/P
Pyrimethamine
Sulfadoxine
⫺0.08
⫺0.32
0.30
1
0.43
0.35
0.17
3
3 47
0.40 0.48
0.64 0.82
⫺0.83** (NS) ⫺0.38
4
88 89
0.86** (NS) ⫺0.14
0.57 0.38
⫺0.47 ⫺0.45
5
57 25 95 50
⫺0.91* (NS) ⫺0.83** (NS) 0.80 0.52
6
91 24
0.01 0.18
0.24 0.27
7
6 94
1.13*** (**) 0.50
1.28*** (***) 0.50
8
4
⫺0.09
⫺0.15
0.29
9
99 52
⫺0.21
⫺0.39
0.85* (NS)
10
39
0.91
0.73
0.07
11
43 15 29 93 8
0.18 ⫺0.09 0.09 0.70 0.17
1.02 ⫺0.61 ⫺0.13 0.50 ⫺0.01
⫺0.44 0.27 ⫺0.01 ⫺0.13 ⫺0.34
12
27
0.25
0.18
⫺0.50
13
7 37 45 19
⫺0.10 ⫺0.39 ⫺0.67* (NS) ⫺1.00** (NS)
⫺0.38 ⫺0.72 ⫺0.55 ⫺1.17*** (**)
⫺0.13 0.98* (NS) 0.27 1.15*** (**)
14
23 59 98
⫺0.57 0.56 0.18
⫺0.93*** (NS) 0.29 0.15
0.75** (NS) ⫺0.64 ⫺0.57* (NS)
⫺1.07** (NS) ⫺0.62 0.78 0.79* (NS)
1.07** (NS) 0.30 0.28 ⫺0.52 0.15 0.00 ⫺1.08*** (***) ⫺0.05
a Mean difference in log parasitaemia between progeny clones inheriting the AS(50S/P) versus the AJ allele of each marker. The difference between the means provides an estimate of the phenotypic effect of inheriting either allele. Positive values indicate that the AS(50S/P) marker allele is segregating with resistance, while negative values indicate that the AJ allele is segregating with resistance. The statistical significance of the difference between the two means was determined by (i) Student’s t test alone and (ii) Student’s t test followed by Holm-Bonferroni correction (shown in parentheses) (13). The analysis could not be carried out for marker 52 because no progeny inherited the AS(50S/P) allele of this locus. *, P ⬍ 0.05; **, P ⬍ 0.01; ***, P ⬍ 0.001; NS, not significant.
sistance to sulfadoxine. Sensitivity to sulfadoxine was invariably associated with resistance to pyrimethamine and to S/P. However, resistance to all three drugs was found in many clones, including various combinations of high and low levels of resistance to S/P and sulfadoxine which differed from that of the parent AS(50S/P). A notable finding was that S/P-resistant clones could be either resistant or sensitive to sulfadoxine, showing that these two types of resistance had an independent genetic basis. QTL analysis using 31 chromosome-specific markers has
2487
confirmed that mutant P. chabaudi dhfr (chromosome 7) causes pyrimethamine resistance, as expected from the earlier studies (9, 11, 27, 45). The analysis has also shown that this gene, or one or more genes closely linked to it, is a major determinant of S/P resistance and sulfadoxine sensitivity. This finding is consistent with the results reported for P. falciparum by Wang et al. (40), in which an unidentified gene closely linked to P. falciparum dhfr on chromosome 4 was linked to folate utilization and sulfadoxine response in vitro. In addition, our analysis suggests that a second locus for S/P resistance could be present on chromosome 4 and that loci for sulfadoxine resistance could be present on chromosomes 5 and 13. Further, despite the fact that the AJ parent was sensitive to both drugs, AJ alleles of loci on chromosomes 5 and 13 were associated with both S/P and pyrimethamine resistance among the cross progeny. The corresponding AS(50S/P) alleles were linked to sulfadoxine resistance, although AS(50S/P) has a low level of resistance to this drug. AJ has a slightly higher growth rate than AS(50S/P), and it is possible that genes on chromosomes 5 and 13 conferring growth rate advantages inherited from the AJ parent may have contributed to the level of resistance observed. It is also noteworthy that P. chabaudi dhpspppk is located on chromosome 13 (6). While we have ruled out mutation in P. chabaudi dhps-pppk as a cause of S/P and sulfadoxine resistance, it is quite possible that a linked gene is involved. Also, while we have detected no evidence of P. chabaudi dhps-pppk or dhfr-ts gene duplication in this work by pulsed-field electrophoresis of chromosomes of sensitive and mutant parasites (data not shown), it would be of interest to determine whether there are differences in the levels of expression of P. chabaudi dhps in resistant and sensitive parasites which could be due to upregulation of expression of this gene. A major problem with genome mapping in P. chabaudi is the present lack of chromosome-specific markers. Our ability to detect loci other than those closely linked to a genetic marker was not optimal (5). Due to the limited number of markers tested, false-positive associations could be expected. Our analysis of the inheritance of 31 markers for linkage to three drug-resistance phenotypes resulted in 93 marker contrast analyses. Procedures to correct the genome-wide thresholds to account for the analysis of multiple markers and drug resistance traits resulted in loss of power. Nonetheless, the application of the Holm procedure to our t test analysis confirmed the strong likelihood that genes on chromosomes 7 and 13 are determinants of resistance to the drugs studied here. S/P-resistant parasites have been selected in other rodent malarias, i.e., those caused by Plasmodium berghei and Plasmodium yoelii, by using a continuous low-dose method, but in these instances the resistant phenotypes were not stable once drug pressure was removed (19, 20, 25). The fact that in the present study stable resistance in AS(50S/P) was obtained after a single course of drug selection suggests that a mutation in a single gene, or at most a few genes, was the cause of the novel phenotype. Therefore, the selection of S/P resistance in P. chabaudi observed in the present study could reflect the evolution of S/P resistance in P. falciparum, which arose quickly in areas where pyrimethamine resistance and, most probably, mutations in P. falciparum dhfr were already widespread (26). We have shown here that the S/P resistance of AS(50S/P) is not due to mutation events in P. chabaudi dhfr additional to
2488
HAYTON ET AL.
Asn-106 or in P. chabaudi dhps. This does not, of course, rule out the possibility that mutations in these genes may be involved in other S/P-resistant mutants of P. chabaudi. While multiple mutations in both P. falciparum dhfr and P. falciparum dhps appear to confer S/P resistance in many P. falciparum isolates (15, 22, 23, 28, 29, 39), there are certainly exceptions which cannot be attributed to, for example, acquired immunity eliminating resistant forms. As a general rule, the resistance of a microorganism to a drug can be expected to have a variety of genetic causes, especially when resistance has arisen multiple times. Alternative mechanisms of S/P resistance, independent of mutations in P. falciparum dhps, may have been overlooked in P. falciparum due to the difficulties associated with studying the S/P response of this species in vitro. The ease with which S/P-resistant mutants can be selected in P. chabaudi, as well as the similarities in the basic biology of this species to that of P. falciparum, make it a very good model for studying the genetics of this as well as other types of malaria drug resistance (7).
ANTIMICROB. AGENTS CHEMOTHER.
16. 17.
18. 19. 20. 21. 22.
ACKNOWLEDGMENTS
23.
We thank Mike Barrett and Chris Janssen for the gift of the initial HindIII P. chabaudi Vectorette II library, Richard Fawcett and Ronnie Mooney for technical assistance, and Margaret Mackinnon for helpful advice. This work was supported by the Medical Research Council of the United Kingdom.
24.
25.
REFERENCES 1. Arnold, C., and I. J. Hodgson. 1991. Vectorette PCR: a novel approach to genome walking. PCR Methods Appl. 1:39–42. 2. Basco, L. K., R. Tahar, and P. Ringwald. 1998. Molecular basis of in vivo resistance to sulfadoxine-pyrimethamine in African adult patients infected with Plasmodium falciparum malaria parasites. Antimicrob. Agents Chemother. 42:1811–1814. 3. Bloland, P. B., E. M. Lackritz, P. N. Kazembe, J. B. O. Were, R. Steketee, and C. C. Campbell. 1993. Beyond chloroquine: implications of drug resistance for evaluating malaria therapy efficacy and treatment policy in Africa. J. Infect. Dis. 167:932–937. 4. Brooks, D. R., P. Wang, M. Read, W. M. Watkins, P. F. G. Sims, and J. E. Hyde. 1994. Sequence variation of the hydroxymethyldihydropterin pyrophosphate kinase:dihydropteroate synthase gene in lines of the human malaria parasite, Plasmodium falciparum, with differing resistance to sulfadoxine. Eur. J. Biochem. 224:397–405. 5. Carlton, J., M. Mackinnon, and D. Walliker. 1998. A chloroquine resistance locus in the rodent malaria parasite Plasmodium chabaudi. Mol. Biochem. Parasitol. 93:57–72. 6. Carlton, J. M. R., R. Vinkenoog, A. P. Waters, and D. Walliker. 1998. Gene synteny in species of Plasmodium. Mol. Biochem. Parasitol. 93:285–294. 7. Carlton, J. M. R., K. Hayton, P. V. L. Cravo, and D. Walliker. 2001. Of mice and malaria mutants: unravelling the genetics of drug resistance using rodent malaria models. Trends Parasitol. 17:236–242. 8. Carter, R., and D. Walliker. 1975. New observations on the malaria parasites of rodents of the Central African Republic Plasmodium vinckei petteri sbsp. nov. and Plasmodium chabaudi Landau, 1965. Ann. Trop. Med. Parasitol. 69:187–196. 9. Cheng, Q., and A. Saul. 1994. The dihydrofolate reductase domain of rodent malarias: point mutations and pyrimethamine resistance. Mol. Biochem. Parasitol. 65:361–363. 10. Chulay, J. D., W. M. Watkins, and D. G. Sixsmith. 1984. Synergistic antimalarial activity of pyrimethamine and sulfadoxine against Plasmodium falciparum in vitro. Am. J. Trop. Med. Hyg. 33:325–330. 11. Cowman, A. F., and A. M. Lew. 1990. Chromosomal rearrangements and point mutations in the DHFR-TS gene of Plasmodium chabaudi under antifolate selection. Mol. Biochem. Parasitol. 42:21–30. 12. Dallas, W. S., J. E. Gowen, P. H. Ray, M. J. Cox, and I. K. Dev. 1992. Cloning, sequencing, and enhanced expression of the dihydropteroate synthase gene of Escherichia coli MC4100. J. Bacteriol. 174:5961–5970. 13. Holm, S. 1979. A simple sequentially rejective multiple test procedure. Scand. J. Statist. 6:65–70. 14. Jacobs, R. L. 1964. Role of p-aminobenzoic acid in Plasmodium berghei infection in the mouse. Exp. Parasitol. 15:213–225. 15. Kublin, J. G., F. K. Dzinjalamala, D. D. Kamwendo, E. M. Malkin, J. F. Cortese, L. M. Martino, R. A. G. Mukadam, S. J. Rogerson, A. G. Lescano,
26. 27.
28. 29.
30. 31.
32.
33. 34.
35.
36.
37. 38.
M. E. Molyneux, P. A. Winstanley, P. Chimpeni, T. E. Taylor, and C. V. Plowe. 2002. Molecular markers for failure of sulfadoxine-pyrimethamine and chlorproguanil-dapsone treatment of Plasmodium falciparum malaria. J. Infect. Dis. 185:380–388. Lander, E. S., and D. Botstein. 1989. Mapping Mendelian factors underlying quantitative traits using RFLP linkage maps. Genetics 121:185–199. Lopez, P., M. Espinosa, B. Greenberg, and S. A. Lacks. 1987. Sulfonamide resistance in Streptococcus pneumoniae: DNA sequence of the gene encoding dihydropteroate synthase and characterization of the enzyme. J. Bacteriol. 169:4320–4326. Lynch, M., and B. Walsh. 1998. Genetic analysis of quantitative traits. Sinauer Associates, Inc. Sunderland, Mass. Merkli, B., and R. W. Richle. 1980. Studies on the resistance to single and combined antimalarials in the Plasmodium berghei mouse model. Acta Trop. 37:228–231. Merkli, B., and R. Richle. 1983. Plasmodium berghei: diet and drug dosage regimes influencing selection of drug-resistant parasites in mice. Exp. Parasitol. 55:372–376. Ndounga, M., L. K. Basco, and P. Ringwald. 2001. Evaluation of a new sulfadoxine sensitivity assay in vitro for field isolates of Plasmodium falciparum. Trans. R. Soc. Trop. Med. Hyg. 95:55–57. Nzila, A. M., E. K. Mberu, J. Sulo, H. Dayo, P. A. Winstanley, C. H. Sibley, and W. M. Watkins. 2000. Towards an understanding of the mechanism of pyrimethamine-sulfadoxine resistance in Plasmodium falciparum: genotyping of dihydrofolate reductase and dihydropteroate synthase of Kenyan parasites. Antimicrob. Agents Chemother. 44:991–996. Omar, S. A., I. S. Adagu, and D. C. Warhurst. 2001. Can pretreatment screening for dhps and dhfr point mutations in Plasmodium falciparum infections be used to predict sulfadoxine-pyrimethamine treatment failure? Trans. R. Soc. Trop. Med. Hyg. 95:315–319. Pashley, T., F. Volpe, M. Pudney, J. E. Hyde, P. F. H. Sims, and C. J. Delves. 1997. Isolation and molecular characterization of the bifunctional hydroxymethyldihydropterin pyrophosphokinase-dihydropteroate synthase gene from Toxoplasma gondii. Mol. Biochem. Parasitol. 86:37–47. Peters, W., and B. L. Robinson. 1984. The chemotherapy of rodent malaria. XXXV. Further studies on the retardation of drug resistance by the use of a triple combination of mefloquine, pyrimethamine and sulfadoxine in mice infected with Plasmodium berghei and ⬘Plasmodium berghei NS.’ Ann. Trop. Med. Parasitol. 78:459–466. Peters, W. 1987. Chemotherapy and drug resistance in malaria, 2nd ed. Academic Press, London, United Kingdom. Peterson, D. S., D. Walliker, and T. E. Wellems. 1988. Evidence that a point mutation in dihydrofolate reductase-thymidylate synthase confers resistance to pyrimethamine in falciparum malaria. Proc. Natl. Acad. Sci. USA 85: 9114–9118. Plowe, C. V., J. G. Kublin, and O. K. Doumbo. 1998. Plasmodium falciparum dihydrofolate reductase and dihydropteroate synthase mutations: epidemiology and role in clinical resistance. Drug Resistance Updates 1:389–396. Sibley, C. H., J. E. Hyde, P. F. G. Sims, C. V. Plowe, J. G. Kublin, E. K. Mberu, A. F. Cowman, P. A. Winstanley, W. M. Watkins, and A. M. Nzila. 2001. Pyrimethamine-sulfadoxine resistance in Plasmodium falciparum: what next? Trends Parasitol. 17: 582–588. Siriwaraporn, W. 1998. Dihydrofolate reductase and antifolate resistance in malaria. Drug Resistance Updates 1:397–406. Slock, J., D. P. Stahly, C.-H. Han, E. W. Six, and I. P. Crawford. 1990. An apparent Bacillus subtilis folic acid biosynthetic operon containing pab, an amphibolic trpG gene, a third gene required for synthesis of para-aminobenzoic acid, and the dihydropteroate synthase gene. J. Bacteriol. 172:7211– 7226. Snounou, G., S. Viriyakosol, G. Holmquist, K. Berzins, P. Perlmann, and K. N. Brown. 1988. Cloning of genomic fragment from Plasmodium chabaudi expressing a 105 kilodalton antigen epitope. Mol. Biochem. Parasitol. 28: 153–162. Triglia, T., and A. F. Cowman. 1994. Primary structure and expression of the dihydropteroate synthase gene of Plasmodium falciparum. Proc. Natl. Acad. Sci. USA 91:7149–7153. Triglia, T., J. G. T. Menting, C. Wilson, and A. F. Cowman. 1997. Mutations in dihydropteroate synthase are responsible for sulfone and sulfonamide resistance in Plasmodium falciparum. Proc. Natl. Acad. Sci. USA 94:13944– 13949. Triglia, T., P. Wang, P. F. G. Sims, J. E. Hyde, and A. F. Cowman. 1998. Allelic exchange at the endogenous genomic locus in Plasmodium falciparum proves the role of dihydropteroate synthase in sulfadoxine-resistant malaria. EMBO J. 17:3807–3815. Volpe, F., M. Dyer, J. G. Scaife, G. Darby, D. K. Stammers, and C. J. Delves. 1992. The multifunctional folic acid synthesis fas gene of Pneumocystis carinii appears to encode dihydropteroate synthase and hydroxymethyldihydropterin pyrophosphokinase. Gene 112:213–218. Walliker, D., R. Carter, and S. Morgan. 1973. Genetic recombination in Plasmodium berghei. Parasitology 66:309–320. Walliker, D., R. Carter, and A. Sanderson. 1975. Genetic studies on Plas-
VOL. 46, 2002
39.
40.
41. 42.
SULFADOXINE-PYRIMETHAMINE RESISTANCE IN MALARIA
modium chabaudi: recombination between enzyme markers. Parasitology 70:19–24. Wang, P., C.-S. Lee, R. Bayoumi, A. Djimde, O. Doumbo, G. Swedberg, L. Duc Dao, H. Mshinda, M. Tanner, W. M. Watkins, P. F. G. Sims, and J. E. Hyde. 1997. Resistance to antifolates in Plasmodium falciparum monitored by sequence analysis of dihydropteroate synthase and dihydrofolate reductase alleles in a large number of field samples of diverse origins. Mol. Biochem. Parasitol. 89:161–177. Wang, P., M. Read, P. F. G. Sims, and J. E. Hyde. 1997. Sulfadoxine resistance in the human malaria parasite Plasmodium falciparum is determined by mutations in the dihydropteroate synthase and an additional factor associated with folate utilization. Mol. Microbiol. 23:979–986. Wang, P., P. F. G. Sims, and J. E. Hyde. 1997. A modified in vitro sulfadoxine susceptibility assay for Plasmodium falciparum suitable for investigating Fansidar resistance. Parasitology 115:223–230. Watkins, W. M., D. G. Sixsmith, J. D. Chulay, and H. C. Spencer. 1985.
2489
Antagonism of sulfadoxine and pyrimethamine antimalarial activity in vitro by p-aminobenzoic acid, p-aminobenzoylglutamic acid and folic acid. Mol. Biochem. Parasitol. 14:55–61. 43. Watkins, W. M., E. K. Mberu, P. A. Winstanley, and C. V. Plowe. 1997. The efficacy of antifolate antimalarial combinations in Africa: a predictive model based on pharmacodynamic and pharmacokinetic analyses. Parasitol. Today 13:459–464. 44. White, N. J., F. Nosten, S. Looareesuwan, W. M. Watkins, K. Marsh, R. W. Snow, G. Kokwaro, J. Ouma, T. T. Hien, M. E. Molyneux, T. E. Taylor, C. I. Newbold, T. K. Ruebush II, M. Danis, B. M. Greenwood, R. M. Anderson, and P. L. Olliaro. 1999. Averting a malaria disaster. Lancet 353:1965–1967. 45. Wu, Y., L. A. Kirkman, and T. E. Wellems. 1996. Transformation of Plasmodium falciparum malaria parasites by homologous integration of plasmids that confer resistance to pyrimethamine. Proc. Natl. Acad. Sci. USA 93: 1130–1134.