Emergence of Plasmodium falciparum strains with artemisinin partial resistance in East Africa and the Horn of Africa: is there a need to panic?

The emergence and spread of artemisinin partial resistance in East and Horn of Africa is alarming. However, artemisinin-based combination therapy (ACT) generally remains efficacious for the treatment of falciparum malaria. The emergence of partial artemisinin resistance does not currently meet the criteria to initiate change on treatment guidelines nor affect ACT routine procurement and distribution. It is high time for scientists and transitional researchers to be more critical and vigilant on further changes so that national programmes will be able to make informed decisions as well as remain alert and prepared for any change that may be required in the future.

Microsatellites reveal high polymorphism and high potential for use in anti-malarial efficacy studies in areas with different transmission intensities in mainland Tanzania.

BACKGROUND: Tanzania is currently implementing therapeutic efficacy studies (TES) in areas of varying malaria transmission intensities as per the World Health Organization (WHO) recommendations. In TES, distinguishing reinfection from recrudescence is critical for the determination of anti-malarial efficacy. Recently, the WHO recommended genotyping polymorphic coding genes, merozoite surface proteins 1 and 2 (msp1 and msp2), and replacing the glutamate-rich protein (glurp) gene with one of the highly polymorphic microsatellites in Plasmodium falciparum to adjust the efficacy of antimalarials in TES. This study assessed the polymorphisms of six neutral microsatellite markers and their potential use in TES, which is routinely performed in Tanzania. METHODS: Plasmodium falciparum samples were obtained from four TES sentinel sites, Kibaha (Pwani), Mkuzi (Tanga), Mlimba (Morogoro) and Ujiji (Kigoma), between April and September 2016. Parasite genomic DNA was extracted from dried blood spots on filter papers using commercial kits. Genotyping was done using six microsatellites (Poly-α, PfPK2, TA1, C3M69, C2M34 and M2490) by capillary method, and the data were analysed to determine the extent of their polymorphisms and genetic diversity at the four sites. RESULTS: Overall, 83 (88.3%) of the 94 samples were successfully genotyped (with positive results for ≥ 50.0% of the markers), and > 50.0% of the samples (range = 47.6-59.1%) were polyclonal, with a mean multiplicity of infection (MOI) ranging from 1.68 to 1.88 among the four sites. There was high genetic diversity but limited variability among the four sites based on mean allelic richness (RS = 7.48, range = 7.27-8.03, for an adjusted minimum sample size of 18 per site) and mean expected heterozygosity (He = 0.83, range = 0.80-0.85). Cluster analysis of haplotypes using STRUCTURE, principal component analysis, and pairwise genetic differentiation (FST) did not reveal population structure or clustering of parasites according to geographic origin. Of the six markers, Poly-α was the most polymorphic, followed by C2M34, TA1 and C3M69, while M2490 was the least polymorphic. CONCLUSION: Microsatellite genotyping revealed high polyclonality and genetic diversity but no significant population structure. Poly-α, C2M34, TA1 and C3M69 were the most polymorphic markers, and Poly-α alone or with any of the other three markers could be adopted for use in TES in Tanzania.

Temporal and spatial analysis of Plasmodium falciparum genomics reveals patterns of parasite connectivity in a low-transmission district in Southern Province, Zambia.

BACKGROUND: Understanding temporal and spatial dynamics of malaria transmission will help to inform effective interventions and strategies in regions approaching elimination. Parasite genomics are increasingly used to monitor epidemiologic trends, including assessing residual transmission across seasons and importation of malaria into these regions. METHODS: In a low and seasonal transmission setting of southern Zambia, a total of 441 Plasmodium falciparum samples collected from 8 neighbouring health centres between 2012 and 2018 were genotyped using molecular inversion probes (MIPs n = 1793) targeting a total of 1832 neutral and geographically informative SNPs distributed across the parasite genome. After filtering for quality and missingness, 302 samples and 1410 SNPs were retained and used for downstream population genomic analyses. RESULTS: The analyses revealed most (67%, n = 202) infections harboured one clone (monogenomic) with some variation at local level suggesting low, but heterogenous malaria transmission. Relatedness identity-by-descent (IBD) analysis revealed variable distribution of IBD segments across the genome and 6% of pairs were highly-related (IBD ≥ 0.25). Some of the highly-related parasite populations persisted across multiple seasons, suggesting that persistence of malaria in this low-transmission region is fueled by parasites "seeding" across the dry season. For recent years, clusters of clonal parasites were identified that were dissimilar to the general parasite population, suggesting parasite populations were increasingly fragmented at small spatial scales due to intensified control efforts. Clustering analysis using PCA and t-SNE showed a lack of substantial parasite population structure. CONCLUSION: Leveraging both genomic and epidemiological data provided comprehensive picture of fluctuations in parasite populations in this pre-elimination setting of southern Zambia over 7 years.

Performance evaluation of RDT, light microscopy, and PET-PCR for detecting Plasmodium falciparum malaria infections in the 2018 Zambia National Malaria Indicator Survey.

BACKGROUND: Zambia continues to advance on the path to elimination with significant reductions in malaria morbidity and mortality. Crucial components that have contributed to progress thus far and are necessary for achieving the national malaria elimination goals include properly identifying and treating all malaria cases through accurate diagnosis. This study sought to compare and assess the diagnostic performance of Rapid Diagnostic Tests (RDT) and Light Microscopy (LM) with photo-induced electron transfer polymerase chain reaction (PET-PCR) as the gold standard using 2018 Malaria Indicator Survey (MIS) data across Zambia to better understand diagnostic accuracy metrics and how these vary across a transmission gradient. METHODS: Cross-sectional samples collected in a nationally representative survey from 7 provinces in Zambia were tested for the presence of malaria parasites by light microscopy (LM), rapid diagnostic test (RDT) and the gold standard PET-PCR. Diagnostic performance was assessed including sensitivity, specificity, negative- and positive-predictive values across a wide malaria transmission spectrum. Diagnostic accuracy metrics were measured, and statistically significant differences were calculated between test methods for different outcome variables. RESULTS: From the individuals included in the MIS, the overall prevalence of Plasmodium falciparum malaria was 32.9% by RDT, 19.4% by LM, and 23.2% by PET-PCR. Herein, RDT and LM diagnostic performance was compared against gold standard PET-PCR with LM displaying a higher diagnostic accuracy than RDTs (91.3% vs. 84.6% respectively) across the transmission spectrum in Zambia. However, the performance of both diagnostics was significantly reduced in low parasitaemia samples. Consistent with previous studies, RDT diagnostic accuracy was predominantly affected by a high rate of false positives. CONCLUSIONS: RDTs and LM both perform well across a range of transmission intensities within their respective target applications, i.e., in the community, for the former, where ease of use and speed of result is critical, and at the health facility, for the latter, where accuracy is prioritized. However, the performance of both diagnostic methods is adversely affected by low parasitaemia infections. As Zambia moves towards elimination more sensitive tools may be required to identify the last cases.

Monitoring Plasmodium falciparum and Plasmodium vivax using microsatellite markers indicates limited changes in population structure after substantial transmission decline in Papua New Guinea.

Monitoring the genetic structure of pathogen populations may be an economical and sensitive approach to quantify the impact of control on transmission dynamics, highlighting the need for a better understanding of changes in population genetic parameters as transmission declines. Here we describe the first population genetic analysis of two major human malaria parasites, Plasmodium falciparum (Pf) and Plasmodium vivax (Pv), following nationwide distribution of long-lasting insecticide-treated nets (LLINs) in Papua New Guinea (PNG). Parasite isolates from pre- (2005-2006) and post-LLIN (2010-2014) were genotyped using microsatellite markers. Despite parasite prevalence declining substantially (East Sepik Province: Pf = 54.9%-8.5%, Pv = 35.7%-5.6%, Madang Province: Pf = 38.0%-9.0%, Pv: 31.8%-19.7%), genetically diverse and intermixing parasite populations remained. Pf diversity declined modestly post-LLIN relative to pre-LLIN (East Sepik: Rs  = 7.1-6.4, HE  = 0.77-0.71; Madang: Rs  = 8.2-6.1, HE  = 0.79-0.71). Unexpectedly, population structure present in pre-LLIN populations was lost post-LLIN, suggesting that more frequent human movement between provinces may have contributed to higher gene flow. Pv prevalence initially declined but increased again in one province, yet diversity remained high throughout the study period (East Sepik: Rs  = 11.4-9.3, HE  = 0.83-0.80; Madang: Rs  = 12.2-14.5, HE  = 0.85-0.88). Although genetic differentiation values increased between provinces over time, no significant population structure was observed at any time point. For both species, a decline in multiple infections and increasing clonal transmission and significant multilocus linkage disequilibrium post-LLIN were positive indicators of impact on the parasite population using microsatellite markers. These parameters may be useful adjuncts to traditional epidemiological tools in the early stages of transmission reduction.

SNP barcodes provide higher resolution than microsatellite markers to measure Plasmodium vivax population genetics.

BACKGROUND: Genomic surveillance of malaria parasite populations has the potential to inform control strategies and to monitor the impact of interventions. Barcodes comprising large numbers of single nucleotide polymorphism (SNP) markers are accurate and efficient genotyping tools, however may need to be tailored to specific malaria transmission settings, since 'universal' barcodes can lack resolution at the local scale. A SNP barcode was developed that captures the diversity and structure of Plasmodium vivax populations of Papua New Guinea (PNG) for research and surveillance. METHODS: Using 20 high-quality P. vivax genome sequences from PNG, a total of 178 evenly spaced neutral SNPs were selected for development of an amplicon sequencing assay combining a series of multiplex PCRs and sequencing on the Illumina MiSeq platform. For initial testing, 20 SNPs were amplified in a small number of mono- and polyclonal P. vivax infections. The full barcode was then validated by genotyping and population genetic analyses of 94 P. vivax isolates collected between 2012 and 2014 from four distinct catchment areas on the highly endemic north coast of PNG. Diversity and population structure determined from the SNP barcode data was then benchmarked against that of ten microsatellite markers used in previous population genetics studies. RESULTS: From a total of 28,934,460 reads generated from the MiSeq Illumina run, 87% mapped to the PvSalI reference genome with deep coverage (median = 563, range 56-7586) per locus across genotyped samples. Of 178 SNPs assayed, 146 produced high-quality genotypes (minimum coverage = 56X) in more than 85% of P. vivax isolates. No amplification bias was introduced due to either polyclonal infection or whole genome amplification (WGA) of samples before genotyping. Compared to the microsatellite panels, the SNP barcode revealed greater variability in genetic diversity between populations and geographical population structure. The SNP barcode also enabled assignment of genotypes according to their geographic origins with a significant association between genetic distance and geographic distance at the sub-provincial level. CONCLUSIONS: High-throughput SNP barcoding can be used to map variation of malaria transmission dynamics at sub-national resolution. The low cost per sample and genotyping strategy makes the transfer of this technology to field settings highly feasible.

A Worldwide Map of Plasmodium falciparum K13-Propeller Polymorphisms.

BACKGROUND: Recent gains in reducing the global burden of malaria are threatened by the emergence of Plasmodium falciparum resistance to artemisinins. The discovery that mutations in portions of a P. falciparum gene encoding kelch (K13)-propeller domains are the major determinant of resistance has provided opportunities for monitoring such resistance on a global scale. METHODS: We analyzed the K13-propeller sequence polymorphism in 14,037 samples collected in 59 countries in which malaria is endemic. Most of the samples (84.5%) were obtained from patients who were treated at sentinel sites used for nationwide surveillance of antimalarial resistance. We evaluated the emergence and dissemination of mutations by haplotyping neighboring loci. RESULTS: We identified 108 nonsynonymous K13 mutations, which showed marked geographic disparity in their frequency and distribution. In Asia, 36.5% of the K13 mutations were distributed within two areas--one in Cambodia, Vietnam, and Laos and the other in western Thailand, Myanmar, and China--with no overlap. In Africa, we observed a broad array of rare nonsynonymous mutations that were not associated with delayed parasite clearance. The gene-edited Dd2 transgenic line with the A578S mutation, which expresses the most frequently observed African allele, was found to be susceptible to artemisinin in vitro on a ring-stage survival assay. CONCLUSIONS: No evidence of artemisinin resistance was found outside Southeast Asia and China, where resistance-associated K13 mutations were confined. The common African A578S allele was not associated with clinical or in vitro resistance to artemisinin, and many African mutations appear to be neutral. (Funded by Institut Pasteur Paris and others.).

GnRH and prostaglandin-based synchronization protocols as alternatives to progestogen-based treatments in sheep.

The study investigated, for cycling sheep, synchronizing protocols simultaneously to the standard "P" protocol using progestogens priming with intravaginal devices and gonadotropin. In November 2014, 90 adult Menz ewes were assigned to either the "P" protocol, "PGF" treatment where oestrus and ovulation were synchronized using two injections of prostaglandin 11 days apart or a "GnRH" treatment where the ewes had their oestrus and ovulation synchronized with GnRH (day 0)-prostaglandin (day 6)-GnRH (day 9) sequence. The ewes were naturally mated at the induced oestrus and the following 36 days. Plasma progesterone revealed that 92% of the ewes were ovulating before synchronization and all, except one, ovulated in response to the applied treatments. All "P" ewes exhibited oestrus during the 96-hr period after the end of the treatments in comparison with only 79.3% and 73.3% for "PGF" and "GnRH" ewes, respectively (p < .05). Onset and duration of oestrus were affected by the hormonal treatment (p < .05); "GnRH" ewes showed oestrus earliest and had the shortest oestrous duration. Lambing rate from mating at the induced oestrus was lower for "P" than for "PGF" ewes (55.6% and 79.3%, respectively; p < .05). The same trait was also lower for "P" than for "PGF" and "GnRH" ewes (70.4%, 89.7% and 86.7%, respectively; p < .05) following the 36-day mating period. Prostaglandin and GnRH analogue-based protocols are promising alternatives for both controlled natural mating and fixed insemination of Menz sheep after the rainy season when most animals are spontaneously cycling.

Diagnostic performance of rapid diagnostic tests for the diagnosis of malaria at public health facilities in north-west Ethiopia.

OBJECTIVE: To assess the performance of RDTs against nested polymerase chain reaction (nPCR) for the diagnosis of malaria in public health facilities in north-western Ethiopia. METHODS: Cross-sectional study at public health facilities in North Gondar, Ethiopia, of 359 febrile patients with signs and symptoms consistent with malaria. Finger prick blood samples were collected for testing in a P. falciparum/pan-malaria RDTs and for molecular analysis. Sensitivity, specificity and predictive values were determined for the RDTs using nPCR as reference diagnostic method. Kappa value was determined to demonstrate the consistency of the results between the diagnostic tools. RESULTS: By RDTs, 22.28% (80/359) of patients tested positive for malaria, and by nPCR, 27.02% (97/359) did. In nPCR, 1.67% (6/359) and 0.28% (1/359) samples were positive for P. ovale and P. malariae, which had almost all tested negative in the RDTs. The sensitivity, specificity, positive and negative predictive values of RDTs for the diagnosis of malaria were 62.9%, 92.7%, 76.3% and 87.1%, respectively, with 0.589 measurement agreement between RDTs and nPCR. The sensitivity and specificity of RDTs for P. falciparum identification only were 70.8% and 95.2%, and 65.2% and 93.1% for P. vivax. CONCLUSION: Although RDTs are commonly used at health posts in resource-limited environments, their sensitivity and specificity for the detection and species identification of Plasmodium parasites were poor compared to nPCR, suggesting caution in interpreting RDTs results. Particularly, in the light of expanded efforts to eliminate malaria in the country, more sensitive diagnostic procedures will be needed.

Evolutionary relationships of the critically endangered frog Ericabatrachus baleensis Largen, 1991 with notes on incorporating previously unsampled taxa into large-scale phylogenetic analyses.

BACKGROUND: The phylogenetic relationships of many taxa remain poorly known because of a lack of appropriate data and/or analyses. Despite substantial recent advances, amphibian phylogeny remains poorly resolved in many instances. The phylogenetic relationships of the Ethiopian endemic monotypic genus Ericabatrachus has been addressed thus far only with phenotypic data and remains contentious. RESULTS: We obtained fresh samples of the now rare and Critically Endangered Ericabatrachus baleensis and generated DNA sequences for two mitochondrial and four nuclear genes. Analyses of these new data using de novo and constrained-tree phylogenetic reconstructions strongly support a close relationship between Ericabatrachus and Petropedetes, and allow us to reject previously proposed alternative hypotheses of a close relationship with cacosternines or Phrynobatrachus. CONCLUSIONS: We discuss the implications of our results for the taxonomy, biogeography and conservation of E. baleensis, and suggest a two-tiered approach to the inclusion and analyses of new data in order to assess the phylogenetic relationships of previously unsampled taxa. Such approaches will be important in the future given the increasing availability of relevant mega-alignments and potential framework phylogenies.

Antimicrobial resistance trends in clinical Escherichia coli and Klebsiella pneumoniae in Ethiopia.

Background: Clinicians rely on local antimicrobial resistance pattern data to guide empiric treatment for seriously ill patients when culture and antimicrobial susceptibility testing results are not immediately available. Objective: This study aimed to analyse 5-year trends in antimicrobial resistance profiles of Escherichia coli and Klebsiella pneumoniae isolates. Methods: Bacteriology reports from 2017 to 2021 at the Ethiopian Public Health Institute were analysed retrospectively. Isolates were identified using either the VITEK 2 Compact system, the BD Phoenix M50 instrument, or conventional biochemical tests. Antimicrobial susceptibility testing was conducted using either the Kirby-Bauer disk diffusion method or the VITEK 2 Compact system and BD Phoenix M50 systems available at the time of testing. The Cochran Armitage trend test was employed to test the significance of antimicrobial resistance trends over time. P-values less than 0.05 were considered statistically significant. Results: Of the 5382 bacteriology reports examined, 458 (9%) were on E. coli and 266 (5%) were on K. pneumoniae. Both K. pneumoniae (88%) and E. coli (65%) demonstrated high resistance to extended-spectrum cephalosporins. However, both K. pneumoniae (14%) and E. coli (5%) showed lower rates of resistance to carbapenems compared to other antimicrobials. In K. pneumoniae, resistance to carbapenems (from 0% to 38%; p < 0.001) and ciprofloxacin (from 41% to 90%; p < 0.001) increased significantly between 2017 and 2021. Conclusion: Both organisms showed very high resistance to broad-spectrum antibiotics. Additionally, K. pneumoniae demonstrated a statistically significant rise in ciprofloxacin and carbapenem resistance. What this study adds: This study emphasises the significance of regular reporting of local antimicrobial resistance patterns as this information can guide appropriate empiric therapy and efforts to address antimicrobial resistance issues.

Genomics reveals heterogeneous Plasmodium falciparum transmission and selection signals in Zambia.

BACKGROUND: Genomic surveillance is crucial for monitoring malaria transmission and understanding parasite adaptation to interventions. Zambia lacks prior nationwide efforts in malaria genomic surveillance among African countries. METHODS: We conducted genomic surveillance of Plasmodium falciparum parasites from the 2018 Malaria Indicator Survey in Zambia, a nationally representative household survey of children under five years of age. We whole-genome sequenced and analyzed 241 P. falciparum genomes from regions with varying levels of malaria transmission across Zambia and estimated genetic metrics that are informative about transmission intensity, genetic relatedness between parasites, and selection. RESULTS: We provide genomic evidence of widespread within-host polygenomic infections, regardless of epidemiological characteristics, underscoring the extensive and ongoing endemic malaria transmission in Zambia. Our analysis reveals country-level clustering of parasites from Zambia and neighboring regions, with distinct separation in West Africa. Within Zambia, identity by descent (IBD) relatedness analysis uncovers local spatial clustering and rare cases of long-distance sharing of closely related parasite pairs. Genomic regions with large shared IBD segments and strong positive selection signatures implicate genes involved in sulfadoxine-pyrimethamine and artemisinin combination therapies drug resistance, but no signature related to chloroquine resistance. Furthermore, differences in selection signatures, including drug resistance loci, are observed between eastern and western Zambian parasite populations, suggesting variable transmission intensity and ongoing drug pressure. CONCLUSIONS: Our findings enhance our understanding of nationwide P. falciparum transmission in Zambia, establishing a baseline for analyzing parasite genetic metrics as they vary over time and space. These insights highlight the urgency of strengthening malaria control programs and surveillance of antimalarial drug resistance.

Malaria is caused by a parasite that is spread to humans via mosquito bites. It is a leading cause of death in children under five years old in sub-Saharan Africa. Analysis of the malaria parasite's complete set of DNA (its genome) can help us to understand transmission of the disease and how this changes in response to different strategies to control the disease. We analyzed the genomes of malaria parasites from children across Zambia. Our study revealed that 77% of children harbored multiple parasite strains, which suggests that local transmission (transmission between people within the same local area) is high. Genetic evidence for long-distance transmission was rarer. Furthermore, our findings suggest parasites are evolving in response to antimalarial drugs. Our study enhances our understanding of malaria dynamics in Zambia and may help to inform strategies for improved surveillance and control.

High-Throughput Genotyping of Plasmodium vivax in the Peruvian Amazon via Molecular Inversion Probes.

Plasmodium vivax transmission occurs throughout the tropics and is an emerging threat in areas of Plasmodium falciparum decline, causing relapse infections that complicate treatment and control. Targeted sequencing for P. falciparum has been widely deployed to detect population structure and the geographic spread of antimalarial and diagnostic resistance. However, there are fewer such tools for P. vivax . Leveraging global variation data, we designed four molecular inversion probe (MIP) genotyping panels targeting geographically differentiating SNPs, neutral SNPs, putative antimalarial resistance genes, and vaccine candidate genes. We deployed these MIP panels on 866 infections from the Peruvian Amazon and identified transmission networks with clonality (IBD>0.99), copy number variation in Pvdbp and multiple Pvrbps , fixation of putative antimalarial resistance, and balancing selection in 13 vaccine candidate genes. Our MIP panels are the broadest genotyping panel currently available and are poised for successful deployment in other regions of P. vivax transmission.

Ex vivo susceptibility to antimalarial drugs and polymorphisms in drug resistance genes of African Plasmodium falciparum, 2016-2023: a genotype-phenotype association study.

Background: Given the altered responses to both artemisinins and lumefantrine in Eastern Africa, monitoring antimalarial drug resistance in all African countries is paramount. Methods: We measured the susceptibility to six antimalarials using ex vivo growth inhibition assays (IC50) for a total of 805 Plasmodium falciparum isolates obtained from travelers returning to France (2016-2023), mainly from West and Central Africa. Isolates were sequenced using molecular inversion probes (MIPs) targeting fourteen drug resistance genes across the parasite genome. Findings: Ex vivo susceptibility to several drugs has significantly decreased in 2019-2023 versus 2016-2018 parasite samples: lumefantrine (median IC50: 23·0 nM [IQR: 14·4-35·1] in 2019-2023 versus 13·9 nM [8·42-21·7] in 2016-2018, p<0·0001), monodesethylamodiaquine (35·4 [21·2-51·1] versus 20·3 nM [15·4-33·1], p<0·0001), and marginally piperaquine (20·5 [16·5-26·2] versus 18.0 [14·2-22·4] nM, p<0·0001). Only four isolates carried a validated pfkelch13 mutation. Multiple mutations in pfcrt and one in pfmdr1 (N86Y) were significantly associated with altered susceptibility to multiple drugs. The susceptibility to lumefantrine was altered by pfcrt and pfmdr1 mutations in an additive manner, with the wild-type haplotype (pfcrt K76-pfmdr1 N86) exhibiting the least susceptibility. Interpretation: Our study on P. falciparum isolates from West and Central Africa indicates a low prevalence of molecular markers of artemisinin resistance but a significant decrease in susceptibility to the partner drugs that have been the most widely used since a decade -lumefantrine and amodiaquine. These phenotypic changes likely mark parasite adaptation to sustained drug pressure and call for intensifying the monitoring of antimalarial drug resistance in Africa. Funding: This work was supported by the French Ministry of Health (grant to the French National Malaria Reference Center) and by the Agence Nationale de la Recherche (ANR-17-CE15-0013-03 to JC). JAB was supported by NIH R01AI139520. JR postdoctoral fellowship was funded by Institut de Recherche pour le Développement.

National genomic profiling of Plasmodium falciparum antimalarial resistance in Zambian children participating in the 2018 Malaria Indicator Survey.

The emergence of antimalarial drug resistance is a major threat to malaria control and elimination. Using whole genome sequencing of 282 P. falciparum samples collected during the 2018 Zambia National Malaria Indicator Survey, we determined the prevalence and spatial distribution of known and candidate antimalarial drug resistance mutations. High levels of genotypic resistance were found across Zambia to pyrimethamine, with over 94% (n=266) of samples having the Pfdhfr triple mutant (N51 I , C59 R , and S108 N ), and sulfadoxine, with over 84% (n=238) having the Pfdhps double mutant (A437 G and K540 E ). In northern Zambia, 5.3% (n=15) of samples also harbored the Pfdhps A581 G mutation. Although 29 mutations were identified in Pfkelch13 , these mutations were present at low frequency (<2.5%), and only three were WHO-validated artemisinin partial resistance mutations: P441 L (n=1, 0.35%), V568 M (n=2, 0.7%) and R622 T (n=1, 0.35%). Notably, 91 (32%) of samples carried the E431 K mutation in the Pfatpase6 gene, which is associated with artemisinin resistance. No specimens carried any known mutations associated with chloroquine resistance in the Pfcrt gene (codons 72-76). P. falciparum strains circulating in Zambia were highly resistant to sulfadoxine and pyrimethamine but remained susceptible to chloroquine and artemisinin. Despite this encouraging finding, early genetic signs of developing artemisinin resistance highlight the urgent need for continued vigilance and expanded routine genomic surveillance to monitor these changes.

Temporal genomic analysis of Plasmodium falciparum reveals increased prevalence of mutations associated with delayed clearance following treatment with artemisinin-lumefantrine in Choma District, Southern Province, Zambia.

The emergence of antimalarial drug resistance is an impediment to malaria control and elimination in Africa. Analysis of temporal trends in molecular markers of resistance is critical to inform policy makers and guide malaria treatment guidelines. In a low and seasonal transmission region of southern Zambia, we successfully genotyped 85.5% (389/455) of Plasmodium falciparum samples collected between 2013-2018 from 8 spatially clustered health centres using molecular inversion probes (MIPs) targeting key drug resistance genes. Aside from one sample carrying K13 R622I, none of the isolates carried other World Health Organization-validated or candidate artemisinin partial resistance (ART-R) mutations in K13. However, 13% (CI, 9.6-17.2) of isolates had the AP2MU S160N mutation, which has been associated with delayed clearance following artemisinin combination therapy in Africa. This mutation increased in prevalence between 2015-2018 and bears a genomic signature of selection. During this time period, there was an increase in the MDR1 NFD haplotype that is associated with reduced susceptibility to lumefantrine. Sulfadoxine-pyrimethamine polymorphisms were near fixation. While validated ART-R mutations are rare, a mutation associated with slow parasite clearance in Africa appears to be under selection in southern Zambia.

Genomics reveals heterogeneous Plasmodium falciparum transmission and population differentiation in Zambia and bordering countries.

Genomic surveillance plays a critical role in monitoring malaria transmission and understanding how the parasite adapts in response to interventions. We conducted genomic surveillance of malaria by sequencing 241 Plasmodium falciparum genomes from regions with varying levels of malaria transmission across Zambia. We found genomic evidence of high levels of within-host polygenomic infections, regardless of epidemiological characteristics, underscoring the extensive and ongoing endemic malaria transmission in the country. We identified country-level clustering of parasites from Zambia and neighboring countries, and distinct clustering of parasites from West Africa. Within Zambia, our identity by descent (IBD) relatedness analysis uncovered spatial clustering of closely related parasite pairs at the local level and rare cases of long-distance sharing. Genomic regions with large shared IBD segments and strong positive selection signatures identified genes involved in sulfadoxine-pyrimethamine and artemisinin combination therapies drug resistance, but no signature related to chloroquine resistance. Together, our findings enhance our understanding of P. falciparum transmission nationwide in Zambia and highlight the urgency of strengthening malaria control programs and surveillance of antimalarial drug resistance.

Evidence of artemisinin partial resistance in North-western Tanzania: clinical and drug resistance markers study.

Background: Artemisinin-based combination therapies (ACTs) are the recommended antimalarial drugs for the treatment of uncomplicated malaria. The recent emergence of artemisinin partial resistance (ART-R) in Rwanda, Uganda and Eritrea is of great concern. In Tanzania, a nationwide molecular malaria surveillance in 2021 showed a high prevalence of the Kelch13 (K13) 561H mutation in Plasmodium falciparum from the north-western region, close to the border with Rwanda and Uganda. This study was conducted in 2022 to evaluate the efficacy of artemether-lumefantrine (AL) and artesunate-amodiaquine (ASAQ) for the treatment of uncomplicated falciparum malaria and to confirm the presence of ART-R in Tanzania. Methods: This single-arm study evaluated the efficacy of AL and ASAQ in eligible children aged six months to 10 years at Bukangara Dispensary in Karagwe District, Kagera Region. Clinical and parasitological responses were monitored for 28 days according to standard WHO protocol. Mutations in K13 gene and extended haplotypes with these mutations were analysed using Sanger and whole genome sequencing data, respectively. Findings: 176 children (88 in each AL and ASAQ group) were enrolled and all achieved the defined outcomes. PCR-corrected adequate clinical and parasitological response (ACPR) was 98.3% (95% CI: 90.8-100) and 100.0% (95% CI: 95.8-100) for AL and ASAQ, respectively. Parasitaemia on day 3 was observed in 11/88 (12.5%) and 17/88 (19.3%) in the AL and ASAQ groups, respectively. The half-life of parasitaemia was significantly higher (>6.5 hrs) in patients with parasitaemia on day 3 and/or mutations in K13 gene at enrolment. Most patients with parasitaemia on day 3 (8/11 = 72.7% in the AL group and 10/17 = 58.8% in the ASAQ group) had 561H mutation at enrolment. The parasites with K13 mutations were not similar to those from south-east Asia and Rwanda, but had the same core haplotype of a new 561H haplotype reported in Kagera in 2021. Interpretation: These findings confirm the presence of ART-R in Tanzania. A context-specific strategy to respond to artemisinin partial resistance is urgently needed. Although both AL and ASAQ showed high efficacy, increased vigilance for reduced efficacy of these ACTs and detection of ART-R in other parts of the country is critical.

Ex vivo susceptibility to antimalarial drugs and polymorphisms in drug resistance genes of African Plasmodium falciparum, 2016-2023: a genotype-phenotype association study.

Background: Given the altered responses to both artemisinins and lumefantrine in Eastern Africa, monitoring antimalarial drug resistance in all African countries is paramount. Methods: We measured the susceptibility to six antimalarials using ex vivo growth inhibition assays (IC 50 ) for a total of 805 Plasmodium falciparum isolates obtained from travelers returning to France (2016-2023), mainly from West and Central Africa. Isolates were sequenced using molecular inversion probes (MIPs) targeting fourteen drug resistance genes across the parasite genome. Findings: Ex vivo susceptibility to several drugs has significantly decreased in 2019-2023 versus 2016-2018 parasite samples: lumefantrine (median IC 50 : 23·0 nM [IQR: 14·4-35·1] in 2019-2023 versus 13·9 nM [8·42-21·7] in 2016-2018, p<0·0001), monodesethylamodiaquine (35·4 [21·2-51·1] versus 20·3 nM [15·4-33·1], p<0·0001), and marginally piperaquine (20·5 [16·5-26·2] versus 18.0 [14·2-22·4] nM, p<0·0001). Only four isolates carried a validated pfkelch13 mutation. Multiple mutations in pfcrt and one in pfmdr1 (N86Y) were significantly associated with altered susceptibility to multiple drugs. The susceptibility to lumefantrine was altered by pfcrt and pfmdr1 mutations in an additive manner, with the wild-type haplotype ( pfcrt K76- pfmdr1 N86) exhibiting the least susceptibility. Interpretation: Our study on P. falciparum isolates from West and Central Africa indicates a low prevalence of molecular markers of artemisinin resistance but a significant decrease in susceptibility to the partner drugs that have been the most widely used since a decade -lumefantrine and amodiaquine. These phenotypic changes likely mark parasite adaptation to sustained drug pressure and call for intensifying the monitoring of antimalarial drug resistance in Africa.

Prevalence of mutations associated with artemisinin partial resistance and sulfadoxine-pyrimethamine resistance in 13 regions in Tanzania in 2021: a cross-sectional survey.

BACKGROUND: The emergence of the artemisinin partial resistance (ART-R) mutation in the Plasmodium falciparum kelch13 gene (k13), Arg561His, in Rwanda and the regional presence of polymorphisms affecting sulfadoxine-pyrimethamine have raised concern in neighbouring Tanzania. The goal of this study was to assess the status of antimalarial resistance in Tanzania, with a focus on the border with Rwanda, to understand the distribution of the Arg561His mutation, partner drug resistance, and resistance to chemoprevention drugs. METHODS: In this cross-sectional survey, capillary dried blood spots were collected from malaria positive asymptomatic individuals in the community and symptomatic individuals in health facilities aged 6 months and older, in 13 regions of mainland Tanzania from Jan 31 to June 26, 2021. Exclusion criteria included residence of the areas other than the target sites, presenting to the health facility for care and treatment of conditions other than malaria, and not providing informed consent. Samples were assessed for antimalarial resistance polymorphisms and genetic relatedness using molecular inversion probes targeting P falciparum and short-read whole-genome sequencing. The primary outcome was the prevalence of molecular markers of antimalarial resistance at the region level, as well as at the district level in Kagera, a region in the northwest of the country at the border with Rwanda. FINDINGS: 6855 (88·1%) of 7782 capillary dried blood spot samples collected were successfully genotyped. The overall prevalence of k13 Arg561His in Kagera was 7·7% (90% CI 6·0-9·4; 50 of 649), with the highest prevalence in the districts near the Rwandan border (22·8% [31 of 136] in Karagwe, 14·4% [17 of 118]) in Kyerwa, and 1·4% [two of 144] in Ngara). k13 Arg561His was uncommon in the other regions. Haplotype analysis suggested that some of these parasites are related to isolates collected in Rwanda in 2015, supporting regional spread of Arg561His. However, a novel k13 Arg561His haplotype was observed, potentially indicating a second origin in the region. Other validated k13 resistance polymorphisms (one Arg622Ile and two Ala675Val isolates) were also identified. A region of prevalent dihydrofolate reductase Ile164Leu mutation, associated with sulfadoxine-pyrimethamine resistance, was also identified in Kagera (15·2% [12·6-17·8%]; 80 of 526). The mutant crt Lys76Thr mutation, associated with chloroquine and amodiaquine resistance, was uncommon, occurring only in 75 of 2861 genotyped isolates, whereases the wild-type mdr1 Asn86Tyr allele, associated with reduced sensitivity to lumefantrine, was found in 99·7% (3819 of 3830) of samples countrywide. INTERPRETATION: These findings show that the k13 Arg561His mutation is common in northwest Tanzania and that multiple emergences of ART-R, similar as to what was seen in southeast Asia, have occurred. Mutations associated with high levels of sulfadoxine-pyrimethamine resistance are common. These results raise concerns about the long-term efficacy of artemisinin and chemoprevention antimalarials in the region. Understanding how multiple emergences interact with drivers of regional spread is essential for combating ART-R in Africa. FUNDING: This study was funded by the Bill & Melinda Gates Foundation and the National Institutes of Health.