Sulfadoxine-pyrimethamine (SP) is used for prevention of malaria in pregnant women in Angola. We sequenced the Plasmodium falciparum dihydrofolate reductase (pfdhfr) and dihydropteroate synthase (pfdhps) genes, implicated in SP resistance, in samples collected during a 2019 study of artemisinin-based combination therapy efficacy in Benguela, Lunda Sul, and Zaire provinces. A total of 90 day 0 and day of failure samples were individually sequenced, while 508 day 0 samples from participants without recurrent parasitemia were pooled after DNA extraction into 61 pools. The N51I, C59R, and S108N pfdhfr mutations and A437G pfdhps mutations were present at high proportions in all provinces (weighted allele frequencies, 62% to 100%). The K540E pfdhps mutation was present at lower proportions (10% to 14%). The A581G pfdhps mutation was only observed in Zaire, at a 4.6% estimated prevalence. The I431V and A613S mutations were also only observed in Zaire, at a prevalence of 2.8% to 2.9%. The most common (27% to 66%) reconstructed haplotype in all three provinces was the canonical quadruple pfdhfr pfdhps mutant. The canonical quintuple mutant was absent in Lunda Sul and Benguela and present in 7.9% of samples in Zaire. A single canonical sextuple (2.6%) mutant was observed in Zaire Province. Proportions of the pfdhps K540E and A581G mutations were well below the World Health Organization thresholds for meaningful SP resistance (prevalence of 95% for K540E and 10% for A581G). Samples from therapeutic efficacy studies represent a convenient source of samples for monitoring SP resistance markers.

Sequencing large numbers of individual samples is often needed for countrywide antimalarial drug resistance surveillance. Pooling DNA from several individual samples is an alternative cost and time saving approach for providing allele frequency (AF) estimates at a population level. Using 100 individual patient DNA samples of dried blood spots from a 2017 nationwide drug resistance surveillance study in Haiti, we compared codon coverage of drug resistance-conferring mutations in four Plasmodium falciparum genes (crt, dhps, dhfr, and mdr1), for the same deep sequenced samples run individually and pooled. Samples with similar real-time PCR cycle threshold (Ct) values (+/- 1.0 Ct value) were combined with ten samples per pool. The sequencing success for samples in pools were higher at a lower parasite density than the individual samples sequence method. The median codon coverage for drug resistance-associated mutations in all four genes were greater than 3-fold higher in the pooled samples than in individual samples. The overall codon coverage distribution for pooled samples was wider than the individual samples. The sample pools with < 40 parasites/μL blood showed more discordance in AF calls for dhfr and mdr1 between the individual and pooled samples. This discordance in AF estimation may be due to low amounts of parasite DNA, which could lead to variable PCR amplification efficiencies. Grouping samples with an estimated ≥ 40 parasites/μL blood prior to pooling and deep sequencing yielded the expected population level AF. Pooling DNA samples based on estimates of > 40 parasites/μL prior to deep sequencing can be used for rapid genotyping of a large number of samples for these four genes and possibly other drug resistant markers in population-based studies. As Haiti is a low malaria transmission country with very few mixed infections and continued chloroquine sensitivity, the pooled sequencing approach can be used for routine national molecular surveillance of resistant parasites.