Molecular assessment of kelch13 non-synonymous mutations in Plasmodium falciparum isolates from Central African Republic (2017–2019)
Malaria Journal volume 19, Article number: 191 (2020)
Over the last decade, artemisinin-based combination therapy (ACT) has contributed substantially to the decrease in malaria-related morbidity and mortality. The emergence of Plasmodium falciparum parasites resistant to artemisinin derivatives in Southeast Asia and the risk of their spread or of local emergence in sub-Saharan Africa are a major threat to public health. This study thus set out to estimate the proportion of P. falciparum isolates, with Pfkelch13 gene mutations associated with artemisinin resistance previously detected in Southeast Asia.
Blood samples were collected in two sites of Bangui, the capital of the Central African Republic (CAR) from 2017 to 2019. DNA was extracted and nested PCR were carried out to detect Plasmodium species and mutations in the propeller domain of the Pfkelch13 gene for P. falciparum samples.
A total of 255 P. falciparum samples were analysed. Plasmodium ovale DNA was found in four samples (1.57%, 4/255). Among the 187 samples with interpretable Pfkelch13 sequences, four samples presented a mutation (2.1%, 4/187), including one non-synonymous mutation (Y653N) (0.5%, 1/187). This mutation has never been described as associated with artemisinin resistance in Southeast Asia and its in vitro phenotype is unknown.
This preliminary study indicates the absence of Pfkelch13 mutant associated with artemisinin resistance in Bangui. However, this limited study needs to be extended by collecting samples across the whole country along with the evaluation of in vitro and in vivo phenotype profiles of Pfkelch13 mutant parasites to estimate the risk of artemisinin resistance in the CAR.
The introduction of artemisinin-based combination therapy (ACT) as the first-line treatment for Plasmodium falciparum malaria worldwide has greatly helped to reduce the morbidity and the mortality due to malaria from 2010 (251 million cases and 585,000 deaths) to 2017 (228 million cases and 405,000 deaths) [1, 2]. The global implementation of ACT has thus inspired hope, but soon became an issue of concern due to the emergence of artemisinin-resistant parasites in 2006–2007 in the Greater Mekong subregion, and the risk of spreading of these resistant strains to sub-Saharan Africa where malaria transmission and burden are high [3, 4]. Since then, many sub-Saharan African countries have enhanced their surveillance system to assess the efficacy of ACT in clinical drug efficacy studies and to detect artemisinin resistant parasites by in vitro susceptibility testing or Pfkelch13 genotyping [5, 6].
In the Central African Republic (CAR), ACT has been used since 2006 as first- and second-line treatments . Set up in October 2016 and based on the WHO Global Technical Strategy for Malaria 2016–2030, the national CAR policies for malaria control emphasize surveillance of clinical efficacy, in vitro susceptibility testing and screening for molecular markers associated with anti-malarial drug resistance. Nonetheless, there are very few reliable data available. Two studies on in vitro sensitivity tests conducted in Bangui in 2004 (before ACT was introduced) and in 2014 showed that 100% of the circulating P. falciparum strains were sensitive to the main artemisinin-based drug combinations [8, 9]. A clinical study assessing the efficacy of ACT was conducted in Bangui in 2010 and showed a high level of clinical and parasitological response rate estimated to 100% for the artemether–lumefantrine (AL) and artesunate–amodiaquine (ASAQ) . The sole study assessing the prevalence of mutant Pfkelch13 parasites was conducted in 2014 in Bangui and showed a frequency of 4.5% non-synonymous mutations (Q468R, W470Stop, K480R, L505S, Y519C, S522C, N537D, G545E, I552M, W565Stop, V566I, S577P, A578S, F583L, V589I, G591D, E606G, E612G, Q633R, I640V, D641G), all not validated to confer artemisinin resistance .
In addition to the lack of data on artemisinin resistance, since 2012, the CAR has been experiencing unprecedented social unrest and political instability, leading to the arrival of thousands of expatriate civilians and foreign troops, some of whom come from countries located in known multidrug-resistant malaria areas. The presence in the CAR of troops from Cambodia and Thailand, two countries located in the area where artemisinin-resistance have emerged, and from Bangladesh, Bhutan and Nepal, countries that neighbour this epicentre of emergence, may lead to the potential spread of multidrug-resistant P. falciparum in Africa in general, and to the CAR in particular.
This concern is further heightened by studies (Ménard, pers. commun.) that have shown a proportion of 0.9% asymptomatic carriers in the Cambodian military personnel sent to Africa for peacekeeping missions (18 P. falciparum infections among the 1950 military personnel tested between June 2014 and June 2017). Moreover, there are a multiplicity of factors that may favour the emergence, the introduction and—more importantly, the selection—of drug-resistant malaria strains, particularly those resistant to artemisinin derivatives. These factors include population migration within the CAR or emigration to neighbouring countries, the deplorable living conditions in displaced persons and refugee camps, the illicit trafficking of fake anti-malarial drugs in local markets. Since 2014, several single non-synonymous mutations in the propeller domain of the Pfkelch13 gene have been associated with resistance to artemisinin derivatives, defined clinically by a delayed clearance of P. falciparum during the three-day course of ACT or by the increase in the number of parasites (ring stages) resistant to a pulse of 700 nM of DHA as expressed in the ring-stage survival assay (RSA) . Since then, in vivo clinical drug efficacy study and in vitro parasite susceptibility testing along with screening for specific mutations in the Pfkelch13 gene are the recommended approach for the surveillance of the efficacy of ACT . However, in vivo and in vitro approaches have some practical issues regarding their elaborate protocols and the follow-up of patients for more than 1 month post-infection. Therefore, screening for mutations to detect potential resistance markers may be a useful, efficient alternative. This study thus set out to investigate kelch13 polymorphism in P. falciparum isolates collected in Bangui, the capital of the CAR.
Study site and period
The study was carried out in Bangui, where malaria transmission is holoendemic, on samples taken from symptomatic malaria patients at two health centres during two different periods: (1) at the Institut Pasteur in Bangui (IPB) between September 2017 and February 2018 and (2) at the Bangui Paediatric Complex (BPC) between November 2018 and March 2019. Both centres are located in the 1st district of Bangui, but patients come from all districts.
Study population and sampling
The study population was made up of patients visiting the IPB and the BPC for malarial symptoms. Patients included in the study showed positive Giemsa-stained thick blood smears, and blood samples and demographic data were available for them. All blood samples were collected in EDTA blood collection tubes. Some of the sampled blood was then spotted on filter paper for further analyses.
Plasmodium species identification
DNA was extracted using the Chelex-100 method on the dried blood spot samples . Extracted DNA was used to identify Plasmodium species using the technique described by Singh et al. . The targeted 18S (SSU) rRNA gene common to all four Plasmodium species was amplified using a specific primer pair (PCR1) and then species-specific primers were used to screen for each individual Plasmodium species (PCR2) (Table 1). After migration, PCR products were observed under UV light and the bands were compared with the positive controls for species identification.
Detection of mutations in the Pfkelch13 gene
A portion of the Pfkelch13 gene was amplified from the DNA extract on confirmed P. falciparum samples using the method described by Ariey et al. . Briefly, amplification of the Pfkelch13-propeller domain (codons 440–680, 720 bp) was performed as following. Five μl DNA was amplified with 0.25 μM each primer, 0.2 mM dNTP, 2.5 mM MgCl2, and 1.25 U Taq DNA polymerase (Solis Biodyne, Estonia), in 25 µl volume using the following cycling program: 15 min at 95 °C, then 35 cycles of 30 s at 95 °C, 2 min at 58 °C, 2 min at 72 °C, and final extension 10 min at 72 °C. For the nested PCR, 5 μl of primary PCR products were amplified under the same conditions, except for annealing and extension (1 min). PCR products were detected using 2% agarose gel electrophoresis and ethidium bromide staining. Double strand sequencing of PCR products was performed by Eurofins (Germany). Sequences were analysed with the CLC Main Workbench 20 software. The 3D7 strain of P. falciparum (PF3D7_1343700) was used as the reference sequence, to identify polymorphism. All electropherograms were visualized to detect isolates with mixed alleles that were considered to be mutated for the purpose of mutation-frequency estimation. The quality control was assessed by including three blinded quality-control samples in each 96-well sequencing plate.
Data processing and statistical analyses
The data were compiled in a Microsoft Excel spreadsheet (ver. MS Office 2010). The data were analysed using descriptive statistics (mean, frequency standard deviation and the confidence intervals).
Demographic characteristics of the study population
The mean age of the patients included in the study was 9.17 years (range: 2 to 71 years) at IPB and 3.75 years (range: 0.12 to 15 years) at the BPC. The difference of the mean age of the populations seeking anti-malarial treatment between the two sites was due that BPC is a pediatric hospital where only patients under 16 years of age are accepted while IPB is opened to the overall population. The sex ratio (M/F) was 1.5 at the IPB and 1.15 at the CPB. The mean parasite density was five times higher in patients recruited at the IPB than those at the CPB (Fig. 1).
Prevalence of Plasmodium species
Of the positive smear samples, 255 could be PCR-amplified using the Plasmodium primers: 100% (255/255) were positive for P. falciparum and 1.57% (4/255) for P. ovale. At the IPB, P. falciparum prevalence was 100% (100/100) and P. ovale prevalence was 3% (3/100). At the BPC, P. falciparum prevalence was 100% (155/155) and P. ovale prevalence was 0.65% (1/155). The other two Plasmodium species (Plasmodium malariae and Plasmodium vivax) were not observed.
Frequency of mutations in the Pfkelch13 gene
A total of 192 amplicons were sequenced to screen for mutations. Among the 187 interpretable sequences, four mutations (2.14%) were detected, including one non-synonymous mutation (Y653N, 0.54%) and three synonymous mutations (C469C, D464D, A627A, 1.6%) (GenBank Accession Numbers: MT434117-MT434120). The Y653N mutation was observed in a sample collected in 2019 at the CPB (Table 2). The Y653N mutation is located on blade 5 of the Pfkelch13 propeller domain (Fig. 2). The model predicts that the substitution from aromatic (Y) to polar, non-charged (N) residue at 653 position.
The present study provides recent data on the Plasmodium species circulating in Bangui, CAR, as well as on the presence of parasites with Pfkelch13 mutations.
First, the molecular detection of Plasmodium DNA in studied isolates showed that both P. falciparum and P. ovale are circulating in Bangui. Although P. falciparum was found in all malaria-positive cases, this species was associated with P. ovale in 1.57% of cases. This figure is higher than that of previous observations. A study carried out in 2010 in Bangui estimated a prevalence of 0.3% P. ovale . Of note, a co-infection of P. ovale, P. falciparum and P. malariae was observed in Rouen, France in 2017 in an imported malaria case in two children from the CAR . In the other imported malaria study, 4 cases (4/200, i.e. a prevalence of 2%) of P. ovale were diagnosed among Peruvian peace-keepers deployed in support of United Nations operations in the CAR from 2016 to 2017 after to return to Peru .
Since 2006–2007, resistance to artemisinin has emerged in Southeast Asia, along the Thai–Cambodian border. This resistance, first identified as an increase in parasite clearance times after treatment with artesunate monotherapy or after ACT, is now better understood. It involves early ring-stages that resist treatment by ceasing to grow when exposed to the drug. This phenomenon has been demonstrated by the development of a new in vitro test called the ring-stage survival assay (RSA). Resistant parasites show a proportion of > 1% of parasites that survive after 72 h compared with susceptible parasites (< 1%). Since 2014, these two phenotypes (clinical and in vitro) have been clearly associated with the presence of several non-synonymous mutations in the propeller domain of the Pfkelch13 gene. The two main hotspots of emergence are located in the Greater Mekong subregion in Southeast Asia, where the parasites carrying the C580Y or the F446I mutations are now dominant . Other single point mutations (N458Y, Y493H, R539T, I543T, M476I, P553L and R561H) have also been validated as conferring resistance to artemisinin . In addition, some mutations (P441L, G449A, C469F, P527H, N537I, G538V, V568G, P574L, F673I and A675V) are candidates suspected to be associated with artemisinin resistance . In sub-Saharan Africa, validated or candidate mutations associated with resistance (R539T, P574L) have been observed in Angola, Equatorial Guinea and in Rwanda, whereas mutations potentially associated with artemisinin resistance (M476I) have been detected in Senegal [20,21,22,23]. A particular case was the observation of a local mutation (M579I) in Equatorial Guinea . This mutation was shown to be associated with artemisinin resistance, while the A578S mutation, common in Africa, was not found to be associated with artemisinin resistance . Here, analyses of samples collected in Bangui between 2017 and 2019 demonstrate the presence of one non-synonymous mutation (Y653N) and three synonymous mutations (frequency of 2.1%). This frequency of Pfkelch13 mutants is similar to those observed in neighbouring countries (Brazzaville Congo, certain regions of Cameroon), which showed frequencies of 1.57% and 2.9%, respectively [26, 27], confirming the absence of artemisinin resistance in Central Africa . The only previous data on Pfkelch13 polymorphism performed in 2014 in the CAR (K13 artemisinin resistance multicenter resistance, KARMA study), revealed a 4.5% prevalence of non-synonymous mutations . The difference or the fluctuation in the frequency of non-synonymous mutations in the CAR in 2014 and that of Bangui in 2017–2019 is likely due to the representativeness of the samples tested, but also to the fact that mutants appear at random and disappear probably because they do not have a selective advantage compared with wild strains . The previously validated or candidate resistance-associated mutations were not detected in our study. Interestingly, the A578S codon, common in Africa, was not observed either. The sole non-synonymous mutation detected here has never been reported in other African countries. It remains to be seen if this mutation is associated with artemisinin resistance. It is possible that local mutant strains resistant to artemisinin can emerge in addition to the risk of spread of the resistant Southeast Asian strain, as observed previously with strains resistant to chloroquine or sulfadoxine–pyrimethamine.
The dynamics of resistance and its emergence are likely to be complex, particularly due to the interactions specific to each world (sub) region . Similar studies in Cameroon and Nigeria have shown large differences according to region and time period, which makes it difficult to determine the spatio-temporal dynamics on polymorphism [6, 30]. Some countries in this African subregion (Congo, Gabon, and DRC) have revealed the presence of parasites polymorphic for the A578S allele, others have detected novel alleles and still others, particularly in Benin, have not detected any polymorphism [31, 32].
If the Y653N mutation increases in frequency in subsequent studies, it will then be necessary to assess its in vivo impact on the parasite clearance half-life in patients treated with ACT and its susceptibility in vitro to confirm or disprove its association with artemisinin resistance.
This study demonstrates the presence of a strain carrying a non-synonymous mutation in Bangui. Neither non-synonymous mutations involved in previously demonstrated in vitro and in vivo resistance nor the candidate resistance mutations were detected. The A578S mutant, although not associated with resistance but frequently found in Africa, was not observed among the tested samples. The novel mutation detected here is yet another mutation to add to the limited list of non-synonymous mutations detected in Africa. Additional studies with samples collected throughout the CAR are needed to confirm this narrow polymorphism profile.
Availability of data and materials
The database of this study is available from the corresponding author upon request.
Artemisinin-based combination therapy
Central African Republic
Institut Pasteur in Bangui
World Health Organization
Paediatric Complex in Bangui
- Pfkelch13 :
Plasmodium falciparum kelch 13 gene
Ring-stage Survival Assay
Cibulskis ER, Alonso P, Aponte J, Aregawi M, Barrette A, Bergeron L. Malaria: global progress 2000–2015 and future challenges. Infect Dis Poverty. 2016;5:61.
WHO. World malaria report. Geneva: World Health Organization; 2019.
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.
Hassett MR, Roepe PD. Origin and spread of evolving artemisinin-resistant Plasmodium falciparum malarial parasites in Southeast Asia. Am J Trop Med Hyg. 2019;101:1204–11.
Balikagala B, Mita T, Ikeda M, Sakurai M, Yatsushiro S, Takahashi N, et al. Absence of in vivo selection for K13 mutations after artemether–lumefantrine treatment in Uganda. Malar J. 2017;16:23.
Djaman JA, Olefongo D, Ako AB, Roman J, Ngane VF, Basco LK, Tahar R. Molecular epidemiology of malaria in Cameroon and Côte d’Ivoire. XXXI. Kelch 13 propeller sequences in Plasmodium falciparum isolates before and after implementation of artemisinin-based combination therapy. Am J Trop Med Hyg. 2017;97:222–4.
Manirakiza A, Njuimo SP, Le Faou A, Malvy D, Millet P. Availability of antimalarial drugs and evaluation of the attitude and practices for the treatment of uncomplicated malaria in Bangui, Central African Republic. J Trop Med. 2010;2010:510834.
Menard D, Djalle D, Manirakiza A, Yapou F, Siadoua V, Sana S, et al. Drug-resistant malaria in Bangui, Central African Republic: an in vitro assessment. Am J Trop Med Hyg. 2005;73:239–43.
Javelle E, Madamet M, Gaillard T, Velut G, Surcouf C, Michel R, et al. Delayed onset of Plasmodium falciparum malaria after doxycycline prophylaxis in a soldier returning from the Central African Republic. Antimicrob Agents Chemother. 2016;60:2592–3.
Nambei WS, Lango Yaya E, Pounguinza S, Achonduh O, Bogon A, Lengande R, et al. Efficacy and safety of antimalarial combinations for treatment of uncomplicated malaria in children in Bangui, Central African Republic. Med Sante Trop. 2013;23:313–9 (in French).
Ménard D, Khim N, Beghain J, Adegnika AA, Shafiul-Alam M, Amodu O. A worldwide map of Plasmodium falciparum K13-propeller polymorphisms. N Engl J Med. 2016;374:2453–64.
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.
WHO. Global technical strategy for malaria 2016–2030. Geneva: World Health Organization; 2016.
Walsh PS, Metzger DA, Higuchi R. Chelex 100 as a medium for simple extraction of DNA for PCR-based typing from forensic material. Biotechniques. 1991;10:506–13.
Singh B, Bobogare A, Cox-Singh J, Snounou G, Abdullah MS, Rahman HA. A genus- and species-specific nested polymerase chain reaction malaria detection assay for epidemiologic studies. Am J Trop Med Hyg. 1999;60:687–92.
Djallé D, Gody JC, Moyen JM, Tekpa G, Ipero J, Madji N, et al. Performance of Paracheck™-Pf, SD Bioline malaria Ag-Pf and SD Bioline malaria Ag-Pf/pan for diagnosis of falciparum malaria in the Central African Republic. BMC Infect Dis. 2014;14:109.
Bichara C, Flahaut P, Costa D, Bienvenu AL, Picot S, Gargala G. Cryptic Plasmodium ovale concurrent with mixed Plasmodium falciparum and Plasmodium malariae infection in two children from Central African Republic. Malar J. 2017;16:339.
Guerra RI, Ore M, Valdivia HO, Bishop DK, Ramos M, Mores CN, Campbell WR. A cluster of the first reported Plasmodium ovale spp. infections in Peru occuring among returning UN peace-keepers, a review of epidemiology, prevention and diagnostic challenges in non-endemic regions. Malar J. 2019;18:176.
WHO. Artemisinin resistance and artemisinin-based combination therapy efficacy: status report. Geneva: World Health Organization; 2018.
Yang C, Zhang H, Zhou R, Qian D, Liu Y, Zhao Y, et al. Polymorphisms of Plasmodium falciparum k13-propeller gene among migrant workers returning to Henan Province, China from Africa. BMC Infect Dis. 2017;17:560.
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.
Ouattara A, Kone A, Adams M, Fofana B, Maiga AW, Hampton S. Polymorphisms in the K13-propeller gene in artemisinin-susceptible Plasmodium falciparum parasites from Bougoula-Hameau and Bandiagara, Mali. Am J Trop Med Hyg. 2015;92:1202–6.
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.
Lu F, Culleton R, Zhang M, Ramaprasad A, von Seidlein L, Zhou H, et al. Emergence of indigenous artemisinin-resistant Plasmodium falciparum in Africa. N Engl J Med. 2017;376:991–3.
Mishra N, Prajapati SK, Kaitholia K, Bharti RS, Srivastava B, Phookan S, et al. Surveillance of artemisinin resistance in Plasmodium falciparum in India using the kelch13 molecular marker. Antimicrob Agents Chemother. 2015;59:2548–53.
Mayengue PI, Niama RF, Kouhounina Batsimba D, Malonga-Massanga A, Louzolo I, Loukabou Bongolo NC, et al. No polymorphisms in K13-propeller gene associated with artemisinin resistance in Plasmodium falciparum isolated from Brazzaville, Republic of Congo. BMC Infect Dis. 2018;18:538.
Eboumbou Moukoko CE, Huang F, Nsango SE, Kojom Foko LP, Ebong SB, Epee Eboumbou P, et al. K-13 propeller gene polymorphisms isolated between 2014 and 2017 from Cameroonian Plasmodium falciparum malaria patients. PLoS ONE. 2019;14:e0221895.
Day T, Huijben S, Read AF. Is selection relevant in the evolutionary emergence of drug resistance? Trends Microbiol. 2015;23:126–33.
Hastings IM. Complex dynamics and stability of resistance to antimalarial drugs. Parasitology. 2006;132:615–24.
Oboh MA, Ndiaye D, Antony HA, Badiane AS, Singh US, Ali NA, et al. Status of artemisinin resistance in malaria parasite Plasmodium falciparum from molecular analyses of the kelch13 gene in Southwestern Nigeria. Biomed Res Int. 2018;2018:2305062.
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.
Ogouyèmi-Hounto A, Damien G, Deme AB, Ndam NT, Assohou C, Tchonlin D, et al. Lack of artemisinin resistance in Plasmodium falciparum in northwest Benin after 10 years of use of artemisinin-based combination therapy. Parasite. 2016;23:28.
We gratefully acknowledge the participation of staff of Bangui Pediatric Complex (BPC).
This work was financed by the Institut Pasteur de Bangui and the Malaria Genetics and Resistance Unit at Institut Pasteur, Paris.
Ethics approval and consent to participate
The study protocol was accepted by malaria experts from the National Malaria Control Programme, in the absence of the institutional and national ethics committee that did not exist at the start of the study. The consent of each adult and the parents of child were obtained before inclusion.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Nzoumbou-Boko, R., Panté-Wockama, CB.G., Ngoagoni, C. et al. Molecular assessment of kelch13 non-synonymous mutations in Plasmodium falciparum isolates from Central African Republic (2017–2019). Malar J 19, 191 (2020). https://doi.org/10.1186/s12936-020-03264-y
- Plasmodium falciparum
- Antimalarial drug resistance
- Central African Republic