Skip to main content


Association of mutations in the Plasmodium falciparum Kelch13 gene (Pf3D7_1343700) with parasite clearance rates after artemisinin-based treatmentsa WWARN individual patient data meta-analysis



Plasmodium falciparum infections with slow parasite clearance following artemisinin-based therapies are widespread in the Greater Mekong Subregion. A molecular marker of the slow clearance phenotype has been identified: single genetic changes within the propeller region of the Kelch13 protein (pfk13; Pf3D7_1343700). Global searches have identified almost 200 different non-synonymous mutant pfk13 genotypes. Most mutations occur at low prevalence and have uncertain functional significance. To characterize the impact of different pfk13 mutations on parasite clearance, we conducted an individual patient data meta-analysis of the associations between parasite clearance half-life (PC1/2) and pfk13 genotype based on a large set of individual patient records from Asia and Africa.


A systematic literature review following the PRISMA protocol was conducted to identify studies published between 2000 and 2017 which included frequent parasite counts and pfk13 genotyping. Four databases (Ovid Medline, PubMed, Ovid Embase, and Web of Science Core Collection) were searched. Eighteen studies (15 from Asia, 2 from Africa, and one multicenter study with sites on both continents) met inclusion criteria and were shared. Associations between the log transformed PC1/2 values and pfk13 genotype were assessed using multivariable regression models with random effects for study site.


Both the pfk13 genotypes and the PC1/2 were available from 3250 (95%) patients (n = 3012 from Asia (93%), n = 238 from Africa (7%)). Among Asian isolates, all pfk13 propeller region mutant alleles observed in five or more specific isolates were associated with a 1.5- to 2.7-fold longer geometric mean PC1/2 compared to the PC1/2 of wild type isolates (all p ≤ 0.002). In addition, mutant allele E252Q located in the P. falciparum region of pfk13 was associated with 1.5-fold (95%CI 1.4–1.6) longer PC1/2. None of the isolates from four countries in Africa showed a significant difference between the PC1/2 of parasites with or without pfk13 propeller region mutations.

Previously, the association of six pfk13 propeller mutant alleles with delayed parasite clearance had been confirmed. This analysis demonstrates that 15 additional pfk13 alleles are associated strongly with the slow-clearing phenotype in Southeast Asia.


Pooled analysis associated 20 pfk13 propeller region mutant alleles with the slow clearance phenotype, including 15 mutations not confirmed previously.


Artemisinin-based combination therapies (ACTs) have been the WHO recommended treatment for uncomplicated Plasmodium falciparum malaria since 2006. Studies conducted in 2006–2007 reported that P. falciparum parasites in northwest Cambodia had reduced in vivo susceptibility to artemisinins [1, 2]. This was manifest as delayed clearance of parasites from the blood of patients treated with artesunate monotherapy or ACTs [2]. In 2014, parasites selected in vitro for delayed response to artemisinins were found to have a  mutation in the kelch13 gene (pfk13, P. falciparum 3D7_1343700) [3, 4], a locus consistent with genome association studies that had strongly associated a region of chromosome 13 with slow parasite clearance [5, 6]. Subsequent clinical studies in western Cambodia showed that delayed parasite clearance in clinical trials was associated with several non-synonymous mutations in pfk13. In particular, mutations in the distinctive propeller region of the Kelch 13 protein (codons 441–726) were associated with slow parasite clearance and subsequently with reduced artemisinin susceptibility in in vitro studies assessing susceptibility of the ring-stage parasites [3, 7, 8]. Recently, the structure of the propeller region was solved and made public by the Structural Genetics Consortium [9], making it possible to link particular genetic changes with their position in the molecule. A central role of pfk13 propeller mutations in mediating ring-stage resistance was confirmed by demonstration that parasites engineered to contain the mutations showed ring-stage resistance to artemisinin whereas parasites with the wild type allele were sensitive [10].

As pfk13 propeller sequences were determined in P. falciparum from the Greater Mekong Subregion (GMS), other Asian sites [6, 11,12,13,14,15,16,17,18,19,20,21,22], and Africa, Latin America, and Oceania [23,24,25,26]; almost 200 different non-synonymous mutations have been identified (see [27] and supplementary tables 6 and 7 in [28]).

The parasite clearance half-life (PC1/2), which measures the slope of the log-linear component of the parasite clearance curve, has become established as the best in vivo metric of P. falciparum artemisinin susceptibility [29,30,31]. Although sporadic parasite isolates with pfk13 propeller mutations and occasional patients with slow parasite clearance have been identified in many parts of the malaria-endemic world, pfk13 mutant parasites isolated from patients with the slow-clearing in vivo phenotype have been demonstrated only in a circumscribed region of Southeast Asia, initially in western Cambodia and currently extending to parts of Thailand, Vietnam, Lao PDR, Myanmar, and Yunnan Province, China. Within this region, there has been clear evidence of selection and spread of successful artemisinin-resistant parasite lineages [32]. As expected, artemisinin resistance has led to further selection of ACT partner drug resistance [33,34,35]. Outside the GMS, extensive molecular analyses have shown no evidence that parasites carrying pfk13 propeller SNPs are under directional selection by artemisinins [23, 25, 27, 28, 36]. Taken together, available studies generally support the conclusion that certain pfk13 propeller mutations mediate delayed parasite clearance following treatment with artemisinin derivatives in Southeast Asia, but that the pfk13 propeller region mutants observed outside the GMS do not show an association with slow parasite clearance [37, 38]. Exceptions include recent molecular reports of independent emergence of parasites that carry pfk13 propeller mutant alleles observed previously in the GMS in Guyana (C580Y) [39] and Rwanda (P574L and A675V) [40]. Accurate measurement of the parasite clearance phenotype requires frequent timed blood sampling [41]. This is often impractical and means that many surveys only collect information on the pfk13 genotypes in a study area. It is therefore critical to utilize optimally those studies that do include both genotype and in vivo phenotype. The WHO has provided guidelines for associating P. falciparum parasite mutations with slow parasite clearance in vivo [42]. Currently, 6 pfk13 propeller region mutant alleles have been validated in this way, 8 additional alleles have been described as “associated” with the phenotype, and 18 other alleles have been isolated in such low numbers that their phenotype could not be evaluated. It is also clear that some non-synonymous mutations in the pfk13 propeller region (notably A578S) are not associated with artemisinin resistance. The low prevalence of pfk13 propeller alleles outside the Greater Mekong Subregion has limited our ability to confirm the roles of additional genotypes in conferring the delayed clearance phenotype.

We have used meta-analysis of individual patient data from published and unpublished studies to compare standardized parasite clearance half-life estimates and pfk13 propeller region genotypes. This analysis capitalized on the power of the large data set to allow assessment of the importance of a broader range of propeller region alleles in mediating slow parasite clearance in vivo.


Data acquisition

A systematic literature review following the PRISMA protocol [43] was conducted to identify studies published between 2000 and 2017 which included frequent parasite measurements and pfk13 genotyping (last run on 31 January 2018). Four databases (Ovid Medline, PubMed, Ovid Embase, and Web of Science Core Collection) were searched. The search terms and conditions are available in Additional file 1.

At the abstract screening stage, we excluded studies that did not include treatment with artemisinin derivatives, pfk13 genotyping or parasite measurements as well as review and correspondence articles. At full-text screening, we excluded additional studies that measured parasites less frequently than twice a day. Data were actively requested only from studies that reported parasites that carried pfk13 non-synonymous mutations. However, we received two studies that found no parasites that carried these pfk13 mutations and these were analyzed as part of the wild type data set.

Studies that identified non-synonymous mutations in pfk13 were eligible for inclusion in this analysis if individual patient files also met the following criteria: (a) patients were treated with either an ACT or an artemisinin monotherapy; (b) parasitemia was measured in the first days of treatment at least every 12 h, until a negative count or at least until day 3, allowing the standardized calculation of the PC1/2 [29, 30]; (c) the dosing protocol was available; and (d) the weight of each individual patient was recorded. The published papers and data sources are listed in detail in Additional file 2: Table S1 [2, 11, 15, 16, 19, 44,45,46,47,48,49,50,51].

Individual study protocols were available for all trials included, either from the publication or as a metafile submitted with the primary data. Individual patient data from eligible studies were shared, collated, and standardized under the protocols of the WWARN data platform using previously described methodology [52]. Study reports generated from the formatted datasets were sent back to investigators for validation or clarification. Methodologies used in studies regarding parasitemia sampling and molecular analyses are presented in Additional file 2: Table S1.

Statistical analysis

The statistical analysis plan [53] was developed before analysis. All non-synonymous mutations in the pfk13 gene identified in the studies were included in the analysis. Isolates without reported mutations were assumed to be wild type in assessing relationships between parasite genotype and PC1/2. Isolates with a mixed genotype at any nucleotide within the pfk13 coding region (wild type/mutation or two non-synonymous mutations) were excluded from the analysis.

The PC1/2 is defined as the time in hours needed for the parasite density to decline by 50% during the log-linear phase of parasite clearance. PC1/2 was calculated using the WWARN parasite clearance estimator tool [41]. The goodness of fit of parasite clearance models was evaluated for each individual patient parasitemia-time profile used to estimate the PC1/2.

Profiles that satisfied the following criteria (i.e., provided biologically or statistically plausible results) were included in the analysis: (a) standard deviations of residuals < 2, (b) number of data points used to fit the linear part of the curve > 2, (c) duration of lag phase < 12 h, (d) pseudo R2 statistics ≥ 0.8. Additionally, patients who withdrew or had a record of inadequate dosing were excluded. The log transformed half-life metric was modeled for all pfk13 mutant alleles in all studies with information from individual patients on age, initial parasitemia, ACT treatment, and artesunate dose as covariates. The method by which the dose was calculated is documented in the statistical analysis plan. Random effects for study site were used to account for heterogeneity between studies. Residuals were examined for normality and for systematic deviations from the model.

The differences in PC1/2 between infections with P. falciparum parasites bearing a specific pfk13 propeller mutant allele and those with wild type parasites were assessed by the Wald test. The fold change in geometric mean of PC1/2 of infections with pfk13 mutant parasites compared to wild type isolates from the same sites; xPC1/2 was calculated as an exponent of the difference of the corresponding regression coefficients.

In order to determine a PC1/2 threshold value that defined slow parasite clearance, we divided Asian data into two populations: rapid clearing and slow clearing. The slow-clearing population was defined as all isolates with mutations associated with a significant increase in PC1/2 values in this analysis, while the fast-clearing population included all other isolates. The PC1/2 value that corresponds to the 95th percentile of the fast-clearing population (i.e., a value x such that the probability that PC1/2 > x is less than 0.05) was selected as the cutoff for infections with “slow clearing” parasites. Risk of bias in individual studies was assessed based on frequency of parasite counting, molecular methodology, and number of patients excluded because of missing data or unsatisfactory fit of the model for PC1/2 estimation (for details, see Additional file 1). Data from studies/sites that reported results very different from all of the others in the same region were included in the analysis, and a sensitivity analysis was conducted after excluding Tra Lang, Vietnam (study site ID 23; study ID 8), and Pyin Oo Lwin, Myanmar (study site ID 15; study ID 13).


Literature search

The literature search identified 146 articles, with 9 satisfying our inclusion criteria and describing data from 14 studies (Additional file 1). Seven additional studies were contributed in response to the proposal for this study group on the WWARN website. Consequently, we requested individual patient data from 21 studies and received data from 18 studies.

Description of the study sites

There were 16 studies from Asian sites that reported both PC1/2 estimates and pfk13 propeller region genotypes from Myanmar, Thailand, Cambodia, Lao Peoples Democratic Republic (PDR), Vietnam, Bangladesh, and Yunnan province of China bordering Myanmar. Among these 3179 patients, 3012 had sufficient data for PC1/2 estimation and 2631 met the inclusion criteria for evaluation of the association of the pfk13 genotype and the parasite clearance phenotype. Of the Asian isolates, nearly half (43%; 1142/2631) were from four study sites of the Shoklo Malaria Research Unit (SMRU) along the Western Thailand/Eastern Myanmar Border; these isolates are identified as a group called Thai Western Border.

There were 238 isolates with genotype and phenotype information contributed from five sites in Africa; one site in Nigeria (n = 31), the Democratic Republic of Congo—DRC (n = 119), and Tanzania (n = 39); and two sites in Madagascar (n = 49). A total of 204 patients’ data met the inclusion criteria for evaluation of the association of the pfk13 genotype and the parasite clearance phenotype (Table 1, Figs. 1 and 2).

Table 1 Summary of genotyping results by study site. Samples with mixed genotyping are excluded. Details of mutations are listed in Additional file 2: Table S2A (codons 1–440) and S2B (codons 441–726), and mixed genotypes are listed in Additional file 2: Table S3. Results are sorted by continent, country, site, and year of study
Fig. 1

Study site locations. The numbers in the pins correspond to the study sites listed in Table 1

Fig. 2

Study profile. Definitions for specific exclusions are listed at the right of the figure detailing the number of isolates included in each analysis. Unsatisfactory fit was defined as pseudo-R2 statistics < 0.8. Insufficient parasite data includes patients with too few observations to fit the model and patients with only daily counts

Figure 1 depicts the locations of the study sites. Additional file 3: Figure S2 summarizes the numbers of individual data sets that met all criteria for inclusion.

Patient characteristics and the numbers of patients treated with artesunate or artemisinin combination therapies are summarized in Table 2. The median age of patients in Asia was 21 years, and 77% were 12 years of age or older; African patients had a median age of 4.7 years, and only 6% were older than 11 years.

Table 2 Baseline characteristics and treatment administered. Only patients with genotyping and PC1/2 results are presented

Molecular analyses of the pfk13 alleles

Samples were collected between 2001 and 2014. Among the 45 individual study sites, 33 included data from the entire pfk13 gene (codons 1–726); 12 studies reported sequence of only the propeller region codons (441–726). A total of 202 and 1254 isolates with molecular information on mutants in codons 1–440 (P. falciparum-specific region) and codons 441–726 (the propeller region) were recorded respectively. The pfk13 genotype of parasites from 3424 patients was determined; mutant alleles were identified in 1455 samples. The prevalence of mutant codons observed in each study is summarized in Table 1.

Plasmodium-specific region (codons 1–440)

Data were available for the whole pfk13 molecule in two African study sites: Kinshasa, DRC (study ID 13; site ID 24), and Ilorin, Nigeria (study ID 13; site ID 36). Ilorin had the higher prevalence of isolates with mutant codons in the P. falciparum domain (14/36, 39%), compared with Kinshasa (33/117, 28%). In 76% (37/49) of these isolates, the mutations were at codon 189; 35 carried the K189T mutation, 12 carried a synonymous K189K, and 2 carried a K189N codon (Additional file 2: Table S2A).

Complete pfk13 DNA sequence was determined in 18 Asian study sites, including all of the studies from the Thai Western Border. Mutant codons in the 5′ region of the gene were rare in most Asian-derived samples. The K189T mutant codon was observed in samples from Bangladesh (study 7; site 1: n = 1/21, 5%), (study 16; site 2: n = 4/55, 7%) and two sites in Myanmar ((study 16; site 16: n = 3/74, 4%), study 12; site 17: n = 1/71, (1.4%)) (Additional file 2: Table S2A).

Many of the parasites from the Thai Western Border (site 20) shared an allele with a codon change at position 252 from glutamic acid to glutamine (E252Q). Parasites with this allele were observed in study 1 (n = 1/7, 14%), study 4 (n = 68/950, 7%), and study 16 (n = 14/116, 12%). Parasites with this genotype were also observed in Shwe Kyin in 2011–2013 (study 16, site 16), where they comprised 10/74 (14%) of the isolates tested (Additional file 2: Table S2A), and this allele has been observed near the Myanmar-China border, as well [14]. In the relatively conserved stem domain (350–440), mutants were observed only transiently in parasites from Shwe Kyin and the Thai Western Border.

pfk13 propeller region (codons 441–726)

The sequence of the pfk13 propeller region was determined in all 45 data sets (Table 1, Additional file 2: Table S2B). In the African isolates, the only mutant codons identified in the propeller region were in 3/117 isolates from Kinshasa, each with a single mutation (codons S522C, A578S, or Q613L) and in 4/44 isolates from Brickaville, Madagascar, all with mutant codon A578S which is relatively common in Africa and has been reported in Thailand and Bangladesh [27]. The Asian data set was dominated by large numbers of patients from two regions, the Thai Western Border (n = 1291/3087 (42%)) and northwestern Cambodia (n = 847/3087 (27%)), with smaller numbers from Vietnam, Myanmar, the Myanmar-China Border, Bangladesh, and Lao PDR. The detailed data are summarized in Additional file 2: Table S2B. As has been published, the Thai Western Border, Thai, Cambodian, Laotian, and Vietnamese data sets contained numerous isolates with the pfk13 C580Y allele [11,12,13,14, 18, 19, 21, 22, 28, 32, 34, 50, 54, 55]. In contrast, isolates with the Y493H and R539T mutant codons were common among Cambodian isolates, but uncommon in samples from Thailand with the exception of those from Srisaket, a province that borders Cambodia.

The F446I mutant predominated in data sets from Myanmar and the China/Myanmar border (n = 83/141, (59%) [14, 16, 27, 51]. The mutant codons in the propeller region identified among the African isolates, S522C, S613E, and the very common A578S “African” allele, were absent from the Asian data set.

Association of pfk13 mutations and parasite clearance half-life

Among 3329 patients with pfk13 genotyping results that reported only a single mutation in the propeller region, PC1/2 values were estimated in 3156 (95%) patients. In 165 patients, PC1/2could not be estimated because the data were too sparse, the initial parasitemia was too low (< 100 parasites per microliter) (n = 3) or the last recorded parasitemia was too high (> 10,000 parasites per microliter) (n = 6). Overall, 90% (2834 of 3156) of PC1/2 estimates were classified as suitable for analysis (using criteria specified prospectively in the methods, and Fig. 2).

Defining the groups of slow- and fast-clearing parasites

Two populations of parasites were identified among the Asian parasites: fast-clearing (loge PC1/2 mean, 1.0 SD, 0.41) and slow-clearing including all isolates with mutations listed in Tables 3 and 7 (logePC1/2 mean 1.9; SD 0.19). These groups had corresponding half-life geometric means of 2.7 h and 6.7 h, respectively. A PC1/2 equal to 5.5 h corresponds to the 95th percentile of the fast-clearing parasite population, and the probability of observing PC1/2 > 5.5 h in a single isolate coming from a sensitive population is equal to 0.043. On this basis, a median value of PC1/2 > 5.5 h was used to define parasites associated with slow clearance.

Table 3 Comparison of PC1/2 between patients with WT and a pfk13 NS-mutation in Asia. Metric is the fold increase of PC1/2 in mutants vs. wild type. Only mutations with n ≥ 5 are shown. With only two exceptions (K189T and K438N), all patients that carried parasites with the indicated pfk13 mutant codon showed a PC1/2 significantly different from those that carried WT parasites

African isolates

Summary statistics of PC1/2 estimates and their association with pfk13 mutations are presented in Table 5 and Additional file 3: Figure S3A for wild type (median PC1/2 2.2 h, range 0.7–6.3) and mutant K189T (median PC1/2 2.1 h, range 0.8–7.1). Additional file 3: Figure S3B shows individual PC1/2 values for all other mutations observed in African isolates.

The PC1/2 values in wild type isolates from the five African sites were all well below the cutoff of 5.5 h, but varied significantly from one another (p = 0.004, Additional file 3: Figure S1). Compared to wild type isolates in Kinshasa, DRC, an area of high perennial transmission, the PC1/2 in wild type isolates was 1.36 times longer (95%CI 1.08–1.70; p = 0.008) in Brickaville, Madagascar, 1.46 times longer (95%CI 1.18–1.8; p = 0.001) in Ilorin, Nigeria; and not significantly different in Fukayosi, Tanzania (p = 0.117), or Ankazobe, Madagascar (p = 0.199) (Additional file 3: Figure S1).

Most pfk13 mutations were rare in the African data sets, being observed in fewer than four isolates (Figure 3). Among a total of seven isolates with a propeller mutation (S522C, A578S, or Q613L), four met the inclusion criteria for evaluation of the association of the pfk13 genotype and the parasite clearance phenotype. All of the seven isolates had PC1/2 < 2.8 h and exhibited no increase in PC1/2 compared to infections with wild type parasites. In fact, the geometric mean of PC1/2 was lower in isolates with propeller mutations compared to wild type isolates at the same study sites (p = 0.491 and p = 0.096 when seven or four isolates were included). In this dataset, only the A578S allele of the propeller mutations observed in Africa was also observed in Asia, but in only three isolates. The median half-life for mutations anywhere in the pfk13 gene ranged from 1.7 to 2.7 h (Additional file 3: Figure S3A and S3B).

Fig. 3

Distribution of parasite clearance half-life by individual wild type or pfk13 mutant codons in isolates from Africa. Left panel (a) shows violin plots for wild type or mutant codons with five or more isolates. The number of individual isolates tested is at the right of each violin plot. The right panel (b) shows dot plots for mutant codons with < 5 isolates. The red line shows a half-life of 5.5 h. The median is shown as a green circle, the red bar corresponds to the interquartile range, and the curve represents kernel estimate of the density function

Codon K189T/N was present in 33 isolates (31 were K189T and 2 were K189N), sufficient numbers to assess the overall parasite clearance; these infections had a median PC1/2 of 2.1 h (range 0.8–7.1) similar to that of infections with wild type parasites; 2.2 h (range 0.7–6.3) (Fig. 3a).

Asian isolates

Among Asian isolates, 1440 were wild type and 1190 carried a single pfk13 mutation. The prevalence of isolates with a PC1/2 > 5.5 h and the median PC1/2 both varied depending on the pfk13 mutation (Table 5, Fig. 4).

Fig. 4

Distribution of parasite clearance half-life by pfk13 mutant codons in isolates from Asia. Left panel (a) shows violin plots for mutant codons with five or more isolates. The number of individual isolates tested is at the right of each violin plot. The right panel (b) shows dot plots for mutant codons with < 5 isolates. The red line shows a half-life of 5.5 h. The median is shown as a green circle, the red bar corresponds to the interquartile range, and the curve represents kernel estimate of the density function. * 22 P527H and 1 P527L with PC1/2 = 5.8h; ** 6 P667T and 1 P667R without valid PC1/2 measurement, so not considered

The 32 studies with infections caused by P. falciparum isolates that carried a non-synonymous change within codons 1–440 were analyzed, and all 45 studies with sequence information on the propeller regions were analyzed.

Variance in PC1/2 between genotypes was significantly larger than within genotypes (F40, 2554 = 48.43, p < 0.001, after adjusting for study site). All infections with pfk13 mutants in the propeller region that were reported in five or more isolates were associated significantly with a geometric mean PC1/2 greater than that in infections with wild type parasites (p < 0.001 for all but codon A481V, in which p = 0.002). The fold increase in geometric mean of PC1/2 (xPC1/2), compared to the PC1/2 of wild type parasites at the same sites, ranged from 1.5 to 2.7 (Table 3). Parasites with mutations F446I, G449A, A481V, P527L, N537I, G538V, and P574L had a fold increase between 1.5 and 1.9, and all others were between 2.0 and 2.7. Neither isolates with the K189T mutation nor K438N showed any increase in PC1/2 compared with wild type parasites (p = 0.644 and 0.237, respectively). However, isolates with the E252Q mutation (which is also outside the propeller region) from both the Western Thai Border and Shwe Kyin, Myanmar, had a significant increase in PC1/2, 1.5-fold (xPC1/2 = 1.5; 95%CI 1.4–1.6) compared with wild type parasites (p < 0.001).

For infections caused by parasites that carried only wild type pfk13, the distribution of PC1/2 values varied significantly between locations (Additional file 3: Figure S2, p < 0.001); one site in Myanmar, Pyin Oo Lwin, and one in Vietnam, Tra Lang, had a median PC1/2 of wild type isolates that was above the cutoff of 5.5 h and above the PC1/2 observed for wild types in other sites (p < 0.001). At both sites, no significant differences were observed in PC1/2 between K13 wild type and K13 mutant isolates. The corresponding factor for change in PC1/2 in mutant isolates was estimated as 0.91 (0.74–1.12) p = 0.367 for Tra Lang, Vietnam, and 0.90 (0.73–1.10) p = 0.274 for Myanmar, Pyin Oo Lwin. Four mutant pfk13 alleles were associated with significant differences in PC1/2 between different study sites: E252Q (p < 0.001), F446I (p = 0.005), M476I (p = 0.014), Y493H (p = 0.029) and C580Y (p < 0.001) (Table 4, Additional file 3: Figures S3-S7).

Table 4 Comparison of PC1/2 for pfk13 NS mutations between study sites in Asia. Only mutations with isolates available from at least five patients from at least two sites are shown

Between 2009 and 2014, the prevalence of the C580Y allele in parasite populations on the Thailand western border increased [18]. This change was accompanied by a progressive increase in the median PC1/2 of infections with pfk13 C580Y mutants, from a median of 5.4 h in 2009 to 7.2 h in 2014 (Additional file 3: Figure S8). A linear trend of increase in the geometric mean of PC1/2 by 5.0% (95%CI 2.3–7.9) each year (xPC1/2 = 1.05; 95% CI 1.02–1.08) was observed (p < 0.001), and the linear trend was not affected by differences in patient treatment or artemisinin derivative dose.

Influence of patient characteristics and antimalarial treatment drugs on PC1/2 in Asian isolates

Differences between PC1/2 values for parasites with the same pfpk13 genotype were not significant when adjusted for patients’ initial parasitemia, artesunate dose, age, and study site. However, within the data set, patients whose parasite clearance was assessed had been treated with a range of artemisinin derivatives and doses and ACT partner drugs (Table 2). Groups with sufficient numbers of isolates (WT, F446I, Y493H, R539T, I543T, P553L, P574L, C580Y) were compared to assess the effects of the drugs administered on the half-life of parasite clearance. No differences in parasite clearance half-life among these different drug treatments were observed with one exception: the comparison between treatment with artesunate alone and artesunate/mefloquine in isolates with a WT genotype. Infections in patients treated initially with artesunate + mefloquine had a geometric mean PC1/2 0.8 times (95%CI 0.7–0.9, p < 0.001) shorter than those treated with artesunate alone.

Correspondence of PC1/2 > 5.5 h and day 3 parasitemia

Collecting frequent parasite counts in patients several times per day is not always feasible, and so recording the proportion of patients that remain parasitemic at day 3 after treatment has been proposed as a simple more practical warning signal for the possible slow clearance phenotype [29, 56]. To assess the correspondence between these two approaches, the proportion of patients with day 3 parasitemia in these data sets was compared with the proportion of isolates with a PC1/2 > 5.5 h (Table 5, Additional file 3: Figure S9). Use of the day 3 positivity test yielded an overall 11% (226/1975) false positive (FP) rate and 14% (107/772) false negative (FN) rate when compared to the use of the cutoff of PC1/2 > 5.5 h.

Table 5 Summary of PC1/2 by mutant codon and region showing proportion of isolates from a study site with PC1/2 > 5.5 h and proportion of Day 3 positive isolates. Mutant codons represented by fewer than 5 isolates are indicated in bold

However, since the initial parasitemia is a major determinant of time to clear microscopy detectable parasites after treatment, the proportion of FP and FN varies across a range of initial parasitemia values with 51% FN for parasite numbers less than 10,000/μL and 29% FP for parasite numbers greater than 300,000/μL (Table 6).

Table 6 Utility of day 3 parasite positivity for detection of isolates with long PC1/2 (> 5.5 h). False positive and false negative outcomes are estimated over a range of initial parasitemia values

Association of additional pfk13 mutant alleles with slow parasite clearance

The current pooled analysis has allowed definition of additional pfk13 genotypes that are associated strongly with slow parasite clearance (Tables 3 and 7). Six pfk13 genotypes have already been associated with slow clearance by WHO, and 11 alleles were listed as candidates [57]. These have all now been demonstrated to be strongly associated with slow clearance (p < 0.001 except allele A481V, p = 0.002). In addition, 6 alleles that previously had been observed infrequently were present in high enough numbers in the analysis data set to allow assessment of the associated PC1/2; these alleles also showed a strong association with prolonged parasite clearance. All of the associated mutant alleles had single changes in the pfk13 propeller region, with one exception, allele E252Q, located in the N-terminal Plasmodium-specific region of the protein. Two alleles with a smaller effect on parasite clearance, E252Q and F446I, could now be assessed with confidence since a larger sample size was available (n = 703 and n = 382, respectively). Parasites with the pfk13 E252Q mutation had a median PC1/2 equal to 4.4 h and 11% had a PC1/2 > 5.5 h, those that carried a F446I mutation had a median PC1/2 of 5.4 h and a PC1/2 > 5.5 h was recorded in 47% of the patients.

Table 7 List of new pfk13 mutant alleles strongly associated with slow parasite clearance

Risk of bias

The risk of bias in individual studies was considered low for majority of studies (Additional file 1).

Sensitivity analysis performed after exclusion of the only study (study ID 8, site ID 23) with moderate risk of bias due to genotyping methods shows the same results as analysis on the full dataset (Table 3). Only the mutant codon I543T was reported at this study/site. Additionally, one study site (study ID 13, site ID 15) presented significantly different PC1/2 results from other sites for parasites with no pfk13 mutant codons (Additional file 3: Figure S1). The mutant codon P574L predominated in that site. We tested the sensitivity of the association of codon P574L with slow clearance by including and excluding the site ID 15 data (Table 3 and Table 7). The association of the P574L codon with slow parasite clearance was highly significant in both analyses (p ≤ 0.001).

Data from three apparently relevant published studies were not shared [3, 22, 38]. However, the association of particular mutations with slow clearing parasites coincides with our analysis. In Southeast Asia, alleles Y493H, 539T, I543T, P553L, V568G, and C580Y in Thuy-Nhien et al. [22] and C580Y, Y493H, and R539T in Ariey et al. [3] were associated with slow clearance, and in Mali, the 9 pfk13 mutant parasites identified by Ouattara et al. were not associated with slow clearance [38].


This is the largest evaluation of the association of different pfk13 mutations with slow parasite clearance following antimalarial treatment with the artemisinin derivatives. Our study confirms that many, but not all (notably A578S), mutations in the propeller region of the pfk13 gene are associated with slow parasite clearance. It also shows that a mutation outside the propeller region, E252Q which had been proposed to confer artemisinin resistance, is strongly associated with slow parasite clearance. This is also consistent with the increased prevalence for several years on the northwest Thailand-Myanmar border of parasites that carried the E252Q mutation, presumably as a result of selection, before the ascendance of the C580Y genotype conferring a more extreme phenotype [18].

Several factors affect the rates of parasite clearance following administration of artemisinin derivatives and would have contributed to the inter-individual and between-site differences observed in this study. In infections with pfk13 wild type parasites, treatment with an artemisinin derivative results in accelerated ring-form clearance, and this is reflected in a steep slope of the parasitemia-time profile and a derived PC1/2 which is usually well below 5.5 h. Background immunity, or its surrogate, age has a significant additional effect, further accelerating parasite clearance. Thus, immunity as reflected by antibody concentrations can have a relatively small but significant effect [58]. Another important contributor to the half-life is the stage of parasite development at the initial presentation of the patient. This variable, too, could disproportionately affect clearance of pfk13 mutant parasites compared to wild type infections [59] but should only affect inter-individual and not inter-site differences. In artemisinin combination treatments, the partner drug also makes a small contribution to the rate of parasite clearance, which in this study, was significant only for the mefloquine combination, confirming previous findings [60]. It is important to note that despite modest site-specific differences, the clearance half-life derived for each pf13 propeller alleles in the pooled data set was at least 1.5-fold higher than the half-life of the wild type parasites, and for 14 of the 20 propeller mutant alleles, the ratio of mutant to wild type PC1/2 was between 2.0 and 2.7.

The prominent role of pfk13 propeller mutants in the selective response to artemisinin drug pressure has been recently demonstrated by the rapid increase in parasites that carry a genotype with a particular version of the C580Y allele. These parasites originated in northwestern Cambodian and Thai foci and spread recently to new areas of Cambodia and Vietnam [22, 32]. On a more local scale, in this study, we found that the median PC1/2 of parasites that already carried a C580Y allele of pfk13 increased from 5.4 h in 2009 to 7.2 h in 2014 in the western border region of Thailand. More widely, among the 14 pfk13 propeller mutants newly associated with slow clearance, 4 of the 14 mutant alleles associated with slow clearance were present in at least 8% of the isolates identified in the corresponding parasite population (range 8–31%). Ten of the newly identified pfk13 mutant alleles were present at low prevalence (range 1–6%) in the sites where they have been reported and even among the validated alleles, prevalence was below 10% for four of the six alleles. Because the pooled data set contains isolates collected from many different sites and over a time span of 14 years, no trends in overall prevalence can be inferred in this analysis. However, the value of following temporal changes in the prevalence of propeller mutants in local sites is clear, as has been demonstrated for parasites that carry the C580Y allele in Thailand, Cambodia, and Vietnam [22, 32].

Detailed temporal studies pfk13 mutant allele prevalence both in the GMS and in other worldwide sites represent a powerful tool for assessing whether a particular pfk13 mutant population is increasing, a signal that further studies of the parasite phenotype may be warranted in that location. Temporal studies will be particularly important for the seven newly identified propeller mutant alleles in our data set that were observed at prevalences of 4% or lower.

In contrast, the pfk13 mutant parasites in all five African sites remained at very low prevalence, generally below 3%, and no evidence of slow-clearing parasites or selection for mutant parasites has been identified. The low prevalence of pfk13 mutant parasites in the African studies in our data set depends on relatively small numbers of parasites (e.g., 31, 26, 29). However, the generality of this observation is supported by extensive studies in many sites in Africa including some assessments of parasite clearance in vivo [37, 38, 61,62,63], protection against parasite exposure to artemisinins of cultured parasites in vitro [10, 62, 64], and widespread molecular surveillance of pfk13 propeller region sequence ([23,24,25,26, 55, 65,66,67,68] and see the WWARN K13 surveyor and WHO drug resistance threats map for details). These many reports also support the conclusions of our study: in African sites, pfk13 propeller mutants are diverse but rare, with no evidence for selection even where artemisinin-based antimalarials have been intensively used. Moreover, the pfk13 propeller mutant alleles commonly observed in Asia have rarely been observed in Africa [40, 69]. In addition, only two African sites report a somewhat higher prevalence of parasites with pfk13 propeller mutations, 14% (18/130) in Mali [23] and an increase from 3% (1/31) in 2007 to 27% (7/26) in 2014 in Grande Comoros [70] and 17% (5/29) mutants in individuals from the neighboring island of Mayotte [71].

The pfk13 gene has been identified as essential to P. falciparum in its blood stream stage [72], and it is likely that parasites dependent upon a mutant K13 protein would incur a considerable fitness disadvantage [64]. In this situation, parasite populations subjected to intensive use of artemisinin-based antimalarials would be expected to be under intense selection for pfk13 mutant parasites that carry additional genetic changes that could compensate for the fitness cost as has been observed in other parasite populations under drug selection [73]. There is increasing evidence that such adaptive responses have evolved in Southeast Asian P. falciparum populations [74]. For example, intensive use of the ACT dihydroartemisinin-piperaquine in Cambodia and Thailand has apparently exerted further selection on the parasites that carry the pfk13 C580Y “spreading” genotype, selecting parasites that also carry an increased copy number of plasmepsin 2–3 that confers resistance to the piperaquine partner drug [33,34,35]. In addition to selection and spread of pfk13 mutant parasites, independent populations that carry new versions of the common C580Y or one of many novel pfk13 alleles have emerged in many other sites within the GMS region [12,13,14,15, 21, 22, 28, 54].

Even before the pfk13 gene was identified as a potential molecular marker of reduced artemisinin response, genome-wide association studies identified several regions of the P. falciparum genome associated with slow parasite clearance on chromosomes 10 and 14, in addition to chromosome 13 where the pfk13 gene is located [5, 6, 75]. Whole genome analyses have also identified genetic changes associated both with pfk13 propeller mutants and with slow clearance in independent, low diversity populations in Cambodia [76, 77] and in Western Thailand [78]. Parasites that carry these additional genetic changes are present in many populations in the western regions of Southeast Asia [21, 28, 79]. Moreover, the importance of other genetic changes for the expression of the slow clearance phenotype has been observed in vitro. Introduction of a pfk13 mutant allele into a cloned, fast-clearing wild type parasite of recent Cambodian origin conferred more protection against dihydroartemisinin in vitro than when the same allele was expressed in a wild type parasite line of African origin [10]. These observations support the interpretation that the independent emergence of parasite populations with new pfk13 alleles is facilitated in Southeast Asia by the presence of local parasite populations with genomes already adapted to support pk13 mutant parasites when they arise. Exposure to artemisinins could then quite rapidly select emergence and selection of a local population of pfk13 mutant parasites.

These specific adaptive genomic changes that are common in Southeast Asian parasites have not been observed in African isolates, but comprehensive analyses have not yet been undertaken [28]. The slow clearance phenotype is extremely complex, dependent on changes in the overall stress response to artemisinins (see [80, 81] for comprehensive reviews) and there may be many other combinations of adaptive genomic changes that could support the fitness of pfk13 mutant parasites. In addition, it is clear that genetic changes other than pfk13 mutations can confer diminished response to artemisinins in vitro and in vivo [82,83,84,85] and surveillance plans will need to remain open to the possibility of these alternative genetic strategies.

This meta-analysis still depends on detailed multiple quantitative assessments of the parasite count, not an approach that is likely to be feasible for routine surveillance in many study sites. This study highlights the limitations in sensitivity and specificity of the standard assessment of parasite clearance on day 3 as a phenotypic signal of delayed clearance. There is no question that simpler methods for assessment of parasite responses to artemisinins will be needed in all endemic areas, with particular attention to areas where malaria transmission is low and drug pressure will be most likely to select resistance.

Our large data set allowed identification of significant association of pfk13 parasite genotype and PC1/2 phenotype even when the sample size of parasites with a particular allele was very small. For example, four different mutant alleles were observed in fewer than eight patient isolates. Even in these cases, the analysis allowed a strong genotype-phenotype association to be identified. Among these rare mutant alleles, the ratio of mutant to wild type PC1/2 ranged from 1.9 to 2.7, ten of the mutant alleles are present in more than 10% of the parasites in the tested isolates, and all but four were observed in more than one study site. This robust outcome demonstrates that analyses of the PC1/2 in patients can allow identification of mutant pfk13 alleles associated with slow-clearing parasites even when malaria transmission is low and few patients can be assessed. This suggests that some researchers in areas outside the GMS may be able to assess the PC1/2 phenotype in vivo [11, 37, 38, 61] or parasite responses to artemisinins in vitro [10, 62, 85,86,87] as a screen for artemisinin response as some groups have already done.

Finally, 18 data sets were included in this meta-analysis of individual patient data sets, including 4 trials that had not yet been published. We believe that this demonstrates that the engagement of research groups to share individual patient data can provide a unique opportunity for advancement of public health. This approach can combine and facilitate collaborative analysis of critical data sets, defining outcomes that cannot be achieved by only extracting outcomes from multiple single publications. Such meta-analyses can provide evidence for policy makers, while not precluding subsequent individual study publications.

There are limitations of this meta-analysis. Notably, the uneven representation of isolates derived from different sites of origin between Africa and Asia and within Southeast Asia and the imbalance of longitudinal data among the sites. The extensive long-term record of clinical outcomes of artemisinin treatment from patients on the Thai Western Border contributed many more isolates to the available data set than any other site (1291/3087, 42%). However, this rich sample set is itself diverse, contributing information on 31 of the 44 mutant alleles observed in the Asian data set. Moreover, the patients studied were not from a very limited area, but came from sites situated over an extended region of the northwest Thai/eastern Myanmar border. The diversity of propeller mutant genotypes supports the assumption that this large number of isolates from the Thai Western Border is representative of the Eastern part of the whole region that includes Thailand, Cambodia, Lao PDR, and Vietnam.


The slow clearance phenotype that has evolved in P. falciparum parasite populations in response to artemisinin challenge is complex and can involve at least one non-synonymous change in the P. falciparum-specific region of the pfk13 gene, but most frequently, specific single mutations in the propeller region. These different propeller mutations may also compromise parasite fitness, so genomic studies that further define the other genomic changes in parasite populations that support pfk13 mutants to overcome fitness deficits will be needed.

Despite this complexity of the phenotype, the appearance and selection of mutations in the pfk13 propeller region are valuable markers for surveillance of diminished artemisinin responsiveness in parasite populations. This meta-analysis demonstrates that 14 pfk13 propeller mutants, in addition to the 6 previously validated, are associated with prolonged parasite clearance and could be considered for implementation as molecular markers in the Greater Mekong Subregion, as a potential early signal of slow parasite clearance.

The confirmation of the apparent absence outside of Southeast Asia of common pfk13 alleles associated with slow clearance is welcome. However, testing of a more practical and sensitive protocol than persistence of parasites on day 3 that can be a first signal of prolonged parasite clearance is urgently needed. Such an approach should also include vigilance for diminished artemisinin responses signaled by increasing prevalence of parasites with mutant alleles of pfk13 or of other genetic changes that may also diminish parasite responses to artemisinins.



Artemisinin combination therapies




Artesunate + artemisinin combination therapy










Data Management and Statistical Analysis Plan


False negative


False positive


Greater Mekong Subregion




Oxford Tropical Research Ethics Committee

PC1/2 :

Parasite clearance half-life


Parasite clearance estimator tool


Ring-stage survival assay


Tracking resistance to artemisinin collaboration


World Health Organization


  1. 1.

    Noedl H, Se Y, Schaecher K, Smith BL, Socheat D, Fukuda MM. Evidence of artemisinin-resistant malaria in western Cambodia. New Eng J Med. 2008;359(24):2619–20.

  2. 2.

    Dondorp AM, Nosten F, Yi P, Das D, Phyo AP, Tarning J, Lwin KM, Ariey F, Hanpithakpong W, Lee SJ, et al. Artemisinin resistance in Plasmodium falciparum malaria. New Eng J Med. 2009;361(5):455–67.

  3. 3.

    Ariey F, Witkowski B, Amaratunga C, Beghain J, Langlois A-C, Khim N, Kim S, Duru V, Bouchier C, Ma L, et al. A molecular marker of artemisinin-resistant Plasmodium falciparum malaria. Nature. 2014;505(7481):50–5.

  4. 4.

    Plowe CV. Malaria: resistance nailed. Nature. 2014;505(7481):30–1.

  5. 5.

    Cheeseman IH, Miller BA, Nair S, Nkhoma S, Tan A, Tan JC, Al Saai S, Phyo AP, Moo CL, Lwin KM, et al. A major genome region underlying artemisinin resistance in malaria. Science. 2012;336(6077):79–82.

  6. 6.

    Takala-Harrison S, Clark TG, Jacob CG, Cummings MP, Miotto O, Dondorp AM, Fukuda MM, Nosten F, Noedl H, Imwong M, et al. Genetic loci associated with delayed clearance of Plasmodium falciparum following artemisinin treatment in Southeast Asia. Proc Nat Acad Sci U S A. 2013;110(1):240–5.

  7. 7.

    Witkowski B, Amaratunga C, Khim N, Sreng S, Chim P, Kim S, Lim P, Mao S, Sopha C, Sam B, et al. Novel phenotypic assays for the detection of artemisinin-resistant Plasmodium falciparum malaria in Cambodia: in-vitro and ex-vivo drug-response studies. Lancet Infect Dis. 2013;13(12):1043–9.

  8. 8.

    Chotivanich K, Tripura R, Das D, Yi P, Day NPJ, Pukrittayakamee S, Chuor CM, Socheat D, Dondorp AM, White NJ. Laboratory detection of artemisinin-resistant Plasmodium falciparum. Antimicrob Agents Chemother. 2014;58(6):3157–61.

  9. 9.

    4YY8: crystal structure analysis of Kelch protein from Plasmodium falciparum [] Accessed 25 June 2018.

  10. 10.

    Straimer J, Gnadig NF, Witkowski B, Amaratunga C, Duru V, Ramadani AP, Dacheux M, Khim N, Zhang L, Lam S, et al. Drug resistance. K13-propeller mutations confer artemisinin resistance in Plasmodium falciparum clinical isolates. Science. 2015;347(6220):428–31.

  11. 11.

    Ashley EA, Dhorda M, Fairhurst RM, Amaratunga C, Lim P, Suon S, Sreng S, Anderson JM, Mao S, Sam B, et al. Spread of artemisinin resistance in Plasmodium falciparum malaria. New Eng J Med. 2014;371(5):411–23.

  12. 12.

    Nyunt MH, Hlaing T, Oo HW, Tin-Oo LL, Phway HP, Wang B, Zaw NN, Han SS, Tun T, San KK, et al. Molecular assessment of artemisinin resistance markers, polymorphisms in the k13 propeller, and a multidrug-resistance gene in the eastern and western border areas of Myanmar. Clin Infect Dis. 2015;60(8):1208–15.

  13. 13.

    Tun KM, Imwong M, Lwin KM, Win AA, Hlaing TM, Hlaing T, Lin K, Kyaw MP, Plewes K, Faiz MA, et al. Spread of artemisinin-resistant Plasmodium falciparum in Myanmar: a cross-sectional survey of the K13 molecular marker. Lancet Infect Dis. 2015;15(4):415–21.

  14. 14.

    Wang Z, Wang Y, Cabrera M, Zhang Y, Gupta B, Wu Y, Kemirembe K, Hu Y, Liang X, Brashear A, et al. Artemisinin resistance at the China-Myanmar border and association with mutations in the K13 propeller gene. Antimicrob Agents Chemother. 2015;59(11):6952–9.

  15. 15.

    Takala-Harrison S, Jacob CG, Arze C, Cummings MP, Silva JC, Dondorp AM, Fukuda MM, Hien TT, Mayxay M, Noedl H, et al. Independent emergence of artemisinin resistance mutations among Plasmodium falciparum in Southeast Asia. J Infect Dis. 2015;211(5):670–9.

  16. 16.

    Huang F, Takala-Harrison S, Jacob CG, Liu H, Sun X, Yang H, Nyunt MM, Adams M, Zhou S, Xia Z, et al. A single mutation in K13 predominates in southern China and is associated with delayed clearance of Plasmodium falciparum following artemisinin treatment. J Infect Dis. 2015;212(10):1629–35.

  17. 17.

    Ye R, Hu D, Zhang Y, Huang Y, Sun X, Wang J, Chen X, Zhou H, Zhang D, Mungthin M, et al. Distinctive origin of artemisinin-resistant Plasmodium falciparum on the China-Myanmar border. SciRep. 2016;6:20100.

  18. 18.

    Phyo AP, Ashley EA, Anderson TJC, Bozdech Z, Carrara VI, Sriprawat K, Nair S, White MM, Dziekan J, Ling C, et al. Declining efficacy of artemisinin combination therapy against P. falciparum malaria on the Thai–Myanmar Border (2003–2013): the role of parasite genetic factors. Clin Infect Dis. 2016;63(6):784–91.

  19. 19.

    Bonnington CA, Phyo AP, Ashley EA, Imwong M, Sriprawat K, Parker DM, Proux S, White NJ, Nosten F. Plasmodium falciparum Kelch 13 mutations and treatment response in patients in Hpa-Pun District, Northern Kayin State, Myanmar. Malar J. 2017;16(1):480.

  20. 20.

    Myint MK, Rasmussen C, Thi A, Bustos D, Ringwald P, Lin K. Therapeutic efficacy and artemisinin resistance in northern Myanmar: evidence from in vivo and molecular marker studies. Malar J. 2017;16(1):143.

  21. 21.

    Nyunt MH, Soe MT, Myint HW, Oo HW, Aye MM, Han SS, Zaw NN, Cho C, Aung PZ, Kyaw KT, et al. Clinical and molecular surveillance of artemisinin resistant falciparum malaria in Myanmar (2009-2013). Malar J. 2017;16(1):333.

  22. 22.

    Thuy-Nhien N, Tuyen NK, Tong NT, Vy NT, Thanh NV, Van HT, Huong-Thu P, Quang HH, Boni MF, Dolecek C, et al. K13 propeller mutations in Plasmodium falciparum populations in regions of malaria endemicity in Vietnam from 2009 to 2016. Antimicrob Agents Chemother. 2017;61(4):e01578-01516.

  23. 23.

    Taylor SM, Parobek CM, DeConti DK, Kayentao K, Coulibaly SO, Greenwood BM, Tagbor H, Williams J, Bojang K, Njie F, et al. Absence of putative artemisinin resistance mutations among Plasmodium falciparum in sub-Saharan Africa: a molecular epidemiologic study. J Infect Dis. 2015;211(5):680–8.

  24. 24.

    Kamau E, Campino S, Amenga-Etego L, Drury E, Ishengoma D, Johnson K, Mumba D, Kekre M, Yavo W, Mead D, et al. K13-propeller polymorphisms in Plasmodium falciparum parasites from sub-Saharan Africa. J Infect Dis. 2015;211(8):1352–5.

  25. 25.

    Talundzic E, Ndiaye YD, Deme AB, Olsen C, Patel DS, Biliya S, Daniels R, Vannberg FO, Volkman SK, Udhayakumar V, et al. Molecular epidemiology of Plasmodium falciparum kelch13 mutations in senegal determined by using targeted amplicon deep sequencing. Antimicrob Agents Chemother. 2017;61(3):e02116.

  26. 26.

    Mvumbi DM, Bobanga TL, Kayembe J-MN, Mvumbi GL, Situakibanza HN-T, Benoit-Vical F, Melin P, De Mol P, Hayette M-P. Molecular surveillance of Plasmodium falciparum resistance to artemisinin-based combination therapies in the Democratic Republic of Congo. PLoS One. 2017;12(6):e0179142.

  27. 27.

    Menard D, Khim N, Beghain J, Adegnika AA, Shafiul-Alam M, Amodu O, Rahim-Awab G, Barnadas C, Berry A, Boum Y, et al. A worldwide map of Plasmodium falciparum K13-propeller polymorphisms. New Eng J Med. 2016;374(25):2453–64.

  28. 28.

    MalariaGEN. Genomic epidemiology of artemisinin resistant malaria. eLife. 2016;5:e08714.

  29. 29.

    Stepniewska K, Ashley E, Lee SJ, Anstey N, Barnes KI, Binh TQ, D'Alessandro U, Day NP, de Vries PJ, Dorsey G, et al. In vivo parasitological measures of artemisinin susceptibility. J Infect Dis. 2010;201(4):570–9.

  30. 30.

    Flegg JA, Guerin PJ, White NJ, Stepniewska K. Standardizing the measurement of parasite clearance in falciparum malaria: the parasite clearance estimator. Malar J. 2011;10:339.

  31. 31.

    Nkhoma SC, Stepniewska K, Nair S, Phyo AP, McGready R, Nosten F, Anderson TJ. Genetic evaluation of the performance of malaria parasite clearance rate metrics. J Infect Dis. 2013;208(2):346–50.

  32. 32.

    Imwong M, Suwannasin K, Kunasol C, Sutawong K, Mayxay M, Rekol H, Smithuis FM, Hlaing TM, Tun KM, van der Pluijm RW, et al. The spread of artemisinin-resistant Plasmodium falciparum in the Greater Mekong subregion: a molecular epidemiology observational study. Lancet Infect Dis. 2017;17(5):491–7.

  33. 33.

    Amato R, Lim P, Miotto O, Amaratunga C, Dek D, Pearson RD, Almagro-Garcia J, Neal AT, Sreng S, Suon S, et al. Genetic markers associated with dihydroartemisinin-piperaquine failure in Plasmodium falciparum malaria in Cambodia: a genotype-phenotype association study. Lancet Infect Dis. 2017;17(2):164–73.

  34. 34.

    Imwong M, Hien TT, Thuy-Nhien NT, Dondorp AM, White NJ. Spread of a single multidrug resistant malaria parasite lineage (PfPailin) to Vietnam. Lancet Infect Dis. 2017;17(10):1022–3.

  35. 35.

    Amato R, Pearson RD, Almagro-Garcia J, Amaratunga C, Lim P, Suon S, Sreng S, Drury E, Stalker J, Miotto O, et al. Origins of the current outbreak of multidrug-resistant malaria in Southeast Asia: a retrospective genetic study. Lancet Infect Dis. 2018;18(3):337–45.

  36. 36.

    Conrad MD, Bigira V, Kapisi J, Muhindo M, Kamya MR, Havlir DV, Dorsey G, Rosenthal PJ. Polymorphisms in K13 and falcipain-2 associated with artemisinin resistance are not prevalent in Plasmodium falciparum isolated from Ugandan children. PLoS One. 2014;9(8):e105690.

  37. 37.

    Maiga AW, Fofana B, Sagara I, Dembele D, Dara A, Traore OB, Toure S, Sanogo K, Dama S, Sidibe B, et al. No evidence of delayed parasite clearance after oral artesunate treatment of uncomplicated falciparum malaria in Mali. Am J Trop Med Hyg. 2012;87(1):23–8.

  38. 38.

    Ouattara A, Kone A, Adams M, Fofana B, Maiga AW, Hampton S, Coulibaly D, Thera MA, Diallo N, Dara A, et al. Polymorphisms in the K13-propeller gene in artemisinin-susceptible Plasmodium falciparum parasites from Bougoula-Hameau and Bandiagara, Mali. Am J Trop Med Hye. 2015;92(6):1202–6.

  39. 39.

    Chenet SM, Akinyi Okoth S, Huber CS, Chandrabose J, Lucchi NW, Talundzic E, Krishnalall K, Ceron N, Musset L, Macedo de Oliveira A, et al. Independent emergence of the Plasmodium falciparum Kelch propeller domain mutant allele C580Y in Guyana. J Infect Dis. 2016;213(9):1472–5.

  40. 40.

    Tacoli C, Gai PP, Bayingana C, Sifft K, Geus D, Ndoli J, Sendegeya A, Gahutu JB, Mockenhaupt FP. Artemisinin resistance-associated K13 polymorphisms of Plasmodium falciparum in southern Rwanda, 2010-2015. Am J Trop Med Hyg. 2016;95(5):1090–3.

  41. 41.

    Flegg JA, Guerin PJ, Nosten F, Ashley EA, Phyo AP, Dondorp AM, Fairhurst RM, Socheat D, Borrmann S, Bjorkman A, et al. Optimal sampling designs for estimation of Plasmodium falciparum clearance rates in patients treated with artemisinin derivatives. Malar J. 2013;12:411.

  42. 42.

    Artemisinin and artemisinin-based combination therapy resistance [] Accessed 12 June 2018.

  43. 43.

    Stewart LA, Clarke M, Rovers M, Riley RD, Simmonds M, Stewart G, Tierney JF, Group P-ID. Preferred reporting items for systematic review and meta-analyses of individual participant data: the PRISMA-IPD statement. JAMA. 2015;313(16):1657–65.

  44. 44.

    Bethell D, Se Y, Lon C, Tyner S, Saunders D, Sriwichai S, Darapiseth S, Teja-Isavadharm P, Khemawoot P, Schaecher K, et al. Artesunate dose escalation for the treatment of uncomplicated malaria in a region of reported artemisinin resistance: a randomized clinical trial. PLoS One. 2011;6(5):e19283.

  45. 45.

    Das D, Tripura R, Phyo AP, Lwin KM, Tarning J, Lee SJ, Hanpithakpong W, Stepniewska K, Menard D, Ringwald P, et al. Effect of high-dose or split-dose artesunate on parasite clearance in artemisinin-resistant falciparum malaria. Clin Infect Dis. 2013;56(5):e48–58.

  46. 46.

    Hien T, Thuy-Nhien N, Phu N, Boni M, Thanh N, Nha-Ca N, Thai L, Thai C, Toi P, Thuan P, et al. In vivo susceptibility of Plasmodium falciparum to artesunate in Binh Phuoc Province, Vietnam. Malar J. 2012;11(1):355.

  47. 47.

    Mayxay M, Khanthavong M, Chanthongthip O, Imwong M, Lee S, Stepniewska K, Soonthornsata B, Pongvongsa T, Phompida S, Hongvanthong B, et al. No evidence for spread of Plasmodium falciparum artemisinin resistance to Savannakhet Province, southern Laos. Am J Trop Med Hyg. 2012;86(3):403–8.

  48. 48.

    Phyo AP, Nkhoma S, Stepniewska K, Ashley EA, Nair S, McGready R, ler Moo C, Al-Saai S, Dondorp AM, Lwin KM, et al. Emergence of artemisinin-resistant malaria on the western border of Thailand: a longitudinal study. Lancet. 2012;379(9830):1960–6.

  49. 49.

    Starzengruber P, Swoboda P, Fuehrer HP, Khan WA, Hofecker V, Siedl A, Fally M, Graf O, Teja-Isavadharm P, Haque R, et al. Current status of artemisinin-resistant falciparum malaria in South Asia: a randomized controlled artesunate monotherapy trial in Bangladesh. PLoS One. 2012;7(12):e52236.

  50. 50.

    Thriemer K, Hong NV, Rosanas-Urgell A, Phuc BQ, Ha do M, Pockele E, Guetens P, Van NV, Duong TT, Amambua-Ngwa A et al: Delayed parasite clearance after treatment with dihydroartemisinin-piperaquine in Plasmodium falciparum malaria patients in central Vietnam. Antimicrob Agents Chemother 2014, 58(12):7049–7055.

  51. 51.

    Tun KM, Jeeyapant A, Imwong M, Thein M, Aung SSM, Hlaing TM, Yuentrakul P, Promnarate C, Dhorda M, Woodrow CJ, et al. Parasite clearance rates in Upper Myanmar indicate a distinctive artemisinin resistance phenotype: a therapeutic efficacy study. Malar J. 2016;15(1):185.

  52. 52.

    Clinical module: data management and statistical analysis plan version 1.2. [] Accessed 3 June 2018.

  53. 53.

    Correlation between K13 mutations and phenotype study group Statistical Analysis Plan Version 1.0. [] Accessed 12 Jan 2018.

  54. 54.

    Talundzic E, Okoth SA, Congpuong K, Plucinski MM, Morton L, Goldman IF, Kachur PS, Wongsrichanalai C, Satimai W, Barnwell JW, et al. Selection and spread of artemisinin-resistant alleles in Thailand prior to the global artemisinin resistance containment campaign. PLoS Path. 2015;11(4):e1004789.

  55. 55.

    Kobasa T, Talundzic E, Sug-aram R, Boondat P, Goldman IF, Lucchi NW, Dharmarak P, Sintasath D, Fukuda M, Whistler T, et al. Emergence and spread of kelch13 mutations associated with artemisinin resistance in Plasmodium falciparum parasites in 12 Thai provinces from 2007 to 2016. Antimicrob Agents Chemother. 2018;62(4):e02141-02117.

  56. 56.

    WHO, editor. Methods and techniques for assessing exposure to antimalarial drugs in clinical field studies. Geneva: World Health Organization; 2011.

  57. 57.

    WHO: artemisinin resistance and artemisinin-based combination therapy efficacy. [] accessed Nov 26, 2018.

  58. 58.

    Ataide R, Ashley EA, Powell R, Chan J-A, Malloy MJ, O’Flaherty K, Takashima E, Langer C, Tsuboi T, Dondorp AM, et al. Host immunity to Plasmodium falciparum and the assessment of emerging artemisinin resistance in a multinational cohort. Proc Natl Acad Sci U S A. 2017;114(13):3515-20.

  59. 59.

    White N. The parasite clearance curve. Malar J. 2011;10(1):278.

  60. 60.

    Alin MH, Ashton M, Kihamia CM, Mtey GJ, Bjorkman A. Clinical efficacy and pharmacokinetics of artemisinin monotherapy and in combination with mefloquine in patients with falciparum malaria. Br J Clin Pharmacol. 1996;41(6):587–92.

  61. 61.

    Lopera-Mesa TM, Doumbia S, Chiang S, Zeituni AE, Konate DS, Doumbouya M, Keita AS, Stepniewska K, Traore K, Diakite SA, et al. Plasmodium falciparum clearance rates in response to artesunate in Malian children with malaria: effect of acquired immunity. J Infect Dis. 2013;207(11):1655–63.

  62. 62.

    Cooper RA, Conrad MD, Watson QD, Huezo SJ, Ninsiima H, Tumwebaze P, Nsobya SL, Rosenthal PJ. Lack of artemisinin resistance in Plasmodium falciparum in Uganda based on parasitological and molecular assays. Antimicrob Agents Chemother. 2015;59(8):5061–4.

  63. 63.

    Sissoko S, Kone A, Fofana B, Sangaré C, Dembele D, Toure S, Sanogo K, Diakite H, Toure S, Doumbia D, et al. PF3D7_1343700 Kelch propeller (K13-propeller) polymorphisms and artesunate monotherapy efficacy in uncomplicated malaria treatment in Mali. BMJ Global Health. 2017;2(Suppl 2):A18.

  64. 64.

    Straimer J, Gnädig NF, Stokes BH, Ehrenberger M, Crane AA, Fidock DA. Plasmodium falciparum K13 mutations differentially impact ozonide susceptibility and parasite fitness in vitro. mBio. 2017;8(2):e00303-17.

  65. 65.

    Menegon M, Nurahmed AM, Talha AA, Nour BYM, Severini C. Molecular surveillance of antimalarial drug resistance related genes in Plasmodium falciparum isolates from Eritrea. Acta Trop. 2016;157:158–61.

  66. 66.

    Muwanguzi J, Henriques G, Sawa P, Bousema T, Sutherland CJ, Beshir KB. Lack of K13 mutations in Plasmodium falciparum persisting after artemisinin combination therapy treatment of Kenyan children. Malar J. 2016;15(1):36.

  67. 67.

    Guerra M, Neres R, Salgueiro P, Mendes C, Ndong-Mabale N, Berzosa P, de Sousa B, Arez AP. Plasmodium falciparum genetic diversity in continental Equatorial Guinea before and after introduction of artemisinin-based combination therapy. Antimicrob Agents Chemother. 2017;61(1):e00303-17.

  68. 68.

    Ljolje D, Dimbu PR, Kelley J, Goldman I, Nace D, Macaia A, Halsey ES, Ringwald P, Fortes F, Udhayakumar V, et al. Prevalence of molecular markers of artemisinin and lumefantrine resistance among patients with uncomplicated Plasmodium falciparum malaria in three provinces in Angola, 2015. Malar J. 2018;17(1):84.

  69. 69.

    Ocan M, Bwanga F, Okeng A, Katabazi F, Kigozi E, Kyobe S, Ogwal-Okeng J, Obua C. Prevalence of K13-propeller gene polymorphisms among Plasmodium falciparum parasites isolated from adult symptomatic patients in northern Uganda. BMC Infect Dis. 2016;16(1):428.

  70. 70.

    Huang B, Deng C, Yang T, Xue L, Wang Q, Huang S, Su XZ, Liu Y, Zheng S, Guan Y, et al. Polymorphisms of the artemisinin resistant marker (K13) in Plasmodium falciparum parasite populations of Grande Comore Island 10 years after artemisinin combination therapy. Parasit Vectors. 2015;8:634.

  71. 71.

    Torrentino-Madamet M, Collet L, Lepère JF, Benoit N, Amalvict R, Ménard D, Pradines B. K13-propeller polymorphisms in Plasmodium falciparum isolates from patients in Mayotte in 2013 and 2014. Antimicrob Agents Chemother. 2015;59(12):7878–81.

  72. 72.

    Zhang M, Wang C, Otto TD, Oberstaller J, Liao X, Adapa SR, Udenze K, Bronner IF, Casandra D, Mayho M, et al. Uncovering the essential genes of the human malaria parasite Plasmodium falciparum by saturation mutagenesis. Science. 2018;360(6388): eaap7847.

  73. 73.

    Jiang H, Patel JJ, Yi M, Mu J, Ding J, Stephens R, Cooper RA, Ferdig MT, Su X-Z. Genome-wide compensatory changes accompany drug-selected mutations in the Plasmodium falciparum crt gene. PLoS One. 2008;3(6):e2484.

  74. 74.

    Kümpornsin K, Modchang C, Heinberg A, Ekland EH, Jirawatcharadech P, Chobson P, Suwanakitti N, Chaotheing S, Wilairat P, Deitsch KW, et al. Origin of robustness in generating drug-resistant malaria parasites. Mol Biol Evol. 2014;31(7):1649–60.

  75. 75.

    Cheeseman IH, McDew-White M, Phyo AP, Sriprawat K, Nosten F, Anderson TJC. Pooled sequencing and rare variant association tests for identifying the determinants of emerging drug resistance in malaria parasites. Mol Biol Evol. 2015;32(4):1080–90.

  76. 76.

    Miotto O, Amato R, Ashley EA, MacInnis B, Almagro-Garcia J, Amaratunga C, Lim P, Mead D, Oyola SO, Dhorda M, et al. Genetic architecture of artemisinin-resistant Plasmodium falciparum. Nat Genet. 2015;47(3):226–34.

  77. 77.

    Dwivedi A, Reynes C, Kuehn A, Roche DB, Khim N, Hebrard M, Milanesi S, Rivals E, Frutos R, Menard D, et al. Functional analysis of Plasmodium falciparum subpopulations associated with artemisinin resistance in Cambodia. Malar J. 2017;16(1):493.

  78. 78.

    Anderson TJ, Nair S, McDew-White M, Cheeseman IH, Nkhoma S, Bilgic F, McGready R, Ashley E, Pyae Phyo A, White NJ, et al. Population parameters underlying an ongoing soft sweep in southeast Asian malaria parasites. Mol Biol Evol. 2017;34(1):131–44.

  79. 79.

    Cerqueira GC, Cheeseman IH, Schaffner SF, Nair S, McDew-White M, Phyo AP, Ashley EA, Melnikov A, Rogov P, Birren BW, et al. Longitudinal genomic surveillance of Plasmodium falciparum malaria parasites reveals complex genomic architecture of emerging artemisinin resistance. Genome Biol. 2017;18(1):78.

  80. 80.

    Tilley L, Straimer J, Gnadig NF, Ralph SA, Fidock DA. Artemisinin action and resistance in Plasmodium falciparum. Trends Parasitol. 2016;32(9):682–96.

  81. 81.

    Mok S, Ashley EA, Ferreira PE, Zhu L, Lin Z, Yeo T, Chotivanich K, Imwong M, Pukrittayakamee S, Dhorda M, et al. Drug resistance. Population transcriptomics of human malaria parasites reveals the mechanism of artemisinin resistance. Science. 2015;347(6220):431–5.

  82. 82.

    Mbengue A, Bhattacharjee S, Pandharkar T, Liu H, Estiu G, Stahelin RV, Rizk SS, Njimoh DL, Ryan Y, Chotivanich K, et al. A molecular mechanism of artemisinin resistance in Plasmodium falciparum malaria. Nature. 2015;520(7549):683–7.

  83. 83.

    Mukherjee A, Bopp S, Magistrado P, Wong W, Daniels R, Demas A, Schaffner S, Amaratunga C, Lim P, Dhorda M, et al. Artemisinin resistance without pfkelch13 mutations in Plasmodium falciparum isolates from Cambodia. Malar J. 2017;16(1):195.

  84. 84.

    Rocamora F, Zhu L, Liong KY, Dondorp A, Miotto O, Mok S, Bozdech Z. Oxidative stress and protein damage responses mediate artemisinin resistance in malaria parasites. PLoS Path. 2018;14(3):e1006930.

  85. 85.

    Siddiqui FA, Cabrera M, Wang M, Brashear A, Kemirembe K, Wang Z, Miao J, Chookajorn T, Yang Z, Cao Y, et al. Plasmodium falciparum Falcipain-2a polymorphisms in Southeast Asia and their association with artemisinin resistance. J Infect Dis. 2018;218(3):434–42.

  86. 86.

    Straimer J, Gnadig NF, Stokes BH, Ehrenberger M, Crane AA, Fidock DA. Plasmodium falciparum K13 mutations differentially impact ozonide susceptibility and parasite fitness in vitro. mBio. 2017;8(2):e00303-17.

  87. 87.

    Ikeda M, Kaneko M, Tachibana SI, Balikagala B, Sakurai-Yatsushiro M, Yatsushiro S, Takahashi N, Yamauchi M, Sekihara M, Hashimoto M, et al. Artemisinin-resistant Plasmodium falciparum with high survival rates, Uganda, 2014-2016. Emerg Infect Dis. 2018;24(4):718–26.

Download references


We appreciate the careful review of the manuscript by Pascal Ringwald and the advice and contributions of Paul Newton and the Structural Genomics Consortium who have publically shared the data on the structure of the propeller region of the K13 protein. The members of the WorldWide Antimalarial Resistance Network (WWARN) K13 Genotype-Phenotype Correlation Study Group are the authors of this paper:

Chanaki Amaratunga

Laboratory of Malaria and Vector Research, Division of Intramural Research, National Institute of Allergy and Infectious Diseases, National Institutes of Health, Rockville, Maryland, USA

Voahangy Hanitriniaina Andrianaranjaka

Malaria Research Unit, Institut Pasteur de Madagascar, Antananarivo, Madagascar, and Faculté des Sciences, Université d’Antananarivo, Antananarivo, Madagascar

Elizabeth Ashley

Myanmar-Oxford Clinical Research Unit (MOCRU), Yangon, Myanmar, and Centre for Tropical Medicine and Global Health, University of Oxford, UK

Delia Bethell

Armed Forces Research Institute of Medical Sciences, Bangkok, Thailand

Anders Björkman

Department of Molecular, Tumor, and Cell Biology, Karolinska Institutet, Stockholm, Sweden

Craig A. Bonnington

Shoklo Malaria Research Unit, Mae Sot, Thailand

Roland A. Cooper

Department of Natural Sciences and Mathematics, Dominican University of California, San Rafael, CA, USA

Mehul Dhorda

WWARN, Centre for Tropical Medicine, Nuffield Department of Clinical Medicine,University of Oxford, Oxford, UK

Arjen Dondorp

Mahidol-Oxford Research Unit, Faculty of Tropical Medicine, Mahidol University, Bangkok, Thailand, and WWARN, Centre for Tropical Medicine, Nuffield Department of Clinical Medicine, University of Oxford, Oxford, UK

Annette Erhart

Department of Public Health, ITM Antwerp, Belgium, and Institute of Tropical Medicine, MRC Unit The Gambia, Fajara, Gambia

Rick M. Fairhurst

Laboratory of Malaria and Vector Research, Division of Intramural Research, National Institute of Allergy and Infectious Diseases, National Institutes of Health, Rockville, Maryland, USA

Abul Faiz

Dev Care Foundation, Dhaka, Bangladesh

Caterina Fanello

Centre for Tropical Medicine, Nuffield Department of Clinical Medicine,University of Oxford, Oxford, UK, and Mahidol-Oxford Research Unit, Bangkok, Thailand

Mark M. Fukuda

Armed Forces Research Institute of Medical Sciences, Bangkok, Thailand

Philippe Guérin

WWARN, Centre for Tropical Medicine, Nuffield Department of Clinical Medicine, University of Oxford, Oxford, UK

Rob Hooft van Huijsduijnen

Medicines for Malaria Venture, Meyrin, Geneva, Switzerland

Tran Tinh Hien

Centre for Tropical Medicine Oxford University, Clinical Research Unit, Ho Chi Minh City, Vietnam

NV Hong

National Institute of Malariology, Parasitology and Entomology, Hanoi, Vietnam

Ye Htut

Department of Medical Research, Lower Myanmar, Yangon, Myanmar

Fang Huang

National Institute of Parasitic Diseases, Chinese Center for Disease Control and Prevention, Shanghai, China

Georgina Humphreys

WWARN, Centre for Tropical Medicine, Nuffield Department of Clinical Medicine,University of Oxford, Oxford, UK

Mallika Imwong

Department of Molecular Tropical Medicine and Genetics Faculty of Tropical Medicine, Mahidol University, Bangkok, Thailand and Mahidol–Oxford Tropical Medicine Research Unit, Faculty of Tropical Medicine, Mahidol University, Bangkok, Thailand

Kalynn Kennon

WWARN, Centre for Tropical Medicine, Nuffield Department of Clinical Medicine,University of Oxford, Oxford,UK

Pharath Lim

Laboratory of Malaria and Vector Research, Division of Intramural Research, National Institute of Allergy and Infectious Diseases, National Institutes of Health, Rockville, Maryland, USA

Khin Lin

Department of Medical Research, (PYIN OO LWIN BRANCH), Myanmar

Chanthap Lon

Armed Forces Research Institute of Medical Sciences, Bangkok, Thailand 

Andreas Mårtensson

Department of Women’s and Children’s Health, International Maternal and Child health (IMCH) Uppsala University, Uppsala, Sweden

Mayfong Mayxay

Lao-Oxford-Mahosot Hospital-Wellcome Trust Research Unit (LOMWRU), Vientiane, Lao PDR, and Faculty of Postgraduate Studies, University of Health Sciences, Ministry of Health, Vientiane, Lao PDR, and Centre for Tropical Medicine and Global Health, Nuffield Department of Medicine, Churchill Hospital, Oxford, UK

Olugbenga Mokuolu

Department of Paediatrics and Child Health, College of Health Sciences, University of Ilorin, Ilorin, Nigeria, and Centre for Malaria and Other Tropical Diseases Care, University of Ilorin Teaching Hospital, Ilorin, Nigeria

Ulrika Morris

Department of Molecular, Tumor, and Cell Biology, Karolinska Institutet, Stockholm, Sweden

Billy E. Ngasala

Department of Parasitology and Medical Entomology, Muhimbili University of Health and Allied Sciences, Dar es Salaam, Tanzania

Alfred Amambua-Ngwa

Institute of Tropical Medicine, MRC Unit The Gambia, Fajara, Gambia

Harald Noedl

Institute of Specific Prophylaxis and Tropical Medicine, Medical University of Vienna, Vienna, Austria

François Nosten

Shoklo Malaria Research Unit, Mae Sot, Thailand, Mahidol-Oxford Tropical Medicine Research Unit, Faculty of Tropical Medicine, Mahidol University, Bangkok, Thailand, and Centre for Tropical Medicine, Nuffield Department of Clinical Medicine, University of Oxford, Oxford, UK

Marie Onyamboko

Kinshasa School of Public Health, Kinshasa, DRC and Mahidol-Oxford Research Unit, Bangkok, Thailand

Aung Pyae Phyo

Shoklo Malaria Research Unit, Mae Sot, Thailand, Mahidol-Oxford Tropical Medicine Research Unit, Faculty of Tropical Medicine, Mahidol University, Bangkok, Thailand and Centre for Tropical Medicine, Nuffield Department of Clinical Medicine, University of Oxford, Oxford, UK

Christopher V. Plowe

Duke Global Health Institute, Duke University, Durham, NC, USA

Sasithon Pukrittayakamee

Department of Clinical Tropical Medicine, Mahidol University, Bangkok, Thailand and Royal Society of Thailand, Sanam Suea Pa, Khet Dusit, Bangkok, Thailand

Milijaona Randrianarivelojosia

Malaria Research Unit, Institut Pasteur de Madagascar, Antananarivo, Madagascar and Faculté des Sciences, Université de Toliara, Toliara, Madagascar

Philip J. Rosenthal

Department of Medicine and Division of HIV, Infectious Diseases, and Global Medicine, University of California, San Francisco, San Francisco, CA, USA

David L. Saunders

Armed Forces Research Institute of Regenerative Medicine, Bangkok, Thailand and US Army Medical Materiel Development Activity, Fort Detrick, Maryland, USA

Carol Hopkins Sibley

Department of Genome Sciences, University of Washington, Seattle, WA, USA and WWARN, Centre for Tropical Medicine, Nuffield Department of Clinical Medicine, University of Oxford, Oxford, UK

Frank Smithuis

Myanmar Oxford Clinical Research Unit, Yangon, Myanmar

Michele D. Spring

Department of Immunology and Medicine, Armed Forces Research Institute of Medical Sciences, Bangkok, Thailand

Paul Sondo

WWARN, Centre for Tropical Medicine, Nuffield Department of Clinical Medicine,University of Oxford, Oxford, UK and Clinical Research Unit of Nanoro (CRUN), Burkina Faso

Sokunthea Sreng

National Center for Parasitology, Entomology and Malaria Control, Phnom Penh, Cambodia

Peter Starzengruber

Institute of Specific Prophylaxis and Tropical Medicine, Medical University of Vienna, Vienna, Austria and Division of Clinical Microbiology, Department of Laboratory Medicine, Medical University of Vienna, Vienna, Austria

Kasia Stepniewska

WWARN, Centre for Tropical Medicine and Global Health, University of Oxford, Oxford, UK

Seila Suon

National Center for Parasitology, Entomology and Malaria Control, Phnom Penh, Cambodia

Shannon Takala-Harrison

Division of Malaria Research, Institute for Global Health, University of Maryland School of Medicine, Baltimore, MD, USA

Kamala Thriemer

Institute of Tropical Medicine, Antwerp, Belgium and Menzies School of Health Research, Darwin, Australia

Nguyen Thuy-Nhien

Centre for Tropical Medicine Oxford University Clinical Research Unit. Ho Chi Minh City, Vietnam

Kyaw Myo Tun

Defense Services Medical Academy, Yangon, Myanmar and Myanmar-Oxford Clinical Research Unit, Yangon, Myanmar

Nicholas J. White

Mahidol-Oxford Research Unit, Faculty of Tropical Medicine, Mahidol University, Bangkok, Thailand and Centre for Tropical Medicine, Nuffield Department of Clinical Medicine, University of Oxford, Oxford, UK

Charles Woodrow

Mahidol-Oxford Research Unit, Faculty of Tropical Medicine, Mahidol University, Bangkok,Thailand and Centre for Tropical Medicine, Nuffield Department of Clinical Medicine, University of Oxford, Oxford, UK


This pooled analysis was funded by grants from the ExxonMobil Foundation and the Bill & Melinda Gates Foundation. Funders had no input on the planning, analysis, or presentation of the study. Information on funding support for the primary data collection and analysis can be found in the respective referenced manuscripts; these publications are also listed in Additional file 2: Table S1.

Data access

The data are available for access via the WorldWide Antimalarial Resistance Network ( Requests for access will be reviewed by a Data Access Committee to ensure that use of data protects the interests of the participants and researchers according to the WWARN principles of data sharing ( Requests can be submitted by email to via the Data Access Form available at WWARN is registered with the Registry of Research Data Repositories (

Author information

KS, CHS, PG, NJW, MD, EA, and FN conceived/designed the meta-analysis. NJW, EA, AB, MD, MI, AE, RAC, PJR, AAN, RMF, AM, UM, HN, FN, APP, MR, DLS, KT, and KMT conceived/designed the original experiments. CHS and KS wrote the first draft. KS, CHS, PG, NJW, PJR, AB, EA, CVP, MD, RAC, CA, KT, DB, RMF, AE, FH, RHH, STH, KK, and P Sondo contributed to the writing of the manuscript. EA, AB, AAN, RMF, AM, UM, HN, FN, MI, APP, MR, DLS, KT, KMT, CA, VHA, DB, CAB, MD, AE, AF, CF, MMF, RHH, TTH, NVH, YH, FH, CL, MM, OM, BEN, MO, CVP, SP, MDS, P Starzengruber, STH, NTN, CW, S Sreng, and S Suon contributed the reagents, data, or materials. EA, AB, RMF, AM, UM, FN, APP, MR, DLS, KT, KMT, CA, VHA, DB, CAB, AE, AF, CF, MMF, NVH, YH, FH, CL, MM, OM, BEN, MO, SP, MDS, NTN, S Sreng, and S Suon enrolled patients/performed experiments. EA, AAN, MI, RMF, CA, PL, FN, FH, APP, NTN, and CW analyzed the data in original studies. KS, CHS, PG, KK, P Sondo, and GH performed the IPD meta-analysis. All authors read and approved the manuscript.

Ethics declarations

Ethics approval and consent to participate

All data included in this analysis were obtained from studies conducted with ethical approvals from the country of origin. Use of existing data which are fully anonymized and which researchers cannot trace back to identifiable individuals does not require the review of the Ethics Committee under the guidelines of the Oxford Central University Research Ethics Committee.

Consent for publication

Not required

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Details of the systematic literature review. (PDF 213 kb)

Additional file 2:

Supplementary Tables S1-S3. (XLSX 66 kb)

Additional file 3:

Supplementary Figures S1- S9. (PDF 1168 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Amaratunga, C., Andrianaranjaka, V.H., Ashley, E. et al. Association of mutations in the Plasmodium falciparum Kelch13 gene (Pf3D7_1343700) with parasite clearance rates after artemisinin-based treatmentsa WWARN individual patient data meta-analysis. BMC Med 17, 1 (2019) doi:10.1186/s12916-018-1207-3

Download citation