- Open Access
Surveillance for sulfadoxine-pyrimethamine resistant malaria parasites in the Lake and Southern Zones, Tanzania, using pooling and next-generation sequencing
Malaria Journal volume 16, Article number: 236 (2017)
Malaria in pregnancy (MiP) remains a major public health challenge in areas of high malaria transmission. Intermittent preventive treatment in pregnancy (IPTp) with sulfadoxine-pyrimethamine (SP) is recommended to prevent the adverse consequences of MiP. The effectiveness of SP for IPTp may be reduced in areas where the dhps581 mutation (a key marker of high level SP resistance) is found; this mutation was previously reported to be common in the Tanga Region of northern Tanzania, but there are limited data from other areas. The frequency of molecular markers of SP resistance was investigated in malaria parasites from febrile patients at health centres (HC) in seven regions comprising the Lake and Southern Zones of mainland Tanzania as part of the ongoing efforts to generate national-wide data of SP resistance.
A cross-sectional survey was conducted in the outpatient departments of 14 HCs in seven regions from April to June, 2015. 1750 dried blood spot (DBS) samples were collected (117 to 160 per facility) from consenting patients with positive rapid diagnostic tests for malaria, and no recent (within past 2 months) exposure to SP or related drugs. DNA was extracted from the DBS, pooled by HC, and underwent pooled targeted amplicon deep sequencing to yield estimates of mutated parasite allele frequency at each locus of interest.
The dhps540 mutation was common across all 14 sites, ranging from 55 to 98.4% of sequences obtained. Frequency of the dhps581 mutation ranged from 0 to 2.4%, except at Kayanga HC (Kagera Region, Lake Zone) where 24.9% of sequences obtained were mutated. The dhfr164 mutation was detected only at Kanyanga HC (0.06%).
By pooling DNA extracts, the allele frequency of mutations in 14 sites could be directly determined on a single deep-sequencing run. The dhps540 mutant was very common at all locations. Surprisingly, the dhps581 was common at one health center, but rare in all the others, suggesting that there is geographic micro-heterogeneity in mutant distribution and that accurate surveillance requires inclusion of multiple sites. A better understanding of the effect of the dhps581 mutant on the efficacy of IPTp-SP is needed.
Malaria in pregnancy (MiP) remains a major public health challenge in areas of high malaria transmission. In pregnant women, malaria can cause mild to severe maternal anaemia, particularly in primigravid women, and placental infection can interfere with the maternal-fetal exchange of nutrients and oxygen, leading to preterm delivery and low birth weight, and consequently increasing neonatal mortality [1, 2]. In order to prevent the adverse consequences of MiP, the World Health Organization (WHO) recommends the use of intermittent preventative treatment in pregnancy (IPTp) with sulfadoxine-pyrimethamine (SP)—a full treatment dose administered to pregnant women during routine ANC visits in the 2nd and 3rd trimesters. IPTp has a protective efficacy of nearly 25% against low birth weight and 21% against neonatal mortality .
Plasmodium falciparum resistance to SP results from an ordered accumulation of mutations in two genes, namely P. falciparum dihydrofolate reductase (Pfdhfr) and P. falciparum dihydropteroate synthase (Pfdhps) that code for enzymes targeted by sulfadoxine and pyrimethamine, respectively. Resistance increases with the number of mutant alleles. The presence of five mutant alleles, the dhfr/dhps “quintuple mutant”, including the dhfr substitutions N51I, C59R, and S108N and the dhps substitutions A437G and K540E, are associated with a very high rate of failure when SP is used for the treatment of uncomplicated falciparum malaria . However, SP remains effective for IPTp, even where the prevalence of the quintuple mutant is high, and thus its use continues to be recommended by the WHO [5, 6]. The efficacy of SP for IPTp appears to be compromised in the presence of a sixth mutation at dhps581 (sextuple mutant; A581G); infection with parasites harbouring the sextuple mutant has been associated with increased placental parasitaemia and inflammation and failure of IPTp-SP to improve birth weight [7,8,9].
Previous studies have reported that the dhps 581 mutation is common in Tanga Region of Tanzania, with up to 57% of parasites harbouring the sextuple mutant ; however, there are limited data from other areas of Tanzania. In light of the possibility that IPTp-SP may provide no benefit to women in areas with a high prevalence of the sextuple mutant, it is critical to define the extent of this mutation. This study aimed to investigate the prevalence of the SP resistance mutations dhps K540E (a surrogate marker for quintuple mutant ) and dhps A581G  and dhfr I164L  (surrogate markers for sextuple mutant parasites, with A581G being the more widely reported) in the parasite population in the Lake Zone (Mwanza, Geita, Mara, and Kagera Regions) and Southern Zone (Lindi, Mtwara, and Ruvuma Regions) of Tanzania.
From April to June, 2015, a cross-sectional survey was conducted in 14 health facilities (HFs) in seven regions of mainland Tanzania: all four regions in the Lake Zone (Mwanza, Geita, Mara, and Kagera) and all three regions in the Southern Zone (Lindi, Mtwara, and Ruvuma Regions) (Fig. 1). These regions were selected as they had high prevalence of malaria in the 2011/12 Tanzania HIV and Malaria Indicator Survey (THMIS) . In each region, two health facilities were conveniently sampled from facilities with laboratory and where staff had been previously trained on collection of malaria rapid diagnostic tests (RDTs) for other projects to determine whether the dhpsA581G mutation was present in these regions.
It was estimated that a sample size of 100 malaria positive patients from each facility would produce a two-sided 95% confidence interval of 1.6–11.3%, assuming the true proportion of patients infected with a given mutation was 5%. To account for false positive results that may arise by use of RDTs, the target sample size was inflated by 20% to approximately 120 samples per facility.
All patients age ≥6 months presenting to the outpatient department who were diagnosed with malaria by RDT (Malaria Ag Pf/Pan, SD Bioline®) or blood smear and no self-reported history of recent exposure to SP or trimethoprim-sulfamethoxazole (within the past 2 months if not pregnant, or during the pregnancy, if pregnant) were consented. In each health facility, consenting eligible patients were consecutively enrolled until the required number of DBS samples (120) were obtained.
Sample collection procedures
Two blood spots were collected on filter paper (Whatman no. 3, GE Healthcare Life Sciences) at the time the blood sample was collected for malaria testing; if the RDT was negative or the patient did not consent or otherwise did not meet inclusion criteria, the filter paper was destroyed (Fig. 2). Blood spots from consenting individuals were air dried for 3–4 h on a flat, nonabsorbent surface, then sealed in plastic bags with a desiccant, and stored at room temperature at the study health facility. The samples from all 14 health facilities were packaged in a water-proof container, and shipped to the reference laboratory at the University of North Carolina (UNC) for pooled deep sequencing analyses at three genetic loci associated with sulfadoxine-pyrimethamine (SP) resistance, dhfr164, dhps540, and dhps581 [8, 12,13,14].
Optimizing DNA extraction
120 Dried Blood Spots (DBS) from the Kharumwa health center were used in a preliminary experiment to determine if DNA from batched DBS could be extracted reliably. Each subject’s DBS was punched 3 times. One of these punches was added to a single well of a 96-well plate. The other punches were pooled so that wells would contain either three or ten punches from individuals. The plate was extracted using the Chelex method using overnight incubation at 4 °C on a shaking platform at 700 RPM (TPM2 shaker, Sarstedt, Nümbrecht, Germany) . After extraction, 3, 9 and 30 ml aliquots of supernatant were removed from the n = 1, n = 3, and n = 10 extracts, respectively, and combined with aliquots from similarly extracted samples: 1×: aliquots of all 120 individually extracted samples; 3×: aliquots of all 40 of the 3-punch extracts; and 10× aliquots of all 12 of the 10-punch extracts.
Nested PCR amplification was done on each of the three DNA pools for dhps and dhfr using primers shown in Table 1 based on a previous report . PCR conditions were 95 °C × 2 min, followed by 40 rounds of [95 °C × 30 s → 45 °C × 30 s → 72 °C × 1 min] and a final round of 72 °C for 10 m. Individual indexed libraries for each gene were made for each DBS pooling condition using the NEB Next Ultra DNA Library Prep Kit (New England Biolabs, Ipswich, USA) and pooled in equimolar fashion. Pooled libraries were sequenced on an Illumina MiSeq using 2 × 300 bp reads at the UNC High Throughput Sequencing Facility. Allele frequencies for single nucleotide polymorphisms (SNPs) within the demultiplexed libraries were determined using methods previously described [17, 18].
Sequencing reads were aligned using bowtie2 with default settings to the P. falciparum 3D7 reference strain and the minor allele frequency (MAF) at all positions across the amplicon was determined using a revised version of our previously described Minor Allele Catcher (MAC) . Bases and reads were excluded from analysis if they had a mapping quality less than or equal to 10, a read length less than or equal to 200 bp, a base quality of less than or equal to 20, a base depth of less than or equal to 5000, if SNPs occurred in the first or last twenty-five bases of a read, or if 30% or fewer of the reads were on the forward or reverse strand.
Extraction and analysis of DBS from 13 other health centres
Based on the results of this pilot (Table 2), DNA from all of the remaining samples was extracted in batches of 10 by health centre. The DNA extracts from each of the health centres were then pooled yielding one amplicon pool for each site (13 pools). Library preps and analysis were performed as above. Data from the Kharumwa health centre, 10-punch extract, was re-run a second time to generate the 14-health center dataset.
A total of 3575 febrile patients were screened, and 1770 were found to be positive for malaria, yielding an overall positivity rate of 49.5% (Table 2). In this pilot work to determine the extraction efficiency of individual versus pooled extraction using the 120 participants from Kharumwa health centre, the proportion of mutant alleles varied only slightly as a result of batching DBS prior to extraction, with 9.4–10.7% alleles mutated at dhps540, 2.8–3.9% at dhps581, and none at dhfr164 (Table 3). As these differences in allele frequencies were considered negligible, DBS from all other health centers were extracted in pools of 10.
A total of 1750 dried blood spot (DBS) samples were collected (117–160 samples per facility) and used in the final study analysis. Deep sequencing resulted in an average of 341,420 reads at dhfr (range 208,769–644,927) and 168,260 reads at dhps (range 49,419–313,645). The library from Kharumwa was re-sequenced during this project. Based upon this repeat sequencing (Table 4) and the 3 sequences in Table 3, the allele frequencies ±standard deviations for dhps540 and dhps581 mutants were 90.0 ± 2.1% and 3.2 ± 0.7%, respectively. The dhps540 mutation was common across all 14 sites, with allele frequency ranging from 55 to 98.4%, with higher allele frequency at sites in Lake Zone compared to Southern Zone (Table 4). Frequency of the dhps581 mutation ranged from 0 to 2.4%, with the exception of Kayanga health centre (Kagera Region, Lake Zone) where 24.9% of sequences were mutated (Fig. 1). The dhfr164 mutation was detected only at Kanyanga health centre (0.6%).
In a survey of parasites from patients at 14 health facilities in the Lake and Southern Zones of Tanzania, the dhpsK540E mutation was very common [in 10 facilities (71%) more than 85% of alleles at dhpsK540E were mutated, and nowhere were fewer than 55% mutated], while the dhpsA581G mutation remained rare and focal, with frequency greater than 2.4% in only one facility, Kayanga Health Facility, where nearly 25% of alleles carried the dhpsA581G mutation. The quintuple mutation (represented by the dhpsK540E) has been reported widely across Tanzania, with prevalence ranging from 64 to 98% in one recent report assessing seven regions  and 77 to 95% in another . Previous reports from Tanga Region have found a high prevalence of parasites harbouring the dhpsA581G mutation (44% in Korogwe , 51% in Muheza , and 57% in Bondo ). Although the Lake and Southern Zones of Tanzania are not immediately adjacent to Tanga Region, the absence of the dhpsA581G mutation in the majority of sites highlights that there may be considerable geographic micro-heterogeneity. This is supported by data from Kavishe et al., who similarly found the dhpsA581G to be present at high proportion in only Tanga and Kagera Regions, though the prevalence of mutants at their sites were higher than the allele frequencies reported here . This is a reassuring finding, suggesting that IPTp-SP retains efficacy in the majority of Tanzania, but highlighting the need for monitoring in multiple geographic sites.
The pooling and sequencing methodology presented here is a cost-effective alternative to individual allele-specific PCR. Using this 2-step pooling methodology, it was possible to perform 90% fewer DNA extractions and 99% fewer sequencing runs than if samples from each participant had been sequenced individually, saving a substantial amount of both time and money. The individual extraction and bi-directional Sanger sequencing of 1750 samples would cost on the order of $19,000. Here the same data has been compiled for under $4000.
Another advantage of the pooling method is that it allows direct calculation of the allele frequency, rather than prevalence. That is, it calculates the percentage of the parasite population bearing the mutation rather than the prevalence of individuals bearing parasites with mutated alleles. From an evolutionary point of view, the allele frequency of a mutation is more important than prevalence of the mutation because it approximates the likelihood that a mosquito will become infected with a mutant parasite after exposure to a given population.
Allele frequency and prevalence can be different when mixed infections are present and when individuals have infections with varying levels of parasitaemia. Traditionally, prevalence data have to be transformed by a complex equation to yield predicted allele frequencies . However, in general, the allele frequency approximates the prevalence of major strain in a human population, barring any biases [16, 21].
This study has a number of limitations. Samples were collected only from the Lake and Southern Zones, neither of which are immediately adjacent to Tanga region, where the highest prevalence of the dhpsA581G have been reported, thus more studies in those areas are needed to better define the extent of the mutant. With regard to the pooling methodology employed here, while much more cost-effective and time-saving than traditional PCR, it is not possible to trace a parasite strain back to a participant. Also, some haplotypes tend to be amplified better than others, leading to PCR amplification bias. However, newer barcoding methods, such as primer ID , compensate for amplification bias and could allow identification of individuals (although DNA would still have to be extracted from individual DBS).
While overall, these data are reassuring with respect to the efficacy of IPTp-SP, there are select areas with a high prevalence of the sextuple mutant, where IPTp-SP may no longer provide a useful benefit against malaria. Spread of this sextuple mutant will threaten the usefulness of SP for IPTp. Given the fact that the quintuple mutant is already found throughout Tanzania, and the sextuple mutant has been found in high prevalence in several sites, continued surveillance in multiple sites, particularly in and around Tanga and Kagera, is warranted to monitor for the spread of the sextuple mutant. The pooling technique presented here provides a highly efficient and cost effective means to screen many samples from multiple sites.
Although the quintuple mutant was very common at all sampled facilities, dhpsA581G remains geographically restricted, suggesting that IPTp-SP remains effective in most of Tanzania. However, additional surveillance, particularly in and around Tanga and Kagera Region is warranted. This can be achieved efficiently through sequencing pooled samples as described here. Finally, a better understanding of the effect of the dhps581 mutant on the efficacy of IPTp-SP is needed.
dried blood spot on filter paper
intermittent preventive treatment in pregnancy
indoor residual spray
malaria focal person
malaria in pregnancy
malaria rapid diagnostic test
World Health Organization
Rogerson SJ, Mwapasa V, Meshnick SR. Malaria in pregnancy: linking immunity and pathogenesis to prevention. Am J Trop Med Hyg. 2007;77:14–22.
Desai M, ter Kuile FO, Nosten F, McGready R, Asamoa K, Brabin B, et al. Epidemiology and burden of malaria in pregnancy. Lancet Infect Dis. 2007;7:93–104.
Eisele TP, Larsen DA, Anglewicz PA, Keating J, Yukich J, Bennett A, et al. Malaria prevention in pregnancy, birthweight, and neonatal mortality: a meta-analysis of 32 national cross-sectional datasets in Africa. Lancet Infect Dis. 2012;12:942–9.
Kublin JG, Dzinjalamala FK, Kamwendo DD, Malkin EM, Cortese JF, Martino LM, et al. Molecular markers for failure of sulfadoxine-pyrimethamine and chlorproguanil-dapsone treatment of Plasmodium falciparum malaria. J Infect Dis. 2002;185:380–8.
Gutman J, Mwandama D, Wiegand RE, Ali D, Mathanga DP, Skarbinski J. Effectiveness of intermittent preventive treatment with sulfadoxine-pyrimethamine during pregnancy on maternal and birth outcomes in Machinga district, Malawi. J Infect Dis. 2013;208:907–16.
Kayentao K, Garner P, van Eijk AM, Naidoo I, Roper C, Mulokozi A, et al. Intermittent preventive therapy for malaria during pregnancy using 2 vs 3 or more doses of sulfadoxine-pyrimethamine and risk of low birth weight in Africa: systematic review and meta-analysis. JAMA. 2013;309:594–604.
Minja DTR, Schmiegelow C, Mmbando B, Boström S, Oesterholt M, Magistrado P, et al. Plasmodium falciparum mutant haplotype infection during pregnancy associated with reduced birthweight, Tanzania. Emerg Infect Dis. 2013;19:1446–54.
Harrington WE, Mutabingwa TK, Muehlenbachs A, Sorensen B, Bolla MC, Fried M, et al. Competitive facilitation of drug-resistant Plasmodium falciparum malaria parasites in pregnant women who receive preventive treatment. Proc Natl Acad Sci USA. 2009;106:9027–32.
Gutman J, Kalilani L, Taylor S, Zhou Z, Wiegand RE, Thwai KL, 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. 2015;211:1997–2005.
Kavishe RA, Kaaya RD, Nag S, Krogsgaard C, Notland JG, Kavishe AA, et al. Molecular monitoring of Plasmodium falciparum super-resistance to sulfadoxine–pyrimethamine in Tanzania. Malar J. 2016;15:335.
Tanzania Commission for AIDS (TACAIDS), Zanzibar AIDS Commission (ZAC), National Bureau of, Statistics (NBS), Office of the Chief Government Statistician (OCGS), and ICF International. Tanzania HIV/AIDS and Malaria Indicator Survey 2011–12 [Internet]. Dar es Salaam, Tanzania: TACAIDS, ZAC, NBS, OCGS, and ICF International; 2013. http://www.nbs.go.tz/index.php?option=com_content&view=article&id=353:2011-12-tanzania-hivaids-and-malaria-indicator-survey-thmis-report&catid=55:this&Itemid=145.
Plowe CV, Kublin JG, Doumbo OK. P. falciparum dihydrofolate reductase and dihydropteroate synthase mutations: epidemiology and role in clinical resistance to antifolates. Drug Resist Updates. 1998;1:389–96.
Gesase S, Gosling RD, Hashim R, Ord R, Naidoo I, Madebe R, et al. High resistance of Plasmodium falciparum to sulphadoxine/pyrimethamine in northern Tanzania and the emergence of dhps resistance mutation at Codon 581. PLoS ONE. 2009;4:e4569.
Naidoo I, Roper C. Mapping, “partially resistant”, “fully resistant”, and “super resistant” malaria. Trends Parasitol. 2013;29:505–15.
Taylor SM, Antonia AL, Harrington WE, Goheen MM, Mwapasa V, Chaluluka E, et al. Independent lineages of highly sulfadoxine-resistant Plasmodium falciparum haplotypes, eastern Africa. Emerg Infect Dis. 2014;20:1140–8.
Taylor SM, Parobek CM, Aragam N, Ngasala BE, Mårtensson A, Meshnick SR, et al. Pooled deep sequencing of Plasmodium falciparum isolates: an efficient and scalable tool to quantify prevailing malaria drug-resistance genotypes. J Infect Dis. 2013;208:1998–2006.
Taylor SM, Parobek CM, DeConti DK, Kayentao K, Coulibaly SO, Greenwood BM, et al. Absence of putative artemisinin resistance mutations among Plasmodium falciparum in sub-Saharan Africa: a molecular epidemiologic study. J Infect Dis. 2015;211:680–8.
Brazeau NF, Hathaway N, Parobek CM, Lin JT, Bailey JA, Lon C, et al. Longitudinal pooled deep sequencing of the Plasmodium vivax K12 Kelch gene in Cambodia reveals a lack of selection by artemisinin. Am J Trop Med Hyg. 2016;95:1409–12.
Baraka V, Ishengoma DS, Fransis F, Minja DTR, Madebe RA, Ngatunga D, et al. High-level Plasmodium falciparum sulfadoxine-pyrimethamine resistance with the concomitant occurrence of septuple haplotype in Tanzania. Malar J. 2015;14:439.
Hastings IM, Smith TA. MalHaploFreq: a computer programme for estimating malaria haplotype frequencies from blood samples. Malar J. 2008;7:130.
Aragam NR, Thayer KM, Nge N, Hoffman I, Martinson F, Kamwendo D, et al. Diversity of T cell epitopes in Plasmodium falciparum circumsporozoite protein likely due to protein-protein interactions. PLoS ONE. 2013;8:e62427.
Jabara CB, Jones CD, Roach J, Anderson JA, Swanstrom R. Accurate sampling and deep sequencing of the HIV-1 protease gene using a primer ID. Proc Natl Acad Sci USA. 2011;108:20166–71.
JG, JMN, SM, OMM, DSI, RAW designed the study. DSI, SM, OMM, RAW, SL, NK, CK, and LAP oversaw the collection of the samples. SMD, KLT, NJH, JAB, JJJ, and SRM analysed the samples. JG, JMN, and SRM drafted the manuscript which all authors edited and approved. All authors read and approved the final manuscript.
We are grateful to the staff of the HCs, NIMR Tanga team (Ezekiel Malecela and August Nyaki), and RTI staff who implemented the study and district malaria focal persons (DMFP) in Lake and Southern Zones for supervising the study.
The authors declare that they have no competing interests.
Availability of data and materials
The datasets generated during the current study are available from the corresponding author on reasonable request.
The findings and conclusions presented in this manuscript are those of the authors and do not necessarily reflect the official position of the US President’s Malaria Initiative, US Agency for International Development, or US Centers for Disease Control and Prevention.
Ethical approval and consent to participate
The study protocol was approved by the Medical Research Coordinating Committee (MRCC) of the Tanzania National Health Institute for Medical Research (NIMR) [reference NIMR/HQ/R.8a/Vol.IX/1924]. The US Centers for Disease Control and Prevention (CDC) human research protection office reviewed the protocol and granted a non-engaged determination for CDC’s role in this project. Informed consent was obtained from all patients/parents or guardians of patients (for minors) participating in the study.
This study was made possible through support provided to the Tanzania Vector Control Scale-up Project (Cooperative Agreement 621-A-00-10-00015-00) by the President’s Malaria Initiative (PMI) via the US Agency for International Development (USAID) and an inter-agency agreement with the US Centers for Disease Control and Prevention (CDC).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Ngondi, J.M., Ishengoma, D.S., Doctor, S.M. et al. Surveillance for sulfadoxine-pyrimethamine resistant malaria parasites in the Lake and Southern Zones, Tanzania, using pooling and next-generation sequencing. Malar J 16, 236 (2017) doi:10.1186/s12936-017-1886-9
- Plasmodium falciparum