Plasmodium falciparum K13 mutations in Africa and Asia impact artemisinin resistance and parasite fitness

  1. Barbara H Stokes
  2. Satish K Dhingra
  3. Kelly Rubiano
  4. Sachel Mok
  5. Judith Straimer
  6. Nina F Gnädig
  7. Ioanna Deni
  8. Kyra A Schindler
  9. Jade R Bath
  10. Kurt E Ward
  11. Josefine Striepen
  12. Tomas Yeo
  13. Leila S Ross
  14. Eric Legrand
  15. Frédéric Ariey
  16. Clark H Cunningham
  17. Issa M Souleymane
  18. Adama Gansané
  19. Romaric Nzoumbou-Boko
  20. Claudette Ndayikunda
  21. Abdunoor M Kabanywanyi
  22. Aline Uwimana
  23. Samuel J Smith
  24. Olimatou Kolley
  25. Mathieu Ndounga
  26. Marian Warsame
  27. Rithea Leang
  28. François Nosten
  29. Timothy JC Anderson
  30. Philip J Rosenthal
  31. Didier Ménard
  32. David A Fidock  Is a corresponding author
  1. Department of Microbiology and Immunology, Columbia University Irving Medical Center, United States
  2. Department of Microbiology and Immunology, University of Otago, New Zealand
  3. Malaria Genetics and Resistance Unit, Institut Pasteur, INSERM U1201, CNRS ERL9195, France
  4. Institut Cochin, INSERM U1016, Université Paris Descartes, France
  5. Department of Genetics, University of North Carolina at Chapel Hill, United States
  6. Programme National de Lutte Contre le Paludisme au Tchad, Chad
  7. Centre National de Recherche et de Formation sur le Paludisme, Burkina Faso
  8. Laboratoire de Parasitologie, Institut Pasteur de Bangui, Central African Republic
  9. University Teaching Hospital of Kamenge, Burundi
  10. Ifakara Health Institute, United Republic of Tanzania
  11. Malaria and Other Parasitic Diseases Division, Rwanda Biomedical Centre, Rwanda
  12. National Malaria Control Program, Sierra Leone
  13. National Malaria Control Program, Gambia
  14. Programme National de Lutte Contre le Paludisme, Democratic Republic of the Congo
  15. School of Public Health and Community Medicine, University of Gothenburg, Sweden
  16. National Center for Parasitology, Entomology & Malaria Control, Cambodia
  17. Shoklo Malaria Research Unit, Mahidol-Oxford Tropical Medicine Research Unit, Faculty of Tropical Medicine, Mahidol University, Thailand
  18. Centre for Tropical Medicine and Global Health, Nuffield Department of Medicine, University of Oxford, United Kingdom
  19. Texas Biomedical Research Institute, United States
  20. Department of Medicine, University of California, San Francisco, United States
  21. Division of Infectious Diseases, Department of Medicine, Columbia University Irving Medical Center, United States

Abstract

The emergence of mutant K13-mediated artemisinin (ART) resistance in Plasmodium falciparum malaria parasites has led to widespread treatment failures across Southeast Asia. In Africa, K13-propeller genotyping confirms the emergence of the R561H mutation in Rwanda and highlights the continuing dominance of wild-type K13 elsewhere. Using gene editing, we show that R561H, along with C580Y and M579I, confer elevated in vitro ART resistance in some African strains, contrasting with minimal changes in ART susceptibility in others. C580Y and M579I cause substantial fitness costs, which may slow their dissemination in high-transmission settings, in contrast with R561H that in African 3D7 parasites is fitness neutral. In Cambodia, K13 genotyping highlights the increasing spatio-temporal dominance of C580Y. Editing multiple K13 mutations into a panel of Southeast Asian strains reveals that only the R561H variant yields ART resistance comparable to C580Y. In Asian Dd2 parasites C580Y shows no fitness cost, in contrast with most other K13 mutations tested, including R561H. Editing of point mutations in ferredoxin or mdr2, earlier associated with resistance, has no impact on ART susceptibility or parasite fitness. These data underline the complex interplay between K13 mutations, parasite survival, growth and genetic background in contributing to the spread of ART resistance.

Introduction

Despite recent advances in chemotherapeutics, diagnostics and vector control measures, malaria continues to exert a significant impact on human health (Hanboonkunupakarn and White, 2020). In 2019, cases were estimated at 229 million, resulting in 409,000 fatal outcomes, primarily in Sub-Saharan Africa as a result of Plasmodium falciparum infection (WHO, 2020). This situation is predicted to worsen as a result of the ongoing SARS-CoV-2 pandemic that has compromised malaria treatment and prevention measures (Sherrard-Smith et al., 2020). In the absence of an effective licensed malaria vaccine, control and elimination strategies are critically reliant on the continued clinical efficacy of first-line artemisinin-based combination therapies (ACTs) (White et al., 2014). These ACTs pair fast-acting artemisinin (ART) derivatives with partner drugs such as lumefantrine, amodiaquine, mefloquine, or piperaquine (PPQ). ART derivatives can reduce the biomass of drug-sensitive parasites by up to 10,000-fold within 48 hr (the duration of one intra-erythrocytic developmental cycle); however, these derivatives are rapidly metabolized in vivo. Longer-lasting, albeit slower-acting, partner drugs are co-administered to reduce the selective pressure for ART resistance and to clear residual parasitemias (Eastman and Fidock, 2009).

P. falciparum resistance to ART derivatives has now swept across Southeast (SE) Asia, having first emerged a decade ago in western Cambodia (Dondorp et al., 2009; Noedl et al., 2009; Ariey et al., 2014; Imwong et al., 2020). Clinically, ART resistance manifests as delayed clearance of circulating asexual blood stage parasites following treatment with an ACT, but does not result in treatment failure as long as the partner drug remains effective. The accepted threshold for resistance is a parasite clearance half-life (the time required for the peripheral blood parasite density to decrease by 50%) of >5.5 hr. Sensitive parasites are typically cleared in <2–3 hr (WHO, 2019). Resistance can also be evidenced as parasite-positive blood smears on day three post initiation of treatment. In vitro, ART resistance manifests as increased survival of tightly synchronized early ring-stage parasites (0–3 hr post invasion) exposed to a 6 hr pulse of 700 nM dihydroartemisinin (DHA, the active metabolite of all ARTs used clinically) in the ring-stage survival assay (RSA) (Witkowski et al., 2013; Ariey et al., 2014). Recently, ART-resistant strains have also acquired resistance to PPQ, which is widely used in SE Asia as a partner drug in combination with DHA (Wicht et al., 2020). Failure rates following DHA-PPQ treatment now exceed 50% in parts of Cambodia, Thailand, and Vietnam (van der Pluijm et al., 2019).

In vitro selections, supported by clinical epidemiological data, have demonstrated that ART resistance is primarily determined by mutations in the beta-propeller domain of the P. falciparum Kelch protein K13, also known as Kelch13 (Ariey et al., 2014; Ashley et al., 2014; MalariaGEN Plasmodium falciparum Community Project, 2016; Ménard et al., 2016; Siddiqui et al., 2020). Recent evidence suggests that these mutations result in reduced endocytosis of host-derived hemoglobin and thereby decrease the release of Fe2+-heme that serves to activate ART, thus reducing the drug's potency (Yang et al., 2019; Birnbaum et al., 2020). Mutations in other genes including ferredoxin (fd) and multidrug resistance protein 2 (mdr2) have also been associated with ART resistance in K13 mutant parasites, leading to the suggestion that they may contribute to a multigenic basis of resistance and/or parasite fitness, or serve as genetic markers of founder populations (Miotto et al., 2015).

In SE Asia, the most prevalent K13 mutation is C580Y, which associates with delayed clearance in vivo (Ariey et al., 2014; Ashley et al., 2014; MalariaGEN Plasmodium falciparum Community Project, 2016; Ménard et al., 2016; Imwong et al., 2017). This mutation also mediates ART resistance in vitro, as demonstrated by RSAs performed on gene-edited parasites (Ghorbal et al., 2014; Straimer et al., 2015; Straimer et al., 2017; Mathieu et al., 2020; Uwimana et al., 2020). Other studies have documented the emergence of nearly 200 other K13 mutations, both in SE Asia and in other malaria-endemic regions, including the Guiana Shield and the western Pacific (MalariaGEN Plasmodium falciparum Community Project, 2016; Ménard et al., 2016; Das et al., 2019; WWARN K13 Genotype-Phenotype Study Group, 2019; Mathieu et al., 2020; Miotto et al., 2020). Aside from C580Y, however, only a handful of other K13 mutations (N458Y, M476I, Y493H, R539T, I543T, and R561H) have been validated by gene-editing experiments as conferring ART resistance in vitro (Straimer et al., 2015; Siddiqui et al., 2020). Nonetheless, multiple mutations in this gene have been associated with the clinical delayed clearance phenotype and have been proposed as candidate markers of ART resistance (WWARN K13 Genotype-Phenotype Study Group, 2019; WHO, 2019).

Here, we define the role of a panel of K13 mutations identified in patient isolates, and address the key question of whether these mutations can confer resistance in African and Asian strains. We include the K13 R561H mutation, first associated with delayed parasite clearance in SE Asia (Ashley et al., 2014; Phyo et al., 2016), and very recently identified at up to 13% prevalence in certain districts in Rwanda (Uwimana et al., 2020; Bergmann et al., 2021; Uwimana et al., 2021). We also assess the impact of the parasite genetic background on in vitro phenotypes, including mutations in ferredoxin and mdr2 that were earlier associated with resistance (Miotto et al., 2015). Our results show that K13 mutations can impart ART resistance across multiple Asian and African strains, at levels that vary widely depending on the mutation and the parasite genetic background. Compared with K13 mutant Asian parasites, we observed stronger in vitro fitness costs in most K13 edited African strains, which might predict a slower dissemination of ART resistance in high-transmission African settings. Nonetheless, our data highlight the threat of the R561H mutation emerging in Rwanda, which confers elevated RSA survival and a minimal fitness cost in African 3D7 parasites.

Results

Non-synonymous K13 mutations are present at low frequencies in Africa

To examine the status of K13 mutations across Africa, we analyzed K13 beta-propeller domain sequences in 3257 isolates from eleven malaria-endemic African countries, including The Gambia, Sierra Leone, and Burkina Faso in West Africa; Chad, Central African Republic, Republic of the Congo, and Equatorial Guinea in Central Africa; and Burundi, Tanzania, Rwanda, and Somalia in East Africa. Samples were collected between 2011 and 2019, with most countries sampled across multiple years. A total of 1038 (32%) samples originated from The Gambia, Republic of the Congo, or Burundi, and have not been previously reported, whereas the remaining samples, including those from Rwanda, have been published (Figure 1—source data 1; Supplementary file 1).

Of all samples, 98% (3179) were K13 wild-type, that is, they matched the 3D7 (African) reference sequence or harbored a synonymous (non-coding) mutation. For individual countries, the percentage of K13 wild-type samples ranged from 95% to 100% (Figure 1; Figure 1—source data 1). In total, we identified 35 unique non-synonymous mutations in K13. Of these, only two have been validated as resistance mediators in the SE Asian Dd2 strain: the M476I mutation initially identified from long-term ART selection studies, and the R561H mutation observed in SE Asia and now in Rwanda (Ariey et al., 2014; Straimer et al., 2015; Uwimana et al., 2020).

Frequency and distribution of K13 alleles in eleven African countries.

Map of Africa with pie charts representing the proportions of sequenced samples per country that harbor the K13 wild-type sequence (3D7 reference), the R561H variant (the most commonly identified mutation, unique to Rwanda; see inset), or another less frequent non-synonymous K13 mutation. Sample sizes and years of sample collection are indicated. Mutations and numbers of African samples sequenced per country, and prior citations as appropriate, are listed in Figure 1—source data 1.

Figure 1—source data 1

Distribution of K13 alleles over time in African countries (2011–2019).

https://cdn.elifesciences.org/articles/66277/elife-66277-fig1-data1-v3.xlsx

Of the 35 non-synonymous mutations, only two were present in ≥ten samples: R561H (n=20, found only in Rwanda, sampled from 2012 to 2015; Uwimana et al., 2020), and A578S (n=10; observed in four African countries across multiple years). Previously, A578S was shown not to confer in vitro resistance in Dd2 (Ménard et al., 2016). In the set of 927 genotyped Rwandan isolates, R561H accounted for 44% of mutant samples and 2% of all samples (Figure 1 inset).

K13 R561H, M579I, and C580Y mutations can confer in vitro artemisinin resistance in African parasites

To test whether the K13 R561H mutation could mediate ART resistance in African strains, we developed a CRISPR/Cas9-mediated K13 editing strategy (Supplementary file 2) to introduce this mutation into 3D7 and F32 parasites. On the basis of whole-genome sequence analysis of African isolates, 3D7 was recently shown to segregate phylogenetically with parasites from Rwanda (Uwimana et al., 2020). F32 was derived from an isolate from Tanzania (Witkowski et al., 2010). We also tested the C580Y mutation that predominates in SE Asia, as well as the M579I mutation earlier identified in a P. falciparum-infected migrant worker in Equatorial Guinea who displayed delayed parasite clearance following ACT treatment (Lu et al., 2017). The positions of these residues are highlighted in the K13 beta-propeller domain structure shown in Supplementary file 3. For 3D7, F32 and other lines used for this study, geographic origins and genotypes at drug resistance loci are described in Table 1 and Supplementary file 4. All parental lines were cloned by limiting dilution prior to transfection. Edited parasites were identified by PCR and Sanger sequencing, and cloned. These and other edited parasite lines used herein are described in Supplementary file 5.

Table 1
Plasmodium falciparum lines employed herein.
ParasiteOriginYearK13Resistance
3D7WTAfrica1981WT--
F32WTTanzania1982WT--
UG659WTUganda2007WTCQ, SP
UG815WTUganda2008WTCQ, SP
Dd2WTIndochina1980WTCQ, MFQ, SP
Cam3.IIWTCambodia2010WTCQ, SP
CamWTC580YCambodia2010C580YART, CQ, SP
RF7C580YCambodia2012C580YART, CQ, PPQ, SP
Thai1WTThailand2003WTCQ, SP
Thai2WTThailand2004WTCQ, MFQ, SP
Thai3WTThailand2003WTCQ, SP
Thai4WTThailand2003WTCQ, SP
Thai5WTThailand2011WTCQ, SP
Thai6E252QThailand2008E252QART (low), CQ, MFQ, SP
Thai7E252QThailand2010E252QART (low), CQ, MFQ, SP
  1. Parasite superscripts refer to the K13 sequence.

    ART, artemisinin; CQ, chloroquine; MFQ, mefloquine; PPQ, piperaquine; SP, sulfadoxine/pyrimethamine; WT, wild type.

RSAs, used to measure in vitro ART susceptibility, revealed a wide range of mean survival values for K13 mutant lines. For 3D7 parasites, the highest RSA survival rates were observed with 3D7R561H parasites, which averaged 6.6% RSA survival. For the 3D7M579I and 3D7C580Y lines, mean RSA survival rates were both 4.8%, a three- to fourfold increase relative to the 3D7WT line. No elevated RSA survival was seen in a 3D7 control line (3D7ctrl) that expresses only the silent shield mutations used at the guide RNA (gRNA) cut site (Figure 2A; Figure 2—source data 1). Western blots performed on tightly synchronized ring-stage parasites revealed an ~30% reduction in K13 protein expression levels in these three K13 mutant lines relative to the parental 3D7WT line (Figure 2—figure supplement 1; Figure 2—figure supplement 1—source data 1).

Figure 2 with 1 supplement see all
Gene-edited mutant K13 African parasites display variable levels of RSA survival.

(A–D) RSA survival rates for (A) 3D7 (Africa), (B) F32 (Tanzania), (C) UG659 (Uganda), or (D) UG815 (Uganda) K13 wild-type parental lines and CRISPR/Cas9-edited K13 R561H, M579I, or C580Y mutant clones. Unedited parental lines are described in Table 1 and Supplementary file 4. For 3D7, we also included a K13 wild-type control (ctrl) line harboring silent shield mutations at the K13 gRNA cut site. Results show the percentage of early ring-stage parasites (0–3 hr post invasion) that survived a 6 hr pulse of 700 nM DHA, relative to DMSO-treated parasites assayed in parallel. Percent survival values are shown as means ± SEM (detailed in Figure 2—source data 1). Results were obtained from three to eight independent experiments, each performed in duplicate. p Values were determined by unpaired t tests and were calculated for K13 mutant lines relative to their isogenic wild-type lines. ** p<0.01; *** p<0.001; **** p<0.0001.

Figure 2—source data 1

Ring-stage survival (RSA) assay data for K13 edited African parasites and controls.

https://cdn.elifesciences.org/articles/66277/elife-66277-fig2-data1-v3.xlsx

Interestingly, for F32 parasites, the introduction of K13 mutations yielded no significant increases in RSA survival, with survival rates in the range of 0.3% to 0.5% for lines expressing R561H, M579I, C580Y, or wild-type K13. (Figure 2B). Previously we reported that introduction of M476I into F32 parasites resulted in a modest gain of resistance (mean survival of 1.7%), while this same mutation conferred RSA survival levels of ~10% in edited Dd2 parasites (Straimer et al., 2015). These data suggest that while K13 mutations differ substantially in their impact on ART susceptibility, there is an equally notable contribution of the parasite genetic background.

We next introduced M579I and C580Y into the cloned Ugandan isolates UG659 and UG815. Editing of both mutations into UG659 yielded moderate RSA survival rates (means of 6.3% and 4.7% for UG659M579I or UG659C580Y respectively, vs 1.0% for UG659WT; Figure 2C). These values resembled our results for 3D7. Strikingly, introducing K13 M579I or C580Y into UG815 yielded the highest rates of in vitro resistance, with mean survival levels reaching ~12% in both UG815M579I and UG815C580Y. These results were confirmed in a second independent clone of UG815M579I (Figure 2D). M579I and C580Y also conferred equivalent levels of resistance in edited Dd2 parasites (RSA survival rates of 4.0% and 4.7%, respectively; Figure 2—source data 1). These data show that mutant K13-mediated ART resistance in African parasites can be achieved in some strains at levels comparable to or above those seen in SE Asian parasites.

K13 C580Y, M579I, and R561H mutations are associated with variable in vitro fitness costs in African parasites

To examine the relationship between resistance and fitness in African parasites harboring K13 mutations, we developed an in vitro fitness assay that uses quantitative real-time PCR (qPCR) for allelic discrimination. Assays were conducted by pairing K13 wild-type lines (3D7, F32, UG659, and UG815) with their isogenic K13 edited R561H, M579I, or C580Y counterparts.

Assays were initiated with tightly synchronized trophozoites that were mixed in 1:1 ratios of wild-type to mutant isogenic parasites, and cultures were maintained over a period of 36 days (~18 generations of asexual blood stage growth). Cultures were sampled every four days for genomic DNA (gDNA) preparation and qPCR analysis. TaqMan probes specific to the K13 wild-type or mutant (R561H, M579I, or C580Y) alleles were used to quantify the proportion of each allele.

Results showed that the K13 M579I and C580Y mutations each conferred significant fitness defects in most strains tested, with the proportions of K13 mutant lines declining over time. For both mutations, the largest reductions were observed with edited 3D7 or UG815 parasites. In contrast, these mutations exerted a minimal impact on fitness in UG659. For R561H, we observed no impact on fitness in 3D7 parasites, although in F32 this mutation exerted a fitness defect similar to M579I and C580Y (Figure 3A–D; Figure 3—source data 1). From these data, we calculated the fitness cost, which represents the percent reduction in growth rate per 48 hr generation of a test line compared to its wild-type isogenic comparator. These costs ranged from <1% to 12% per generation across mutations and lines, with the lowest costs observed in 3D7R561H and UG659C580Y parasites, and the greatest costs observed in the K13 edited UG815 lines (Figure 3E). Comparing data across these four African strains revealed that high RSA survival rates were generally accompanied by high fitness costs, and conversely that low fitness costs were associated with low survival rates. Exceptions were the 3D7R561H and UG659C580Y lines that showed moderate resistance with little to no apparent fitness costs (Figure 3F).

K13 mutations cause differential impacts on in vitro growth rates across gene-edited African strains.

(A–D) Percentages of mutant alleles relative to the wild-type allele over time in (A) 3D7, (B) F32, (C) UG659, and (D) UG815 parasite cultures in which K13 mutant clones were co-cultured at 1:1 starting ratios with isogenic K13 wild-type controls over a period of 36 days. Results, shown as means ± SEM, were obtained from two to five independent experiments, each performed in duplicate. Values are provided in Figure 3—source data 1. (E) The percent reduction in growth rate per 48 hr generation, termed the fitness cost, is presented as mean ± SEM for each mutant line relative to its isogenic wild-type comparator. (F) Fitness costs for mutant lines and isogenic wild-type comparators plotted relative to RSA survival values for the same lines.

Figure 3—source data 1

Fitness assay data for K13 edited African parasite lines and controls.

https://cdn.elifesciences.org/articles/66277/elife-66277-fig3-data1-v3.xlsx

The K13 C580Y mutation has swept rapidly across Cambodia, displacing other K13 variants

We next examined the spatio-temporal distribution of K13 alleles in Cambodia, the epicenter of ART resistance in SE Asia. In total, we analyzed K13 propeller domain sequences from 3327 parasite isolates collected from fourteen Cambodian provinces in the western, northern, eastern, and southern regions (Figure 4—figure supplement 1). Samples were collected between 2001 and 2017, except for the southern region where sample collection was initiated in 2010. A total of 1412 samples (42%) were obtained and sequenced during the period from 2015–2017 and have not previously been published. Earlier samples were reported in Ariey et al., 2014; Ménard et al., 2016. In sum, 19 nonsynonymous polymorphisms in K13 were identified across all regions and years. Of these, only three were present in >10 samples: Y493H (n=83), R539T (n=87), and C580Y (n=1915). Each of these mutations was previously shown to confer ART resistance in vitro (Straimer et al., 2015). Rarer mutations included A418V, I543T, P553L, R561H, P574L, and D584V (Figure 4; Figure 4—source data 1).

Figure 4 with 1 supplement see all
The K13 C580Y allele has progressively outcompeted all other alleles in Cambodia.

(A–D) Stacked bar charts representing the percentages of sequenced samples expressing the K13 wild-type allele or individual variants, calculated based on the total number of samples (listed in parentheses) for a given period. Sample collection was segregated into four regions in Cambodia (detailed in Figure 4—figure supplement 1). All K13 mutant samples harbored a single non-synonymous nucleotide polymorphism. Mutations and numbers of Cambodian samples sequenced per region/year, including prior citations as appropriate, are listed in Figure 4—source data 1.

Figure 4—source data 1

Distribution of K13 alleles over time in Cambodia (2001–2017).

https://cdn.elifesciences.org/articles/66277/elife-66277-fig4-data1-v3.xlsx

This analysis revealed a significant proportion of K13 wild-type parasites in the early 2000s, particularly in northern and eastern Cambodia, where 96% of isolates in 2001–2002 were wild type (Figure 4). In western Cambodia, where ART resistance first emerged (Dondorp et al., 2009; Noedl et al., 2009), the wild-type allele percentage in 2001–2002 had already fallen to 56%. This is striking given that delayed parasite clearance following ACT or artesunate treatment was first documented in 2008–2009 (Noedl et al., 2008; Noedl et al., 2009).

In all four regions, the frequency of the wild-type allele declined substantially over time and the diversity of mutant alleles contracted, with nearly all wild-type and non-K13 C580Y mutant parasites being replaced by parasites harboring the C580Y mutation (Figure 4). This effect was particularly pronounced in the western and the southern regions, where the prevalence of C580Y in 2016–2017 effectively attained 100%, increasing from 22% and 58% in the initial sample sets, respectively (Figure 4A,D). In northern and eastern Cambodia, C580Y also outcompeted all other mutant alleles; however, 19–25% of parasites remained K13 wild type in 2016–2017 (Figure 4B,C). These data show rapid dissemination of K13 C580Y across Cambodia.

SE Asian K13 mutations associated with delayed parasite clearance differ substantially in their ability to confer artemisinin resistance in vitro

Given that most K13 polymorphisms present in the field have not previously been characterized in vitro, we selected a set of mutations to test by gene editing, namely E252Q, F446I, P553L, R561H, and P574L. The positions of these residues are highlighted in Supplementary file 3. F446I is the predominant mutation in Myanmar (Imwong et al., 2020). P553L, R561H, and P574L have each been shown to have multiple independent origins throughout SE Asia (Ménard et al., 2016) and were identified at low frequencies in our sequencing study in Cambodia (Figure 4). Lastly, the E252Q mutation was formerly prevalent on the Thai-Myanmar border, and, despite its location upstream of the beta-propeller domain, has been associated with delayed parasite clearance in vivo (Anderson et al., 2017; Cerqueira et al., 2017; WWARN K13 Genotype-Phenotype Study Group, 2019).

Zinc-finger nuclease- or CRISPR/Cas9-based gene-edited lines expressing K13 E252Q, F4461, P553L, R561H, or P574L were generated in Dd2 or Cam3.II lines expressing wild-type K13 (Dd2WT or Cam3.IIWT) and recombinant parasites were cloned. Early ring-stage parasites were then assayed for their ART susceptibility using the RSA. For comparison, we included published Dd2 and Cam3.II lines expressing either K13 C580Y (Dd2C580Y and Cam3.IIC580Y) or R539T (Dd2R539T and the original parental line Cam3.IIR539T) (Straimer et al., 2015), as well as control lines expressing only the guide-specific silent shield mutations (Dd2ctrl and Cam3.IIctrl).

Both the P553L and R561H mutations yielded mean RSA survival rates comparable to C580Y (4.6% or 4.3% RSA survival for Dd2P553L or Dd2R561H, respectively, vs 4.7% for Dd2C580Y; Figure 5A; Figure 5—source data 1). F446I and P574L showed only modest increases in survival relative to the wild-type parental line (2.0% and 2.1% for Dd2F446I and Dd2P574L, respectively, vs 0.6% for Dd2WT). No change in RSA survival relative to Dd2WT was observed for the Dd2E252Q line. The resistant benchmark Dd2R539T showed a mean RSA survival level of 20.0%, consistent with earlier reports of this mutation conferring high-grade ART resistance in vitro (Straimer et al., 2015; Straimer et al., 2017).

Figure 5 with 1 supplement see all
Southeast Asian K13 mutations yield elevated RSA survival and minor impacts on in vitro growth in gene-edited parasite lines.

(A, B) RSA survival rates for Dd2 (Indochina) and Cam3.II (Cambodia) P. falciparum parasites expressing wild-type or mutant K13. Gene-edited parasites were generated using CRISPR/Cas9 or zinc-finger nucleases. Control (ctrl) lines express silent shield mutations at the K13 gRNA cut site. Parental lines are described in Table 1 and Supplementary file 4. Results show the percentages of early ring-stage parasites (0–3 hr post invasion) that survived a 6 hr pulse of 700 nM DHA, relative to DMSO-treated parasites processed in parallel. Percent survival values are shown as means ± SEM (detailed in Figure 5—source data 1). Results were obtained from three to thirteen independent experiments, each performed in duplicate. p Values were determined by unpaired t tests and were calculated for mutant lines relative to the isogenic line expressing wild-type K13. *** p<0.001; **** p<0.0001. (C) Percent reductions in growth rate per 48 hr generation, expressed as fitness costs, for Dd2 mutant lines relative to the Dd2WT line. Fitness costs were determined by co-culturing the Dd2eGFP reporter line with either the Dd2 K13 wild-type parental line (Dd2WT) or gene-edited K13 mutant lines. Co-cultures were maintained for 20 days and percentages of eGFP+ parasites were determined by flow cytometry (see Figure 5—source data 2 and Figure 5—figure supplement 1). Fitness costs were initially calculated relative to the Dd2eGFP reporter line (Figure 5—figure supplement 1) and then normalized to the Dd2WT line. Mean ± SEM values were obtained from three independent experiments, each performed in triplicate. (D) Fitness costs for K13 mutant lines relative to the Dd2WT line plotted against their corresponding RSA survival values.

Figure 5—source data 1

Ring-stage survival (RSA) assay data for K13 edited SE Asian parasites and controls (Dd2 and Cam3.II strains).

https://cdn.elifesciences.org/articles/66277/elife-66277-fig5-data1-v3.xlsx
Figure 5—source data 2

Fitness assay data (percent eGFP+ parasites) for K13 edited Dd2 parasites and parental control.

https://cdn.elifesciences.org/articles/66277/elife-66277-fig5-data2-v3.xlsx

In contrast to Dd2, editing of the F446I, P553L, and P574L mutations into Cambodian Cam3.IIWT parasites did not result in statistically significant increases in survival rates relative to the K13 wild-type line, in part because the background survival rate of Cam3.IIWT was higher than for Dd2WT. All survival rates were <3%, contrasting with the Cam3.IIR539T parental strain that expresses the R539T mutation (20.4% mean survival; Figure 5B; Figure 5—source data 1). The E252Q mutation did not result in elevated RSA survival in the Cam3.II background, a result also observed with Dd2. Nonetheless, ART resistance was apparent upon introduction of the R561H mutation into Cam3.IIWT parasites, with mean survival rates for Cam3.IIR561H exceeding those for the Cam3.IIC580Y line (13.2% vs 10.0%, respectively). No elevated survival was seen in the Cam3.IIctrl line expressing only the silent shield mutations used at the gRNA cut site.

SE Asian K13 mutations do not impart a significant fitness impact on Dd2 parasites

Prior studies with isogenic gene-edited SE Asian lines have shown that certain K13 mutations can exert fitness costs, as demonstrated by reduced intra-erythrocytic asexual blood stage parasite growth (Straimer et al., 2017; Nair et al., 2018). To determine the fitness impact of the K13 mutations described above, we used an eGFP-based parasite competitive growth assay (Ross et al., 2018). For these experiments, Dd2E252Q, Dd2F446I, Dd2P553L, Dd2R561H, or Dd2P574L parasites were co-cultured with a K13 wild-type eGFP-positive (eGFP+) Dd2 reporter line at starting ratios of 1:1. Proportions of eGFP+ parasites were then assessed every two days by flow cytometry. As controls, we also included the Dd2WT, Dd2bsm, and Dd2C580Y lines. These data provided evidence of a minimal growth impact with the F446I, P553L, and C580Y mutations. In contrast, E252Q, R561H, and P574L resulted in greater fitness costs when compared to Dd2WT parasites (Figure 5C; Figure 5—figure supplement 1; Figure 5—source data 2). Both the C580Y and P553L mutations yielded elevated RSA survival and minimal fitness costs in the Dd2 strain, providing optimal traits for dissemination (Figure 5D). We note that all fitness costs in K13 edited Dd2 parasites were considerably lower than those observed in the majority of the K13 edited African lines described above (Figure 3).

Strain-dependent genetic background differences significantly impact RSA survival rates in culture-adapted Thai isolates

Given the earlier abundance of the R561H and E252Q alleles in border regions of Thailand and Myanmar, we next tested the impact of introducing these mutations into five K13 wild-type Thai isolates (Thai1-5). For comparison, we also edited C580Y into several of these same isolates. These studies revealed a major contribution of the parasite genetic background in dictating the level of mutant K13-mediated ART resistance, as exemplified by the C580Y mutation, which yielded mean survival rates ranging from 2.1% to 15.4% in edited parasites. Trends observed for individual mutations were maintained across strains, with the R561H mutation consistently yielding moderate to high in vitro resistance, at or above the level of C580Y. Consistent with our Dd2 results, E252Q edited parasites did not display significant increases in survival rates relative to isogenic K13 wild-type lines (Figure 6A–E; Figure 6—source data 1).

Thai isolates expressing mutant K13 display variable RSA survival rates.

RSA survival rates for (A–E) K13 edited Thai isolates and (F) K13 E252Q unedited Thai lines, shown as means ± SEM (detailed in Figure 6—source data 1). Results were obtained from three to seven independent experiments, each performed in duplicate. p Values were determined by unpaired t tests and were calculated for mutant lines relative to the isogenic line expressing wild-type K13. * p<0.05; ** p<0.01; *** p<0.001.

Figure 6—source data 1

Ring-stage survival (RSA) assay data for K13 edited Thai parasites and controls.

https://cdn.elifesciences.org/articles/66277/elife-66277-fig6-data1-v3.xlsx

We also profiled two unedited culture-adapted Thai isolates (Thai6E252Q and Thai7E252Q), which express the K13 E252Q mutation that occurs upstream of the propeller domain. Notably, both lines exhibited mean RSA survival rates above the 1% threshold for ART sensitivity (2.7% for Thai6E252Q and 5.1% for Thai7E252Q; Figure 6F). These data suggest that additional genetic factors present in these two Thai isolates are required for E252Q to manifest ART resistance.

Mutations in P. falciparum multidrug resistance protein 2 and ferredoxin do not modulate resistance to artemisinin or parasite fitness in vitro

In a prior genome-wide association study of SE Asian parasites, K13-mediated ART resistance was associated with the D193Y and T484I mutations in the ferredoxin (fd) and multidrug resistance protein 2 (mdr2) genes, respectively (Miotto et al., 2015). To directly test the impact of these mutations in parasite resistance and fitness, we applied CRISPR/Cas9 editing (Supplementary file 6) to revert the fd D193Y and mdr2 T484I mutations to their wild-type sequences. These experiments were performed in the Cambodian K13 C580Y strains RF7C580Y and Cam3.IIC580Y. Isogenic RF7C580Y parasites expressing either the mutant or wild-type fd residue at position 193 showed no difference in RSA survival rates, either at 700 nM (averaging ~27%), or across a range of DHA concentrations down to 1.4 nM (Figure 7A,C; Figure 7—figure supplement 1; Figure 7—source data 1). Editing fd D193Y into the recombinant CamWTC580Y line that expresses K13 C580Y (Straimer et al., 2015) also had no impact on RSA survival (mean RSA survival rate of 11% vs 12.5%). Likewise, Cam3.IIC580Y parasites maintained the same rates of in vitro RSA survival (means ~19–22%) irrespective of their mdr2 allele. Silent shield mutations had no impact on RSA survival rates for either fd or mdr2. eGFP-based fitness assays initiated at different starting ratios of the Dd2 eGFP+ reporter line and fd or mdr2 edited RF7C580Y or Cam3.IIC580Y lines revealed no changes in growth rates for mutant lines compared with their wild-type controls (Figure 7B,D; Figure 7—figure supplement 1; Figure 7—source data 2 and 3). These data suggest that the fd D193Y and mdr2 T484I mutations may be markers of ART-resistant founder populations, but themselves do not contribute directly to ART resistance or augment parasite fitness.

Figure 7 with 1 supplement see all
Ferredoxin (fd) and multidrug resistance protein 2 (mdr2) mutations do not impact RSA survival or in vitro growth rates in K13 C580Y parasites.

RSA survival rates for (A) RF7C580Y parasite lines expressing the fd variant D193Y (parent), this variant plus silent shield mutations (edited control), or fd D193 (edited revertant), and (C) Cam3.IIC580Y parasite lines expressing the mdr2 variant T484I (parent), this variant plus silent shield mutations (edited control), or mdr2 T484 (edited revertant). Parental lines are described in Table 1 and Supplementary file 4. Mean ± SEM survival rates were generated from three independent experiments, each performed in duplicate. (B, D) In vitro eGFP-based fitness assays performed with (B) fd and (D) mdr2 RF7C580Y or Cam3.IIC580Y edited lines, respectively. Competitive growth assays were seeded with individual lines plus the Dd2eGFP+ reporter line at starting ratios of 10:1. Results show percentages of eGFP+ parasites over time. Co-cultures were maintained over a period of 24 days (fd edited lines) or 30 days (mdr2 edited lines), and percentages of eGFP+ parasites were determined by flow cytometry. Results were obtained from two to three independent experiments, each performed in triplicate, and are shown as means ± SEM. All values are provided in Figure 7—source data 1, 2, 3.

Figure 7—source data 1

Ring-stage survival (RSA) assay data for fd and mdr2 edited parasites and controls.

https://cdn.elifesciences.org/articles/66277/elife-66277-fig7-data1-v3.xlsx
Figure 7—source data 2

Fitness assay data (percent eGFP+ parasites) for RF7 fd edited parasites and parental control.

https://cdn.elifesciences.org/articles/66277/elife-66277-fig7-data2-v3.xlsx
Figure 7—source data 3

Fitness assay data (percent eGFP+ parasites) for Cam3.II mdr2 edited parasites and parental control.

https://cdn.elifesciences.org/articles/66277/elife-66277-fig7-data3-v3.xlsx

Discussion

Mutant K13-mediated ART resistance has substantially compromised the efficacy of antimalarial treatments across SE Asia, and the relatively high prevalence of the R561H variant that has recently been associated with delayed clearance in Rwanda highlights the risk of ART resistance emerging and spreading in sub-Saharan Africa (Uwimana et al., 2020; Bergmann et al., 2021; Uwimana et al., 2021). Using gene editing and phenotypic analyses, we provide definitive evidence that the K13 R561H, M579I and C580Y mutations can confer in vitro ART resistance in several African strains. In vitro resistance, as defined using the RSA, was comparable between gene-edited K13 R561H 3D7 parasites (originating from or near Rwanda) and C580Y Dd2 and Cam3.II parasites (from SE Asia). Further investigations into edited African 3D7 parasites showed that these mutations also resulted in an ~30% decrease in K13 protein levels, consistent with earlier studies into the mechanistic basis of mutant K13-mediated ART resistance (Birnbaum et al., 2017; Siddiqui et al., 2017; Yang et al., 2019; Gnädig et al., 2020; Mok et al., 2021). We also observed that K13 mutant African strains differed widely in their RSA survival rates. As an example, when introduced into the Tanzanian F32 strain, the C580Y mutation yielded a 0.3% RSA survival rate, contrasting with 11.8% survival in the Ugandan UG815 strain. These data suggest that F32 parasites lack additional genetic determinants that are required for mutant K13 to confer ART resistance. Collectively, our results provide evidence that certain African strains present no major biological obstacle to becoming ART resistant in vitro upon acquiring K13 mutations. Further gene editing experiments are merited to extend these studies to additional African strains, and to incorporate other variants such as C469Y and A675V that are increasing in prevalence in Uganda (Asua et al., 2021).

Our mixed culture competition assays with African parasites revealed substantial fitness costs with the K13 C580Y mutation in three of the four strains tested (UG659 was the exception). The largest growth defect was observed with the edited UG815 C580Y line, which also yielded the highest level of ART resistance in vitro. These data suggest that K13 C580Y may not easily take hold in Africa where, unlike in SE Asia, infections are often highly polyclonal, generating intra-host competition that impacts a strain’s ability to succeed at the population level. In addition, individuals in highly-endemic African settings generally have high levels of acquired immunity, potentially minimizing infection by relatively unfit parasites, and often have asymptomatic infections that go untreated and are thus less subject to selective drug pressure compared with individuals in SE Asia (Eastman and Fidock, 2009). This situation recalls the history of chloroquine use in Africa, where fitness costs caused by mutations in the primary resistance determinant PfCRT resulted in the rapid resurgence of wild-type parasites following the implementation of other first-line antimalarial therapies (Kublin et al., 2003; Laufer et al., 2006; Ord et al., 2007; Frosch et al., 2014).

An even greater fitness cost was observed with the M579I mutation, earlier detected in an infection acquired in Equatorial Guinea with evidence of in vivo ART resistance (Lu et al., 2017), but which was notably absent in all 3257 African samples reported herein. In contrast, we observed no evident fitness cost in 3D7 parasites expressing the R561H variant, which might help contribute to its increasing prevalence in Rwanda. While our Rwandan samples from 2012 to 2015 observed this mutation at 2% prevalence, samples collected by others in 2018 and 2019 identified this mutation at 12–13% prevalence (Bergmann et al., 2021; Uwimana et al., 2021). One of these reports included evidence associating R561H with delayed parasite clearance in patients treated with the ACT artemether-lumefantrine (Uwimana et al., 2021). These recent data heighten the concern that mutant K13 might be taking hold in certain areas in Africa where it may begin to compromise ACT efficacy.

In Cambodia, our spatio-temporal analysis of K13 sequence diversity highlights the initial emergence of C580Y in the western provinces, and its progressive replacement of other variants in the country. Interestingly, this mutation was already at high prevalence in western Cambodia several years before the first published reports of delayed parasite clearance in ART-treated patients (Noedl et al., 2008; Dondorp et al., 2009; Ariey et al., 2014). The success of this mutation in Cambodia, and elsewhere in the eastern Greater Mekong subregion (Imwong et al., 2020), cannot be explained by resistance alone, as we previously reported that the less common R539T and I543T variants conferred greater ART resistance in vitro (Straimer et al., 2015). Similarly, we now report that the R561H and P553L mutations yield ART resistance at levels comparable to C580Y in Dd2 parasites. In contrast, low-level resistance was observed with F446I, which has nonetheless spread across Myanmar (Imwong et al., 2020). In a separate recent gene editing study, F446I yielded no significant in vitro resistance in 3D7 parasites and was fitness neutral (Siddiqui et al., 2020), consistent with our findings for this mutation in edited Dd2 parasites.

Our studies into the impact of K13 mutations on in vitro growth in Asian Dd2 parasites provide evidence that the C580Y mutation generally exerts less of a fitness cost relative to other K13 variants, as measured in K13 edited parasites co-cultured with an eGFP+ reporter line. A notable exception was P553L, which compared with C580Y was similarly fitness neutral and showed similar RSA values. P553L has nonetheless proven far less successful in its regional dissemination compared with C580Y (Ménard et al., 2016). These data suggest that additional factors have contributed to the success of C580Y in sweeping across SE Asia. These might include specific genetic backgrounds that have favored the dissemination of C580Y parasites, possibly resulting in enhanced transmission potential (Witmer et al., 2020), or ACT use that favored the selection of partner drug resistance in these parasite backgrounds (van der Pluijm et al., 2019). In terms of growth rates in our isogenic Dd2 lines, the most detrimental impacts were observed with E252Q and R561H, which earlier predominated near the Thailand-Myanmar border region, but were later overtaken by C580Y (Phyo et al., 2016). In our study, C580Y produced an optimal combination of no measurable fitness cost and relatively high RSA survival rates in Dd2 parasites. In a prior independent study, however, R561H showed slightly improved fitness relative to C580Y in paired isogenic parasites from Thailand (generated in the NHP4302 strain), providing further evidence that both fitness and resistance are strain-dependent (Nair et al., 2018).

Further research is required to define secondary genetic determinants that could augment mutant K13-mediated ART resistance, and to explore other potential mediators of resistance. Proposed candidates have included fd, mdr2, ap-2μ, ubp1, and pfcoronin, which have earlier been associated with P. falciparum ART susceptibility (Demas et al., 2018; Henrici et al., 2019; Sutherland et al., 2021). Our data argue against a direct role for mutations in fd and mdr2 in the strains tested herein. We also observed no evident association between the genotypes of pfcrt, mdr1, arps10, ap-2μ, or ubp1 and the degree to which mutant K13 conferred ART resistance in vitro in our set of African or Asian strains (Supplementary file 4). Mutations associated with enhanced DNA repair mechanisms have also been observed in ART-resistant SE Asian parasites, supporting the idea that mutant K13 parasites may have an improved ability to repair ART-mediated DNA damage (Xiong et al., 2020). Additional studies are merited to investigate whether these DNA repair mutations may provide a favorable background for the development of ART resistance.

At the population level, we note that P. falciparum genomic structures in Africa tend to be far more diverse than in the epicenter of resistance in Cambodia, where parasite strains are highly sub-structured into a few lineages that can readily maintain complex genetic traits (Amato et al., 2018). A requirement to transmit mutant K13 and additional determinants of resistance in African malaria-endemic settings, where genetic outcrossing is the norm, would predict that ART resistance will spread more gradually on this continent than in SE Asia. It is nonetheless possible that secondary determinants will allow some African strains to offset fitness costs associated with mutant K13, or otherwise augment K13-mediated ART resistance. Identifying such determinants could be possible using genome-wide association studies or genetic crosses between ART-resistant and sensitive African parasites in the human liver-chimeric mouse model of P. falciparum infection (Vaughan et al., 2015; Amambua-Ngwa et al., 2019). Reduced transmission rates in areas of Africa where malaria is declining, leading to lower levels of immunity, may also benefit the emergence and dissemination of mutant K13 (Conrad and Rosenthal, 2019).

Another impediment to the dissemination of ART resistance in Africa is the continued potent activity of lumefantrine, the partner drug in the first line treatment artemether-lumefantrine (Conrad and Rosenthal, 2019). This situation contrasts with SE Asia where ART-resistant parasites have also developed high-level resistance to the partner drug PPQ, with widespread treatment failures enabling the dissemination of multidrug-resistant strains (Conrad and Rosenthal, 2019; van der Pluijm et al., 2019). While the genotyping data presented herein and other recent molecular surveillance studies reveal a low prevalence of mutant K13 in Africa (Kayiba et al., 2021; Schmedes et al., 2021), the emergence and spread of the R561H variant in Rwanda is cause for significant concern. These data call for continuous continent-wide monitoring of the emergence and spread of mutant K13 in Africa, and for studies into whether its emergence in Rwanda is a harbinger of subsequent partner drug resistance and ACT treatment failure.

Materials and methods

Key resources table
Reagent type (species) or resourceDesignationSource or referenceIdentifiersAdditional information
Gene (Plasmodium falciparum 3D7 strain)Kelch13 (K13)PlasmoDBPF3D7_1343700
Gene (Plasmodium falciparum 3D7 strain)Ferredoxin (fd)PlasmoDBPF3D7_1318100
Gene (Plasmodium falciparum 3D7 strain)Multidrug resistance protein 2 (mdr2)PlasmoDBPF3D7_1447900
Strain, strain background
(Plasmodium falciparum)
3D7 clone A10 (3D7WT)D. Goldberg, Washington University School of Medicine, St. Louis, MO, USAsee Table 1 and Supplementary file 4 for additional details on all P. falciparum strains employed herein
Strain, strain background
(Plasmodium falciparum)
F32-TEM (F32WT)F. Benoit-Vical, Université de Toulouse, Toulouse, France
Ariey et al., 2014
Strain, strain background
(Plasmodium falciparum)
UG659 (UG659WT)P. Rosenthal, University of California, San Francisco, CA, USA
Strain, strain background
(Plasmodium falciparum)
UG815 (UG815WT)P. Rosenthal, University of California, San Francisco, CA, USA
Strain, strain background
(Plasmodium falciparum)
Dd2 (Dd2WT)The Malaria Research and Reference Reagent Resource Center (MR4), BEI ResourcesMRA-156
Strain, strain background
(Plasmodium falciparum)
Cam3.II (Cam3.IIR539T)R. Fairhurst, NIAID, NIH, Bethesda, MD, USA
Straimer et al., 2015
PH0306-C
Strain, strain background
(Plasmodium falciparum)
CamWTR. Fairhurst, NIAID, NIH, Bethesda, MD, USA
Straimer et al., 2015
PH0164-C
Strain, strain background
(Plasmodium falciparum)
RF7
(RF7C580Y)
R. Fairhurst, NIAID, NIH, Bethesda, MD, USA
Ross et al., 2018
PH1008-C
Strain, strain background
(Plasmodium falciparum)
Thai1WTT. Anderson, Texas Biomedical Research Institute, San Antonio, TX, USATA32A2A4
Strain, strain background
(Plasmodium falciparum)
Thai2WTT. Anderson, Texas Biomedical Research Institute, San Antonio, TX, USATA50A2B2
Strain, strain background
(Plasmodium falciparum)
Thai3WTT. Anderson, Texas Biomedical Research Institute, San Antonio, TX, USATA85R1
Strain, strain background
(Plasmodium falciparum)
Thai4WTT. Anderson, Texas Biomedical Research Institute, San Antonio, TX, USATA86A3
Strain, strain background
(Plasmodium falciparum)
Thai5WTT. Anderson, Texas Biomedical Research Institute, San Antonio, TX, USANHP-01334-6B
Strain, strain background
(Plasmodium falciparum)
Thai6E252QT. Anderson, Texas Biomedical Research Institute, San Antonio, TX, USANHP4076
Strain, strain background
(Plasmodium falciparum)
Thai7E252QT. Anderson, Texas Biomedical Research Institute, San Antonio, TX, USANHP4673
Strain, strain background
(Escherichia coli)
HST08TakaraCat. #636766Stellar Competent Cells
Genetic reagent
(Plasmodium falciparum)
Transgenic parasite linesThis study and Straimer et al., 2015See Supplementary file 5Available from D. Fidock upon request
Commercial assay or kitIn-Fusion HD Cloning Plus kitTakaraCat. #638909
Commercial assay or kitQuantiFast Multiplex PCR KitQiagenCat. #204654
Sequence-based reagentsOligonucleotidesThis studySee Supplementary file 7
Recombinant DNA reagentsPlasmidsThis studySee Supplementary file 8Available from D. Fidock upon request
Sequence-based reagentsqPCR primers and probesThis studySee Supplementary file 9
AntibodyAnti-K13 (P. falciparum)
(Mouse monoclonal)
I. Trakht, Columbia University Medical Center, New York, NY, USA
Gnädig et al., 2020
Antibody clone E9
WB (1:1000)
AntibodyAnti-ERD2 (P. falciparum)
(Rabbit polyclonal)
MR4, BEI ResourcesMRA-1WB (1:1000)
AntibodyStarBright Blue 700 goat anti-mouseBio-Rad12004158WB (1:200)
AntibodyStarBright Blue 520 goat anti-rabbitBio-Rad12005869WB (1:1000)
Other4–20% Criterion TGX Precast Protein GelBio-Rad5671093Used with recommended buffers, also purchased from Bio-Rad
Chemical compound, drugCarbenicillin disodium saltSigmaC1389
Chemical compound, drugWR99210Jacobus Pharmaceuticals
Chemical compound, drugDihydroartemisinin (DHA)SigmaD7439
Software, algorithmGraphPad Prism Version 9GraphPad Software, San Diego, CA, USAgraphpad.com
Software, algorithmImageJ softwareNIH, Bethesda, MD, USAimagej.nih.gov

Sample collection and K13 genotyping

Request a detailed protocol

Samples were obtained as blood-spot filter papers from patients seeking treatment at sites involved in national surveys of antimalarial drug resistance, or patients enrolled in therapeutic efficacy studies, or asymptomatic participants enrolled in surveillance programs. Collection details for African and Cambodian samples are provided in Figure 1—source data 1 and Figure 4—source data 1, respectively. Samples were processed at the Pasteur Institute in Paris or the Pasteur Institute in Cambodia, as detailed in Supplementary file 1. These investigators vouch for the accuracy and completeness of the molecular data. DNA was extracted from dried blood spots using QIAmp Mini kits, as described previously (Ménard et al., 2016). A nested PCR was performed on each sample to amplify the K13 propeller domain sequence, corresponding to codons 440–680. PCR products were sequenced using internal primers and electropherograms were analyzed on both strands using the Pf3D7_1343700 3D7 sequence as the wild-type reference. Quality controls included adding six blinded quality-control samples to each 96-well sequencing plate prepared from samples from each in-country partner, and independently retesting randomly selected blood samples. Isolates with mixed alleles were considered to be mutant for the purposes of estimating mutation frequencies.

P. falciparum parasite in vitro culture

Request a detailed protocol

P. falciparum asexual blood-stage parasites were cultured in human erythrocytes at 3% hematocrit in RPMI-1640 medium supplemented with 2 mM L-glutamine, 50 mg/L hypoxanthine, 25 mM HEPES, 0.21% NaHCO3, 10 mg/L gentamycin, and 0.5% w/v Albumax II (Invitrogen). Parasites were maintained at 37°C in 5% O2, 5% CO2, and 90% N2. The geographic origins and years of culture adaptation for lines employed herein are described in Supplementary file 4. Parasite lines were authenticated by genotyping resistance genes and were screened by PCR for Mycoplasma every 3–6 months.

Whole-genome sequencing of parental lines

Request a detailed protocol

To define the genome sequences of P. falciparum lines used for transfection, we lysed parasites in 0.05% saponin, washed them with 1×PBS, and purified genomic DNA (gDNA) using the QIAamp DNA Blood Midi Kit (Qiagen). DNA concentrations were quantified by NanoDrop (Thermo Scientific) and Qubit (Invitrogen) prior to sequencing. 200 ng of gDNA was used to prepare sequencing libraries using the Illumina Nextera DNA Flex library prep kit with dual indices. Samples were multiplexed and sequenced on an Illumina MiSeq to obtain 300 bp paired-end reads at an average of 50× depth of coverage. Sequence reads were aligned to the P. falciparum 3D7 reference genome (PlasmoDB version 36) using Burrow-Wheeler Alignment. PCR duplicates and unmapped reads were filtered out using Samtools and Picard. Reads were realigned around indels using GATK RealignerTargetCreator, and base quality scores were recalibrated using GATK BaseRecalibrator. GATK HaplotypeCaller (version 3.8) was used to identify all single nucleotide polymorphisms (SNPs). SNPs were filtered based on quality scores (variant quality as function of depth QD >1.5, mapping quality >40, min base quality score >18) and read depth (>5) to obtain high-quality SNPs, which were annotated using snpEFF. Integrated Genome Viewer was used to visually verify the presence of SNPs. BIC-Seq was used to check for copy number variations using the Bayesian statistical model (Xi et al., 2011). Copy number variations in highly polymorphic surface antigens and multi-gene families were removed as these are prone to stochastic changes during in vitro culture.

Whole-genome sequencing data were used to determine the genotypes of the antimalarial drug resistance loci pfcrt, mdr1, dhfr, and dhps (Haldar et al., 2018). We also genotyped fd, arps10, mdr2, ubp1, and ap-2μ, which were previously associated with ART resistance (Henriques et al., 2014; Miotto et al., 2015; Cerqueira et al., 2017; Adams et al., 2018). These results are described in Supplementary file 4.

Cloning of K13, fd, and mdr2 plasmids

Request a detailed protocol

Zinc-finger nuclease-meditated editing of select mutations in the K13 locus was performed as previously described (Straimer et al., 2015). CRISPR/Cas9 editing of K13 mutations was achieved using the pDC2-cam-coSpCas9-U6-gRNA-hdhfr all-in-one plasmid that contains a P. falciparum codon-optimized Cas9 sequence, a human dihydrofolate reductase (hdhfr) gene expression cassette (conferring resistance to WR99210) and restriction enzyme insertion sites for the guide RNA (gRNA) and donor template (White et al., 2019). A K13 propeller domain-specific gRNA was introduced into this vector at the BbsI restriction sites using the primer pair p1+p2 (Supplementary file 7) using T4 DNA ligase (New England BioLabs). Oligos were phosphorylated and annealed prior to cloning. A donor template consisting of a 1.5 kb region of the K13 coding region including the entire propeller domain was amplified using the primer pair p3+p4 and cloned into the pGEM T-easy vector system (Promega). This donor sequence was subjected to site-directed mutagenesis in the pGEM vector to introduce silent shield mutations at the Cas9 cleavage site using the primer pair p5+p6, and to introduce allele-specific mutations using primer pairs (p7 to p20). K13 donor sequences were amplified from the pGEM vector using the primer pair p21+p22 and sub-cloned into the pDC2-cam-coSpCas9-U6-gRNA-hdhfr plasmid at the EcoRI and AatII restriction sites by In-Fusion Cloning (Takara). Final plasmids (see Supplementary file 8) were sequenced using primers p23 to p25. A schematic showing the method of K13 plasmid construction can be found in Supplementary file 2. Both our customized zinc-finger nuclease and CRISPR/Cas9 approaches generated the desired amino acid substitutions without the genomic integration of any plasmid sequences or any additional amino acid changes in the K13 locus, and thus provide fully comparable data.

CRISPR/Cas9 editing of fd and mdr2 was performed using a separate all-in-one plasmid, pDC2-cam-Cas9-U6-gRNA-hdhfr, generated prior to the development of the codon-optimized version used above for K13 (Lim et al., 2016). Cloning was performed as for K13, except for gRNA cloning that was performed using In-Fusion cloning (Takara) rather than T4 ligase. Cloning of gRNAs was performed using primer pair p29+p30 for fd and p42+p43 for mdr2. Donor templates were amplified and cloned into the final vector using the primer pairs p31+p32 for fd and p44+p45 for mdr2. Site-directed mutagenesis was performed using the allele-specific primer pairs p33+p34 or p35+p36 for fd, and p46+p47 or p48+p49 for mdr2. All final plasmids (for both fd and mdr2, see Supplementary file 8) were sequenced using the primer pair p37+p38 (Supplementary file 7). Schematic representations of final plasmids are shown in Supplementary file 6.

Generation of K13, fd, and mdr2 gene-edited parasite lines

Request a detailed protocol

Gene-edited lines were generated by electroporating ring-stage parasites at 5–10% parasitemia with 50 μg of purified circular plasmid DNA resuspended in Cytomix. Transfected parasites were selected by culturing in the presence of WR99210 (Jacobus Pharmaceuticals) for six days post electroporation. Parental lines harboring 2–3 mutations in the P. falciparum dhfr gene were exposed to 2.5 nM WR99210, while parasites harboring four dhfr mutations were selected under 10 nM WR99210 (see Supplementary file 4 for dhfr genotypes of transfected lines). Parasite cultures were monitored for recrudescence by microscopy for up six weeks post electroporation. To test for successful editing, the K13 locus was amplified directly from parasitized whole blood using the primer pair p26+p27 (Supplementary file 7) and the MyTaq Blood-PCR Kit (Bioline). Primer pairs p39+p40 and p50+p51 were used to amplify fd and mdr2, respectively. PCR products were submitted for Sanger sequencing using the PCR primers as well as primer p28 in the case of K13, p41 (fd) or p52 (mdr2). Bulk-transfected cultures showing evidence of editing by Sanger sequencing were cloned by limiting dilution. All gene-edited transgenic lines generated herein are described in Supplementary file 5.

Parasite synchronization, ring-stage survival assays (RSAs), and flow cytometry

Request a detailed protocol

Synchronized parasite cultures were obtained by exposing predominantly ring-stage cultures to 5% D-Sorbitol (Sigma) for 15 min at 37°C to remove mature parasites. After 36 hr of subsequent culture, multinucleated schizonts were purified over a density gradient consisting of 75% Percoll (Sigma). Purified schizonts were incubated with fresh RBCs for 3 hr, and early rings (0–3 hr post invasion; hpi) were treated with 5% D-Sorbitol to remove remaining schizonts.

In vitro RSAs were conducted as previously described, with minor adaptations (Straimer et al., 2015). Briefly, tightly synchronized 0–3 hpi rings were exposed to a pharmacologically-relevant dose of 700 nM DHA or 0.1% dimethyl sulfoxide (DMSO; vehicle control) for 6 hr at 1% parasitemia and 2% hematocrit, washed three times with RPMI medium to remove drug, transferred to fresh 96-well plates, and cultured for an additional 66 hr in drug-free medium. Removal of media and resuspension of parasite cultures was performed on a Freedom Evo 100 liquid-handling instrument (Tecan). Parasitemias were measured at 72 hr by flow cytometry (see below) with at least 100,000 events captured per sample. Parasite survival was expressed as the percentage value of the parasitemia in DHA-treated samples divided by the parasitemia in DMSO-treated samples processed in parallel. We considered any RSA mean survival rates <2% to be ART sensitive.

Flow cytometry was performed on an BD Accuri C6 Plus cytometer with a HyperCyt plate sampling attachment (IntelliCyt), or on an iQue Screener Plus cytometer (Sartorius). Cells were stained with 1×SYBR Green (ThermoFisher) and 100 nM MitoTracker DeepRed (ThermoFisher) for 30 min and diluted in 1×PBS prior to sampling. Percent parasitemia was determined as the percentage of MitoTracker-positive and SYBR Green-positive cells.

Western blot analysis of K13 expression levels in edited lines

Request a detailed protocol

Western blots were performed with lysates from tightly synchronized rings harvested 0–6 hr post invasion. Parasite cultures were washed twice in ice-cold 1× phosphate-buffered saline (PBS), and parasites were isolated by treatment with 0.05% saponin in PBS. Released parasites were lysed in 4% SDS, 0.5% Triton X-100 and 0.5% PBS supplemented with 1× protease inhibitors (Halt Protease Inhibitors Cocktail, ThermoFisher). Samples were centrifuged at 14,000 rpm for 10 min to pellet cellular debris. Supernatants were collected and protein concentrations were determined using the DC protein assay kit (Bio-Rad). Laemmli Sample Buffer (Bio-Rad) was added to lysates and samples were denatured at 90°C for 10 min. Proteins were electrophoresed on precast 4–20% Tris-Glycine gels (Bio-Rad) and transferred onto nitrocellulose membranes. Western blots were probed with a 1:1000 dilution of primary antibodies to K13 (Gnädig et al., 2020) or the loading control ERD2 (BEI Resources), followed by a 1:200 dilution of fluorescent StarBright secondary antibodies (Bio-Rad). Western blots were imaged on a ChemiDoc system (Bio-Rad) and band intensities quantified using ImageJ.

TaqMan allelic discrimination real-time (quantitative) PCR-based fitness assays

Request a detailed protocol

Fitness assays with African K13 edited parasite lines were performed by co-culturing isogenic wild-type unedited and mutant edited parasites in 1:1 ratios. Assays were initiated with tightly synchronized trophozoites. Final culture volumes were 3 mL. Cultures were maintained in 12-well plates and monitored every four days over a period of 36 days (18 generations) by harvesting at each time point a fraction of each co-culture for saponin lysis. gDNA was then extracted using the QIAamp DNA Blood Mini Kit (Qiagen). The percentage of the wild-type or mutant allele in each sample was determined in TaqMan allelic discrimination real-time PCR assays. TaqMan primers (forward and reverse) and TaqMan fluorescence-labeled minor groove binder probes (FAM or HEX, Eurofins) are described in Supplementary file 9. Probes were designed to specifically detect the K13 R561H, M579I, or C580Y propeller mutations. The efficiency and sensitivity of the TaqMan primers was assessed using standard curves comprising 10-fold serially diluted templates ranging from 10 ng to 0.001 ng. Robustness was demonstrated by high efficiency (88–95%) and R2 values (0.98–1.00). The quantitative accuracy in genotype calling was assessed by performing multiplex qPCR assays using mixtures of wild-type and mutant plasmids in fixed ratios (0:100, 20:80, 40:60, 50:50, 60:40, 80:20, 100:0). Triplicate data points clustered tightly, indicating high reproducibility in the data across the fitted curve (R2 = 0.89–0.91).

Purified gDNA from fitness co-cultures was subsequently amplified and labeled using the primers and probes described in Supplementary file 9. qPCR reactions for each sample were run in triplicate. 20 μL reactions consisted of 1×QuantiFAST reaction mix containing ROX reference dye (Qiagen), 0.66 µM forward and reverse primers, 0.16 µM FAM-MGB and HEX-MGB TaqMan probes, and 10 ng genomic DNA. Amplification and detection of fluorescence were carried out on a QuantStudio3 qPCR machine (Applied Biosystems) using the genotyping assay mode. Cycling conditions were as follows: 30 s at 60°C; 5 min at 95°C; and 40 cycles of 30 s at 95°C and 1 min at 60°C for primer annealing and extension. Every assay was run with positive controls (wild-type or mutant plasmids at different fixed ratios). No-template negative controls (water) in triplicates were processed in parallel. Rn, the fluorescence of the FAM or HEX probe, was normalized to the fluorescence signal of the ROX reporter dye. Background-normalized fluorescence (Rn minus baseline, or ΔRn) was calculated as a function of cycle number.

To determine the wild-type or mutant allele frequency in each sample, we first confirmed the presence of the allele by only retaining values where the threshold cycle (Ct) of the sample was less than the no-template control by at least three cycles. Next, we subtracted the ΔRn of the samples from the background ΔRn of the no-template negative control. We subsequently normalized the fluorescence to 100% using the positive control plasmids to obtain the percentage of the wild-type and mutant alleles for each sample. The final percentage of the mutant allele was defined as the average of these two values: the normalized percentage of the mutant allele, and 100% minus the normalized percentage of the wild-type allele.

eGFP-based fitness assays

Request a detailed protocol

Fitness assays with Dd2, RF7C580Y, and Cam3.IIC580Y K13, fd, or mdr2 edited parasite lines were performed using mixed culture competition assays with an eGFP+ Dd2 reporter line (Ross et al., 2018). This reporter line uses a calmodulin (cam) promoter sequence to express high levels of GFP and includes hdhfr and blasticidin S-deaminase expression cassettes. This line was earlier reported to have a reduced growth rate relative to parental non-recombinant Dd2, presumably at least in part because of its high levels of GFP expression (Ross et al., 2018; Dhingra et al., 2019). With our Dd2 parasites, K13 edited lines were co-cultured in 1:1 ratios with the reporter line. This ratio was adjusted to 10:1 or 100:1 for fd edited RF7C580Y and mdr2 edited Cam3.IIC580Y parasites relative to the eGFP line, given the slower rate of growth with RF7 and Cam3.II parasites. Fitness assays were initiated with tightly synchronized trophozoites in 96-well plates with 200 μL culture volumes. Percentages of eGFP+ parasites were monitored by flow cytometry every two days over a period of 20 days (10 generations). Flow cytometry was performed as written above, except that only 100 nM MitoTracker DeepRed staining was used to detect total parasitemias, since SYBR Green and eGFP fluoresce in the same channel.

Fitness costs

Request a detailed protocol

The fitness cost associated with a line expressing a given K13 mutation was calculated relative to its isogenic wild-type counterpart using the following equation:

P'=P((1-x)n)

where P’ is equal to the parasitemia at the assay endpoint, P is equal to the parasitemia on day 0, n is equal to the number of generations from the assay start to finish, and x is equal to the fitness cost. This equation assumes 100% growth for the wild-type comparator line. For qPCR and GFP-based fitness assays, days 36 and 20 were set as the assay endpoints, resulting in the number of parasite generations (n) being set to 18 and 10, respectively.

Data availability

All data generated or analysed during this study are included in the manuscript and supporting files. Source data files have been provided for Figures 1-7.

References

Article and author information

Author details

  1. Barbara H Stokes

    Department of Microbiology and Immunology, Columbia University Irving Medical Center, New York, United States
    Contribution
    Conceptualization, Data curation, Formal analysis, Validation, Investigation, Visualization, Methodology, Writing - original draft, Writing - review and editing
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-9519-487X
  2. Satish K Dhingra

    Department of Microbiology and Immunology, Columbia University Irving Medical Center, New York, United States
    Contribution
    Conceptualization, Formal analysis, Investigation, Methodology
    Competing interests
    No competing interests declared
  3. Kelly Rubiano

    Department of Microbiology and Immunology, Columbia University Irving Medical Center, New York, United States
    Contribution
    Data curation, Investigation
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-4149-4030
  4. Sachel Mok

    Department of Microbiology and Immunology, Columbia University Irving Medical Center, New York, United States
    Contribution
    Formal analysis, Investigation, Methodology
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-9605-0154
  5. Judith Straimer

    Department of Microbiology and Immunology, Columbia University Irving Medical Center, New York, United States
    Contribution
    Conceptualization, Investigation
    Competing interests
    No competing interests declared
  6. Nina F Gnädig

    Department of Microbiology and Immunology, Columbia University Irving Medical Center, New York, United States
    Contribution
    Investigation
    Competing interests
    No competing interests declared
  7. Ioanna Deni

    Department of Microbiology and Immunology, Columbia University Irving Medical Center, New York, United States
    Contribution
    Investigation
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-5266-8243
  8. Kyra A Schindler

    Department of Microbiology and Immunology, Columbia University Irving Medical Center, New York, United States
    Contribution
    Investigation
    Competing interests
    No competing interests declared
  9. Jade R Bath

    Department of Microbiology and Immunology, Columbia University Irving Medical Center, New York, United States
    Contribution
    Investigation
    Competing interests
    No competing interests declared
  10. Kurt E Ward

    1. Department of Microbiology and Immunology, Columbia University Irving Medical Center, New York, United States
    2. Department of Microbiology and Immunology, University of Otago, Dunedin, New Zealand
    Contribution
    Investigation
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-3614-436X
  11. Josefine Striepen

    Department of Microbiology and Immunology, Columbia University Irving Medical Center, New York, United States
    Contribution
    Investigation
    Competing interests
    No competing interests declared
  12. Tomas Yeo

    Department of Microbiology and Immunology, Columbia University Irving Medical Center, New York, United States
    Contribution
    Data curation, Software, Formal analysis, Methodology
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-2923-6060
  13. Leila S Ross

    Department of Microbiology and Immunology, Columbia University Irving Medical Center, New York, United States
    Contribution
    Investigation
    Competing interests
    No competing interests declared
  14. Eric Legrand

    Malaria Genetics and Resistance Unit, Institut Pasteur, INSERM U1201, CNRS ERL9195, Paris, France
    Contribution
    Investigation
    Competing interests
    No competing interests declared
  15. Frédéric Ariey

    Institut Cochin, INSERM U1016, Université Paris Descartes, Paris, France
    Contribution
    Software, Methodology
    Competing interests
    No competing interests declared
  16. Clark H Cunningham

    Department of Genetics, University of North Carolina at Chapel Hill, Chapel Hill, United States
    Contribution
    Investigation
    Competing interests
    No competing interests declared
  17. Issa M Souleymane

    Programme National de Lutte Contre le Paludisme au Tchad, Ndjamena, Chad
    Contribution
    Resources, Investigation
    Competing interests
    No competing interests declared
  18. Adama Gansané

    Centre National de Recherche et de Formation sur le Paludisme, Ouagadougou, Burkina Faso
    Contribution
    Resources, Investigation
    Competing interests
    No competing interests declared
  19. Romaric Nzoumbou-Boko

    Laboratoire de Parasitologie, Institut Pasteur de Bangui, Bangui, Central African Republic
    Contribution
    Resources, Investigation
    Competing interests
    No competing interests declared
  20. Claudette Ndayikunda

    University Teaching Hospital of Kamenge, Bujumbura, Burundi
    Contribution
    Resources, Investigation
    Competing interests
    No competing interests declared
  21. Abdunoor M Kabanywanyi

    Ifakara Health Institute, Dar es Salaam, United Republic of Tanzania
    Contribution
    Resources, Investigation
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-2980-0429
  22. Aline Uwimana

    Malaria and Other Parasitic Diseases Division, Rwanda Biomedical Centre, Kigali, Rwanda
    Contribution
    Resources, Investigation
    Competing interests
    No competing interests declared
  23. Samuel J Smith

    National Malaria Control Program, Freetown, Sierra Leone
    Contribution
    Resources, Investigation
    Competing interests
    No competing interests declared
  24. Olimatou Kolley

    National Malaria Control Program, Banjul, Gambia
    Contribution
    Resources, Investigation
    Competing interests
    No competing interests declared
  25. Mathieu Ndounga

    Programme National de Lutte Contre le Paludisme, Brazzaville, Democratic Republic of the Congo
    Contribution
    Resources, Investigation
    Competing interests
    No competing interests declared
  26. Marian Warsame

    School of Public Health and Community Medicine, University of Gothenburg, Gothenburg, Sweden
    Contribution
    Resources, Investigation
    Competing interests
    No competing interests declared
  27. Rithea Leang

    National Center for Parasitology, Entomology & Malaria Control, Phnom Penh, Cambodia
    Contribution
    Formal analysis, Investigation
    Competing interests
    No competing interests declared
  28. François Nosten

    1. Shoklo Malaria Research Unit, Mahidol-Oxford Tropical Medicine Research Unit, Faculty of Tropical Medicine, Mahidol University, Mae Sot, Thailand
    2. Centre for Tropical Medicine and Global Health, Nuffield Department of Medicine, University of Oxford, Oxford, United Kingdom
    Contribution
    Resources
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-7951-0745
  29. Timothy JC Anderson

    Texas Biomedical Research Institute, San Antonio, United States
    Contribution
    Resources
    Competing interests
    No competing interests declared
  30. Philip J Rosenthal

    Department of Medicine, University of California, San Francisco, San Francisco, United States
    Contribution
    Resources, Writing - review and editing
    Competing interests
    No competing interests declared
  31. Didier Ménard

    Malaria Genetics and Resistance Unit, Institut Pasteur, INSERM U1201, CNRS ERL9195, Paris, France
    Contribution
    Conceptualization, Resources, Supervision, Funding acquisition, Writing - review and editing
    Competing interests
    No competing interests declared
  32. David A Fidock

    Division of Infectious Diseases, Department of Medicine, Columbia University Irving Medical Center, New York, United States
    Contribution
    Conceptualization, Supervision, Funding acquisition, Visualization, Writing - original draft, Project administration, Writing - review and editing
    For correspondence
    df2260@cumc.columbia.edu
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-6753-8938

Funding

National Institute of Allergy and Infectious Diseases (R01 AI109023)

  • David A Fidock

U.S. Department of Defense (W81XWH1910086)

  • David A Fidock

Bill and Melinda Gates Foundation (OPP1201387)

  • David A Fidock

Wellcome Trust (106698)

  • François Nosten

National Institute of Allergy and Infectious Diseases (R37 AI048071)

  • Timothy Anderson

National Institute of Allergy and Infectious Diseases (T32 AI106711)

  • David A Fidock

World Health Organization

  • Didier Menard
  • David A Fidock

National Institute of Allergy and Infectious Diseases (R01 AI075045)

  • Philip J Rosenthal

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

We thank Dr. Pascal Ringwald (World Health Organization) for his support and feedback. DAF gratefully acknowledges the US National Institutes of Health (R01 AI109023), the Department of Defense (W81XWH1910086) and the Bill and Melinda Gates Foundation (OPP1201387) for their financial support. BHS was funded in part by T32 AI106711 (PD: D Fidock). SM is a recipient of a Human Frontiers of Science Program Long-Term Fellowship. CHC was supported in part by the NIH (R01 AI121558; PI: Jonathan Juliano). FN is supported by the Wellcome Trust of Great Britain (Grant ID: 106698). TJCA acknowledges funding support from the NIH (R37 AI048071). PJR received funding support from the NIH (R01 AI075045). DAF and DM gratefully acknowledge the World Health Organization for their funding. We thank the following individuals for their kind help with the K13 genotyped samples – Chad: Ali S Djiddi, Mahamat S I Diar, Kodbessé Boulotigam, Mbanga Djimadoum, Hamit M Alio, Mahamat M H Taisso, Issa A Haggar; Burkina Faso: TES 2017–2018 team and the US President’s Malaria Initiative through the Improving Malaria Care Project as the funding agency for the study in Burkina Faso, Chris-Boris G Panté-Wockama; Burundi: Dismas Baza; Tanzania: Mwaka Kakolwa, Celine Mandara, Tanzania TES coordination team for the Ministry of Health; Sierra Leone: Anitta R Y Kamara, Foday Sahr, Mohamed Samai; The Gambia: Balla Kandeh, Joseph Okebe, Serign J Ceesay, Baboucarr Babou, Emily Jagne, Alsan Jobe; Congo: Brice S Pembet, Jean M Youndouka; Somalia: Jamal Ghilan Hefzullah Amran, Abdillahi Mohamed Hassan, Abdikarim Hussein Hassan and Ali Abdulrahman; Rwanda: extended TES team for the Malaria and Other Parasitic Diseases Division, Rwanda Biomedical Centre.

Ethics

Human subjects: Health care facilities were in charge of collecting anonymized P. falciparum positive cases. Identification of individuals cannot be established. The studies were approved by ethics committees listed in Supplementary file 1. We note that the sponsors had no role in the study design or in the collection or analysis of the data. There was no confidentiality agreement between the sponsors and the investigators.

Copyright

© 2021, Stokes et al.

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 4,139
    views
  • 758
    downloads
  • 100
    citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Barbara H Stokes
  2. Satish K Dhingra
  3. Kelly Rubiano
  4. Sachel Mok
  5. Judith Straimer
  6. Nina F Gnädig
  7. Ioanna Deni
  8. Kyra A Schindler
  9. Jade R Bath
  10. Kurt E Ward
  11. Josefine Striepen
  12. Tomas Yeo
  13. Leila S Ross
  14. Eric Legrand
  15. Frédéric Ariey
  16. Clark H Cunningham
  17. Issa M Souleymane
  18. Adama Gansané
  19. Romaric Nzoumbou-Boko
  20. Claudette Ndayikunda
  21. Abdunoor M Kabanywanyi
  22. Aline Uwimana
  23. Samuel J Smith
  24. Olimatou Kolley
  25. Mathieu Ndounga
  26. Marian Warsame
  27. Rithea Leang
  28. François Nosten
  29. Timothy JC Anderson
  30. Philip J Rosenthal
  31. Didier Ménard
  32. David A Fidock
(2021)
Plasmodium falciparum K13 mutations in Africa and Asia impact artemisinin resistance and parasite fitness
eLife 10:e66277.
https://doi.org/10.7554/eLife.66277

Share this article

https://doi.org/10.7554/eLife.66277

Further reading

    1. Epidemiology and Global Health
    2. Genetics and Genomics
    Tianyu Zhao, Hui Li ... Li Chen
    Research Article

    Alzheimer’s disease (AD) is a complex degenerative disease of the central nervous system, and elucidating its pathogenesis remains challenging. In this study, we used the inverse-variance weighted (IVW) model as the major analysis method to perform hypothesis-free Mendelian randomization (MR) analysis on the data from MRC IEU OpenGWAS (18,097 exposure traits and 16 AD outcome traits), and conducted sensitivity analysis with six models, to assess the robustness of the IVW results, to identify various classes of risk or protective factors for AD, early-onset AD, and late-onset AD. We generated 400,274 data entries in total, among which the major analysis method of the IVW model consists of 73,129 records with 4840 exposure traits, which fall into 10 categories: Disease, Medical laboratory science, Imaging, Anthropometric, Treatment, Molecular trait, Gut microbiota, Past history, Family history, and Lifestyle trait. More importantly, a freely accessed online platform called MRAD (https://gwasmrad.com/mrad/) has been developed using the Shiny package with MR analysis results. Additionally, novel potential AD therapeutic targets (CD33, TBCA, VPS29, GNAI3, PSME1) are identified, among which CD33 was positively associated with the main outcome traits of AD, as well as with both EOAD and LOAD. TBCA and VPS29 were negatively associated with the main outcome traits of AD, as well as with both EOAD and LOAD. GNAI3 and PSME1 were negatively associated with the main outcome traits of AD, as well as with LOAD, but had no significant causal association with EOAD. The findings of our research advance our understanding of the etiology of AD.

    1. Epidemiology and Global Health
    Xiaoning Wang, Jinxiang Zhao ... Dong Liu
    Research Article

    Artificially sweetened beverages containing noncaloric monosaccharides were suggested as healthier alternatives to sugar-sweetened beverages. Nevertheless, the potential detrimental effects of these noncaloric monosaccharides on blood vessel function remain inadequately understood. We have established a zebrafish model that exhibits significant excessive angiogenesis induced by high glucose, resembling the hyperangiogenic characteristics observed in proliferative diabetic retinopathy (PDR). Utilizing this model, we observed that glucose and noncaloric monosaccharides could induce excessive formation of blood vessels, especially intersegmental vessels (ISVs). The excessively branched vessels were observed to be formed by ectopic activation of quiescent endothelial cells (ECs) into tip cells. Single-cell transcriptomic sequencing analysis of the ECs in the embryos exposed to high glucose revealed an augmented ratio of capillary ECs, proliferating ECs, and a series of upregulated proangiogenic genes. Further analysis and experiments validated that reduced foxo1a mediated the excessive angiogenesis induced by monosaccharides via upregulating the expression of marcksl1a. This study has provided new evidence showing the negative effects of noncaloric monosaccharides on the vascular system and the underlying mechanisms.