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

Antimalarial drug resistance has historically arisen through convergent de novo mutations in Plasmodium falciparum parasite populations in Southeast Asia and South America. For the past decade in Southeast Asia, artemisinins, the core component of first-line antimalarial therapies, have experienced delayed parasite clearance associated with several pfk13 mutations, primarily C580Y. We report that mutant pfk13 has emerged independently in Guyana, with genome analysis indicating an evolutionary origin distinct from Southeast Asia. Pfk13 C580Y parasites were observed in 1.6% (14/854) of samples collected in Guyana in 2016–2017. Introducing pfk13 C580Y or R539T mutations by gene editing into local parasites conferred high levels of in vitro artemisinin resistance. In vitro growth competition assays revealed a fitness cost associated with these pfk13 variants, potentially explaining why these resistance alleles have not increased in frequency more quickly in South America. These data place local malaria control efforts at risk in the Guiana Shield.

The origin of Plasmodium falciparum in South America is controversial. Some studies suggest a recent introduction during the European colonizations and the transatlantic slave trade. Other evidence—archeological and genetic—suggests a much older origin. We collected and analyzed P. falciparum isolates from different regions of the world, encompassing the distribution range of the parasite, including populations from sub-Saharan Africa, the Middle East, Southeast Asia, and South America. Analyses of microsatellite and SNP polymorphisms show that the populations of P. falciparum in South America are subdivided in two main genetic clusters (northern and southern). Phylogenetic analyses, as well as Approximate Bayesian Computation methods suggest independent introductions of the two clusters from African sources. Our estimates of divergence time between the South American populations and their likely sources favor a likely introduction from Africa during the transatlantic slave trade.

Background Plasmodium vivax is the most prevalent human malaria parasite, causing serious public health problems in malaria-endemic countries. Until recently the Duffy-negative blood group phenotype was considered to confer resistance to vivax malaria for most African ethnicities. We and others have reported that P. vivax strains in African countries from Madagascar to Mauritania display capacity to cause clinical vivax malaria in Duffy-negative people. New insights must now explain Duffy-independent P. vivax invasion of human erythrocytes.Methods/Principal FindingsThrough recent whole genome sequencing we obtained ≥70× coverage of the P. vivax genome from five field-isolates, resulting in ≥93% of the Sal I reference sequenced at coverage greater than 20×. Combined with sequences from one additional Malagasy field isolate and from five monkey-adapted strains, we describe here identification of DNA sequence rearrangements in the P. vivax genome, including discovery of a duplication of the P. vivax Duffy binding protein (PvDBP) gene. A survey of Malagasy patients infected with P. vivax showed that the PvDBP duplication was present in numerous locations in Madagascar and found in over 50% of infected patients evaluated. Extended geographic surveys showed that the PvDBP duplication was detected frequently in vivax patients living in East Africa and in some residents of non-African P. vivax-endemic countries. Additionally, the PvDBP duplication was observed in travelers seeking treatment of vivax malaria upon returning home. PvDBP duplication prevalence was highest in west-central Madagascar sites where the highest frequencies of P. vivax-infected, Duffy-negative people were reported.Conclusions/SignificanceThe highly conserved nature of the sequence involved in the PvDBP duplication suggests that it has occurred in a recent evolutionary time frame. These data suggest that PvDBP, a merozoite surface protein involved in red cell adhesion is rapidly evolving, possibly in response to constraints imposed by erythrocyte Duffy negativity in some human populations.

Suspected artemisinin resistance in Plasmodium falciparum can be explored by examining polymorphisms in the Kelch (PfK13) propeller domain. Sequencing of PfK13 and other gene resistance markers was performed on 98 samples from Guyana. Five of these samples carried the C580Y allele in the PfK13 propeller domain, with flanking microsatellite profiles different from those observed in Southeast Asia. These molecular data demonstrate independent emergence of the C580Y K13 mutant allele in Guyana, where resistance alleles to previously used drugs are fixed. Therefore, in Guyana and neighboring countries, continued molecular surveillance and periodic assessment of the therapeutic efficacy of artemisinin-based combination therapy are warranted.