- Open Access
Observing in real time the evolution of artemisinin resistance in Plasmodium falciparum
BMC Medicine volume 13, Article number: 67 (2015)
Simple genetic changes that correlate with drug resistance are used routinely to identify resistant pathogens. These “molecular markers” have usually been defined long after the phenotype of resistance was noted. The molecular changes at the “end game” reflect a long and complex evolution of genetic changes, but once a solidly resistant set of changes assembles under drug selection, that genotype is likely to become fixed, and resistant pathogens will spread widely.
Artemisinins are currently used worldwide to treat malaria caused by Plasmodium falciparum, but parasite response has diminished rapidly in the Mekong region of Southeast Asia. Should artemisinins lose potency completely and this effect spread worldwide, effective malaria treatment would be almost impossible. The full range of modern methods has been applied to define rapidly the genetic changes responsible. Changes associated with artemisinin resistance are complex and seem to be evolving rapidly, especially in Southeast Asia. This is a rare chance to observe the early stages in evolution of resistance, and to develop strategies to reverse or mitigate the trend and to protect these key medicines.
There may be no better scientific place to study an evolutionary arms race than the interaction between pathogens and the drugs needed to treat them. The clear importance to public health, the relatively rapid time frame of pathogen responses, and the richness of the data available provide both context and motivation for these studies . The interaction of Plasmodium falciparum, the parasite that causes the most dangerous form of malaria, with antimalarial drugs provides an outstanding example. In 2008, it was first noted that P. falciparum in western Cambodia were evolving resistance to the artemisinin drugs, the cornerstone of current global antimalarial policy [2-4]. A paper published recently in Nature Genetics provides new insights into the genetic complexity and population structure of P. falciparum in an area where the resistance is evolving rapidly .
Artemisinins are one principal component of artemisinin combination therapies (ACTs). In the combination, the artemisinin component clears the parasite rapidly from the patient’s blood, but depends upon a long-acting partner drug to dispatch the remaining parasite biomass. One of the earliest indications that the artemisinin drugs are compromised is a delay in the initial clearance of the parasites in the first days after treatment [2,3]. The rate of parasite clearance is a challenging metric to quantify, but further studies were then launched to determine the extent of the resistant phenotype in the Mekong region [6-11] and Africa . A collaborative effort by the research community developed a practical tool to assess the rate of parasite clearance reproducibly between studies , and a laboratory-based assay that correlates with the slow clearance phenotype was developed and validated [14-16]. In 2011 an extensive multicenter trial of artemisinin efficacy in 13 sites in the Greater Mekong region and 2 in Africa was launched, the Tracking Resistance to Artemisinin Collaboration (TRAC) .
Molecular markers of resistance
Previous studies defined molecular markers for resistance to antimalarials long after the resistant parasites were widespread , but this time, modern genomic and analytical tools were available, and identification of a simple molecular marker to track the slow-clearing parasites was an immediate goal. The phenotype was shown to depend largely on the genotype of the parasite , and combining data from the field studies, slow-clearing parasites from Cambodia [10,20] and western Thailand  were shown to have some broad genomic regions in common. In late 2014, a molecular signature in what is called the “propeller region” of the highly conserved Kelch 13 gene (K13, PF3D7_1343700), was shown to correlate strongly with the slow-clearing phenotype in Cambodian parasites . Rapidly, the K13 genotypes of parasites from the TRAC study and earlier field studies in the Southeast Asian region were determined [17,23,24], and other studies are being added rapidly . More than 30 different K13 mutant alleles have been identified; each carries a single nucleotide polymorphism (SNP) that changes one amino acid in the propeller region of the gene. So far in the Mekong region, most of the parasites with these mutant alleles correlate with slow clearance [17,22-24].
A few of these mutant alleles are very common, observed in several locations, but others are rare, identified in a single location or parasite population. When the DNA flanking the coding region of K13 was examined, two very common alleles had spread locally, but these same K13 alleles had also emerged independently on different genetic backgrounds . Miotto and colleagues determined full genome sequences of parasites collected in the TRAC study and compared the genetic signatures of slow- and fast-clearing parasites . As expected, the slow-clearing parasites carried a variety of K13 mutant alleles; analysis of the extended haplotypes surrounding the various K13 alleles confirmed that both common and rare alleles have arisen recently and independently in many different parasite populations.
Miotto and colleagues employed a genome-wide association study (GWAS) to identify genes associated with the slow clearance phenotype. This approach verified that the presence of any K13 propeller mutant allele was by far the strongest signal associated with slow clearance, as had been noted earlier . However, they also identified four other loci on three different chromosomes that showed very strong association; all had prior connections with antimalarial resistance [fd (ferredoxin), arps10 (apicoplast ribosomal protein S10), mdr2 (multidrug resistance protein 2), and crt (chloroquine resistance transporter)]. When the carriage of “any K13 allele” was included as a covariate in the GWAS, these other loci made only modest contributions to the association, so this correlation with slow clearance largely reflects their population-based relationship with K13 alleles, not a functional role in artemisinin resistance.
Genetic complexity and artemisinin use
In Cambodia and many other locations, artemisinins have been used and misused over an extended period, and malaria transmission has declined rapidly since 2008. Under these circumstances, it is not surprising that some of these parasite populations showed a reduced genetic complexity characteristic of founder populations . To understand this observation more clearly, the authors defined seven founder populations that each carried a prominent K13 mutant allele. Despite their overall genetic difference, these populations frequently shared the same alleles of the fd, arps10, mdr2, and crt loci originally identified in the GWAS of all parasites under study. Analysis of these genetic relationships demonstrated that the auxiliary alleles must have evolved in a population ancestral to the founder populations, and the different K13 mutant alleles were selected later in each founder population, but on this common genetic background. Viewed in this light, the auxiliary alleles constitute a “permissive genome” common among parasites in the Mekong basin, a foundation on which K13 mutants not only occur, but survive, and often rise to high prevalence under selection by artemisinins.
This work from the Greater Mekong region makes it clear that the evolution of resistance to artemisinin is still evolving rapidly. Molecular surveys of the prevalence of K13 mutant alleles in Africa and India also suggest that the changes are at a very early stage there. The propeller region of K13 is highly conserved [26,27], and slow parasite clearance has not been observed in either area [12,17,28,29], so mutations in the propeller of the K13 gene were expected to be rare, or absent. Instead, K13 mutants were identified at low prevalence (<5%) in almost all African locations examined [28,30-32]; 25 K13 mutant alleles were identified, 21 that were novel and 4 that had been identified previously in the Mekong. Four sites in India each had a single isolate carrying a different K13 mutant allele; and two of these were novel .
More than 50 different K13 mutant alleles have now been identified, and many more are being reported as geographic surveys are widened. In the Mekong region, the correlation of K13 mutants with slow clearance is strong, and appears to depend upon the presence of the permissive genomes identified by Miotto and colleagues. In other parts of the world, it is not yet known whether the presence of a K13 mutant allele alone is enough to identify slow-clearing parasites. New tools may help to answer this important question. For example, when artemisinin-sensitive parasites were engineered to express a common K13 mutant allele from Cambodia, they were protected against artemisinin exposure in the laboratory [33,34]; testing the African or Indian alleles in this system could be a first step in determining the importance of these novel alleles to the parasite response to artemisinins. In addition, earlier studies have shown that artemisinin treatment increases oxidative stresses in the parasites  and slow-clearing parasites survive the treatment by retarding their progress through the life cycle and activating transcription of several families of genes that protect against these stresses [36,37]. These transcriptional changes are likely to be another manifestation of the “permissive genome” observed in the Mekong parasites. Other genetic changes may be required to support the emergence of K13 mutants in new selective environments, but sets of genes that are adapted to those locations could be identified using the approach described by Miotto and colleagues.
All of these studies have involved collaborations among a very large international group of researchers, combining expertise in clinical management, epidemiology, and genomics, and these studies provide an expanding understanding of the mechanism of resistance to the artemisinins. In the past, molecular markers of resistance have been identified only at the culmination of a long selection process, and a few alleles sufficed to identify resistant parasite populations over wide geographic areas [38-42]. At this stage, it is still too early to be sure whether the K13 molecular signature can be used alone as a valid marker for the slow clearance phenotype worldwide. A few parasite populations carrying common K13 alleles may prove to be best adapted to the selection pressures, increase, and slowly spread into neighboring areas, while most parasites with rare alleles may be transitory, disappearing in the evolutionary struggle . Whatever the outcome, the tools are being assembled to answer that important question as quickly as possible.
From the perspective of public health, patients treated with an ACT still recover in most of the world, so complete failure of ACT treatment is not as widespread as the K13 genotypes in the Mekong region and certainly not in other areas. However, with continued artemisinin pressure, other genetic changes may evolve and render these parasites even more artemisinin resistant. Even more worrying, the efficacy of ACTs requires both components, and when artemisinin loses potency, the partner drug is exposed to far greater selective pressure. When that happens, clinical efficacy of the ACT can diminish rapidly . To short circuit this threat, elimination of P. falciparum parasites in the Mekong region must be a top priority. These new insights on the evolution of artemisinin resistance must also be turned into practical tools to detect artemisinin-resistant parasites in all areas, so that public health measures can be mobilized to contain their emergence or spread long before the clinical efficacy of ACTs is exhausted.
artemisinin combination therapies
- arps10 :
apicoplast ribosomal protein S10 gene, PF3D7_1460900.1
- crt :
chloroquine resistance transporter gene, PF3D7_0709000
- fd :
ferredoxin gene, PF3D7_1318100
genome-wide association study
Kelch 13 gene, PF3D7_1343700
- mdr2 :
multidrug resistance protein 2 gene, PF3D7_1447900
single nucleotide polymorphism
Tracking Resistance to Artemisinin Collaboration
Davies SC, Grant J, Catchpole M. The drugs don’t work - a global threat. Penguin Books, Ltd.; 2013.
Noedl H, Se Y, Schaecher K, Smith BL, Socheat D, Fukuda MM. Evidence of artemisinin-resistant malaria in western Cambodia. N Engl J Med. 2008;359:2619–20.
Dondorp AM, Nosten F, Yi P, Das D, Phyo AP, Tarning J, et al. Artemisinin resistance in Plasmodium falciparum malaria. N Engl J Med. 2009;361:455–67.
World Health Organization. Global plan for containment of artemisinin resistance. http://www.who.int/malaria/publications/atoz/9789241500838/en/index.html. Accessed 22 Feb 2015.
Miotto O, Amato R, Ashley EA, MacInnis B, Almagro-Garcia J, Amaratunga C, et al. Genetic architecture of artemisinin-resistant Plasmodium falciparum. Nat Genet. 2015. doi:10.1038/ng.3189.
Bethell D, Se Y, Lon C, Tyner S, Saunders D, Sriwichai S, 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:e19283.
Starzengruber P, Swoboda P, Fuehrer HP, Khan WA, Hofecker V, Siedl A, et al. Current status of artemisinin-resistant falciparum malaria in South Asia: a randomized controlled artesunate monotherapy trial in Bangladesh. PLoS One. 2012;7:e52236.
Amaratunga C, Sreng S, Suon S, Phelps ES, Stepniewska K, Lim P, et al. Artemisinin-resistant Plasmodium falciparum in Pursat province, western Cambodia: a parasite clearance rate study. Lancet Infect Dis. 2012;12:851–8.
Phyo AP, Nkhoma S, Stepniewska K, Ashley EA, Nair S, McGready R, et al. Emergence of artemisinin-resistant malaria on the western border of Thailand: a longitudinal study. Lancet. 2012;379:1960–6.
Takala-Harrison S, Clark TG, Jacob CG, Cummings MP, Miotto O, Dondorp AM, et al. Genetic loci associated with delayed clearance of Plasmodium falciparum following artemisinin treatment in Southeast Asia. Proc Natl Acad Sci U S A. 2013;110:240–5.
Kyaw MP, Nyunt MH, Chit K, Aye MM, Aye KH, Aye MM, et al. Reduced susceptibility of Plasmodium falciparum to artesunate in southern Myanmar. PLoS One. 2013;8:e57689.
Maiga AW, Fofana B, Sagara I, Dembele D, Dara A, Traore OB, 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:23–8.
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.
Witkowski B, Lelievre J, Lopez Barragan MJ, Laurent V, Su X-z, Berry A, et al. Increased tolerance to artemisinin in Plasmodium falciparum is mediated by a quiescence mechanism. Antimicrob Agents Chemother. 2010;54:1872–7.
Witkowski B, Khim N, Chim P, Kim S, Ke S, Kloeung N, et al. Reduced artemisinin susceptibility of Plasmodium falciparum ring stages in western Cambodia. Antimicrob Agents Chemother. 2013;57:914–23.
Amaratunga C, Witkowski B, Khim N, Menard D, Fairhurst RM. Artemisinin resistance in Plasmodium falciparum. Lancet Infect Dis. 2014;14:449–50.
Ashley EA, Dhorda M, Fairhurst RM, Amaratunga C, Lim P, Suon S, et al. Spread of artemisinin resistance in Plasmodium falciparum malaria. N Engl J Med. 2014;371:411–23.
Sibley CH, Price RN. Monitoring antimalarial drug resistance: applying lessons learned from the past in a fast-moving present. Int J Parasitol Drugs Drug Resist. 2012;2:126–33.
Anderson TJ, Nair S, Nkhoma S, Williams JT, Imwong M, Yi P, et al. High heritability of malaria parasite clearance rate indicates a genetic basis for artemisinin resistance in western Cambodia. J Infect Dis. 2010;201:1326–30.
Miotto O, Almagro-Garcia J, Manske M, Macinnis B, Campino S, Rockett KA, et al. Multiple populations of artemisinin-resistant Plasmodium falciparum in Cambodia. Nat Genet. 2013;45:648–55.
Cheeseman IH, Miller BA, Nair S, Nkhoma S, Tan A, Tan JC, et al. A major genome region underlying artemisinin resistance in malaria. Science. 2012;336:79–82.
Ariey F, Witkowski B, Amaratunga C, Beghain J, Langlois A-C, Khim N, et al. A molecular marker of artemisinin-resistant Plasmodium falciparum malaria. Nature. 2014;505:50–5.
Takala-Harrison S, Jacob CG, Arze C, Cummings MP, Silva JC, Dondorp AM, et al. Independent emergence of artemisinin resistance mutations among Plasmodium falciparum in Southeast Asia. J Infect Dis. 2015;211:670–9.
Thriemer K, Hong NV, Rosanas-Urgell A, Phuc BQ, do Ha M, Pockele E, et al. Delayed parasite clearance after treatment with dihydroartemisinin-piperaquine in Plasmodium falciparum malaria patients in central Vietnam. Antimicrob Agents Chemother. 2014;58:7049–55.
Tun KM, Imwong M, Lwin KM, Win AA, Hlaing TM, Lin K, et al.Westward expansion of artemisinin resistant P. falciparum across Myanmar. Lancet Infect Dis. 2015. doi: 10.1016/S1473-3099(15)70032-0.
Preston MD, Assefa SA, Ocholla H, Sutherland CJ, Borrmann S, Nzila A, et al. PlasmoView: a Web-based resource to visualise global Plasmodium falciparum genomic variation. J Infect Dis. 2014;209:1808–15.
Plasmodium falciparum genetic variations (SNPs) & allele frequencies. http://explorercat.sanger.ac.uk/ExplorerCat-pgv/pub/executePlugin.action?pluginName=PopGen+Explorer&selectedCatalogId=1. Accessed 22 Feb 2015.
Conrad MD, Bigira V, Kapisi J, Muhindo M, Kamya MR, Havlir DV, et al. Polymorphisms in K13 and falcipain-2 associated with artemisinin resistance are not prevalent in Plasmodium falciparum isolated from Ugandan children. PLoS One. 2014;9:e105690.
Mishra N, Prajapati SK, Kaitholia K, Bharti RS, Srivastava B, Phookan S, et al. Surveillance for artemisinin resistance in Plasmodium falciparum in India using the kelch13 molecular marker. Antimicrob Agents Chemother. 2015. doi: 10.1128/AAC.04632-14.
Kamau E, Campino S, Amenga-Etego L, Drury E, Ishengoma D, Johnson K, et al. K13-propeller polymorphisms in Plasmodium falciparum parasites from sub-Saharan Africa. J Infect Dis. 2014.doi: 10.1093/infdis/jiu608.
Taylor SM, Parobek CM, DeConti DK, Kayentao K, Coulibaly SO, Greenwood BM, et al. Absence of putative artemisinin resistance mutations among Plasmodium falciparum in sub-Saharan Africa: a molecular epidemiologic study. J Infect Dis. 2015;211:680–8.
Torrentino-Madamet M, Fall B, Benoit N, Camara C, Amalvict R, Fall M, et al. Limited polymorphisms in k13 gene in Plasmodium falciparum isolates from Dakar, Senegal in 2012-2013. Malar J. 2014;13:472.
Ghorbal M, Gorman M, Macpherson CR, Martins RM, Scherf A, Lopez-Rubio JJ. Genome editing in the human malaria parasite Plasmodium falciparum using the CRISPR-Cas9 system. Nat Biotechnol. 2014;32:819–21.
Straimer J, Gnadig NF, Witkowski B, Amaratunga C, Duru V, Ramadani AP, et al. Drug resistance. K13-propeller mutations confer artemisinin resistance in Plasmodium falciparum clinical isolates. Science. 2015;347:428–31.
Klonis N, Creek DJ, Tilley L. Iron and heme metabolism in Plasmodium falciparum and the mechanism of action of artemisinins. Curr Opin Microbiol. 2013;16:722–7.
Klonis N, Xie SC, McCaw JM, Crespo-Ortiz MP, Zaloumis SG, Simpson JA, et al. Altered temporal response of malaria parasites determines differential sensitivity to artemisinin. Proc Natl Acad Sci U S A. 2013;110:5157–62.
Mok S, Imwong M, Mackinnon MJ, Sim J, Ramadoss R, Yi P, et al. Artemisinin resistance in Plasmodium falciparum is associated with an altered temporal pattern of transcription. BMC Genomics. 2011;12:391.
Cortese JF, Caraballo A, Contreras CE, Plowe CV. Origin and dissemination of Plasmodium falciparum drug-resistance mutations in South America. J Infect Dis. 2002;186:999–1006.
Roper C, Pearce R, Nair S, Sharp B, Nosten F, Anderson T. Intercontinental spread of pyrimethamine-resistant malaria. Science. 2004;305:1124.
Ariey F, Fandeur T, Durand R, Randrianarivelojosia M, Jambou R, Legrand E, et al. Invasion of Africa by a single pfcrt allele of South East Asian type. Malar J. 2006;5:34.
Plowe CV. The evolution of drug-resistant malaria. Trans R Soc Trop Med Hyg. 2009;103:S11–4.
Roper C, Alifrangis M, Ariey F, Talisuna A, Menard D, Mercereau-Puijalon O, et al. Molecular surveillance for artemisinin resistance in Africa. Lancet Infect Dis. 2014;14:668–70.
Menard D, Ariey F. Towards real-time monitoring of artemisinin resistance. J Infect Dis. 2015: doi:10.1016/s1473-3099(15)70046-0.
This work was completed in the context of the WorldWide Antimalarial Resistance Network (WWARN), with funding from the Bill & Melinda Gates Foundation and the ExxonMobile Foundation. Neither organization has reviewed or influenced the content or conclusions of the article.
The author acknowledges the colleagues from many groups who continue to work together to rapidly define the molecular bases that underlie the responses of P. falciparum to artemisinins and to use that information to sustain the efficacy of artemisinins.
The author declares that she has no competing interests.
The author is Professor Emerita of Genome Sciences at the University of Washington and the Scientific Director of WWARN.
About this article
Cite this article
Sibley, C.H. Observing in real time the evolution of artemisinin resistance in Plasmodium falciparum. BMC Med 13, 67 (2015). https://doi.org/10.1186/s12916-015-0316-5
- Plasmodium falciparum
- Artemisinin Resistance
- Molecular Markers
- Genomic Analysis