Analysis of Plasmodium falciparum diversity in natural infections by deep sequencing.

Malaria elimination strategies require surveillance of the parasite population for genetic changes that demand a public health response, such as new forms of drug resistance. Here we describe methods for the large-scale analysis of genetic variation in Plasmodium falciparum by deep sequencing of parasite DNA obtained from the blood of patients with malaria, either directly or after short-term culture. Analysis of 86,158 exonic single nucleotide polymorphisms that passed genotyping quality control in 227 samples from Africa, Asia and Oceania provides genome-wide estimates of allele frequency distribution, population structure and linkage disequilibrium. By comparing the genetic diversity of individual infections with that of the local parasite population, we derive a metric of within-host diversity that is related to the level of inbreeding in the population. An open-access web application has been established for the exploration of regional differences in allele frequency and of highly differentiated loci in the P. falciparum genome.

The MSPDBL2 codon 591 polymorphism is associated with lumefantrine in vitro drug responses in Plasmodium falciparum isolates from Kilifi, Kenya.

The mechanisms of drug resistance development in the Plasmodium falciparum parasite to lumefantrine (LUM), commonly used in combination with artemisinin, are still unclear. We assessed the polymorphisms of Pfmspdbl2 for associations with LUM activity in a Kenyan population. MSPDBL2 codon 591S was associated with reduced susceptibility to LUM (P = 0.04). The high frequency of Pfmspdbl2 codon 591S in Kenya may be driven by the widespread use of lumefantrine in artemisinin combination therapy (Coartem).

Repeat polymorphisms in the low-complexity regions of Plasmodium falciparum ABC transporters and associations with in vitro antimalarial responses.

The Plasmodium falciparum genome is rich in regions of low amino acid complexity which evolve with few constraints on size. To explore the extent of diversity in these loci, we sequenced repeat regions in pfmdr1, pfmdr5, pfmdr6, pfmrp2, and the antigenic locus pfmsp8 in laboratory and cultured-adapted clinical isolates. We further assessed associations between the repeats and parasite in vitro responses to 7 antimalarials to determine possible adaptive roles of these repeats in drug tolerance. Our results show extensive repeat variations in the reference and clinical isolates in all loci. We also observed a modest increase in dihydroartemisinin activity in parasites harboring the pfmdr1 sequence profile 7-2-10 (reflecting the number of asparagine repeats, number of aspartate repeats, and number of asparagine repeats in the final series of the gene product) (P = 0.0321) and reduced sensitivity to chloroquine, mefloquine, quinine, and dihydroartemisinin in those with the 7-2-11 profile (P = 0.0051, 0.0068, 0.0011, and 0.0052, respectively). Interestingly, we noted an inverse association between two drugs whereby isolates with 6 asparagine repeats encoded by pfmdr6 were significantly more susceptible to piperaquine than those with 8 (P = 0.0057). Against lumefantrine, those with 8 repeats were, however, more sensitive (P = 0.0144). In pfmrp2, the 7-DNNNTS/NNNNTS (number of DNNNTS or NNNNTS motifs; underlining indicates dimorphism) repeat group was significantly associated with a higher lumefantrine 50% inhibitory concentration (IC50) (P = 0.008) than in those without. No associations were observed with pfmsp8. These results hint at the probable utility of some repeat conformations as markers of in vitro antimalarial response; hence, biochemical functional studies to ascertain their role in P. falciparum are required.

Antimalarial activity of isoquine against Kenyan Plasmodium falciparum clinical isolates and association with polymorphisms in pfcrt and pfmdr1 genes.

BACKGROUND: The use of amodiaquine in prophylaxis is associated with serious toxicity, resulting from its metabolic conversion into a reactive quinone-imine metabolite by the hepatic cytochrome P450. To circumvent this toxicity, several amodiaquine analogues that lack the potential to form a quinone-imine derivative, while retaining antimalarial activity, have been designed. Isoquine is one of these promising molecules that has already reached Phase I clinical trials in humans. METHODS: We analysed the in vitro activity of isoquine against 62 Plasmodium falciparum isolates collected in Kenya and the association of this activity with polymorphisms in pfcrt and pfmdr1 genes. RESULTS: The median concentration of isoquine that inhibited 50% of parasite growth (IC50) was 9 nM, compared with 56 nM chloroquine, 8 nM amodiaquine, 10 nM desethylamodiaquine, 69 nM lumefantrine and 1 nM dihydroartemisinin. Isoquine activity was correlated with polymorphisms in pfcrt at codon 76, but not in pfmdr1 at codon 86. CONCLUSIONS: The high activity of isoquine against field isolates, including chloroquine-resistant isolates, with IC50 <10 nM, warrants its further development as an antimalarial.

Baseline in vitro activities of the antimalarials pyronaridine and methylene blue against Plasmodium falciparum isolates from Kenya.

We have analyzed the in vitro activities of pyronaridine and methylene blue against 59 Plasmodium falciparum isolates from Kenya in association with polymorphisms in Pfcrt (codon 76), Pfmdr1 (codon 86), and Pfnhe (full sequence). The median inhibitory concentrations that kill 50% of parasites were 13.5 and 3.3 nM for pyronaridine and methylene blue, respectively. Their activities were not associated with polymorphisms in these genes. The drugs' high in vitro activities indicate that they would be efficacious against Kenyan isolates in vivo.

In vitro activities of quinine and other antimalarials and pfnhe polymorphisms in Plasmodium isolates from Kenya.

Resistance to the amino alcohol quinine has been associated with polymorphisms in pfnhe, a sodium hydrogen exchanger. We investigated the role of this gene in quinine resistance in vitro in isolates from Kenya. We analyzed pfnhe whole-gene polymorphisms, using capillary sequencing, and pfcrt at codon 76 (pfcrt-76) and pfmdr1 at codon 86 (pfmdr1-86), using PCR-enzyme restriction methodology, in 29 isolates from Kilifi, Kenya, for association with the in vitro activities of quinine and 2 amino alcohols, mefloquine and halofantrine. In vitro activity was assessed as the drug concentration that inhibits 50% of parasite growth (IC50). The median IC50s of quinine, halofantrine, and mefloquine were 92, 22, and 18 nM, respectively. The presence of 2 DNNND repeats in microsatellite ms4760 of pfnhe was associated with reduced susceptibility to quinine (60 versus 227 nM for 1 and 2 repeats, respectively; P<0.05), while 3 repeats were associated with restoration of susceptibility. The decrease in susceptibility conferred by the 2 DNNND repeats was more pronounced in parasites harboring the pfmdr1-86 mutation. No association was found between susceptibility to quinine and the pfcrt-76 mutation or between susceptibility to mefloquine or halofantrine and the pfnhe gene and the pfcrt-76 and pfmdr1-86 mutations. Using previously published data on the in vitro activities of chloroquine, lumefantrine, piperaquine, and dihydroartemisinin, we investigated the association of their activities with pfnhe polymorphism. With the exception of a modulation of the activity of lumefantrine by a mutation at position 1437, pfnhe did not modulate their activities. Two DNNND repeats combined with the pfmdr1-86 mutation could be used as an indicator of reduced susceptibility to quinine.

In vitro activity of antifolate and polymorphism in dihydrofolate reductase of Plasmodium falciparum isolates from the Kenyan coast: emergence of parasites with Ile-164-Leu mutation.

We have analyzed the activities of the antifolates pyrimethamine (PM), chlorcycloguanil (CCG), WR99210, trimethoprim (TMP), methotrexate (MTX), and trimetrexate (TMX) against Kenyan Plasmodium falciparum isolates adapted in vitro for long-term culture. We have also assessed the relationship between these drug activities and mutations in dihydrofolate reductase (dhfr), a domain of the gene associated with antifolate resistance. As expected, WR99210 was the most potent drug, with a median 50% inhibitory concentration (IC50) of <0.075 nM, followed by TMX, with a median IC50 of 30 nM. The median IC50 of CCG was 37.80 nM, and that of MTX was 83.60 nM. PM and TMP were the least active drugs, with median IC50s of 733.26 nM and 29,656.04 nM, respectively. We analyzed parasite dhfr genotypes by the PCR-enzyme restriction technique. No wild-type dhfr parasite was found. Twenty-four of 33 parasites were triple mutants (mutations at codons 108, 51, and 59), and only 8/33 were double mutants (mutations at codons 108 and 51 or at codons 108 and 59). IC50s were 2.1-fold (PM) and 3.6-fold (TMP) higher in triple than in double mutants, though these differences were not statistically significant. Interestingly, we have identified a parasite harboring a mutation at codon 164 (Ile-164-Leu) in addition to mutations at codons 108, 51, and 59. This quadruple mutant parasite had the highest TMP IC50 and was in the upper 10th percentile against PM and CCG. We confirmed the presence of this mutation by sequencing. Thus, TMX and MTX are potent against P. falciparum, and quadruple mutants are now emerging in Africa.

In vitro activities of piperaquine, lumefantrine, and dihydroartemisinin in Kenyan Plasmodium falciparum isolates and polymorphisms in pfcrt and pfmdr1.

We have analyzed the in vitro chemosensitivity profiles of 115 Kenyan isolates for chloroquine (CQ), piperaquine, lumefantrine (LM), and dihydroartemisinin in association with polymorphisms in pfcrt at codon 76 and pfmdr1 at codon 86, as well as with variations of the copy number of pfmdr1. The median drug concentrations that inhibit 50% of parasite growth (IC(50)s) were 41 nM (interquartile range [IQR], 18 to 73 nM), 50 nM (IQR, 29 to 96 nM), 32 nM (IQR, 17 to 46 nM), and 2 nM (IQR, 1 to 3 nM) for CQ, LM, piperaquine, and dihydroartemisinin, respectively. The activity of CQ correlated inversely with that of LM (r(2) = -0.26; P = 0.02). Interestingly, parasites for which LM IC(50)s were higher were wild type for pfcrt-76 and pfmdr1-86. All isolates had one pfmdr1 copy. Thus, the decrease in LM activity is associated with the selection of wild-type pfcrt-76 and pfmdr1-86 parasites, a feature that accounts for the inverse relationship between CQ and LM. Therefore, the use of LM-artemether is likely to lead to the selection of more CQ-susceptible parasites.

Chloroquine resistance before and after its withdrawal in Kenya.

BACKGROUND: The spread of resistance to chloroquine (CQ) led to its withdrawal from use in most countries in sub-Saharan Africa in the 1990s. In Malawi, this withdrawal was followed by a rapid reduction in the frequency of resistance to the point where the drug is now considered to be effective once again, just nine years after its withdrawal. In this report, the polymorphisms of markers associated with CQ-resistance against Plasmodium falciparum isolates from coastal Kenya (Kilifi) were investigated, from 1993, prior to the withdrawal of CQ, to 2006, seven years after its withdrawal. Changes to those that occurred in the dihydrofolate reductase gene (dhfr) that confers resistance to the replacement drug, pyrimethamine/sulphadoxine were also compared. METHODS: Mutations associated with CQ resistance, at codons 76 of pfcrt, at 86 of pfmdr1, and at codons 51, 59 and 164 of dhfr were analysed using PCR-restriction enzyme methods. In total, 406, 240 and 323 isolates were genotyped for pfcrt-76, pfmdr1-86 and dhfr, respectively. RESULTS: From 1993 to 2006, the frequency of the pfcrt-76 mutant significantly decreased from around 95% to 60%, while the frequency of pfmdr1-86 did not decline, remaining around 75%. Though the frequency of dhfr mutants was already high (around 80%) at the start of the study, this frequency increased to above 95% during the study period. Mutation at codon 164 of dhfr was analysed in 2006 samples, and none of them had this mutation. CONCLUSION: In accord with the study in Malawi, a reduction in resistance to CQ following official withdrawal in 1999 was found, but unlike Malawi, the decline of resistance to CQ in Kilifi was much slower. It is estimated that, at current rates of decline, it will take 13 more years for the clinical efficacy of CQ to be restored in Kilifi. In addition, CQ resistance was declining before the drug's official withdrawal, suggesting that, prior to the official ban, the use of CQ had decreased, probably due to its poor clinical effectiveness.

A genome wide association study of Plasmodium falciparum susceptibility to 22 antimalarial drugs in Kenya.

BACKGROUND: Drug resistance remains a chief concern for malaria control. In order to determine the genetic markers of drug resistant parasites, we tested the genome-wide associations (GWA) of sequence-based genotypes from 35 Kenyan P. falciparum parasites with the activities of 22 antimalarial drugs. METHODS AND PRINCIPAL FINDINGS: Parasites isolated from children with acute febrile malaria were adapted to culture, and sensitivity was determined by in vitro growth in the presence of anti-malarial drugs. Parasites were genotyped using whole genome sequencing techniques. Associations between 6250 single nucleotide polymorphisms (SNPs) and resistance to individual anti-malarial agents were determined, with false discovery rate adjustment for multiple hypothesis testing. We identified expected associations in the pfcrt region with chloroquine (CQ) activity, and other novel loci associated with amodiaquine, quinazoline, and quinine activities. Signals for CQ and primaquine (PQ) overlap in and around pfcrt, and interestingly the phenotypes are inversely related for these two drugs. We catalog the variation in dhfr, dhps, mdr1, nhe, and crt, including novel SNPs, and confirm the presence of a dhfr-164L quadruple mutant in coastal Kenya. Mutations implicated in sulfadoxine-pyrimethamine resistance are at or near fixation in this sample set. CONCLUSIONS/SIGNIFICANCE: Sequence-based GWA studies are powerful tools for phenotypic association tests. Using this approach on falciparum parasites from coastal Kenya we identified known and previously unreported genes associated with phenotypic resistance to anti-malarial drugs, and observe in high-resolution haplotype visualizations a possible signature of an inverse selective relationship between CQ and PQ.

Population genetic analysis of Plasmodium falciparum parasites using a customized Illumina GoldenGate genotyping assay.

The diversity in the Plasmodium falciparum genome can be used to explore parasite population dynamics, with practical applications to malaria control. The ability to identify the geographic origin and trace the migratory patterns of parasites with clinically important phenotypes such as drug resistance is particularly relevant. With increasing single-nucleotide polymorphism (SNP) discovery from ongoing Plasmodium genome sequencing projects, a demand for high SNP and sample throughput genotyping platforms for large-scale population genetic studies is required. Low parasitaemias and multiple clone infections present a number of challenges to genotyping P. falciparum. We addressed some of these issues using a custom 384-SNP Illumina GoldenGate assay on P. falciparum DNA from laboratory clones (long-term cultured adapted parasite clones), short-term cultured parasite isolates and clinical (non-cultured isolates) samples from East and West Africa, Southeast Asia and Oceania. Eighty percent of the SNPs (n = 306) produced reliable genotype calls on samples containing as little as 2 ng of total genomic DNA and on whole genome amplified DNA. Analysis of artificial mixtures of laboratory clones demonstrated high genotype calling specificity and moderate sensitivity to call minor frequency alleles. Clear resolution of geographically distinct populations was demonstrated using Principal Components Analysis (PCA), and global patterns of population genetic diversity were consistent with previous reports. These results validate the utility of the platform in performing population genetic studies of P. falciparum.

Characteristics of Plasmodium falciparum dhfr haplotypes that confer pyrimethamine resistance, Kilifi, Kenya, 1987--2006.

Resistance to the antimalarial drug sulfadoxine-pyrimethamine (SP) emerged in Plasmodium falciparum from Asia in the 1960s and subsequently spread to Africa. It is not known whether alleles that confer SP resistance also arose independently in Africa. We defined the coding region and microsatellite haplotypes of dhfr alleles in P. falciparum collected in Kilifi, Kenya, during 1987--2006, which spans the period when SP was first introduced. Isolates that carried a double-mutant or triple-mutant dhfr allele were detected at a low frequency, even during 1987--1988. Each of 2 double mutants carried a unique haplotype, and both were related to wild-type haplotypes from the same population. The number of isolates that carried a triple-mutant dhfr allele increased rapidly after introduction of SP and shared the haplotype of the triple mutant derived form Asia. We observed no triple-mutant alleles with haplotypes related to those of the Africa-derived wild-type and double-mutant alleles.