Plasmodium falciparum resistance to artemisinin derivatives in southeast Asia threatens malaria control and elimination activities worldwide. To monitor the spread of artemisinin resistance, a molecular marker is urgently needed. Here, using whole-genome sequencing of an artemisinin-resistant parasite line from Africa and clinical parasite isolates from Cambodia, we associate mutations in the PF3D7_1343700 kelch propeller domain (‘K13-propeller’) with artemisinin resistance in vitro and in vivo. Mutant K13-propeller alleles cluster in Cambodian provinces where resistance is prevalent, and the increasing frequency of a dominant mutant K13-propeller allele correlates with the recent spread of resistance in western Cambodia. Strong correlations between the presence of a mutant allele, in vitro parasite survival rates and in vivo parasite clearance rates indicate that K13-propeller mutations are important determinants of artemisinin resistance. K13-propeller polymorphism constitutes a useful molecular marker for large-scale surveillance efforts to contain artemisinin resistance in the Greater Mekong Subregion and prevent its global spread.

The emergence of artemisinin resistance in Southeast Asia imperils efforts to reduce the global malaria burden. We genetically modified the Plasmodium falciparum K13 locus using zinc-finger nucleases and measured ring-stage survival rates after drug exposure in vitro; these rates correlate with parasite clearance half-lives in artemisinin-treated patients. With isolates from Cambodia, where resistance first emerged, survival rates decreased from 13 to 49% to 0.3 to 2.4% after the removal of K13 mutations. Conversely, survival rates in wild-type parasites increased from ≤0.6% to 2 to 29% after the insertion of K13 mutations. These mutations conferred elevated resistance to recent Cambodian isolates compared with that of reference lines, suggesting a contemporary contribution of additional genetic factors. Our data provide a conclusive rationale for worldwide K13-propeller sequencing to identify and eliminate artemisinin-resistant parasites.

We report a large multicenter genome-wide association study of Plasmodium falciparum resistance to artemisinin, the frontline antimalarial drug. Across 15 locations in Southeast Asia, we identified at least 20 mutations in kelch13 (PF3D7_1343700) affecting the encoded propeller and BTB/POZ domains, which were associated with a slow parasite clearance rate after treatment with artemisinin derivatives. Nonsynonymous polymorphisms in fd (ferredoxin), arps10 (apicoplast ribosomal protein S10), mdr2 (multidrug resistance protein 2) and crt (chloroquine resistance transporter) also showed strong associations with artemisinin resistance. Analysis of the fine structure of the parasite population showed that the fd, arps10, mdr2 and crt polymorphisms are markers of a genetic background on which kelch13 mutations are particularly likely to arise and that they correlate with the contemporary geographical boundaries and population frequencies of artemisinin resistance. These findings indicate that the risk of new resistance-causing mutations emerging is determined by specific predisposing genetic factors in the underlying parasite population.

Summary Background Artemisinin resistance in Plasmodium falciparum lengthens parasite clearance half-life during artemisinin monotherapy or artemisinin-based combination therapy. Absence of in-vitro and ex-vivo correlates of artemisinin resistance hinders study of this phenotype. We aimed to assess whether an in-vitro ring-stage survival assay (RSA) can identify culture-adapted P falciparum isolates from patients with slow-clearing or fast-clearing infections, to investigate the stage-dependent susceptibility of parasites to dihydroartemisinin in the in-vitro RSA, and to assess whether an ex-vivo RSA can identify artemisinin-resistant P falciparum infections. Methods We culture-adapted parasites from patients with long and short parasite clearance half-lives from a study done in Pursat, Cambodia, in 2010 (registered with ClinicalTrials.gov, number NCT00341003) and used novel in-vitro survival assays to explore the stage-dependent susceptibility of slow-clearing and fast-clearing parasites to dihydroartemisinin. In 2012, we implemented the RSA in prospective parasite clearance studies in Pursat, Preah Vihear, and Ratanakiri, Cambodia (NCT01736319), to measure the ex-vivo responses of parasites from patients with malaria. Continuous variables were compared with the Mann-Whitney U test. Correlations were analysed with the Spearman correlation test. Findings In-vitro survival rates of culture-adapted parasites from 13 slow-clearing and 13 fast-clearing infections differed significantly when assays were done on 0–3 h ring-stage parasites (10·88% vs 0·23%; p=0·007). Ex-vivo survival rates significantly correlated with in-vivo parasite clearance half-lives (n=30, r=0·74, 95% CI 0·50–0·87; p<0·0001). Interpretation The in-vitro RSA of 0–3 h ring-stage parasites provides a platform for the molecular characterisation of artemisinin resistance. The ex-vivo RSA can be easily implemented where surveillance for artemisinin resistance is needed. Funding Institut Pasteur du Cambodge and the Intramural Research Program, NIAID, NIH.

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.

Artemisinin resistance in Plasmodium falciparum threatens global efforts to control and eliminate malaria. Polymorphisms in the kelch domain–carrying protein K13 are associated with artemisinin resistance, but the underlying molecular mechanisms are unknown. We analyzed the in vivo transcriptomes of 1043 P. falciparum isolates from patients with acute malaria and found that artemisinin resistance is associated with increased expression of unfolded protein response (UPR) pathways involving the major PROSC and TRiC chaperone complexes. Artemisinin-resistant parasites also exhibit decelerated progression through the first part of the asexual intraerythrocytic development cycle. These findings suggest that artemisinin-resistant parasites remain in a state of decelerated development at the young ring stage, whereas their up-regulated UPR pathways mitigate protein damage caused by artemisinin. The expression profiles of UPR-related genes also associate with the geographical origin of parasite isolates, further suggesting their role in emerging artemisinin resistance in the Greater Mekong Subregion.

We describe an analysis of genome variation in 825 Plasmodium falciparum samples from Asia and Africa that reveals an unusual pattern of parasite population structure at the epicentre of artemisinin resistance in western Cambodia. Within this relatively small geographical area we have discovered several distinct but apparently sympatric parasite subpopulations with extremely high levels of genetic differentiation. Of particular interest are three subpopulations, all associated with clinical resistance to artemisinin, which have skewed allele frequency spectra and remarkably high levels of haplotype homozygosity, indicative of founder effects and recent population expansion. We provide a catalogue of SNPs that show high levels of differentiation in the artemisinin-resistant subpopulations, including codon variants in various transporter proteins and DNA mismatch repair proteins. These data provide a population genetic framework for investigating the biological origins of artemisinin resistance and for defining molecular markers to assist its elimination.

Summary Background As the prevalence of artemisinin-resistant Plasmodium falciparum malaria increases in the Greater Mekong Subregion (GMS), emerging resistance to partner drugs in artemisinin combination therapies (ACTs) seriously threatens global efforts to treat and eliminate this disease. Molecular markers for ACT failure are urgently needed to monitor the spread of partner drug resistance, and to recommend alternative treatments in Southeast Asia and beyond. Methods We performed a genome-wide association study (GWAS) of 297 P. falciparum isolates from Cambodia to investigate the relationship of 11,630 exonic single-nucleotide polymorphisms (SNPs) and 43 copy number variations (CNVs) with in-vitro piperaquine 50% inhibitory concentrations (IC50s), and tested whether these genetic variants are markers of dihydroartemisinin-piperaquine failures. We then performed a survival analysis of 133 patients to determine whether candidate molecular markers predicted parasite recrudescence following dihydroartemisinin-piperaquine treatment. Findings Piperaquine IC50s increased significantly from 2011 to 2013 in 3 Cambodian provinces. Genome-wide analysis of SNPs identified a chromosome 13 region that associates with elevated piperaquine IC50s. A nonsynonymous SNP (encoding a Glu415Gly substitution) in this region, within a gene encoding an exonuclease, associates with parasite recrudescence following dihydroartemisinin-piperaquine treatment. Genome-wide analysis of CNVs revealed that a single copy of the mdr1 gene on chromosome 5 and a novel amplification of the plasmepsin II and plasmepsin III genes on chromosome 14 also associate with elevated piperaquine IC50s. After adjusting for covariates, both exo-E415G and plasmepsin II-III markers significantly associate with decreased treatment efficacy (0.38 and 0.41 survival rates, respectively). Interpretation The exo-E415G SNP and plasmepsin II-III amplification are markers of piperaquine resistance and dihydroartemisinin-piperaquine failures in Cambodia, and can help monitor the spread of these phenotypes into GMS countries, and elucidate the mechanism of piperaquine resistance. Since plasmepsins are involved in the parasite’s haemoglobin-to-haemozoin conversion pathway, targeted by related antimalarials, plasmepsin II-III amplification likely mediates piperaquine resistance. Funding Intramural Research Program of the US National Institute of Allergy and Infectious Diseases, National Institutes of Health; Wellcome Trust; Bill and Melinda Gates Foundation; Medical Research Council; and UK Department for International Development.