The recent emergence of artemisinin resistance in the Greater Mekong Subregion poses a major threat to the global effort to control malaria. Tracking the spread and evolution of artemisinin-resistant parasites is critical in aiding efforts to contain the spread of resistance. A total of 417 patient samples from the year 2007, collected during malaria surveillance studies across ten provinces in Thailand, were genotyped for the candidate Plasmodium falciparum molecular marker of artemisinin resistance K13. Parasite genotypes were examined for K13 propeller mutations associated with artemisinin resistance, signatures of positive selection, and for evidence of whether artemisinin-resistant alleles arose independently across Thailand. A total of seven K13 mutant alleles were found (N458Y, R539T, E556D, P574L, R575K, C580Y, S621F). Notably, the R575K and S621F mutations have previously not been reported in Thailand. The most prevalent artemisinin resistance-associated K13 mutation, C580Y, carried two distinct haplotype profiles that were separated based on geography, along the Thai-Cambodia and Thai-Myanmar borders. It appears these two haplotypes may have independent evolutionary origins. In summary, parasites with K13 propeller mutations associated with artemisinin resistance were widely present along the Thai-Cambodia and Thai-Myanmar borders prior to the implementation of the artemisinin resistance containment project in the region.

The Thailand-Cambodia border is the epicenter for drug-resistant falciparum malaria. Previous studies have shown that chloroquine (CQ) and pyrimethamine resistance originated in this region and eventually spread to other Asian countries and Africa. However, there is a dearth in understanding the origin and evolution of dhps alleles associated with sulfadoxine resistance. The present study was designed to reveal the origin(s) of sulfadoxine resistance in Cambodia and its evolutionary relationship to African and South American dhps alleles. We sequenced 234 Cambodian Plasmodium falciparum isolates for the dhps codons S436A/F, A437G, K540E, A581G and A613S/T implicated in sulfadoxine resistance. We also genotyped 10 microsatellite loci around dhps to determine the genetic backgrounds of various alleles and compared them with the backgrounds of alleles prevalent in Africa and South America. In addition to previously known highly-resistant triple mutant dhps alleles SGEGA and AGEAA (codons 436, 437, 540, 581, 613 are sequentially indicated), a large proportion of the isolates (19.3%) contained a 540N mutation in association with 437G/581G yielding a previously unreported triple mutant allele, SGNGA. Microsatellite data strongly suggest the strength of selection was greater on triple mutant dhps alleles followed by the double and single mutants. We provide evidence for at least three independent origins for the double mutants, one each for the SGKGA, AGKAA and SGEAA alleles. Our data suggest that the triple mutant allele SGEGA and the novel allele SGNGA have common origin on the SGKGA background, whereas the AGEAA triple mutant was derived from AGKAA on multiple, albeit limited, genetic backgrounds. The SGEAA did not share haplotypes with any of the triple mutants. Comparative analysis of the microsatellite haplotypes flanking dhps alleles from Cambodia, Kenya, Cameroon and Venezuela revealed an independent origin of sulfadoxine resistant alleles in each of these regions.