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. 1,2 Here we describe methods for 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 SNPs 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 exploration of regional differences in allele frequency and of highly differentiated loci in the P. falciparum genome.

Summary Background Artemisinin-resistant Plasmodium falciparum has been reported in Pailin, western Cambodia, detected as a slow parasite clearance rate in vivo. Emergence of this phenotype in western Thailand and possibly elsewhere threatens to compromise the effectiveness of all artemisinin-based combination therapies. Parasite genetics is associated with parasite clearance rate but does not account for all variation. We investigated contributions of both parasite genetics and host factors to the artemisinin-resistance phenotype in Pursat, western Cambodia. Methods Between June 19 and Nov 28, 2009, and June 26 and Dec 6, 2010, we enrolled patients aged 10 years or older with uncomplicated falciparum malaria, a density of asexual parasites of at least 10 000 per μL of whole blood, no symptoms or signs of severe malaria, no other cause of febrile illness, and no chronic illness. We gave participants 4 mg/kg artesunate at 0, 24, and 48 h, 15 mg/kg mefloquine at 72 h, and 10 mg/kg mefloquine at 96 h. We assessed parasite density on thick blood films every 6 h until undetectable. The parasite clearance half-life was calculated from the parasite clearance curve. We genotyped parasites with 18 microsatellite markers and patients for haemoglobin E, α-thalassaemia, and a mutation of G6PD, which encodes glucose-6-phosphate dehydrogenase. To account for the possible effects of acquired immunity on half-life, we used three surrogates for increased likelihood of exposure to P falciparum: age, sex, and place of residence. This study is registered with ClinicalTrials.gov, number NCT00341003. Findings We assessed 3504 individuals from all six districts of Pursat province seeking treatment for malaria symptoms. We enrolled 168 patients with falciparum malaria who met inclusion criteria. The geometric mean half-life was 5.85 h (95% CI 5.54–6.18) in Pursat, similar to that reported in Pailin (p=0.109). We identified two genetically different parasite clone groups: parasite group 1 (PG1) and parasite group 2 (PG2). Non-significant increases in parasite clearance half-life were seen in patients with haemoglobin E (0.55 h; p=0.078), those of male sex (0.96 h; p=0.064), and in 2010 (0.68 h; p=0.068); PG1 was associated with a significant increase (0.79 h; p=0.033). The mean parasite heritability of half-life was 0.40 (SD 0.17). Interpretation Heritable artemisinin resistance is established in a second Cambodian province. To accurately identify parasites that are intrinsically susceptible or resistant to artemisinins, future studies should explore the effect of erythrocyte polymorphisms and specific immune responses on half-life variation. Funding Division of Intramural Research, National Institute of Allergy and Infectious Diseases, National Institutes of Health.

Antimalarial drugs impose strong pressure on Plasmodium falciparum parasites and leave signatures of selection in the parasite genome 1,2. Search for signals of selection may lead to genes encoding drug or immune targets 3. The lack of high-throughput genotyping methods, inadequate knowledge of parasite population history, and time-consuming adaptations of parasites to in vitro culture have hampered genome-wide association studies (GWAS) of parasite traits. Here we report genotyping of DNA from 189 culture-adapted P. falciparum parasites using a custom-built array with thousands of single nucleotide polymorphisms (SNPs). Population structure, variation in recombination rate, and loci under recent positive selection were detected. Parasite half maximum inhibitory concentrations (IC50) to seven antimalarial drugs were obtained and used in GWAS to identify genes associated with drug responses. The SNP array and genome-wide parameters provide valuable tools and information for new advances in P. falciparum genetics.