Sub-Saharan Africa is courting the risk of artemisinin resistance (ARTr) emerging in Plasmodium falciparum malaria parasites. Current molecular surveillance efforts for ARTr have been built on the utility of P. falciparum kelch13 (pfk13) validated molecular markers. However, whether these molecular markers will serve the purpose of early detection of artemisinin-resistant parasites in Ghana is hinged on a pfk13 dependent evolution. Here, we tested the hypothesis that the background pfk13 genome may be present before the pfk13 ARTr-conferring variant(s) is selected and that signatures of balancing selection on these genomic loci may serve as an early warning signal of ARTr. We analyzed 12 198 single nucleotide polymorphisms (SNPs) in Ghanaian clinical isolates in the Pf3K MalariaGEN dataset that passed a stringent filtering regimen. We identified signatures of balancing selection in 2 genes (phosphatidylinositol 4-kinase and chloroquine resistance transporter) previously reported as background loci for ARTr. These genes showed statistically significant and high positive values for Tajima’s D, Fu and Li’s F, and Fu and Li’s D. This indicates that the biodiversity required to establish a pfk13 background genome may have been primed in clinical isolates of P. falciparum from Ghana as of 2010. Despite the absence of ARTr in Ghana to date, our finding supports the current use of pfk13 for molecular surveillance of ARTr in Ghana and highlights the potential utility of monitoring malaria parasite populations for balancing selection in ARTr precursor background genes as early warning molecular signatures for the emergence of ARTr.

The molecular determinants of Plasmodium falciparum artemisinin resistance are the single nucleotide polymorphisms in the parasite’s kelch propeller domain, pfk13. Validated and candidate markers are under surveillance in malaria endemic countries using artemisinin-based combination therapy. However, pfk13 mutations which may confer parasite artemisinin resistance in Africa remains elusive. It has therefore become imperative to report all observed pfk13 gene polymorphisms in malaria therapeutic efficacy studies for functional characterization. We herein report all novel pfk13 mutations observed only in the Ghanaian parasite population. In all, 977 archived samples from children aged 12 years and below with uncomplicated malaria from 2007 to 2017 were used. PCR/Sanger sequencing analysis revealed 78% (763/977) of the samples analyzed were wild type (WT) for pfk13 gene. Of the 214 (22%) mutants, 78 were novel mutations observed only in Ghana. The novel SNPs include R404G, P413H, N458D/H/I, C473W/S, R529I, M579T/Y, C580R/V, D584L, N585H/I, Q661G/L. Some of the mutations were sites and ecological zones specific. There was low nucleotide diversity and purifying selection at the pfk13 locus in Ghanaian parasite population. With increasing drug pressure and its consequent parasite resistance, documenting these mutations as baseline data is crucial for future molecular surveillance of P. falciparum resistance to artemisinin in Ghana.

Plasmodium falciparum malaria is still an important disease in sub-Saharan Africa (sSA). Great strides have been made in its control spear-headed by artemisinin (ART)-based combination therapies (ACTs). However, concerns about the imminent spread of ART-resistant (ARTr) malaria parasites to sSA threaten gains already made. Attempts to mitigate this risk have highlighted the need to discover novel P. falciparum drug targets. Therefore, studies to deepen our understanding of the biology of P. falciparum are needed. The role of extracellular vesicles (EVs) in the biology of malaria parasites is not fully understood. Recently, the ART resistance-associated transcriptional profile has been reported to involve several biological processes connected to vesicular trafficking, proteotoxic stress, erythrocyte remodelling, and mitochondrial metabolism. We explored a role for EVs in developing the P. falciparum ARTr phenotype using bulk RNA sequencing of unsynchronized parasite cultures under untreated, 0.1% dimethyl sulfoxide and 700nM dihydroartemisinin treated conditions for six hours. As pathway and gene ontology analysis is limited in its curated knowledge repertoire on EVs biogenesis in P. falciparum, we used a modular (gene set) analysis approach to explore whether an EVs biogenesis module is associated with the ARTr phenotype in P. falciparum. We first generated well-defined EVs modules of interest and used statistical tools to determine differences in their expression among the parasite and treatment conditions. Then we used gene set enrichment analysis to determine the strength of the association between each EVs module of interest and the ARTr phenotype. This transcriptome-module phenotype association study (TMPAS) represents a well-powered approach to making meaningful discoveries out of bulk gene expression data. We identified four EVs module of interest and report that one module representing gene sets with correlated expression to PF3D7_1441800 – involved with EVs biogenesis in P. falciparum - is associated with the ARTr phenotype (R539T_DHA_treated versus R539T_untreated: normalized enrichment score (NES) = 1.1830174, FDR q-value < 0.25; C580R_DHA_treated versus C580R_untreated: NES = 1.2457103, FDR q-value < 0.25). PF3D7_1441800 has been reported to reduce EVs production when knocked out in P. falciparum. Altogether, our findings suggest a role for EVs in developing ART resistance and warrant further studies interrogating this association.

Molecular markers give evidence of Plasmodium falciparum resistance to antimalarials. We determined the spatial and temporal distribution of antimalarial drug resistance markers in the P. falciparum population in Ghana using targeted amplicon deep sequencing. We contextualized our findings within the framework that duration of drug exposure and transmission intensity drive trends in molecular markers. About 87% (901/1,037) of samples collected from children ≤9 years with uncomplicated malaria from 2018 to 2023 were sequenced. There was a decreasing trend in the amodiaquine/lumefantrine pfmdr1 N F SND haplotype in the Coastal zone. For chloroquine, the pfcrt CV IET haplotype showed a decreasing trend in all zones. For sulfadoxine-pyrimethamine, the pfdhps S G KAA haplotype showed increasing trends in all zones in contrast with the pfdhps AG KAA, pfdhfr IRN , and quintuple IRN + AG KAA haplotypes. Finally, validated markers of artemisinin resistance P441L (1/709), M476I (1/709), N537I (1/709), A481V (1/709), P574L (1/709), C469Y (2/709), P553L (2/709), R561H (2/709), A578S (4/709), and A675V (5/709) were identified. Zonal differences in trends for pfmdr1 N F SND haplotype are likely due to lower amodiaquine drug pressure in the Coastal zone. Policy review of amodiaquine/lumefantrine antimalarial use might help reduce the increasing trends seen in the Forest and Savannah zones. The declining prevalence of the pfcrt mutant CV IET after two decades of chloroquine disuse in Ghana suggests the return of chloroquine-sensitive parasites. Therapeutic efficacy studies must be done to verify this observation. The opposing trends of sulfadoxine-pyrimethamine molecular markers highlight concern for its use in malaria prophylaxis in pregnancy and children. Molecular surveillance remains vital to mitigating the risk of artemisinin-resistant parasites evolving in malaria-endemic regions.