Regular monitoring of P. falciparum resistance molecular markers is a key factor in the early detection of emerging resistance and achieving malaria elimination. In this study, pfmdr1 and pfk13 genes were studied for their implication in P. falciparum resistance to antimalarial drugs, particularly amodiaquine, lumefantrine, and artemisinin derivatives. The limitation of this study was the lack of data to investigate the temporal comparison in the same site to identify subnational trends in prevalence in different geographic areas.
A selection for pfmdr 86N wild-type alleles was observed as a result of using AS + AQ and AL for many years. This selection is commonly associated with lumefantrine tolerance (Raman et al. 2019). Indeed, pfmdr1 polymorphism is usually studied because of its involvement in P. falciparum’s resistance to many antimalarial drugs, particularly those used as partner drugs of artemisinin derivatives in AS + AQ and AL (amodiaquine and lumefantrine). Single-nucleotide polymorphisms (SNPs) at codons 86 (N → Y), 184 (Y → F), and 1246 (D → Y) of pfmdr1 are among the most studied. Monitoring of NFD and YYY haplotypes is useful for detecting an early decrease in susceptibility to lumefantrine and amodiaquine, respectively (Venkatesan et al. 2014). Usually, after the change from chloroquine to ACTs, an increase in parasites carrying mdr86N was found, as reported in this study from 2013 to 2016 (Okell et al. 2018). In the same way, a previous study conducted in 2006 i.e. one year before the official adoption of ACTs as first-line treatment for uncomplicated malaria cases showed a higher value of the mdr186Y allele (47.5%) than ours (unpublished data). In addition, an upward trend of mdr184F mutant-allele prevalence is usually reported after the adoption of ACTs, as observed in the present study from 2013 to 2016 (Okell et al. 2018). The prevalence of this mutation in Abidjan in 2010 was 57% (Trebissou et al. 2014), which was lower than ours. In this study, we did not investigate the mdrD1246Y mutation. A previous study carried out in Côte d'Ivoire on this mutation reported a low incidence of mutant-types (1.6%) [unpublished data]. A scarcity of mutant mdr1246Y has also been found in Burkina Faso (Sondo et al. 2016a). These results point to lumefantrine tolerance in strains circulating through the country, as found in many reports from sub-Saharan African areas that use AL combination (Raman et al. 2019; Sondo et al. 2016; Dama et al. 2017). The pfmdr1 gene results observed in this study could be due to frequent use of AL combination over AS + AQ due to patient complaints following the use of the latter (Azagoh-Kouadio et al. 2017). Moreover, ex vivo susceptibility tests of P. falciparum strains circulating in Côte d’Ivoire should confirm this result. In Mali, a reduced ex vivo susceptibility of P. falciparum was observed after AL treatment (Dama et al. 2017).
Mutations in the pfk13 gene (C580Y, R561H, R539T, I543T, Y493H, M476I, N458Y, Y493H, and recently F446I) (Wang et al. 2018) are molecular markers of artemisinin resistance validated by in vitro and in vivo studies, while P441L, G449A, G538V, P553L, R561H, V568G, P574L, A578S, and A675V are candidate markers (WHO 2017). These validated mutations have been reported in Cambodia and other SEA countries and serve as tools for monitoring ACTs resistance (Ariey et al. 2014). These mutations are absent in this study as well as in a previous study conducted in 2017 in Abidjan (Basco et al. 2017). This result suggests that P. falciparum parasites circulating in Côte d’Ivoire are still sensitive to artemisinin. This is consistent with the findings of previous results from clinical trials at sentinel sites in the country that have shown adequate clinical and parasitological responses of up to 100% (Toure et al. 2014; Yavo et al. 2015; Konaté et al. 2018a). However, recent studies in Rwanda and Tanzania have showed the presence of R561H mutations (Uwimana et al. 2020; Bwire et al. 2020; Bergmann et al. 2021). In addition, phylogenetic analysis of Rwanda strains revealed the expansion of an indigenous R561H lineage and evidence for the de novo emergence of Pfkelch13-mediated artemisinin resistance in Rwanda (Uwimana et al. 2020, 2021). These data highlight the importance of surveilling antimalarial molecular markers in Côte d’Ivoire. Additional mutations were also identified in this study. Indeed, A557S, a non-synonymous mutation already reported in Côte d'Ivoire (Kamau et al. 2014) and Congo by Taylor et al. (Taylor et al. 2015), was found in this study. In addition, the S522C mutation observed in this study has already been described in Africa, particularly West Africa (Togo), Central Africa (Central African Republic, Gabon, Democratic Republic of Congo), and East Africa (Kenya) (Menard et al. 2016). The A578S mutation, a candidate marker for artemisinin resistance (WHO 2017), has the highest prevalence in this study (0.5%). This mutation seems to be the most widespread worldwide. Indeed, it has been reported in a small proportion in Cambodia (Straimer et al. 2015a) and several sub-Saharan African countries, including the Democratic Republic of Congo (Kamau et al. 2014; Mayengue et al. 2018), Uganda (Conrad et al. 2014), Equatorial Guinea (Li et al. 2016), Gabon (Kamau et al. 2014), Ghana (Kamau et al. 2014), Kenya (Kamau et al. 2014), Togo (Dorkenoo et al. 2016), and Mali (Kamau et al. 2014; Dama et al. 2017). In vitro studies have confirmed that this mutation is not responsible for artemisinin resistance (Straimer et al. 2015; de Laurent et al. 2018; Ashley et al. 2014). However, because of its common occurrence in many countries, further characterization, as well as an assessment of its role in in vivo parasite clearance in sub-Saharan African countries is required. Moreover, the mutations P475H, Y482H, G665S, and A676S reported in our study have not been described before. These results indicate that in Africa, parasites show a high polymorphism in pfk13 gene; however, the mutations observed are different from those reported in SEA. These differences seem to be affected by demographic and epidemiological parameters (Joy et al. 2003) rather than differential selective pressures that may be transient, and not necessarily due to artemisinin pressure (de Laurent et al. 2018). Although some candidate marker mutations have been reported in Africa, they were not involved in any case of treatment failure. Thus, the molecular mechanisms underlying artemisinin resistance remain unclear (Takala-Harrison et al. 2013). This phenomenon is likely to involve multiple genetic loci (Ariey et al. 2014b; Miotto et al. 2015), and requires complete genotyping of circulating P. falciparum strains to detect new markers. Further studies, such as site-specific genome editing experiments using zinc finger nucleases (Mayengue et al. 2018), or the CRISPR-Cas9 system (Ghorbal et al. 2014), could shed light on the role of mutations discovered in sub-Saharan Africa, particularly in Côte d'Ivoire, in the resistance of P. falciparum to artemisinin derivatives.