Steven M. Kiara

Next-generation sequencing is used here to analyse Plasmodium falciparum genome variation directly from clinical blood samples, as well as cultured isolates, from Africa, Asia and Oceania. Resistance to the major antimalarial drug artemisinin is emerging in the Plasmodium falciparum parasite across Southeast Asia, and there is concern that the increased deployment of antimalarials in pursuit of disease eradication might simply lead to increased drug resistance. To monitor these risks it is important to survey the parasite population for genetic changes. Next-generation sequencing is used here to analyse P. falciparum genome variation directly from nearly 300 clinical blood samples, and from cultured isolates from Africa, Asia and Oceania. The authors use these data to analyse the diversity of the parasite population across different geographical locations, as well as within-host diversity at the level of the whole genome, and they show how this may be used to estimate inbreeding rates, which are important for the evolution of drug resistance. Malaria elimination strategies require surveillance of the parasite population for genetic changes that demand a public health response, such as new forms of drug resistance1,2. 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.

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.

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.

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.

The folate derivatives folic acid (FA) and folinic acid (FNA) decrease the in vivo and in vitro activities of antifolate drugs in Plasmodium falciparum. However, the effects of 5-methyl-tetrahydrofolate (5-Me-THF) and tetrahydrofolate (THF), the two dominant circulating folate forms in humans, have not been explored yet. We have investigated the effects of FA, FNA, 5-Me-THF, and THF on the in vitro activity of the antimalarial antifolates pyrimethamine and chlorcycloguanil and the anticancer antifolates methotrexate (MTX), aminopterin, and trimetrexate (TMX), against P. falciparum. The results indicate that these anticancers are potent against P. falciparum, with IC50 < 50 nM. 5-Me-THF does not significantly decrease the activity of all tested drugs, and none of the tested folate derivatives significantly decrease the activity of these anticancers. Thus, malaria folate metabolism has features different from those in human, and the exploitation of this difference could lead to the discovery of new drugs to treat malaria. For instance, the combination of 5-Me-THF with a low dose of TMX could be used to treat malaria. In addition, the safety of a low dose of MTX in the treatment of arthritis indicates that this drug could be used alone to treat malaria.