SummaryBackgroundEmergence of artemisinin resistance in southeast Asia poses a serious threat to the global control of Plasmodium falciparum malaria. Discovery of the K13 marker has transformed approaches to the monitoring of artemisinin resistance, allowing introduction of molecular surveillance in remote areas through analysis of DNA. We aimed to assess the spread of artemisinin-resistant P falciparum in Myanmar by determining the relative prevalence of P falciparum parasites carrying K13-propeller mutations.MethodsWe did this cross-sectional survey at malaria treatment centres at 55 sites in ten administrative regions in Myanmar, and in relevant border regions in Thailand and Bangladesh, between January, 2013, and September, 2014. K13 sequences from P falciparum infections were obtained mainly by passive case detection. We entered data into two geostatistical models to produce predictive maps of the estimated prevalence of mutations of the K13 propeller region across Myanmar.FindingsOverall, 371 (39%) of 940 samples carried a K13-propeller mutation. We recorded 26 different mutations, including nine mutations not described previously in southeast Asia. In seven (70%) of the ten administrative regions of Myanmar, the combined K13-mutation prevalence was more than 20%. Geospatial mapping showed that the overall prevalence of K13 mutations exceeded 10% in much of the east and north of the country. In Homalin, Sagaing Region, 25 km from the Indian border, 21 (47%) of 45 parasite samples carried K13-propeller mutations.InterpretationArtemisinin resistance extends across much of Myanmar. We recorded P falciparum parasites carrying K13-propeller mutations at high prevalence next to the northwestern border with India. Appropriate therapeutic regimens should be tested urgently and implemented comprehensively if spread of artemisinin resistance to other regions is to be avoided.FundingWellcome Trust–Mahidol University–Oxford Tropical Medicine Research Programme and the Bill & Melinda Gates Foundation.

BackgroundArtemisinin-based combination therapies are the first line of treatment for Plasmodium falciparum infections worldwide, but artemisinin resistance has risen rapidly in Southeast Asia over the past decade. Mutations in the kelch13 gene have been implicated in this resistance. We used longitudinal genomic surveillance to detect signals in kelch13 and other loci that contribute to artemisinin or partner drug resistance. We retrospectively sequenced the genomes of 194 P. falciparum isolates from five sites in Northwest Thailand, over the period of a rapid increase in the emergence of artemisinin resistance (2001–2014).ResultsWe evaluate statistical metrics for temporal change in the frequency of individual SNPs, assuming that SNPs associated with resistance increase in frequency over this period. After Kelch13-C580Y, the strongest temporal change is seen at a SNP in phosphatidylinositol 4-kinase, which is involved in a pathway recently implicated in artemisinin resistance. Furthermore, other loci exhibit strong temporal signatures which warrant further investigation for involvement in artemisinin resistance evolution. Through genome-wide association analysis we identify a variant in a kelch domain-containing gene on chromosome 10 that may epistatically modulate artemisinin resistance.ConclusionsThis analysis demonstrates the potential of a longitudinal genomic surveillance approach to detect resistance-associated gene loci to improve our mechanistic understanding of how resistance develops. Evidence for additional genomic regions outside of the kelch13 locus associated with artemisinin-resistant parasites may yield new molecular markers for resistance surveillance, which may be useful in efforts to reduce the emergence or spread of artemisinin resistance in African parasite populations.Electronic supplementary materialThe online version of this article (doi:10.1186/s13059-017-1204-4) contains supplementary material, which is available to authorized users.

Multiple kelch13 alleles conferring artemisinin resistance (ART-R) are currently spreading through Southeast Asian malaria parasite populations, providing a unique opportunity to observe an ongoing soft selective sweep, investigate why resistance alleles have evolved multiple times and determine fundamental population genetic parameters for Plasmodium. We sequenced kelch13 (n = 1,876), genotyped 75 flanking SNPs, and measured clearance rate (n = 3,552) in parasite infections from Western Thailand (2001–2014). We describe 32 independent coding mutations including common mutations outside the kelch13 propeller associated with significant reductions in clearance rate. Mutations were first observed in 2003 and rose to 90% by 2014, consistent with a selection coefficient of ∼0.079. ART-R allele diversity rose until 2012 and then dropped as one allele (C580Y) spread to high frequency. The frequency with which adaptive alleles arise is determined by the rate of mutation and the population size. Two factors drive this soft sweep: (1) multiple kelch13 amino-acid mutations confer resistance providing a large mutational target—we estimate the target is 87–163 bp. (2) The population mutation parameter (Θ = 2Neμ) can be estimated from the frequency distribution of ART-R alleles and is ∼5.69, suggesting that short term effective population size is 88 thousand to 1.2 million. This is 52–705 times greater than Ne estimated from fluctuation in allele frequencies, suggesting that we have previously underestimated the capacity for adaptive evolution in Plasmodium. Our central conclusions are that retrospective studies may underestimate the complexity of selective events and the Ne relevant for adaptation for malaria is considerably higher than previously estimated.

Fitness costs are key determinants of whether drug resistance alleles establish and how fast they spread within populations. More than 125 different kelch13 alleles, each containing a different amino acid substitution, have arisen in Southeast Asian malaria parasite (Plasmodium falciparum) populations under artemisinin selection over the past 15 years in a dramatic example of a soft selective event.

Determining the genetic basis of fitness is central to understanding evolution and transmission of microbial pathogens. In human malaria parasites (Plasmodium falciparum), most experimental work on fitness has focused on asexual blood stage parasites, because this stage can be easily cultured, although the transmission of malaria requires both female Anopheles mosquitoes and vertebrate hosts. We explore a powerful approach to identify the genetic determinants of parasite fitness across both invertebrate and vertebrate life-cycle stages of P. falciparum. This combines experimental genetic crosses using humanized mice, with selective whole genome amplification and pooled sequencing to determine genome-wide allele frequencies and identify genomic regions under selection across multiple lifecycle stages. We applied this approach to genetic crosses between artemisinin resistant (ART-R, kelch13-C580Y) and ART-sensitive (ART-S, kelch13-WT) parasites, recently isolated from Southeast Asian patients. Two striking results emerge: we observed (i) a strong genome-wide skew (>80%) towards alleles from the ART-R parent in the mosquito stage, that dropped to ~50% in the blood stage as selfed ART-R parasites were selected against; and (ii) repeatable allele specific skews in blood stage parasites with particularly strong selection (selection coefficient (s) ≤ 0.18/asexual cycle) against alleles from the ART-R parent at loci on chromosome 12 containing MRP2 and chromosome 14 containing ARPS10. This approach robustly identifies selected loci and has strong potential for identifying parasite genes that interact with the mosquito vector or compensatory loci involved in drug resistance.