Skip to main content
Download PDF

Molecular surveillance of Kelch 13 polymorphisms in Plasmodium falciparum isolates from Kenya and Ethiopia

Malaria Journal volume 23, Article number: 36 (2024) Cite this article

Abstract

Background

Timely molecular surveillance of Plasmodium falciparum kelch 13 (k13) gene mutations is essential for monitoring the emergence and stemming the spread of artemisinin resistance. Widespread artemisinin resistance, as observed in Southeast Asia, would reverse significant gains that have been made against the malaria burden in Africa. The purpose of this study was to assess the prevalence of k13 polymorphisms in western Kenya and Ethiopia at sites representing varying transmission intensities between 2018 and 2022.

Methods

Dried blood spot samples collected through ongoing passive surveillance and malaria epidemiological studies, respectively, were investigated. The k13 gene was genotyped in P. falciparum isolates with high parasitaemia: 775 isolates from four sites in western Kenya (Homa Bay, Kakamega, Kisii, and Kombewa) and 319 isolates from five sites across Ethiopia (Arjo, Awash, Gambella, Dire Dawa, and Semera). DNA sequence variation and neutrality were analysed within each study site where mutant alleles were detected.

Results

Sixteen Kelch13 haplotypes were detected in this study. Prevalence of nonsynonymous k13 mutations was low in both western Kenya (25/783, 3.19%) and Ethiopia (5/319, 1.57%) across the study period. Two WHO-validated mutations were detected: A675V in three isolates from Kenya and R622I in four isolates from Ethiopia. Seventeen samples from Kenya carried synonymous mutations (2.17%). No synonymous mutations were detected in Ethiopia. Genetic variation analyses and tests of neutrality further suggest an excess of low frequency polymorphisms in each study site. Fu and Li’s F test statistic in Semera was 0.48 (P > 0.05), suggesting potential population selection of R622I, which appeared at a relatively high frequency (3/22, 13.04%).

Conclusions

This study presents an updated report on the low frequency of k13 mutations in western Kenya and Ethiopia. The WHO-validated R622I mutation, which has previously only been reported along the north-west border of Ethiopia, appeared in four isolates collected from eastern Ethiopia. The rapid expansion of R622I across Ethiopia signals the need for enhanced monitoring of the spread of drug-resistant P. falciparum parasites in East Africa. Although ACT remains currently efficacious in the study areas, continued surveillance is necessary to detect early indicators of artemisinin partial resistance.

Background

Efficacious antimalarial drugs are mainstays of malaria control and elimination programs. Widespread resistance to chloroquine (CQ) and sulfadoxine-pyrimethamine (SP) led to the adoption of artemisinin-based combination therapy (ACT) as the first- and second-line treatments for uncomplicated Plasmodium falciparum malaria and chloroquine-resistant Plasmodium vivax malaria in most malaria-endemic regions at the turn of the twenty-first century. Within less than a decade, resistance to artemisinins and partner drugs emerged in Southeast Asia. Historically, resistance to CQ and SP arose in Southeast Asia and spread to the African continent, leading to increases in malaria-related deaths. Africa disproportionately shoulders the global burden of malaria [1]. The spread of ACT resistance to Africa would result in devastating human health and economic costs, as there is currently no new class of antimalarials that is suitable for immediate, large-scale implementation as an alternative to ACT [2].

Genetic markers that are associated with drug resistance serve as valuable tools for assessing the geographic origins of resistance and tracking the spread of resistance. The identification of nonsynonymous single nucleotide polymorphisms (SNPs) in a gene encoding the Kelch propeller domain on P. falciparum chromosome 13 (PfKelch13, or Pfk13) as molecular markers of artemisinin partial resistance has greatly facilitated global surveillance of artemisinin resistance [3]. Over 200 nonsynonymous SNPs in the k13 gene have been reported worldwide. However, k13 mutations have differential effects on clearance phenotypes, and not all nonsynonymous SNPs indicate artemisinin partial resistance. Extensive clinical and laboratory evidence are required for SNPs to be validated as molecular markers of artemisinin partial resistance [4]. Each newly discovered k13 mutation must exhibit both significant association with delayed parasite clearance and > 1% survival in ring-stage survival assays to be classified as a validated k13 marker by the World Health Organization (i.e., WHO-validated). A mutation that meets only one of these two criteria is considered a candidate resistance marker, while a mutation that meets neither of these criteria is considered neutral [4, 5].

Nonsynonymous k13 mutations are currently relatively rare in Africa, and ACT remains efficacious across the continent [6]. However, partial artemisinin resistance and signs of resistance to ACT partner drugs have been identified, mainly in East Africa [7]. Validated k13 mutations have been detected in Uganda, Rwanda, and Tanzania, raising alarm in neighbouring East African countries [8,9,10]. In both Kenya and Ethiopia, CQ was replaced by SP for the treatment of uncomplicated P. falciparum malaria in 1998, followed by yet another policy change from SP to ACT in 2004. The switch to ACT was widely implemented in Kenya in 2006 when artemether-lumefantrine (AL) was made available in government hospitals, while rollout of AL across Ethiopia began in 2005 [11, 12]. Studies that have been conducted to investigate changes in k13 mutant allele frequency and diversity in Kenya [13,14,15,16,17] and Ethiopia [18,19,20,21] have largely focused on samples collected within the first ten years after the introduction of ACT.

The WHO recommends testing the efficacy of first-line antimalarial treatments at least once every 24 months at sentinel sites to monitor the emergence, importation, and spread of resistant alleles and to ensure effective case management. To this end, dried blood spot samples collected through ongoing passive surveillance and malaria epidemiological studies were leveraged to assess the prevalence of k13 polymorphisms in western Kenya and Ethiopia at sites representing varying transmission intensities between 2018 and 2022.

Methods

Study sites

Samples were collected from four study sites of varying transmission levels in western Kenya (Fig. 1). Kombewa (00°07′N, 34°30′E, 1,260 m above sea level [m a.s.l.]), in the lowland area of the Lake Victoria basin in Kisumu County, is a holoendemic site. Kakamega (00°10′N, 34°45′E, 1,500 m a.s.l.), in Kakamega County, represents the highland mesoendemic area. Kisii (00°35′S, 34°48′E, 1700 m a.s.l.) is a hypoendemic site in Kisii County. Homa Bay (0°31′S, 34°27′E, 1,190 m a.s.l), in southern Homa Bay County, represents the Lake mesoendemic transmission level. Parasite prevalence in Homa Bay County was historically greater than 20%, but indoor residual spraying campaigns conducted by the Kenyan Ministry of Health and the U.S. President’s Malaria Initiative since 2018 have reduced transmission significantly [22, 23]. At all sites, P. falciparum is the primary malaria parasite species. Western Kenya experiences a bimodal pattern of rainfall. January-March is the hot and dry season, April-June is the long rainy season, July-September is the cool and dry season, and October-December is the short rainy season.

Fig. 1
figure 1

Map of study sites in Kenya (Homa Bay, Kakamega, Kisii, and Kombewa) and Ethiopia (Arjo, Awash, Dire Dawa, Gambella, and Semera)

Full size image

Samples from five study sites across Ethiopia (Fig. 1), representing different transmission risk strata, as defined based on historical annual parasite incidence by Ethiopia’s National Malaria Control Strategy, were also analyzed [24]. Dire Dawa (9°36′N, 41°52′E, 1,280 m a.s.l.) in eastern Ethiopia and Arjo (8°42′N, 36°26′E, 1,350 m a.s.l.) in western Ethiopia are in the “very low risk” stratum. Semera (11°47′N, 41°0′E, 430 m a.s.l.) in northeastern Ethiopia and Hawassa (7°3′N, 38°28′E, 1700 m a.s.l.) in the central Ethiopian Rift Valley are in the “moderate risk” stratum. Gambella (7°53′N, 34°30′E, 500 m a.s.l.) in remote western Ethiopia is in the “high risk” stratum. P. falciparum and P. vivax are co-endemic in these sites. Ethiopia has three rainy seasons. The primary rainy season occurs from June-September, and a short rainy season occurs from February-May. Southern Ethiopia experiences a drier and colder third rainy season from October-December [25].

Study design

The primary purpose of this study was to assess the overall prevalence of k13 polymorphisms in study sites established under the National Institutes of Health’s sub-Saharan Africa International Centers of Excellence for Malaria Research (ICEMR) programme. Thus, dried blood spot samples (DBS) collected through ICEMR surveillance and research activities between 2018 and 2022 were investigated. Samples from four ICEMR study sites in Kenya (Homa Bay, Kakamega, Kisii, and Kombewa) and five ICEMR study sites in Ethiopia (Arjo, Awash, Dire Dawa, Gambella, Semera) were analysed (Table 1). Samples were identified through convenience sampling, therefore the number of samples per site and collection year was not balanced. Briefly, the sources of samples used in this study included passive case detection (PCD) at health centers, clinics, and hospitals; seasonal, cross-sectional mass blood surveys (MBS) in randomly selected households; schoolchildren dynamics studies, which seek to describe transmission dynamics among schoolchildren aged 5-18; and a three-year longitudinal cohort study with monthly follow-up. PCD captured cases of malaria that were confirmed via light microscopy or rapid diagnostic tests among community members who sought treatment at health facilities, while the remaining studies captured asymptomatic cases at households or schools within communities. For each study, DBS were collected via finger prick from all community members who consented to participating. DBS were stored with silica gel desiccant at 4 ºC until usage.

Table 1 ICEMR samples collected between 2018 and 2022 used for Pfk13 polymorphism investigation
Full size table

Parasite genotyping

DNA was extracted from DBS using the saponin/ Chelex method [26]. Plasmodium falciparum infections were identified by varATS real-time quantitative PCR (qPCR) on a QuantStudio 3 Real-Time PCR System (Thermo Fisher Scientific, Carlsbad, CA) [27]. Isolates with cycle threshold (Ct) values ≤ 32 were selected to ensure sufficiently high parasitaemia for k13 genotyping. In samples that met this Ct threshold, the k13 propeller domain was amplified by nested PCR, using primers designed by Ariey et al. [3] and PerfeCTa qPCR ToughMix (QuantaBio, Beverly, MA). Cycling conditions and specifications of reaction mixtures are presented in Additional file 1: Table S1. Gel electrophoresis was performed at 120 V (2% agarose gel with ethidium bromide) on 5 µL of each nested PCR product, along with a 100 bp DNA ladder (New England Biolabs, Ipswich, MA), and gels were visualised under UV light to confirm amplification and correct amplicon size. For each successfully amplified sample, enzymatic cleanup on the remaining PCR product was performed using exonuclease I and shrimp alkaline phosphatase (Thermo Fisher Scientific). Purified PCR amplicons were screened for polymorphisms by Sanger sequencing in the forward direction (GENEWIZ, Inc., La Jolla, CA). For all samples with polymorphisms in the initial forward direction screening, fresh aliquots of DNA were used to prepare and clean PCR amplicons, following the steps described above. The purified PCR products were then sequenced in both directions to confirm polymorphisms.

Data analysis

CodonCode Aligner 10.0.3 (CodonCode Corporation, Centerville, MA) was used to trim low-quality ends, assemble sequences, generate contigs, and align sequences. Sequences were exported and analysed using BioEdit 7.2 sequence alignment editor software [28]. Nucleotide and amino acid sequences were compared with the P. falciparum 3D7 chromosome 13 reference strain (GenBank Accession Number, CP017003.1). DNA Sequence Polymorphism (DnaSP) version 6 software [29] was used to first infer phased haplotypes for samples where heterozygous peaks were observed in sequencing chromatograms (e.g., presumed infection with two clones), and then to calculate population-level genetic variation indices, including the proportion of polymorphic loci, number of haplotypes, haplotype diversity, and nucleotide diversity. DnaSP was also used to perform Tajima’s and Fu and Li’s neutrality tests for all sites with mutant alleles (i.e., Homa Bay, Kakamega, and Kombewa in Kenya; Arjo, Dire Dawa, and Semera in Ethiopia). Haplotype sequences were deposited in NCBI GenBank under accession numbers OR571911-OR571926.

Results

A total of 1,094 DBS between 2018 and 2022, comprising 775 DBS from Kenya and 319 DBS from Ethiopia, were successfully sequenced. Eight samples from Kenya showed multiclonal infections, as evidenced by double peaks in single nucleotide positions in the sequencing chromatograms. Thus, an additional eight sequences were obtained by phasing the isolates with heterozygote peaks, yielding 783 sequences from the Kenyan isolates. Sixteen unique haplotypes were identified from a total of 1,102 sequences analysed in this study. k13 polymorphisms were detected in 42 samples from Kenya and five samples from Ethiopia (Table 2). Each mutant sample carried only one SNP in the k13 propeller region, and only one haplotype (Hap_2, carrying the A582V mutation) was shared between the Kenyan and Ethiopian samples (Fig. 2). Among these samples, limited occurrences of two mutations that were added to the list of WHO-validated markers of artemisinin partial resistance in the 2022 World Malaria Report were observed [6]: A675V in Kenya and R622I in Ethiopia.

Table 2 Frequency and distribution of Pfk13 propeller domain haplotypes in samples from Kenya and Ethiopia, 2018-2022
Full size table
Fig. 2
figure 2

Haplotype network of P. falciparum isolates harbouring k13 polymorphisms in study sites in Kenya and Ethiopia. The size of each circle corresponds to the haplotype prevalence. Each color represents the site where the haplotype was found; pie charts comprising multiple colors indicate the presence of the same haplotype in multiple study sites. Hap_1 represents the reference k13 sequence (Pf3D7)

Full size image

The prevalence of nonsynonymous k13 mutations in Kenya was 3.25%. The most common mutations overall among the Kenyan samples were A578S (9/783, 1.15%) and A582V (10/783, 1.28%). A578S appeared in Kakamega in 2019, 2021, and 2022 at increasing frequencies (Additional file 1: Table S3). Three samples collected in 2022 in Homa Bay also harboured the A578S mutation. A582V was observed at low frequencies in Homa Bay in 2018 (1.45%) and again in 2022 (0.64%). A675V occurred only in Kakamega in 2021, at a frequency of 25.00% (3/12). The nonsynonymous mutations A569S, I646L, and P667S were detected in only one sample each. Frequencies of synonymous k13 polymorphisms in Kenya were low, ranging from 0.13-0.50%. Haplotype diversity was low in all four Kenyan study sites, with Hd ranging from 0 in Kisii, where no polymorphisms were detected, to 0.302 in Kakamega, where eight haplotypes were found. Nucleotide diversity was similarly low (πi < 0.00015). Tajima’s D and Fu and Li’s F were negative for each population tested, congruous with the observation of an excess of rare alleles compared to the expected diversity under the assumptions of neutral selection and constant population size. However, both test statistics were only significant at a p-value < 0.05 for Homa Bay and Kombewa (Table 3).

Table 3 Molecular diversity of the Pfk13 propeller gene in Kenya and Ethiopia, 2018-2022
Full size table

The prevalence of nonsynonymous k13 mutations in Ethiopia was 1.57%. One sample collected in Arjo in 2020 harboured the A582V mutation. Four occurrences of R622I were found among the 2022 samples: one sample from Dire Dawa and three samples from Semera. No synonymous k13 mutations were detected in the 319 Ethiopian samples tested (Table 2). Site-level analyses revealed low haplotype diversity (Hd: 0-0.247) and nucleotide diversity (πi < 0.00011). Negative Tajima’s D and Fu and Li’s F statistic values indicate an excess of singletons in Arjo and Dire Dawa (P > 0.05 for all neutrality tests). In Semera, Tajima’s D value was -0.17 (P > 0.05) and Fu and Li’s F value was 0.48 (P > 0.05), suggesting possible balancing selection of the R622I mutation in the study area.

Discussion

In this study, the prevalence of synonymous and nonsynonymous mutations in 775 P. falciparum isolates collected from four sites in western Kenya and 319 P. falciparum isolates collected from five sites across Ethiopia between 2018 and 2022 was investigated. Low frequency of k13 polymorphisms and few occurrences of shared alleles between geographically distinct populations were observed. These findings align with previous studies conducted in western Kenya [15, 16, 20, 30], Ethiopia [20, 21, 31], and across Africa [5, 20, 32, 33] following the implementation of ACT.

Seven unique nonsynonymous mutations were found in this study. A582V was the most common mutation among the Kenyan samples and is the only nonsynonymous mutation that was observed in both Kenya and Ethiopia. The mutation has been reported once among parasites collected across Kilifi County, Kenya between 1994 and 2018 [34]. A578S was the second most common mutation among the Kenyan samples in this study and is the most common k13 polymorphism found across Africa [5]. An increase in the relative frequency of A578S mutants in Kakamega between 2019 and 2022 was observed, in line with a previous study conducted in Kombewa and Kakamega [30]. However, A578S is not associated with clinical or in vitro artemisinin partial resistance [5]. A675V, a recently validated marker of artemisinin partial resistance [6], appeared in three out of the 12 samples tested from Kakamega in 2021; the mutation was not detected at any study site in 2022. A569S, which appeared in only one 2019 sample, was reported in Kenya in 2013 in one sample from Mbita district, Homa Bay County [13] and one sample from Kisumu County [20]. I646L has been previously reported in Asia but has yet to be evaluated by the WHO or the Worldwide Antimalarial Resistance Network [35]. Another nonsynonymous mutation at the same codon, I646T, was detected in Sierra Leone in 2016, but the mutation similarly has not yet been evaluated for resistance [36]. In Rwanda, Straimer et al. reported parasite clearance time > 5 h. in a patient infected with P. falciparum parasites carrying the P667S mutation [37]. However, further investigation is required to confirm the potential impact of P667S on ring-stage survival in vitro.

R622I, which appeared in four PCD samples collected in 2022 in Ethiopia, is another recently validated marker of partial artemisinin resistance [38]. The mutation was first detected in samples from Northwest Ethiopia and increased in prevalence from 2.4% in 2014 to 9.5% in 2017-2018 [39]. Within Ethiopia, R622I has also been reported in samples from districts further north, along the border with Eritrea, and further south, along the border with South Sudan [40]. R622I has also been found in samples from Somalia, Sudan, Mozambique, Zambia, and Eritrea [38, 41]. In Eritrea, which shares a border with northern Ethiopia, the cumulative frequency of R622I in 729 samples collected between 2015 and 2020 was approximately 14% [6]. This study is the first to report cases of R622I in Semera (northeast Ethiopia), further underscoring the pressing need for monitoring clonal expansion of R622I across the country.

Samples were aggregated by study site for molecular diversity analyses in this study due to the limited number of k13 polymorphisms found when samples were stratified by year and infection type. The imbalances in the number of samples analysed by site, year, and infection type reflect not only variations in transmission intensity by site but also differences in the ongoing activities among study sites. Most of the Kenyan samples investigated arose from asymptomatic infections (88.88%, 711/800), whereas most of the Ethiopian samples came from symptomatic infections detected at health facilities (92.16%, 294/319). In 2020, many health facilities were converted to treat patients afflicted with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Thus, the symptomatic samples used in this study were predominantly collected from health facilities in the years preceding the SARS-CoV-2 pandemic. On the other hand, passive surveillance remained largely uninterrupted in the study sites in Ethiopia, as evidenced by the symptomatic (i.e., PCD) samples available for use in this study, whereas community-based malaria surveys (i.e., asymptomatic samples) were pared down due to limited field team capacity and resources during the SARS-CoV-2 pandemic. No significant differences in the frequency or diversity of k13 polymorphisms were observed between symptomatic and asymptomatic infections (Additional file 1: Tables S3, 4) when samples derived from both infection types were available for analysis, suggesting that the nature of infections did not clearly bias overall findings. However, the sample sizes available for comparison were small. Additionally, other studies have found that the extent to which symptomatic cases represent the asymptomatic parasite reservoir may vary with changes in transmission intensity [42, 43]. Therefore, spatiotemporally matched screening of both asymptomatic and symptomatic infections is necessary to evaluate potential sampling bias in molecular surveillance of artemisinin resistance.

The low number of mutations precluded statistical comparisons between years or study sites. Nevertheless, one key strength of this study was the optimization of available samples collected through various study sites and designs, which helped reduce confounding due to factors such as age and transmission intensity among the samples from Kenya. Several studies reporting low frequency of k13 polymorphisms have focused exclusively on school-aged children [30, 44,45,46]. Despite constituting important drivers of malaria transmission, schoolchildren are frequently overlooked in malaria control programmes [47,48,49,50]. DBS were collected from asymptomatic school children aged 5-18 in Kakamega, Kisii, and Kombewa and from participants spanning the age spectrum in all studies conducted in Homa Bay (Additional file 1: Table S2). Interestingly, in nearly all studies conducted in Homa Bay, infections with high parasitemia (Ct ≤ 32) occurred most frequently in children 15 years old or younger. Furthermore, all 42 Kenyan samples carrying k13 mutations (synonymous and nonsynonymous) occurred in children who were 4-15 years old. No associations between transmission intensity and nonsynonymous k13 mutant allele frequency or diversity were observed. Taken together, these observations suggest that schoolchildren may be an appropriate sentinel group for artemisinin partial resistance surveillance in western Kenya, particularly when resources are limited. On the other hand, the age distribution among the Ethiopian samples tested for k13 mutations more closely matched the age distribution of all enrolled participants within each study from which the samples were sourced. A higher frequency of R622I was observed in Semera compared to Dire Dawa (“moderate” vs. “very low” transmission sites, respectively). However, no k13 mutations were found in Awash, which is similarly located in the “moderate risk” stratum. Surveillance strategies in Ethiopia might benefit from consideration of additional factors, aside from age and transmission intensity, in defining important sentinel groups and risk areas.

The variations in sample sizes also reflect differences in parasite density distributions by study sites and years. The number of samples tested for k13 mutations in each study and year does not reflect the total number of samples collected but rather the number of samples with Ct ≤ 32. As this study was restricted to high density infections, rare mutant alleles present in low density infections may have been missed. Low density infections yield low concentrations of parasite DNA, rendering identification of minor alleles challenging to replicate. Additionally, low density infections often resolve without treatment and thus may have limited clinical relevance [51]. However, these infections may still contribute to malaria transmission and should be explored using more sensitive methods in future studies. Compared to next generation sequencing methods, Sanger sequencing, as used in this study, is less sensitive for detecting minor alleles in polyclonal infections [52]. Although the diversity of the parasite population is likely underestimated in this study, Sanger sequencing remains a highly accurate and rapid method for variant screening, and the use of more sensitive deep sequencing methods would sustain the overall trends reported here (no observations of clonal spread).

Conclusion

An updated report on the low frequency of k13 mutations in several study sites of western Kenya and across Ethiopia is presented in this study, in line with many other studies that have concluded that k13 mutations remain rare in Africa. The k13 validated marker A675V was restricted to three samples from asymptomatic schoolchildren in Kakamega, Kenya in 2019. In Ethiopia, four occurrences of the k13 validated marker R622I were detected among symptomatic 2022 PCD cases from Dire Dawa and Semera, constituting one of the first reports of R622I from eastern Ethiopia. While ACT remains largely efficacious in these study sites in western Kenya and Ethiopia, enhanced molecular surveillance of P. falciparum resistance to artemisinin is necessary to monitor the spread of validated markers of resistance and the emergence of additional mutants of concern.

Availability of data and materials

The dataset generated in this study is available from the corresponding author upon reasonable request.

Abbreviations

Pfk13:

Kelch propeller domain on Plasmodium falciparum chromosome 13

DBS:

Dried blood spot

SNP:

Single nucleotide polymorphism

CQ:

Chloroquine

SP:

Sulfadoxine-pyrimethamine

AL:

Artemether-lumefantrine

ACT:

Artemisinin-based combination therapy

WHO:

World Health Organization

Ct:

Cycle threshold value

m a.s.l.:

Metres above sea level

ICEMR:

International Centers of Excellence for Malaria Research

PCD:

Passive case detection

MBS:

Mass blood surveys

varATS qPCR:

Var gene acidic terminal sequence-based quantitative polymerase chain reaction

References

  1. Phyo AP, Nosten F, Phyo AP, Nosten F. The Artemisinin Resistance in Southeast Asia: An Imminent Global Threat to Malaria Elimination. In: Manguin S, Dev V (eds.); Towards Malaria Elimination - A Leap Forward. IntechOpen; 2018 [cited 2023 Jul 12]. https://www.intechopen.com/chapters/60988

  2. Lubell Y, Dondorp A, Guérin PJ, Drake T, Meek S, Ashley E, et al. Artemisinin resistance—modelling the potential human and economic costs. Malar J. 2014;13:452.

    Article  PubMed  PubMed Central  Google Scholar 

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

    Article  PubMed  Google Scholar 

  4. Global Malaria Programme. Artemisinin and artemisinin-based combination therapy resistance, April 2017. World Health Organization; 2017 Apr https://iris.who.int/bitstream/handle/10665/255213/WHO-HTM-GMP-2017.9-eng.pdf?sequence=1.

  5. 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 

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

    Google Scholar 

  7. Kayiba NK, Yobi DM, Tshibangu-Kabamba E, Tuan VP, Yamaoka Y, Devleesschauwer B, et al. Spatial and molecular mapping of Pfkelch13 gene polymorphism in Africa in the era of emerging Plasmodium falciparum resistance to artemisinin: a systematic review. Lancet Infect Dis. 2021;21:e82-92.

    Article  PubMed  CAS  Google Scholar 

  8. Balikagala B, Fukuda N, Ikeda M, Katuro OT, Tachibana SI, Yamauchi M, et al. Evidence of artemisinin-resistant malaria in Africa. N Engl J Med. 2021;385:1163–71.

    Article  PubMed  CAS  Google Scholar 

  9. Uwimana A, Umulisa N, Venkatesan M, Svigel SS, Zhou Z, Munyaneza T, et al. Association of Plasmodium falciparum kelch13 R561H genotypes with delayed parasite clearance in Rwanda: an open-label, single-arm, multicentre, therapeutic efficacy study. Lancet Infect Dis. 2021;21:1120–8.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  10. Bwire GM, Ngasala B, Mikomangwa WP, Kilonzi M, Kamuhabwa AAR. Detection of mutations associated with artemisinin resistance at k13-propeller gene and a near complete return of chloroquine susceptible falciparum malaria in Southeast of Tanzania. Sci Rep. 2020;10:3500.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  11. National strategic plan for malaria prevention control and elimination in Ethiopia 2011–2015. Addis Ababa, Ethiopia; 2010 [cited 2023 Jul 13]. https://www.nationalplanningcycles.org/sites/default/files/country_docs/Ethiopia/ethiopia_malaria_national_strategic_plan_2011-2015_130810.pdf

  12. Amin AA, Zurovac D, Kangwana BB, Greenfield J, Otieno DN, Akhwale WS, et al. The challenges of changing national malaria drug policy to artemisinin-based combinations in Kenya. Malar J. 2007;6:72.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Isozumi R, Uemura H, Kimata I, Ichinose Y, Logedi J, Omar AH, et al. Novel mutations in K13 propeller gene of artemisinin-resistant Plasmodium falciparum. Emerg Infect Dis. 2015;21:490–2.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  14. Wamae K, Okanda D, Ndwiga L, Osoti V, Kimenyi KM, Abdi AI, et al. No Evidence of Plasmodium falciparum k13 artemisinin resistance-conferring mutations over a 24-year analysis in Coastal Kenya but a near complete reversion to chloroquine-sensitive parasites. Antimicrob Agents Chemother. 2019;63:e01067-e1119.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  15. de Laurent ZR, Chebon LJ, Ingasia LA, Akala HM, Andagalu B, Ochola-Oyier LI, et al. Polymorphisms in the K13 gene in Plasmodium falciparum from different malaria transmission areas of Kenya. Am J Trop Med Hyg. 2018;98:1360–6.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Maniga JN, Samuel M, John O, Rael M, Muchiri JN, Bwogo P, et al. Novel Plasmodium falciparum k13 gene polymorphisms from Kisii County, Kenya during an era of artemisinin-based combination therapy deployment. Malar J. 2023;22:87.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  17. Lucchi NW, Komino F, Okoth SA, Goldman I, Onyona P, Wiegand RE, et al. In vitro and molecular surveillance for antimalarial drug resistance in Plasmodium falciparum parasites in Western Kenya reveals sustained artemisinin sensitivity and increased chloroquine sensitivity. Antimicrob Agents Chemother. 2015;59:7540–7.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  18. Bayih AG, Getnet G, Alemu A, Getie S, Mohon AN, Pillai DR. A unique plasmodium falciparum Kelch 13 gene mutation in Northwest Ethiopia. Am J Trop Med Hyg. 2016;94:132–5.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  19. Lo E, Hemming-Schroeder E, Yewhalaw D, Nguyen J, Kebede E, Zemene E, et al. Transmission dynamics of co-endemic Plasmodium vivax and P. falciparum in Ethiopia and prevalence of antimalarial resistant genotypes. PLoS Negl Trop Dis. 2017;11:0005806.

    Article  Google Scholar 

  20. 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.

    PubMed  CAS  Google Scholar 

  21. Brazeau NF, Assefa A, Mohammed H, Seme H, Tsadik AG, Parr JB, et al. Pooled deep sequencing of drug resistance loci from Plasmodium falciparum parasites across Ethiopia. Am J Trop Med Hyg. 2019;101:1139–43.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  22. Omondi CJ, Otambo WO, Odongo D, Ochwedo KO, Otieno A, Onyango SA, et al. Asymptomatic and submicroscopic Plasmodium infections in an area before and during integrated vector control in Homa Bay, western Kenya. Malar J. 2022;21:272.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  23. National Malaria Control Programme. Kenya Malaria Strategy 2019–2023. Nairobi, Kenya: Ministry of Health; [cited 2020 Oct 23]. http://fountainafrica.org/wp-content/uploads/2020/01/Kenya-Malaria-Strategy-2019-2023.pdf

  24. Federal Ministry of Health. Ethiopia Malaria Elimination Strategic Plan: 2021–2025. Addis Ababa, Ethiopia: 2020.

  25. The World Bank Group. Climate risk country profile: Ethiopia [Internet]. Washington DC, USA: 2021 [cited 2023 Jul 17]. https://climateknowledgeportal.worldbank.org/sites/default/files/2021-05/15463A-WB_Ethiopia%20Country%20Profile-WEB.pdf

  26. Bereczky S, Mårtensson A, Gil JP, Färnert A. Rapid DNA extraction from archive blood spots on filter paper for genotyping of Plasmodium falciparum. Am J Trop Med Hyg. 2005;72:249–51.

    Article  PubMed  Google Scholar 

  27. Zhong D, Hemming-Schroeder E, Wang X, Kibret S, Zhou G, Atieli H, et al. Extensive new Anopheles cryptic species involved in human malaria transmission in western Kenya. Sci Rep. 2020;10:16139.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  28. Hall TA. BioEdit : a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp Ser. 1999;41:95–8.

    CAS  Google Scholar 

  29. Rozas J, Ferrer-Mata A, Sánchez-DelBarrio JC, Guirao-Rico S, Librado P, Ramos-Onsins SE, et al. DnaSP 6: DNA Sequence polymorphism analysis of large data sets. Mol Biol Evol. 2017;34:3299–302.

    Article  PubMed  CAS  Google Scholar 

  30. Hemming-Schroeder E, Umukoro E, Lo E, Fung B, Tomás-Domingo P, Zhou G, et al. Impacts of antimalarial drugs on Plasmodium falciparum drug resistance markers, Western Kenya, 2003–2015. Am J Trop Med Hyg. 2018;98:692–9.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  31. Heuchert A, Abduselam N, Zeynudin A, Eshetu T, Löscher T, Wieser A, et al. Molecular markers of anti-malarial drug resistance in southwest Ethiopia over time: regional surveillance from 2006 to 2013. Malar J. 2015;14:208.

    Article  PubMed  PubMed Central  Google Scholar 

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

  33. Ndwiga L, Kimenyi KM, Wamae K, Osoti V, Akinyi M, Omedo I, et al. A review of the frequencies of Plasmodium falciparum Kelch 13 artemisinin resistance mutations in Africa. Int J Parasitol Drugs Drug Resist. 2021;16:155–61.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  34. Omedo I, Bartilol B, Kimani D, Gonçalves S, Drury E, Rono MK, et al. Spatio-temporal distribution of antimalarial drug resistant gene mutations in a Plasmodium falciparum parasite population from Kilifi, Kenya: a 25-year retrospective study. Wellcome Open Res. 2022;7:45.

    Article  Google Scholar 

  35. Kagoro FM, Barnes KI, Marsh K, Ekapirat N, Mercado CEG, Sinha I, et al. Mapping genetic markers of artemisinin resistance in Plasmodium falciparum malaria in Asia: a systematic review and spatiotemporal analysis. Lancet Microbe. 2022;3:e184–92.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  36. Smith SJ, Kamara ARY, Sahr F, Samai M, Swaray AS, Menard D, et al. Efficacy of artemisinin-based combination therapies and prevalence of molecular markers associated with artemisinin, piperaquine and sulfadoxine-pyrimethamine resistance in Sierra Leone. Acta Trop. 2018;185:363–70.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  37. Straimer J, Gandhi P, Renner KC, Schmitt EK. High prevalence of Plasmodium falciparum K13 mutations in Rwanda is associated with slow parasite clearance after treatment with artemether-lumefantrine. J Infect Dis. 2021;225:1411–4.

    Article  PubMed Central  Google Scholar 

  38. WHO. Report on antimalarial drug efficacy, resistance and response: 10 years of surveillance (2010–2019). Geneva: World Health Organization; 2020.

    Google Scholar 

  39. Alemayehu AA, Castaneda-Mogollon D, Tesfa H, Getie S, Mohon AN, Balasingam N, et al. Expansion of the Plasmodium falciparum Kelch 13 R622I mutation in Northwest Ethiopia [Internet]. Res Square. 2021. https://www.researchsquare.com/article/rs-171038/v1

  40. Fola AA, Feleke SM, Mohammed H, Brhane BG, Hennelly CM, Assefa A, et al. Plasmodium falciparum resistant to artemisinin and diagnostics have emerged in Ethiopia. Nat Microbiol. 2023;8:1911–9.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  41. Wang X, Ruan W, Zhou S, Huang F, Lu Q, Feng X, et al. Molecular surveillance of Pfcrt and k13 propeller polymorphisms of imported Plasmodium falciparum cases to Zhejiang Province, China between 2016 and 2018. Malar J. 2020;19:59.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  42. Pava Z, Handayuni I, Trianty L, Utami RAS, Tirta YK, Puspitasari AM, et al. Passively versus actively detected malaria: similar genetic diversity but different complexity of infection. Am J Trop Med Hyg. 2017;97:1788–96.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  43. Falq G, Van Den Bergh R, De Smet M, Etienne W, Nguon C, Rekol H, et al. Assessing the asymptomatic reservoir and dihydroartemisinin–piperaquine effectiveness in a low transmission setting threatened by artemisinin resistant Plasmodium falciparum. Malar J. 2016;15:446.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Osoti V, Akinyi M, Wamae K, Kimenyi KM, de Laurent Z, Ndwiga L, et al. Targeted amplicon deep sequencing for monitoring antimalarial resistance markers in Western Kenya. Antimicrob Agents Chemother. 2022;66:e01945-e2021.

    Article  PubMed  PubMed Central  Google Scholar 

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

  46. Musyoka KB, Kiiru JN, Aluvaala E, Omondi P, Chege WK, Judah T, et al. Prevalence of mutations in Plasmodium falciparum genes associated with resistance to different antimalarial drugs in Nyando, Kisumu County in Kenya. Infect Genet Evol. 2020;78: 104121.

    Article  PubMed  CAS  Google Scholar 

  47. Coalson JE, Walldorf JA, Cohee LM, Ismail MD, Mathanga D, Cordy RJ, et al. High prevalence of Plasmodium falciparum gametocyte infections in school-age children using molecular detection: patterns and predictors of risk from a cross-sectional study in southern Malawi. Malar J. 2016;15:527.

    Article  PubMed  PubMed Central  Google Scholar 

  48. Bashir IM, Nyakoe N, van der Sande M. Targeting remaining pockets of malaria transmission in Kenya to hasten progress towards national elimination goals: an assessment of prevalence and risk factors in children from the Lake endemic region. Malar J. 2019;18:233.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Rek J, Blanken SL, Okoth J, Ayo D, Onyige I, Musasizi E, et al. Asymptomatic school-aged children are important drivers of malaria transmission in a high endemicity setting in Uganda. J Infect Dis. 2022;226:708–13.

    Article  PubMed  PubMed Central  Google Scholar 

  50. Touray AO, Mobegi VA, Wamunyokoli F, Herren JK. Diversity and multiplicity of P. falciparum infections among asymptomatic school children in Mbita. Western Kenya Sci Rep. 2020;10:5924.

    PubMed  CAS  Google Scholar 

  51. Kanai M, Yeo T, Asua V, Rosenthal PJ, Fidock DA, Mok S. Comparative analysis of Plasmodium falciparum genotyping via SNP detection, microsatellite profiling, and whole-genome sequencing. Antimicrob Agents Chemother. 2022;66:e01163-e1221.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  52. Hartley MA, Hofmann N, Keitel K, Kagoro F, Moniz CA, Mlaganile T, et al. Clinical relevance of low-density Plasmodium falciparum parasitemia in untreated febrile children: a cohort study. PLoS Med. 2020;17:e1003318.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

Download references

Acknowledgements

We sincerely thank the study participants, communities, and hospital staff who participated in the studies at all study sites. This work was made possible by the sub-Saharan Africa ICEMR field teams, study coordinators, and Yan lab members.

Funding

This study was supported by grants U19AI129326 and D43TW001505 from the National Institutes of Health. B.J. is supported by the National Institute of Allergy and Infectious Diseases of the National Institutes of Health under award number F31AI164846. The content presented in this study is solely the responsibility of the authors and does not represent the official views of the NIH.

Author information

Authors and Affiliations

  1. Program in Public Health, University of California Irvine, Irvine, CA, USA

    Brook Jeang, Daibin Zhong, Ming-Chieh Lee & Guiyun Yan

  2. School of Public Health and Community Development, Maseno University, Kisumu, Kenya

    Harrysone Atieli

  3. International Center of Excellence for Malaria Research, Tom Mboya University College, Homa Bay, Kenya

    Harrysone Atieli

  4. School of Medical Laboratory Sciences, Faculty of Health Sciences, Jimma University, Jimma, Ethiopia

    Delenasaw Yewhalaw

  5. Tropical and Infectious Diseases Research Center, Jimma University, Jimma, Ethiopia

    Delenasaw Yewhalaw

Authors
  1. Brook Jeang
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Daibin Zhong
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Ming-Chieh Lee
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Harrysone Atieli
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Delenasaw Yewhalaw
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Guiyun Yan
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

BJ and DZ conducted laboratory analyses and drafted the manuscript with contributions from M-CL and GY. HA, DY, and GY conceived the study and coordinated sample collection. All authors edited and approved the final manuscript.

Corresponding author

Correspondence to Guiyun Yan.

Ethics declarations

Ethics approval and consent to participate

Ethical approval was obtained from the institutional review boards of the University of California, Irvine; Maseno University, Kenya; and Jimma University Institute of Health, Ethiopia. For all studies, approval was first received from the division offices of the Ministry of Health. Written informed consent/ assent was then obtained from all consenting heads of households and for everyone who was willing to participate. All dried blood spot samples were de-identified. For schoolchildren dynamics studies, permission was received from school administrators. Assent was obtained from all students who were willing to participate in the surveys. Individuals who were unwilling to participate were excluded from all studies. Participants harbouring asymptomatic infections were not treated with antimalarials. Symptomatic volunteers were referred to local health facilities for diagnosis and treatment free of charge, in accordance with national policies.

Competing interests

All authors declare that they have 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 file1:

Table S1. Nested Pfk13 PCR reaction mixtures and thermal cycling conditions. Table S2. Age distribution of samples tested for Pfk13 polymorphisms. Table S3. Frequency of synonymous and nonsynonymous Pfk13 polymorphisms in Kenya from 2018 to 2022, stratified by year and infection status. Table S4. Frequency of nonsynonymous Pfk13 polymorphisms in Ethiopia from 2018 to 2022, stratified by year and infection status.

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 http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated 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

Jeang, B., Zhong, D., Lee, MC. et al. Molecular surveillance of Kelch 13 polymorphisms in Plasmodium falciparum isolates from Kenya and Ethiopia. Malar J 23, 36 (2024). https://0-doi-org.brum.beds.ac.uk/10.1186/s12936-023-04812-y

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s12936-023-04812-y

Keywords

Malaria Journal

ISSN: 1475-2875

Contact us