Molecular characterization of Plasmodium falciparum antifolate resistance markers in Thailand between 2008 and 2016.

BACKGROUND: Resistance to anti-malarials is a major threat to the control and elimination of malaria. Sulfadoxine-pyrimethamine (SP) anti-malarial treatment was used as a national policy for treatment of uncomplicated falciparum malaria in Thailand from 1973 to 1990. In order to determine whether withdrawal of this antifolate drug has led to restoration of SP sensitivity, the prevalence of genetic markers of SP resistance was assessed in historical Thai samples. 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 (varphi = 0.213, p-value = 0.001) or pfdhps K540E/N (varphi = 0.399, p-value