Origins and spread of novel genetic variants of sulfadoxine–pyrimethamine resistance in Plasmodium falciparum isolates in Indonesia
Malaria Journalvolume 17, Article number: 475 (2018)
While malaria incidence in Indonesia has decreased threefold in the last decade, more than 200,000 cases were reported in 2016. Different endemicity of Plasmodium falciparum malaria among several islands in Indonesia has been recognized and two unique mutations of P. falciparum dihydropteroate synthase (pfdhps) affecting sulfadoxine–pyrimethamine (SP) resistance were detected from the research of SP efficiency and genotype analysis in South Kalimantan. In this study, geographical distribution and origin of these pfdhps K540T and I588F mutations were analysed.
Malaria parasites DNA from several endemic areas in Indonesia; Sumatera, Java, Kalimantan, Lombok, Sumbawa, Timor, Sulawesi, and Papua islands; in two periods, 2004–2006 and 2009–2012 were subjected for pfdhfr and pfdhps sequence analysis.
Different genotype polymorphisms of pfdhfr and pfdhps were observed in the parasites from various regions in Indonesia and relatively more divergent genotypes were determined from Kalimantan isolates in both 2004–2006 and 2009–2012. The parasites containing K540T mutation were identified in 2004–2006 isolates from East Kalimantan, East Java and Sumbawa as an SGTGA haplotype. The other I588F mutation was also determined in 2004–2006 parasites, isolated from Lombok and Sumbawa islands as an SGEAA(588F) haplotype. The parasites with pfdhfr/pfdhps quintuple or sextuple mutation, a genotype marker of SP resistance, were determined mostly in Kalimantan in both 2004–2006 and 2009–2012.
Analysis of the prevalence and pfdhfr/pfdhps combined genotypes of K540T or I588F mutations suggested that K540T might be origin in Kalimantan Island and I588F in Sumbawa Island and then these were spread to other areas along with people movement. This research indicates regular monitoring of drug efficacy and parasite genotype analysis is important to keep efficiency and prevent the spread of resistance. It is also essential for the latest anti-malarial drug artemisinin-based combination therapy.
Malaria incidence in Indonesia has decreased threefold in the last decade , and there is no evidence for the presence of the parasites resistant to artemisinin-based combination therapy (ACT) . However, about 25% of Indonesia’s population (total population: 261 million) is at risk of malaria, and more than 200,000 positive cases were reported in 2016 . Detection of drug resistant malaria parasites and prevention of these spreading are critical to keep efficacy of the malaria treatment and to obtain elimination of malaria. In this research, origin and distribution of unique mutations of Plasmodium falciparum sulfadoxine–pyrimethamine (SP) resistance in Indonesia were analysed.
SP has been widely used as an anti-malarial drug for treatment of uncomplicated malaria, and for intermittent preventive treatment in vulnerable populations, pregnant women (IPTp) and infants (IPTi) in high malaria transmission areas in Africa [4, 5]. However, emergence and spread of SP resistant P. falciparum has been reported worldwide [6,7,8,9]. In Indonesia, SP was recommended as a second line anti-malarial drug after chloroquine resistance had been determined in 1973, and P. falciparum resistance to SP was reported for the first time in Jayapura (Papua Province) in 1979 [10, 11]. Chloroquine and SP had been used in Indonesia until 2008, when the malaria treatment policy was changed. ACT, using a combination of an artemisinin derivative with another anti-malarial, such as piperaquine, lumefantrine or amodiaquine, is provided as the first line anti-malarial drug for treatment of uncomplicated malaria, and SP is not administered for malaria treatment. However, people, especially in local areas, use SP when they are suffering from malaria, or for chemoprophylaxis , because of some effectiveness, low cost, simple administration as a single oral dose, and fewer side effects. It is important to obtain information about SP resistance in malaria endemic areas in Indonesia.
The mutations in P. falciparum dihydrofolate reductase (pfdhfr) and dihydropteroate synthase (pfdhps) are responsible for pyrimethamine and sulfadoxine resistance, respectively [13,14,15,16]. Stepwise accumulation of point mutations in pfdhfr and pfdhps genes is associated with higher level of resistance to SP in vitro and in vivo [17, 18]. The amino acid substitution at position 108 serine to asparagine or threonine (S108N/T) in pfdhfr is essential for subsequent A16V, N51I, C59R and I164L mutations (underlined bold type indicates the mutant allele), leading to high-level of resistance to cycloguanil or pyrimethamine . Similarly, a single mutation in the pfdhps converting alanine to glycine at amino acid position 437 (A437G), which is normally the first mutation under the sulfadoxine drug pressure, conferred on the parasite a fivefold higher level of drug resistance . Additional mutations K540E and A581G, then S436A/F and A613S/T are associated with increasing resistance to sulfadoxine. A combination of pfdhfr triple (N51I, C59R and S108 N) and pfdhps double (A437G and K540E) mutations collectively form the quintuple mutations, which is strongly associated with in vitro and in vivo SP resistance [17, 20,21,22,23,24].
There are several reports of mutation analysis of pfdhfr and pfdhps genes from different malaria endemic areas in Indonesia. Different ratio of mutation level was reported from sample analysis of several island and district parasites. From the West Papua sample analysis obtained in 1996 to 1999 by Nagesha et al. , C59R and S108N mutations in pfdhfr and A437G in pfdhps were commonly determined in SP resistant parasites. In addition, they reported an additional K540E mutation in pfdhps was observed in more resistant level parasites. Extensive analysis of the parasite genotypes from eight malaria endemic areas were reported by Syafruddin et al. in 2005, representing a broad region of the western and eastern Indonesian archipelago . Polymorphisms in pfdhfr gene at S108N/T, A16V and C59R were frequently identified, in which A16V were observed in association with S108T and these were differently distributed, more common among samples from eastern regions. The polymorphism in pfdhps was less frequent in this report; about 15% of A437G and less than 5% of K540E were detected. Among the Sumba island samples of 2007, less frequent mutations were reported by Asih et al. in 2009 ; about 25% parasites presented S108N and C59R mutations in pfdhfr and only few % of the isolates presented A437G mutant allele of pfdhps. More prevalence of mutations in pfdhfr and pfdhps was observed from a study of SP efficacy and genotype analysis in South Kalimantan in 2009–2010 . More than 90% of the isolates exhibited S108N and C59R mutations in pfdhfr gene, in addition I164L substitution was detected in 30% of the parasites. The alterations in pfdhps were detected at the amino acid positions of A437, K540, A581 and I588 to glycine (97%), glutamine or threonine (36%, 36%), glycine (45%) and phenylalanine (23%), respectively. This result of pfdhfr and pfdhps genotypes is quite unique compared with the previous Indonesia parasites; more prevalence of mutation ratio and detection of novel mutations of pfdhps K540T and I588F.
In this study, pfdhfr and pfdhps sequences of P. falciparum from different malaria endemic regions, mostly eastern part of Indonesia were analysed to investigate polymorphisms of pfdhfr and pfdhps genotypes and predict susceptibility for sulfadoxine–pyrimethamine in malaria parasites in Indonesia, and to obtain information about distribution of the previously identified pfdhps K540T and I588F novel mutations.
Study sites and malaria patients
Malaria parasites collected from the following several endemic areas in Indonesia were analysed in this study: Sumatera (Indragiri Hilir, Merangin), Kalimantan (Paser, Seruyan, Banjar), Java (Pacitan), Lombok (West Lombok), Sumbawa (Sumbawa), Timor (Timor Tengah Selatan), Sulawesi (Gorontalo), and Papua (Jayapura) islands in 2004–2006 and in 2009–2012. Malaria patients were recruited at primary health centers, district hospitals or local field areas and written informed consent were obtained from each participant, or from caretakers if participations were under 12 years old, after explanation of the purpose of the study in local language. Finger prick blood samples were collected on a glass slide for microscopical observation and on a filter paper (Advantec, Toyo Roshi Kaisha, Ltd., Japan) for parasite DNA analysis. Thick smear blood films were stained with Giemsa (Merck, Germany) and examined microscopically for the presence of malaria parasites. Dried blood spot on a filter paper was kept in a small plastic clips prior to parasite DNA extraction. In some field studies, BinaxNOW Malaria diagnosis kit (Binax Inc, Portland, ME, USA) were used to identify P. falciparum malaria patients. Patients with positive malaria diagnosis results were treated with an anti-malarial drug according to national policy. The study protocol was reviewed and approved by the Ethical Committee, Faculty of Medicine, Universitas Airlangga, Surabaya, Indonesia and Institute of Tropical Medicine, Nagasaki University, Nagasaki, Japan.
Parasite DNA analysis
Parasite DNA was extracted from the dried blood spots on filter paper by using QIAamp DNA blood mini kit (Qiagen, Netherlands) and kept at − 30 °C. Plasmodium falciparum samples were selected by nested PCR methodology using species specific primer sets of 18S rRNA genes described in Snounou et al.  and Kimura et al. . Pfdhfr and pfdhps genotypes were determined by sequencing as previously described by Isozumi et al. . using amplified PCR product as a template directly for sequence analysis. Alleles corresponding to amino acid positions at 16, 51, 59, 108 and 164 of the pfdhfr gene and at 436, 437, 540, 581, 588 and 613 of the pfdhps gene were read more carefully, and at least two independent PCR products were prepared for sequence analysis in the case of rare mutations. Previously reported results from South Kalimantan Province  were included in this analysis.
Data were entered in Microsoft Excel and exported to SPSS version 17.0 for analysis. Chi square and Fisher’s exact tests were used, where applicable, to assess the relationship of mutations between two periods of studies. Allele proportions were calculated as the number carrying a certain allele divided by the number of samples with positive PCR outcome.
A total of 622 P. falciparum samples from symptomatic to asymptomatic malaria patients during the two periods, 2004–2006 and 2009–2012, were analysed in this study. The research areas were different levels of malaria endemicity, and presented seasonal differences in number of malaria patients. Same districts in Lombok and Papua islands, the latter is reported as high malaria endemicity, were included for sample analysis in both periods. Of the 384 samples from 622 patients, P. falciparum pfdhfr and pfdhps genes were successfully amplified and analysed both genotypes (Tables 1, 2, 3). The genotypes were obtained from more than 20 samples of Kalimantan, Lombok and Papua isolates in both periods of 2004–2006 and 2009–2012, except Papua in 2004–2006 which was 9 samples.
Pfdhfr and pfdhps genotypes
The pfdhfr and pfdhps genotypes were determined from 384 samples, 205 in 2004–2006 and 179 in 2009–2012. Polymorphisms were identified at amino acid positions of C59, S108 and I164 in pfdhfr to arginine, aspargine and leucine, respectively, but no other substitutions including at A16 and N51 were observed. In pfdhps gene, substitutions were observed at amino acid positions of A437, K540, A581 and I588 to glycine, glutamic acid or threonine, glycine and phenylalanine, respectively. No mutations were detected at the positions of S436 and A613. Comparison of the mutation ratios in 2004–2006 and 2009–2012 samples from Kalimantan, Lombok and Papua is presented in Fig. 1. Similar ratio of mutation was observed in 2004–2006 and 2009–2012 at each amino acid position of pfdhfr in Kalimantan, Lombok and Papua samples, but more mutation ratio was observed at I164L in 2009–2012 than 2004–2006 in Kalimantan. Some different percentages of the mutations were detected in pfdhps gene from Kalimantan, Lombok and Papua samples, and between 2004–2006 and 2009–2012 in the same island parasite samples. Remarkable change was noticed in the pfdhps I588F mutation; new appearance in Kalimantan and Papua samples in 2009–2012 isolates, but disappearance from Lombok samples in 2009–2012.
With the pfdhfr mutations at C59R, S108N and I164L, in total, the wild type pfdhfr allele ANCSI (at positions 16, 51, 59, 108 and 164, respectively) and three mutant alleles, ANCNI, ANRNI and ANRNL were detected in this study (Table 1).
Relatively high ratio of mutation in pfdhps gene was determined at A437 to glycine (Fig. 1), which is regarded as the first and most essential substitution for sulfadoxine resistance . Two types of mutations at amino acid position K540 were detected as AAA (lysine codon) to GAA (glutamic acid) or ACA (threonine). One of this K540T and another unique mutation I588F are recently reported alleles in Indonesia, recognized through previous South Kalimantan sample analysis . The K540T was also reported in Sabah, Malaysia , which is located in the same island. Altogether, the wild type pfdhps haplotype SAKAA (at positions 436, 437, 540, 581 and 613, respectively) and six different mutant alleles SGKAA, SGKGA, SGEAA, SGEAA(588F), SGTGA and SGEGA were identified in this study (Table 2). The unique I588F mutation was detected only in SGEAA(588F) haplotype and K540T was only in SGTGA allele in both 2004–2006 and 2009–2012 samples.
Diverse polymorphisms in different malaria endemic areas
Among the observed four pfdhfr genotypes with simple stepwise accumulation pattern of mutations, ANRNI is the most common pfdhfr allele in most of the research locations both in 2004–2006 and 2009–2012 periods except Seruyan (Middle Kalimantan) in 2009–2012, where the ANRNL is dominant (Fig. 2a). The ANRNL allele is frequently detected in Paser (East Kalimantan) and Pacitan (East Jawa) in 2004–2006 samples, and Banjar (South Kalimantan) in 2009–2012. A remarkable observation is the wild type ANCSI allele found in Jayapura (Papua) samples in 2009–2012 had not been observed in 2004–2006.
Complicated genotype variations were observed in pfdhps allele haplotypes. Similar variant haplotypes were detected from Paser (East Kalimantan) and Pacitan (Easr Java) in 2004–2006 samples, but different haplotype variations were observed in parasite populations from most of the research areas (Fig. 2b). The substitution of K540E was detected as one of the SGEAA, SGEAA(588F) or SGEGA alleles at most of the malaria research areas in 2004–2006 and 2009–2012, except Indragiri Hilir (Riau, Sumatera Island), and Timor Tengah Selatan (Timor Island), where only small number of samples were obtained (Table 2). The unique I588F substitution was detected from the parasites with K540E mutation as the pfdhps SGEAA(588F) allele only in Lombok and Sumbawa Island samples in 2004–2006. It would have been then spreading to the other islands of eastern Indonesia and detected from Banjar in Kalimantan, Gorontalo in Sulawesi and Jayapura in Papua Islands in 2009–2012 (Fig. 2b). This mutant was not detected in Lombok in 2009–2012.
Another unique mutation of K540T was identified as a major pfdhps SGTGA allele in Paser (East Kalimantan) and Pacitan (East Java), and in Sumbawa as a minority in 2004–2006 (Fig. 2b). This was then detected only in Kalimantan Island in 2009–2012, as a major allele of Banjar (South Kalimantan) parasites, and one case of Seruyan (Middle Kalimantan) isolates (Fig. 2b).
Combined genotypes of pfdhfr/pfdhps haplotypes
With a combination of four pfdhfr and seven pfdhps haplotypes, totally 16 different combined pfdhfr/pfdhps genotypes were determined; 13 combined genotypes in 2004–2006 parasites and 16 in 2009–2012 samples (Table 3). Relatively more polymorphisms in the pfdhfr/pfdhps combined genotype were observed in Paser (East Kalimantan), Pacitan (East Java), Sumbawa samples in 2004–2006, and Seruyan (Middle Kalimantan), Banjar (South Kalimantan) and Jayapura (Papua) samples in 2009–2012. Only one ANRNI/SGEAA genotype was observed in 2004–2006 Jayapura samples, and several different combined genotypes were detected in 2009–2012 samples. The originally detected ANRNI/SGEAA remained a major genotype, but introduction of the parasites with the wild type pfdhfr and pfdhps between two research periods had increased number of types in the combined genotype.
Accumulation of mutations in pfdhfr and pfdhps genes enhances parasite resistance level against sulfadoxine–pyrimethamine. The parasites containing more than five mutations in combined pfdhfr/pfdhps genotype have been shown not to respond adequately to SP treatment [8, 22,23,24]. Prevalence of the quadruple ANRNI/SGEAA(588F) genotype, and quintuple or sextuple mutant genotypes in each research area is presented in Fig. 3. The parasites from Paser (East Kalimantan) and Pacitan (East Java) in 2004–2006 and from Seruyan (Middle Kalimantan) and Banjar (South Kalimantan) in 2009–2012 mostly presented SP resistant genotypes, especially the Pacitan parasite samples in 2004–2006 and Seruyan parasites in 2009–2012 were all belonging to the quintuple or sextuple mutants. The parasites from the other areas presented the genotypes containing four or less mutations in combination genotype that tend to respond adequately to SP treatment. Prevalence of parasite ratios harboring five or more mutations between 2004–2006 and 2009–2012 are not significantly different between two periods of time in Kalimantan, Lombok and Papua statistically (Table 4).
In this report, polymorphisms of pfdhfr and pfdhps genes in several malaria endemic areas in Indonesia are analysed in two periods of time, 2004–2006 and 2009–2012. Previously reported amino acid substitutions, important for sulfadoxine–pyrimethamine (SP) resistance in Indonesian parasites, were determined in this analysis. The C59R, S108N and I164L mutations in pfdhfr gene and A437G, K540E/T, A581G and I588F in pfdhps were observed from both 2004–2006 and 2009–2012 samples, but pfdhfr A16V and S108T mutations for cycloguanil resistance were not detected (Fig. 1). In pfdhfr gene, wild type ANCSI and three mutant ANCNI, ANRNI and ANRNL alleles were observed (Table 1). This simple accumulation pattern of mutations in the pfdhfr genotype supports the stepwise selection hypothesis of resistant gene evolution . Wild type SAKAA and six different mutant alleles were determined in pfdhps gene (Table 2). Comparison of wild type and each mutant alleles suggests stepwise accumulation of mutations in the pfdhps genotypes.
The unique K540T and I588F mutations of pfdhps, both of which were detected previously from 2009 to 2010 South Kalimantan sample analysis  and the former was also reported by Lau et al.  from 2010 Sabah, Malaysia parasites, were identified in 2004–2006 parasites. The evidence presenting here that K540T was detected in 2004–2006 as an SGTGA haplotype from several parasites in East Kalimantan, East Java and Sumbawa Island suggests a possibility of Kalimantan Island origin for this mutation. In East Kalimantan, the SGTGA was found from many patients as a combination genotype of either ANRNI/SGTGA or ANRNL/SGTGA, additionally as a mixed infection of both genotypes (Table 3). Meanwhile, the SGTGA existed as the combination genotype of ANRNI/SGTGA only in East Java and Sumbawa parasites.
Another novel mutation I588F was also detected in 2004–2006 samples from Lombok and Sumbawa islands; the parasites of SGEAA both with and without I588F mutation were detected in Sumbawa, whereas SGEAA(588F) and wild type SAKAA were observed in Lombok. It suggests a possibility that the initial mutation of I588F had occurred in SGEAA type parasites in Sumbawa Island, then introduced into its neighbor Lombok Island. Further analysis and comparison of microsatellite loci around the pfdhps haplotypes will provide additional information for understanding the origin and spreading of these K540T and I588F mutant alleles.
In Riau and Jambi Provinces (Sumatera Island), only small number of malaria patients had been detected and these were mostly infected with P. vivax and only a few P. falciparum were observed in both 2004–2006 and 2009–2012 periods. Timor Tengah Selatan (Timor Island) and Gorontalo (Sulawesi Island) were similar situation in 2009–2012. Remarkably pfdhps I588F mutation was observed in Gorontalo as the same combined genotype of ANRNI/SGEAA(588F).
Genotype comparison of 2004–2006 and 2009–2012 parasites in each malaria areas provides additional information (Table 3). In Lombok, parasites of the major combined genotype ANRNI/SAKAA and a few number of minor type were detected in both 2004–2006 and 2009–2012 periods. The parasites might have well responded to SP or other malaria therapies and new SP resistant type parasites had not settled in Lombok. On the other hand, increased number of pfdhfr/pfdhps combined genotype was determined in 2009–2012 parasites from Papua Island compared to 2004–2006 parasites. The genotypes were ANRNI/SGEAA(588F) and new combinations with either pfdhfr or pfdhps wild type. This suggests some parasites with wild type pfdhfr ANCSI and pfdhps SAKAA were introduced as the both wild type combination or either haplotype individually between two periods. The ANRNI/SGEAA had been selected under SP pressure before 2004–2006, and the new parasites with pfdhfr and/or pfdhps wild type alleles and ANRNI/SGEAA(588F) were maintained under reduced SP pressure after DHP was commonly used. ACT was introduced in Papua Province in 2005; first artesunate-amodiaquine (AA) was used experimentally, then changed to dihydroartemisinin–piperaquine (DHP) in 2008 . Wild type pfdhfr and pfdhps genotype detection was reported in the parasites from West Papua and several islands by Syafruddin et al. in 2005 . It suggests the parasites with wild type allele were introduced from a neighbour districts or it might exist in Papua in 2004–2006, however it could not detected because the number of samples was not sufficient.
Parasite populations from different Provinces in Kalimantan Island presented diverse polymorphisms. Analysis of Paser (East Kalimantan) 2004–2006 parasites revealed dominance of the pfdhps K540T mutation (SGTGA haplotype, Fig. 2b) and many cases of mixed genotype infections (Table 3). This suggests high P. falciparum infection rate under strong pressure of SP. The mixed genotype infection was not common in 2009–2012 samples from Middle and South Kalimantan, and the parasites from these Provinces presented characteristic genotype polymorphisms. Most of the parasites from Seruyan (Middle Kalimantan) acquired pfdhfr I164L mutation (ANRNL haplotype, Fig. 2a), and several pfdhps genotypes involving unique pfdhps K540T or I588F mutations were detected in Banjar (South Kalimantan) parasites (Fig. 2b).
High heterogeneity of malaria epidemiology and divergent genetic polymorphisms across islands, districts, and even close neighbour sub-districts in one island are not uncommon in Indonesia . Divergent polymorphisms in pfdhfr/pfdhps genotype of the parasite populations from different Kalimantan provinces were observed in this study from analysis of the parasite samples before introduction of ACT at each research site in these Provinces. AA was implemented in 2006 in East Kalimantan and then DHP treatment was started in 2009. DHP has been first-line treatment therapy against malaria in Middle Kalimantan and South Kalimantan since 2010. However in fact, ACT was introduced gradually form one site to the other, and the parasite samples from Kalimantan analysed in this research were collected before ACT was applied in the research areas. Current situation and comparison of genetic divergence among Kalimantan provinces are important subject to provide information how application of ACT influences malaria parasite populations on genotypes and polymorphisms.
Different prevalence of quintuple or sextuple mutant parasites in pfdhfr/pfdhps combined genotypes were observed (Fig. 3). While the parasites from Kalimantan and Pacitan (East Java) belonged to SP resistant quintuple or sextuple mutant genotype, parasites from the other areas in Indonesia presented four or less mutations in combined genotype that tend to adequate response for SP treatment. In addition, efficacy of SP treatment for P. vivax malaria infection in Indonesia was reported by Asih et al.  recently. These suggest SP could be considered for prevention or treatment of malaria as a single prescription or in combination with artemisinin in Indonesia except in Kalimantan Island. In such a case, regular monitoring of the efficacy and regular genotype analysis are essential to prevent spread of resistance.
Mobility of people among thousands of islands is one of the important factors for malaria control in Indonesia. Some of outbreaks, resurgences of malaria had occurred under the unique circumstances of people migration. Marwoto et al. reported immigrant workers who worked as transmigrants or seasonal workers in malaria endemic areas outside Java Island returned to their home villages brought imported malaria cases [35, 36]. Many migrants and temporal workers have moved from Pacitan (East Java) to several islands historically. Similar pfdhfr/pfdhps genotype polymorphisms in 2004–2006 parasites from Pacitan (East Java) and 2009–2012 Saruyan (Middle Kalimantan) in Table 3 suggests the malaria parasites in these districts could have been transferred along with human migrations.
The Indonesia National Malaria Control Programme desires to eliminate malaria in the whole country by 2030 . Most districts in Java and Bali islands, also several districts in other islands fall under the World Health Organization criteria of elimination and malaria cases have gradually decreased over the last several years in Indonesia .
Considering presence of malaria vector mosquitoes in most of the places used to be malaria endemic, and current situation of frequent human migrations within the country, continuous maintenance of early malaria diagnosis and treatment system is essential for malaria elimination programme. Parasite genotype analysis is helpful to follow malaria transmission and drug treatment efficiencies.
Different polymorphisms of pfdhfr/pfdhps genotypes from malaria endemic areas in Indonesia were observed in both 2004–2006 and 2009–2012 parasite isolates. The unique K540T mutations of pfdhps were detected as an SGTGA haplotype in the parasites from Kalimantan and Sumbawa islands, and East Java. The other pfdhps I588F mutation were observed as an SGEAA(588F) haplotype in Sumbawa and Lombok parasites in 2004–2006 and South Kalimantan, Gorontalo (Sulawesi) and Papua in 2009–2012. Analytical study of the prevalence and pfdhfr/pfdhps combined genotypes of K540T or I588F mutations suggested that K540T might be origin in Kalimantan Island and I588F in Sumbawa Island and then these were brought to other areas together with people movement. This research indicates regular monitoring of drug efficacy and parasite genotype analysis is important to keep efficiency and prevent the spread of resistance. It is also essential for the latest anti-malarial drug artemisinin-based combination therapy (ACT).
Sitohang V, Sariwati E, Fajariyani SB, Hwang D, Kurnia B, Hapsari RK, et al. Malaria elimination in Indonesia: halfway there. Lancet Glob Health. 2018;6:e604–6.
Poespoprodjo JR, Kenangalem E, Wafom J, Chandrawati F, Puspitasari AM, Ley B, et al. Therapeutic response to dihydroartemisinin–piperaquine for P. falciparum and P. vivax nine years after its introduction in Southern Papua, Indonesia. Am J Trop Med Hyg. 2018;98:677–82.
WHO. World Malaria Report 2017. Geneva: World Health Organization; 2017.
WHO. Recommendations on antenatal care for a positive pregnancy experience. Geneva: World Health Organization; 2016.
WHO. Intermittent preventive treatment for infants using sulfadoxine–pyrimethamine (IPTi-SP) for malaria control in Africa: implementation field guide. Geneva: Global Malaria Programme HIV/AIDS, TB, Malaria and Neglected Tropical Diseases, World Health Organization; 2011.
Roper C, Pearce R, Bredenkamp B, Gumede J, Drakeley C, Mosha F, et al. Antifolate antimalarial resistance in southeast Africa: a population-based analysis. Lancet. 2003;361:1174–81.
Mita T, Tanabe K, Kita K. Spread and evolution of Plasmodium falciparum drug resistance. Parasitol Int. 2009;58:201–9.
Muller IB, Hyde JE. Folate metabolism in human malaria parasites—75 years on. Mol Biochem Parasitol. 2013;188:63–77.
Lucchi NW, Okoth SA, Komino F, Onyona P, Goldman IF, Ljolje D, et al. Increasing prevalence of a novel triple-mutant dihydropteroate synthase genotype in Plasmodium falciparum in western Kenya. Antimicrob Agents Chemother. 2015;59:3995–4002.
Hutapea AM. Treatment of malaria quartan and prophylaxis against malaria with combined sulfadoxine–pyrimethamine in Jayapura, Indonesia. Indonesian Bull Health Res (BPK). 1979;7:9–13. http://ejournal.litbang.depkes.go.id/index.php/BPK/article/view/313. Accessed 25 June 2018.
Tjitra E, Gunawan S, Laihad F, Marwoto H, Sulaksono S, Arjoso S, et al. Evaluation of antimalarial drugs in Indonesia, 1981–1995. Bul Penelit Kesehat. 1997;25:25–58.
Laksono RD. Profilaksis Malaria di Perbatasan Indonesia-Timor Leste. Cermin Dunia Kedokteran CDK. 2011;188(38):503–7.
Cowman AF, Morry MJ, Biggs BA, Cross GA, Foote SJ. Amino acid changes linked to pyrimethamine resistance in the dihydrofolate reductase-thymidylate synthase gene of Plasmodium falciparum. Proc Natl Acad Sci USA. 1988;85:9109–13.
Peterson DS, Walliker D, Wellems TE. Evidence that a point mutation in dihydrofolate reductase-thymidylate synthase confers resistance to pyrimethamine in falciparum malaria. Proc Natl Acad Sci USA. 1988;85:9114–8.
Cowman A. The molecular basis of resistance to the sulfones, sulfonamides, and dihydrofolate reductase inhibitors. In: Sherman IW, editor. Malaria: parasite biology, pathogenesis, and protection. Washinton, DC: American Society of Microbiology Press; 1998. p. 317–29.
Triglia T, Wang P, Sims PF, Hyde JE, Cowman AF. Allelic exchange at the endogenous genomic locus in Plasmodium falciparum proves the role of dihydropteroate synthase in sulfadoxine-resistant malaria. EMBO J. 1998;17:3807815.
Gregson A, Plowe CV. Mechanisms of resistance of malaria parasites to antifolates. Pharmacol Rev. 2005;57:117–45.
Mita T, Ohashi J, Venkatesan M, Marma AS, Nakamura M, Plowe CV, et al. Ordered accumulation of mutations conferring resistance to sulfadoxine–pyrimethamine in the Plasmodium falciparum parasite. J Infect Dis. 2014;209:130–9.
Sirawaraporn W, Sathitkul T, Sirawaraporn R, Yuthavong Y, Santi DV. Antifolate-resistant mutants of Plasmodium falciparum dihydrofolate reductase. Proc Natl Acad Sci USA. 1997;94:1124–9.
Sibley CH, Hyde JE, Sims PFG, Plowe CV, Kublin JG, Mberu EK, et al. Pyrimethamine–sulfadoxine resistance in Plasmodium falciparum: what next? Trends Parasitol. 2001;17:582–8.
Omar SA, Adagu IS, Warhurst DC. Can pretreatment screening for dhps and dhfr point mutations in Plasmodium falciparum infections be used to predict sulfadoxine–pyrimethamine treatment failure? Trans R Soc Trop Med Hyg. 2001;95:315–9.
Kublin JG, Dzinjalamala FK, Kamwendo DD, Malkin EM, Cortese JF, Martino LM, et al. Molecular markers for failure of sulfadoxine–pyrimethamine and chlorproguanil-dapsone treatment of Plasmodium falciparum malaria. J Infect Dis. 2002;185:380–8.
Mugittu K, Ndejembi M, Malisa A, Lemnge M, Premji Z, Mwita A, et al. Therapeutic efficacy of sulfadoxine–pyrimethamine and prevalence of resistance markers in Tanzania prior to revision of malaria treatment policy: Plasmodium falciparum dihydrofolate reductase and dihydropteroate synthase mutations in monitoring in vivo resistance. Am J Trop Med Hyg. 2004;71:696–702.
Plowe CV. The evolution of drug-resistant malaria. Trans R Soc Trop Med Hyg. 2009;103(Suppl 1):S11–4.
Nagesha HS, Din S, Casey GJ, Susanti AI, Fryauff DJ, Reeder JC, et al. Mutations in the pfmdr1, dhfr and dhps genes of Plasmodium falciparum are associated with in vivo drug resistance in West Papua, Indonesia. Trans R Soc Trop Med Hyg. 2001;95:43–9.
Syafruddin D, Asih PB, Casey GJ, Maguire J, Baird JK, Nagesha HS, et al. Molecular epidemiology of Plasmodium falciparum resistance to antimalarial drugs in Indonesia. Am J Trop Med Hyg. 2005;72:174–81.
Asih PB, Rogers WO, Susanti AI, Rahmat A, Rozi IE, Kusumaningtyas MA, et al. Seasonal distribution of anti-malarial drug resistance alleles on the island of Sumba, Indonesia. Malar J. 2009;8:222.
Basuki S, Fitriah, Riyanto S, Budiono, Dachlan YP, Uemura H. Two novel mutations of pfdhps K540T and I588F, affecting sulphadoxine–pyrimethamine-resistant response in uncomplicated falciparum malaria at Banjar district, South Kalimantan Province, Indonesia. Malar J. 2014;13:135.
Snounou G, Viriyakosol S, Zhu XP, Jarra W, Pinheiro L, do Rosario VE, et al. High sensitivity of detection of human malaria parasites by the use of nested polymerase chain reaction. Mol Biochem Parasitol. 1993;61:315–20.
Kimura M, Kaneko O, Liu Q, Zhou M, Kawamoto F, Wataya Y, et al. Identification of the four species of human malaria parasites by nested PCR that targets variant sequences in the small subunit rRNA gene. Parasitol Int. 1997;46:91–5.
Isozumi R, Uemura H, Le DD, Truong VH, Nguyen DG, Ha VV, et al. Longitudinal survey of Plasmodium falciparum infection in Vietnam: characteristics of antimalarial resistance and their associated factors. J Clin Microbiol. 2010;48:70–7.
Lau TY, Sylvi M, William T. Mutational analysis of Plasmodium falciparum dihydrofolate reductase and dihydropteroate synthase genes in the interior division of Sabah, Malaysia. Malar J. 2013;12:445.
DitJen PPL, DepKes RI. Gebrak Malaria: Pedoman Penatalaksanaan Kasus Malaria di Indonesia. Jakarta: Ministry of Health; 2008.
Asih PB, Marantina SS, Nababan R, Lobo NF, Rozi IE, Sumarto W, et al. Distribution of Plasmodium vivax pvdhfr and pvdhps alleles and their association with sulfadoxine–pyrimethamine treatment outcomes in Indonesia. Malar J. 2015;14:365.
Marwoto HA, Sulaksono STE. Peningkatan kasus malaria di Pulau Jawa, Kepulauan Seribu dan Lampung. Media Litbang Kesehatan. 2003, 13.
Marwoto HA. Situasi Malaria Jawa Timur 1989–2002. Media Penelitian dan Pengembangan Kesehatan. 2007;17:29–41.
Ministry of Health Republic of Indonesia. Summary of National Malaria Control Program Strategic Plan 2015–2019. 2014. http://static1.1.sqspcdn.com/static/f/471029/26502872/1441109745877/Summary+of. Accessed 22 June 2018.
Kementerian Kesehatan Republik Indonesia. Profil Kesehatan Indonesia Tahun 2016. 2017. http://www.depkes.go.id/resources/download/pusdatin/profil-kesehatan-indonesia/Profil-Kesehatan-Indonesia-2016.pdf. Accessed 25 June 2018.
SB and HU prepared the protocol, summarized and interpreted data, wrote the paper. SB, F and HU analysed the samples. F, PMR, K, PA, DS, SR, I and BA coordinated field studies and collected data at the field locations. B developed protocol. YPD and HK reviewed and discussed the protocol, results, interpretation and the manuscript. All authors read and approved the final manuscript.
We thank to Prof. Dr. P.J. van den Broek, Department of Infectious Diseases, Leiden University Medical Centre, for critical advice and comments on this manuscript. We are grateful to all staff members of Health Organizations and Primary Health Centers, who were involved in investigation of malaria, sampling, and management of research resources that enabled this study.
The authors declare that they have no competing interests.
Availability of data and materials
Data presented in this study are available from the corresponding author upon reasonable request.
Consent for publication
Ethics approval and consent to participate
This study was reviewed and approved by the Ethical Committee, Faculty of Medicine, Universitas Airlangga, Surabaya, Indonesia (No. 003/EC/KEPK/FKUA/2011), and Institute of Tropical Medicine, Nagasaki University, Nagasaki, Japan (12021084).
This work was supported by JSPS KAKENHI Grant Number JP16H05817 to HU, and a research grant from Ministry of Education and Culture, Republic of Indonesia 018/SP2H/LT/DRPM/II/2016 to SB. This research was accomplished with agreement on Academic Corporation between Nagasaki University and Universitas Airlangga, and supported by JSPS-DG-RSTHE 2016 Bilateral Joint Research Program.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.