Skip to main content

Assessment of genetic polymorphisms associated with malaria antifolate resistance among the population of Libreville, Gabon



Gabon is a malaria-threatened country with a stable and hyperendemic transmission of Plasmodium falciparum monoinfection. Malaria drug resistance is widely spread in many endemic countries around the world, including Gabon. The molecular surveillance of drug resistance to antifolates and artemisinin-based combination therapy (ACT) is one of the strategies for combating malaria. As Plasmodium parasites continue to develop resistance to currently available anti-malarial drugs, this study evaluated the frequency of the polymorphisms and genetic diversity associated with this phenomenon among the parasites isolates in Gabon.


To assess the spread of resistant haplotypes among the malaria-infected population of Libreville, single nucleotide polymorphisms linked to sulfadoxine–pyrimethamine (SP) and artemisinin drugs resistance were screened for P. falciparum dihydrofolate reductase (Pfdhfr), P. falciparum dihydropteroate synthase (Pfdhps), and P. falciparum kelch 13-propeller domain (Pfk13) point mutations.


The analysis of 70 malaria-positive patient samples screened for polymorphism showed 92.65% (n = 63) mutants vs. 7.35% (n = 5) wild parasite population in Pfdhfr, with high prevalence mutations at S108N(88.24%, n = 60), N51I(85.29%, n = 58), C59R(79.41%, n = 54); however, I164L(2.94%, n = 2) showed low frequency mutation. No wild haplotype existed for Pfdhps, and there were no mutations at the K540E, A581G, and A613T/S positions. However, the mutation rate at A437G(93.38%, n = 62) was the highest, followed by S436A/F(15.38%, n = 10). A higher frequency of quadruple IRNI–SGKAA (69.84%) than quintuple IRNI–(A/F)GKAA (7.94%) mutations was observed in the PfdhfrPfdhps combination. Furthermore, none of the mutations associated with ACT resistance, especially those commonly found in Africa, were observed in Pfk13.


High polymorphism frequencies of Pfdhfr and Pfdhps genes were observed, with alternative alanine/phenylalanine mutation at S436A/F (7.69%, n = 5) for the first time. Similar to that of other areas of the country, the patterns of multiple polymorphisms were consistent with selection owing to drug pressure. Although there was no evidence of a medication failure haplotype in the studied population, ACT drug efficacy should be regularly monitored in Libreville, Gabon.


Malaria remains one of the most infectious and deadly diseases worldwide, caused by the Plasmodium parasite. The World Health Organization (WHO) reported about 247 million cases from 84 malaria endemic countries and 619,000 estimated malaria deaths in 2021 [1]. Most malaria deaths were reported in the WHO African region, with almost 76% of the total deaths recorded in children under 5 years old [1, 2]. In addition to the four African countries accounting for the highest malaria rates, malaria epidemiology in Gabon, a country in Central Africa, is characterized by a stable and hyperendemic transmission of > 90% Plasmodium falciparum monoinfection along with mixed infection of various species in a single individual [1,2,3,4]. In endemic areas, the most vulnerable target population for malaria disease includes children and pregnant women [1, 5, 6]. However, the malaria burden persists among residents of perennial transmission zones, where asymptomatic carriers are reportedly vast parasite reservoirs and are often adults because of acquired immunity within the exposure time and age [7].

Malaria antifolate drug resistance in the population of Gabon was previously described in 1995, which has gradually increased with the persistent use of these drugs [8]. Following the official implementation of artemisinin-based combination therapy (ACT) under the WHO’s recommendation to reduce the risk of drug resistance in 2003, the first-line treatment for nonsevere P. falciparum malaria has been artesunate + amodiaquine (AS + AQ) or artemether–lumefantrine (AL), and sulfadoxine–pyrimethamine (SP) has been the intermittent preventive treatment (IPT) for pregnant women [6, 9, 10]. Moreover, various studies have evidenced that the drug failure of SP is linked to point mutations, including N51I, C59R, S108N, and I164L in dihydrofolate reductase (dhfr) and S436A/F, A437G, K540E, A581G, and A613T/S in dihydropteroate synthase (dhps) [11,12,13]. The combined quintuple mutation (I51R59N108G437E540) in the Plasmodium parasite has been strongly associated with reduced parasite clearance ability, and SP showed reduced efficacy as an IPT in pregnancy [12, 14]. Since the introduction of these artemisinin-based drug combinations, polymorphisms in dhfr and dhps, which are linked to their resistance, have been continuously spreading across the nation in different regions [15,16,17,18].

Recently, mutations at the Kelch 13-propeller domain were identified to play a key role in delaying parasite clearance following ACT, especially at the C580Y, Y493H, R539T, and I543T locations, in a Southeast Asia population [19, 20]. Several studies conducted in Africa did not report these specific Asian Pfk13 mutations strongly associated with artemisinin resistance [21,22,23,24]. However, patients carried other common non-synonymous mutations, such as A578S [21, 25, 26]. These observations raise real concerns regarding the efficacy of the drugs used in combination [27, 28]. Nevertheless, Pfk13 mutation-associated artemisinin resistance in the propeller domain of the parasites is yet to be confirmed among the field samples in Gabon [18, 29]. Therefore, this study aimed to (i) screen the status of circulating dhfr, dhps, and K13 haplotypes according to the polymorphisms associated with SP and ACT resistance and (ii) estimate the frequency of the mutations (single to quintuple) of each gene and in the combination dhfrdhps genes in the Plasmodium parasites obtained from the population in Libreville, Gabon.


Study area and study design

Gabon is located in the sub-Saharan region along the Atlantic coast of Central Africa where malaria is endemic throughout the country. Libreville is the capital city of Gabon and home to one third of the country's population, located in the northwestern province of Estuaire bordering the Komo River (Fig. 1).

Fig. 1
figure 1

Localization of the study site, Libreville in Gabon

In Gabon, malaria is one of the leading causes of consultation and hospitalization, and most P. falciparum prevalence statistics are obtained from clinical data [10, 30]. Over the last couple of decades, fluctuations in the decrease or increase of malaria prevalence reported in Libreville’s urban area were 24–42% [30,31,32,33]. Moreover, there is an increased prevalence of P. falciparum asymptomatic carriage in the population and a scarcity of documentation [30, 32, 34, 35]. This study included asymptomatic and symptomatic patients who were thought to be potential carriers of malaria parasites. From the total number of malaria-confirmed carriers, only 70 samples were included in the resistance profile study. However, the final sample size (N) for each gene marker investigated was determined by considering the total number of observations obtained through sequencing and validating them twice. Thus, the frequency distribution of a haplotype was calculated by dividing the frequency of each occurrence (n) by the total number of observations (N) and multiplying by 100 to obtain a percentage value.

Study population and ethics approval

From June to July 2019, patients attending the Army Hospital (HIABO) who presented with or without malaria-related symptoms, including fever, general body weakness, chills, and headache, and who were asked by the physicians to take a malaria test were randomly enrolled according to the regulations of the ethics committee. Each patient or caregiver of the child was asked to complete and sign a written consent form before the recruitment. Thus, the study population included 94 women and 74 men between the ages of 5 months and 81 years. The study protocol was reviewed and approved by the National Ethics Committee of Gabon (PROT N°030/2018/PR/SG/CNE).

Laboratory procedures for malaria diagnosis

A rapid diagnostic test was performed using ABON™ Plus Malaria Pf/Pan Rapid Test Device (Biopharm, China), and regardless of the RDT results, blood samples of the patients were collected using dry Whatman filter paper. The spotted blood samples were dried, and each sample was placed in a single desiccant sachet for storage. Genomic DNA was extracted from the dried blood spots using DNeasy® Blood & Tissue Kit (QIAGEN, Germany) according to the manufacturer’s instructions. To confirm whether the patients had a monoinfection or multi-infection of the Plasmodium species, all the samples were amplified via conventional polymerase chain reaction (PCR) assays targeting the small subunit ribosomal RNA genes (18s) according to Singh’s protocol [36]. Only the samples showing monoinfection with P. falciparum were considered for the analysis of the resistance markers.

Detection of Pfdhfr, Pfdhps, and Pfk13mutations via nested PCR assay

In agreement with previous studies [19, 37,38,39,40], the protocols for evaluating the potential point mutations that induce anti-malarial drug resistance in Pfdhfr, Pfdhps, and Pfk13 were slightly modified. The primer pairs used for the primary and secondary PCR amplifications are listed in Additional file 1: Table S1. The mutations examined in the current study are as follows; four mutations at the codons N51I, C59R, S108N, and I164L for Pfdhfr, five mutations at the codons S436A/F, A437G, K540E, A581G, and A613T/S for Pfdhps, and six mutations at the codons Y493H, R539T, I543T, A578S, C580Y, and V589I for Pfk13.

Both primary and secondary PCR were conducted with 50 µL reaction volume containing ~ 25 ng genomic DNA, and 0.5 µL primary PCR products was used as the DNA template for the secondary reaction. All the amplification reactions were conducted using Biometra Tone 96 thermal cycler under the following conditions. For Pfdhfr, both amplification reactions were initiated with denaturation at 94 °C for 3 min followed by 40 cycles of denaturation, annealing, and extension at 94 °C/1 min, 48 °C/2 min, and 72 °C/1 min, respectively, with a final extension at 72 °C/10 min for the primary PCR. Furthermore, secondary PCR involved 94 °C/30s, 50 °C/20s, and 72 °C/1 min, with 5 min of final extension at 72 °C. The cycling conditions for Pfdhps were the same, except the annealing temperatures of 52 °C/2 min and 53 °C/20s for the primary and secondary PCR amplifications, respectively. However, PfK13 cycling was performed as described by Ariey et al. [19]. The purified genomic DNA from P. falciparum 3D7 clone and distilled water were used as positive and negative controls, respectively. The purified PCR amplicons were sent for sequencing to Macrogen (Korea) where Sanger sequencing of the reaction products was performed using ABI PRISM 3730XL Analyzer (96 capillary type). Further, the obtained sequences data were analysed using Clustal Omega for multiple sequence alignment (


Distribution of polymorphisms in dhfr, dhps, and kelch 13

Of the 168 samples, only 43.45% (n = 73) were positive for P. falciparum infection as determined via PCR targeting 18sRNA. However, because of the very poor DNA quality of 3 samples, 70 positive samples were preserved to analyse the resistance profile harboured among them. Among the 70 P. falciparum isolates, 68, 65, and 58 samples were successfully sequenced for Pfdhfr, Pfdhps, and Pfk13, respectively. After comparing the major codons mainly responsible for inducing anti-malarial drug resistance, the frequency results of the mutant alleles were obtained (Table 1). Percentages of Pfdhfr mutations were detected in 88.24% (n = 60) for the codon S108N, 85.29% (n = 58) for the codon N51I, 79.41% (n = 54) for the codon C59R, and 2.94% (n = 2) for the codon I164L. Mutations were not identified at the codons K540E, A581G, and A613T/S in Pfdhps. However, mutant alleles were detected in 93.38% (n = 62) for the codon A437G and 15.38% (n = 10) for the codon S436A/F in Pfdhps. Moreover, none of the known mutations for Pfk13 that induce artemisinin resistance were observed in the analysed samples (n = 58).

Table 1 Frequency distribution of Pfdhfr, Pfdhps, and PfK13 point mutations vs. wild-type in P. falciparum isolates from Libreville

Prevalence of Pfdhfr, Pfdhps, andPfdhfr–Pfdhps combined haplotypes

To determine the distribution of multiple mutations in the population of Libreville, haplotype analysis was performed (Table 2). Pfdhfr revealed a high mutation rate (92.65%, n = 63) against the wild-type parasites (7.35%, n = 5). Triple mutations were more frequently detected (77.95%, n = 53) than double mutations and single mutations, which were equally detected (7.35%, n = 5). The predominant IRN haplotype (76.48%, n = 52), with the N51I, C59R, and S108N codons was responsible for this finding. There was no wild-type for Pfdhps among the genotyped samples, and the single mutation was more predominant (89.25%, n = 58) than the double mutation (10.76%, n = 7). SGKAA was the highest single mutant haplotype (84.62%, n = 55) mainly due to the A437G mutation point. The analysis of the overall mutation prevalence of the combined Pfdhfr–Pfdhps genotypes included 63 fully sequenced samples, identifying 13 haplotypes. The quadruple mutation (69.84%, n = 44) was present in majority of the isolates, followed by double, triple, and quintuple mutations that were equivalent to 7.94% (n = 5) and single mutations (6.35%, n = 4). The common quadruple genotype was primarily from the combination of the N51I, C59R, S108N, and A437G mutant genotypes (63.49%, n = 40), whereas the quintuple haplotype combined the N51I, C59R, S108N–S436A/F, and A437G mutations (7.94%, n = 5).

Table 2 Frequency of haplotypes of Pfdhfr, Pfdhps, and combined Pfdhfr–Pfdhps in P. falciparum isolates

Short trends review of the resistance

Monitoring drug resistance over time remains crucial for controlling malaria. To assess the constant spread of resistant haplotypes and evaluate its trends, the polymorphism patterns observed in the study population (100%, n = 63) were compared with those reported in previous studies (Table 3). The studied population in the Bakoumba city of Gabon harbored ~ 64.3% IRN-G (quadruple) vs. 1.23% IRN-AG (quintuple) mutations before ACT implementation [41]. Following the free use of ACT among the population, compared with the quadruple mutations, the quintuple mutations were almost absent between 2005 and 2014 in certain studied areas, such as Oyem [42], Koulamoutou, Lastrouville, and Franceville [18]. However, quintuple sets of mutations were found in 2005 in Lambaréné (8%) and Libreville (17.9–22%) [15, 17, 43]. For years, malaria parasites have constantly maintained high prevalence rates of multiple sets of quadruple mutations [42,43,44,45]. The data showed that the quintuple mutation rates increased while a wide range of single nucleotide polymorphism (SNP) point mutations also simultaneously increased [41,42,43, 45].

Table 3 Short literature review of the distribution of the combined dhfr/dhps mutations in P. falciparum in the population of Gabon after introducing ACT


The trend of malaria transmission in Libreville, the capital city of Gabon, exhibits an extensively dynamic pattern throughout the rainy season and high heterogeneity according to district areas [46]. Since approximately two decades, SP administration has been limited to intermittent preventive treatment in pregnancy (IPTp)-SPs for pregnant women in Gabon, and a few studies have been conducted on pregnant women between 2005 and 2006 and 2005 vs. 2011 to determine the levels of SP-based ACT drug resistance in the population of Libreville [6, 15, 43, 47]. Additionally, a 2015–2016 study involving the adult population reported a combination of data associated with the cross resistance of SP and trimethoprim–sulfamethoxazole (cotrimazole, CTX) drugs in Oyem and Koulamoutou, including Libreville [34]. The molecular basis of the resistance of the parasites against the antifolates drugs has been clearly associated with mutations in the dhfr and dhps genes of P. falciparum [11]. Following the introduction of antifolates drugs and implementation of ACT in Gabon, a notable rise in parasite polymorphisms has been reported [15, 16]. Consequently, the current prevalence of these drug-resistant haplotypes among the population of Libreville, Gabon was investigated.

In vitro and in vivo studies reported that the S108N mutation in Pfdhfr was sufficient to confer low-level resistance toward pyrimethamine. Additive mutations at N51I and C59R caused high-level resistance toward the latter drug in combination with the mutation at I164L [17, 41]. The present study revealed high mutations rates at these positions (Table 1), which induce a moderate-to-high level resistance toward pyrimethamine. During the early use of antifolates in Gabon, the Pfdhfr triple IRN mutation rates were low before the year 2000 at Franceville and Lambaréné (13.58% and 34%, respectively) [8, 16]. A drastic increase in the mutation rates of up to 100% was observed between 2011 and 2013 at Franceville [18]. Although no wild phenotype was observed among the Lambaréné, Koulamoutou, and Franceville populations in 2014 [18], only 7.35% of the wild-type phenotype was detected among the isolates in the present investigation, which was close to that of the population of Fougamou [45]. However, the I164L mutation was detected in two isolates (2.94%), and no mutations were detected in the Pfdhfr codons A16V and C50R nor the codon S108T, which is associated with cycloguanil resistance, in the present study [16, 41]. Since the deployment of SP, a high incidence of triple mutations has been reported across the country [15, 17, 18, 34], with the triple mutant IRNI haplotype (76.48%) maintaining dominance over the single and double mutations of the studied population (7.35% each) over the years, as shown in this present study. This was consistent with earlier findings at Bakoumba (72%) [41], Oyem (72–92%) [42], and Lambaréné (92%) [17]. Similar to previously reported high levels of resistance across the country, these resistance patterns occurred because the selected parasites are constantly spreading under the influence of drugs.

Similar patterns involving high triple mutant prevalence were observed in the population of other African countries, such as Cameroon, Congo, Equatorial Guinea, Senegal, and Tanzania [23, 44, 48,49,50]. On comparing with the high multiple mutations sets documented in Nigeria [51] and Kenya [52], where various SNPs increased the double, triple, and quadruple mutation sets in Pfdhfr, likewise NRNL, IRSI, ICNI, and ICNL multiple mutation sets were expressed among the current studied population. In addition, contrary to the I164U finding in individuals with malaria infection living with human immunodeficiency virus, I164L mutation, which is rarely reported in Central Africa [34, 49, 53] and somewhat reported in East and South Africa and Asia [54,55,56], was detected in two isolates of the Gabon population of the current study. Thus, these findings indicate the importance of regularly checking the status of the polymorphisms that induce drug resistance.

Regarding Pfdhps, no wild phenotype was detected in the studied population. Notably, after the 2005–2007 period, Pfdhps wild-type was not detected in the screened populations [16,17,18, 45], and similarly, a shift from the high selection of the S436A mutation to that of the A437G mutation following the full implementation of IPTs-SP and ACT was reported [8, 16,17,18]. Herein, the frequency of the A437G single mutation was higher than that of the alternative alanine/phenylalanine mutation (S436 (A/F)), as shown in the results (93.38% vs. 15.38%, respectively). In contrast to the findings of the isolates from Kenya, the present results did not detect S436H mutation [52], instead S436F mutation was detected for the first time in five samples among the current studied population of Gabon, similar to that in the population of other countries, such as Sierra Leone, Kenya, and Vietnam [11, 57]. Moreover, it was particularly detected as a FGKAA double mutation (10.67%) (Table 2) in the current study population, whereas no triple or quadruple mutations were observed in dhps.

Herein, Pfdhfr and Pfdhps mutations were commonly identified as multiple mutations than single mutations in the isolates. However, within 10 years following the introduction of ACT and IPTp–SP and its daily use as prescription drugs, a number of investigations highlighted two- to three-fold rates of multiple SNPs among the parasite population owing to high levels of triple and quadruple mutations. However, quintuple or higher combined mutations were not recorded or were scarce [9, 15, 18, 42]. Herein, no sextuple or septuplet genotypes were noted in combined dhfr–dhps; however, the quintuple mutations were due to the combinations of IRNI–AGKAA and IRNI–FGKAA. These results highlight the fact that within the time of use, the number of haplotypes increases under the drug pressure. An in vitro study reported that the A437G mutant of Pfdhps exhibited a lesser degree of tolerance toward sulfadoxine than the double A437G–K540E mutant of Pfdhps. Therefore, quadrupleIRNI-G mutants exhibit a less noxious effect in parasites than quintuple mutants associated A437G–K540E mutations [58]. Herein, no mutations were detected at the 540 position (Table 1), and neither the quintupleIRNI–SGEAA nor the sextupleIRNL–SGEAA mutations were linked to SP drug failure (Table 2). These outcomes support the continued effectiveness of SP.

The emergence of artemisinin resistance has been reported in Myanmar, Vietnam, Laos, and China [19,20,21]. In fact, a wide range of Pfk13-propeller domain mutations linked to artemisinin resistance played a role in parasite clearance [19, 59]. However, A578S mutation was the most widespread in sub-Saharan Africa [52]. ACT is administered as the first-line treatment for non-severe malarial infection in Gabon from 2003 [47], and evidence of malarial morbidity decline in Libreville up to 2008 was reported following free access to and the huge distribution of ACT [31]. Among all the investigated samples, no mutations that have previously been linked to artemisinin resistance were noted. Furthermore, few non-synonymous mutations were identified in some areas of Gabon. The results of this study are consistent with previous findings [18, 45]. The current study included a small number of patients from one healthcare facility; thus, it may be necessary to expand the screening population to include patients from various health centers. This will potentially enhance the representative study scale of the situation of genetic polymorphisms associated with resistance to these drugs in the population of Libreville, Gabon.


High frequencies of Pfdhfr and Pfdhps mutant haplotypes among the studies population were observed. The constant rise of mutations associated with SP drug resistance more than 20 years following its use for treating uncomplicated malaria cases that predated the implementation of IPTp–SP is alarming. However, despite several investigations, the haplotypes associated with drug failure are yet to be established. The constant increase of drugs resistant haplotypes sets across time and place over the country, constitute a permanent danger leading to drugs failure. Thus, to control the epidemiology accompanied by the administration of the treatment drugs, regular monitoring of drug efficacy must be mandated across the country.

Availability of data and materials

All data generated or analysed during this study are included in this published article.



Artemisinin-based combination therapy






Dihydrofolate reductase


Dihydropteroate synthase


Intermittent preventive treatment in pregnancy

Pf k13:

Plasmodium falciparum kelch 13




World Health Organization


  1. WHO. World malaria report 2022. Geneva: World Health Organization; 2022.

    Google Scholar 

  2. WHO. Regional data and trends briefing kit. Geneva: World malaria report; 2022.

    Google Scholar 

  3. Lalremruata A, Jeyaraj S, Engleitner T, Joanny F, Lang A, Belard S, et al. Species and genotype diversity of Plasmodium in malaria patients from Gabon analysed by next generation sequencing. Malar J. 2017;16:398.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Woldearegai TG, Lalremruata A, Nguyen TT, Gmeiner M, Veletzky L, Tazemda-Kuitsouc GB, et al. Characterization of Plasmodium infections among inhabitants of rural areas in Gabon. Sci Rep. 2019;9:9784.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Maghendji-Nzondo S, Nzoughe H, Lemamy GJ, Kouna LC, Pegha-Moukandja I, Lekoulou F, et al. Prevalence of malaria, prevention measures, and main clinical features in febrile children admitted to the Franceville Regional Hospital, Gabon. Parasite. 2016;23:32.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Bouyou-Akotet MK, Nzenze-Afene S, Ngoungou EB, Kendjo E, Owono-Medang M, Lekana-Douki JB, et al. Burden of malaria during pregnancy at the time of IPTp/SP implementation in Gabon. Am J Trop Med Hyg. 2010;82:202–9.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Dal-Bianco MP, Koster KB, Kombila UD, Kun JF, Grobusch MP, Ngoma GM, et al. High prevalence of asymptomatic Plasmodium falciparum infection in gabonese adults. Am J Trop Med Hyg. 2007;77:939–42.

    Article  CAS  PubMed  Google Scholar 

  8. Kun JF, Lehman LG, Lell B, Schmidt-Ott R, Kremsner PG. Low-dose treatment with sulfadoxine-pyrimethamine combinations selects for drug-resistant Plasmodium falciparum strains. Antimicrob Agents Chemother. 1999;43:2205–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Ramharter M, Schuster K, Bouyou-Akotet MK, Adegnika AA, Schmits K, Mombo-Ngoma G, et al. Malaria in pregnancy before and after the implementation of a national IPTp program in Gabon. Am J Trop Med Hyg. 2007;77:418–22.

    Article  PubMed  Google Scholar 

  10. Mawili-Mboumba DP, Bouyou Akotet MK, Kendjo E, Nzamba J, Medang MO, Mbina JR, et al. Increase in malaria prevalence and age of at risk population in different areas of Gabon. Malar J. 2013;12:3.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Wang P, Lee CS, Bayoumi R, Djimde A, Doumbo O, Swedberg G, et al. Resistance to antifolates in Plasmodium falciparum monitored by sequence analysis of dihydropteroate synthetase and dihydrofolate reductase alleles in a large number of field samples of diverse origins. Mol Biochem Parasitol. 1997;89:161–77.

    Article  CAS  PubMed  Google Scholar 

  12. Gesase S, Gosling RD, Hashim R, Ord R, Naidoo I, Madebe R, et al. High resistance of Plasmodium falciparum to sulphadoxine/pyrimethamine in northern Tanzania and the emergence of dhps resistance mutation at Codon 581. PLoS ONE. 2009;4:e4569.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Kayode AT, Ajogbasile FV, Akano K, Uwanibe JN, Oluniyi PE, Eromon PJ, et al. Polymorphisms in Plasmodium falciparum dihydropteroate synthetase and dihydrofolate reductase genes in nigerian children with uncomplicated malaria using high-resolution melting technique. Sci Rep. 2021;11:471.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Desai M, Gutman J, Taylor SM, Wiegand RE, Khairallah C, Kayentao K, et al. Impact of sulfadoxine-pyrimethamine resistance on effectiveness of intermittent preventive therapy for malaria in pregnancy at clearing infections and preventing low birth weight. Clin Infect Dis. 2016;62:323–33.

    Article  CAS  PubMed  Google Scholar 

  15. Bouyou-Akotet MK, Mawili-Mboumba DP, Tchantchou Tde D, Kombila M. High prevalence of sulfadoxine/pyrimethamine-resistant alleles of Plasmodium falciparum isolates in pregnant women at the time of introduction of intermittent preventive treatment with sulfadoxine/pyrimethamine in Gabon. J Antimicrob Chemother. 2010;65:438–41.

    Article  CAS  PubMed  Google Scholar 

  16. Mawili-Mboumba DP, Ekala MT, Lekoulou F, Ntoumi F. Molecular analysis of DHFR and DHPS genes in P. falciparum clinical isolates from the Haut–Ogooué region in Gabon. Acta Trop. 2001;78:231–40.

    Article  CAS  PubMed  Google Scholar 

  17. Mombo-Ngoma G, Oyakhirome S, Ord R, Gabor JJ, Greutélaers KC, Profanter K, et al. High prevalence of dhfr triple mutant and correlation with high rates of sulphadoxine-pyrimethamine treatment failures in vivo in gabonese children. Malar J. 2011;10:123.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Voumbo-Matoumona DF, Kouna LC, Madamet M, Maghendji-Nzondo S, Pradines B, Lekana-Douki JB. Prevalence of Plasmodium falciparum antimalarial drug resistance genes in Southeastern Gabon from 2011 to 2014. Infect Drug Resist. 2018;11:1329–38.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. 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  Google Scholar 

  20. 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  CAS  PubMed  Google Scholar 

  21. 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  Google Scholar 

  22. Taylor SM, Parobek CM, DeConti DK, Kayentao K, Coulibaly SO, Greenwood BM, et al. Absence of putative artemisinin resistance mutations among Plasmodium falciparum in sub-saharan Africa: a molecular epidemiologic study. J Infect Dis. 2015;211:680–8.

    Article  CAS  PubMed  Google Scholar 

  23. Boussaroque A, Fall B, Madamet M, Camara C, Benoit N, Fall M, et al. Emergence of mutations in the K13 propeller gene of Plasmodium falciparum isolates from Dakar, Senegal, in 2013–2014. Antimicrob Agents Chemother. 2016;60:624–7.

    Article  CAS  PubMed  Google Scholar 

  24. Muwanguzi J, Henriques G, Sawa P, Bousema T, Sutherland CJ, Beshir KB. Lack of K13 mutations in Plasmodium falciparum persisting after artemisinin combination therapy treatment of kenyan children. Malar J. 2016;15:36.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Torrentino-Madamet M, Fall B, Benoit N, Camara C, Amalvict R, Fall M, et al. Limited polymorphisms in k13 gene in Plasmodium falciparum isolates from Dakar, Senegal in 2012–2013. Malar J. 2014;13:472.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Maïga-Ascofaré O, May J. Is the A578S single-nucleotide polymorphism in K13-propeller a marker of emerging resistance to artemisinin among Plasmodium falciparum in Africa? J Infect Dis. 2016;213:165–6.

    Article  PubMed  Google Scholar 

  27. Roper C, Alifrangis M, Ariey F, Talisuna A, Menard D, Mercereau-Puijalon O, et al. Molecular surveillance for artemisinin resistance in Africa. Lancet Infect Dis. 2014;14:668–70.

    Article  PubMed  Google Scholar 

  28. Sutherland CJ, Lansdell P, Sanders M, Muwanguzi J, van Schalkwyk DA, Kaur H, et al. pfk13-Independent treatment failure in four imported cases of Plasmodium falciparum malaria treated with artemether-lumefantrine in the United Kingdom. Antimicrob Agents Chemother. 2017;61:e02382–16.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Kamau E, Campino S, Amenga-Etego L, Drury E, Ishengoma D, Johnson K, et al. K13-propeller polymorphisms in Plasmodium falciparum parasites from sub-saharan Africa. J Infect Dis. 2015;211:1352–5.

    CAS  PubMed  Google Scholar 

  30. Bouyou-Akotet MK, Offouga CL, Mawili-Mboumba DP, Essola L, Madoungou B, Kombila M. Falciparum malaria as an emerging cause of fever in adults living in Gabon, Central Africa. Biomed Res Int. 2014;2014:351281.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Bouyou-Akotet MK, Mawili-Mboumba DP, Kendjo E, Mabika-Mamfoumbi M, Ngoungou EB, Dzeing-Ella A, et al. Evidence of decline of malaria in the general hospital of Libreville, Gabon from 2000 to 2008. Malar J. 2009;8:300.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Mbang Nguema OA, Obiang Ndong GP, M’Bondoukwe NP, Ndong Ngomo JM, Koumba Lengongo JV, Pongui Ngondza B, et al. Prevalence of asymptomatic malaria in urban and semi-urban areas in 2016 and 2019, to Gabon. Eur J Applied Sci. 2022;10:62–72.

    Google Scholar 

  33. Bouyou-Akotet MK, Nzenze-Afene S, Mawili-Mboumba DP, Owono-Medang M, Guiyedi V, Kombila M. [Trends in the prevalence of malaria and anemia at delivery in Libreville from 1995 to 2011](in French). Sante. 2011;21:199–203.

    PubMed  Google Scholar 

  34. Koumba Lengongo JV, Ndiaye YD, Tshibola Mbuyi ML, Ndong Ngomo JM, Ndiaye D, Bouyou Akotet MK, et al. Increased frequency of Pfdhps A581G mutation in Plasmodium falciparum isolates from gabonese HIV-infected individuals. Malar Res Treat. 2019;2019:9523259.

    PubMed  PubMed Central  Google Scholar 

  35. Pegha Moukandja I, Biteghe Bi Essone JC, Sagara I, Kassa Kassa RF, Ondzaga J, Lékana Douki JB, et al. Marked rise in the prevalence of asymptomatic Plasmodium falciparum infection in rural Gabon. PLoS ONE. 2016;11:e0153899.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Snounou G, Viriyakosol S, Jarra W, Thaithong S, Brown KN. Identification of the four human malaria parasite species in field samples by the polymerase chain reaction and detection of a high prevalence of mixed infections. Mol Biochem Parasitol. 1993;58:283–92.

    Article  CAS  PubMed  Google Scholar 

  37. Pearce RJ, Drakeley C, Chandramohan D, Mosha F, Roper C. Molecular determination of point mutation haplotypes in the dihydrofolate reductase and dihydropteroate synthase of Plasmodium falciparum in three districts of northern Tanzania. Antimicrob Agents Chemother. 2003;47:1347–54.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Parola P, Pradines B, Simon F, Carlotti MP, Minodier P, Ranjeva MP, et al. Antimalarial drug susceptibility and point mutations associated with drug resistance in 248 Plasmodium falciparum isolates imported from Comoros to Marseille, France in 2004 2006. Am J Trop Med Hyg. 2007;77:431–7.

    Article  CAS  PubMed  Google Scholar 

  39. Boussaroque A, Fall B, Madamet M, Wade KA, Fall M, Nakoulima A, et al. Prevalence of anti-malarial resistance genes in Dakar, Senegal from 2013 to 2014. Malar J. 2016;15:347.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Torrentino-Madamet M, Collet L, Lepère JF, Benoit N, Amalvict R, Ménard D, et al. K13-Propeller polymorphisms in Plasmodium falciparum isolates from patients in Mayotte in 2013 and 2014. Antimicrob Agents Chemother. 2015;59:7878–81.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Aubouy A, Jafari S, Huart V, Migot-Nabias F, Mayombo J, Durand R, et al. DHFR and DHPS genotypes of Plasmodium falciparum isolates from Gabon correlate with in vitro activity of pyrimethamine and cycloguanil, but not with sulfadoxine-pyrimethamine treatment efficacy. J Antimicrob Chemother. 2003;52:43–9.

    Article  CAS  PubMed  Google Scholar 

  42. Ndong Ngomo J-M, Mawili-Mboumba DP, M’Bondoukwe NP, Nikiéma Ndong Ella R, Bouyou Akotet MK. Increased prevalence of mutant allele Pfdhps 437G and pfdhfr triple mutation in Plasmodium falciparum isolates from a rural area of Gabon, three years after the change of malaria treatment policy. Malar Res Treat. 2016;2016:9694372–72.

    PubMed  PubMed Central  Google Scholar 

  43. Bouyou-Akotet MK, Tshibola ML, Mawili-Mboumba DP, Nzong J, Bahamontes-Rosa N, Tsoumbou-Bakana G, et al. Frequencies of dhfr/dhps multiple mutations and Plasmodium falciparum submicroscopic gametocyte carriage in gabonese pregnant women following IPTp-SP implementation. Acta Parasitol. 2015;60:218–25.

    Article  CAS  PubMed  Google Scholar 

  44. Voumbo-Matoumona DF, Akiana J, Madamet M, Kouna LC, Lekana-Douki JB, Pradines B. High prevalence of Plasmodium falciparum antimalarial drug resistance markers in isolates from asymptomatic patients from the Republic of the Congo between 2010 and 2015. J Glob Antimicrob Resist. 2018;14:277–83.

    Article  PubMed  Google Scholar 

  45. Boukoumba FM, Lekana-Douki JB, Matsiegui PB, Moukodoum DN, Adegnika AA, Oyegue-Liabagui SL. High prevalence of genotypes associated with sulfadoxine/pyrimethamine resistance in the rural area of Fougamou, Gabon. J Glob Antimicrob Resist. 2021;25:181–6.

    Article  CAS  PubMed  Google Scholar 

  46. Mourou JR, Coffinet T, Jarjaval F, Cotteaux C, Pradines E, Godefroy L, et al. Malaria transmission in Libreville: results of a one year survey. Malar J. 2012;11:40.

    Article  PubMed  PubMed Central  Google Scholar 

  47. Ministère de la Santé. Rapport de l’atelier national de consensus sur les perspectives thérapeutiques du paludisme. Programme national de lutte contre le paludisme, Gabon. 2003: 13–14.

  48. Alifrangis M, Lusingu JP, Mmbando B, Dalgaard MB, Vestergaard LS, Ishengoma D, et al. Five-year surveillance of molecular markers of Plasmodium falciparum antimalarial drug resistance in Korogwe District, Tanzania: accumulation of the 581G mutation in the P. falciparum dihydropteroate synthase gene. Am J Trop Med Hyg. 2009;80:523–7.

    Article  CAS  PubMed  Google Scholar 

  49. Tuedom AGB, Sarah-Matio EM, Moukoko CEE, Feufack-Donfack BL, Maffo CN, Bayibeki AN, et al. Antimalarial drug resistance in the Central and Adamawa regions of Cameroon: prevalence of mutations in P. falciparum crt, Pfmdr1, Pfdhfr and Pfdhps genes. PLoS ONE. 2021;16:e0256343.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Guerra M, Neres R, Salgueiro P, Mendes C, Ndong-Mabale N, Berzosa P, et al. Plasmodium falciparum genetic diversity in Continental Equatorial Guinea before and after introduction of artemisinin-based combination therapy. Antimicrob Agents Chemother. 2017;61:e02556–15.

    Article  CAS  PubMed  Google Scholar 

  51. Quan H, Igbasi U, Oyibo W, Omilabu S, Chen S-B, Shen H-M, et al. High multiple mutations of Plasmodium falciparum-resistant genotypes to sulphadoxine-pyrimethamine in Lagos, Nigeria. Infect Dis Poverty. 2020;9:91–1.

    Article  PubMed  PubMed Central  Google Scholar 

  52. Pacheco MA, Schneider KA, Cheng Q, Munde EO, Ndege C, Onyango C, et al. Changes in the frequencies of Plasmodium falciparum dhps and dhfr drug-resistant mutations in children from western Kenya from 2005 to 2018: the rise of Pfdhps S436H. Malar J. 2020;19:378.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Berzosa P, Esteban-Cantos A, García L, González V, Navarro M, Fernández T, et al. Profile of molecular mutations in pfdhfr, pfdhps, pfmdr1, and pfcrt genes of Plasmodium falciparum related to resistance to different anti-malarial drugs in the Bata District (Equatorial Guinea). Malar J. 2017;16:28.

    Article  PubMed  PubMed Central  Google Scholar 

  54. Lynch CA, Pearce R, Pota H, Egwang C, Egwang T, Bhasin A, et al. Travel and the emergence of high-level drug resistance in Plasmodium falciparum in southwest Uganda: results from a population-based study. Malar J. 2017;16:150.

    Article  PubMed  PubMed Central  Google Scholar 

  55. Kaingona-Daniel EP, Gomes LR, Gama BE, Almeida-de-Oliveira NK, Fortes F, Ménard D, et al. Low-grade sulfadoxine-pyrimethamine resistance in Plasmodium falciparum parasites from Lubango, Angola. Malar J. 2016;15:309.

    Article  PubMed  PubMed Central  Google Scholar 

  56. Basuki S, Fitriah, Risamasu PM, Kasmijati, Ariami P, Riyanto S, et al. Origins and spread of novel genetic variants of sulfadoxine-pyrimethamine resistance in Plasmodium falciparum isolates in Indonesia. Malar J. 2018;17:475.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Iriemenam NC, Shah M, Gatei W, van Eijk AM, Ayisi J, Kariuki S, et al. Temporal trends of sulphadoxine-pyrimethamine (SP) drug-resistance molecular markers in Plasmodium falciparum parasites from pregnant women in western Kenya. Malar J. 2012;11:134.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  58. Griffin JT, Cairns M, Ghani AC, Roper C, Schellenberg D, Carneiro I, et al. Protective efficacy of intermittent preventive treatment of malaria in infants (IPTi) using sulfadoxine-pyrimethamine and parasite resistance. PLoS ONE. 2010;5:e12618.

    Article  PubMed  PubMed Central  Google Scholar 

  59. WWARN K13 Genotye-Phenotype Study Group. Association of mutations in the Plasmodium falciparum Kelch13 gene (Pf3D7_1343700) with parasite clearance rates after artemisinin-based treatments-a WWARN individual patient data meta-analysis. BMC Med. 2019;17:1.

    Article  Google Scholar 

Download references


We thank all the members of the two teams for their invaluable contributions, especially Litchangou Fred, Bignoumba Corliane, Angomo Désiré Wasslin, and Ngondet Mpoumboue Igor Eudes at the Service de Santé Militaire, for their assistance and technical support. We are also grateful to all the participants for their approval.


This research was financially supported by the Basic Science Research Program (NRF-2019R1C1C1002170), the framework of the international cooperation program (2022K2A9A1A01098057, FY2022), and Brain Pool program (NRF-2018H1D3A1A02074759) funded by the Ministry of Science, ICT and Future Planning through the National Research Foundation of Korea (NRF).

Author information

Authors and Affiliations



SDDB, BAI and YKG conceived the study. BAI supervised the fieldwork in Gabon. LYM and FLA conducted the samples collection and administrative part. SDDB, ZM, and JAM investigated. SL, DIL, YH and YKG analyzed the data. SDDB and YKG wrote a first draft to the published version of the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Youn-Kyoung Goo.

Ethics declarations

Ethics approval and consent to participate

The study was reviewed and approved by the National Ethics Committee of Gabon (PROT N°030/2018/PR/SG/CNE). Written informed consent was obtained from each participant.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

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

Supplementary information

Additional file 1:  Table S1

. Primer pairs used for theprimary and secondary amplification of drug resistance genes.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Dinzouna-Boutamba, SD., Iroungou, B.A., Akombi, F.L. et al. Assessment of genetic polymorphisms associated with malaria antifolate resistance among the population of Libreville, Gabon. Malar J 22, 183 (2023).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: