Genetic polymorphisms in P. falciparum can be used to indicate the parasite’s susceptibility to antimalarial drugs as well as its geographical origin. Both of these factors are key to monitoring development and spread of antimalarial drug resistance. In this study, we combine multiplex PCR, custom designed dual indexing and Miseq sequencing for high throughput SNP-profiling of 457 malaria infections from Guinea-Bissau, at the cost of 10 USD per sample. By amplifying and sequencing 15 genetic fragments, we cover 20 resistance-conferring SNPs occurring in pfcrt, pfmdr1, pfdhfr, pfdhps, as well as the entire length of pfK13, and the mitochondrial barcode for parasite origin. SNPs of interest were sequenced with an average depth of 2,043 reads, and bases were called for the various SNP-positions with a p-value below 0.05, for 89.8–100% of samples. The SNP data indicates that artemisinin resistance-conferring SNPs in pfK13 are absent from the studied area of Guinea-Bissau, while the pfmdr1 86 N allele is found at a high prevalence. The mitochondrial barcodes are unanimous and accommodate a West African origin of the parasites. With this method, very reliable high throughput surveillance of antimalarial drug resistance becomes more affordable than ever before.

Antimalarial drug resistance is a major constraint for malaria control and elimination efforts. Artemisinin-based combination therapy is now the mainstay for malaria treatment. However, delayed parasite clearance following treatment with artemisinin derivatives has now spread in the Greater Mekong Sub region and may emerge or spread to other malaria endemic regions. This spread is of great concern for malaria control programmes, as no alternatives to artemisinin-based combination therapies are expected to be available in the near future. There is a need to strengthen surveillance systems for early detection and response to the antimalarial drug resistance threat. Current surveillance is mainly done through therapeutic efficacy studies; however these studies are complex and both time- and resource-intensive. For multiple common antimalarials, parasite drug resistance has been correlated with specific genetic mutations, and the molecular markers associated with antimalarial drug resistance offer a simple and powerful tool to monitor the emergence and spread of resistant parasites. Different techniques to analyse molecular markers associated with antimalarial drug resistance are available, each with advantages and disadvantages. However, procedures are not adequately harmonized to facilitate comparisons between sites. Here we describe the target product profiles for tests to analyse molecular markers associated with antimalarial drug resistance, discuss how use of current techniques can be standardised, and identify the requirements for an ideal product that would allow malaria endemic countries to provide useful spatial and temporal information on the spread of resistance.

BackgroundSulfadoxine–pyrimethamine (SP) is recommended for prophylactic treatment of malaria in pregnancy while artemisinin combination therapy is the recommended first-line anti-malarial treatment. Selection of SP resistance is ongoing since SP is readily available in health facilities and in private drug shops in sub-Saharan Africa. This study reports on the prevalence and distribution of Pfdhps mutations A540E and A581G in Tanzania. When found together, these mutations confer high-level SP resistance (sometimes referred to as ‘super-resistance’), which is associated with loss in protective efficacy of SP-IPTp.MethodsDNA samples were extracted from malaria-positive blood samples on filter paper, used malaria rapid diagnostic test strips and whole blood collected from eight sites in seven administrative regions of Tanzania. PCR–RFLP and SSOP-ELISA techniques were used to genotype the A540E and A581G Pfdhps. Data were analysed using SPSS version 18 while Chi square and/or Fischer Exact tests were used to compare prevalence between regions.ResultsA high inter-regional variation of Pfdhps-540E was observed (χ2 = 76.8, p < 0.001). High inter-regional variation of 581G was observed (FE = 85.3, p < 0.001). Both Tanga and Kagera were found to have the highest levels of SP resistance. A high prevalence of Pfdhps-581G was observed in Tanga (56.6 %) in northeastern Tanzania and in Kagera (20.4 %) in northwestern Tanzania and the 540–581 EG haplotype was found at 54.5 and 19.4 %, respectively. Pfdhps-581G was not detected in Pwani and Lindi regions located south of Tanga region.ConclusionsSelection of SP super-resistant Pfdhps A581G is highest in northern Tanzania. Variation in distribution of SP resistance is observed across the country: northeastern Tanga region and northwestern Kagera region have highest prevalence of SP super-resistance markers, while in Pwani and Lindi in the southeast the prevalence of super-resistance was zero. More studies should be conducted to understand the factors underlying the remarkable heterogeneity in SP resistance in the country.