National genomic profiling of Plasmodium falciparum antimalarial resistance in Zambian children participating in the 2018 Malaria Indicator Survey
==================================================================================================================================================

* Abebe A. Fola
* Ilinca I. Ciubotariu
* Jack Dorman
* Mulenga C. Mwenda
* Brenda Mambwe
* Conceptor Mulube
* Rachael Kasaro
* Moonga B. Hawela
* Busiku Hamainza
* John M. Miller
* Jeffrey A. Bailey
* William J. Moss
* Daniel J. Bridges
* Giovanna Carpi

## Abstract

The emergence of antimalarial drug resistance is a major threat to malaria control and elimination. Using whole genome sequencing of 282 *P. falciparum* samples collected during the 2018 Zambia National Malaria Indicator Survey, we determined the prevalence and spatial distribution of known and candidate antimalarial drug resistance mutations. High levels of genotypic resistance were found across Zambia to pyrimethamine, with over 94% (n=266) of samples having the *Pfdhfr* triple mutant (N51**I**, C59**R**, and S108**N**), and sulfadoxine, with over 84% (n=238) having the *Pfdhps* double mutant (A437**G** and K540**E**). In northern Zambia, 5.3% (n=15) of samples also harbored the *Pfdhps* A581**G** mutation. Although 29 mutations were identified in *Pfkelch13*, these mutations were present at low frequency (<2.5%), and only three were WHO-validated artemisinin partial resistance mutations: P441**L** (n=1, 0.35%), V568**M** (n=2, 0.7%) and R622**T** (n=1, 0.35%). Notably, 91 (32%) of samples carried the E431**K** mutation in the *Pfatpase6* gene, which is associated with artemisinin resistance. No specimens carried any known mutations associated with chloroquine resistance in the *Pfcrt* gene (codons 72-76). *P. falciparum* strains circulating in Zambia were highly resistant to sulfadoxine and pyrimethamine but remained susceptible to chloroquine and artemisinin. Despite this encouraging finding, early genetic signs of developing artemisinin resistance highlight the urgent need for continued vigilance and expanded routine genomic surveillance to monitor these changes.

Keywords
*   Malaria
*   *Plasmodium falciparum*
*   drug resistance
*   genomic epidemiology
*   Zambia

## Introduction

Among all *Plasmodium* species that infect humans, *P. falciparum* is of the greatest significance, accounting for over 95% of malaria deaths.1,2 The World Health Organization (WHO) recommends front line artemisinin combination therapy (ACT), such as artemether-lumefantrine (AL), artesunate-amodiaquine (ASAQ) or dihydroartemisinin-piperaquine (DHAP), in most African countries for treatment of malaria.3 For prevention, sulfadoxine-pyrimethamine (SP) is recommended for intermittent preventive treatment in pregnancy (IPTp) and for infants (IPTi) living in high-transmission areas.3,4 Antimalarial drugs have played a major role in achieving a significant reduction in the malaria burden globally since 2000 but progress towards malaria elimination has stalled, with resurgence in several endemic countries.5,6 A major threat to control and elimination is the emergence and spread of antimalarial drug resistance.

*P. falciparum* parasites have developed resistance to nearly every available antimalarial drug and resistant strains have spread across malaria endemic countries.7 For example, *P. falciparum* resistance to chloroquine (CQ) emerged in 1957 in Thailand, spread to Southeast Asia in the 1970’s, and by 1982 resistance had spread across the entire African continent.8 Several alleles are associated with CQ resistance around codon 72-76 in the *P. falciparum* CQ resistance transporter gene (*Pfcrt*), including the diagnostic K76**T** mutation.9 Similarly, mutations in the enzymes dihydropteroate synthase (*Pfdhps*) (S436**A**, A437**G**, K540**E**, A581**G**, A613**S**) and dihydrofolate reductase (*Pfdhfr*) (N51**I**, C59**R**, S108**N**, I164**L**) are associated with varying degree of resistance to sulfadoxine and pyrimethamine, respectively,10,11 and are widespread in Africa, Asia, South America, and Oceania.12 Combinations of these mutations i.e., triple *Pfdhfr* mutations of N51**I**, C59**R**, and S108**N**, plus double *Pfdhps* mutations of A437**G**, K540**E** (IRNGE - ‘Quintuple mutant’) confer full resistance to SP.4,13 Moreover, parasites that have the additional *Pfdhps* A581**G** mutation (IRNGEG - ‘Sextuple-mutant’) are associated with enhanced SP resistance *in vitro* that contribute to super SP resistance and IPTp failure.14,15

Since 2008, *P. falciparum* parasites resistant to first-line artemisinin (ART) treatments have emerged in Southeast Asia16,17 and spread in the Greater Mekong Subregion (GMS).18 Studies identified point mutations in the beta-propeller domain of *kelch* 13 (*Pfkelch13*) that were associated with reduced susceptibility to ART and its derivatives, manifested by delayed parasite clearance times.19 To classify ART resistant (ART-R) parasites, WHO provided a list of 9 “validated” and 12 “associated/candidate” *kelch13* ART resistance markers.20 Resistance to partner drugs has also been identified, specifically mutations N86**Y**, Y184**F** and D1246**Y** in the *P. falciparum* multidrug resistance gene 1 (*Pfmdr1*), which have been associated with reduced susceptibility to lumefantrine.21 With reports of decreased ACT efficacy and treatment failure in Africa,22–24 and genotypic identification of an increasing prevalence of WHO-validated *kelch13* ART-R marker (R561**H**),25,26 the threat of antimalarial resistance is increasing such that surveillance should be high priority.27,28

In Zambia, despite the continued use of SP for IPTp and AL for the treatment of uncomplicated malaria since 2002 (the first African country to adopt AL as a first-line treatment policy nationwide), only a few small-scale studies in Southern, Western, and Luapula Provinces have investigated the prevalence of CQ and ACT resistance markers.29,30 However, our recent nationwide genomic study,31 aimed at understanding malaria transmission across Zambia and establishing a baseline for parasite genetic metrics, identified strong positive selection signatures in genes involved in SP and ACT resistance that warrants further investigation into the regional prevalence and spatial variations of drug resistance markers. Moreover, genetic background mutations that could augment ART-R and other novel mutations that may confer antimalarial resistance have yet to be thoroughly assessed in Zambia.

To address these important knowledge gaps and support the Zambian National Malaria Elimination Program, we conducted a genomic surveillance study nested within the 2018 Zambia Malaria Indicator Survey (MIS), which comprises nationwide representative samples.31 We generated 282 *P. falciparum* whole genome sequences (WGS) from seven provinces and mined these data for known and candidate antimalarial resistance mutations across 5 key *P. falciparum* drug resistance genes *(Pfdhfr*, *Pfdhps*, *Pfkelch13*, *Pfcrt*, and *Pfmdr1*) and other gene that may contribute to ART-R phenotypes: apicoplast ribosomal protein S10 precursor *(Pfarps10),* multidrug resistance protein 2 *(Pfmdr2),* ferredoxin (*Pffd),* adaptor protein 2 complex subunit mu gene (*Pfap2mu*), ubiquitin carboxyl-terminal hydrolase 1 gene (*Pfubp1*),32,33 and reticulum Ca2+ ATPase (*Pfatp6*).34 This study defines the most complete genetic landscape of antimalarial drug resistance markers in Zambia, allows spatial and temporal trends to b identified, and provides the foundation for future studies.

## Results

### Whole genome sequencing and mining of drug resistance markers

A total of 282 specimens collected during the 2018 Zambia MIS from children younger than five years of age (**Figure 1**), and were whole genome sequenced (**Figure S1**) with high coverage across 13 *P. falciparum* genes associated with antimalarial drug resistance (**Table S1**). Within the open-reading frames of these genes, 489 non-synonymous (NS) mutations (**Table S2** and **Figure S2**) were identified with variable spatial frequencies across Zambia (**Table S3**). The prevalence of key mutations associated with antimalarial drug resistance is shown in **Table 1**.

![Figure 1:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2024/08/06/2024.08.05.24311512/F1.medium.gif)

[Figure 1:](http://medrxiv.org/content/early/2024/08/06/2024.08.05.24311512/F1)

Figure 1: Spatial distribution of samples retained (n=282) for downstream drug resistanc analysis.
**A)** Sample distribution at the cluster level across Zambia. The size of each sample collection cluster (red) is shown in proportion to the cluster sample size. **B)** Sample size per province. Three provinces (Central, Lusaka and Southern) were excluded from provincial prevalence calculations due to low sample size.

View this table:
[Table 1:](http://medrxiv.org/content/early/2024/08/06/2024.08.05.24311512/T1)

Table 1: 
Prevalence of key mutations associated with mono- and multiple antimalarial drug resistance across Zambia in 2018. n = number of samples that carried the mutant allele out of the total sequenced samples (282)

### Chloroquine resistance

The *Pfcrt* K76**T** mutation associated with CQ resistance was not identified in the analyzed samples (**Table 1, Table S2**), indicating reversal of CQ sensitive *P. falciparum* strains across Zambia.31,35 Seven NS mutations not associated with CQ resistance were identified at very low frequencies in other codons across the *Pfcrt* gene (PF3D7_0709000), but 94.0% (265/282) of sequenced samples did not carry any of these mutations (**Figure S3, Table S2**). All sequenced specimens were wild-type (N) at codon 86 in the multidrug resistance transporter *Pfmdr1* gene, which has been shown to enhance resistance to CQ in some genetic backgrounds,36 further supporting reversal to CQ sensitive parasites in Zambia.

### Antifolate resistance

Key mutations (N51**I**, C59**R**, S108**N**) in the *Pfdhfr* gene were identified in 98.2% (277/282), 96.5% (272/282), and 99.6% (281/282) of samples, respectively, with 94.6% (267/282) classified as a *Pfdhfr* triple (**IRN**) mutants (**Table 1 and Figure 2A**). Moreover, key mutations A437**G** and K540**E**, in the *Pfdhps* gene were identified in 93.3% (263/282) and 87.6% (247/282) of sequenced samples, respectively, with 84.0% (237/282) classified as a *Pfdhps* double (**GE**) mutant (**Table 1 and Figure 2C**). We found some degree of spatial heterogeneity for both *Pfdhfr* triple mutants and *Pfdhps* double mutants at the cluster level across the country (**Figure 2B** and **Figure 2D**, respectively).

![Figure 2.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2024/08/06/2024.08.05.24311512/F2.medium.gif)

[Figure 2.](http://medrxiv.org/content/early/2024/08/06/2024.08.05.24311512/F2)

Figure 2. Prevalence of *Pfdhfr* and *Pfdhps* mutations across Zambia.
UpSet plots showing the number of times each combination of mutations was seen for *Pfdhfr* (**A**) and *Pfdhps* (**C**). The lines and circles below the bars represent the different combinations of resistant genotype in individual samples. **B**) Spatial prevalence of pyrimethamine resistant triple mutant (*Pfdhf* 51**I***/*59**R**/108**N**, **IRN**), **D**) Spatial prevalence of sulfadoxine resistant double mutant (*Pfdhp* 437**G***/*540**E**, **GE**). Color code heat map shows prevalence at the district level, while white indicates districts where no samples were available.

Overall, 83.0% (234/282) of sequenced samples were typed as quintuple mutants (*Pfdhfr* IRN, *Pfdhps* GE, **Figure 3A**), a genotype associated with near complete SP resistance (**Table 1**). In the *Pfdhps* gene, the A581**G** mutation was found in 5.3% (15/282) of the analyzed samples with high geographic variation across Zambia and with clustering in the adjacent Luapula and Northern Provinces (**Figure 3B**). There was no evidence for the presence of *Pfdhfr* I164**L** or *Pfdhps* A613**S**/**T**, mutations that can enhance the quintuple mutant SP resistance profile.

![Figure 3.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2024/08/06/2024.08.05.24311512/F3.medium.gif)

[Figure 3.](http://medrxiv.org/content/early/2024/08/06/2024.08.05.24311512/F3)

Figure 3. Cumulative prevalence of *Pfdhfr-dhps* mutations (A) and spatial prevalence of *Pfdhps* A581G (B) across Zambia.
Color code heat map shows prevalence at the district-level, while white indicates districts where no samples were available.

### Partial artemisinin resistance

A total of 29 mutations were identified across the entire *kelch13* gene, of which 8 were located inside the propeller domain (**Figure 4**). Of these, three mutations (P441**L**: n**=**1 from Western Province, V568**M**: n=2 from Luapula Province and R622**T**: n**=**1 from Western Province) were WHO-validated mutations associated with partial artemisinin resistance. Apart from the K189**T** mutation that was found in 20.6% (58/282) of sequenced samples and was reported in other clinical studies, all other mutations were found at low frequency (<2.5%). The A578**S** mutation, which was identified in seven specimens (2.5%), has been commonly reported in other African countries, although this mutation is not associated with ART resistance in vitro and/or delayed parasite clearance. Overall, 35.1% (99/282) of sequenced samples carried one or more *Pfkelch13* mutations with variable prevalence at the provincial level (**Table S3**), suggesting high polymorphism in the *Pfkelch13* gene due to increased ACT pressure in Zambia.

![Figure 4.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2024/08/06/2024.08.05.24311512/F4.medium.gif)

[Figure 4.](http://medrxiv.org/content/early/2024/08/06/2024.08.05.24311512/F4)

Figure 4. Non-synonymous mutations identified across the *kelch13* gene in Zambia.
Colors indicate mutations as outside (blue) or within (orange) the propeller domain.

### Lumefantrine resistance

*P. falciparum* parasite harboring N86 (wild-type), 184F (mutant), and D1246 (wild-type) genotypes in the multi-drug resistance gene 1 (*Pfmdr1*), are associated with decreased sensitivity to lumefantrine. Our analysis revealed that 99.3% (280/282) of the samples carried N86, 55.0% (155/282) of the samples carried 184F, and 98.2% (277/282) carried D1246 (**Figure 5A**). Overall, 53.9% (152/282) of samples carried N86/184**F**/D1246 (N**F**D) genotype with marked spatial variation across districts (**Figure 5B**).

![Figure 5.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2024/08/06/2024.08.05.24311512/F5.medium.gif)

[Figure 5.](http://medrxiv.org/content/early/2024/08/06/2024.08.05.24311512/F5)

Figure 5. Prevalence of *Pfmdr1* mutations associated with decreased sensitivity to lumefantrine.
**A)** Cumulative prevalence of *Pfmdr1* haplotypes, **B)** Spatial prevalence of th *Pfmdr1* 184**F** mutation. Color code heat map shows the prevalence at the district level, while white indicates districts where no samples were available.

### Additional ART resistance-associated mutations

Several mutations in the *Pfatp6* gene, a sarcoplasmic and endoplasmic reticulum Ca2+ ATPase (SERCA)-type protein previously associated with artemisinin resistance were identified (**Table S2**). Most mutations occurred at low frequency except for *Pfatp6* N569**K** (54.6%), E431**K** (32.3%), L402**V** (9.6%), and H243**Y** (3.5%). Of the mutations previously shown to mediate resistance, E431**K** (32.3%) and A623**E** (0.35%) were identified, while S769**N** and L263**E** were not reported (**Figure S4**). No other mutations were identified in *Pfarps10*, *Pffd*, *Pfmdr2*, *Pfpib*, *Pfpp*, *Pfap2mu*, *Pfubp1*, and *Pfcrt*, that could potentially augment artemisinin resistance. Several mutations of variable but generally low frequency (**Table S2**) and variable spatial distribution (**Table S3**) were identified in these genes (**Figure S5**).

## Discussion

Intense control activities, especially the continuous use of effective antimalarial drugs, apply significant selective pressure on the malaria parasite population. As transmission declines, most infected individuals carry single clones promoting a higher rate of inbreeding37 that favors the spread of drug resistance phenotypes when they arise.38,39 The independent emergence or spread of artemisinin resistance in Zambia due to drug selection pressure could significantly increase the malaria burden leading to malaria resurgence and increased morbidity and mortality.40 Close monitoring of the efficacy of available antimalarial drugs and improved surveillance for resistance mutations will be key to inform control strategies such that rapid action can be taken to mitigate the impact and slow or prevent the spread of drug resistant parasites.22,23

This nationally representative spatial analysis of antimalarial drug resistance mutations supports several conclusions that together strongly suggest that individual genes are under markedly different selection pressures. Firstly, there was a noticeable lack of mutations in *Pfcrt* codons 72-76 and a low number of polymorphisms (only seven unique NS mutations) across the gene. Zambia withdrew CQ treatment and therefore drug selective pressure on the *Pfcrt* gene in 2003 allowing reversion to the wild-type. While other studies have previously identified this reversion,35,41 the selection patterns in Zambia lacked a commonly selected region on chromosome 7 (*Pfcrt*), contrasting with parasite populations from other regions. Similarly, we did not observe selection signatures in *Pfaat1*, the second important transporter gene for chloroquine resistance.42 This is consistent with other African countries where CQ withdrawal resulted in declines in CQ resistance alleles and reductions in CQ median IC50 values.43,44 While this finding is encouraging in terms of the potential to reintroduce CQ as part of the chemotherapeutic arsenal in Zambia, preferably as a combination therapy, however additional phenotype-genotype association studies would be needed to confirm the susceptibility of this wild-type strain to CQ.

In contrast to CQ, SP resistance markers were highly prevalent throughout Zambia, suggesting strong selective pressure is maintained through national IPTp implementation and high private sector SP utilization without prescription for self-medication of suspected malaria.45 Pyrimethamine associated resistance in *Pfdhfr* was very high, with 95% of samples having triple mutants (IRN) and codon S108N (99.6%) approaching fixation with negligible spatial variation. A similar picture was observed for sulfadoxine associated *Pfdhps* mutations, with 84% double mutants (codons A437**G** and K540**E**). Overall, 82.9% of samples were *Pfdhfr* and *Pfdhps* quintuple mutants *Pfdhfr-dhps* **(**IRNGE) correlating with full genotypic SP resistance and expected treatment failure. Furthermore, with a concentration in Luapula and Northern Provinces, 5.3% (15/282) of samples also carried the *Pfdhps* A581**G** mutation in addition to the *Pfdhfr-dhps* IRNGE background. This genotype confers extreme SP resistance and is of concern for SP efficacy in Zambia. While some variation between this study and historical data29,46 may be explained by study differences (subject selection, sites, implementation period) and data type (PCR genotyping vs. WGS), overall a marked increase in SP genotypic resistance has occurred. Our findings are consistent with our previous evidence of positive selection for SP markers,31 as well as with other African countries where *Pfdhfr*–*Pfdhps* quintuple mutant prevalence is high, while the sextuple mutant remains rare,45,47 albeit with high spatial heterogeneity.4,48 While overall SP resistance is clearly high in Zambia, two mutations (*Pfdhfr* I164**L**49 and *Pfdhps* A613**S**/**T**47,50) that confer even higher SP resistance were not identified in Zambia.

The high prevalence of SP resistance (>90%) are in children younger than five years) indicates a strong selective pressure has been applied to these genes, even though SP is primarily only used for IPTp. The WHO recommends that countries withdraw SP for IPTp use when the prevalence of *Pfdhps* K540**E** is >95% and *Pfdhps* A581**G** is >10%.11 At 87.6% and 5% respectively, Zambia as a country remains below these thresholds, although some districts e.g., Nchelenge and Mansa in Luapula Province, did exceed them. Based on these findings, SP should continue to be used for IPTp. In contrast, the WHO threshold for SP-based IPTi withdrawal is when K540**E** is >50% and thus, SP would not be recommended to be used for IPTi in Zambia at this time.

Until novel therapies are developed, maintaining the efficacy of ART based treatments is fundamental to global control and elimination efforts. In Zambia, AL has been used as a first line combination antimalarial treatment for uncomplicated *P. falciparum* malaria since 2002.51 Considering the historical use and importance of ART to malaria control in Zambia, increased polymorphisms in the *kelch13* gene and the identification three WHO-validated mutations associated with partial artemisinin resistance suggest an early signal of partial artemisinin resistance in Zambia. Close monitoring of local emerging or spreading *kelch13* mutations (R561**H**, A675**V** and C469**Y**)52–54 that were recently reported from East Africa and confer partial artemisinin resistance is warranted. While not unexpected, it was also encouraging to note that no mutations (*Pfcrt* I356**T**, *Pffd* D193**Y**, *Pfmdr2* T484**I**, *Pfap2mu* S160**N** and *Pfubp1* E1528**D)** associated with ART resistance in Southeast Asia55 were identified. Nevertheless, considering the variation between Asian and African parasite populations,56 it is possible that other Africa-specific mutations may augment ART resistance. Similarly, we must continue to track all mutations in any key genes, irrespective of their genotypic resistance status. For example, 26 *Pfkelch13* mutations were identified, including one (A578**S**) that has been commonly reported in Africa,2,29 the implications of which remains unclear. Unfortunately, while ART appears to be efficacious, the main partner drug lumefantrine does not fare as well, with all but two specimens containing one or more key mutations in *Pfmdr1*. In fact, more than 50% of all specimens carried mdr1 (N**F**D) haplotype. This confirms results provided by other studies performed in the Southern and Western Provinces of Zambia2,29,57 but, as with SP, the trend is that genotypic resistance is increasing. While this may not correlate with ACT clinical treatment failure with AL, it does potentially remove the partner drug from the combination therapy leaving ART exposed as monotherapy. Such an environment would be primed to enable rapid selection and spread of ART resistance irrespective of resistance evolving independently or through an introduction event into Zambia.

In summary, this study support two worrying and two encouraging conclusions with respect to antimalarial drug resistance. Firstly, Zambia has very high, and for some loci almost fixed, resistance to SP. While still under WHO recommended limits, this warrants further SP efficacy studies in pregnancy to assess the drugs’ ability to reduce deleterious maternal and birth outcomes, especially in Luapula and Northern Provinces where WHO frequency thresholds were crossed. Secondly, there are also very high levels of genotypic resistance to lumefantrine, the main ART partner ACT drug used in Zambia. While therapeutic efficacy studies have not identified significant treatment failure several years after these samples were collected, it may be prudent to switch to an alternative ACT in the near future, or at least prepare for a switch should treatment failures occur. Finally, there is some encouraging data, namely that CQ sensitivity has been restored and there is no evidence of ART resistance in Zambia. Together these findings along with the recent evidence of strong positive selection signatures genes involved in sulfadoxine-pyrimethamine and artemisinin combination therapies drug resistance31 highlight the need of sustained surveillance of antimalarial drug resistance across the country. Furthermore, this work underlines the utility of high-quality genomic surveillance, which if performed and acted upon, gives every chance of effective malaria treatment continuing for the foreseeable future despite the constant threat of drug resistance. Without surveillance, resistance will only be detected following treatment failure, at which point options to respond will be limited.

## Methods

### Sample selection and whole-genome sequencing

This work is a secondary analysis of a larger parent study focused on understanding malaria transmission across Zambia and establishing a baseline for parasite genetic metrics.31 Briefly, whole-genome sequencing was performed on 459 *P. falciparum* PET-PCR positive dried blood spots (DBS) collected from children in all ten provinces of Zambia as part of the 2018 Zambia National Malaria Indicator Survey (MIS).58 The samples were processed for genomic DNA extraction as previously described.31, 59 We adopted a 4-plex hybrid capture method using SeqCap EZ custom probes60 to selectively enrich *P. falciparum* genomes prior to sequencing as previously described.31 Genomic libraries and hybridization capture were prepared using a modified Roche/Nimblegen SeqCap EZ Library Short Read protocol and sequenced on Illumina NextSeq 6000 (2 × 101-bp) at Yale Center for Genomic Analysis with a target of 30 million reads per sample. Raw sequence reads are available at the Sequence Read Archive (PRJNA932927).

### Variant identification

We identified *P. falciparum* genomic variation from whole genome sequence data with the pipeline previously described.31,61 Briefly, Illumina raw paired-end reads were aligned to the *P. falciparum* 3D7 reference genome with BWA-MEM 0.7.1762 and removed using Picard Tools 2.20.8. For this study, variant calling was performed only on samples with >30% *P. falciparum* 3D7 reference genome with >5X coverage, resulting in a total of 282 *P. falciparum* samples. Variants were called using GATK v4.1.4.121 following best practices ([https://software.broadinstitute.org/gatk/best-practices](https://software.broadinstitute.org/gatk/best-practices)). We used GATK HaplotypeCaller in GVCF mode to call single-sample variants (ploidy 2 and standard-min-confidence-threshold for calling = 30), followed by GenotypeGVCFs to genotype the parasites. Prior to variant filtering, we scored 1,219,517 SNPs with a VQSLOD >0 across the 282 genomes. The VCF was functionally annotated with SnpEff v4.3 (build 2017-11-24 10:18). Variants removed included those located in telomeric and hypervariable regions ([ftp://ngs.sanger.ac.uk/production/malaria/pf-crosses/1.0/regions-20130225.onebased.txt](ftp://ngs.sanger.ac.uk/production/malaria/pf-crosses/1.0/regions-20130225.onebased.txt)), SNPs with >20% missingness, and minor allele frequency (MAF) >0.02, leaving a total of 27,163 high quality biallelic SNPs.

### Mining drug resistance loci and estimating mutation prevalence

SNPs located in 13 genes associated with antimalarial drug resistance (**Table S1**) were identified using the *VariantAnnotation* R package. All non-synonymous mutations with high read coverage (>30x) were identified and classified as mutant (heterozygous or homozygous mutant) or wild-type (homozygous reference). The reference allele for *Pfdhps* 437 encodes the mutant allele and was therefore re-coded to A437**G** for clarity as the reference carries a mutant allele unlike all other alleles. The prevalence of each mutation was calculated as (p = m/n*100, where p = prevalence, m = number of infections with mutant alleles, n = number of successfully genotyped infections) using R software version 4.2.0. Mutant combinations were plotted and visualized using the UpSet package in R63 and maps were created using the sf package in R.64

### Ethical statement

All study participants, and/or their parents or legal guardians, provided written informed consent, and this study was conducted with the approval of the Biomedical Research Ethics Committee from the University of Zambia (Ref 011-02-18) and from the Zambian National Health Research Authority.

## Data availability

The sequence data for the parent project are available in the NCBI Sequence Read Archive, in BioProject PRJNA932927.

## Supporting information

Supplementary Tables 1-3 [[supplements/311512_file03.xlsx]](pending:yes)

## Funding information

This work was supported by funds to G.C. from the Purdue Department of Biological Sciences. Partial funding was provided by the Bill & Melinda Gates Foundation through a grant to PATH (OPP1134518 / INV-009984 to D.J.B.) The Southern and Central Africa International Center of Excellence for Malaria Research was supported by funding from the National Institute of Allergy and Infectious Diseases (U19AI089680 to W.J.M.).

## Authors’ contributions

G.C., W.J.M and D.J.B. contributed to funding acquisition, project resources and supervision. G.C., A.A.F., W.J.M., and D.J.B., conceived and designed the study. A.A.F., D.J.B. and G.C., coordinated sample selection and curation. M.C.M., B.M., C.M., R.K., M.B.H., B.H., J.M.M. and D.J.B. collected samples and epidemiological data. A.A.F., I.C., and J.D. performed laboratory analysis. A.A.F., D.J.B, G.C. contributed to formal genomic analysis, visualization, interpretation and writing the original draft. All authors contributed to review and editing the manuscript.

## SUPPLEMENTAL MATERIAL

### SUPPLEMENTAL FIGURES

![Figure S1.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2024/08/06/2024.08.05.24311512/F6.medium.gif)

[Figure S1.](http://medrxiv.org/content/early/2024/08/06/2024.08.05.24311512/F6)

Figure S1. 
Read coverage per all identified mutations associated with antimalarial drug resistance across 282 samples. **A**) Density plot showing cumulative read depth and dashed red lines indicate median coverage=45X. **B**) Read coverage per mutation for top 5 known antimalarial drug resistance genes.


![Figure S2.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2024/08/06/2024.08.05.24311512/F7.medium.gif)

[Figure S2.](http://medrxiv.org/content/early/2024/08/06/2024.08.05.24311512/F7)

Figure S2. 
Number of mutations per gene and read depth. **A**) Number of mutations identified across 13 *P. falciparum* genes from 282 samples. **B**) Drug resistance gene size in base pairs.

![Figure S3.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2024/08/06/2024.08.05.24311512/F8.medium.gif)

[Figure S3.](http://medrxiv.org/content/early/2024/08/06/2024.08.05.24311512/F8)

Figure S3. Frequency of mutations in the *Pfcrt* gene.
The intersection size bar plots represent the number of parasites carrying different combinations of antimalarial drug resistanc associated mutations in the *Pfcrt* gene. The circles below the bars represent the different combinations of resistant genotypes in individual samples. The set size represents the number of samples in which individual drug resistance-conferring markers were genotyped. Gene position and nucleotide change for each mutation are shown in Table S1. None = samples carrying synonymous mutations or wild-type alleles

![Figure S4.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2024/08/06/2024.08.05.24311512/F9.medium.gif)

[Figure S4.](http://medrxiv.org/content/early/2024/08/06/2024.08.05.24311512/F9)

Figure S4. Prevalence of 37 mutations in *Pfatpase6* from 282 Zambian samples ordered by genomic loci.
Colours indicate whether the mutations are validated or not as resistance markers for ACT.

![Figure S5.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2024/08/06/2024.08.05.24311512/F10.medium.gif)

[Figure S5.](http://medrxiv.org/content/early/2024/08/06/2024.08.05.24311512/F10)

Figure S5. Prevalence of different mutations across genes associated with partial artemisinin resistance other than *Kelch13* gene.
We identified several mutations across *PfATPase6* protein, apicoplast ribosomal protein S10 precursor gene (*Pfarps10*), (*Pfmdr*) protein 2 codon, AP-2 complex subunit mu gene (*Pfap2-mu*), ubiquitin carboxyl-terminal hydrolase 1 gene (*Pfubp-1*) and Phosphoinositide-binding protein (*Pfpib7*). No mutation observed in genes: *Pfcrt* I356**T**, *Pffd* D193**Y**, *Pfmdr2* T484I, *Pfap2-mu* S160N/T and *Pfubp-1* E1528D – known background mutations associated with ART in SEA. Also, no mutation was observed in the AP-2 complex subunit mu gene (ap2-mu, S160N/T) and the ubiquitin carboxyl-terminal hydrolase 1 gene (ubp-1, E1528D), mutation that is associated with delayed ACT clearance in Africa.

### SUPPLEMENTAL TABLES

*Supplementary tables are compiled into a single file for ease of viewing*

**Table S1. Description of 13 *P. falciparum* drug resistance genes across successfully sequenced samples and included in this analysis and the drugs with which they are associated.**

See uploaded file.

**Table S2. Sequencing coverage across successfully sequenced samples and cumulative prevalence of all mutations across 13 genes in Zambia.** n = number of samples successfully sequenced samples

See uploaded file.

**Table S3. Prevalence of all mutations across 13 genes per region in Zambia.** n = number of samples successfully sequenced samples

See uploaded file.

## Acknowledgments

The authors are grateful to the Zambian communities, particularly the volunteers and their families, for providing samples during the MIS. We would like to thank the staff of the Zambia National Malaria Elimination Centre for their generous support, especially the field researchers who conducted the nationwide survey.

*   Received August 5, 2024.
*   Revision received August 5, 2024.
*   Accepted August 6, 2024.


*   © 2024, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/)

## References

1.  1.Blasco, B., Leroy, D. & Fidock, D. A. Antimalarial drug resistance: linking *Plasmodium falciparum* parasite biology to the clinic. Nat. Med. 23, 917–928 (2017).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nm.4381&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28777791&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F06%2F2024.08.05.24311512.atom) 

2.  2.Menard, D. & Dondorp, A. Antimalarial Drug Resistance: A Threat to Malaria Elimination. Cold Spring Harbor Perspectives in Medicine vol. 7 a025619 Preprint at doi:10.1101/cshperspect.a025619 (2017).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTQ6ImNzaHBlcnNwZWN0bWVkIjtzOjU6InJlc2lkIjtzOjExOiI3LzcvYTAyNTYxOSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDI0LzA4LzA2LzIwMjQuMDguMDUuMjQzMTE1MTIuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

3.  3.Conrad, M. D. & Rosenthal, P. J. Antimalarial drug resistance in Africa: the calm before the storm? Lancet Infect. Dis. 19, e338–e351 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1473-3099(19)30261-0&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31375467&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F06%2F2024.08.05.24311512.atom) 

4.  4.Amimo, F., et al. *Plasmodium falciparum* resistance to sulfadoxine-pyrimethamine in Africa: a systematic analysis of national trends. BMJ Glob Health 5, (2020).
    
    

5.  5.Dhiman, S. Correction to: Are malaria elimination efforts on right track? An analysis of gains achieved and challenges ahead. Infect Dis Poverty 8, 19 (2019).
    
    

6.  6.Makenga, G. et al. Prevalence of malaria parasitaemia in school-aged children and pregnant women in endemic settings of sub-Saharan Africa: A systematic review and meta-analysis. Parasite Epidemiology and Control 11, e00188 (2020).
    
    

7.  7.Thu, A. M., Phyo, A. P., Landier, J., Parker, D. M. & Nosten, F. H. Combating multidrug-resistant *Plasmodium falciparum* malaria. FEBS J. 284, 2569–2578 (2017).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/febs.14127&link_type=DOI) 

8.  8.Wellems, T. E. & Plowe, C. V. Chloroquine-resistant malaria. J. Infect. Dis. 184, 770–776 (2001).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1086/322858&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11517439&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F06%2F2024.08.05.24311512.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000170992000014&link_type=ISI) 

9.  9.Djimdé, A. et al. A Molecular Marker for Chloroquine-Resistant Falciparum Malaria. New England Journal of Medicine vol. 344 257–263 Preprint at doi:10.1056/nejm200101253440403 (2001).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJM200101253440403&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11172152&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F06%2F2024.08.05.24311512.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000166545200003&link_type=ISI) 

10. 10.Kublin, J. G. et al. Molecular markers for failure of sulfadoxine-pyrimethamine and chlorproguanil-dapsone treatment of *Plasmodium falciparum* malaria. J. Infect. Dis. 185, 380–388 (2002).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1086/338566&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11807721&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F06%2F2024.08.05.24311512.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000173458900014&link_type=ISI) 

11. 11.Okell, L. C., Griffin, J. T. & Roper, C. Mapping sulphadoxine-pyrimethamine-resistant *Plasmodium falciparum* malaria in infected humans and in parasite populations in Africa. Sci. Rep. 7, 7389 (2017).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41598-017-06708-9&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F06%2F2024.08.05.24311512.atom) 

12. 12.Abdul-Ghani, R., Farag, H. F. & Allam, A. F. Sulfadoxine-pyrimethamine resistance in *Plasmodium falciparum*: a zoomed image at the molecular level within a geographic context. Acta Trop. 125, 163–190 (2013).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.actatropica.2012.10.013&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23131424&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F06%2F2024.08.05.24311512.atom) 

13. 13.Desai, M. et al. Impact of Sulfadoxine-Pyrimethamine Resistance on Effectiveness of Intermittent Preventive Therapy for Malaria in Pregnancy at Clearing Infections and Preventing Low Birth Weight. Clinical Infectious Diseases vol. 62 323–333 Preprint at doi:10.1093/cid/civ881 (2016).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/cid/civ881&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26486699&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F06%2F2024.08.05.24311512.atom) 

14. 14.van Eijk, A. M. et al. Effect of *Plasmodium falciparum* sulfadoxine-pyrimethamine resistance on the effectiveness of intermittent preventive therapy for malaria in pregnancy in Africa: a systematic review and meta-analysis. The Lancet Infectious Diseases vol. 19 546– 556 Preprint at doi:10.1016/s1473-3099(18)30732-1 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1473-3099(18)30732-1&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F06%2F2024.08.05.24311512.atom) 

15. 15.Gutman, J. et al. The A581G mutation in the gene encoding *Plasmodium falciparum* dihydropteroate synthetase reduces the effectiveness of sulfadoxine-pyrimethamine preventive therapy in Malawian pregnant women. J. Infect. Dis. 211, 1997–2005 (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/infdis/jiu836&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25564249&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F06%2F2024.08.05.24311512.atom) 

16. 16.Phyo, A. P. et al. Emergence of artemisinin-resistant malaria on the western border of Thailand: a longitudinal study. Lancet 379, 1960–1966 (2012).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(12)60484-X&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22484134&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F06%2F2024.08.05.24311512.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000304757700030&link_type=ISI) 

17. 17.Imwong, M. et al. The spread of artemisinin-resistant *Plasmodium falciparum* in the Greater Mekong subregion: a molecular epidemiology observational study. Lancet Infect. Dis. 17, 491–497 (2017).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1473-3099(17)30048-8&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28161569&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F06%2F2024.08.05.24311512.atom) 

18. 18.Roberts, L. Malaria wars. Science 352, 398–402, 404–5 (2016).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIzNTIvNjI4NC80MDQiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyNC8wOC8wNi8yMDI0LjA4LjA1LjI0MzExNTEyLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

19. 19.Miotto, O. et al. Multiple populations of artemisinin-resistant *Plasmodium falciparum* in Cambodia. Nat. Genet. 45, 648–655 (2013).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.2624&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23624527&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F06%2F2024.08.05.24311512.atom) 

20. 20.World Health Organization. Report on Antimalarial Drug Efficacy, Resistance and Response: 10 Years of Surveillance (2010-2019). (World Health Organization, 2020).
    
    

21. 21.Grais, R. F. et al. Molecular markers of resistance to amodiaquine plus sulfadoxine– pyrimethamine in an area with seasonal malaria chemoprevention in south central Niger. Malar. J. 17, 1–9 (2018).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s12936-018-2247-z&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F06%2F2024.08.05.24311512.atom) 

22. 22.Ekland, E. H. & Fidock, D. A. Advances in understanding the genetic basis of antimalarial drug resistance. Curr. Opin. Microbiol. 10, 363–370 (2007).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.mib.2007.07.007&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17709280&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F06%2F2024.08.05.24311512.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000250104100011&link_type=ISI) 

23. 23.Guyant, P. et al. Past and new challenges for malaria control and elimination: the role of operational research for innovation in designing interventions. Malar. J. 14, 279 (2015).
    
    

24. 24.Arya, A., Kojom Foko, L. P., Chaudhry, S., Sharma, A. & Singh, V. Artemisinin-based combination therapy (ACT) and drug resistance molecular markers: A systematic review of clinical studies from two malaria endemic regions – India and sub-Saharan Africa. International Journal for Parasitology: Drugs and Drug Resistance vol. 15 43–56 Preprint at doi:10.1016/j.ijpddr.2020.11.006 (2021).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ijpddr.2020.11.006&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=33556786&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F06%2F2024.08.05.24311512.atom) 

25. 25.Uwimana, A. et al. Emergence and clonal expansion of in vitro artemisinin-resistant *Plasmodium falciparum* kelch13 R561H mutant parasites in Rwanda. Nat. Med. 26, 1602– 1608 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41591-020-1005-2&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32747827&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F06%2F2024.08.05.24311512.atom) 

26. 26.Uwimana, A. et al. Author Correction: Emergence and clonal expansion of in vitro artemisinin-resistant *Plasmodium falciparum* kelch13 R561H mutant parasites in Rwanda. Nat. Med. 27, 1113–1115 (2021).
    
    

27. 27.Takala-Harrison, S. & Laufer, M. K. Antimalarial drug resistance in Africa: key lessons for the future. Ann. N. Y. Acad. Sci. 1342, 62–67 (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/nyas.12766&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25891142&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F06%2F2024.08.05.24311512.atom) 

28. 28.Ndiaye, Y. D. et al. Genetic surveillance for monitoring the impact of drug use on *Plasmodium falciparum* populations. Int. J. Parasitol. Drugs Drug Resist. 17, 12–22 (2021).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ijpddr.2021.07.004&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=34333350&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F06%2F2024.08.05.24311512.atom) 

29. 29.Sitali, L. et al. Surveillance of molecular markers for antimalarial resistance in Zambia: Polymorphism of Pfkelch 13, Pfmdr1 and Pfdhfr/Pfdhps genes. Acta Trop. 212, 105704 (2020).
    
    

30. 30.Ippolito, M. M. et al. Therapeutic Efficacy of Artemether-Lumefantrine for Uncomplicated Falciparum Malaria in Northern Zambia. Am. J. Trop. Med. Hyg. 103, 2224–2232 (2020).
    
    

31. 31.Fola, A. A. et al. Genomics reveals heterogeneous *Plasmodium falciparum* transmission and selection signals in Zambia. Commun. Med. 4, 67 (2024).
    
    

32. 32.Henriques, G. et al. Directional selection at the pfmdr1, pfcrt, pfubp1, and pfap2mu loci of Plasmodium falciparum in Kenyan children treated with ACT. J. Infect. Dis. 210, 2001– 2008 (2014).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/infdis/jiu358&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24994911&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F06%2F2024.08.05.24311512.atom) 

33. 33.Adams, T. et al. Prevalence of *Plasmodium falciparum* delayed clearance associated polymorphisms in adaptor protein complex 2 mu subunit (pfap2mu) and ubiquitin specific protease 1 (pfubp1) genes in Ghanaian isolates. Parasites & Vectors vol. 11 Preprint at doi:10.1186/s13071-018-2762-3 (2018).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s13071-018-2762-3&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29530100&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F06%2F2024.08.05.24311512.atom) 

34. 34.Afonso, A. et al. Malaria parasites can develop stable resistance to artemisinin but lack mutations in candidate genes atp6 (encoding the sarcoplasmic and endoplasmic reticulum Ca2+ ATPase), tctp, mdr1, and cg10. Antimicrob. Agents Chemother. 50, 480–489 (2006).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYWFjIjtzOjU6InJlc2lkIjtzOjg6IjUwLzIvNDgwIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjQvMDgvMDYvMjAyNC4wOC4wNS4yNDMxMTUxMi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

35. 35.Mwanza, S. et al. The return of chloroquine-susceptible *Plasmodium falciparum* malaria in Zambia. Malar. J. 15, 584 (2016).
    
    

36. 36.Shafik, S. H., Richards, S. N., Corry, B. & Martin, R. E. Mechanistic basis for multidrug resistance and collateral drug sensitivity conferred to the malaria parasite by polymorphisms in *PfMDR1* and *PfCRT*. PLoS Biol. 20, e3001616 (2022).
    
    

37. 37.Wirth, D. F. & Alonso, P. L. Malaria: Biology in the Era of Eradication. (Cold Spring Harbor Laboratory Press, 2017).
    
    

38. 38.Carrasquilla, M. et al. Resolving drug selection and migration in an inbred South American *Plasmodium falciparum* population with identity-by-descent analysis. PLoS Pathog. 18, e1010993 (2022).
    
    

39. 39.Auburn, S. & Barry, A. E. Dissecting malaria biology and epidemiology using population genetics and genomics. Int. J. Parasitol. 47, 77–85 (2017).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ijpara.2016.08.006&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27825828&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F06%2F2024.08.05.24311512.atom) 

40. 40.Bødker, R., Kisinza, W., Malima, R., Msangeni, H. & Lindsay, S. Resurgence of Malaria in the Usambara Mountains, Tanzania, An Epidemic of Drug-Resistant Parasites. Global Change and Human Health 1, 134–153 (2000).
    
    

41. 41.Verity, R. et al. The impact of antimalarial resistance on the genetic structure of *Plasmodium falciparum* in the DRC. Nat. Commun. 11, 2107 (2020).
    
    

42. 42.Amambua-Ngwa, A. et al. Chloroquine resistance evolution in *Plasmodium falciparum* is mediated by the putative amino acid transporter AAT1. Nat Microbiol 8, 1213–1226 (2023).
    
    

43. 43.Wamae, K. et al. No evidence of *P. falciparum* K13 artemisinin conferring mutations over a 24-year analysis in Coastal Kenya, but a near complete reversion to chloroquine wild type parasites. Antimicrob. Agents Chemother. 63, (2019).
    
    

44. 44.Lu, F. et al. Return of chloroquine sensitivity to Africa? Surveillance of African *Plasmodium falciparum* chloroquine resistance through malaria imported to China. Parasit. Vectors 10, 355 (2017).
    
    

45. 45.Alifrangis, M. et al. Independent Origin of *Plasmodium falciparum*Antifolate Super-Resistance, Uganda, Tanzania, and Ethiopia. Emerging Infectious Diseases vol. 20 1280– 1286 Preprint at doi:10.3201/eid2008.131897 (2014).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3201/eid2008.131897&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25061906&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F06%2F2024.08.05.24311512.atom) 

46. 46.Siame, M. N. P., Mharakurwa, S., Chipeta, J., Thuma, P. & Michelo, C. High prevalence of dhfr and dhps molecular markers in *Plasmodium falciparum* in pregnant women of Nchelenge district, Northern Zambia. Malar. J. 14, 190 (2015).
    
    

47. 47.Naidoo, I. & Roper, C. Mapping ‘partially resistant’, ‘fully resistant’, and ‘super resistant’ malaria. Trends Parasitol. 29, 505–515 (2013).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.pt.2013.08.002&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24028889&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F06%2F2024.08.05.24311512.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000325589300006&link_type=ISI) 

48. 48.Chaturvedi, R. et al. Geographical spread and structural basis of sulfadoxine-pyrimethamine drug-resistant malaria parasites. Int. J. Parasitol. 51, 505–525 (2021).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ijpara.2020.12.011&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=33775670&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F06%2F2024.08.05.24311512.atom) 

49. 49.Jiang, T. et al. Molecular surveillance of anti-malarial resistance Pfdhfr and Pfdhps polymorphisms in African and Southeast Asia Plasmodium falciparum imported parasites to Wuhan, China. Malaria Journal vol. 19 Preprint at doi:10.1186/s12936-020-03509-w (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s12936-020-03509-w&link_type=DOI) 

50. 50. Nkoli Mandoko, P., et al. Prevalence of *Plasmodium falciparum* parasites resistant to sulfadoxine/pyrimethamine in the Democratic Republic of the Congo: emergence of highly resistant pfdhfr/pfdhps alleles. J. Antimicrob. Chemother. 73, 2704–2715 (2018).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/jac/dky258&link_type=DOI) 

51. 51.Sipilanyambe, N. et al. From chloroquine to artemether-lumefantrine: the process of drug policy change in Zambia. Malar. J. 7, 25 (2008).
    
    

52. 52.Uwimana, A. et al. Association of *Plasmodium falciparum* kelch13 R561H genotypes with delayed parasite clearance in Rwanda: an open-label, single-arm, multicentre, therapeutic efficacy study. Lancet Infect. Dis. 21, 1120–1128 (2021).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1473-3099(21)00142-0&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=33864801&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F06%2F2024.08.05.24311512.atom) 

53. 53.Juliano, J. J. et al. Country wide surveillance reveals prevalent artemisinin partial resistance mutations with evidence for multiple origins and expansion of high level sulfadoxine-pyrimethamine resistance mutations in northwest Tanzania. medRxiv (2023) doi:10.1101/2023.11.07.23298207.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMy4xMS4wNy4yMzI5ODIwN3YyIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjQvMDgvMDYvMjAyNC4wOC4wNS4yNDMxMTUxMi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

54. 54.Young, N. W. et al. High frequency of artemisinin partial resistance mutations in the great lake region revealed through rapid pooled deep sequencing. medRxiv (2024) doi:10.1101/2024.04.29.24306442.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyNC4wNC4yOS4yNDMwNjQ0MnYxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjQvMDgvMDYvMjAyNC4wOC4wNS4yNDMxMTUxMi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

55. 55.Miotto, O. et al. Genetic architecture of artemisinin-resistant *Plasmodium falciparum*. Nat. Genet. 47, 226–234 (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3189&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25599401&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F06%2F2024.08.05.24311512.atom) 

56. 56.Amambua-Ngwa, A. et al. Major subpopulations of *Plasmodium falciparum* in sub-Saharan Africa. Science 365, 813–816 (2019).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIzNjUvNjQ1NS84MTMiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyNC8wOC8wNi8yMDI0LjA4LjA1LjI0MzExNTEyLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

57. 57.Fola, A. A. et al. Temporal genomic analysis of *Plasmodium falciparum* reveals increased prevalence of mutations associated with delayed clearance following treatment with artemisinin-lumefantrine in Choma District, Southern Province, Zambia. medRxiv (2024) doi:10.1101/2024.06.05.24308497.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyNC4wNi4wNS4yNDMwODQ5N3YxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjQvMDgvMDYvMjAyNC4wOC4wNS4yNDMxMTUxMi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

58. 58.Zambia National Malaria Indicator Survey (MIS) 2018. [https://www.path.org/resources/zambia-natl-malaria-indicator-survey-mis-2018/](https://www.path.org/resources/zambia-natl-malaria-indicator-survey-mis-2018/).
    
    

59. 59. A. Fola,  A., Dorman,  J., Levy,  M., Ciubotariu,  I. & Carpi, G. Optimized HT gDNA extraction from dried blood spot using QIAcube HT for malaria genomic epidemiology studies v1. (2020) doi:10.17504/protocols.io.bh69j9h6.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.17504/protocols.io.bh69j9h6&link_type=DOI) 

60. 60.Carpi, G. et al. Whole genome capture of vector-borne pathogens from mixed DNA samples: a case study of Borrelia burgdorferi. BMC Genomics 16, 434 (2015).
    
    

61. 61.Carpi, G., Gorenstein, L., Harkins, T. T., Samadi, M. & Vats, P. A GPU-accelerated compute framework for pathogen genomic variant identification to aid genomic epidemiology of infectious disease: a malaria case study. Brief. Bioinform. 23, (2022).
    
    

62. 62.Li, H. Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. arXiv [q-bio.GN*]* (2013).
    
    

63. 63.Conway, J. R., Lex, A. & Gehlenborg, N. UpSetR: an R package for the visualization of intersecting sets and their properties. Bioinformatics 33, 2938–2940 (2017).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btx364&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28645171&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F06%2F2024.08.05.24311512.atom) 

64. 64.Pebesma, E. Simple features for R: Standardized support for spatial vector data. R J. 10, 439 (2018).