Malaria Journal volume  21, Article number: 286 (2022 ) Cite this article
Plasmodium falciparum resistance to series of anti-malarial drugs is a major challenge in efforts to control and/or eliminate malaria globally. In 1998, following the widespread of chloroquine (CQ) resistant P. falciparum, Ethiopia switched from CQ to sulfadoxine–pyrimethamine (SP) and subsequently in 2004 from SP to artemether–lumefantrine (AL) for the treatment of uncomplicated falciparum malaria. Data on the prevalence of CQ resistance markers after more than two decades of its removal is important to map the selection pressure behind the targets codons of interest. The present study was conducted to determine the prevalence of mutations in Pfcrt K76T and Pfmdr1 N86Y codons among malaria-infected patients from Adama, Olenchiti and Metehara sites of East Shewa zone, Oromia Regional State, Ethiopia.
Finger-prick whole blood samples were collected on 3MM Whatman ® filter papers from a total of 121 microscopically confirmed P. falciparum infected patients. Extraction of parasite DNA was done by Chelex-100 method from dried blood spot (DBS). Genomic DNA template was used to amplify Pfcrt K76T and Pfmdr1 N86Y codons by nested PCR. Nested PCR products were subjected to Artherobacter protophormiae-I (APoI) restriction enzyme digestion to determine mutations at codons 76 and 86 of Pfcrt and Pfmdr1 genes, respectively.
Of 83 P. falciparum isolates successfully genotyped for Pfcrt K76T, 91.6% carried the mutant genotypes (76T). The prevalence of Pfcrt 76T was 95.7%, 92.5% and 84.5% in Adama, Metehara and Olenchiti, respectively. The prevalence of Pfcrt 76T mutations in three of the study sites showed no statistical significance difference (χ2 = 1.895; P = 0.388). On the other hand, of the 80 P. falciparum samples successfully amplified for Pfmdr1, all carried the wild-type genotypes (Pfmdr1 N86).
Although CQ officially has been ceased for the treatment of falciparum malaria for more than two decades in Ethiopia, greater proportions of P. falciparum clinical isolates circulating in the study areas carry the mutant 76T genotypes indicating the presence of indirect CQ pressure in the country. However, the return of Pfmdr1 N86 wild-type allele may be favoured by the use of AL for the treatment of uncomplicated falciparum malaria.
Despite the reduction of malaria-related morbidity and mortality in recent years, there are 241Â million estimated cases and nearly 627 thousand deaths in 2020. The World Health Organization (WHO) African Region alone contributed to 95% and 96% of the estimated global malaria cases and deaths, respectively. Children under 5 years of age are the most risk group, accounting for 77% of all malaria deaths globally in 2020. Plasmodium falciparum is the major causes of disease and deaths in this region [1].
Plasmodium falciparum developed resistance to every class of anti-malarial drugs to date [2]. Chloroquine resistance (CQR) has been associated with mutations in P. falciparum chloroquine resistance transporter (Pfcrt) gene and in P. falciparum multidrug resistance 1 (Pfmdr1) gene; however the former is a strong predictor of CQR [3, 4]. In vitro and in vivo studies indicated that substitutions of charged lysine (K) amino acid by uncharged amino acid threonine (T) at amino acid position 76 (K76T) is a primary marker of CQR [3, 4]. Other mutations at amino acid position 72, 74, 75, 97, 152, 163, 220, 271, 326, 356 and 371 have been also identified in CQR P. falciparum in different geographical regions [4, 5]. CQR strains of P. falciparum originating from CVMNK (CQ sensitive haplotype) are also characterized by CVIET (mostly in Africa) and SVMNT (Asia and South America) haplotypes spanning codons 72–76 of the Pfcrt gene [6, 7].
Five point mutations in Pfmdr1 gene at positions 754, 1049, 3598, 3622 and 4234 result in amino acid changes at codons 86 (N86Y), 184 (Y184F), 1034 (S1034C), 1042 (N1042D) and 1246 (D1246Y), respectively [8]. These mutations have been implicated in increase resistance to CQ, amodiaquine, mefloquine (MQ), quinine (QN) and halofantrine [9]. Out of these mutations N86Y appears to be the most important as it alters the transport activity of the protein and has been implicated as high resistance to CQ [10]. Meanwhile this mutation increases susceptibility to a number of anti-malarial drugs including lumefantrine (LUM), MQ and dihydroartemisinin [9]. Strong linkage disequilibrium was observed between Pfmdr1 N86Y and Pfcrt K76T alleles in CQR isolates of P. falciparum [8]. Amplification of the Pfmdr1 gene also leads to decreased drug sensitivity to QN, LUM and artemisinin [11].
In Ethiopia, CQ treatment failure for P. falciparum and P. vivax was reported in 1996 [12]. Then, due to widespread of CQ resistant P. falciparum, the country changed its treatment policy from CQ to sulfadoxine–pyrimethamine (SP) in 1998 which in turn was replaced by artemether–lumefantrine (AL) for the treatment of uncomplicated falciparum malaria in 2004; however CQ is retained and still continued to be used as a first-line treatment for vivax malaria [6, 13]. Molecular studies from Ethiopia have showed high prevalence of Pfcrt K76T mutation that ranged from 41.5 to 100% [6, 13, 14]. However, only one study has showed the return of CQ sensitive P. falciparum with 84% of Pfcrt K76 sensitive alleles [15]. On the other hand, in Ethiopia, AL pressure might result in directional selection of Pfmdr1 variants N86, 184 F and D1246 due to LUM partner drug pressure [6, 13,14,15]. Therefore, it is important to continue monitoring the prevalence of these markers (Pfcrt and Pfmdr1) that confer resistance to CQ in the hope of possible re-emergence of CQ sensitive P. falciparum as previously observed in Malawi [16]. Due to co-existence of P. falciparum and P. vivax; mixed infections of both parasites are common in Ethiopia [17] that favor the presence of indirect but sustained CQ drug pressure. Therefore, the aim of the present study was to investigate the prevalence of genotypic resistance markers for CQ (Pfcrt K76T and Pfmdr1 N86Y) in P. falciparum isolates from East Shewa zone, Oromia Regional State, Ethiopia.
Cross-sectional health facility based study was conducted at Adama (8° 32′ 31″ N and 39° 16′ 33″ E), Olenchiti (8° 38′ 16″ N and 39° 24′ 19″ E) and Metehara (8° 86′ 95″ N and 39° 92′ 02″ E) health facilities of East Shewa zone, Oromia Regional State, Ethiopia. Adama the capital town of East Shewa zone is found at around 99 km from Addis Ababa, the capital city of Ethiopia. Olenchiti and Metehara towns are found to the East of Adama town at about 25 and 91 km, respectively (Fig. 1). The current study was conducted at purposively selected health facilities due to their physical accessibility to Adama malaria research and diagnostic centre to conduct membrane-feeding assay for another experiment.
Map of the study area. Arc GIS version 10.8
The areas are characterized by low altitude (< 2000 m above sea level), heavy rainfall pattern from June to August followed by major transmission of malaria from September to December and short rains of February and March contributing for minor malaria transmission from April to May [18]. Geographical location of the study areas in Great Rift Valley, at low altitude, with rainfall pattern and average annual temperature of 16–32 °C make the areas favorable environments for breeding of Anopheles mosquitoes, of which Anopheles arabiensis is the major malaria vector [19]. Both P. falciparum and P. vivax co-exist in the area as it is common in different parts of the country.
Blood samples from one hundred twenty-one microscopically confirmed P. falciparum infected patients were collected between October, 2019 and September, 2021. Inclusion criteria include patients who presented with uncomplicated falciparum malaria and history of fever onset since 24 h of clinical examination. Blood samples were collected after obtaining consents from P. falciparum infected patients or their parents or legal guardians in case of children less than 18 years of age. After consent was obtained patients sex, age and residency were recorded into log-book. Finger-prick blood samples were collected onto microscope slides and Whatman-3MM ® filter papers by experienced laboratory technologists. Filter papers were allowed to air dry in horizontal position at room temperature for at least 4 h away from dust, insects and direct sun light. Each filter paper was labelled with patient code and date of collection and kept into zip-locked plastic bags with silica gels to protect the specimens from moisture. These zip-locked plastic bags then packed in a larger plastic bags and stored in refrigerator at 2–8 °C for short period in the field and further stored in deep freezer at − 20 °C at Aklilu Lemma Institute of Pathobiology for long period until DNA analysis for molecular study.
Microscopy diagnosis of both Giemsa-stained thick and thin blood films was conducted according to WHO protocol [20]. At each health facilities slides were read by two experienced laboratory technologists to confirm species identification for P. falciparum. In addition slides re-examinations and parasite counts were performed by two WHO certified malaria microscopists at Adama malaria research and diagnostic centre. Further, asexual stages of the parasite per micro-litre (µl) was determined by counting the number of parasites against 200 white blood cells (WBCs) on a thick blood film assuming a total of standard WBCs count of 8000/µl of blood. Density of the parasite was determined by using the following formula [21]:
For comparison with Pfcrt K76T wild and mutant alleles, parasitaemia was categorized as low (< 1000 parasites/µl of blood), moderate (1000–9999 parasites/µl of blood) and high (≥ 10,000 parasites/µl of blood) [22].
Plasmodium falciparum genomic DNA was extracted from dried blood spot (DBS) by Chelex extraction method as described elsewhere [19]. Extracted DNA was transferred into 0.5 ml tubes and stored at − 20 °C until further use.
A fragment from Pfcrt gene encompassing 76th codon was successively amplified for 83 DNA samples. Nested polymerase chain reaction (PCR) protocol adapted and modified from the work of Mulenga et al. [23]. Briefly, in primary PCR, 20.0 µl reaction mixture consisting of 11.6 µl nuclease free water, 0.2 µl of each forward and reverse primers, 4.0 µl 5X FIREPol® Master Mix (Solis BioDyne) and 4.0 µl of DNA template was used. PCR was conducted in a thermo cycler (VWR) Schmidt, Germany. Amplification was performed under the following conditions: Hot start denaturations at 94oC for 3 min, followed by 30 cycles of 94 °C for 30 s, 56 °C for 30 s and 65 °C for 1 min and final elongation at 65 °C for 3 min. Then, the primary PCR was followed with nested PCR. Twenty microlitre reaction mixture consisting of 14.6 µl nuclease free water, 0.2 µl of each forward and reverse primers (Additional file 1), 4.0 µl 5X FIREPol Master Mix (Solis BioDyne) and 1.0 µl of DNA template was used. The nested PCR was conducted as follows: Hot start denaturation at 94 °C for 5 min, followed by 40 cycles of 94 °C for 30 s, 52 °C for 30 s and 65 °C for 1 min and lastly, final elongation at 65 °C for 3 min.
Pfmdr1 gene was amplified to identify mutation at codon N86Y according to the protocol described by Fontecha et al. [24] with certain modifications. Briefly, both primary and nested PCRs were conducted in a total reaction volume of 18.0 µl consisting of 11.4 µl ddH2O, 1.0 µl of each forward and reverse primers (Additional file 1), 3.6 µl 5X FIREPol Master Mix and 1.0 µl DNA template. In primary PCR, the target gene was amplified starting with hot start denaturation at 95 °C for 3 min, followed by 20 cycles of 95 °C for 30 s, 52 °C for 30 s and 72 °C for 1 min and lastly, final extension at 60 °C for 5 min. Nested PCR amplification was conducted beginning with hot start denaturation at 95 °C for 3 min, followed by 25 cycles of 95 °C for 30 s, 54 °C for 30 s, and 72 °C for 60 s and final extension at 72 °C for 5 min.
For both Pfcrt K76T and Pfmdr1 N86Y codons, PCR conditions were conducted along with nuclease free water (negative control) instead of parasite DNA, 3D7 (CQ sensitive) and K1 (MRA-159G, CQ resistant) clones of P. falciparum genomic DNA (positive controls). The two strains were provided by Ethiopian Public Health Institute and Aklilu Lemma Institute of Pathobiology, respectively. Five microlitre of nested PCR products were resolved in 2.0% and 2.5% agarose, type I gel (Sigma Alorich, USA) stained with ethidium bromide (EtBr) respectively, for Pfcrt K76T and Pfmdr1 N86Y codons. Electrophoresis conditions were conducted at 100 volts for 45 min in a Thermo Scientific apparatus filled with 0.5X TBE buffer solution (BIONEER Corp, Korea). Then, DNA bands were observed through ultra-violet (UV) trans-illuminator (USA). The band sizes were estimated by relating with 100 bp DNA ladder marker (Bionexus). The remaining 15 µl and 13 µl nested PCR products of Pfcrt K76T and Pfmdr1 N86Y, respectively were stored at 2–8 °C until used for genotyping by restriction fragment length polymorphism (RFLP).
Following amplification of the two codons, genotyping by PCR-RFLP was made by digesting the fragments with APoI restriction enzyme (New England Biolabs). Nested PCR products were used without purification procedures as this will represent common situations in malaria settings [25]. Restriction digestion was conducted according to the manufacturer’s protocol. Restriction enzyme digestion was done in 20.0 µl of reaction mixture containing 12.5 µl nuclease free water, 5U ApoI restriction enzyme, 2.0 µl NEBuffer r3.1 and 5 µl nested PCR products. Known positive controls (Pf3D7, wild-type and K1, mutant type) and uncut DNA samples were included in the reactions. The enzyme cuts Pfcrt 76K and Pfmdr1 N86 but not Pfcrt 76T and Pfmdr1 86Y alleles. Cleaved PCR products (bands) of Pfcrt K76T and Pfmdr1 N86Y codons were then visualized with UV trans-illuminator after electrophoresis using 2.0% and 2.5% agarose gel, respectively.
For the subsets of samples that amplifications of Pfcrt and Pfmdr1 genes were not successful, msp1 gene amplification was made. Thirty-eight samples which were microscopically detected as P. falciparum were not amplified particularly for Pfcrt K76T codon. In order to verify the presence of intact genes and quality of DNA for those samples, nested PCR targeting msp-1 of polymorphic region (block-2) was performed as described elsewhere [19].
Data were first entered into Microsoft Office Excel ® 2007 (Microsoft Corporation). Comparisons of prevalence of 76T allele and parasite densities across study sites and associations of K76 and 76T alleles with parasite density were performed using Pearson’s Chi-square test of Statistical Package for Social Science (SPSS) (IBM Inc., Chicago, IL, USA) version 26. Fisher exact test was used to test correlation between Pfcrt 76T alleles and age groups. Statistical significance was set at a P-value < 0.05.
Eighty- three and eighty samples were successfully analysed for Pfcrt K76T and Pfmdr1 N86Y codons, respectively. The patients were aged between 1 and 68 years and the mean age was 23 ± 14.7 years (interquartile range: 13–30). Seventy-percent of the patients were adults older than 15 years of age. Among patients involved in the study 65% and 35% were males and females, respectively. Parasite density ranged from 79 to 157,538 asexual stages/µl of blood, with a parasitaemia geometric mean of 5366 or median 6579 asexual parasites/µl (Table 1).
Average parasites density at Adama, Olenchiti and Metehara sites were 10781.84/µl, 9870.95/µl and 15089.22/µl, respectively; however, the difference is not statistically significant (χ2 = 162.75; P = 0.469) (Fig. 2).
Mean parasites density versus study sites
A total of 121 samples (31 from Adama, 22 from Olenchiti and 68 from Metehara) which were collected from health centres confirmed for P. falciparum isolates were allowed for PCR. Eighty three (68.6%) samples were successfully amplified for Pfcrt K76T codon; 38 samples were not amplified for the codon. From those samples amplified for Pfcrt K76T, 80 were found successfully amplified for Pfmdr1 N86Y codon.
Digestion of the nested PCR products carrying K76 alleles of Pfcrt gene with APoI restriction produced 100 and 45 bps, but remained uncut when 76T mutant allele presents (hence a single band of 145Â bp) while digestion of N86 allele of Pfmdr1 gene with the same restriction enzyme resulted in 239 and 179Â bp fragments (denoting wild-type) (Additional file 2).
Among 83 P. falciparum isolates successfully genotyped for Pfcrt K76T, 91.6% and 8.4% mutant- type allele (76T) and wild-type alleles (K76), respectively were detected (Fig. 3). Accordingly, high prevalence of 76T mutant allele was detected at Adama (95.7%), Olenchiti (84.2%) and Metehara (92.7%) study sites with no statistical significant difference (χ2 = 1.895; P = 0.388). On the other hand low prevalence of wild-type K76 alleles were detected at each health centres (Table 2). Most of them were detected in P. falciparum isolates with moderate parasitaemia level (1000–99,999 parasites/µl). Isolates with 76T allele had a higher mean parasite density of 15,420 parasites/µl while those with K76 allele had a mean parasite density of 7623.14 parasites/µl. However, no significant association was observed between parasite density and allele type (χ2 = 76.53; P = 0.620). Eighty nested PCR products of N86Y codons were analysed by restriction enzyme digestion produced N86 (100%) wild-type allele; no 86Y mutant type-allele was detected (Fig. 3).
Prevalence of wild and mutant alleles for Pfcrt K76T and Pfmdr1 N86Y codons among P. falciparum isolates from Adama, Olenchiti and Metehara study sites
Out of 83 isolates successfully genotyped for Pfcrt K76T codon, 54 and 29 were males and females, respectively. Males were more represented and contributed for higher prevalence (64.5%) of mutant Pfcrt 76T allele. However, the difference was not statistically significant (P = 0.712).
Prevalence of Pfcrt 76T allele was assessed among age groups at Adama, Olenchiti and Metehara sites. High prevalence of Pfcrt 76T alleles was detected among 15–30 age groups with frequencies of 56.3%, 45.5% and 42.1% at Olenchiti, Adama and Metehara sites, respectively (Fig. 4). However, no statistical significance difference of 76T alleles were detected among age groups (Fisher’s exact test, P = 0.69).
Prevalence of Pfcrt 76T CQ-resistant marker among age groups per study sites
Out of thirty-eight (38) samples which were not amplified especially for Pfcrt K76T, 26 (68%) of them were amplified by nested PCR for Pfmsp1 allelic families of polymorphic region (Block-2). However, still amplification was unsuccessful for 12 of them.
To date, the most important challenge to control and/or eliminate malaria globally is because of emerging and wide-spread anti-malarial drug resistance Plasmodium parasites. In this regard a development in drug resistance by P. falciparum to series of anti-malarial drugs is the most challenge in fighting malaria for several years [2]. The current study aimed to investigate the prevalence of Pfcrt K76T and Pfmdr1 N86Y codons, in P. falciparum isolates collected from Adama, Olenchiti and Metehara sites of East Shewa zone, Oromia Regional State, Ethiopia.
In this study, high frequencies of Pfcrt 76T (91.6%) with no significance difference (P = 0.388) was reported from the study areas. Wild-type Pfcrt K76 allele was only detected in 8.4% P. falciparum isolates. The high proportions of Pfcrt 76T allele in this study is in agreement with studies from various parts of Ethiopia [6, 13, 26]. This persistence of high frequency of 76T allele in Ethiopia may be associated with incomplete withdrawal of CQ as it is currently the first-line anti-malarial drug for the treatment of P. vivax which accounts nearly for 40% of malaria cases in Ethiopia. Moreover, the persistence of vivax malaria in the form of relapse during dry season may also expose P. falciparum to CQ pressure [18]. The dominance of P. vivax over P. falciparum in the study area was already documented which might contribute for such high frequency of 76T allele [13, 18].
Another potential probable reason for circulations of high frequency of Pfcrt 76T mutant allele in the study areas may be because of the existence of mixed infections of P. falciparum and P. vivax [17]. Misdiagnosis of malaria [27], such as mixed infections, as P. vivax and treating patients with CQ was commonly encountered during the study period (unpublished observations) and this may expose P. falciparum to indirect CQ pressure selecting for the mutant genotypes.
To the best of our knowledge, in the present study areas there is no base-line data regarding the prevalence of 76T allele at the time when CQ was replaced with AL. However, 2 years later study conducted by Schunk et al. [28] in Southern Ethiopia showed fixation of 76T in parasite isolates in the area. Ten years later, since the implementation of AL, study conducted at Adama by Golassa et al. [6] revealed fixation of 76T mutant allele in parasite population, however; in the current study wild-type K76 allele was recovered in the frequency of 8.4%. In the current low malaria endemic area [29], the mutant allele still may circulate at high frequencies in parasite populations between human host and Anopheles vectors in the study areas. In this low malaria transmission setting, lower rate of sexual recombination parasite genotypes may take place in Anopheles mosquitoes, which might result in more resistant genotypes to be established in the human host. Therefore, malaria transmission involving multiple genotypes is uncommon in low transmission areas. Conversely, in high malaria transmission areas competition between drug sensitive and drug resistant parasites may slow the spread of drug resistant parasites [23, 30]. Study conducted by Hailemeskel et al. [13] which investigated the prevalence of Pfcrt K76T genotype in low and high malaria endemic areas of Ethiopia [29] showed significantly higher prevalence of 76T and K76 alleles in low and high malaria endemic areas, respectively.
In countries like Ethiopia where CQ drug pressure persists for the treatment of vivax malaria, high frequency of the mutant allele is being reported from Yemen [31], Malaysia [22], Nigeria [32], India [33], Bangladesh [34] and Angola [35]. Conversely, progressive decrease in 76T allele and recovery of K76 allele were documented particularly in sub-Saharan African countries such as Malawi [16] Tanzania [36], Kenya [37] and Uganda [38], where AL is being implemented. Though AL is used as first-line treatment for P. falciparum in Ethiopia low recovery of parasites harbouring wild-type K76 was detected in the study areas may be because of reasons already described above.
In this study, associations of Pfcrt K76T mutation with parasite density, age and gender were investigated. Some studies have associated this mutation with parasite density [22, 31]; while others have found no such associations [39]. In the current study we have found no association of these mutations with parasite density. On the other hand, the detection of high prevalence of this mutation among 15–30 age groups in this study may be explained by more exposure of these groups to mosquito bite and CQ pressure than others. However, no significant association was found between K76T mutation and age of patients which is in agreement with findings of Atrosh et al. [22] and Edogun et al. [40]. Similarly, no positive association of this mutation with gender was detected in this study.
In this study, P. falciparum clinical isolates analysed carried the wild-type genotypes for Pfmdr1 N86. The return to the wild-type genotypes was successful in the gene under report after cessation of CQ and use of AL for the treatment of uncomplicated falciparum malaria in Ethiopia. No single mutant Pfmdr1 86Y genotype was detected in the study under report. This finding is comparable with the previous reports from different parts of Ethiopia [6, 13, 15, 41]. Study conducted 2 years later after implementation of AL in Ethiopia showed the dominance of 86Y allele over its wild allele counterpart [28]. The current high prevalence of N86 allele (CQ sensitive) showed the survival advantage of wild-type over the mutant type since the implementation of AL for the treatment of falciparum malaria in Ethiopia as evident the selection of this allele by AL elsewhere [42].
In different African countries after adoption of AL for the treatment of falciparum malaria progressive increase in N86 allele and decline in 86Y allele have been documented [43, 44]. In vivo therapeutic studies have showed that N86 wild-type is more prevalent in recurrent infections following AL treatment [45]. In support of this meta-analysis of thirty-one clinical trials have showed a five-fold risk of parasite recrudescence in patients infected with parasites harbouring wild-type N86 allele following treatment with AL compared to patients infected with 86Y mutant parasites [46]. Similarly, allelic replacement study conducted by Veiga et al. [9] revealed that substitution of N86Y mutation with N86 withstands three-four fold higher LUM (partner drug of artemether) and MQ concentrations compared to 86Y residues. Although there is no clear evidence linking N86 allele with treatment failure, high prevalence of this allele in treatment failures in Angola [47] and parasite tolerance to LUM [9, 42], is a clear warning signal for plausible emergence of resistance against AL.
There was unsuccessful nested-PCR amplification for 38 of the samples; particularly for Pfcrt K76T codon as was the case for other studies [48]. However, amplification of Pfmsp1gene was successful for 26 of them. This indicated that such parasites may have low concentrated intact DNA that could not be amplified. Still, those 12 samples which were not again amplified for Pfmsp1 indicated that they may be microscopically false positive for P. falciparum [27].
The current study is not without limitations. Some of the limitations of the study are small sample size analysed due to logistic constraints, lack of data of copy number variations of the Pfmdr1 -86 as parasites possessing multiply copy of the Pfmdr1 N86 showed decreased susceptibility to MQ, artesunate-amodiaquine and LUM [46] and lack of clinical data on the susceptibility of the parasite to LUM that could be obtained from therapeutic efficacy study.
In conclusion, the present study gives an update on the prevalence of K76T and N86Y alleles in the study areas. A reversal of K76 wild-type allele was detected in non-significant proportions while a fixation of N86 allele was investigated across the study areas. The data obtained from this study calls for continuous monitoring anti-malarial drug resistance against falciparum malaria to maintain efficacy of AL for the treatment of uncomplicated P. falciparum malaria. Based on this assessment we also recommend careful monitoring of CQ usage for P. vivax treatment.
The datasets used and/or analysed for this study are available from the correspondent author on reasonable request.
Plasmodium falciparum chloroquine resistane transporter
Statistical package for social science
WHO. World Malaria Report 2021. Geneva: World Health Organization; 2021.
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Â
Fidock DA, Nomura T, Talley AK, Cooper RA, Dzekunov SM, Ferdig MT, et al. Mutations in the P. falciparum digestive vacuole transmembrane protein Pfcrt and evidence for their role in chloroquine resistance. Mol Cell. 2000;6:861–71.
Article CAS PubMed PubMed Central Google ScholarÂ
Djimde A, Pharm D, Ogobara K, Doumbo MD, Joseph F, et al. A molecular marker for chloroquine resistant falciparum malaria. N Ingl J Med. 2001;344:257–63.
Ibrahim ML, Steenkste N, Khim N, Adam HH, Konate L, Coppee JY, et al. Field-based evidence of fast and global increase of Plasmodium falciparum drug resistance by DNA-microarrays and PCR/RFLP in Niger. Malar J. 2009;8:32.
Article PubMed PubMed Central Google ScholarÂ
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Â
Ikegbunam MN, Nkonganyi CN, Thomas BN, Esimone CO, Velavan TP, Ojurongbe O. Analysis of Plasmodium falciparum Pfcrt and Pfmdr1 genes in parasite isolates from asymptomatic individuals in Southeast Nigeria 11 years after withdrawal of chloroquine. Malar J. 2019;18:343.
Article PubMed PubMed Central Google ScholarÂ
Shrivastava SK, Gupta RK, Mahanta J, Dubey ML. Correlation of molecular markers, Pfmdr1-N86Y and Pfcrt-K76T, with in vitro chloroquine resistant Plasmodium falciparum, isolated in the malaria endemic states of Assam and Arunachal Pradesh, Northeast India. PLoS ONE. 2014;9:e103848.
Article PubMed PubMed Central Google ScholarÂ
Veiga MI, Dhingra SK, Henrich PP, Straimer J, Gna N, Uhlemann A, 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Â
Singh G, Singh R, Urhekar AD, Rane K. Gene sequence polymorphisms mutations in Pfmdr-1 and Pfcrt-o genes of Plasmodium falciparum. Int J Curr Microbiol Appl Sci. 2016;5:451–8.
Cui L, Mharakurwa S, Ndiaye D, Rathod PK, Rosenthal PJ. Antimalarial drug resistance: literature review and activities and findings of the ICEMR Network. Am J Trop Med Hyg. 2015;93:57–68.
Article CAS PubMed PubMed Central Google ScholarÂ
Tulu AN, Webbeg RH, Schellenberg AJS, Bradley DJ. Failure of chloroquine treatment for malaria in the highlands of Ethiopia. Trans R Soc Trop Med Hyg. 1996;44:556–7.
Hailemeskel E, Menberu T, Shumie G, Behaksra S, Chali W, Keffale M, et al. 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Â
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Â
Mekonnen SK, 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 PubMed PubMed Central Google ScholarÂ
Kublin JG, Cortese JF, Njunju EM, Mukadam RAG, 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.
Ketema T, Bacha K, Getahun K, Del Portillo HA, Bassat Q. Plasmodium vivax epidemiology in Ethiopia 2000–2020: a systematic review and meta-analysis. PLoS Negl Trop Dis. 2021;15:e0009781.
Article PubMed PubMed Central Google ScholarÂ
Golassa L, White MT. Population – level estimates of the proportion of Plasmodium vivax blood-stage infections attributable to relapses among febrile patients attending Adama Malaria Diagnostic Centre. Malar J. 2017;16:307.
File T, Chekol T, Solomon G, Dinka H, Golassa L. Detection of high frequency of MAD20 allelic variants of Plasmodium falciparum merozoite surface protein 1 gene from Adama and its surroundings, Oromia, Ethiopia. Malar J. 2021;20:385.
Article CAS PubMed PubMed Central Google ScholarÂ
WHO. Basic Malaria Microscopy. Part I. Learner’s guide, 2nd Edn. Geneva: World Health Organization; 2010.
Nega D, Daniel D, Tefera T, Eshetu T. Prevalence and predictors of asymptomatic malaria parasitemia among pregnant women in the rural surroundings of Arbaminch town. PLoS ONE. 2015;10:e0123630.
Article PubMed PubMed Central Google ScholarÂ
Atrosh WM, Al-Mekhalfi HM, Mahdy MAK, Surin J. The detection of Pfcrt and Pfmdr1 point mutations as molecular markers of chloroquine drug resistance, Pahang, Malaysia. Malar J. 2012;11:251.
Mulenga MC, Sitali L, Ciubotariu II, Hawela MB, Hamainza B. Decreased prevalence of the Plasmodium falciparum Pfcrt K76T and Pfmdr1 N86Y mutations post-chloroquine treatment withdrawal in Katete District, Eastern Zambia. Malar J. 2021;20:329.
Article CAS PubMed PubMed Central Google ScholarÂ
Fontecha G, Pinto A, Archaga O, Betancourth S, Escober L, HenrÃquez J, et al. Assessment of Plasmodium falciparum anti-malarial drug resistance markers in Pfcrt and Pfmdr1 genes in isolates from Honduras and Nicaragua, 2018–2021. Malar J. 2021;20:465.
Article CAS PubMed PubMed Central Google ScholarÂ
Veiga MI, Ferreira PE, Gil PJ. Multiplex PCR-RFLP methods for Pfcrt, Pfmdr1 and Pfdhfr mutations in Plasmodium falciparum. Mol Cell Probes. 2006;20:100–4.
Article CAS PubMed Google ScholarÂ
Mula P, Fernandez-martinez A, Lucio A, De, Ramos JM, Reyes F, Gonzalez V, et al. Detection of high levels of mutations involved in anti-malarial drug resistance in Plasmodium falciparum and Plasmodium vivax at a rural hospital in Southern Ethiopia. Malar J. 2011;10:214.
Article CAS PubMed PubMed Central Google ScholarÂ
Jaleta F, Garoma G, Gerenfes T. Evaluation of malaria microscopy diagnosis performance in public hospitals of Eastern and Central part of Oromia Region, Ethiopia, 2019. Pathol Lab Med Int. 2020;12:1–8.
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Â
MOH. Ethiopia Malaria Elimination Strategic Plan: 2021–2025. Ministry of Health; Addis Ababa, Ethiopia; 2020.
Yobi DM, Kayiba NK, Mvumbi DM, Boreux R, Kabututu PZ, Situakibanza HNT, et al. Molecular surveillance of anti-malarial drug resistance in Democratic Republic of Congo: high variability of chloroquine resistance and lack of amodiaquinoresistance. Malar J. 2020;19:121.
Article CAS PubMed PubMed Central Google ScholarÂ
Al-Mekhlafi AM, Mahdy MAK, Al-Mekhlafi HM, Azazy AA, Fong MY. High frequency of Plasmodium falciparum chloroquine resistance marker (Pfcrt T76 mutation) in Yemen: an urgent need to re-examine malaria drug policy. Parasit Vectors. 2011;4:94.
Article CAS PubMed PubMed Central Google ScholarÂ
Adam R, Mukhtar MM, Abubakar UF, Damudi HA, Muhammad A, Ibrahim SS. Polymorphism analysis of Pfmdr1 and Pfcrt from Plasmodium falciparum isolates in Northwestern Nigeria revealed the major markers associated with antimalarial resistance. Diseases. 2021;9:6.
Article CAS PubMed PubMed Central Google ScholarÂ
Das S, Tripathy S, Chattopadhayay S, Das B, Kar Mahapatra S, Hati AK, et al. Progressive increase in point mutations associates chloroquine resistance: even after withdrawal of chloroquine use in India. IJP Drugs Drug Resist. 2017;7:251–61.
Johora FT, Elahi R, Nima MK, Hossain MS, Rashid H, Kibria MG, et al. Persistence of markers of chloroquine resistance in Plasmodium falciparum in Bangladesh. Am J Trop Med Hyg. 2021;104:276–82.
Article CAS PubMed Google ScholarÂ
Ebel ER, Reis F, Petrov DA, Beleza S. Historical trends and new surveillance of Plasmodium falciparum drug resistance markers in Angola. Malar J. 2021;20:175.
Article CAS PubMed PubMed Central Google ScholarÂ
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 CAS PubMed PubMed Central Google ScholarÂ
Chebore W, Zhou Z, Westercamp N, Otieno K, Shi YP, Sergent SB, et al. Assessment of molecular markers of anti-malarial drug resistance among children participating in a therapeutic efficacy study in Western Kenya. Malar J. 2020;19:291.
Article CAS PubMed PubMed Central Google ScholarÂ
Balikagala B, Sakurai-Yatsushiro M, Tachibana SI, Ikeda M, Yamauchi M, Katuro OT, et al. Recovery and stable persistence of chloroquine sensitivity in Plasmodium falciparum parasites after its discontinued use in Northern Uganda. Malar J. 2020;19:76.
Article CAS PubMed PubMed Central Google ScholarÂ
Afoakwah R, Boampong JN, Egyir-Yawson A, Nwaefuna EK, Verner ON, Asare KK. High prevalence of PfCRT K76T mutation in Plasmodium falciparum isolates in Ghana. Acta Trop. 2014;136:32–6.
Article CAS PubMed Google ScholarÂ
Edogun HA, Daramola GO, Ojerinde AO, Esan CO, Adegbuyi AT. Prevalence of mutant alleles responsible for chloroquine resistance among Plasmodium falciparum isolates in North Central, Nigeria. J Adv Biol Biotechnol. 2019;20:1–8.
Lo E, Hemming-schroeder E, Yewhalaw D, Nguyen J, Kebede E, Zemene E, et al. Transmission dynamics of co-endemic Plasmodium vivax and Plasmodium falciparum in Ethiopia and prevalence of anti-malarial resistant genotypes. PLoS Negl Trop Dis. 2017;11:e0005806.
Article PubMed PubMed Central Google ScholarÂ
Otienoburu SD, Maiga-Ascofare O, Schramm B, Jullien V, Jones JJ, Zolia YM, et al. Selection of Plasmodium falciparum Pfcrt and Pfmdr1 polymorphisms after treatment with artesunate–amodiaquine fixed dose combination or artemether–lumefantrine in Liberia. Malar J. 2016;15:452.
Article PubMed PubMed Central Google ScholarÂ
Asare KK, Africa J, Mbata J, Opoku YK. The emergence of chloroquine-sensitive Plasmodium falciparum is influenced by selected communities in some parts of the Central Region of Ghana. Malar J. 2021;20:447.
Article CAS PubMed PubMed Central Google ScholarÂ
Chidimatembue A, Svigel SS, Mayor A, Aide P, Nhama A, Nhamussua L, et al. Molecular surveillance for polymorphisms associated with artemisinin-based combination therapy resistance in Plasmodium falciparum isolates collected in Mozambique, 2018. Malar J. 2021;20:398.
Article CAS PubMed PubMed Central Google ScholarÂ
Malmberg M, Ferreira PE, Tarning J, Ursing J, Ngasala B, Bjorkman A, et al. Plasmodium falciparum drug resistance phenotype as assessed by patient antimalarial drug levels and its association with Pfmdr1 polymorphisms. J Infect Dis. 2013;207:842–7.
Article CAS PubMed Google ScholarÂ
Venkatesan M, Gadalla NB, Stepniewska K, Dahal P, Nsanzabana C, Price RN, et al. Polymorphisms in Plasmodium falciparum chloroquine resistance transporter and multidrug resistance 1 genes: parasite risk factors that affect treatment outcomes for Plasmodium falciparum malaria after artemether–lumefantrine and artesunate–amodiaquine. Am J Trop Med Hyg. 2014;91:833–43.
Article PubMed PubMed Central Google ScholarÂ
Dimbu PR, Horth R, Candido LM, Ferreira M, Caquece F. Continued low efficacy of artemether–lumefantrine in Angola in 2019. Antimicrob Agents Chemother. 2021;65:e01949-20.
Article PubMed PubMed Central Google ScholarÂ
Adamu A, Jada MS, Haruna HMS, Yakubu BO, Ibrahim MA, Balogun EO, et al. Plasmodium falciparum multidrug resistance gene-1 polymorphisms in Northern Nigeria: implications for the continued use of artemether–lumefantrine in the region. Malar J. 2020;19:439.
Article CAS PubMed PubMed Central Google ScholarÂ
This work is supported by Human, Heredity and Health in Africa (H3Africa) [H3A-18-002]. H3Africa is managed by Science for Africa Foundation (SFA Foundation) in partnership with Wellcome, NIH and AfSHG. We would like to thank participants of this study and laboratory technologists who were involved in collection and examining blood specimen collected from patients attending health facilities. We gratefully acknowledge Oromia Health Bureau for their consent and organizing selected health facilities to collect the data. We would like to extend our gratitude to Aklilu Lemma Institute of Pathobiology who allowed us to use laboratory facilities to perform PCR-RFLP experiments. We also thank Ms. Yergalem G/hiwot and Roza Yoseph for their assistance in laboratory work.
This work was partially supported by the PAMGEN H3 Africa project (H3A/002/18) through Addis Ababa University, Aklilu Lemma Institute of Pathobiology.
Department of Applied Biology, School of Applied Natural Science, Adama Science and Technology University, P. O. Box 1888, Adama, Ethiopia
Jifar Hassen & Hunduma Dinka
Research and Community Service Center, College of Health Science, Defense University, P. O. Box 1419, Bishoftu, Ethiopia
Aklilu Lemma Institute of Pathobiology, Addis Ababa University, P. O. Box 1176, Addis Ababa, Ethiopia
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
JH designed the study, performed experiments, analysed data, drafted and wrote draft manuscript. LG and GS organized molecular laboratory platform. GS reviewed the manuscript. HD supervised the work and reviewed the manuscript. LG conceived the idea, designed the study, supervised the overall work, participated in the interpretation of results and reviewed the manuscript. All authors read and approved the final version of the manuscript.
The study was reviewed and approved by Adama Science and Technology University research and ethical committee (certificate Reference No: RECSoANS/BIO/02/2019). Permission latter was also obtained from Oromia Regional State Health Bureau.
Authors declared that they do not have competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Primers sequences used in nested PCR for the amplification of Pfcrt and Pfmdr1 genes.
Representative gel image showing digested nested PCR products. (A): K76T codon cleaved with APoI restriction enzyme run in 2.0% agarose gel. L = 100 bp DNA ladder marker, lanes 1–9 are samples, lane 10 is uncut control and lane 11 is wild-type control (Pf3D7). (B): N86Y codon cleaved with APoI restriction enzyme run in 2.5% agarose gel. L = 100 bp ladder, Lanes 1–9 are samples, lane 10 is wild-type pos. cont. (3D7), lane 11 is mutant- type pos. cont. (K1), and lane 12 uncut control.
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.
Hassen, J., Alemayehu, G.S., Dinka, H. et al. High prevalence of Pfcrt 76T and Pfmdr1 N86 genotypes in malaria infected patients attending health facilities in East Shewa zone, Oromia Regional State, Ethiopia. Malar J 21, 286 (2022). https://doi.org/10.1186/s12936-022-04304-5
DOI: https://doi.org/10.1186/s12936-022-04304-5
Anyone you share the following link with will be able to read this content:
Sorry, a shareable link is not currently available for this article.
Provided by the Springer Nature SharedIt content-sharing initiative
By using this website, you agree to our Terms and Conditions, California Privacy Statement, Privacy statement and Cookies policy. Manage cookies/Do not sell my data we use in the preference centre.
© 2022 BioMed Central Ltd unless otherwise stated. Part of Springer Nature.