Spatio-temporal distribution of antimalarial drug resistant gene mutations in a Plasmodium falciparum parasite population from Kilifi, Kenya: A 25-year retrospective study

Background: Antimalarial drug resistance is a major obstacle to sustainable malaria control. Here we use amplicon sequencing to describe molecular markers of drug resistance in Plasmodium falciparum parasites from Kilifi county in the coastal region of Kenya over a 25-year period. Methods: We performed P. falciparum amplicon sequencing on 1162 malaria-infected blood samples collected between 1994 and 2018 to identify markers of antimalarial drug resistance in the Pfcrt, Pfdhfr, Pfdhps, Pfmdr1, Pfexo, Pfkelch13, plasmepsin 2/3, Pfarps10, Pffd, and Pfmdr2 genes. We further interrogated parasite population structure using a genetic barcode of 101 drug resistance-unrelated single nucleotide polymorphisms (SNPs) distributed across the genomes of 1245 P. falciparum parasites. Results: Two major changes occurred in the parasite population over the 25 years studied. In 1994, approximately 75% of parasites carried the marker of chloroquine resistance, CVIET. This increased to 100% in 1999 and then declined steadily, reaching 6.7% in 2018. Conversely, the quintuple mutation form of sulfadoxine-pyrimethamine resistance increased from 16.7% in 1994 to 83.6% in 2018. Several non-synonymous mutations were identified in the Kelch13 gene, although none of them are currently associated with artemisinin resistance. We observed a temporal increase in the Pfmdr1 NFD haplotype associated with lumefantrine resistance, but observed no evidence of piperaquine resistance. SNPs in other parts of the genome showed no significant temporal changes despite the marked changes in drug resistance loci over this period. Conclusions: We identified substantial changes in molecular markers of P. falciparum drug resistance over 25 years in coastal Kenya, but no associated changes in the parasite population structure.