Skip to main content

Clinical isolates of uncomplicated falciparum malaria from high and low malaria transmission areas show distinct pfcrt and pfmdr1 polymorphisms in western Ethiopia

Abstract

Background

Pfcrt gene has been associated with chloroquine resistance and the pfmdr1 gene can alter malaria parasite susceptibility to lumefantrine, mefloquine, and chloroquine. In the absence of chloroquine (CQ) and extensive use of artemether–lumefantrine (AL) from 2004 to 2020 to treat uncomplicated falciparum malaria, pfcrt haplotype, and pfmdr1 single nucleotide polymorphisms (SNPs) were determined in two sites of West Ethiopia with a gradient of malaria transmission.

Methods

230 microscopically confirmed P. falciparum isolates were collected from Assosa (high transmission area) and Gida Ayana (low transmission area) sites, of which 225 of them tested positive by PCR. High-Resolution Melting Assay (HRM) was used to determine the prevalence of pfcrt haplotypes and pfmdr1 SNPs. Furthermore, the pfmdr1 gene copy number (CNV) was determined using real-time PCR. A P-value of less or equal to 0.05 was considered significant.

Results

Of the 225 samples, 95.5%, 94.4%, 86.7%, 91.1%, and 94.2% were successfully genotyped with HRM for pfcrt haplotype, pfmdr1-86, pfmdr1-184, pfmdr1-1042 and pfmdr1-1246, respectively. The mutant pfcrt haplotypes were detected among 33.5% (52/155) and 80% (48/60) of isolates collected from the Assosa and Gida Ayana sites, respectively. Plasmodium falciparum with chloroquine-resistant haplotypes was more prevalent in the Gida Ayana area compared with the Assosa area (COR = 8.4, P = 0.00). Pfmdr1-N86Y wild type and 184F mutations were found in 79.8% (166/208) and 73.4% (146/199) samples, respectively. No single mutation was observed at the pfmdr1-1042 locus; however, 89.6% (190/212) of parasites in West Ethiopia carry the wild-type D1246Y variants. Eight pfmdr1 haplotypes at codons N86Y–Y184F–D1246Y were identified with the dominant NFD 61% (122/200). There was no difference in the distribution of pfmdr1 SNPs, haplotypes, and CNV between the two study sites (P > 0.05).

Conclusion

Plasmodium falciparum with the pfcrt wild-type haplotype was prevalent in high malaria transmission site than in low transmission area. The NFD haplotype was the predominant haplotype of the N86Y–Y184F–D1246Y. A continuous investigation is needed to closely monitor the changes in the pfmdr1 SNPs, which are associated with the selection of parasite populations by ACT.

Background

Plasmodium falciparum resistance to chloroquine, sulfadoxine–pyrimethamine; artemisinin and its partner drugs threatens effective malaria treatment [1, 2]. Worryingly, studies in Africa have indicated the emergence of artemisinin-resistant parasites in eastern region of the continent [3, 4]. Mutations in the P. falciparum chloroquine-resistant transporter (Pfcrt) gene, specifically, K76T including mutations in three other amino acids in the adjoining K76 (residues 72, 74 and 75), region are considered to contribute to chloroquine and amodiaquine resistance [5]. After CQ withdrawal, studies have shown that chloroquine resistance isolates have gradually reduced in frequency with the circulation of chloroquine sensitive strains [6] as has been reported from different parts of Africa, 84.1% of the parasites became the wild type in Ethiopia [7], 76.3% of isolates were chloroquine sensitive in Kenya [8] and another study in Kenya also showed 93.3% of the parasites were the wild type in 2018 [9], the wild-type Pfcrt haplotype CVMNK was found in 22 of the 26 isolates in Niger [10], all detected parasites were the wild type in Malawi in 2001 [11], and a significant reduction in Pfcrt 76T from 97.0 to 66.9% in Cameroon [12].

Most previous studies in Ethiopia have shown the fixation of pfcrt 76T mutation despite withdrawal of CQ for treatment of uncomplicated falciparum malaria [13,14,15]. In Ethiopia, the low rate of reversion in CQ-sensitive parasites is associated with the continued use of CQ to treat vivax malaria [16]. Knowledge on the re-emergence of CQ-sensitive parasites would be essential towards the likely reuse of the safe and cheap CQ [17].

Plasmodium falciparum multidrug resistance 1 (pfmdr1) plays a central role in parasite resistance to artemisinin-based combination therapy (ACT) partner drugs, such as mefloquine, lumefantrine and amodiaquine through single nucleotide polymorphismsand/or gene copy number variations [18, 19]. Thus, protection of the partner drug efficacy is crucial for maintaining the effectiveness of ACT [20]. Mutations at positions N86Y, Y184F, S1034, N1042D, and D1246Y of Pfmdr1 are suggested to be involved in altered drug transport from the parasite’s cytosol into the digestive vacuole [21]. N86Y mutation increases parasite susceptibility to the partner drugs lumefantrine and mefloquine and the active artemisinin metabolite dihydroartemisinin and conversely augments resistance to the ACT partner drug amodiaquine and the former first-line agent CQ [19]. Studies from malaria endemic sites have reported a significant increase in the wild-type pfmdr1 N86 and D1246 following ACT [8, 22]. The pfmdr1, NFD haplotype (codons N86, 184F and D1246), has been associated with reduced sensitivity to lumefantrine [23] and a temporal increase in the Pfmdr1 NFD haplotype was reported from different parts of Africa [24, 25].

The dominance of the wide type pfmdr-1 N86 and the mutant 184F alleles have been reported from different parts of Ethiopia; all the parasites carried the wild-type Pfmdr1 N86 [15], 77.3–100% of isolates collected from different sites of Ethiopia carried the Pfmdr1-N86 allele (16), and 98.8% and 100% of the parasites carried the wild type Pfmdr1 N86 and the mutant 184F, respectively [26]. Similarly, a study in South East Ethiopia indicated all the isolates had the wild-type pfmdr-1 D1246 [7].

Pfmdr1 copy number variation (CNV) has been associated with differential susceptibility to anti-malarial, with multiple copies suggesting reduced sensitivity to lumefantrine and mefloquine [27, 28]. In Cambodia, for instance, pfmdr1 CNV was associated with ACT treatment failure [29]. A study conducted in Kenya showed that post-ACT, parasites had significantly higher pfmdr1 CNV compared to pre-ACTs [30]. In addition, a study from Sudan showed that pfmdr1gene amplification before ACT was associated with recurrent infection during follow-up [31]. Studies have shown different proportion of multiple pfmdr1 copies, 5.5% in Mynamar [32], 11% in Suriname [33], 20% in Venezuela [34] and 54.2% in Ethiopia [35].

Measuring the prevalence of pfmdr1 polymorphisms has implications for designing policies such as drug cycling, sequential artemisinin-based combination treatments or multiple first-line therapies [22]. Moreover, monitoring of drug resistance in Plasmodium populations is crucial for malaria control and elimination [2]. This study was initiated to determine mutations in pfcrt and pfmdr1 genes among the clinical isolates of uncomplicated falciparum malaria patients in western Ethiopia's with gradient of malaria transmission.

Methods

Study site and period

This study was carried out from September through December 2020 in two areas of western Ethiopia, Assosa and Anger Gute, with different levels of malaria transmission. From the high transmission area, the samples were collected from Sherkole and Horazhab health centres of Assosa zone in the Benishangul-Gumuz Region of Ethiopia. Sherkole health centre is located in the Sherkole district which is bordered by Sudan in the north. Sherkole town is about 754 km from Addis Ababa. Horazhab health centre is located in Kurmuk district, and it is bordered by Sudan in the north and west. It is about 769 km from Addis Ababa. Sherkole and Kurmuk districts border Eth-Sudan where there is human mobility between the two countries. In addition, there is a Sherkole Refugee camp that inhabitant more than 14,000 Sudan and South Sudan immigrants [36]. So, there is a chance of imported P. falciparum from Sudan and South Sudan that might affect the pfcrt and pfmdr1 polymorphism in those study districts. Sherkole and Kurmuk districts are located in areas along the western borders with Sudan, characterized by high malaria transmission intensity [37, 38]. In Ethiopia, high malaria-risk areas are mainly located in the western lowland areas of the country [39] (Fig. 1).

Fig. 1
figure 1

Ethiopia Malaria risk stratification, 2020 [39]

From the low transmission area, the samples were collected from Anger Gute and warabo health centers located in Gida Ayana district, East Wollega Zone, Oromia Regional State (Fig. 2). In the Anger Gute area malaria transmission is low and stable and P. falciparum infection prevalence among children 2–10 years was < 5% [37]. The incidence of malaria in and around Anger Gute town was 3.43 per 1000 population at risk of the disease, and the malaria trend from 2014 to 2018 indicated nearly unchanged numbers of malaria cases [40].

Fig. 2
figure 2

Study area map, Gida Ayana, Sherkole, and Kurmuk districts, Western Ethiopia

Study design and population

A health facility-based cross-sectional study was conducted. At each study health center, patients with uncomplicated malaria whose age was greater than 6 months, and who were microscopically-confirmed to be infected with P. falciparum, were enrolled.

Sample size

The sample size was determined using a single population formula, a 13.1% prevalence of malaria in the Benishangul-Gumuz regional state [41], a 95% confidence level and a 5% precision. Accordingly, the calculated sample size was 175. With a 23% adjustment for the non-response rate, 216 uncomplicated falciparum malaria patients were included in the study. The numbers of patients included in the study were proportional to the number of confirmed uncomplicated falciparum malaria cases reported from the study health facilities in 2019. Accordingly, a total of 230 P. falciparum isolates, 112, 51, and 67, were collected from Sherkole, Kurmuk, and Gida Ayana districts, respectively.

Microscopic diagnosis

For confirmatory diagnosis of P. falciparum, capillary blood samples were collected and used for preparation of thick and thin blood films for slide microscopy examination. Each slide was stained with 10% Giemsa for 10 min, and 100 fields were examined before designating a negative sample [42]. Once the patients were diagnosed for P. falciparum infection as determined by microscope, they gave finger-prick blood samples to prepare dried blood spots (DBS) on filter paper (Whatman No. 1001 320, International Ltd. Maidstone, England) for molecular analyses. The DBS were individually kept in plastic bags with desiccants until molecular biology analysis. Patients with mixed infections with other Plasmodium species were excluded from the study.

Molecular genotyping

Parasite genomic DNA extraction was done from DBS using the Chelex protocol as earlier described [43]. NanoDrop™ Spectrophotometer did quantification of the DNA for extracted samples. All samples were normalized to 30 ng/µL DNA concentrations. Molecular genotyping was done at Medical Research Council Unit Gambia at the London School of Hygiene and Tropical Medicine.

Plasmodium falciparum detection was performed by var gene acidic terminal sequence (varATS) real-time PCR as previously described [44]. Parasites drug resistance genotyping assays were done using high-resolution melting assay (HRM) with a Light-Cycler®480 real-time polymerase chain reaction (PCR) machine and previously described procedure [45].We evaluated loci of pfcrt associated with chloroquine resistance, Pfcrt haplotype C72/M74/N75/K76 was considered wild-type and any polymorphism at any of these amino acids positions was considered mutant by HRM analysis.

Pfmdr1 codon positions 86, 184, 1042 and 1246 were analysed for SNPs. In addition, the detection of Pfmdr1 CNV was performed using a Bio-Rad CFX96 real-time PCR (Qiagen, Valencia, USA) as per the published protocol [46] and calculated using the formula 2−ΔCt (ΔCt = Ct PfMDR1-Ctβ-tubulin). PfMDR-1 copy number > 1.6 is defined as an amplification of the gene. All PCR-positive samples were examined by HRM across the two drug-resistance loci.

Data analysis

The prevalence of Pfcrt wild or mutant parasites was determined by counting the number of samples observed with wild or mutant by HRM analysis divided by the total number of samples that were successfully genotyped by HRM for Pfcrt locus. Similarly, the prevalence of different SNPs in the pfmdr1 gene was determined by counting the number of samples with a particular allele determined by HRM divided by the total number of samples that were successfully genotyped by HRM in a given codon.

Frequencies of pfmdr1 haplotypes at codons N86Y, Y184F and D1246Y and Pfmdr1copy numbers were determined. Logistic regression analysis was done to show an association between the antimalarial drug resistance markers between the study sites. A P-value of less or equal to 0.05 was considered statistically significant.

Ethical considerations

Ethical clearance was obtained from the Ethiopian National Ethics Review Committee and Addis Ababa University, Aklilu Lemma Institute of Pathobiology, IRB. The relevant regional and district health authorities sought permission to conduct the study at the health facilities. Written informed consent was obtained from adult study participants and a parent or guardian of a child. Written informed assent was also taken from children.

Results

Out of 230 microscopically confirmed P. falciparum cases enrolled in the study, 225 of them tested positive by PCR. Of these PCR positive samples, 72.4% (163/225) and 27.6% (62/225) were collected from Asossa and Gida Ayana sites, respectively. The mean age of the study participants was 17.8 ± 12.7 years and with an age range of 7 months to 75 years old. The ratio of males to females was 1.27: 1.

HRM of pfcrt genes, including residues 72–76, was successfully for 95.5% (215/225) of the samples. Overall, the proportions of parasites carrying the wild-type CVMNK and the mutant haplotypes were 51.6% (111/215) and 46.5% (100/215), respectively. Mixed infections accounted for 1.9% (4/215). A total of 33.5% (52/155) and 80% (48/60) of isolates collected from the Assosa and Gida Ayana sites carry the mutant haplotypes, respectively suggesting that isolates circulating at Gida Ayana area carry the mutant genotypes as compared with the Assosa area (COR = 8.4, P = 0.00) (Table 1).

Table 1 Pfcrt haplotypes distribution and association with study sites in West Ethiopia

For the pfmdr1 gene, 94.4%, 86.7%, 91.1% and 94.2% samples were successfully genotyped at codons 86, 184, 1042 and 1246, respectively. The N86Y wild type was detected among 79.8% (166/208) of the samples, of which 77% (114/148) were from Assosa and 86.7% (52/60) were from Gida Ayana. The Y184F mutation was found in 73.4% (146/199) samples, of which 105 were from Assosa and 41 were from Gida Ayana. None of the isolates carry N1042D mutation in this study. In addition, 89.6% (190/212) of parasites in west Ethiopia carry the wild-type D1246Y variants. There was no significant difference in pfmdr1 SNPs between the two study sites (P > 0.05) (Table 2).

Table 2 pfmdr1 SNPs distribution and association with study sites in western Ethiopia

Pfmdr1 haplotypes were determined by combining SNPs at pfmdr1 codons 86, 184 and 1246 and 89% (200/225) of samples were successfully genotyped at the three pfmdr1 loci. Eight pfmdr1 haplotypes at the N86Y–Y184F–D1246Y were identified among P. falciparum isolates in West Ethiopia. The NFD haplotype was observed in 61% (122/200) of the isolates and the NYD wild-type haplotype was found at a frequency of 17% (34/200) (Fig. 3).

Fig. 3
figure 3

Distribution of pfmdr1 haplotypes at N86Y–Y184F–D1246Y among P. falciparum isolates collected from two sites in west Ethiopia

While the haplotypes NFD, NFY, NYD and YFD were identified in both study areas; the YYD, YFY and YYY haplotypes were found only in clinical isolates collected from the Assosa site. More importantly, there was no statistically significant difference in the distribution of N86Y–Y184F–D1246Y haplotypes between the two study sites (P > 0.05) (Table 3).

Table 3 Distribution of pfmdr1haplotypes at N86Y–Y184F–D1246Y among P. falciparum isolates in West Ethiopia

From a total of 225 samples that were successfully analysed for CNV in the pfmdr1 gene, the proportion of parasites with multiple copies greater than 1.6 was 8.4% (19/225). The maximum copy number detected was 4.8. No statistically significant difference in pfmdr1 CNV between the studies sites (P > 0.05) (Table 4).

Table 4 Distribution of pfmdr1 gene copy number variation among P. falciparumisolates in West Ethiopia

Discussion

Surveillance of anti-malarial molecular markers is important to scan the persistence of known mutations and/or temporal genetic changes in the parasite population exposed to changing drug pressures. This study compares directional changes of pfcrt haplotype and pfmdr1 genes polymorphism among parasites in the western part of Ethiopia that have been exposed to similar drug pressures but found in different levels of malaria transmission settings.

The wild-type pfcrt haplotype was detected among 51.6% (111/215) of P. falciparum isolates. This study showed the gradual return of the wild-type strains in west Ethiopia after more than two decades of removal of CQ pressure. This finding was in contrast to studies conducted in other parts of Ethiopia; South-central Oromia [13], East Shoa and West Arsi [14] and Southern Ethiopia [47] were all isolates were mutant and also with studies in East Shewa zone [15], Gonder [35] and North, south and east part of Ethiopia [48] that shown majority of parasites was mutant. The observed level of variation in a reversal of mutations associated with CQ resistance might be related to the dominance of P. falciparum infection in western Ethiopia compared with high levels of Plasmodium vivax co-endemicity in the country, reflecting the discrepancy of drug pressure in different geographical areas of the country [16].

The proportion of P. falciparum isolates with mutant pfcrt haplotype was higher in the Gida Ayana area compared with Assosa. This variation in the study sites might relate to the difference in the level of malaria transmission between the two sites. The low malaria transmission in the Gida Ayana area might fix the chloroquine-resistant P. falciparum as compared to the high malaria transmission site, Assosa. Spread of P. falciparum resistance malaria has been associated with intensity of malaria transmission [49, 50]. Thus, after withdrawal of chloroquine, the spread of chlorquine sensitive parasites is higher in Asossa where malaria transmission is high compared with Gida Ayana in which transmission is low. On the other hand, the difference might also associate with the different levels of P. vivax co-endemicity between the two study areas and the continued availability of CQ for treating P. vivax [16]. In 2019/2020, the prevalence of P. vivax infection in the study districts was 24.6% (167/680) in the Gida Ayana district [51] compared with 0.1% (74/12,358) and 1.9% (72/3860) in Sherkole and Kurmuk districts of Assosa, respectively [52]. Therefore, the fixation of the CQ-resistant haplotype in the Anger Gute area was likely related to CQ drug pressure used to treat P. vivax infections.

The wild-type Pfmdr1-N86 allele that related to reduce sensitivity to lumefantrine was observed in 79.8% (166/208) of the investigated P. falciparum clinical isolates. This prevalence was in agreement with a past study reported from the study site, Benishangul-Gumuz [16]. Our finding was also consistent with other studies conducted in other parts of Ethiopia; Southeast Ethiopia [7], Adama and Gambella [14]. This study also showed that 73.4% (146/199) of the parasites carried the 184F variant that related to reduced susceptibility to mefloquine and lumefantrine. This finding was comparable to previous reports from North, South and east Ethiopia (48) and Southwest Ethiopia (26). Thus, a successive use of AL treatment has induced the selection of the pfmdr1-N86 and pfmdr1-184F strains that reduce parasite susceptibility to lumefantrine in the study area.

At codon N1042D, no mutation was observed in clinical isolates collected from west Ethiopia. The absence of mutation at this codon was previously reported from Southeast Iran [53], Saudi Araba [22] and Vanuatu and Solomon Islands [54]. However, the finding was in contrast to studies from South America that reported mutation of the N1042D [20, 55], and this SNP contributes to resistance to quinine [56]. Therefore, distribution of pfmdr1-N1042D mutation might differ over different geographic areas due to the discrepancy in the genetic variation of the parasites and the type of stress that each strain had experienced. It also demonstrates the low or no contribution of the N1042D mutation to antimalarial drug pressure in the study area.

The majority (89.6%) of parasites in west Ethiopia carried the wild-type D1246 allele. The observed high level of the wild-type D1246 allele was in line with previous studies conducted in Southeast Ethiopia [27], NorthernUganda [57], and East-Central Gabon [58]. Therefore, the high frequency of D1246 was likely associated with selective pressures of AL treatment in the study area as parasites with genotypes N86, 184F, and D1246 are partially resistant to the treatment [23].

This study identified the presence of eight pfmdr1haplotypes at the N86Y–Y184F–D1246Y circulating in west Ethiopia, with the NFD haplotype dominant (61%) followed by the wild type NYD (13.7%). The dominance of the NFD haplotype in this study was in agreement with previous studies in Africa [9, 23, 59]. Thus, a widespread use of AL as first-line anti-malarial treatment raises the NFD haplotype selection that reduces parasite sensitivity to lumefantrine. Although AL remained efficacious in the study area, persistence of residual submicroscopic P. falciparum parasitaemia following the treatment need continuing surveillance using genetic markers may help track the spread of lumefantrine-resistant parasites [60].

In the present study, 8.4% of isolates harbor multiple copies of pfmdr1. This prevalence was comparable to studies done in Kenya [30], Ghana [59], Mynamar [32] and Suriname [33]. However, it was lower than the amplification of the gene reported in North part of Ethiopia [35], Venezuela [34] and Thailand [61]. The difference might be related to underlying genetic background of the parasites and drug pressure caused by antimalarial drugs used in the study sites.

There were no significant differences in the distributions of pfmdr1 SNPs at codons 86, 184, 1042 and 1246. Similarly, pfmdr1 haplotypes and CNV were not associated with study sites (P > 0.05). Thus, frequencies of the mutations are not related to geographic selection and level of malaria transmission but are associated with rates of selection to AL treatment [23].

Conclusion

There is a gradual regaining of chloroquine-sensitive haplotype, and the return is more in areas of high transmission than low transmission sites in western Ethiopia with principally P. falciparum infections. A high prevalence of the wild-type alleles N86, D1042 and D1246 and of the mutant-type allele 184F was detected in the study sites. The NFD haplotype was the predominant haplotype of the N86Y–Y184F–D1246Y, and 8.4% of the parasite carried multiple copies of the pfmdr1 gene. Continuous investigation is needed to closely monitor the changes in the pfmdr1 SNPs associated with the selection of P. falciparum parasites populations by ACT.

Availability of data and materials

The datasets and analysed results of the study are available from the corresponding author and can be obtained on reasonable request.

Abbreviations

ACT:

Artemisinin-based combination therapy

AL:

Artemether–lumefantrine

CQ:

Chloroquine

Pfcrt :

P. falciparum chloroquine-resistant transporter

CNV:

Copy number variations

DBS:

Dried blood spots

HRM:

High-resolution melting assay

IRB:

Institute Review Board

PCR:

Polymerase chain reaction

pfmdr1 :

P. falciparum multidrug resistance 1

VarATS :

Var gene acidic terminal sequence

SNP:

Single nucleotide polymorphisms

References

  1. WHO. Global report on antimalarial drug efficacy and drug resistance: 2000–2010. [Internet]. Geneva, World Health Organization; 2010. https://reliefweb.int/report/world/global-report-antimalarial-efficacy-and-drug-resistance-2000-2010?.

  2. Menard D, Dondorp A. Antimalarial drug resistance: a threat to malaria elimination. Cold Spring Harb Perspect. 2017;7:a025619.

    Article  Google Scholar 

  3. Uwimana A, Legrand E, Stokes BH, Ndikumana JM, Warsame M, Umulisa N, et al. Emergence and clonal expansion of in vitro artemisinin-resistant Plasmodium falciparum kelch13 R561H mutant parasites in Rwanda. Nat Med. 2020;26:1602–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Tumwebaze PK, Conrad MD, Okitwi M, Orena S, Byaruhanga O, Katairo T, et al. Decreased susceptibility of Plasmodium falciparum to both dihydroartemisinin and lumefantrine in northern Uganda. Nat Commun. 2022;3:6353.

    Article  Google Scholar 

  5. Awasthi G, Krishna BK, Das A. Pfcrt haplotypes and the evolutionary history of chloroquine-resistant Plasmodium falciparum. Mem Inst Oswaldo Cruz. 2012;107:129–34.

    Article  PubMed  Google Scholar 

  6. Pelleau S, Moss EL, Dhingra SK, Volney B, Casteras J, Gabryszewski SJ, et al. Adaptive evolution of malaria parasites in French Guiana: reversal of chloroquine resistance by acquisition of a mutation in pfcrt. Proc Natl Acad Sci USA. 2015;112:11672–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Kebede S, Aseffa A, Berhe N, Teklehaymanot T, Clouse RM, Gebru T, et al. Return of chloroquine-sensitive Plasmodium falciparum parasites and emergence of chloroquine-resistant Plasmodium vivax in Ethiopia. Malar J. 2014;13:244.

    Article  Google Scholar 

  8. Kishoyian G, Njagi EN, Orinda GO, Kimani FT. Chloroquine sensitivity and prevalence of chloroquine-resistant genes pfcrt and pfmdr-1 in Western Kenya after two decades of chloroquine withdrawal. Ann Med Health Sci. 2018;8:331–5.

    Google Scholar 

  9. 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:1–21.

    Article  Google Scholar 

  10. Salissou A, Zamanka H, Binze BB, Rivière T, Tichit M, Ibrahim ML, et al. Low prevalence of pfcrt resistance alleles among patients with uncomplicated falciparum malaria in Niger six years after chloroquine withdrawal. Malar Res Treat. 2014;2014:614190.

    PubMed  PubMed Central  Google Scholar 

  11. Kublin JG, Cortese JF, Njunju EM, Mukadam RA, Wirima JJ, Kazembe PN, et al. Reemergence of chloroquine-sensitive Plasmodium falciparum malaria after cessation of chloroquine use in Malawi. J Infect Dis. 2003;187:1870–5.

    Article  PubMed  Google Scholar 

  12. Moyeh MN, Njimoh DL, Evehe MS, Ali IM, Nji AM, Nkafu DN, et al. Effects of drug policy changes on evolution of molecular markers of Plasmodium falciparum resistance to chloroquine, amodiaquine, and sulphadoxine-pyrimethamine in the South West Region of Cameroon. Malar Res Treat. 2018;2018:7071383.

    PubMed  PubMed Central  Google Scholar 

  13. Golassa L, Enweji N, Erko B, Aseffa A, Swedberg G. High prevalence of pfcrt-CVIET haplotype in isolates from asymptomatic and symptomatic patients in south-central Oromia. Ethiopia Malar J. 2014;13:120.

    Article  PubMed  Google Scholar 

  14. Golassa L, Kamugisha E, Ishengoma DS, Baraka V, Shayo A, Baliraine FN, et al. Identification of large variation in pfcrt, pfmdr-1 and pfubp-1 markers in Plasmodium falciparum isolates from Ethiopia and Tanzania. Malar J. 2015;14:264.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Hassen J, Solomon G, Dinka H, Golassa L. High prevalence of Pfcrt76T and Pfmdr1 N86 genotypes in malaria infected patients attending health facilities in East Shewa zone, Oromia Regional State, Ethiopia. Malar J. 2022;21:286.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Hailemeskel E, Menberu T, Shumie G, Behaksra S, Chali W, Keffale M. Prevalence of Plasmodium falciparum Pfcrt and Pfmdr1 alleles in settings with different levels of Plasmodium vivax co-endemicity in Ethiopia. Int J Parasitol Drugs Drug Resist. 2019;11:8–12.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Roux AT, Maharaj L, Oyegoke O, Akoniyon OP, Adeleke MA, Maharaj R, et al. Chloroquine and sulfadoxine–pyrimethamine resistance in sub-Saharan Africa: a review. Front Genet. 2021;12:668574.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Qidwai T. Exploration of copy number variation in genes related to anti-malarial drug resistance in Plasmodium falciparum. Gene. 2020;736:144414.

    Article  CAS  PubMed  Google Scholar 

  19. Arya A, Kojom LP, Chaudhry S, Sharma A, Singh V. Artemisinin-based combination therapy (ACT) and drug resistance molecular markers: a systematic review of clinical studies from two malaria endemic regions India and sub-Saharan Africa. Int J Parasitol Drugs Drug Resist. 2021;15:43–56.

    Article  CAS  PubMed  Google Scholar 

  20. Ibraheem ZO, AbdMajid RA, Noor SM, Sedik HM, Basir R. Role of different Pfcrt and Pfmdr-1 mutations in conferring resistance to antimalaria drugs in Plasmodium falciparum. Malar Res Treat. 2014;2014:950424.

    PubMed  PubMed Central  Google Scholar 

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

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Al-Mekhlafi HM, Madkhali AM, Ahmed A, Abdulhaq AA, Atroosh WM, Ghzwani AH, et al. Polymorphism analysis of pfmdr1gene in Plasmodium falciparum isolates 11 years post-adoption of artemisinin-based combination therapy in Saudi Arabia. Sci Rep. 2022;12:517.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Okell LC, Reiter LM, Ebbe LS, Baraka V, Bisanzio D, Watson OJ, et al. Emerging implications of policies on malaria treatment: genetic changes in the Pfmdr-1 gene affecting susceptibility to artemether–lumefantrine and artesunate–amodiaquine in Africa. BMJ Glob Health. 2018;3:e000999.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Achieng AO, Muiruri P, Ingasia LA, Opot BH, Juma DW, Yeda R, et al. Temporal trends in prevalence of Plasmodium falciparum molecular markers selected for by artemether lumefantrine treatment in pre-ACT and post-ACT parasites in western Kenya. Int J Parasitol Drugs Drug Resist. 2015;5:92–9.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Thomsen TT, Madsen LB, Hansson HH, Tomás EV, Charlwood D, Bygbjerg IC, et al. Rapid selection of Plasmodium falciparum chloroquine resistance transporter gene and multidrug resistance gene-1 haplotypes associated with past chloroquine and present artemether–lumefantrine use in Inhambane District, Southern Mozambique. Am J Trop Med Hyg. 2013;88:536–41.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Heuchert A, Abduselam N, Zeynudin A, Eshetu T, Loscher 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 

  27. Shah NK, Alker AP, Sem R, Susanti AI, Muth S, Maguire JD, et al. Molecular surveillance for multidrug-resistant Plasmodium falciparum. Cambodia Emerg Infect Dis. 2008;14:1637–40.

    Article  PubMed  Google Scholar 

  28. Srimuang K, Miotto O, Lim P, Fairhurst RM, Kwiatkowski DP, Woodrow CJ, et al. Analysis of anti-malarial resistance markers in pfmdr1 and pfcrt across Southeast Asia in the tracking resistance to artemisinin collaboration. Malar J. 2016;15:541.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Lim P, Alker AP, Khim N, Shah NK, Incardona S, Doung S, et al. Pfmdr1 copy number and arteminisin derivatives combination therapy failure in falciparum malaria in Cambodia. Malar J. 2009;8:11.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Ngalah BS, Ingasia LA, Cheruiyot AC, Chebon LJ, Juma DW, Muiruri P, et al. Analysis of major genome loci underlying artemisinin resistance and pfmdr1 copy number in pre- and post-ACTs in Western Kenya. Sci Rep. 2013;5:8308.

    Article  Google Scholar 

  31. Gadalla NB, Adam I, Elzaki S, Bashi S, Mukhtar I, Oguike M, et al. Increased pfmdr1 copy number and sequence polymorphisms in Plasmodium falciparum isolates from Sudanese malaria patients treated with artemether–lumefantrine. Antimicrob Agents Chemother. 2011;55:5408–11.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Win AA, Imwong M, Kyaw MP, Woodrow CJ, Chotivanich K, Hanboonkunupakarn B, et al. K13 mutations and pfmdr1 copy number variation in Plasmodium falciparum malaria in Myanmar. Malar J. 2016;15:110.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Labadie-Bracho M, Adhin MR. Increased pfmdr1 copy number in Plasmodium falciparum isolates from Suriname. Trop Med Int Health. 2013;18:796–9.

    Article  CAS  PubMed  Google Scholar 

  34. Pacheco C, Moreno J, Herrera F. A high number of pfmdr1 gene copies in P. falciparum from Venezuela. Parasitol Res. 2019;118:3085–9.

    Article  PubMed  Google Scholar 

  35. Tajebe A, Aemero M, Francis K, Magoma G. Identification of chloroquine resistance Pfcrt-K76T and determination of Pfmdr1-N86Y copy number by SYBR Green I qPCR. Asian Pac J Trop Biomed. 2015;5:208–20.

    Article  CAS  Google Scholar 

  36. UNHCR. Ethiopia Operation. Sherkole camp profile. [Internet]; 2020. https://data.unhcr.org/en/documents/details/83700.

  37. Ministry of Health, Federal Democratic Republic of Ethiopia. An epidemiological profile of malaria in Ethiopia [Internet]. Addis Ababa, 2014. https://web-archive.lshtm.ac.uk/www.linkmalaria.org/sites/link/files/content/country/profiles/Ethiopia%20Epi%20Report%20(240314).pdf.

  38. Taffese HS, Hemming-Schroeder E, Koepfli C, Tesfaye G, Lee MC, Kazura J, et al. Malaria epidemiology and interventions in Ethiopia from 2001 to 2016. Infect Dis Poverty. 2018;7:103.

    Article  PubMed  PubMed Central  Google Scholar 

  39. FMOH. Ethiopia Malaria elimination strategic plan: 2021–2025.

  40. Dufera M, Kenea O, Tadele G. Malaria incidence and associated risk factors in and around Anger Gute Town. Western Ethiopia Res Square. 2020. https://0-doi-org.brum.beds.ac.uk/10.21203/rs.3.rs-29458/v1.

    Article  Google Scholar 

  41. Assefa A. The third Ethiopian Malaria indicator survey 2015. EMIS-2015.

  42. WHO. Basic malaria microscopy. Geneva: World Health Organization; 2010.

  43. Simon N, Shallat J, Williams C, Harrington WE. Optimization of Chelex 100 resin-based extraction of genomic DNA from dried blood spots. Biol Methods Protoc. 2020;5:bpaa009.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Hofmann N, Mwingira F, Shekalaghe S, Robinson LJ, Mueller I, Felger I. Ultra-sensitive detection of Plasmodium falciparum by amplification of multi-copy subtelomeric targets. PLoS Med. 2015;12:e1001788.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Daniels R, Ndiaye D, Wall M, McKinney J, Séne PD, Sabeti PC, et al. Rapid, field-deployable method for genotyping and discovery of single-nucleotide polymorphisms associated with drug resistance in Plasmodium falciparum. Antimicrob Agents Chemother. 2012;56:2976–86.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Nyang H. Copy number variation assay: ASSAY-MDEE-025. MRC Unit the Gambia at LSHTM.2021; Version: 1.0.

  47. Schunk M, Kumma WP, Miranda IB, Osman ME, Roewer S, Alano A, et al. High prevalence of drug-resistance mutations in Plasmodium falciparum and Plasmodium vivax in southern Ethiopia. Malar J. 2006;5:54.

    Article  PubMed  PubMed Central  Google Scholar 

  48. 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:e0005806.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Talisuna AO, Okello PE, Erhart A, Coosemans M, Alessandro UD. Intensity of malaria transmission and the spread of Plasmodium falciparum—resistant malaria: a review of epidemiologic field evidence. Am J Trop Med Hyg. 2007;77:170–80.

    Article  PubMed  Google Scholar 

  50. Masserey T, Lee T, Golumbeanu M, Shattock AJ, Kelly SL, Hastings IM, et al. The influence of biological, epidemiological, and treatment factors on the establishment and spread of drug resistant Plasmodium falciparum. Elife. 2022;11:e77634.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Gida Ayana District Health Office. 2019/2020 Annual malaria report.

  52. Benishangul-Gumuz Regional Health Bureau. 2019/2020 Annual malaria report.

  53. Jalousian F, Dalimi A, Samiee SM, Ghaffarifar F, Soleymanloo F, Naghizadeh R. Mutation in pfmdr1 gene in chloroquine-resistant Plasmodium falciparum isolates, Southeast Iran. Int J Infect Dis. 2008;12:630–4.

    Article  CAS  PubMed  Google Scholar 

  54. Gresty KJ, Gray K, Bobogare A, Taleo G, Hii J, Wini L, et al. Genetic mutations in pfcrt and pfmdr1 at the time of artemisinin combination therapy introduction in South Pacific islands of Vanuatu and Solomon Islands. Malar J. 2014;13:406.

    Article  PubMed  PubMed Central  Google Scholar 

  55. Plummer WB, Lexley LP, Carrington CV. Pfcrt and pfmdr1 alleles associated with Chloroquine resistance in Plasmodium falciparum from Guyana, South America. Mem Inst Oswaldo Cruz. 2004;99:389–92.

    Article  CAS  Google Scholar 

  56. Sidhu AB, Valderramos SG, Fidock DA. Pfmdr1 mutations contribute to quinine resistance and enhance mefloquine and artemisinin sensitivity in Plasmodium falciparum. Mol Microbiol. 2005;57:913–26.

    Article  CAS  PubMed  Google Scholar 

  57. Achol E, Ochaya S, Malinga GM, Edema H, Echodu R. High prevalence of Pfmdr1 N86 and D1246 genotypes detected among febrile malaria outpatients attending Lira Regional Referral Hospital, Northern Uganda. BMC Res Notes. 2019;12:235.

    Article  PubMed  PubMed Central  Google Scholar 

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

  59. Duah NO, Matrevi SA, Souza DK, Binnah DD, Tamakloe MM, Opoku VS, et al. Increased pfmdr1 gene copy number and the decline in pfcrt and pfmdr1 resistance alleles in Ghanaian Plasmodium falciparum isolates after the change of anti-malarial drug treatment policy. Malar J. 2013;12:377.

    Article  PubMed  PubMed Central  Google Scholar 

  60. Tadele G, Jaiteh FK, Oboh M, Oriero E, Dugassa S, Amambua-Ngwa A, et al. Persistence of residual submicroscopic P. falciparum parasitemia following treatment of artemether–lumefantrine in Ethio-Sudan Border, Western Ethiopia. Antimicrob Agents Chemother. 2022;66:e0000222.

    Article  PubMed  Google Scholar 

  61. Chaisatit C, Sai-Ngam P, Thaloengsok S, Sriwichai S, Jongsakul K, Fukuda M, et al. Molecular detection of mutations in the propeller domain of kelch 13 and pfmdr1 copy number variation in Plasmodium falciparum isolates from Thailand collected from 2002 to 2007. Am J Trop Med Hyg. 2021;105:1093–6.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

The authors would like to thank all P. falciparum patients for their willingness to participate in the study. We would also like to thank Medical Research Council Unit the Gambia, at the London School of Hygiene and Tropical Medicine, Gambia for the permission to use the molecular biology laboratory. This work was supported by the PAMGEN H3 Africa project (H3A/002/18) and Addis Ababa University, Post Graduate studies.

Author information

Authors and Affiliations

Authors

Contributions

Conceptualization: Dr. LG and GT. Supervision of the overall work: Dr. LG, Dr. AA-N, and Dr. SD. Methodology: GT. Laboratory work: GT, AJ, EO, and MO. Write up of the original draft of the manuscript: GT. Reviewing and editing the manuscript: Dr. LG, Dr. AA-N, Dr. SD, EO, and MO. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Lemu Golassa.

Ethics declarations

Ethics approval and consent to participate

Ethical approval was obtained from the Ethiopian National Ethic review committee, Addis Ababa, Ethiopia (No. MoSHE 04/246/66) and Institutional review board of Aklilu Lemma Institute of Pathobiology, Addis Ababa University (No. ALIPB IRB/19/2012/20). Written informed consent was obtained from adult study participants and a parent or guardian of a child. Written informed assent was also taken from children.

Consent for publication

Not applicable.

Competing interests

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

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

Tadele, G., Jawara, A., Oboh, M. et al. Clinical isolates of uncomplicated falciparum malaria from high and low malaria transmission areas show distinct pfcrt and pfmdr1 polymorphisms in western Ethiopia. Malar J 22, 171 (2023). https://0-doi-org.brum.beds.ac.uk/10.1186/s12936-023-04602-6

Download citation

  • Received:

  • Accepted:

  • Published:

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

Keywords