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.

BACKGROUND Recent gains in reducing the global burden of malaria are threatened by the emergence of Plasmodium falciparum resistance to artemisinins. The discovery that mutations in portions of a P. falciparum gene encoding kelch (K13)–propeller domains are the major determinant of resistance has provided opportunities for monitoring such resistance on a global scale. METHODS We analyzed the K13-propeller sequence polymorphism in 14,037 samples collected in 59 countries in which malaria is endemic. Most of the samples (84.5%) were obtained from patients who were treated at sentinel sites used for nationwide surveillance of antimalarial resistance. We evaluated the emergence and dissemination of mutations by haplotyping neighboring loci. RESULTS We identified 108 nonsynonymous K13 mutations, which showed marked geographic disparity in their frequency and distribution. In Asia, 36.5% of the K13 mutations were distributed within two areas — one in Cambodia, Vietnam, and Laos and the other in western Thailand, Myanmar, and China — with no overlap. In Africa, we observed a broad array of rare nonsynonymous mutations that were not associated with delayed parasite clearance. The gene-edited Dd2 transgenic line with the A578S mutation, which expresses the most frequently observed African allele, was found to be susceptible to artemisinin in vitro on a ring-stage survival assay. CONCLUSIONS No evidence of artemisinin resistance was found outside Southeast Asia and China, where resistance-associated K13 mutations were confined. The common African A578S allele was not associated with clinical or in vitro resistance to artemisinin, and many African mutations appear to be neutral.

SummaryBackgroundWestern Cambodia is the epicentre of Plasmodium falciparum multidrug resistance and is facing high rates of dihydroartemisinin–piperaquine treatment failures. Genetic tools to detect the multidrug-resistant parasites are needed. Artemisinin resistance can be tracked using the K13 molecular marker, but no marker exists for piperaquine resistance. We aimed to identify genetic markers of piperaquine resistance and study their association with dihydroartemisinin–piperaquine treatment failures.MethodsWe obtained blood samples from Cambodian patients infected with P falciparum and treated with dihydroartemisinin–piperaquine. Patients were followed up for 42 days during the years 2009–15. We established in-vitro and ex-vivo susceptibility profiles for a subset using piperaquine survival assays. We determined whole-genome sequences by Illumina paired-reads sequencing, copy number variations by qPCR, RNA concentrations by qRT-PCR, and protein concentrations by immunoblotting. Fisher’s exact and non-parametric Wilcoxon rank-sum tests were used to identify significant differences in single-nucleotide polymorphisms or copy number variants, respectively, for differential distribution between piperaquine-resistant and piperaquine-sensitive parasite lines.FindingsWhole-genome exon sequence analysis of 31 culture-adapted parasite lines associated amplification of the plasmepsin 2–plasmepsin 3 gene cluster with in-vitro piperaquine resistance. Ex-vivo piperaquine survival assay profiles of 134 isolates correlated with plasmepsin 2 gene copy number. In 725 patients treated with dihydroartemisinin–piperaquine, multicopy plasmepsin 2 in the sample collected before treatment was associated with an adjusted hazard ratio (aHR) for treatment failure of 20·4 (95% CI 9·1–45·5, p<0·0001). Multicopy plasmepsin 2 predicted dihydroartemisinin–piperaquine failures with 0·94 (95% CI 0·88–0·98) sensitivity and 0·77 (0·74–0·81) specificity. Analysis of samples collected across the country from 2002 to 2015 showed that the geographical and temporal increase of the proportion of multicopy plasmepsin 2 parasites was highly correlated with increasing dihydroartemisinin–piperaquine treatment failure rates (r=0·89 [95% CI 0·77–0·95], p<0·0001, Spearman’s coefficient of rank correlation). Dihydroartemisinin–piperaquine efficacy at day 42 fell below 90% when the proportion of multicopy plasmepsin 2 parasites exceeded 22%.InterpretationPiperaquine resistance in Cambodia is strongly associated with amplification of plasmepsin 2–3, encoding haemoglobin-digesting proteases, regardless of the location. Multicopy plasmepsin 2 constitutes a surrogate molecular marker to track piperaquine resistance. A molecular toolkit combining plasmepsin 2 with K13 and mdr1 monitoring should provide timely information for antimalarial treatment and containment policies.FundingInstitut Pasteur in Cambodia, Institut Pasteur Paris, National Institutes of Health, WHO, Agence Nationale de la Recherche, Investissement d’Avenir programme, Laboratoire d’Excellence In...

fThe declining efficacy of artemisinin derivatives against Plasmodium falciparum in western Cambodia is a major concern. The knowledge gap in the understanding of the mechanisms involved hampers designing monitoring tools. Here, we culture-adapted 20 isolates from Pailin and Ratanakiri (areas of artemisinin resistance and susceptibility in western and eastern Cambodia, respectively) and studied their in vitro response to dihydroartemisinin. No significant difference between the two sets of isolates was observed in the classical isotopic test. However, a 6-h pulse exposure to 700 nM dihydroartemisinin (ring-stage survival assay -RSA]) revealed a clear-cut geographic dichotomy. The survival rate of exposed ring-stage parasites (ring stages) was 17-fold higher in isolates from Pailin (median, 13.5%) than in those from Ratanakiri (median, 0.8%), while exposed mature stages were equally and highly susceptible (0.6% and 0.7%, respectively). Ring stages survived drug exposure by cell cycle arrest and resumed growth upon drug withdrawal. The reduced susceptibility to artemisinin in Pailin appears to be associated with an altered in vitro phenotype of ring stages from Pailin in the RSA.