Skip directly to site content Skip directly to page options Skip directly to A-Z link Skip directly to A-Z link Skip directly to A-Z link
Volume 28, Number 4—April 2022

Molecular Surveillance for Imported Antimicrobial Resistant Plasmodium falciparum, Ontario, Canada

Ruwandi Kariyawasam, Rachel Lau, Eric Shao, Katherine Tan, Adrienne Showler, Filip Ralevski, Samir N. Patel, and Andrea K. BoggildComments to Author 
Author affiliations: University of Alberta, Edmonton, Alberta, Canada (R. Kariyawasam); Alberta Precision Labs-Public Health, Edmonton (R. Kariyawasam); Public Health Ontario Laboratory, Toronto, Ontario, Canada (R. Lau, F. Ralevski, S.N. Patel); University of Western Ontario, London, Ontario (E. Shao); University of Toronto, Toronto (K. Tan, A. Showler, S.N. Patel, A.K. Boggild); Georgetown University, Washington, DC, USA (A. Showler); Toronto General Hospital, Toronto (A. Showler, A.K. Boggild)

Cite This Article


Single-nucleotide polymorphisms at several loci have been correlated with Plasmodium falciparum drug resistance. We examined the prevalence of resistance markers in P. falciparum from imported malaria cases in Canada during 3 time periods, 2008–2009, 2013–2014, and 2017–2018. We evaluated single-nucleotide polymorphisms at atpase6 (pfATPase6), pfcrt (chloroquine resistance transporter), cytb (cytochrome b), dhfr (dihydrofolate reductase), dhps (dihydropteroate synthetase), mdr1 (multidrug resistance protein) and mdr1 copy number, and kelch13 (kelch protein gene on chromosome 13). Over time, we observed increasing mutant genotypes for dhfr S108N and dhps A613T and decreasing mutant genotypes for mdr1 N86Y, D1246Y, pfcrt K76T, and pfcrt 74–75; we identified no kelch13 mutations. We observed fewer mutations indicative of chloroquine resistance over time, which may reflect reduced chloroquine pressure in specimens from travelers to Africa. Mutations conferring proguanil resistance increased over time. Minor genotypes confirm the heterogeneous nature of infection and may affect treatment success.

Malaria remains the deadliest vectorborne infectious disease worldwide (1). Plasmodium spp., most commonly P. falciparum, are responsible for ≈229 million cases and ≈500,000 deaths from malaria annually (2). Although malaria incidence and death have decreased over the past decade, emerging antimalarial drug resistance, fueled by counterfeiting, overuse, and underdosing, threatens control and elimination efforts (2). P. falciparum resistance mutations have resulted in waning efficacy in multiple antimalarial classes including artemisinins, quinolines, and antifolates (3). Increased international travel and climate change are exacerbating the spread of malarial vectors (4), making drug-resistant P. falciparum malaria imported from endemic regions, in particular Africa and Southeast Asia, a growing concern (5). Surveillance of antimalarial drug resistance, particularly among P. falciparum specimens, is crucial both to inform universal treatment guidelines and to identify global patterns of emerging resistance. We aimed to identify the prevalence of several resistance markers, including genes that confer resistance to chloroquine, mefloquine, atovaquone/proguanil, and artemisinins, and to quantify the copy number of multidrug resistance genes (pfmdr1) in P. falciparum isolates from malaria cases imported to Ontario, Canada, over a 10-year period.



From the malaria biobank at the Public Health Ontario Laboratory in Toronto, Ontario, Canada, we retrieved unique surplus whole-blood clinical specimens containing P. falciparum from the years 2008–2009, 2013–2014, and 2017–2018. We confirmed the P. falciparum specimens from the biobank as monoinfections using multiplex real-time quantitative PCR (qPCR) (1,6,7) after clinical testing, which included Giemsa-stained thick and thin blood film examination by certified medical lab technologists and rapid diagnostic test using the Abbott BinaxNOW malaria kit ( To eliminate the possibility of recrudescent infections, we excluded specimens from patients who had positive blood smears in the year before our enrollment period.

Extracting DNA from Whole Blood and Detecting Parasite DNA by Real-Time qPCR

We extracted DNA from surplus blood (0.4 mL) using the KingFisher Pure DNA Blood Kit (Thermo Fisher Scientific, and stored it in 50-μL aliquots: one at −20°C for working stock and another at −80°C for storage stock. We used included BEI Resources W2 MRA-157 ( and ATCC 3D7 ( samples as P. falciparum–positive controls. We tested samples using 2 species-specific duplex real-time qPCRs (P. falciparum and P. vivax; P. malariae and P. ovale) (6,7). Samples were 25 μL total volume, including 12.5 μL Taqman Universal MasterMix (Thermo Fisher) and 5 μL template DNA. To perform all reactions, we used the Applied Biosystems 7500 Fast Real-Time PCR system (Thermo Fisher) for 2 min at 50°C and 10 min at 95°C, followed by 45 cycles of 15 s at 95°C and 1 min at 60°C.


We analyzed 23 different single-nucleotide polymorphisms (SNPs) across 6 genes: atpase6 (PfATPase6; SNPs A623E, S769N), cytb (cytochrome b; SNP Y268NSC), dhfr (dihydrofolate reductase; SNPs A16V, C50R, N51I, C59R, S108N, I164L), dhps (dihydropteroate synthetase; SNPs S436FA, A437G, K540E, A581G, A613TS), mdr1 (multidrug resistance protein; SNPs N86Y, Y184F, S1034CTR, N1042D, D1246Y), and pfcrt (chloroquine resistance transporter; SNPs K76T, N75E, M74I, C72S). We performed 16 PCR reactions according to the Pyromark PCR Kit protocol (QIAGEN, using 400 nmol/L–concentration primers published elsewhere (Table 1). We performed all reactions on ABI Veriti Thermal Cyclers (Thermo Fisher) with adjusted annealing temperatures of 44°C for dhfr 164 and 63.8°C for pfcrt. We ran 3 μL of PCR product on a 1% agarose gel at 100 V for 30 min to ensure band presence before pyrosequencing.

We performed pyrosequencing using the Pyromark Q24 Vacuum Workstation and Pyromark Q24 Pyrosequencer (QIAGEN) according to manufacturer protocols. We designed 21 pyrosequencing assays using PSQ assay design software version 1.0.6 ( (Table 1). To determine run-to-run consistency and reproducibility, we calculated the intra-assay coefficient of variation by running the cytb assay on 3 samples in quintuplicate. We analyzed samples using Pyromark Q24 version 1.0.10 software (QIAGEN).

Sanger Sequencing

We interrogated kelch13 (kelch protein gene on chromosome 13) for 20 SNPs according to Institut Pasteur protocols using primers published elsewhere (8,9) (Table 1). To validate the results obtained from pyrosequencing, we performed a quality control check of the gene targets using the first 10% of cases from each time period and controls. We designed sequencing primers with Primer 3 (10,11) and performed PCR using primers published elsewhere (8,9,12) that target drug resistance genes according to the Accuprime Pfx SuperMix protocol (Thermo Fisher) (Table 1); we adjusted the annealing temperature to 60°C and cycles to 40 for cytb and 45 for all mdr1 and dhfr assays. We ran 3 μL of PCR product on a 1% agarose gel at 100 V for 30 min to ensure band presence before sequencing.

We performed sequencing reactions using the Big Dye Terminator v3.1 Cycle Sequencing Kit (Thermo Fisher) with a volume of 20 μL, including 1 μL of PCR product, 2 μL of Big Dye ready reaction mix, 3 μL of Thermo Fisher sequencing buffer, and 2 μL of 10 μmol/L of primer. We used a Bio-Rad C1000 Thermal Cycler ( for 1 min at 96°C, 25 cycles of 10 s at 96°C, 5 s at 50°C, and 4 min at 60°C. We purified product using an in-house isopropanol protocol including incubating at room temperature using 15 μL of 100% isopropanol for 15 min, centrifuging at 3,000 × g for 30 min at room temperature, centrifuging at 2,000 × g for 10 min after 35 μL of 70% isopropanol had been added, then inverting the plate and centrifuging for 700 × g at room temperature for 1 min. Last, we added 20 μL of Hi-Di Formamide (Thermo Fisher) to reconstitute each sample and sequenced using an Applied Biosystems 3730xl DNA Analyzer (Thermo Fisher). We standardized data using a Thermo Fisher sequencing analysis program and used a BioEdit sequence alignment editor version 7.2.5 ( to analyze the sequence for nucleotide changes at each identified SNP location.

Pfmdr1 Copy Number Analysis

We tested specimens for pfmdr1 copy number using real-time qPCR with primers and probes published elsewhere (12), using a total volume of 25 μL, including 12.5 μL Taqman Universal MasterMix (Thermo Fisher) and 5 μL template DNA. We performed reactions on an Applied Biosystems 7500 Fast Real-Time PCR system (Thermo Fisher) using the same conditions used to detect parasite DNA. We ran the process in triplicate for specimens and used the 2-ΔΔCT method to calculate relative expression (13).

Statistical Analyses

We considered a genotype mutant dominant if the mutant allelic frequency within the specimen was >50% and wildtype dominant if <50%. We calculated median and range pfmdr1 relative copy numbers for each time period and compared all categorical variables (proportions of mutant versus wildtype genotypes) using the Fisher exact test. We calculated median and range for mutant allelic frequencies within specimens and presented them descriptively across time periods to inform surveillance for resistance mutations over time. We compared differences in median pfmdr1 relative copy numbers across time periods using the Kruskal-Wallis H test and set significance at p <0.05. We analyzed data using Stata version 13 software (StataCorp LLC,


Over the enrollment period, we retrieved a total of 574 primary blood specimens from the malaria biobank containing >1 Plasmodium spp. determined using real-time qPCR; 566 (98.6%) specimens contained just 1 species of Plasmodium and the other 8 (1.4%) contained 2 Plasmodium species. Of specimens with single-species infections, 365/574 (64%) were positive only for P. falciparum, 142/574 (25%) only for P. vivax, 46/574 (8%) only for P. ovale, and 13/574 (2.3%) only for P. malariae. Of the specimens with 2 Plasmodium species, 3/574 (0.5%) were P. falciparum–P. ovale co-infections and 5/574 (0.9%) were P. ovale–P. vivax co-infections. Only 243/365 (67%) specimens with P. falciparum monoinfection were identified as unique cases, 75 (31%) occurring during July 2008–June 2009, 79 (33%) during July 2013–June 2014, and 89 (37%) during July 2017–June 2018. Of those 243 P. falciparum–only specimens from unique cases, 169 (70%) were from male patients, 66 (27%) from female patients, and 8 (3%) had no reported sex; there was no difference in proportions of sex across time periods (p = 0.47) (Table 2). Mean age was 39.2 (range 3–88) years, which did not differ between time periods (p = 0.61) (Table 2). Mean parasitemia was 0.3% (<0.01%–24.0%) across all 3 time periods, and we saw no significant differences over time (p = 0.10) (Table 2). A total of 186 (77%) patients had documented travel history, 81 (33%) in West Africa, the most common region, and 40 (16%) in Nigeria, the most common country (Table 2). Only 5/243 (2%) unique P. falciparum specimens were from travelers to Southeast Asia, and 1 (0.4%) was from a traveler to the Caribbean (Table 2).

Prevalence of Mutations in Mutant Dominant Genotypes

All 243 unique P. falciparum­–only specimens contained >1 resistance mutation. We detected the highest prevalence of mutant genotypes in dhfr; 212 (91%) cases demonstrated triple codon mutation N51I, C59R, and S108N (Appendix Table 1). dhfr SNP analysis revealed an increase in S108N mutant genotypes from 89% in 2008–2009 to 100% in 2017–2018 (Appendix Table 1). pfcrt SNP analysis revealed a gradual decline in K76T mutant genotype from 57% in 2008–2009 to 33% in 2017–2018 (p = 0.011), coupled with a decrease in M74I and N75E from 52% in 2008–2009 to 27% in 2017–2018 (p = 0.010; Appendix Table 2). dhps A613T SNP analysis revealed an overall increase in mutant genotypes, from 12% in 2008–2009 to 28% in 2017–2018 (p = 0.024; Appendix Table 3). mdr1 N86Y and D1246Y analysis revealed a decline in mutant genotypes from 2008–2009 and 2017–2018 (p<0.003; Appendix Table 4). We observed no difference in mutant genotype populations across time periods for dhps K540E, A581G, A437G, S436A, S436F mdr1 Y184F, S1034T, or atpase6 A623E (Appendix Tables 3–5). Last, we observed no mutant genotype populations for cytb Y268NSC across time periods (Appendix Table 6).

Mutant Allelic Frequency within Wildtype Dominant Specimens

We presented summaries of mutant allelic frequencies within wild-type–dominant specimens (i.e., with <50% within-specimen mutation prevalence) for dhfr (Appendix Table 1), pfcrt (Appendix Table 2), dhps (Appendix Table 3), mdr1 and mdr1 copy number (Appendix Table 4), atpase6 (Appendix Table 5), and cytb (Appendix Table 6). Many minor changes in allelic frequencies of questionable importance occurred across time periods (Appendix Tables 1–6). Notable findings include that pfcrt K76T and C72S mutant allele populations among wild-type infections declined significantly from 2008–2009 to 2017–2018 (Appendix Table 2). Conversely, mutant allele frequency within mdr1 D1246Y increased from 7% in 2008–2009 to 16% in 2017–2018, despite a decline in overall proportion of mutant-dominant genotypes (i.e., with >50% within-specimen mutation frequency) from 17.6% to 3.5% across time periods (Appendix Table 4). Frequency of any allelic mutant among wildtype infections at cytb Y268N increased from 0% in 2008–09 and 2013–2014 to 3% in 2017–2018 (Appendix Table 6).

mdr1 Copy Number and kelch13 Mutations

Using the 3D7 comparator reference strain, which has a pfmdr1 copy number of 1, we found an increase in median pfmdr1 copy number from 1.1 (range 0.83–1.4) in 2008–2009 to 1.9 (range 0.73–5.4) in 2017–2018 (p<0.001) (Appendix Table 4). kelch13 SNP analysis of >20 codons revealed no mutations for all analyzed specimens, despite 5 specimens with confirmed P. falciparum being from Southeast Asia (data not shown). We observed 2 silent mutations at codon 553 and 561 in the 2017–2018 specimens (sequences deposited to GenBank under accession nos. OM489472 and OM489473); however, travel history was not recorded for those patients.


We analyzed 243 unique cases of P. falciparum imported to Canada using pyrosequencing assays based on previously reported genetic markers known to confer drug resistance. Our results have provided stakeholders with a resource for tracking antimalarial resistance over time. Conventional SNP analysis involves Sanger sequencing and PCR–restriction fragment-length polymorphism, but over the past decade pyrosequencing has added utility for detecting drug-resistance markers (14). Now, pyrosequencing has been partially replaced by more powerful techniques, such as targeted next-generation sequencing, which remain accessible mostly just in reference or research laboratories because of bioinformatics requirements. Our analysis of P. falciparum markers of genotypic resistance using pyrosequencing provides evidence of a quantitative advantage over Sanger sequencing because allelic frequencies within wild-type populations can be monitored over time to capture emerging molecular resistance. We observed significant decreases in mutant genotypes for chloroquine resistance genes and an increase over time in the proportion of mutant genotypes for dhfr, the gene conferring resistance to proguanil. We also observed significant increases in markers of antifolate drug resistance correlating with increases in mutant allelic frequencies within wild-type populations.

The emergence of antimicrobial-resistant P. falciparum began in 1957 when the first cases of resistance to chloroquine were observed along the Cambodia–Thailand border (15). pfcrt encodes for a transmembrane protein in the digestive vacuoles of P. falciparum parasites. Mutations of the pfcrt gene, specifically at positions 72, 74, 75, and 76, confer resistance to chloroquine (1618). pfmdr1 mutation N86Y combined with increased pfmdr1 copy numbers, both of which have been documented in many parts of Southeast Asia and recently in sub-Saharan Africa, is also implicated in chloroquine resistance (12,17,1921). Our data document a decline over time in pfcrt mutants from positions 72 and 76 in imported cases of P. falciparum, indicating molecular reversion in chloroquine resistance at these codons. This phenomenon may reflect the reduction in chloroquine drug pressure on specimens from West Africa across the 10-year study period, during which oral artemisinin combination therapy has supplanted chloroquine use.

Drugs targeting folate metabolism, including pyrimethamine, sulfadoxine, and proguanil, are common components of fixed-dose combination antimalarials and inhibit purine and pyrimidine biosynthetic pathways (3). Point mutations in the dhfr and dhps genes conferring sulfadoxine and pyrimethamine resistance were identified as early as 1988 in Kenya. dhps mutations S436AF, A437G, A581G, and A613TS have been identified in sulfonamide-resistant specimens primarily found in Southeast Asia (2123). Higher levels of resistance have been identified in specimens carrying multiple mutations conferring a synergistic effect, especially double mutation S436FA and A613TS (23). Overall, we documented an increase in all mutant dhps genotypes except K540E. Increases in mutant genotypes coincided with overuse of chloroquine alternatives, which have been heavily deployed in sub-Saharan Africa.

dhfr S108N is the essential driving force associated with resistance to pyrimethamine (3,24). Our data highlight evidence of antifolate resistance over time, which may relate to increasing reliance on antifolate drugs in the face of chloroquine resistance. Additional mutations including dhfr N51I, C59R, and I164 have been associated with more pyrimethamine-resistant specimens (3,24,25). We did not identify increases in either N51I or C59R mutations and detected only one I164L specimen in the 2008–2009 time frame, with no subsequent occurrences.

Our findings are limited by a small sample set and incomplete travel history, but this limitation may suggest a decreased prevalence of I164L mutation in sub-Saharan Africa. Of interest, the dhfr A16V mutation generally occurs in conjunction with S108N, which is known to cause cycloguanil resistance (26). However, our data demonstrate no such correlation in the specimens we analyzed. Analyses by many groups over the years have led to the conclusion that greater numbers of point mutations in dhfr confer greater resistance, and cross-resistance between pyrimethamine and cycloguanil is a possibility, whereas S108N is the necessary first mutation responsible for pyrimethamine resistance (3). Our results coincide with previous literature resulting in data on a wave of antifolate resistance markers in specimens of P. falciparum imported to Ontario.

Artemisinin derivatives, the most powerful antimalarials used clinically, have been the main focus of concern surrounding recent efforts to mitigate P. falciparum antimicrobial resistance. After the initial surge of chloroquine resistance, which was quickly followed by a diminution of antifolate effectiveness, the World Health Organization recommended use of artemisinin-based combination therapies as first-line drugs of choice. However, similar to the experience when generations of antimalarial drugs have been heavily used in other endemic areas, many countries in Africa have seen strains of multi­drug-resistant P. falciparum appeared within 1 year of introducing artemisinin-based combination therapy formulations, occurring as early as 2003 (27). In Southeast Asia, artemisinin resistance has been associated with mutations in the propeller region of P. falciparum kelch13 (28). At present, there is no documented large-scale spread of kelch13 mutations in Africa; however, kelch13 mutations have been documented in specimens from Bangladesh, suggesting a westward migration of mutations from their origins in the Thailand–Myanmar region (2935).

Recently, there have been reports of kelch13-independent treatment failure with artemether/lumefantrine in P. falciparum cases imported to the United Kingdom by persons with travel history to Angola, Liberia, and Uganda (30). Despite novel mutations, including G112E, outside of the propeller-encoding domain in some patients, no reduced artemisinin susceptibility has been reported in the Greater Mekong region of Southeast Asia (8,9,28,30,34). Our data from specimens predominantly from West Africa correlate with these and other findings corroborating an absence of detectable mutations in 20 SNPs identified in the kelch13 gene between codon 440 and the 3′ end of the coding region. On the other hand, pfATP6, which is also known as pfSERCA or pfATPase6, is a calcium ATPase gene involved in calcium ion transport. Point mutations at A623E and S769N in pfATP6 could be associated with artemisinin resistance. Increased prevalence of these point mutations may be responsible for the rise of resistance in areas such as Cambodia, French Guiana, and Senegal, where artemisinin use is uncontrolled (27). Our data suggest that artemisinin-resistant P. falciparum infection is unlikely to occur in travelers returning to Canada, given that only 1 atpase6 A623E mutant was detected in both the 2008–2009 and 2013–2014 time frames, and none for the atpase6 S769N and kelch13 mutations.

The fixed-dose combination atovaquone/proguanil, an antimalarial prophylactic popular among travelers, has been recommended more broadly in response to potential resistance to other antimicrobial classes, including artemisinins (36). However, documented cases with cytb mutations have resulted in treatment failure associated with resistance to atovaquone. We documented an exceptionally high prevalence of P. falciparum specimens (83%) with triple codon dhfr 51, 59, 108 mutation, which confer resistance to proguanil, leaving the partner drug, atovaquone, as the sole antimalarial effective in clearing P. falciparum (4 of every 5 cases in our setting). Any resistance to atovaquone, conferred by cytb268 mutations, or impaired atovaquone absorption can lead to treatment failure (37,38); resistance as a result of cytb Y268NSC has been identified in travelers returning from Kenya, Angola, and many parts of East Africa (22,3941). We observed no cytb-dominant mutants in our study, but in some specimens, the mutant allele frequency in the wild-type specimens were close to 20%, suggesting ≈1/5 of the infection could be comprised of mutant strains. Whether having a minor population of allelic mutants instead of dominant mutant allelic frequency equating to an overall mutant genotype can lead to treatment failure or recrudescence is unknown, but nevertheless, that condition may be an early indicator of emerging atovaquone resistance. Identifying subpopulations of resistant mutant strains can help determine appropriate treatment strategies and follow-up; prevalence of resistant strains as low as 10% can double in vitro 50% inhibitory concentration to chloroquine (42). This possibility reiterates the value of pyrosequencing and targeted next generation sequencing, which enable quantification of minor allelic populations.

Our conclusions may be limited by lack of complete data on travel history, history of malaria infection, and treatment strategies and outcomes among study participants. Subtle changes in mutant allelic frequency noted across time periods may be attributable to travel proclivities, period-to-period natural variation in expression, or a true change in prevalence. The technical limitations of developing a pyrosequencing assay targeting positions 74 and 75 only enabled qualitative analysis for pfcrt double mutations at those positions. A technical limitation of pyrosequencing includes diminished accuracy of determining proportions <5%. Although kelch13 remains the only gene target analyzed via Sanger sequencing, minor alleles in mixed populations cannot be detected by Sanger sequencing; thus, pyrosequencing is a more useful and reliable sequencing technology for SNP analysis. Although pyrosequencing has largely been replaced by more effective techniques, such as targeted next-generation sequencing and possibly whole-genome sequencing, those newer techniques are limited by cost and the requirement for bioinformatics infrastructure and expertise, which generally confine their application to large reference or research laboratories. Future studies involving prospective recruitment of P. falciparum–infected travelers, with complete demographic, clinical, and treatment outcome data and use of targeted next-generation sequencing to identify other SNPs, risk factors, or confounders of resistance, could fill these knowledge gaps.

Surveillance of SNPs can serve as a sentinel of clinical resistance. Triple codon dhfr mutants, for example, were detected years before observation of increased sulfadoxine/pyrimethamine treatment failures in Kenya (43). Robust surveillance analyses over time of known markers conferring drug resistance in P. falciparum imported by travelers are lacking. We observed a high percentage of molecular resistance to proguanil and none to atovaquone. We observed an increase in imported antifolate-resistant of P. falciparum strains, including a decrease in chloroquine-resistant strains and virtually no signs of artemisinin resistance because of mutations at kelch13. Importation of resistant P. falciparum mono-infections highlights the importance of developing better surveillance tools to monitor drug resistance patterns based on time and source region. Such surveillance would also inform clinical decision-making and serve as a resource for tracking antimalarial resistance that clinicians could monitor for possible genetic markers of reduced drug efficacy.

Dr. Kariyawasam is a clinical microbiology fellow at the University of Alberta, Edmonton, Alberta, Canada. She completed her PhD on American tegumentary leishmaniasis in the laboratory of Dr. Boggild and has a longstanding interest in other imported vectorborne diseases.



We thank MR4 for providing us with Plasmodium small subunit ribosomal RNA clones contributed by Peter A. Zimmerman.

This study was funded by the Project Initiation Fund of Public Health Ontario. A.K.B. is supported as a clinician scientist by the Departments of Medicine at the University Health Network and University of Toronto.

Author contributions: A.K.B. conceived the study; A.K.B., R.K. and R.L., contributed to study design, data collection, and analysis and interpretation, and were primarily responsible for writing the manuscript. E.S., K.T., F.R., and S.N.P. contributed to data collection and critical appraisal and revision of the manuscript. A.S. contributed to data analysis and interpretation and appraisal and revision of the manuscript.



  1. Phuong  M, Lau  R, Ralevski  F, Boggild  AK. Sequence-based optimization of a quantitative real-time PCR assay for detection of Plasmodium ovale and Plasmodium malariae. J Clin Microbiol. 2014;52:106873. DOIPubMedGoogle Scholar
  2. World Health Organization. Fact Sheet: World Malaria Report 2020 [cited 2022 Feb 27].
  3. Gregson  A, Plowe  CV. Mechanisms of resistance of malaria parasites to antifolates. Pharmacol Rev. 2005;57:11745. DOIPubMedGoogle Scholar
  4. Paul  REL, Diallo  M, Brey  PT. Mosquitoes and transmission of malaria parasites - not just vectors. Malar J. 2004;3:39. DOIPubMedGoogle Scholar
  5. McCarthy  AE, Morgan  C, Prematunge  C, Geduld  J. Severe malaria in Canada, 2001-2013. Malar J. 2015;14:151. DOIPubMedGoogle Scholar
  6. Phuong  M, Lau  R, Ralevski  F, Boggild  AK. Survival analysis of diagnostic assays in Plasmodium falciparum malaria. Malar J. 2015;14:350. DOIPubMedGoogle Scholar
  7. Shokoples  SE, Ndao  M, Kowalewska-Grochowska  K, Yanow  SK. Multiplexed real-time PCR assay for discrimination of Plasmodium species with improved sensitivity for mixed infections. J Clin Microbiol. 2009;47:97580. DOIPubMedGoogle Scholar
  8. Menard  D, Ariey  F. PCR and sequencing for genotyping of candidate Plasmodium falciparum artemisinin resistance SNPs in the Kelch 13 gene. Paris: Institut Pasteur du Cambodge; 2013.
  9. Ariey  F, Witkowski  B, Amaratunga  C, Beghain  J, Langlois  AC, Khim  N, et al. A molecular marker of artemisinin-resistant Plasmodium falciparum malaria. Nature. 2014;505:505. DOIPubMedGoogle Scholar
  10. Koressaar  T, Remm  M. Enhancements and modifications of primer design program Primer3. Bioinformatics. 2007;23:128991. DOIPubMedGoogle Scholar
  11. Untergasser  A, Cutcutache  I, Koressaar  T, Ye  J, Faircloth  BC, Remm  M, et al. Primer3—new capabilities and interfaces. Nucleic Acids Res. 2012;40:e115. DOIPubMedGoogle Scholar
  12. Price  RN, Uhlemann  AC, Brockman  A, McGready  R, Ashley  E, Phaipun  L, et al. Mefloquine resistance in Plasmodium falciparum and increased pfmdr1 gene copy number. Lancet. 2004;364:43847. DOIPubMedGoogle Scholar
  13. Livak  KJ, Schmittgen  TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Δ Δ C(T)) Method. Methods. 2001;25:4028. DOIPubMedGoogle Scholar
  14. Lulitanond  V, Intapan  PM, Tantrawatpan  C, Sankuntaw  N, Sanpool  O, Janwan  P, et al. Molecular markers for detection and differentiation of Plasmodium falciparum and Plasmodium vivax in human blood samples by pyrosequencing. J Clin Microbiol. 2012;50:14557. DOIPubMedGoogle Scholar
  15. Dondorp  AM, Yeung  S, White  L, Nguon  C, Day  NPJ, Socheat  D, et al. Artemisinin resistance: current status and scenarios for containment. Nat Rev Microbiol. 2010;8:27280. DOIPubMedGoogle Scholar
  16. Djimdé  A, Doumbo  OK, Cortese  JF, Kayentao  K, Doumbo  S, Diourté  Y, et al. A molecular marker for chloroquine-resistant falciparum malaria. N Engl J Med. 2001;344:25763. DOIPubMedGoogle Scholar
  17. Olukosi  YA, Oyebola  MK, Ajibaye  O, Orok  BA, Aina  OO, Agomo  CO, et al. Persistence of markers of chloroquine resistance among P. falciparum isolates recovered from two Nigerian communities. MWJ. 2014;5:3.
  18. Mekonnen  SK, Aseffa  A, Berhe  N, Teklehaymanot  T, Clouse  RM, Gebru  T, et al. Return of chloroquine-sensitive Plasmodium falciparum parasites and emergence of chloroquine-resistant Plasmodium vivax in Ethiopia. Malar J. 2014;13:244. DOIPubMedGoogle Scholar
  19. Price  RN, Cassar  C, Brockman  A, Duraisingh  M, van Vugt  M, White  NJ, et al. The pfmdr1 gene is associated with a multidrug-resistant phenotype in Plasmodium falciparum from the western border of Thailand. Antimicrob Agents Chemother. 1999;43:29439. DOIPubMedGoogle Scholar
  20. Ngo  T, Duraisingh  M, Reed  M, Hipgrave  D, Biggs  B, Cowman  AF. Analysis of pfcrt, pfmdr1, dhfr, and dhps mutations and drug sensitivities in Plasmodium falciparum isolates from patients in Vietnam before and after treatment with artemisinin. Am J Trop Med Hyg. 2003;68:3506. DOIPubMedGoogle Scholar
  21. Zalis  MG, Pang  L, Silveira  MS, Milhous  WK, Wirth  DF. Characterization of Plasmodium falciparum isolated from the Amazon region of Brazil: evidence for quinine resistance. Am J Trop Med Hyg. 1998;58:6307. DOIPubMedGoogle Scholar
  22. Triglia  T, Cowman  AF. Primary structure and expression of the dihydropteroate synthetase gene of Plasmodium falciparum. Proc Natl Acad Sci U S A. 1994;91:714953. DOIPubMedGoogle Scholar
  23. Brooks  DR, Wang  P, Read  M, Watkins  WM, Sims  PF, Hyde  JE. Sequence variation of the hydroxymethyldihydropterin pyrophosphokinase: dihydropteroate synthase gene in lines of the human malaria parasite, Plasmodium falciparum, with differing resistance to sulfadoxine. Eur J Biochem. 1994;224:397405. DOIPubMedGoogle Scholar
  24. Mharakurwa  S, Kumwenda  T, Mkulama  MA, Musapa  M, Chishimba  S, Shiff  CJ, et al. Malaria antifolate resistance with contrasting Plasmodium falciparum dihydrofolate reductase (DHFR) polymorphisms in humans and Anopheles mosquitoes. Proc Natl Acad Sci U S A. 2011;108:18796801. DOIPubMedGoogle Scholar
  25. Peterson  DS, Walliker  D, Wellems  TE. Evidence that a point mutation in dihydrofolate reductase-thymidylate synthase confers resistance to pyrimethamine in falciparum malaria. Proc Natl Acad Sci U S A. 1988;85:91148. DOIPubMedGoogle Scholar
  26. Foote  SJ, Galatis  D, Cowman  AF. Amino acids in the dihydrofolate reductase-thymidylate synthase gene of Plasmodium falciparum involved in cycloguanil resistance differ from those involved in pyrimethamine resistance. Proc Natl Acad Sci U S A. 1990;87:30147. DOIPubMedGoogle Scholar
  27. Jambou  R, Legrand  E, Niang  M, Khim  N, Lim  P, Volney  B, et al. Resistance of Plasmodium falciparum field isolates to in-vitro artemether and point mutations of the SERCA-type PfATPase6. Lancet. 2005;366:19603. DOIPubMedGoogle Scholar
  28. Ashley  EA, Dhorda  M, Fairhurst  RM, Amaratunga  C, Lim  P, Suon  S, et al.; Tracking Resistance to Artemisinin Collaboration (TRAC). Spread of artemisinin resistance in Plasmodium falciparum malaria. N Engl J Med. 2014;371:41123. DOIPubMedGoogle Scholar
  29. Torrentino-Madamet  M, Fall  B, Benoit  N, Camara  C, Amalvict  R, Fall  M, et al. Limited polymorphisms in k13 gene in Plasmodium falciparum isolates from Dakar, Senegal in 2012-2013. Malar J. 2014;13:472. DOIPubMedGoogle Scholar
  30. Sutherland  CJ, Lansdell  P, Sanders  M, Muwanguzi  J, van Schalkwyk  DA, Kaur  H, et al. Pfk13-independent treatment failure in four imported cases of Plasmodium falciparum malaria treated with artemether-lumefantrine in the United Kingdom. Antimicrob Agents Chemother. 2017;61:e0238216. DOIPubMedGoogle Scholar
  31. Mohon  AN, Alam  MS, Bayih  AG, Folefoc  A, Shahinas  D, Haque  R, et al. Mutations in Plasmodium falciparum K13 propeller gene from Bangladesh (2009-2013). Malar J. 2014;13:431. DOIPubMedGoogle Scholar
  32. Tun  KM, Imwong  M, Lwin  KM, Win  AA, Hlaing  TM, Hlaing  T, et al. Spread of artemisinin-resistant Plasmodium falciparum in Myanmar: a cross-sectional survey of the K13 molecular marker. Lancet Infect Dis. 2015;15:41521. DOIPubMedGoogle Scholar
  33. Cooper  RA, Conrad  MD, Watson  QD, Huezo  SJ, Ninsiima  H, Tumwebaze  P, et al. Lack of artemisinin resistance in Plasmodium falciparum in Uganda based on parasitological and molecular assays. Antimicrob Agents Chemother. 2015;59:50614. DOIPubMedGoogle Scholar
  34. Ménard  D, Khim  N, Beghain  J, Adegnika  AA, Shafiul-Alam  M, Amodu  O, et al.; KARMA Consortium. KARMA Consortium. KARMA Consortium. A worldwide map of Plasmodium falciparum k13-propeller polymorphisms. N Engl J Med. 2016;374:245364. DOIPubMedGoogle Scholar
  35. Daily  JP. k13-propeller mutations and malaria resistance. N Engl J Med. 2016;374:24923. DOIPubMedGoogle Scholar
  36. Schwartz  E, Bujanover  S, Kain  KC. Genetic confirmation of atovaquone-proguanil-resistant Plasmodium falciparum malaria acquired by a nonimmune traveler to East Africa. Clin Infect Dis. 2003;37:4501. DOIPubMedGoogle Scholar
  37. Musset  L, Bouchaud  O, Matheron  S, Massias  L, Le Bras  J. Clinical atovaquone-proguanil resistance of Plasmodium falciparum associated with cytochrome b codon 268 mutations. Microbes Infect. 2006;8:2599604. DOIPubMedGoogle Scholar
  38. David  KP, Alifrangis  M, Salanti  A, Vestergaard  LS, Rønn  A, Bygbjerg  IB. Atovaquone/proguanil resistance in Africa: a case report. Scand J Infect Dis. 2003;35:8978. DOIPubMedGoogle Scholar
  39. Pimentel  S, Nogueira  F, Benchimol  C, Quinhentos  V, Bom  J, Varandas  L, et al. Detection of atovaquone-proguanil resistance conferring mutations in Plasmodium falciparum cytochrome b gene in Luanda, Angola. Malar J. 2006;5:30. DOIPubMedGoogle Scholar
  40. Savini  H, Bogreau  H, Bertaux  L, Bouchiba  H, Kraemer  P, Parzy  D, et al. First case of emergence of atovaquone-proguanil resistance in Plasmodium falciparum during treatment in a traveler in Comoros. Antimicrob Agents Chemother. 2008;52:22834. DOIPubMedGoogle Scholar
  41. Wichmann  O, Muehlen  M, Gruss  H, Mockenhaupt  FP, Suttorp  N, Jelinek  T. Malarone treatment failure not associated with previously described mutations in the cytochrome b gene. Malar J. 2004;3:14. DOIPubMedGoogle Scholar
  42. Liu  S, Mu  J, Jiang  H, Su  XZ. Effects of Plasmodium falciparum mixed infections on in vitro antimalarial drug tests and genotyping. Am J Trop Med Hyg. 2008;79:17884. DOIPubMedGoogle Scholar
  43. Sibley  CH, Ringwald  P. A database of antimalarial drug resistance. Malar J. 2006;5:48. DOIPubMedGoogle Scholar




Cite This Article

DOI: 10.3201/eid2804.210533

Original Publication Date: March 08, 2022

Table of Contents – Volume 28, Number 4—April 2022

EID Search Options
presentation_01 Advanced Article Search – Search articles by author and/or keyword.
presentation_01 Articles by Country Search – Search articles by the topic country.
presentation_01 Article Type Search – Search articles by article type and issue.



Please use the form below to submit correspondence to the authors or contact them at the following address:

Andrea K. Boggild, Tropical Disease Unit, Toronto General Hospital, 200 Elizabeth St, 13EN-218, Toronto, ON M5G 2C4, Canada

Send To

10000 character(s) remaining.


Page created: March 08, 2022
Page updated: March 19, 2022
Page reviewed: March 19, 2022
The conclusions, findings, and opinions expressed by authors contributing to this journal do not necessarily reflect the official position of the U.S. Department of Health and Human Services, the Public Health Service, the Centers for Disease Control and Prevention, or the authors' affiliated institutions. Use of trade names is for identification only and does not imply endorsement by any of the groups named above.