Abstract
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.
•Illumina MiSeq deep amplicon sequencing was used to identify pyrimethamine resistance associated SNPs in Plasmodium vivax.•Pyrimethamine resistance was identified in field populations.•Positive selection pressure for resistance in human patients was investigated.•Hard and soft selective sweep models were used for the emergence of resistance.•Both single and multiple recurrent resistance mutations were spread by migration.