METHODOLOGY/PRINCIPAL FINDINGS: Anopheles spp. were sampled using human landing catch (HLC) method at Paradason village in Kudat district of Sabah. The collected Anopheles were identified morphologically and then subjected to total DNA extraction and polymerase chain reaction (PCR) to detect Plasmodium parasites in the mosquitoes. Identification of Plasmodium spp. was confirmed by sequencing the SSU rRNA gene with species specific primers. MEGA4 software was then used to analyse the SSU rRNA sequences and bulid the phylogenetic tree for inferring the relationship between simian malaria parasites in Sabah. PCR results showed that only 1.61% (23/1,425) of the screened An. balabacensis were infected with one or two of the five simian Plasmodium spp. found in Sabah, viz. Plasmodium coatneyi, P. inui, P. fieldi, P. cynomolgi and P. knowlesi. Sequence analysis of SSU rRNA of Plasmodium isolates showed high percentage of identity within the same Plasmodium sp. group. The phylogenetic tree based on the consensus sequences of P. knowlesi showed 99.7%-100.0% nucleotide identity among the isolates from An. balabacensis, human patients and a long-tailed macaque from the same locality.
CONCLUSIONS/SIGNIFICANCE: This is the first study showing high molecular identity between the P. knowlesi isolates from An. balabacensis, human patients and a long-tailed macaque in Sabah. The other common simian Plasmodium spp. found in long-tailed macaques and also detected in An. balabacensis were P. coatneyi, P. inui, P. fieldi and P. cynomolgi. The high percentage identity of nucleotide sequences between the P. knowlesi isolates from the long-tailed macaque, An. balabacensis and human patients suggests a close genetic relationship between the parasites from these hosts.
METHODS: Plasmodium falciparum DNA was collected from the Thailand-Myanmar, Thailand-Malaysia and Thailand-Cambodia borders during 2008-2016 (N = 233). Semi-nested PCR and nucleotide sequencing were used to assess mutations in Plasmodium falciparum dihydrofolate reductase (pfdhfr), P. falciparum dihydropteroate synthase (pfdhps). Gene amplification of Plasmodium falcipaurm GTP cyclohydrolase-1 (pfgch1) was assessed by quantitative real-time PCR. The association between pfdhfr/pfdhps mutations and pfgch1 copy numbers were evaluated.
RESULTS: Mutations in pfdhfr/pfdhsp and pfgch1 copy number fluctuated overtime through the study period. Altogether, 14 unique pfdhfr-pdfhps haplotypes collectively containing quadruple to octuple mutations were identified. High variation in pfdhfr-pfdhps haplotypes and a high proportion of pfgch1 multiple copy number (51% (73/146)) were observed on the Thailand-Myanmar border compared to other parts of Thailand. Overall, the prevalence of septuple mutations was observed for pfdhfr-pfdhps haplotypes. In particular, the prevalence of pfdhfr-pfdhps, septuple mutation was observed in the Thailand-Myanmar (50%, 73/146) and Thailand-Cambodia (65%, 26/40) border. In Thailand-Malaysia border, majority of the pfdhfr-pfdhps haplotypes transaction from quadruple (90%, 9/10) to quintuple (65%, 24/37) during 2008-2016. Within the pfdhfr-pfdhps haplotypes, during 2008-2013 the pfdhps A/S436F mutation was observed only in Thailand-Myanmar border (9%, 10/107), while it was not identified later. In general, significant correlation was observed between the prevalence of pfdhfr I164L (ϕ = 0.213, p-value = 0.001) or pfdhps K540E/N (ϕ = 0.399, p-value ≤ 0.001) mutations and pfgch1 gene amplification.
CONCLUSIONS: Despite withdrawal of SP as anti-malarial treatment for 17 years, the border regions of Thailand continue to display high prevalence of antifolate and anti-sulfonamide resistance markers in falciparum malaria. Significant association between pfgch1 amplification and pfdhfr (I164L) or pfdhps (K540E) resistance markers were observed, suggesting a compensatory mutation.