Hitchhiking and selective sweeps of Plasmodium falciparum sulfadoxine and pyrimethamine resistance alleles in a population from central Africa

A computational method for predicting the most likely evolutionary trajectories in the stepwise accumulation of resistance mutations

Pathogen evolution of drug resistance often occurs in a stepwise manner via the accumulation of multiple mutations that in combination have a non-additive impact on fitness, a phenomenon known as epistasis. The evolution of resistance via the accumulation of point mutations in the DHFR genes of Plasmodium falciparum (Pf) and Plasmodium vivax (Pv) has been studied extensively and multiple studies have shown epistatic interactions between these mutations determine the accessible evolutionary trajectories to highly resistant multiple mutations. Here, we simulated these evolutionary trajectories using a model of molecular evolution, parameterised using Rosetta Flex ddG predictions, where selection acts to reduce the target-drug binding affinity. We observe strong agreement with pathways determined using experimentally measured IC50 values of pyrimethamine binding, which suggests binding affinity is strongly predictive of resistance and epistasis in binding affinity strongly influences the order of fixation of resistance mutations. We also infer pathways directly from the frequency of mutations found in isolate data, and observe remarkable agreement with the most likely pathways predicted by our mechanistic model, as well as those determined experimentally. This suggests mutation frequency data can be used to intuitively infer evolutionary pathways, provided sufficient sampling of the population.

Assessment of Plasmodium falciparum drug resistance associated molecular markers in Mandla, Madhya Pradesh, India

BACKGROUND

Resistance against artemisinin-based combination therapy is one of the challenges to malaria control and elimination globally. Mutations in different genes (Pfdhfr, Pfdhps, Pfk-13 and Pfmdr1) confer resistance to artesunate and sulfadoxine-pyrimethamine (AS + SP) were analysed from Mandla district, Madhya Pradesh, to assess the effectiveness of the current treatment regimen against uncomplicated Plasmodium falciparum.

METHODS

Dried blood spots were collected during the active fever survey and mass screening and treatment activities as part of the Malaria Elimination Demonstration Project (MEDP) from 2019 to 2020. Isolated DNA samples were used to amplify the Pfdhfr, Pfdhps, Pfk13 and Pfmdr1 genes using nested PCR and sequenced for mutation analysis using the Sanger sequencing method.

RESULTS

A total of 393 samples were subjected to PCR amplification, sequencing and sequence analysis; 199, 215, 235, and 141 samples were successfully sequenced for Pfdhfr, Pfdhps, Pfk13, Pfmdr1, respectively. Analysis revealed that the 53.3% double mutation (C59R, S108N) in Pfdhfr, 89.3% single mutation (G437A) in Pfdhps, 13.5% single mutants (N86Y), and 51.1% synonymous mutations in Pfmdr1 in the study area. Five different non-synonymous and two synonymous point mutations found in Pfk13, which were not associated to artemisinin resistance.

CONCLUSION

The study has found that mutations linked to SP resistance are increasing in frequency, which may reduce the effectiveness of this drug as a future partner in artemisinin-based combinations. No evidence of mutations linked to artemisinin resistance in Pfk13 was found, suggesting that parasites are sensitive to artemisinin derivatives in the study area. These findings are a baseline for routine molecular surveillance to proactively identify the emergence and spread of artemisinin-resistant parasites.

Malaria Genomics, Vaccine Development, and Microbiome

Recent advances in malaria genetics and genomics have transformed many aspects of malaria research in areas of molecular evolution, epidemiology, transmission, host-parasite interaction, drug resistance, pathogenicity, and vaccine development. Here, in addition to introducing some background information on malaria parasite biology, parasite genetics/genomics, and genotyping methods, we discuss some applications of genetic and genomic approaches in vaccine development and in studying interactions with microbiota. Genetic and genomic data can be used to search for novel vaccine targets, design an effective vaccine strategy, identify protective antigens in a whole-organism vaccine, and evaluate the efficacy of a vaccine. Microbiota has been shown to influence disease outcomes and vaccine efficacy; studying the effects of microbiota in pathogenicity and immunity may provide information for disease control. Malaria genetics and genomics will continue to contribute greatly to many fields of malaria research.

A maximum-likelihood method to estimate haplotype frequencies and prevalence alongside multiplicity of infection from SNP data

The introduction of genomic methods facilitated standardized molecular disease surveillance. For instance, SNP barcodes in Plasmodium vivax and Plasmodium falciparum malaria allows the characterization of haplotypes, their frequencies and prevalence to reveal temporal and spatial transmission patterns. A confounding factor is the presence of multiple genetically distinct pathogen variants within the same infection, known as multiplicity of infection (MOI). Disregarding ambiguous information, as usually done in ad-hoc approaches, leads to less confident and biased estimates. We introduce a statistical framework to obtain maximum-likelihood estimates (MLE) of haplotype frequencies and prevalence alongside MOI from malaria SNP data, i.e., multiple biallelic marker loci. The number of model parameters increases geometrically with the number of genetic markers considered and no closed-form solution exists for the MLE. Therefore, the MLE needs to be derived numerically. We use the Expectation-Maximization (EM) algorithm to derive the maximum-likelihood estimates, an efficient and easy-to-implement algorithm that yields a numerically stable solution. We also derive expressions for haplotype prevalence based on either all or just the unambiguous genetic information and compare both approaches. The latter corresponds to a biased ad-hoc estimate of prevalence. We assess the performance of our estimator by systematic numerical simulations assuming realistic sample sizes and various scenarios of transmission intensity. For reasonable sample sizes, and number of loci, the method has little bias. As an example, we apply the method to a dataset from Cameroon on sulfadoxine-pyrimethamine resistance in P. falciparum malaria. The method is not confined to malaria and can be applied to any infectious disease with similar transmission behavior. An easy-to-use implementation of the method as an R-script is provided.

Genetic characterisation of the Theileria annulata cytochrome b locus and its impact on buparvaquone resistance in bovine

Control of tropical theileriosis, caused by the apicomplexan Theileria annulata, depends on the use of a single drug, buparvaquone, the efficacy of which is compromised by the emergence of resistance. The present study was undertaken to improve understanding of the role of mutations conferring buparvaquone resistance in T. annulata, and the effects of selection pressures on their emergence and spread. First, we investigated genetic characteristics of the cytochrome b locus associated with buparvaquone resistance in 10 susceptible and 7 resistant T. annulata isolates. The 129G (GGC) mutation was found in the Q₀₁ binding pocket and 253S (TCT) and 262S (TCA) mutations were identified within the Q₀₂ binding pocket. Next, we examined field isolates and identified cytochrome b mutations 129G (GGC), 253S (TCT) and 262S (TCA) in 21/75 buffalo-derived and 19/119 cattle-derived T. annulata isolates, providing evidence of positive selection pressure. Both hard and soft selective sweeps were identified, with striking differences between isolates. For example, 19 buffalo-derived and 7 cattle-derived isolates contained 129G (GGC) and 253S (TCT) resistance haplotypes at a high frequency, implying the emergence of resistance by a single mutation. Two buffalo-derived and 12 cattle-derived isolates contained equally high frequencies of 129G (GGC), 253S (TCT), 129G (GGC)/253S (TCT) and 262S (TCA) resistance haplotypes, implying the emergence of resistance by pre-existing or recurrent mutations. Phylogenetic analysis further revealed that 9 and 21 unique haplotypes in buffalo and cattle-derived isolates were present in a single lineage, suggesting a single origin. We propose that animal migration between farms is an important factor in the spread of buparvaquone resistance in endemic regions of Pakistan. The overall outcomes will be useful in understanding how drug resistance emerges and spreads, and this information will help design strategies to optimise the use and lifespan of the single most drug use to control tropical theileriosis.

Anti-folate quintuple mutations in Plasmodium falciparum asymptomatic infections in Yaoundé, Cameroon

A major challenge in the fight to effectively control malaria is the emergence of resistant parasite to drugs used in therapy as well as for chemoprevention. In this study, single nucleotide polymorphisms (SNPs) associated with Plasmodium falciparum resistance to sulfadoxine-pyrimethamine (SP), one of the partner drugs in artemisinin-based therapies (ACTs) were studied in asymptomatic P. falciparum isolates from Cameroon. Dried Blood spots were collected from children with asymptomatic malaria enrolled during a household survey. The P. falciparum dihydrofolate reductase (Pfdhfr), dihydropteroate synthase (Pfdhps) and Kelch 13 genes were amplified and point mutations in these gene sequences were analyzed by sequencing. Among a total of 234 samples collected, 51 showed parasitaemia after microscopic examination of which 47 were P. falciparum mono-infections. Molecular analysis revealed 97.3% of mutant alleles at codons 51I, 59R and 108 N in Pfdhfr gene. In Pfdhps gene the most common mutation was 437G (83.3%); followed by 436A (47.6%) and 436F (28.6%). The association of mutations in the two genes (dhfr + dhps) showed 11 different haplotypes including three sextuple mutants (IRNI + AGKGA, IRNI + AAKGS, IRNI + AGKAS) and one septuple mutant (IRNI + AGKGS). For K13 gene no SNPs were seen in the studied asymptomatic malaria samples. The findings revealed presence of SP-resistant alleles in asymptomatic infected individuals with presence of sextuples and septuple SNPs. This emphasizes that regular profiling of antimalarial drugs resistance markers in such population is essential for malaria control and elimination programmes.

Population-level genome-wide STR discovery and validation for population structure and genetic diversity assessment of Plasmodium species

Short tandem repeats (STRs) are highly informative genetic markers that have been used extensively in population genetics analysis. They are an important source of genetic diversity and can also have functional impact. Despite the availability of bioinformatic methods that permit large-scale genome-wide genotyping of STRs from whole genome sequencing data, they have not previously been applied to sequencing data from large collections of malaria parasite field samples. Here, we have genotyped STRs using HipSTR in more than 3,000 Plasmodium falciparum and 174 Plasmodium vivax published whole-genome sequence data from samples collected across the globe. High levels of noise and variability in the resultant callset necessitated the development of a novel method for quality control of STR genotype calls. A set of high-quality STR loci (6,768 from P. falciparum and 3,496 from P. vivax) were used to study Plasmodium genetic diversity, population structures and genomic signatures of selection and these were compared to genome-wide single nucleotide polymorphism (SNP) genotyping data. In addition, the genome-wide information about genetic variation and other characteristics of STRs in P. falciparum and P. vivax have been available in an interactive web-based R Shiny application PlasmoSTR (https://github.com/bahlolab/PlasmoSTR).

Polymorphism of Antifolate Drug Resistance in Plasmodium vivax From Local Residents and Migrant Workers Returned From the China-Myanmar Border

Drug-resistant Plasmodium vivax malaria impedes efforts to control, eliminate, and ultimately eradicate malaria in Southeast Asia. P. vivax resistance to antifolate drugs derives from point mutations in specific parasite genes, including the dihydropteroate synthase (pvdhps), dihydrofolate reductase (pvdhfr), and GTP cyclohydrolase I (pvgch1) genes. This study aims to investigate the prevalence and spread of drug resistance markers in P. vivax populating the China-Myanmar border. Blood samples were collected from symptomatic patients with acute P. vivax infection. Samples with single-clone P. vivax infections were sequenced for pvdhps and pvdhfr genes and genotyped for 6 flanking microsatellite markers. Copy number variation in the pvgch1 gene was also examined. Polymorphisms were observed in six different codons of the pvdhps gene (382, 383, 512, 549, 553, and 571) and six different codons of the pvdhfr gene (13, 57, 58, 61, 99, 117) in two study sites. The quadruple mutant haplotypes 57I/L/58R/61M/117T of pvdhfr gene were the most common (comprising 76% of cases in Myitsone and 43.7% of case in Laiza). The double mutant haplotype 383G/553G of pvdhps gene was also prevalent at each site (40.8% and 31%). Microsatellites flanking the pvdhfr gene differentiated clinical samples from wild type and quadruple mutant genotypes (F_ST= 0.259-0.3036), as would be expected for a locus undergoing positive selection. The lack of copy number variation of pvgch1 suggests that SP-resistant P. vivax may harbor alternative mechanisms to secure sufficient folate.

Targeted deep amplicon sequencing of antimalarial resistance markers in Plasmodium falciparum isolates from Cameroon

BACKGROUND

Recent studies showed the first emergence of the R561H artemisinin-associated resistance marker in Africa, which highlights the importance of continued molecular surveillance to assess the selection and spread of this and other drug resistance markers in the region.

METHOD

In this study, we used targeted amplicon deep sequencing of 116 isolates collected in two areas of Cameroon to genotype the major drug resistance genes, k13, crt, mdr1, dhfr, and dhps, and the cytochrome b gene (cytb) in Plasmodium falciparum.

RESULTS

No confirmed or associated artemisinin resistance markers were observed in Pfk13. In comparison, both major and minor alleles associated with drug resistance were found in Pfcrt, Pfmdr1, Pfdhfr, and Pfdhps. Notably, a high frequency of other nonsynonymous mutations was observed across all the genes, except for Pfcytb, suggesting continued selection pressure.

CONCLUSIONS

The results from this study supported the continued use of artemisinin-based combination therapy and administration of sulfadoxine-pyrimethamine for intermittent preventive therapy in pregnant women, and for seasonal chemoprevention in these study sites in Cameroon.

Association of Plasmodium falciparum kelch13 R561H genotypes with delayed parasite clearance in Rwanda: an open-label, single-arm, multicentre, therapeutic efficacy study

BACKGROUND

Partial artemisinin resistance is suspected if delayed parasite clearance (ie, persistence of parasitaemia on day 3 after treatment initiation) is observed. Validated markers of artemisinin partial resistance in southeast Asia, Plasmodium falciparum kelch13 (Pfkelch13) R561H and P574L, have been reported in Rwanda but no association with parasite clearance has been observed. We aimed to establish the efficacy of artemether-lumefantrine and genetic characterisation of Pfkelch13 alleles and their association with treatment outcomes.

METHODS

This open-label, single-arm, multicentre, therapeutic efficacy study was done in 2018 in three Rwandan sites: Masaka, Rukara, and Bugarama. Children aged 6-59 months with P falciparum monoinfection and fever were eligible and treated with a 3-day course of artemether-lumefantrine. Treatment response was monitored for 28 days using weekly microscopy screenings of blood samples for P falciparum. Mutations in Pfkelch13 and P falciparum multidrug resistance-1 (Pfmdr1) genes were characterised in parasites collected from enrolled participants. Analysis of flanking microsatellites surrounding Pfkelch13 was done to define the origins of the R561H mutations. The primary endpoint was PCR-corrected parasitological cure on day 28, as per WHO protocol.

FINDINGS

228 participants were enrolled and 224 (98·2%) reached the study endpoint. PCR-corrected efficacies were 97·0% (95% CI 88-100) in Masaka, 93·8% (85-98) in Rukara, and 97·2% (91-100) in Bugarama. Pfkelch13 R561H mutations were present in 28 (13%) of 218 pre-treatment samples and P574L mutations were present in two (1%) pre-treatment samples. 217 (90%) of the 240 Pfmdr1 haplotypes observed in the pretreatment samples, had either the NFD (N86Y, Y184F, D1246Y) or NYD haplotype. Eight (16%) of 51 participants in Masaka and 12 (15%) of 82 participants in Rukara were microscopically positive 3 days after treatment initiation, which was associated with pre-treatment presence of Pfkelch13 R561H in Masaka (p=0·0005). Genetic analysis of Pfkelch13 R561H mutations suggest their common ancestry and local origin in Rwanda.

INTERPRETATION

We confirm evidence of emerging artemisinin partial resistance in Rwanda. Although artemether-lumefantrine remains efficacious, vigilance for decreasing efficacy, further characterisation of artemisinin partial resistance, and evaluation of additional antimalarials in Rwanda should be considered.

FUNDING

The US President's Malaria Initiative.

TRANSLATION

For the French translation of the abstract see Supplementary Materials section.

Phylogenetic analysis suggests single and multiple origins of dihydrofolate reductase mutations in Plasmodium vivax

Pyrimethamine was first introduced for the treatment of malaria in Asia and Africa during the early 1980s, replacing chloroquine, and has become the first line of drugs in many countries. In recent years, development of pyrimethamine resistance in Plasmodium vivax has become a barrier to effective malaria control strategies. Here, we describe the use of meta-barcoded deep amplicon sequencing technology to assess the evolutionary origin of pyrimethamine resistance by analysing the flanking region of dihydrofolate reductase (dhfr) locus. The genetic modelling suggests that 58R and 173L single mutants and 58R/117N double mutants are present on a single lineage; suggesting a single origin of these mutations. The triple mutants (57L/58R/117N, 58R/61M/117N and 58R/117N/173L) share the lineage of 58R/117N, suggesting a common origin. In contrast, the 117N mutant is present on two separate lineages suggesting that there are multiple origins of this mutation. We characterised the allele frequency of the P. vivax dhfr locus. Our results support the view that the single mutation of 117N and double mutations of 58R/117N arise commonly, whereas the single mutation of 173L and triple mutations of 57L/58R/117N, 58R/61M/117N and 58R/117N/173L are less common. Our work will help to inform mitigation strategies for pyrimethamine resistance in P. vivax.

Drug resistance markers within an evolving efficacy of anti-malarial drugs in Cameroon: a systematic review and meta-analysis (1998-2020)

BACKGROUND

Malaria remains highly endemic in Cameroon. The rapid emergence and spread of drug resistance was responsible for the change from monotherapies to artemisinin-based combinations. This systematic review and meta-analysis aimed to determine the prevalence and distribution of Plasmodium falciparum drug resistance markers within an evolving efficacy of anti-malarial drugs in Cameroon from January 1998 to August 2020.

METHODS

The PRISMA-P and PRISMA statements were adopted in the inclusion of studies on single nucleotide polymorphisms (SNPs) of P. falciparum anti-malarial drug resistance genes (Pfcrt, Pfmdr1, Pfdhfr, Pfdhps, Pfatp6, Pfcytb and Pfk13). The heterogeneity of the included studies was evaluated using the Cochran's Q and I2 statistics. The random effects model was used as standard in the determination of heterogeneity between studies.

RESULTS

Out of the 902 records screened, 48 studies were included in this aggregated meta-analysis of molecular data. A total of 18,706 SNPs of the anti-malarial drug resistance genes were genotyped from 47,382 samples which yielded a pooled prevalence of 35.4% (95% CI 29.1-42.3%). Between 1998 and 2020, there was significant decline (P < 0.0001 for all) in key mutants including Pfcrt 76 T (79.9%-43.0%), Pfmdr1 86Y (82.7%-30.5%), Pfdhfr 51I (72.2%-66.9%), Pfdhfr 59R (76.5%-67.8%), Pfdhfr 108 N (80.8%-67.6%). The only exception was Pfdhps 437G which increased over time (30.4%-46.9%, P < 0.0001) and Pfdhps 540E that remained largely unchanged (0.0%-0.4%, P = 0.201). Exploring mutant haplotypes, the study observed a significant increase in the prevalence of Pfcrt CVIET mixed quintuple haplotype from 57.1% in 1998 to 57.9% in 2020 (P < 0.0001). In addition, within the same study period, there was no significant change in the triple Pfdhfr IRN mutant haplotype (66.2% to 67.3%, P = 0.427). The Pfk13 amino acid polymorphisms associated with artemisinin resistance were not detected.

CONCLUSIONS

This review reported an overall decline in the prevalence of P. falciparum gene mutations conferring resistance to 4-aminoquinolines and amino alcohols for a period over two decades. Resistance to artemisinins measured by the presence of SNPs in the Pfk13 gene does not seem to be a problem in Cameroon. Systematic review registration PROSPERO CRD42020162620.

Host-Malaria Parasite Interactions and Impacts on Mutual Evolution

Malaria is the most deadly parasitic disease, affecting hundreds of millions of people worldwide. Malaria parasites have been associated with their hosts for millions of years. During the long history of host-parasite co-evolution, both parasites and hosts have applied pressure on each other through complex host-parasite molecular interactions. Whereas the hosts activate various immune mechanisms to remove parasites during an infection, the parasites attempt to evade host immunity by diversifying their genome and switching expression of targets of the host immune system. Human intervention to control the disease such as antimalarial drugs and vaccination can greatly alter parasite population dynamics and evolution, particularly the massive applications of antimalarial drugs in recent human history. Vaccination is likely the best method to prevent the disease; however, a partially protective vaccine may have unwanted consequences that require further investigation. Studies of host-parasite interactions and co-evolution will provide important information for designing safe and effective vaccines and for preventing drug resistance. In this essay, we will discuss some interesting molecules involved in host-parasite interactions, including important parasite antigens. We also discuss subjects relevant to drug and vaccine development and some approaches for studying host-parasite interactions.

Plasmodium Genomics and Genetics: New Insights into Malaria Pathogenesis, Drug Resistance, Epidemiology, and Evolution

Protozoan Plasmodium parasites are the causative agents of malaria, a deadly disease that continues to afflict hundreds of millions of people every year. Infections with malaria parasites can be asymptomatic, with mild or severe symptoms, or fatal, depending on many factors such as parasite virulence and host immune status. Malaria can be treated with various drugs, with artemisinin-based combination therapies (ACTs) being the first-line choice. Recent advances in genetics and genomics of malaria parasites have contributed greatly to our understanding of parasite population dynamics, transmission, drug responses, and pathogenesis. However, knowledge gaps in parasite biology and host-parasite interactions still remain. Parasites resistant to multiple antimalarial drugs have emerged, while advanced clinical trials have shown partial efficacy for one available vaccine. Here we discuss genetic and genomic studies of Plasmodium biology, host-parasite interactions, population structures, mosquito infectivity, antigenic variation, and targets for treatment and immunization. Knowledge from these studies will advance our understanding of malaria pathogenesis, epidemiology, and evolution and will support work to discover and develop new medicines and vaccines.

Selective sweep and phylogenetic models for the emergence and spread of pyrimethamine resistance mutations in Plasmodium vivax

Pyrimethamine resistance is a major concern for the control of human haemoprotozoa, especially Plasmodium species. Currently, there is little understanding of how pyrimethamine resistance developed in Plasmodium vivax in the natural field conditions. Here, we present for the first time evidence of positive selection pressure on a dihydrofolate reductase locus and its consequences on the emergence and the spread of pyrimethamine resistance in P. vivax in the Punjab province of Pakistan. First, we examined the dihydrofolate reductase locus in 38 P. vivax isolates to look for evidence of positive selection pressure in human patients. The S58R (AGA)/S117N (AAC) double mutation was most common, being detected in 10/38 isolates. Single mutation S117N (AAC), I173L (CTT) and S58R (AGA) SNPs were detected in 8/38, 2/38 and 1/38 isolates, respectively. The F57L/I (TTA/ATA) and T61M (ATG) SNPs were not detected in any isolates examined. Although both soft and hard selective sweeps have occurred with striking differences between isolates, there was a predominance of hard sweeps. A single resistance haplotype was present at high frequency in 9/14 isolates, providing a strong evidence for single emergence of resistance by the single mutation, characteristics of hard selective sweeps. In contrast, 5/14 isolates carried multiple resistance haplotypes at high frequencies, providing an evidence of the emergence of resistance by recurrent mutations, characteristics of soft selective sweeps. Our phylogenetic relationship analysis suggests that S58R (AGA)/S117N (AAC) and S117N (AAC) mutations arose multiple times from a single origin and spread to multiple different cities in the Punjab province through gene flow. Interestingly, the I173L (CTT) mutation was present on a single haplotype, suggesting that it arises rarely and has not spread between cities. Our work shows the need for responsible use of existing and new antimicrobial drugs and their combinations, control the movement of infected patients and mosquito vector control strategies.

Plasmodium genomics: an approach for learning about and ending human malaria

Malaria causes high levels of morbidity and mortality in human beings worldwide. According to the World Health Organization (WHO), about half a million people die of this disease each year. Malaria is caused by six species of parasites belonging to the Plasmodium genus: P. falciparum, P. knowlesi, P. vivax, P. malariae, P. ovale curtisi, and P. ovale wallikeri. Currently, malaria is being kept under control with varying levels of elimination success in different countries. The development of new molecular tools as well as the use of next-generation sequencing (NGS) technologies and novel bioinformatic approaches has improved our knowledge of malarial epidemiology, diagnosis, treatment, vaccine development, and surveillance strategies. In this work, the genetics and genomics of human malarias have been analyzed. Since the first P. falciparum genome was sequenced in 2002, various population-level genetic and genomic surveys, together with transcriptomic and proteomic studies, have shown the importance of molecular approaches in supporting malaria elimination.

Plasmodium falciparum Genetic Diversity in Continental Equatorial Guinea before and after Introduction of Artemisinin-Based Combination Therapy

Efforts to control malaria may affect malaria parasite genetic variability and drug resistance, the latter of which is associated with genetic events that promote mechanisms to escape drug action. The worldwide spread of drug resistance has been a major obstacle to controlling Plasmodium falciparum malaria, and thus the study of the origin and spread of associated mutations may provide some insights into the prevention of its emergence. This study reports an analysis of P. falciparum genetic diversity, focusing on antimalarial resistance-associated molecular markers in two socioeconomically different villages in mainland Equatorial Guinea. The present study took place 8 years after a previous one, allowing the analysis of results before and after the introduction of an artemisinin-based combination therapy (ACT), i.e., artesunate plus amodiaquine. Genetic diversity was assessed by analysis of the Pfmsp2 gene and neutral microsatellite loci. Pfdhps and Pfdhfr alleles associated with sulfadoxine-pyrimethamine (SP) resistance and flanking microsatellite loci were investigated, and the prevalences of drug resistance-associated point mutations of the Pfcrt, Pfmdr1, Pfdhfr, and Pfdhps genes were estimated. Further, to monitor the use of ACT, we provide the baseline prevalences of K13 propeller mutations and Pfmdr1 copy numbers. After 8 years, noticeable differences occurred in the distribution of genotypes conferring resistance to chloroquine and SP, and the spread of mutated genotypes differed according to the setting. Regarding artemisinin resistance, although mutations reported as being linked to artemisinin resistance were not present at the time, several single nucleotide polymorphisms (SNPs) were observed in the K13 gene, suggesting that closer monitoring should be maintained to prevent the possible spread of artemisinin resistance in Africa.

Whole Genome Sequencing of Field Isolates Reveals Extensive Genetic Diversity in Plasmodium vivax from Colombia

Plasmodium vivax is the most prevalent malarial species in South America and exerts a substantial burden on the populations it affects. The control and eventual elimination of P. vivax are global health priorities. Genomic research contributes to this objective by improving our understanding of the biology of P. vivax and through the development of new genetic markers that can be used to monitor efforts to reduce malaria transmission. Here we analyze whole-genome data from eight field samples from a region in Cordóba, Colombia where malaria is endemic. We find considerable genetic diversity within this population, a result that contrasts with earlier studies suggesting that P. vivax had limited diversity in the Americas. We also identify a selective sweep around a substitution known to confer resistance to sulphadoxine-pyrimethamine (SP). This is the first observation of a selective sweep for SP resistance in this species. These results indicate that P. vivax has been exposed to SP pressure even when the drug is not in use as a first line treatment for patients afflicted by this parasite. We identify multiple non-synonymous substitutions in three other genes known to be involved with drug resistance in Plasmodium species. Finally, we found extensive microsatellite polymorphisms. Using this information we developed 18 polymorphic and easy to score microsatellite loci that can be used in epidemiological investigations in South America.

Although P. vivax is not as deadly as the more widely studied P. falciparum, it remains a pressing global health problem. Here we report the results of a whole-genome study of P. vivax from Cordóba, Colombia, in South America. This parasite is the most prevalent in this region. We show that the parasite population is genetically diverse, which is contrary to expectations from earlier studies from the Americas. We also find molecular evidence that resistance to an anti-malarial drug has arisen recently in this region. This selective sweep indicates that the parasite has been exposed to a drug that is not used as first-line treatment for this malaria parasite. In addition to extensive single nucleotide and microsatellite polymorphism, we report 18 new genetic loci that might be helpful for fine-scale studies of this species in the Americas.

Prevalence of Plasmodium falciparum resistance markers to sulfadoxine-pyrimethamine among pregnant women receiving intermittent preventive treatment for malaria in Uganda

The aim of this study was to assess the prevalence of mutations in Plasmodium falciparum dihydrofolate reductase (Pfdhfr) and dihydropteroate synthase (Pfdhps) genes among pregnant women using sulfadoxine-pyrimethamine (SP) as an intermittent preventive treatment (IPTp). A molecular epidemiological study of P. falciparum parasite resistance markers to SP was conducted from August 2010 to February 2012 in Mukono district in central Uganda. DNA was extracted from 413 P. falciparum-positive samples. Real-time PCR, followed by melting curve analysis, was used to characterize point mutations in the Pfdhfr and Pfdhps genes that are associated with SP resistance. The prevalence of the single-nucleotide mutations in Pfdhfr at codons 51I, 59R, and 108N and in Pfdhps at codons 437G and 540E was high (>98%), reaching 100% fixation after one dose of SP, while the prevalence of 581G was 3.3% at baseline, reaching 12.5% after one dose of SP. At baseline, the prevalence of Pfdhfr and Pfdhps quintuple mutations was 89%, whereas the sextuple mutations (including 581G) were not prevalent (3.9%), reaching 16.7% after one dose of SP. However, the numbers of infections at follow-up visits were small, and hence there was insufficient statistical power to test whether there was a true rise in the prevalence of this allele. The overall high frequency of Pfdhfr and Pfdhps quintuple mutations throughout pregnancy excluded further analyses of possible associations between certain haplotypes and the risk of lower birth weight and anemia. However, women infected with P. falciparum had 1.3-g/dl-lower hemoglobin levels (P = 0.001) and delivered babies with a 400-g-lower birth weight (P = 0.001) compared to nonparasitemic women. Despite this, 44 women who were P. falciparum positive at baseline became negative after one or two doses of SP (i.e., 50.5%), implying that SP-IPTp still has some efficacy. P. falciparum resistance markers to SP are high in this population, whereas P. falciparum infection was associated with poor birth outcomes.

Prevalence of Plasmodium falciparum parasites resistant to sulfadoxine/pyrimethamine in pregnant women in Yaoundé, Cameroon: emergence of highly resistant pfdhfr/pfdhps alleles

OBJECTIVES

To determine, 6 years after the adoption of intermittent preventive treatment of pregnant women with sulfadoxine/pyrimethamine (IPTp-SP) in Cameroon, (i) the polymorphism and prevalence of Plasmodium falciparum dihydrofolate reductase (pfdhfr) and dihydropteroate synthase (pfdhps) gene mutations associated with sulfadoxine/pyrimethamine resistance and (ii) the consequences of sulfadoxine/pyrimethamine use in the selection of pfdhfr/pfdhps alleles.

METHODS

pfdhfr and pfdhps genes from P. falciparum isolates collected in Yaoundé (Cameroon) from pregnant women with symptomatic malaria before taking IPTp-SP [SP- group (control) (n = 51)] or afterwards [SP+ group (n = 49)] were sequenced.

RESULTS

The pfdhfr N51I, C59R, S108N triple mutant had a prevalence close to 100% (96/100) and no mutations at codons 50 and 164 were detected in either of the groups. The most frequent pfdhps mutation was A437G with a prevalence of 76.5% (39/51) in the SP- group, which was significantly higher in pregnant women who took sulfadoxine/pyrimethamine [95.9% (47/49)] (P = 0.012). Our study confirmed the presence of the pfdhps K540E mutation in Cameroon, but it remained rare. The prevalence of pfdhps A581G and A613S mutations had increased [5.9% (3/51) and 11.8% (6/51) in the control group, respectively] since the last studies in 2005. Surprisingly, the new pfdhps I431V mutation was detected, at a prevalence of 9.8% (5/51), and was found to be associated with other pfdhfr/pfdhps alleles to form an octuple N51I, C59R, S108N/I431V, S436A, A437G, A581G, A613S mutant.

CONCLUSIONS

Significant changes were found in pfdhps polymorphism. In particular, we observed several parasites carrying eight mutations in pfdhfr/pfdhps genes, which are very susceptible to having a high level of resistance to sulfadoxine/pyrimethamine.

Increasing prevalence of a novel triple-mutant dihydropteroate synthase genotype in Plasmodium falciparum in western Kenya

The molecular basis of sulfadoxine-pyrimethamine (SP) resistance lies in a combination of single-nucleotide polymorphisms (SNPs) in two genes coding for Plasmodium falciparum dihydrofolate reductase (Pfdhfr) and P. falciparum dihydropteroate synthase (Pfdhps), targeted by pyrimethamine and sulfadoxine, respectively. The continued use of SP for intermittent preventive treatment in pregnant women in many African countries, despite SP's discontinuation as a first-line antimalarial treatment option due to high levels of drug resistance, may further increase the prevalence of SP-resistant parasites and/or lead to the selection of new mutations. An antimalarial drug resistance surveillance study was conducted in western Kenya between 2010 and 2013. A total of 203 clinical samples from children with uncomplicated malaria were genotyped for SNPs associated with SP resistance. The prevalence of the triple-mutant Pfdhfr C50 I51R59N108: I164 genotype and the double-mutant Pfdhps S436 G437E540: A581A613 genotype was high. Two triple-mutant Pfdhps genotypes, S436 G437E540G581: A613 and H436G437E540: A581A613, were found, with the latter thus far being uniquely found in western Kenya. The prevalence of the S436 G437E540G581: A613 genotype was low. However, a steady increase in the prevalence of the Pfdhps triple-mutant H436G437E540: A581A613 genotype has been observed since its appearance in early 2000. Isolates with these genotypes shared substantial microsatellite haplotypes with the most common double-mutant allele, suggesting that this triple-mutant allele may have evolved locally. Overall, these findings show that the prevalence of the H436G437E540: A581A613 triple mutant may be increasing in this population and could compromise the efficacy of SP for intermittent preventive treatment in pregnant women if it increases the resistance threshold further.

Evolutionary biology and the avoidance of antimicrobial resistance

Evolutionary biologists have largely left the search for solutions to the drug resistance crisis to biomedical scientists, physicians, veterinarians and public health specialists. We believe this is because the vast majority of professional evolutionary biologists consider the evolutionary science of drug resistance to be conceptually uninteresting. Using malaria as case study, we argue that it is not. We review examples of evolutionary thinking that challenge various fallacies dominating antimalarial therapy, and discuss open problems that need evolutionary insight. These problems are unlikely to be resolved by biomedical scientists ungrounded in evolutionary biology. Involvement by evolutionary biologists in the science of drug resistance requires no intellectual compromises: the problems are as conceptually challenging as they are important.

Human-induced evolution caught in action: SNP-array reveals rapid amphi-atlantic spread of pesticide resistance in the salmon ecotoparasite Lepeophtheirus salmonis

Background

The salmon louse, Lepeophtheirus salmonis, is an ectoparasite of salmonids that causes huge economic losses in salmon farming, and has also been causatively linked with declines of wild salmonid populations. Lice control on farms is reliant upon a few groups of pesticides that have all shown time-limited efficiency due to resistance development. However, to date, this example of human-induced evolution is poorly documented at the population level due to the lack of molecular tools. As such, important evolutionary and management questions, linked to the development and dispersal of pesticide resistance in this parasite, remain unanswered. Here, we introduce the first Single Nucleotide Polymorphism (SNP) array for the salmon louse, which includes 6000 markers, and present a population genomic scan using this array on 576 lice from twelve farms distributed across the North Atlantic.

Results

Our results support the hypothesis of a single panmictic population of lice in the Atlantic, and importantly, revealed very strong selective sweeps on linkage groups 1 and 5. These sweeps included candidate genes potentially connected to pesticide resistance. After genotyping a further 576 lice from 12 full sibling families, a genome-wide association analysis established a highly significant association between the major sweep on linkage group 5 and resistance to emamectin benzoate, the most widely used pesticide in salmonid aquaculture for more than a decade.

Conclusions

The analysis of conserved haplotypes across samples from the Atlantic strongly suggests that emamectin benzoate resistance developed at a single source, and rapidly spread across the Atlantic within the period 1999 when the chemical was first introduced, to 2010 when samples for the present study were obtained. These results provide unique insights into the development and spread of pesticide resistance in the marine environment, and identify a small genomic region strongly linked to emamectin benzoate resistance. Finally, these results have highly significant implications for the way pesticide resistance is considered and managed within the aquaculture industry.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-937) contains supplementary material, which is available to authorized users.

A likelihood approach to estimate the number of co-infections

The number of co-infections of a pathogen (multiplicity of infection or MOI) is a relevant parameter in epidemiology as it relates to transmission intensity. Notably, such quantities can be built into a metric in the context of disease control and prevention. Having applications to malaria in mind, we develop here a maximum-likelihood (ML) framework to estimate the quantities of interest at low computational and no additional costs to study designs or data collection. We show how the ML estimate for the quantities of interest and corresponding confidence-regions are obtained from multiple genetic loci. Assuming specifically that infections are rare and independent events, the number of infections per host follows a conditional Poisson distribution. Under this assumption, we show that a unique ML estimate for the parameter () describing MOI exists which is found by a simple recursion. Moreover, we provide explicit formulas for asymptotic confidence intervals, and show that profile-likelihood-based confidence intervals exist, which are found by a simple two-dimensional recursion. Based on the confidence intervals we provide alternative statistical tests for the MOI parameter. Finally, we illustrate the methods on three malaria data sets. The statistical framework however is not limited to malaria.

Changes in malaria parasite drug resistance in an endemic population over a 25-year period with resulting genomic evidence of selection

Background. Analysis of genome-wide polymorphism in many organisms has potential to identify genes under recent selection. However, data on historical allele frequency changes are rarely available for direct confirmation.

Methods. We genotyped single nucleotide polymorphisms (SNPs) in 4 Plasmodium falciparum drug resistance genes in 668 archived parasite-positive blood samples of a Gambian population between 1984 and 2008. This covered a period before antimalarial resistance was detected locally, through subsequent failure of multiple drugs until introduction of artemisinin combination therapy. We separately performed genome-wide sequence analysis of 52 clinical isolates from 2008 to prospect for loci under recent directional selection.

Results. Resistance alleles increased from very low frequencies, peaking in 2000 for chloroquine resistance-associated crt and mdr1 genes and at the end of the survey period for dhfr and dhps genes respectively associated with pyrimethamine and sulfadoxine resistance. Temporal changes fit a model incorporating likely selection coefficients over the period. Three of the drug resistance loci were in the top 4 regions under strong selection implicated by the genome-wide analysis.

Conclusions. Genome-wide polymorphism analysis of an endemic population sample robustly identifies loci with detailed documentation of recent selection, demonstrating power to prospectively detect emerging drug resistance genes.

Surveillance of travellers: an additional tool for tracking antimalarial drug resistance in endemic countries

INTRODUCTION

There are growing concerns about the emergence of resistance to artemisinin-based combination therapies (ACTs). Since the widespread adoption of ACTs, there has been a decrease in the systematic surveillance of antimalarial drug resistance in many malaria-endemic countries. The aim of this work was to test whether data on travellers returning from Africa with malaria could serve as an additional surveillance system of local information sources for the emergence of drug resistance in endemic-countries.

METHODOLOGY

Data were collected from travellers with symptomatic Plasmodium falciparum malaria returning from Senegal (n = 1,993), Mali (n = 2,372), Cote d'Ivoire (n = 4,778) or Cameroon (n = 3,272) and recorded in the French Malaria Reference Centre during the period 1996-2011. Temporal trends of the proportion of parasite isolates that carried the mutant genotype, pfcrt 76T, a marker of resistance to chloroquine (CQ) and pfdhfr 108N, a marker of resistance to pyrimethamine, were compared for travellers and within-country surveys that were identified through a literature review in PubMed. The in vitro response to CQ was also compared between these two groups for parasites from Senegal.

RESULTS

The trends in the proportion of parasites that carried pfcrt 76T, and pfdhfr 108N, were compared for parasites from travellers and patients within-country using the slopes of the curves over time; no significant differences in the trends were found for any of the 4 countries. These results were supported by in vitro analysis of parasites from the field in Senegal and travellers returning to France, where the trends were also not significantly different.

CONCLUSION

The results have not shown different trends in resistance between parasites derived from travellers or from parasites within-country. This work highlights the value of an international database of drug responses in travellers as an additional tool to assess the emergence of drug resistance in endemic areas where information is limited.

Distinctive origin and spread route of pyrimethamine-resistant Plasmodium falciparum in southern China

Southeast Asia (the Thailand-Cambodia border) has been considered the primal epicenter for most antimalarial drug resistance; however, numerous molecular epidemiological studies have successively reported multiple independent origins of sulfadoxine-pyrimethamine (SP) resistance-associated Plasmodium falciparum dhfr (pfdhfr) and pfdhps alleles in other areas. To better understand the origin and evolutionary pathway of the SP resistance in Southeast Asia, a total of 374 P. falciparum field isolates from the Yunnan-Burma border and Hainan Island in southern China have been collected for comprehensive investigations on the mutation patterns of the pfdhfr/pfdhps genes as well as their microsatellite haplotypes. By comparative analysis of single-nucleotide polymorphism (SNP) genotyping and flanking microsatellite haplotypes, we reveal a unique origin of pyrimethamine-resistant mutations in Pfdhfr gene in Hainan Island and an oriented spread route of the pyrimethamine resistance from the Thailand-Cambodia border into the Hainan area, which reflects the geographical traits and SP administration histories in the two geographically independent areas. Moreover, genetic linkages between the high-level SP resistance-conferring pfdhfr/pfdhps alleles have been established in the isolates from the Yunnan-Burma border, raising the concern of a genetic basis in adopting combination chemotherapies against falciparum malaria.

Selection of pfdhfr/pfdhps alleles and declining artesunate/sulphadoxine-pyrimethamine efficacy against Plasmodium falciparum eight years after deployment in eastern Sudan

Background

Artesunate/sulphadoxine-pyrimethamine (AS/SP) has been the first-line treatment for falciparum malaria in Sudan since 2004. The impact of this combination on anti-malarial resistance-associated molecular markers has not been investigated. In this study, an evaluation of the efficacy and prevalence of drug resistance alleles (pfcrt, pfmdr1, pfdhfr and pfdhps) eight years after the adoption of AS/SP in eastern Sudan is reported.

Methods

A 28-day follow-up efficacy trial of AS/SP was conducted in eastern Sudan during the 2012 transmission season. Blood smears were collected from patients on days 0, 1, 2, 3, 7, 14, 21 and 28. Blood spots on filter paper were obtained pre-treatment and on the day the patient was parasite positive by microscopy. Genotyping of alleles was performed by qPCR (pfcrt 72–76 and pfmdr1 copy number) and direct sequencing of pfmdr1, pfdhfr and pfdhps.

Results

Sixty-three patients out of 68 (93%) completed the 28-day follow-up, adequate clinical, and parasitological response occurred in 90.5% and 85.3% of the patients in the per-protocol and intent-to-treat analyses, respectively. PCR corrected per-protocol efficacy was 93.7%. The enrolment prevalence of pfcrt-CVMNK was 30.2% and pfmdr1-N86 was 40.3%. The pfmdr1 haplotype NFD occurred in 32.8% of pre-treatment samples and was significantly higher than previous reports (Fisher’s exact p = 0.0001). The pfdhfr-51I/108N combination occurred in all sequenced isolates and 59R was observed in a single individual. pfdhps substitutions 436A, 437G, 540E, 581G and 613S were observed at 7.8, 77.3, 76.9%, 33.8% and 0.0%, respectively. Treatment failures were associated with the pfdhps haplotype SGEGA at these five codons (OR 7.3; 95% CI 0.65 - 368; p = 0.048).

Conclusion

The decrease of CQR associated genotypes reflects the formal policy of complete removal of CQ in Sudan. However, the frequency of markers associated with SP failure is increasing in this study area and may be contributing to the treatment efficacy falling below 90%. Further monitoring of AS/SP efficacy and of post-treatment selection of pfdhfr and pfdhps alleles in vivo is required to inform future treatment guidelines.

Population genetics analysis during the elimination process of Plasmodium falciparum in Djibouti

Background

Case management of imported malaria within the context of malaria pre-elimination is increasingly considered to be relevant because of the risk of resurgence. The assessment of malaria importation would provide key data i) to select countries with propitious conditions for pre-elimination phase and ii) to predict its feasibility. Recently, a sero-prevalence study in Djibouti indicated low malaria prevalence, which is propitious for the implementation of pre-elimination, but data on the extent of malaria importation remain unknown.

Methods

Djiboutian plasmodial populations were analysed over an eleven-year period (1998, 1999, 2002 and 2009). The risk of malaria importation was indirectly assessed by using plasmodial population parameters. Based on 5 microsatellite markers, expected heterozygosity (H.e.), multiplicity of infection, pairwise Fst index, multiple correspondence analysis and individual genetic relationship were determined. The prevalence of single nucleotide polymorphisms associated with pyrimethamine resistance was also determined.

Results

Data indicated a significant decline in genetic diversity (0.51, 0.59, 0.51 and 0 in 1998, 1999, 2002 and 2009, respectively) over the study period, which is inconsistent with the level of malaria importation described in a previous study. This suggested that Djiboutian malaria situation may have benefited from the decline of malaria prevalence that occurred in neighbouring countries, in particular in Ethiopia. The high Fst indices derived from plasmodial populations from one study period to another (0.12 between 1999 and 2002, and 0.43 between 2002 and 2009) suggested a random sampling of parasites, probably imported from neighbouring countries, leading to oligo-clonal expansion of few different strains during each transmission season. Nevertheless, similar genotypes observed during the study period suggested recurrent migrations and imported malaria.

Conclusion

In the present study, the extent of genetic diversity was used to assess the risk of malaria importation in the low malaria transmission setting of Djibouti. The molecular approach highlights i) the evolution of Djiboutian plasmodial population profiles that are consistent and compatible with Djiboutian pre-elimination goals and ii) the necessity to implement the monitoring of plasmodial populations and interventions at the regional scale in the Horn of Africa to ensure higher efficiency of malaria control and elimination.

Sulfadoxine-pyrimethamine resistance in Plasmodium falciparum: a zoomed image at the molecular level within a geographic context

Antimalarial chemotherapy is one of the main pillars in the prevention and control of malaria. Following widespread resistance of Plasmodium falciparum to chloroquine, sulfadoxine-pyrimethamine came to the scene as an alternative to the cheap and well-tolerated chloroquine. However, widespread resistance to sulfadoxine-pyrimethamine has been documented. In vivo efficacy tests are the gold standard for assessing drug resistance and treatment failure. However, they have many disadvantages, such as influence of host immunity and drug pharmacokinetics. In vitro tests of antimalarial drug efficacy also have many technical difficulties. Molecular markers of resistance have emerged as epidemiologic tools to investigate antimalarial drug resistance even before becoming clinically evident. Mutations in P. falciparum dihydrofolate reductase and dihydrofolate synthase have been extensively studied as molecular markers for resistance to pyrimethamine and sulfadoxine, respectively. This review highlights the resistance of P. falciparum at the molecular level presenting both supporting and opposing studies on the utility of molecular markers.

Plasmodium falciparum sulfadoxine resistance is geographically and genetically clustered within the DR Congo

Understanding the spatial clustering of Plasmodium falciparum populations can assist efforts to contain drug-resistant parasites and maintain the efficacy of future drugs. We sequenced single nucleotide polymorphisms (SNPs) in the dihydropteroate synthase gene (dhps) associated with sulfadoxine resistance and 5 microsatellite loci flanking dhps in order to investigate the genetic backgrounds, genetic relatedness, and geographic clustering of falciparum parasites in the Democratic Republic of the Congo (DRC). Resistant haplotypes were clustered into subpopulations: one in the northeast DRC, and the other in the balance of the DRC. Network and clonal lineage analyses of the flanking microsatellites indicate that geographically-distinct mutant dhps haplotypes derive from separate lineages. The DRC is therefore a watershed for haplotypes associated with sulfadoxine resistance. Given the importance of central Africa as a corridor for the spread of antimalarial resistance, the identification of the mechanisms of this transit can inform future policies to contain drug-resistant parasite strains.

Local population structure of Plasmodium: impact on malaria control and elimination

BACKGROUND

Regardless of the growing interest in detecting population structures in malarial parasites, there have been limited discussions on how to use this concept in control programmes. In such context, the effects of the parasite population structures will depend on interventions' spatial or temporal scales. This investigation explores the problem of identifying genetic markers, in this case microsatellites, to unveil Plasmodium genetic structures that could affect decisions in the context of elimination. The study was performed in a low-transmission area, which offers a good proxy to better understand problems associated with surveillance at the final stages of malaria elimination.

METHODS

Plasmodium vivax samples collected in Tumeremo, Venezuela, between March 2003 and November 2004 were analysed. Since Plasmodium falciparum also circulates in many low endemic areas, P. falciparum samples from the same locality and time period were included for comparison. Plasmodium vivax samples were assayed for an original set of 25 microsatellites and P. falciparum samples were assayed for 12 microsatellites.

RESULTS

Not all microsatellite loci assayed offered reliable local data. A complex temporal-cluster dynamics is found in both P. vivax and P. falciparum. Such dynamics affect the numbers and the type of microsatellites required for identifying individual parasites or parasite clusters when performing cross-sectional studies. The minimum number of microsatellites required to differentiate circulating P. vivax clusters differs from the minimum number of hyper-variable microsatellites required to distinguish individuals within these clusters. Regardless the extended number of microsatellites used in P. vivax, it was not possible to separate all individual infections.

CONCLUSIONS

Molecular surveillance has great potential; however, it requires preliminary local studies in order to properly interpret the emerging patterns in the context of elimination. Clonal expansions and clusters turnovers need to be taken into account when using molecular markers. Those affect the number and type of microsatellite markers, as well as, the expected genetic patterns in the context of operational investigations. By considering the local dynamics, elimination programmes could cost-effectively use molecular markers. However, population level studies need to consider the local limitations of a given set of loci in terms of providing epidemiologically relevant information.

Source of drug resistant Plasmodium falciparum in a potential malaria elimination site in Saudi Arabia

A major challenge to the success of malaria control program in Saudi Arabia is the high influx of expatriates and holy visitors from malaria endemic countries. In the present study we examined whether drug resistant parasite genotypes reported in Jazan region, southwest of Saudi Arabia are imported or developed locally. We examined 178 Plasmodium falciparum isolates for alleles of dihydropteroate synthase (dhps) and dihydrofolate reductase (dhfr), associated with Sulfadoxine-Pyrimethamine (SP) resistance, and three microsatellites flanking each gene. In addition, we examined a neutral polymorphic gene (Pfg377). We compared the dhfr and dhps haplotypes in Jazan, using network analysis, to an existing similar data set of 94 P. falciparum isolates from eastern Sudan. In Jazan, double mutant dhfr allele (51I, 108N) occurred with a prevalence of 33%. The vast majority (99%) of dhps were wild-type alleles. The mean expected heterozygosity (H(e)) of microsatellites around mutant dhfr alleles (H(e)=0.312; n=60) was lower (P ≤ 0.05) than that around the wild-type allele (H(e)=0.834; n=116). Also, the mutant dhfr isolates showed high H(e) for dhps (H(e)=0.80) and the non-drug resistance locus Pfg377 (H(e)=0.63) indicative of selection for mutant dhfr only. The predominant double mutant dhfr haplotype in Jazan (73%), was prevalent among P. falciparum in east Africa. Network analysis suggests the mutant haplotype of dhfr gene was possibly introduced into Jazan from East Africa. The absence of mutations in dhps as well as triple mutant dhfr haplotype associated with SP failure support the current use of SP as a partner with artesunate as a first line therapy in Saudi Arabia. However, the close relationship between the major mutant dhfr haplotype in Sudan and Saudi Arabia, favour the hypothesis of recent migration as a source of the major resistant dhfr lineage. Thus, regular monitoring of the dhfr and dhps haplotypes is of high priority to guard possible importation of high level SP resistant lineages.

Differences in selective pressure on dhps and dhfr drug resistant mutations in western Kenya

Background

Understanding the origin and spread of mutations associated with drug resistance, especially in the context of combination therapy, will help guide strategies to halt and prevent the emergence of resistance. Unfortunately, studies have assessed these complex processes when resistance is already highly prevalent. Even further, information on the evolutionary dynamics leading to multidrug-resistant parasites is scattered and limited to areas with low or seasonal malaria transmission. This study describes the dynamics of strong selection for mutations conferring resistance against sulphadoxine-pyrimethamine (SP), a combination therapy, in western Kenya between 1992 and 1999, just before SP became first-line therapy (1999). Importantly, the study is based on longitudinal data, which allows for a comprehensive analysis that contrasts with previous cross-sectional studies carried out in other endemic regions.

Methods

This study used 236 blood samples collected between 1992 and 1999 in the Asembo Bay area of Kenya. Pyrosequencing was used to determine the alleles of dihydrofolate reductase (dhfr) and dihydropterote synthase (dhps) genes. Microsatellite alleles spanning 138 kb around dhfr and dhps, as well as, neutral markers spanning approximately 100 kb on chromosomes 2 and 3 were characterized.

Results

By 1992, the South-Asian dhfr triple mutant was already spreading, albeit in low frequency, in this holoendemic Kenyan population, prior to the use of SP as a first-line therapy. Additionally, dhfr triple mutant alleles that originated independently from the predominant Southeast Asian lineage were present in the sample set. Likewise, dhps double mutants were already present as early as 1992. There is evidence for soft selective sweeps of two dhfr mutant alleles and the possible emergence of a selective sweep of double mutant dhps alleles between 1992 and 1997. The longitudinal structure of the dataset allowed estimation of selection pressures on various dhfr and dhps mutants relative to each other based on a theoretical model tailored to P. falciparum. The data indicate that drug selection acted differently on the resistant alleles of dhfr and dhps, as evidenced by fitness differences. Thus a combination drug therapy such as SP, by itself, does not appear to select for "multidrug"-resistant parasites in areas with high recombination rate.

Conclusions

The complexity of these observations emphasizes the importance of population-based studies to evaluate the effects of strong drug selection on Plasmodium falciparum populations.

Identification of pyrimethamine- and chloroquine-resistant Plasmodium falciparum in Africa between 1984 and 1998: genotyping of archive blood samples

BACKGROUND

Understanding the geographical distribution of drug resistance of Plasmodium falciparum is important for the effective treatment of malaria. Drug resistance has previously been inferred mainly from records of clinical resistance. However, clinical resistance is not always consistent with the parasite's genetic resistance. Thus, molecular identification of the parasite's drug resistance is required. In Africa, clinical resistance to pyrimethamine (Pyr) and chloroquine (CQ) was evident before 1980 but few studies investigating the genetic resistance to these drugs were conducted before the late 1990s. In this study, genotyping of genes involved in resistance to Pyr and CQ was performed using archive blood samples from Africa between 1984 and 1998.

METHODS

Parasite DNA was extracted from P. falciparum-infected blood smears collected from travellers returning to Japan from Africa between 1984 and 1998. Genotypes of the dihydrofolate reductase gene (dhfr) and CQ-resistance transporter gene (pfcrt) were determined by polymerase chain reaction amplification and sequencing.

RESULTS

Genotyping of dhfr and pfcrt was successful in 59 and 80 samples, respectively. One wild-type and seven mutant dhfr genotypes were identified. Three dhfr genotypes lacking the S108N mutation (NRSI, ICSI, IRSI; amino acids at positions 51, 59, 108, and 164 with mutations underlined) were highly prevalent before 1994 but reduced after 1995, accompanied by an increase in genotypes with the S108N mutation. The dhfr IRNI genotype was first identified in Nigeria in 1991 in the present samples, and its frequency gradually increased. However, two double mutants (ICNI and NRNI), the latter of which was exclusively found in West Africa, were more frequent than the IRNI genotype. Only two pfcrt genotypes were found, the wild-type and a Southeast Asian type (CVIET; amino acids at positions 72-76 with mutations underlined). The CVIET genotype was already present as early as 1984 in Tanzania and Nigeria, and appeared throughout Africa between 1984 and 1998.

CONCLUSIONS

This study is the first to report the molecular identification of Pyr- and CQ-resistant genotypes of P. falciparum in Africa before 1990. Genotyping of dhfr and pfcrt using archive samples has revealed new aspects of the evolutionary history of Pyr- and CQ-resistant parasites in Africa.

The evolution of pyrimethamine resistant dhfr in Plasmodium falciparum of south-eastern Tanzania: comparing selection under SP alone vs SP+artesunate combination

Background

Sulphadoxine-pyrimethamine (SP) resistance is now widespread throughout east and southern Africa and artemisinin compounds in combination with synthetic drugs (ACT) are recommended as replacement treatments by the World Health Organization (WHO). As well as high cure rates, ACT has been shown to slow the development of resistance to the partner drug in areas of low to moderate transmission. This study looked for evidence of protection of the partner drug in a high transmission African context. The evaluation was part of large combination therapy pilot implementation programme in Tanzania, the Interdisciplinary Monitoring Programme for Antimalarial Combination Therapy (IMPACT-TZ)

Methods

The growth of resistant dhfr in a parasite population where SP Monotherapy was the first-line treatment was measured for four years (2002-2006), and compared with the development of resistant dhfr in a neighbouring population where SP + artesunate (SP+AS) was used as the first-line treatment during the same interval. The effect of the differing treatment regimes on the emergence of resistance was addressed in three ways. First, by looking at the rate of increase in frequency of pre-existing mutant dhfr alleles under monotherapy and combination therapy. Second, by examining whether de-novo mutant alleles emerged under either treatment. Finally, by measuring diversity at three dhfr flanking microsatellite loci upstream of the dhfr gene.

Results

The reduction in SP selection pressure resulting from the adoption of ACT slowed the rate of increase in the frequency of the triple mutant resistant dhfr allele. Comparing between the two populations, the higher levels of genetic diversity in sequence flanking the dhfr triple mutant allele in the population where the ACT regimen had been used indicates the reduction in SP selection pressure arising from combination therapy.

Conclusion

The study demonstrated that, alleles containing two mutations at the dhfr have arisen at least four times independently while those containing triple mutant dhfr arose only once, and were found carrying a single unique Asian-type flanking sequence, which apparently drives the spread of pyrimethamine resistance associated dhfr alleles in east Africa. SP+AS is not recommended for use in areas where SP cure rates are less than 80% but this study reports an observed principle of combination protection from an area where pyrimethamine resistance was already high.

The evolution of drug resistance and the curious orthodoxy of aggressive chemotherapy

The evolution of drug-resistant pathogens is a major challenge for 21st century medicine. Drug use practices vigorously advocated as resistance management tools by professional bodies, public health agencies, and medical schools represent some of humankind's largest attempts to manage evolution. It is our contention that these practices have poor theoretical and empirical justification for a broad spectrum of diseases. For instance, rapid elimination of pathogens can reduce the probability that de novo resistance mutations occur. This idea often motivates the medical orthodoxy that patients should complete drug courses even when they no longer feel sick. Yet "radical pathogen cure" maximizes the evolutionary advantage of any resistant pathogens that are present. It could promote the very evolution it is intended to retard. The guiding principle should be to impose no more selection than is absolutely necessary. We illustrate these arguments in the context of malaria; they likely apply to a wide range of infections as well as cancer and public health insecticides. Intuition is unreliable even in simple evolutionary contexts; in a social milieu where in-host competition can radically alter the fitness costs and benefits of resistance, expert opinion will be insufficient. An evidence-based approach to resistance management is required.

Multiple origins of Plasmodium falciparum dihydropteroate synthetase mutant alleles associated with sulfadoxine resistance in India

With the spread of chloroquine (CQ)-resistant malaria in India, sulfadoxine-pyrimethamine (SP) alone or in combination with artesunate is used as an alternative antimalarial drug. Due to continuous drug pressure, the Plasmodium falciparum parasite is exhibiting resistance to antifolates because of mutations in candidate genes dihydrofolate reductase (dhfr) and dihydropteroate synthetase (dhps). Our earlier study on flanking microsatellite markers of dhfr mutant alleles from India had shown a single origin of the pyrimethamine resistance and some minor haplotypes which shared haplotypes with Southeast Asian (Thailand) strains. In the present study, we have analyzed 193 of these Indian P. falciparum isolates for 15 microsatellite loci around dhps to investigate the genetic lineages of the mutant dhps alleles in different parts of the country. Eighty-one of these samples had mutant dhps alleles, of which 62 were from Andaman and Nicobar Islands and the remaining 19 were from mainland India. Of 112 isolates with a wild-type dhps allele, 109 were from mainland India and only 3 were from Andaman and Nicobar Islands. Consistent with the model of selection, the mean expected heterozygosity (H(e)) around mutant dhps alleles (H(e) = 0.55; n = 81) associated with sulfadoxine resistance was lower (P ≤ 0.05) than the mean H(e) around the wild-type dhps allele (H(e) = 0.80; n = 112). There was more genetic diversity in flanking microsatellites of dhps than dhfr among these isolates, which confirms the assertion that dhps mutations are at a very early stage of fixation in the parasite population. Microsatellite haplotypes around various mutant dhps alleles suggest that the resistant dhps alleles have multiple independent origins in India, especially in Andaman and Nicobar Islands. Determining the genetic lineages of the resistant dhps alleles on Andaman and Nicobar Islands and mainland India is significant, given the role of Asia in the intercontinental spread of chloroquine- and pyrimethamine-resistant parasites in the past.

Selective sweeps and genetic lineages of Plasmodium falciparum drug -resistant alleles in Ghana

BACKGROUND

In 2005, Ghana adopted artemisinin-based combination therapy (ACT) for primary treatment of falciparum malaria. A comprehensive study of the drug-resistance-associated mutations and their genetic lineages will lead to a better understanding of the evolution of antimalarial drug resistance in this region.

METHODS

The pfcrt, pfmdr1, dhps, and dhfr mutations associated with chloroquine (CQ) and sulfadoxine-pyrimethamine (SP) resistance and the microsatellite loci flanking these genes were genotyped in Plasmodium falciparum isolates from Ghana.

RESULTS

The prevalence of mutations associated with both CQ and SP resistance was high in Ghana. However, we observed a decrease in prevalence of the pfcrt K76T mutation in northern Ghana after the change in drug policy from CQ to ACT. Analysis of genetic diversity and differentiation at microsatellite loci flanking all 4 genes indicated that they have been under strong selection, because of CQ and SP use. The triple-mutant pfcrt and dhfr alleles in Ghana were derived from Southeast Asia, whereas the double-mutant dhfr, dhps, and pfmdr1 alleles were of African lineage.

CONCLUSION

Because of the possible role of pfmdr1 in amodiaquine and mefloquine resistance, demonstrating selection on pfmdr1 and defining lineages of resistant alleles in an African population holds great importance.

Whole genome re-sequencing identifies a mutation in an ABC transporter (mdr2) in a Plasmodium chabaudi clone with altered susceptibility to antifolate drugs

In malaria parasites, mutations in two genes of folate biosynthesis encoding dihydrofolate reductase (dhfr) and dihydropteroate synthase (dhps) modify responses to antifolate therapies which target these enzymes. However, the involvement of other genes which modify the availability of exogenous folate, for example, has been proposed. Here, we used short-read whole-genome re-sequencing to determine the mutations in a clone of the rodent malaria parasite, Plasmodium chabaudi, which has altered susceptibility to both sulphadoxine and pyrimethamine. This clone bears a previously identified S106N mutation in dhfr and no mutation in dhps. Instead, three additional point mutations in genes on chromosomes 2, 13 and 14 were identified. The mutated gene on chromosome 13 (mdr2 K392Q) encodes an ABC transporter. Because Quantitative Trait Locus analysis previously indicated an association of genetic markers on chromosome 13 with responses to individual and combined antifolates, MDR2 is proposed to modulate antifolate responses, possibly mediated by the transport of folate intermediates.

Anti-folate drug resistance in Africa: meta-analysis of reported dihydrofolate reductase (dhfr) and dihydropteroate synthase (dhps) mutant genotype frequencies in African Plasmodium falciparum parasite populations

BACKGROUND

Mutations in the dihydrofolate reductase (dhfr) and dihydropteroate synthase (dhps) genes of Plasmodium falciparum are associated with resistance to anti-folate drugs, most notably sulphadoxine-pyrimethamine (SP). Molecular studies document the prevalence of these mutations in parasite populations across the African continent. However, there is no systematic review examining the collective epidemiological significance of these studies. This meta-analysis attempts to: 1) summarize genotype frequency data that are critical for molecular surveillance of anti-folate resistance and 2) identify the specific challenges facing the development of future molecular databases.

METHODS

This review consists of 220 studies published prior to 2009 that report the frequency of select dhfr and dhps mutations in 31 African countries. Maps were created to summarize the location and prevalence of the highly resistant dhfr triple mutant (N51I, C59R, S108N) genotype and dhps double mutant (A437G and K540E) genotype in Africa. A hierarchical mixed effects logistic regression was used to examine the influence of various factors on reported mutant genotype frequency. These factors include: year and location of study, age and clinical status of sampled population, and reporting conventions for mixed genotype data.

RESULTS

A database consisting of dhfr and dhps mutant genotype frequencies from all African studies that met selection criteria was created for this analysis. The map illustrates particularly high prevalence of both the dhfr triple and dhps double mutant genotypes along the Kenya-Tanzania border and Malawi. The regression model shows a statistically significant increase in the prevalence of both the dhfr triple and dhps double mutant genotypes in Africa.

CONCLUSION

Increasing prevalence of the dhfr triple mutant and dhps double mutant genotypes in Africa are consistent with the loss of efficacy of SP for treatment of clinical malaria in most parts of this continent. Continued assessment of the effectiveness of SP for the treatment of clinical malaria and intermittent preventive treatment in pregnancy is needed. The creation of a centralized resistance data network, such as the one proposed by the WorldWide Antimalarial Resistance Network (WWARN), will become a valuable resource for planning timely actions to combat drug resistant malaria.

Approximations for the hitchhiking effect caused by the evolution of antimalarial-drug resistance

An analytically feasible, deterministic model for the spread of drug resistance among human malaria parasites, which incorporates all characteristics of the complex malaria-transmission cycle was introduced by Schneider and Kim (Theor. Popul Biol, 2010). The model accounts for the fact that only a fraction of infected hosts receive drug treatment and that hosts can be co-infected by differently many parasites. Furthermore, the model also incorporates host heterogeneity. Antimalarial-drug resistance is assumed to be caused by a single locus with two alleles-a sensitive one and a resistance one. The most important result for this model is that an analytical solution for the frequencies of a linked neutral biallelic locus exists. However, the exact solution does not admit an explicit form, and cannot straightforwardly be interpreted in terms of the model parameters. Here, we establish simple approximations for the equilibrium frequency at the neutral locus. Under the assumption that the resistant allele is initially rare-the biologically most relevant assumption in this context-and that recombination is weak, the approximations become similar to the approximations in the standard hitchhiking model. However, there are crucial differences. In particular, because of the high degree of selfing among malaria parasites in their sexual phase, a genome-wide reduction of relative heterozygosity occurs if selection is sufficiently strong. It turns out that the approximations are accurate even if the recombination rates are not small and the resistant allele is initially not very rare. The main advantage of our approximations is that they are easy to interpret in terms of model parameters. Moreover, they allow to make predictions of the size of the valley of reduced heterozygosity around the selected locus for given model parameters. Reversely, for a given reduction of heterozygosity, it is possible to identify the corresponding parameters. Moreover, we will show that incorporating host heterogeneity leads to an increased hitchhiking effect.

An analytical model for genetic hitchhiking in the evolution of antimalarial drug resistance

We analytically study a deterministic model for the spread of drug resistance among human malaria parasites. The model incorporates all major characteristics of the complex malaria transmission cycle and accounts for the fact that only a fraction alpha of infected hosts receive drug treatment. Furthermore, the model incorporates that hosts can be co-infected. The number m of parasites co-infecting a host is either a constant or, more generally, follows a given frequency distribution. Although the model is formulated in a multilocus setup, for our results we assume that drug resistance is caused by a single locus with two alleles - a sensitive one and a resistant one. We assume that the resistant allele has a selective advantage only in treated hosts and pays metabolic costs, which causes this allele to be deleterious in untreated hosts. We provide necessary and sufficient conditions for the fixation of the resistant allele. Moreover, provided the resistant allele will sweep through the population, we derive a formula for the time until it reaches a given frequency and in particular for the time until quasi-fixation. Furthermore, we establish an analytical solution for allele frequency changes at a linked neutral biallelic locus due to the rapid increase in frequency of the resistant allele. Our solution describes a local reduction in heterozygosity among parasite chromosomes around the resistant allele, the effect commonly referred to as the hitchhiking effect, as a function of alpha and m. The result therefore allows the investigation of selective sweep patterns under specific demographic settings. We find that the hitchhiking effect is similar but different from the standard model of genetic hitchhiking that assumes random mating and homogeneous selection. In particular, the process of recombination and selection cannot be decoupled. We further explain why standard hitchhiking theory cannot be applied to drug resistance in malaria. Furthermore, we will show that a genome-wide reduction in relative heterozygosity can occur provided a fraction of hosts is infected by a single parasite haplotype. Finally, we show how to incorporate host heterogeneity, and generalize our results to this biologically more realistic case.

Tracing the origins and signatures of selection of antifolate resistance in island populations of Plasmodium falciparum

Background

Resistance of the malaria parasite Plasmodium falciparum to sulfadoxine-pyrimethamine (SP) has evolved worldwide. In the archipelago of São Tomé and Principe (STP), West Africa, although SP resistance is highly prevalent the drug is still in use in particular circumstances. To address the evolutionary origins of SP resistance in these islands, we genotyped point mutations at P. falciparum dhfr and dhps genes and analysed microsatellites flanking those genes.

Methods

Blood samples were collected in July and December 2004 in three localities of São Tomé Island and one in Principe Island. Species-specific nested-PCR was used to identify P. falciparum infected samples. Subsequently, SNPs at the dhfr and dhps genes were identified through PCR-RFLP. Isolates were also analysed for three microsatellite loci flanking the dhfr gene, three loci flanking dhps and four loci located at putative neutral genomic regions.

Results

An increase of resistance-associated mutations at dhfr and dhps was observed, in particular for the dhfr/dhps quintuple mutant, associated with clinical SP failure. Analysis of flanking microsatellites suggests multiple independent introductions for dhfr and dhps mutant haplotypes, possibly from West Africa. A reduced genetic diversity and increased differentiation at flanking microsatellites when compared to neutral loci is consistent with a selective sweep for resistant alleles at both loci.

Conclusions

This study provides additional evidence for the crucial role of gene flow and drug selective pressures in the rapid spread of SP resistance in P. falciparum populations, from only a few mutation events giving rise to resistance-associated mutants. It also highlights the importance of human migration in the spread of drug resistant malaria parasites, as the distance between the islands and mainland is not consistent with mosquito-mediated parasite dispersal.

Origins and spread of pfdhfr mutant alleles in Plasmodium falciparum

The emergence and spread of Plasmodium falciparum parasite resistant to sulfadoxine and pyrimethamine (SP) poses a serious public health problem. Resistance is caused by point mutations in dihydrofolate reductase (pfdhfr) and dihydropteroate synthase (pfdhps), the two key enzymes in the folate biosynthetic pathway. The use of microsatellite markers flanking pfdhfr has recently shown that the invasion of limited resistant lineages may explain the widespread SP resistance in many endemic regions. In Africa, however, multiple indigenous origins of pfdhfr triple mutants have been demonstrated. More new independent lineages and routes of geographical spread of resistance may be found by further molecular evolutionary analyses using samples from various endemic regions. Here, I review recent studies about the history of SP usage and the evolution and spread of resistant lineages while addressing the technical issue of microsatellite analysis.

Transmission and cross-mating of high-level resistance Plasmodium falciparum dihydrofolate reductase haplotypes in The Gambia

A high-level pyrimethamine resistance Plasmodium falciparum lineage with triple dihydrofolate reductase (dhfr) mutations prevails across Africa. However, additional minority lineages were seen. We examined transmission success of mutant dhfr haplotypes among 22 children in The Gambia and 60 infected Anopheles gambiae mosquitoes fed on their blood. Additional polymorphic genes of the gametocyte-specific protein (pfg377) and merozoite surface protein-1 (MSP-1) were examined. Similarities were seen between pfg377 and MSP-1 alleles in children and mosquitoes and evidence of cross-mating between different parasite genotypes was seen in some infected mosquitoes, reflecting high transmission success of existing clones. With regard to dhfr, 16 haplotypes were seen among the children: 2 carried double mutations and 14 carried triple mutations. However, only nine haplotypes, all with triple mutations, were detected among mosquitoes. A single triple-mutant dhfr haplotype, similar to that in other countries in Africa, predominated among children (42%) and mosquitoes (60%), supporting the hypothesis of migration of this haplotype across Africa. However, evidence of cross-mating between the above haplotypes signifies the role of local evolution.

Origin and evolution of sulfadoxine resistant Plasmodium falciparum

The Thailand-Cambodia border is the epicenter for drug-resistant falciparum malaria. Previous studies have shown that chloroquine (CQ) and pyrimethamine resistance originated in this region and eventually spread to other Asian countries and Africa. However, there is a dearth in understanding the origin and evolution of dhps alleles associated with sulfadoxine resistance. The present study was designed to reveal the origin(s) of sulfadoxine resistance in Cambodia and its evolutionary relationship to African and South American dhps alleles. We sequenced 234 Cambodian Plasmodium falciparum isolates for the dhps codons S436A/F, A437G, K540E, A581G and A613S/T implicated in sulfadoxine resistance. We also genotyped 10 microsatellite loci around dhps to determine the genetic backgrounds of various alleles and compared them with the backgrounds of alleles prevalent in Africa and South America. In addition to previously known highly-resistant triple mutant dhps alleles SGEGA and AGEAA (codons 436, 437, 540, 581, 613 are sequentially indicated), a large proportion of the isolates (19.3%) contained a 540N mutation in association with 437G/581G yielding a previously unreported triple mutant allele, SGNGA. Microsatellite data strongly suggest the strength of selection was greater on triple mutant dhps alleles followed by the double and single mutants. We provide evidence for at least three independent origins for the double mutants, one each for the SGKGA, AGKAA and SGEAA alleles. Our data suggest that the triple mutant allele SGEGA and the novel allele SGNGA have common origin on the SGKGA background, whereas the AGEAA triple mutant was derived from AGKAA on multiple, albeit limited, genetic backgrounds. The SGEAA did not share haplotypes with any of the triple mutants. Comparative analysis of the microsatellite haplotypes flanking dhps alleles from Cambodia, Kenya, Cameroon and Venezuela revealed an independent origin of sulfadoxine resistant alleles in each of these regions.

Widespread resistance to chloroquine (CQ) and sulfadoxine-pyrimethamine (SP), the two least expensive and widely available antimalarial drugs, has become a major global public health challenge. It is known that point mutations in Plasmodium falciparum crt, dhfr and dhps genes contribute to resistance to CQ, pyrimethamine and sulfadoxine, respectively. CQ and pyrimethamine resistance spread to Africa and Asia from a few founding mutant lineages originating from the Thailand-Cambodia border. Here, we define the origins of dhps alleles in Cambodia and their relationships to African and South American counterparts. Three different triple mutant alleles including a novel allele comprised of 437G, 540N, and 581G mutations (S₄₃₆G₄₃₇N₅₄₀G₅₈₁A₆₁₃) were found in Cambodia as opposed to a single triple mutant allele in South America and two common double mutant alleles in Africa. Microsatellite data suggest strong selection operating on triple mutant alleles as compared to double and single mutants in Cambodia. We report three major independent origins for the double mutants and at least two independent origins for the highly resistant triple mutant dhps alleles in Cambodia. We also show that the resistant dhps alleles in Africa and South America have distinct origins from Cambodia. These results suggest that the evolution and spread of sulfadoxine resistance is different from CQ and pyrimethamine resistance.

Origins of the recent emergence of Plasmodium falciparum pyrimethamine resistance alleles in Madagascar

The combination of sulfadoxine-pyrimethamine is recommended for use as intermittent preventive treatment of malaria during pregnancy and is deployed in Africa. The emergence and the spread of resistant parasites are major threats to such an intervention. We have characterized the Plasmodium falciparum dhfr (pfdhfr) haplotypes and flanking microsatellites in 322 P. falciparum isolates collected from the Comoros Islands and Madagascar. One hundred fifty-six (48.4%) carried the wild-type pfdhfr allele, 19 (5.9%) carried the S108N single-mutation allele, 30 (9.3%) carried the I164L single-mutation allele, 114 (35.4%) carried the N51I/C59R/S108N triple-mutation allele, and 3 (1.0%) carried the N51I/C59R/S108N/I164L quadruple-mutation allele. Microsatellite analysis showed the introduction from the Comoros Islands of the ancestral pfdhfr triple mutant allele of Asian origin and its spread in Madagascar. Evidence for the emergence on multiple occasions of the I164L single-mutation pfdhfr allele in Madagascar was also obtained. Thus, the conditions required to generate mutants with quadruple mutations are met in Madagascar, representing a serious threat to current drug policy.