Skip to main content
Download PDF

Phenotypic and genotypic characterization of Thai isolates of Plasmodium falciparum after an artemisinin resistance containment project

Malaria Journal volume 17, Article number: 197 (2018) Cite this article

Abstract

Background

In Thailand, artemisinin-based combination therapy (ACT) has been used to treat uncomplicated falciparum malaria since 1995. Unfortunately, artemisinin resistance has been reported from Thailand and other Southeast Asian countries since 2003. Malarone®, a combination of atovaquone–proguanil (ATQ–PG), has been used to cease artemisinin pressure in some areas along Thai–Cambodia border, as part of an artemisinin resistance containment project since 2009. This study aimed to determine genotypes and phenotypes of Plasmodium falciparum isolates collected from the Thai–Cambodia border after the artemisinin resistance containment project compared with those collected before.

Results

One hundred and nine of P. falciparum isolates collected from Thai–Cambodia border from Chanthaburi and Trat provinces during 1988–2016 were used in this study. Of these, 58 isolates were collected after the containment. These parasite isolates were characterized for in vitro antimalarial sensitivities including chloroquine (CQ), quinine (QN), mefloquine (MQ), piperaquine (PPQ), artesunate (AS), dihydroartemisinin (DHA), ATQ and PG and genetic markers for drug resistance including the Kelch13 (k13), Plasmodium falciparum chloroquine resistance transporter (pfcrt), P. falciparum multidrug resistance 1 (pfmdr1) and cytochrome b (cytb) genes. Mean CQ, QN, MQ, PPQ and AS IC50s of the parasite isolates collected from 2009 to 2016 exhibited significantly higher than those of parasites collected before 2009. Approximately 57% exhibited in vitro MQ resistance. Approximately 94% of the isolates collected from 2009 to 2016 contained the pfmdr1 184F allele. Mutations of the k13 gene were detected in approximately 90% of the parasites collected from 2009 to 2016 which were significantly higher than the parasite isolates collected before. No ATQ-resistant genotype and phenotype of P. falciparum were found among the isolates collected after the containment project.

Conclusions

Although the containment project had been implemented in this area, the expansion of artemisinin-resistant parasites did not decline. In addition, reduced sensitivity of the partner drugs of ACT including MQ and PPQ were identified.

Background

Malaria is one of the important parasitic diseases threatening human beings over many decades. Emergence and spread of drug resistance is an important cause of morbidity and mortality in malaria. Artemisinin derivatives are the most potent drugs against multidrug-resistant Plasmodium falciparum. Artemisinin-based combination therapy (ACT) has been recommended by the World Health Organization (WHO) to use as the first-line treatment for multidrug-resistant falciparum malaria [1]. Artesunate-mefloquine (AS-MQ) had been used in Thailand to treat uncomplicated falciparum malaria since 1995 due to the emergence of MQ resistance [2]. Later AS–MQ was used in other Southeast Asian countries including Cambodia [2]. Unfortunately, the efficacy of AS–MQ declined along the Thai-Cambodian border after a few years of implementation [3, 4]. Increased AS–MQ failure rates observed in Thailand and Cambodia were usually associated with MQ resistance [5, 6]. Recently, a few studies have shown evidence of artemisinin resistance in P. falciparum, defined as delayed parasite clearance, i.e., presence of parasitaemia on day 3 following treatment with AS monotherapy or ACT [7, 8]. Subsequently, artemisinin resistance has emerged independently and spread to many areas of the Greater Mekong Subregion (GMS) [9, 10]. To determine the situation of artemisinin resistance, a novel ring survival assays (RSA) (i.e., in vitro or ex vivo RSA) has been used to represent the delayed parasite clearance phenotype [11]. Ariey et al. have identified mutations of the Kelch13 (k13) gene as molecular markers for artemisinin resistance [12]. Mutations in the k13 gene were correlated with delayed parasite clearance and also increased RSA (0–3 h) rate [12, 13]. The mutations were associated with artemisinin resistance in Cambodia and other countries in the GMS [12, 14,15,16].

Due to the emergence of artemisinin resistance in some areas along the Thai-Cambodian border, an artemisinin resistance containment project has been launched to cease the spread of artemisinin resistance since 2009 [17]. The combination of AS–MQ has been replaced by a fixed-dose combination of atovaquone–proguanil (ATQ–PG) (Malarone®) to reduce artemisinin pressure in these areas. ATQ inhibits the mitochondrial electron transport chain at the bc1 complex [18, 19]. Mutations in the cytochrome b (cytb) gene resulted in amino acid changes at codon 268, exchanging tyrosine for serine (Y268S) or, less frequently, asparagine (Y268 N) conferred ATQ resistance and treatment failure [20, 21]. Monitoring in vitro drug sensitivities and also molecular markers of P. falciparum isolates is essential to detect the emergence and spread of drug resistance and provide valuable information for a rational drug use policy. Little is known about the characteristics of P. falciparum after the artemisinin resistance containment project. This study was aimed to determine phenotypes and genotypes of P. falciparum isolated from the Thai–Cambodian border after the artemisinin resistance containment project compared with those collected before.

Methods

Parasite isolates and cultivation

One hundred and nine P. falciparum adapted isolates from along the Thai-Cambodia border of Chanthaburi and Trat Provinces from 1988 to 2016, were used in this study. The research protocol was reviewed and approved by the Ethics Committee of the Royal Thai Army Medical Department. Of these, 51 frozen adapted isolates were collected between 1988 and 2005 from Chanthaburi province. Thirty-six isolates was collected from 1988 to 1993, while fifteen isolates were collected between 2003 and 2005. Fifty-eight newly adapted isolates were collected from the border of Trat province from 2009 to 2016, when the artemisinin resistance containment project was implemented. No clinical profiles and outcomes of the patients were recorded. All isolates were cultured using the modification method of Trager and Jensen [22]. Parasite isolates were cultured in culture flasks containing medium RPMI 1640 (Gibco®), 10% human AB serum and human erythrocytes (O+). Culture flasks were maintained under an atmosphere of 90% N2, 5% O2, and 5% CO2 and incubated at 37 °C.

In vitro drug sensitivity assay

Newly adapted isolates were usually maintained in the culture for 2–3 cycles before in vitro sensitivity assays. Plasmodium falciparum isolates collected before and after 2009 were randomly measured for the IC50 at the same time to minimize batch effects. In addition, the reference laboratory strains, K1 and 3D7 were used as a control in each batch of experiment. Drug sensitivities of P. falciparum isolates to chloroquine (CQ), quinine (QN), MQ, piperaquine (PPQ), AS, dihydroartemisinin (DHA), ATQ and PG were determined by measuring [3H] hypoxanthine incorporated in parasite nucleic acids using the modified technique of Desjadins et al. [23]. The GRAFIT® Program (Erithacus Software Limited, UK) was used to determine inhibitory concentration 50% (IC50).

Characterization of the pfcrt, pfmdr1, cytb and k13 genes

Plasmodium falciparum DNA was extracted using the Chelex-resin method [24]. PCR–RFLP techniques were performed to detect pfcrt mutation at codon 76, pfmdr1 mutation at codons 86, 184, 1034, 1042 and 1246 and cytb mutation at codon 268 as previously described [20, 25, 26]. The pfmdr1 gene copy number was detected by TaqMan real-time PCR (CFX96™ Real-Time PCR Detection System; Bio-Rad Laboratories, Inc.) [27]. The K1 and DD2 clone containing 1 and 4 pfmdr1 copies, respectively, were used as the reference DNA sample. The amplification of pfmdr1 and β-tubulin was performed in triplicate and relative pfmdr1 copy number was determined [27].

The k13 gene was amplified using the method of Ariey et al. [12]. PCR products were performed and visualized using 2% agarose gel electrophoresis and SYBR® Safe DNA Gel Stain (Thermo Fisher Scientific). The PCR products were purified using the QIAquick® Gel Extraction Kit and sequenced by U2Bio Co. Ltd. (Seoul, South Korea). The sequences were aligned against the P. falciparum 3D7 strain (PF3D7_1343700, PlasmoDB Release 28) using BioEdit v7.2.5.

Statistical analysis

Data were analysed by STATA/MP, Version 12. Each IC50 value represented the mean of at least three independent experiments. Each in vitro drug sensitivity experiment was performed in triplicate. Normally distributed IC50 data were assessed by the Kolmogorov–Smirnov test. Differences of the mean IC50 and copy number of the pfmdr1 gene among groups were analysed by Independent t test and One-way ANOVA. Post hoc test (Scheffe) for multiple comparisons was used to test for differences between the two groups. Associations between genotypes and P. falciparum from different areas were analysed by Chi square or Fisher’s exact test and level of significance was set at a p value of < 0.05.

Results

In vitro drug sensitivities of Plasmodium falciparum

In vitro drug sensitivities of P. falciparum isolates collected before and after 2009 are presented in Table 1. Mean CQ, QN, MQ, PPQ and AS IC50s of the parasite isolates collected from 2009 to 2016 showed significantly higher value than those of parasites collected before 2009. In contrast, mean DHA IC50s of the parasites before 2009 and 2009–2016 were not significantly different. In the present study, only ATQ and PG IC50s of the parasites collected from 2009 to 2016 were shown.

Table 1 In vitro drug sensitivity of P. falciparum isolated from along the Thai–Cambodia border collected before 2009 and from 2009 to 2016
Full size table

Table 2 shows the percentage of resistant phenotypes of P. falciparum isolated along the Thai–Cambodian border in different years. Cut-off points for in vitro anti-malarial resistance were used as previously described [28,29,30,31,32]. In vitro CQ, MQ and QN resistance were defined when the IC50 ≥ 25 nM, > 30 nM and > 800 nM, respectively [28,29,30]. However, according to Pradines et al., AS and DHA resistance have been categorized when the IC50 was more than 10.5 nM [31]. In vitro ATQ resistance was defined when the parasites exhibited the IC50 of more than 1900 nM [32]. The cut-off point of PPQ resistance has not been defined yet. Nearly all parasite isolates exhibited in vitro CQ resistance. Parasite isolates collected in 2009 contained significantly higher MQ-resistant phenotype (56.9%) than the parasites before 2009 (22.0%) (p < 0.001, Fisher’s Exact test). No parasite isolates showing QN, AS, DHA or ATQ resistance were identified.

Table 2 The percentage of resistant phenotypes of P. falciparum isolated from along the Thai–Cambodia border collected before 2009 and from 2009 to2016
Full size table

Polymorphisms of pfcrt, pfmdr1, cytb and k13 genes

Plasmodium falciparum isolates collected in the years before and after the artemisinin resistance containment project from Chanthaburi and Trat Provinces were characterized for the polymorphisms of the pfcrt, pfmdr1, cytb and k13 genes (Table 3). Of 104 isolates, only two contained the wild-type pfcrt gene. Approximately 89% of these parasites had the pfmdr1 184F allele. The copy number of pfmdr1 had approximately 1 copy number (1.1 ± 0.6, n = 102) as shown in Table 3. The wild-type cytb gene at the codon 268 was identified in all parasite isolates both before and after the artemisinin containment project.

Table 3 Resistant gene polymorphisms of the parasite isolates collected before 2009 and from 2009 to 2016
Full size table

Single nucleotide polymorphisms (SNPs) on the k13 gene of the parasite isolates were determined by PCR amplification and sequencing. The percentage of the k13 mutation(s) was significantly higher among the parasites collected from 2009 to 2016 (44/52, 84.6%) compared with those collected before 2009 (9/49, 18.4%) (p < 0.001, Fisher’s Exact test). Seven different SNPs were identified including G436S, F483L, Y493H, G538V, R539T, V568G and C580Y (Table 3). The most common SNP found among these isolates were C580Y (38.6%). The k13-R539T allele, another common SNP in the Great Mekong Subregion, was identified in 8.9% of these isolates. Two mutations on the k13 gene were identified in one isolate collected from 2009 to 2016 which were combinations of C580Y/F483L. The percentage of both C580Y and R539T had significantly increased (p = 0.019 and p < 0.001, respectively, Fisher’s Exact test) among the parasites collected from 2009 to 2016.

The parasite isolates were categorized in seven groups according to their genotype of the pfmdr1 and K13 genes (Table 4), i.e., (I) the pfmdr1 86N allele with no mutation on the k13 gene, (II) the pfmdr1 184F allele with no mutation on the k13 gene, (III) the pfmdr1 184F allele with k13 580Y alleles (IV) the pfmdr1 184F allele with k13 539T alleles (V) the pfmdr1 184F + 1042D alleles with no mutation on the k13 gene (VI) the pfmdr1 184F + 1034C + 1042D alleles with no mutation on the k13 gene and (VII) others. Before 2009, the most prevalent parasites were those containing the pfmdr1 184F allele with no mutation on the k13 gene (20/49, 40.8%). Parasites with the pfmdr1 86Y allele or the pfmdr1 184F + 1034C + 1042N alleles with no mutation on the k13 gene were also identified before 2009, but not among the parasites collected from 2009 to 2016. Before 2009, the category VII contained one isolate with no mutation on both genes, with the pfmdr1 184F + 1042D alleles and k13 580Y alleles and with the pfmdr1 184F and k13 483L alleles were identified. Approximately 57% of parasites collected from 2009 to 2016 contained the pfmdr1 184F and k13 580Y alleles. In addition, 13% of these parasites exhibited the pfmdr1 184F and k13 539T alleles. In the group collected from 2009 to 2016, category VII consisted of one isolate with no mutation on both genes, two isolates with the k13 580Y allele but no mutation on the pfmdr1 gene, one isolate with the pfmdr1 184F + 1034C and the k13 580Y alleles, one isolate with the pfmdr1 184F + 1034C and the K1 539T alleles. In addition, the pfmdr1 184F allele with the k13 436S, 493H and 538 V alleles was found in one isolate each.

Table 4 Haplotypes of Plasmodium falciparum isolates collected before 2009 and from 2009 to 2016
Full size table

Correlation of anti-malarial drug sensitivities and genetic polymorphisms of Plasmodium falciparum

Table 5 shows in vitro anti-malarial sensitivities among the parasite isolates with different pfmdr1 genotypes. The parasites containing the pfmdr1 86Y or 1034C or 1042D alleles exhibited significantly increased MQ sensitivity (p < 0.001, Independent t test) while, those containing the pfmdr1 184F allele showed significantly reduced MQ sensitivity (p = 0.043, Independent t test). The parasites with ≤ 1 copy of the pfmdr1 gene showed higher MQ IC50 than those with more than one copy (p < 0.018, Independent t test). The parasites containing the pfmdr1 86Y allele exhibited significantly increased QN sensitivity compared with the wild-type counterpart (p = 0.008, Independent t test). The parasites with the pfmdr1 184F allele exhibited approximately twice less susceptible to QN than the parasites with the pfmdr1 184Y allele (p = 0.003, Independent t test). The parasites containing the pfmdr1 N1042 allele or the k13 580Y alleles exhibited significantly higher PPQ IC50 compared with the others (p = 0.019 and p = 0.013, respectively, Independent t test). The parasites containing the k13 539T alleles exhibited significant higher ARS IC50 compared with the others (p = 0.017, Independent t test).

Table 5 Correlation of antimalarial drug sensitivities and pfmdr1 and K13 gene of P. falciparum collected along the Thai–Cambodian border
Full size table

Table 6 shows the comparisons of in vitro anti-malarial sensitivities of P. falciparum isolates with different pfmdr1 and k13 genotypes. Only six groups (I–VI) were compared because varied genotypes with a small number were added in group VII. Significant differences were observed in the MQ and PPQ IC50s among these six groups (p = 0.002 and 0.013, respectively, One-way ANOVA). Post Hoc analysis using Scheffe test showed no significant difference between any pair of these genotypes which may have been due to a lower number of parasites in some subgroups.

Table 6 In vitro antimalarial sensitivities of different haplotype subgroups of P. falciparum
Full size table

Discussion

Due to the emergence of artemisinin resistance along the Thai–Cambodia border especially in Chanthaburi and Trat Provinces, the artemisinin resistance containment project was launched in 2009 by the Ministry of Public Health, Thailand [17]. Reduction of artemisinin pressure was one of objectives in this project by replacing AS–MQ with a fixed dose combination of ATQ–PG. In this study, both phenotypes and genotypes of P. falciparum isolates collected after the artemisinin resistance containment project, 2009–2016 were compared with the isolates collected before 2009. The parasites collected from 2009 to 2016 showed significantly higher CQ, QN, MQ, PPQ and AS IC50s compared with the parasites collected before 2009. The IC50s of these Thai isolates were in the same range as those reported in Cambodian isolates in 2013 [33]. Previously described cut-off points for in vitro anti-malarial resistance were used to determine the parasite’s resistant phenotypes [28,29,30,31,32]. Of 101, only one isolate collected before 2009 exhibited CQ sensitive. No CQ-sensitive isolate was detected among the parasites collected from 2009 to 2016. Although CQ was not used to treat falciparum malaria, Thai isolates of P. falciparum remain resistant to CQ. This may be due to vivax malaria sharing similar endemic areas with falciparum malaria. Thus, CQ, the first-line treatment for vivax malaria could cause a drug pressure for P. falciparum as well. The cure rate of MQ has rapidly declined soon after using as a monotherapy to treat falciparum malaria in 1991 [34]. Because MQ is a long half-life drug, drug pressure could cause the emergence of MQ resistance. As a result, AS–MQ combination was used as the first-line treatment of uncomplicated falciparum malaria since 1995. The parasites collected from 2009 to 2016 showed significantly increased MQ IC50 compared with the parasites before 2009. In addition, approximately 57% of the isolates collected after 2009 exhibited in vitro MQ resistance. Increased MQ resistance after 2009 may be influenced by the delayed parasite clearance phenotype of P. falciparum against AS in these areas. Slow parasite clearance causes more parasites to be exposed to the partner drug, i.e., MQ, increasing the risk of resistance selection of the partner drug, which in turn increases the risk of treatment failure. In the present study, the parasites with the pfmdr1 184F allele showed a significantly higher MQ IC50 than others.

A few cases of ATQ–PG treatment failure have been reported. Treatment failure of ATQ–PG was due to ATQ resistance and has been linked to point mutations in the target gene, the cytb gene [20, 32, 35]. Determination of the phenotypes and genotypes related to ATQ–PG response in Thai isolates of P. falciparum will be useful for rational drug use. According to Musset and colleagues (2006) [32], the cut-off point for in vitro ATQ resistance was the IC50 > 1900 nM. None of parasite isolates in this study exhibited ATQ resistance. In addition, they contained no mutations in the cytb 268 codon, molecular markers for ATQ resistance. The present results are similar to our survey in 2008 showing no evidence of ATQ resistance in Thai isolates of P. falciparum collected from both Thai–Cambodian and Thai–Myanmar borders [36]. A recent study in Cambodia also showed similar results indicating that P. falciparum isolated from Western Cambodia remained sensitive to ATQ in vitro and showed no point mutations in the cytb gene [33]. Recently, a successful cure of a multidrug-resistant falciparum case after artemisinin-based and QN-based treatment failure was reported in a subject that traveled to Cambodia [37]. These results suggest that a fixed-dose combination of ATQ–PG could be used the artemisinin-resistant areas with careful monitoring.

In this study, the IC50 of AS but not DHA was increased among the parasites isolated from 2009 to 2016. However, using the IC50 of > 10.5 nM as the cut-off point for in vitro AS and DHA resistance [31], no parasite exhibiting AS and DHA resistance was collected in the year before and after the artemisinin resistance containment project. Although no evidence exists of full artemisinin resistance, partial artemisinin resistance defined by delayed parasite clearance following treatment with an AS monotherapy or with an ACT is widespread in the Great Mekong Subregion [38]. To date, more than 200 nonsynonymous mutations in the k13 gene have been reported. Several mutations in the k13 gene were associated with delayed parasite clearance in vivo and in vitro including N458Y, Y493H, R539T, I543T, R561H and C580Y [38]. In the GMS, mutations in the k13 gene have spread and are distinctly reported according to their geographical areas [38, 39]. In the eastern GMS including Thai–Cambodia border, C580Y, R539T, Y493H, I543T, and P553L were commonly identified with the domination of C580Y. In the present study, 7 SNPs were identified including G436S, F483L, Y493H, G538 V, R539T, V568G, and C580Y. The most common SNPs of the isolates collected from 2009 to 2016 were C580Y (63.5%) and R539T (15.4%). Ring survival assay was not performed in the present study, however, reduced in vitro AS sensitivity was identified by in vitro sensitivity assay in the parasites with the k13 539T allele compared with others. For the newly identified SNPs including G436S and F483L, validation as a resistance marker will be required.

After 2009, ATQ–PG has been used to reduce artemisinin pressure in this area. In the absence of drug pressure, some resistant parasites might be less fit than their sensitive counterparts [40]. However, after the artemisinin resistance containment project, k13 mutations had increased significantly from 18.4 to 84.6%. Of these mutations, the k13 C580Y allele is increasing and replacing other haplotypes along the Thai-Cambodia border indicating a selective sweep in these areas. A study of parasites collected in 2007 found that 50% (11/22) of parasites from Chanthaburi and Trat Provinces contained the k13 mutations [41]. Of these, 45.5% (10/22) contained the k13 580Y allele indicating that parasites with the k13 580Y allele spread widely before the artemisinin resistance containment project. Recent studies indicated that the parasites with the k13 580Y allele arose in western Cambodia and then spread to other countries in the western Great Mekong Subregion including Thailand, Lao PDR and Vietnam [12, 41,42,43]. Thus, sensitive parasites might not compete with the main haplotype, the k13 580Y allele. One other factor that might influence the spread of artemisinin resistance in these areas is cross-board migration. Approximately one half of malaria cases in Thailand were foreign migrant workers [2]. As part of the artemisinin resistance containment project in Thailand, AS–MQ has been replaced by ATQ–PG to reduce artemisinin pressure in Chanthaburi and Trat provinces, Thailand which has been implemented since 2009. However, ACT remains the first-line treatment for uncomplicated falciparum malaria in Cambodia [38]. Thus, artemisinin pressure in these areas might not be effectively reduced. In addition, different policies and implementation of primaquine, as a P. falciparum gametocytocide, may influence the spreading of artemisinin resistance along Thai–Cambodian border [44, 45].

DHA-PPQ was used as the first-line drug for multidrug-resistant falciparum malaria in Cambodia [46]. Unfortunately, treatment failure of DHA-PPQ was promptly reported possibly due to the existing resistant parasites because PPQ monotherapy was used in Cambodia in the 1990s [47]. A few studies showed a link between the pfmdr1 copy number and PPQ sensitivity [48,49,50]. However, no significant difference of PPQ IC50 between the parasites with one and more than pfmdr1 copy number was found in the present study (18.8 ± 6.3 and 22.9 ± 13.3, p = 0.078, Independent t test) similar to the recent study using parasites collected from both Thai–Myanmar and Thai–Cambodian border areas [51]. Genetic markers for PPQ resistance including nonsynonymous SNP encoding a Glu415Gly mutation in a putative exonuclease (exo-E415G) and plasmepsin 2–3 amplification have been identified [52, 53]. Both in vitro PPQ sensitivity and ring survival assay were used to identify the association between PPQ resistance and these genes. Treatment failure of DHA-PPQ in Cambodia has been associated with parasites containing the k13 mutations and multiple plasmepsin 2 copy [53]. Recent studies have shown that parasites with the k13 580Y allele and plasmepsin 2 amplification have emerged and spread widely in the western Mekong Basin Subregion causing DHA-PPQ treatment failure [42, 43]. In the present study, the parasites containing the k13 580Y alleles exhibited significant higher PPQ IC50 compared with the others. This could be explained by parasites with the 580Y allele acquiring reduced PPQ susceptibility in this area. Unfortunately, other genetic markers for PPQ resistance including exo-E415G and plasmepsin 2–3 amplification were not determined in our study.

Conclusion

ATQ–PG, one of the non-ACT combinations, might be suitable to treat uncomplicated falciparum malaria in multidrug-resistant areas. Reduced ATQ sensitivity and mutation in the target gene has not been identified after the artemisinin resistance containment project. However, reduced artemisinin pressure using this combination might not be accomplished. Parasites with the k13 mutations, particularly the C580Y mutation, have increased over the years, even after the artemisinin resistance containment project. In addition, reduced sensitivity of the partner drugs of ACT including MQ and PPQ has been shown. New combinations which overcome these resistant phenotypes and genotypes should be carefully selected.

References

  1. World Health Organization. World malaria report 2012. Geneva: World Health Organization; 2012.

    Book  Google Scholar 

  2. World Health Organization. Global report on antimalarial drug efficacy and drug resistance 2000–2010. Geneva: World Health Organization; 2010.

    Google Scholar 

  3. Vijaykadga S, Rojanawatsirivej C, Cholpol S, Phoungmanee D, Nakavej A, Wongsrichanalai C. In vivo sensitivity monitoring of mefloquine monotherapy and artesunate-mefloquine combinations for the treatment of uncomplicated falciparum malaria in Thailand in 2003. Trop Med Int Health. 2006;11:211–9.

    Article  PubMed  CAS  Google Scholar 

  4. Denis MB, Tsuyuoka R, Poravuth Y, Narann TS, Seila S, Lim C, et al. Surveillance of the efficacy of artesunate and mefloquine combination for the treatment of uncomplicated falciparum malaria in Cambodia. Trop Med Int Health. 2006;11:1360–6.

    Article  PubMed  CAS  Google Scholar 

  5. Price RN, Uhlemann AC, Brockman A, McGready R, Ashley E, Phaipun L, et al. Mefloquine resistance in Plasmodium falciparum and increased pfmdr1 gene copy number. Lancet. 2004;364:438–47.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  6. Rogers WO, Sem R, Tero T, Chim P, Lim P, Muth S, et al. Failure of artesunate-mefloquine combination therapy for uncomplicated Plasmodium falciparum malaria in southern Cambodia. Malar J. 2009;8:10.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  7. Noedl H, Se Y, Schaecher K, Smith BL, Socheat D, Fukuda MM, et al. Evidence of artemisinin-resistant malaria in western Cambodia. N Engl J Med. 2008;359:2619–20.

    Article  PubMed  CAS  Google Scholar 

  8. Dondorp AM, Nosten F, Yi P, Das D, Phyo AP, Tarning J, et al. Artemisinin resistance in Plasmodium falciparum malaria. N Engl J Med. 2009;361:455–67.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  9. Phyo AP, Nkhoma S, Stepniewska K, Ashley EA, Nair S, McGready R, et al. Emergence of artemisinin-resistant malaria on the western border of Thailand: a longitudinal study. Lancet. 2012;379:1960–6.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Ashley EA, Dhorda M, Fairhurst RM, Amaratunga C, Lim P, Suon S, et al. Tracking resistance to artemisinin Collaboration (TRAC): spread of artemisinin resistance in Plasmodium falciparum malaria. N Engl J Med. 2014;371:411–23.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  11. Witkowski B, Amaratunga C, Khim N, Sreng S, Chim P, Kim S, et al. Novel phenotypic assays for the detection of artemisinin-resistant Plasmodium falciparum malaria in Cambodia: in vitro and ex vivo drug-response studies. Lancet Infect Dis. 2013;13:1043–9.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  12. Ariey F, Witkowski B, Amaratunga C, Beghain J, Langlois AC, Khim N, et al. A molecular marker of artemisinin-resistant Plasmodium falciparum malaria. Nature. 2014;505:50–5.

    Article  PubMed  CAS  Google Scholar 

  13. Straimer J, Gnädig NF, Witkowski B, Amaratunga C, Duru V, Ramadani AP, et al. Drug resistance. K13-propeller mutations confer artemisinin resistance in Plasmodium falciparum clinical isolates. Science. 2015;347:428–31.

    Article  PubMed  CAS  Google Scholar 

  14. Takala-Harrison S, Jacob CG, Arze C, Cummings MP, Silva JC, Dondorp AM, et al. Independent emergence of artemisinin resistance mutations among Plasmodium falciparum in Southeast Asia. J Infect Dis. 2015;211:670–9.

    Article  PubMed  CAS  Google Scholar 

  15. Tun KM, Imwong M, Lwin KM, Win AA, Hlaing TM, Hlaing T, et al. Spread of artemisinin-resistant Plasmodium falciparum in Myanmar: a cross-sectional survey of the K13 molecular marker. Lancet Infect Dis. 2015;15:415–21.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  16. Ye R, Hu D, Zhang Y, Huang Y, Sun X, Wang J, et al. Distinctive origin of artemisinin-resistant Plasmodium falciparum on the China–Myanmar border. Sci Rep. 2016;6:20100.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  17. Satimai W, Sudathip P, Vijaykadga S, Khamsiriwatchara A, Sawang S, Potithavoranan T, et al. Artemisinin resistance containment project in Thailand. II: responses to mefloquine–artesunate combination therapy among falciparum malaria patients in provinces bordering Cambodia. Malar J. 2012;11:300.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  18. Srivastava IK, Rottenberg H, Vaidya AB. Atovaquone, a broad spectrum antiparasitic drug, collapses mitochondrial membrane potential in a malarial parasite. J Biol Chem. 1997;272:3961–6.

    Article  PubMed  CAS  Google Scholar 

  19. Kessl JJ, Ha KH, Merritt AK, Lange BB, Hill P, Meunier B, et al. Cytochrome b mutations that modify the ubiquinol-binding pocket of the cytochrome bc1 complex and confer anti-malarial drug resistance in Saccharomyces cerevisiae. J Biol Chem. 2005;280:17142–8.

    Article  PubMed  CAS  Google Scholar 

  20. Schwöbel B, Alifrangis M, Salanti A, Jelinek T. Different mutation patterns of atovaquone resistance to Plasmodium falciparum in vitro and in vivo: rapid detection of codon 268 polymorphisms in the cytochrome b as potential in vivo resistance marker. Malar J. 2003;2:5.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Korsinczky M, Chen N, Kotecka B, Saul A, Rieckmann K, Cheng Q. Mutations in Plasmodium falciparum cytochrome b that are associated with atovaquone resistance are located at a putative drug-binding site. Antimicrob Agents Chemother. 2000;44:2100–8.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  22. Trager W, Jensen JB. Human malaria parasites in continuous culture. Science. 1976;193:673–5.

    Article  PubMed  CAS  Google Scholar 

  23. Desjardins RE, Canfield CJ, Haynes JD, Chulay JD. Quantitative assessment of antimalarial activity in vitro by a semiautomated microdilution technique. Antimicrob Agents Chemother. 1979;16:710–8.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  24. Wooden J, Gould EE, Paull AT, Sibley CH. Plasmodium falciparum: a simple polymerase chain reaction method for differentiating strains. Exp Parasitol. 1992;75:207–12.

    Article  PubMed  CAS  Google Scholar 

  25. Djimdé A, Doumbo OK, Cortese JF, Kayentao K, Doumbo S, Diourté Y, et al. A molecular marker for chloroquine-resistant falciparum malaria. N Engl J Med. 2001;344:257–63.

    Article  PubMed  Google Scholar 

  26. Duraisingh MT, Roper C, Walliker D, Warhurst DC. Increased sensitivity to the antimalarials mefloquine and artemisinin is conferred by mutations in the pfmdr1 gene of Plasmodium falciparum. Mol Microbiol. 2000;36:955–61.

    Article  PubMed  CAS  Google Scholar 

  27. Price RN, Uhlemann AC, Brockman A, McGready R, Ashley E, Phaipun L, et al. Mefloquine resistance in Plasmodium falciparum and increased pfmdr1 gene copy number. Lancet. 2004;364:438–47.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  28. Bickii J, Basco LK, Ringwald P. Assessment of three in vitro tests and an in vivo test for chloroquine resistance in Plasmodium falciparum clinical isolates. J Clin Microbiol. 1998;36:243–7.

    PubMed  PubMed Central  CAS  Google Scholar 

  29. Basco LK, Le Bras J. In vitro susceptibility of Cambodian isolates of Plasmodium falciparum to halofantrine, pyronaridine and artemisinin derivatives. Ann Trop Med Parasitol. 1994;88:137–44.

    Article  PubMed  CAS  Google Scholar 

  30. Hatin I, Trape JF, Legros F, Bauchet J, Le Bras J. Susceptibility of Plasmodium falciparum strains to mefloquine in an urban area in Senegal. Bull World Health Organ. 1992;70:363–7.

    PubMed  PubMed Central  CAS  Google Scholar 

  31. Pradines B, Bertaux L, Pomares C, Delaunay P, Marty P. Reduced in vitro susceptibility to artemisinin derivatives associated with multi-resistance in a traveller returning from South-East Asia. Malar J. 2011;10:268.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  32. Musset L, Bouchaud O, Matheron S, Massias L, Le Bras J. Clinical atovaquone-proguanil resistance of Plasmodium falciparum associated with cytochrome b codon 268 mutations. Microbes Infect. 2006;8:2599–604.

    Article  PubMed  CAS  Google Scholar 

  33. Saunders DL, Chaorattanakawee S, Gosi P, Lanteri C, Somethy S, Kuntawunginn W, et al. Atovaquone-proguanil remains a potential stopgap therapy for multidrug-resistant Plasmodium falciparum in areas along the Thai–Cambodian border. Antimicrob Agents Chemother. 2015;60:1896–8.

    Article  PubMed  CAS  Google Scholar 

  34. Nosten F, ter Kuile F, Chongsuphajaisiddhi T, Luxemburger C, Webster HK, Edstein M, et al. Mefloquine-resistant falciparum malaria on the Thai–Burmese border. Lancet. 1991;337:1140–3.

    Article  PubMed  CAS  Google Scholar 

  35. Fivelman QL, Butcher GA, Adagu IS, Warhurst DC, Pasvol G. Malarone treatment failure and in vitro confirmation of resistance of Plasmodium falciparum isolate from Lagos, Nigeria. Malar J. 2002;1:1.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Khositnithikul R, Tan-Ariya P, Mungthin M. In vitro atovaquone/proguanil susceptibility and characterization of the cytochrome b gene of Plasmodium falciparum from different endemic regions of Thailand. Malar J. 2008;7:23.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  37. Dell’Acqua R, Fabrizio C, Di Gennaro F, Lo Caputo S, Saracino A, Menegon M, et al. An intricate case of multidrug resistant Plasmodium falciparum isolate imported from Cambodia. Malar J. 2017;16:149.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  38. World Health Organization. Artemisinin and artemisinin-based combination therapy resistance. Geneva: World Health Organization; 2017.

    Google Scholar 

  39. Ménard D, Khim N, Beghain J, Adegnika AA, Shafiul-Alam M, Amodu O, et al. A worldwide map of Plasmodium falciparum K13-propeller polymorphisms. N Engl J Med. 2016;374:2453–64.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  40. Anderson TJ, Roper C. The origins and spread of antimalarial drug resistance: lessons for policy makers. Acta Trop. 2005;94:269–80.

    Article  PubMed  Google Scholar 

  41. Talundzic E, Okoth SA, Congpuong K, Plucinski MM, Morton L, Goldman IF, et al. Selection and spread of artemisinin-resistant alleles in Thailand prior to the global artemisinin resistance containment campaign. PLoS Pathog. 2015;11:e1004789.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  42. Imwong M, Suwannasin K, Kunasol C, Sutawong K, Mayxay M, Rekol H, et al. The spread of artemisinin-resistant Plasmodium falciparum in the Greater Mekong subregion: a molecular epidemiology observational study. Lancet Infect Dis. 2017;17:491–7.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Imwong M, Hien TT, Thuy-Nhien NT, Dondorp AM, White NJ. Spread of a single multidrug resistant malaria parasite lineage (PfPailin) to Vietnam. Lancet Infect Dis. 2017;17:1022–3.

    Article  PubMed  Google Scholar 

  44. World Health Organization. Single dose primaquine as a gametocyte in Plasmodium falciparum malaria. Updated WHO Policy Recommendation (October 2012). Geneva: Global Malaria Programme, World Health Organization; 2012.

    Google Scholar 

  45. Kitchakarn S, Lek D, Thol S, Hok C, Saejeng A, Huy R, et al. Implementation of G6PD testing and primaquine for P. vivax radical cure: operational perspectives from Thailand and Cambodia. WHO South East Asia. J Public Health. 2017;6:60–8.

    Google Scholar 

  46. Leang R, Barrette A, Bouth DM, Menard D, Abdur R, Duong S, et al. Efficacy of dihydroartemisinin-piperaquine for treatment of uncomplicated Plasmodium falciparum and Plasmodium vivax in Cambodia, 2008–2010. Antimicrob Agents Chemother. 2013;57:818–26.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  47. Davis TM, Hung TY, Sim IK, Karunajeewa HA, Ilett KF. Piperaquine: a resurgent antimalarial drug. Drugs. 2005;65:75–87.

    Article  PubMed  CAS  Google Scholar 

  48. Eastman RT, Dharia NV, Winzeler EA, Fidock DA. Piperaquine resistance is associated with a copy number variation on chromosome 5 in drug-pressured Plasmodium falciparum parasites. Antimicrob Agents Chemother. 2011;55:3908–16.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  49. Veiga MI, Ferreira PE, Malmberg M, Jörnhagen L, Björkman A, Nosten F, et al. Pfmdr1 amplification is related to increased Plasmodium falciparum in vitro sensitivity to the bisquinoline piperaquine. Antimicrob Agents Chemother. 2012;56:3615–9.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  50. Veiga MI, Dhingra SK, Henrich PP, Straimer J, Gnädig N, Uhlemann AC, et al. Globally prevalent PfMDR1 mutations modulate Plasmodium falciparum susceptibility to artemisinin-based combination therapies. Nat Commun. 2016;7:11553.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  51. Mungthin M, Watanatanasup E, Sitthichot N, Suwandittakul N, Khositnithikul R, Ward SA. Influence of the pfmdr1 gene on in vitro sensitivities of piperaquine in Thai isolates of Plasmodium falciparum. Am J Trop Med Hyg. 2017;96:624–9.

    PubMed  PubMed Central  CAS  Google Scholar 

  52. Amato R, Lim P, Miotto O, Amaratunga C, Dek D, Pearson RD, et al. Genetic markers associated with dihydroartemisinin-piperaquine failure in Plasmodium falciparum malaria in Cambodia: a genotype-phenotype association study. Lancet Infect Dis. 2017;17:164–73.

    Article  PubMed  CAS  Google Scholar 

  53. Witkowski B, Duru V, Khim N, Ross LS, Saintpierre B, Beghain J, et al. A surrogate marker of piperaquine-resistant Plasmodium falciparum malaria: a phenotype-genotype association study. Lancet Infect Dis. 2017;17:174–83.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

Download references

Authors’ contributions

TT, PT, and MM conceived of the study, participated in the design and coordination of the study and performed the statistical analysis. TT, NSu, TR and KM performed molecular analysis. TT, PJ, JT and NSi carried out the in vitro cultivation and sensitivity test. All authors read and approved the final manuscript.

Acknowledgements

This study was financially supported by the Health System Research Institute/National Science and Technology Development Agency (P-13-50112) and the Phramongkutklao Research Fund.

Competing interests

The authors declare that they have no competing interests.

Availability of data and materials

The data and detailed protocol can be made available upon request from the corresponding author.

Consent for publication

Not applicable.

Ethics approval and consent to participate

The research protocol was reviewed and approved by the Ethics Committee of the Royal Thai Army Medical Department.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information

Authors and Affiliations

  1. Department of Microbiology, Faculty of Science, Mahidol University, Bangkok, Thailand

    Thunyapit Thita, Pimrat Jadsri, Jarupatr Thamkhantho, Kittiya Mahotorn & Peerapan Tan-ariya

  2. Department of Parasitology, Phramongkutklao College of Medicine, Bangkok, Thailand

    Toon Ruang-areerate, Nantana Suwandittakul, Naruemon Sitthichot & Mathirut Mungthin

Authors
  1. Thunyapit Thita
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Pimrat Jadsri
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Jarupatr Thamkhantho
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Toon Ruang-areerate
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Nantana Suwandittakul
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Naruemon Sitthichot
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Kittiya Mahotorn
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Peerapan Tan-ariya
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Mathirut Mungthin
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding author

Correspondence to Mathirut Mungthin.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Thita, T., Jadsri, P., Thamkhantho, J. et al. Phenotypic and genotypic characterization of Thai isolates of Plasmodium falciparum after an artemisinin resistance containment project. Malar J 17, 197 (2018). https://doi.org/10.1186/s12936-018-2347-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12936-018-2347-9

Keywords

Malaria Journal

ISSN: 1475-2875

Contact us