Genomic epidemiology of artemisinin resistant malaria

An overview of artemisinin-resistant malaria and associated Pfk13 gene mutations in Central Africa

Malaria caused by Plasmodium falciparum is one of the deadliest and most common tropical infectious diseases. However, the emergence of artemisinin drug resistance associated with the parasite's Pfk13 gene, threatens the public health of individual countries as well as current efforts to reduce malaria burdens globally. It is of concern that artemisinin-resistant parasites may be selected or have already emerged in Africa. This narrative review aims to evaluate the published evidence concerning validated, candidate, and novel Pfk13 polymorphisms in ten Central African countries. Results show that four validated non-synonymous polymorphisms (M476I, R539T, P553L, and P574L), directly associated with a delayed therapy response, have been reported in the region. Also, two Pfk13 polymorphisms associated to artemisinin resistance but not validated (C469F and P527H) have been reported. Furthermore, several non-validated mutations have been observed in Central Africa, and one allele A578S, is commonly found in different countries, although additional molecular and biochemical studies are needed to investigate whether those mutations alter artemisinin effects. This information is discussed in the context of biochemical and genetic aspects of Pfk13, and related to the regional malaria epidemiology of Central African countries.

A Plasmodium falciparum genetic cross reveals the contributions of pfcrt and plasmepsin II/III to piperaquine drug resistance

Piperaquine (PPQ) is widely used in combination with dihydroartemisinin as a first-line treatment against malaria. Multiple genetic drivers of PPQ resistance have been reported, including mutations in the Plasmodium falciparum chloroquine resistance transporter (pfcrt) and increased copies of plasmepsin II/III (pm2/3). We generated a cross between a Cambodia-derived multidrug-resistant KEL1/PLA1 lineage isolate (KH004) and a drug-susceptible Malawian parasite (Mal31). Mal31 harbors a wild-type (3D7-like) pfcrt allele and a single copy of pm2/3, while KH004 has a chloroquine-resistant (Dd2-like) pfcrt allele with an additional G367C substitution and multiple copies of pm2/3. We recovered 104 unique recombinant parasites and examined a targeted set of progeny representing all possible combinations of variants at pfcrt and pm2/3. We performed a detailed analysis of competitive fitness and a range of PPQ susceptibility phenotypes with these progenies, including PPQ survival assay, area under the dose response curve, and a limited point IC50. We find that inheritance of the KH004 pfcrt allele is required for reduced PPQ sensitivity, whereas copy number variation in pm2/3 further decreases susceptibility but does not confer resistance in the absence of additional mutations in pfcrt. A deep investigation of genotype-phenotype relationships demonstrates that progeny clones from experimental crosses can be used to understand the relative contributions of pfcrt, pm2/3, and parasite genetic background to a range of PPQ-related traits. Additionally, we find that the resistance phenotype associated with parasites inheriting the G367C substitution in pfcrt is consistent with previously validated PPQ resistance mutations in this transporter.IMPORTANCEResistance to piperaquine, used in combination with dihydroartemisinin, has emerged in Cambodia and threatens to spread to other malaria-endemic regions. Understanding the causal mutations of drug resistance and their impact on parasite fitness is critical for surveillance and intervention and can also reveal new avenues to limiting the evolution and spread of drug resistance. An experimental genetic cross is a powerful tool for pinpointing the genetic determinants of key drug resistance and fitness phenotypes and has the distinct advantage of quantifying the effects of naturally evolved genetic variation. Our study was strengthened since the full range of copies of KH004 pm2/3 was inherited among the progeny clones, allowing us to directly test the role of the pm2/3 copy number on resistance-related phenotypes in the context of a unique pfcrt allele. Our multigene model suggests an important role for both loci in the evolution of this multidrug-resistant parasite lineage.

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.

Impact of different mutations on Kelch13 protein levels, ART resistance, and fitness cost in Plasmodium falciparum parasites

Reduced susceptibility to ART, the first-line treatment against malaria, is common in South East Asia (SEA). It is associated with point mutations, mostly in kelch13 (k13) but also in other genes, like ubp1. K13 and its compartment neighbors (KICs), including UBP1, are involved in endocytosis of host cell cytosol. We tested 135 mutations in KICs but none conferred ART resistance. Double mutations of k13C580Y with k13R539T or k13C580Y with ubp1R3138H, did also not increase resistance. In contrast, k13C580Y parasites subjected to consecutive RSAs did, but the k13 sequence was not altered. Using isogenic parasites with different k13 mutations, we found correlations between K13 protein amount, resistance, and fitness cost. Titration of K13 and KIC7 indicated that the cellular levels of these proteins determined resistance through the rate of endocytosis. While fitness cost of k13 mutations correlated with ART resistance, ubp1R3138H caused a disproportionately higher fitness cost.

IMPORTANCE

Parasites with lowered sensitivity to artemisinin-based drugs are becoming widespread. However, even in these "resistant" parasites not all parasites survive treatment. We found that the proportion of surviving parasites correlates with the fitness cost of resistance-inducing mutations which might indicate that the growth disadvantages prevents resistance levels where all parasites survive treatment. We also found that combining two common resistance mutations did not increase resistance levels. However, selection through repeated ART-exposure did, even-though the known resistance genes, including k13, were not further altered, suggesting other causes of increased resistance. We also observed a disproportionally high fitness cost of a resistance mutation in resistance gene ubp1. Such high fitness costs may explain why mutations in ubp1 and other genes functioning in the same pathway as k13 are rare. This highlights that k13 mutations are unique in their ability to cause resistance at a comparably low fitness cost.

Diversity and selection analyses identify transmission-blocking antigens as the optimal vaccine candidates in Plasmodium falciparum

Background

A highly effective vaccine for malaria remains an elusive target, at least in part due to the under-appreciated natural parasite variation. This study aimed to investigate genetic and structural variation, and immune selection of leading malaria vaccine candidates across the Plasmodium falciparum's life cycle.

Methods

We analyzed 325 P. falciparum whole genome sequences from Zambia, in addition to 791 genomes from five other African countries available in the MalariaGEN Pf3k Rdatabase. Ten vaccine antigens spanning three life-history stages were examined for genetic and structural variations, using population genetics measures, haplotype network analysis, and 3D structure selection analysis.

Findings

Among the ten antigens analyzed, only three in the transmission-blocking vaccine category display P. falciparum 3D7 as the dominant haplotype. The antigens AMA1, CSP, MSP119 and CelTOS, are much more diverse than the other antigens, and their epitope regions are under moderate to strong balancing selection. In contrast, Rh5, a blood stage antigen, displays low diversity yet slightly stronger immune selection in the merozoite-blocking epitope region. Except for CelTOS, the transmission-blocking antigens Pfs25, Pfs48/45, Pfs230, Pfs47, and Pfs28 exhibit minimal diversity and no immune selection in epitopes that induce strain-transcending antibodies, suggesting potential effectiveness of 3D7-based vaccines in blocking transmission.

Interpretations

These findings offer valuable insights into the selection of optimal vaccine candidates against P. falciparum. Based on our results, we recommend prioritizing conserved merozoite antigens and transmission-blocking antigens. Combining these antigens in multi-stage approaches may be particularly promising for malaria vaccine development initiatives.

Funding

Purdue Department of Biological Sciences; Puskas Memorial Fellowship; National Institute of Allergy and Infectious Diseases (U19AI089680).

Genome Analysis of Plasmodium falciparum: A Preliminary Observation - Sierra Leone, 2022-2023

What is already known about this topic?

Sierra Leone, with a gross domestic product (GDP) per capita below $300 and significant poverty, ranks among the world's least developed countries (LDCs). Despite its modest population of 8.6 million, the nation reports approximately 2.6 million malaria cases annually. Previously, there has been no reporting on the malaria genome data from this country.

What is added by this report?

In this study, we present the first reported whole-genome sequence analysis of 19 high parasite-density Plasmodium falciparum isolates from Sierra Leone, providing insights into the genomic epidemiology of this high-prevalence area. We found a high degree of relatedness among infections and substantial genetic diversity, consistent with the gradual reduction in overall case numbers. Moreover, our whole-genome analysis revealed that, beyond drug-resistance genes, gene families related to blood cell invasion, immune evasion, and others are undergoing directional selection. This suggests that the population in Sierra Leone has developed a relatively strong acquired immunity.

What are the implications for public health practice?

The genomic data not only facilitate the creation of single nucleotide polymorphism barcodes for case tracking but also enable the analysis of evolving transmission dynamics and selection pressures. Additionally, the samples from Sierra Leone exhibited higher selective pressures on resistance genes compared to those from Asia, a trend not commonly observed in other African samples. This suggests that less stringent healthcare systems and inconsistent treatment strategies can subject parasites to increased drug pressure, thereby accelerating the development of resistant strains.

Role of seasonal importation and genetic drift on selection for drug-resistant genotypes of Plasmodium falciparum in high-transmission settings

Historically Plasmodium falciparum has followed a pattern of drug resistance first appearing in low-transmission settings before spreading to high-transmission settings. Several features of low-transmission regions are hypothesized as explanations: higher chance of symptoms and treatment seeking, better treatment access, less within-host competition among clones and lower rates of recombination. Here, we test whether importation of drug-resistant parasites is more likely to lead to successful emergence and establishment in low-transmission or high-transmission periods of the same epidemiological setting, using a spatial, individual-based stochastic model of malaria and drug-resistance evolution calibrated for Burkina Faso. Upon controlling for the timing of importation of drug-resistant genotypes and examination of key model variables, we found that drug-resistant genotypes imported during the low-transmission season were (i) more susceptible to stochastic extinction due to the action of genetic drift, and (ii) more likely to lead to establishment of drug resistance when parasites are able to survive early stochastic loss due to drift. This implies that rare importation events are more likely to lead to establishment if they occur during a high-transmission season, but that constant importation (e.g. neighbouring countries with high levels of resistance) may produce a greater risk during low-transmission periods.

Putative molecular markers of Plasmodium falciparum resistance to antimalarial drugs in malaria parasites from Ghana

Introduction

Antimalarial drugs including artemisinin-based combination therapy (ACT) regimens and sulphadoxine-pyrimethamine (SP) are used in Ghana for malaria therapeutics and prophylaxis respectively. The genetic basis of Plasmodium falciparum development of drug resistance involves single nucleotide polymorphisms in genes encoding proteins for multiple cellular and metabolic processes. The prevalence of single nucleotide polymorphisms in nine P. falciparum genes linked to ACT and SP resistance in the malaria parasite population was determined.

Methods

Archived filter paper blood blot samples from patients aged 9 years and below with uncomplicated malaria reporting at 10 sentinel sites located in three ecological zones for the Malaria Therapeutic Efficacy Studies were used. The samples used were collected from 2007-2018 malaria transmission seasons and mutations in the genes were detected using PCR and Sanger sequencing.

Results

In all 1,142 samples were used for the study. For falcipain-2 gene (pffp2), Sanger sequencing was successful for 872 samples and were further analysed. The prevalence of the mutants was 45% (392/872) with pffp2 markers V51I and S59F occurring in 15.0% (128/872) and 3.0% (26/872) of the samples respectively. Prevalence of other P. falciparum gene mutations: coronin (pfcoronin) was 44.8% (37/90); cysteine desulfurase (pfnfs) was 73.9% (68/92); apicoplast ribosomal protein S10 (pfarps10) was 36.8% (35/95); ferredoxin (pffd) was 8.8% (8/91); multidrug resistance protein-1 (pfmrp1) was 95.2.0% (80/84); multidrug resistance protein-2 (pfmrp2) was 91.4% (32/35); dihydrofolate reductase (pfdhfr) was 99.0% (84/85); dihydropteroate synthase (pfdhps) was 72% (68/95).

Discussion

The observation of numerous mutations in these genes of interest in the Ghanaian isolates, some of which have been implicated in delayed parasite clearance is of great interest. The presence of these genotypes may account for the decline in the efficacies of ACT regimens being used to treat uncomplicated malaria in the country. The need for continuous monitoring of these genetic markers to give first-hand information on parasite susceptibility to antimalarial drugs to inform policy makers and stakeholders in malaria elimination in the country is further discussed.

Flexible and cost-effective genomic surveillance of P. falciparum malaria with targeted nanopore sequencing

Genomic surveillance of Plasmodium falciparum malaria can provide policy-relevant information about antimalarial drug resistance, diagnostic test failure, and the evolution of vaccine targets. Yet the large and low complexity genome of P. falciparum complicates the development of genomic methods, while resource constraints in malaria endemic regions can limit their deployment. Here, we demonstrate an approach for targeted nanopore sequencing of P. falciparum from dried blood spots (DBS) that enables cost-effective genomic surveillance of malaria in low-resource settings. We release software that facilitates flexible design of amplicon sequencing panels and use this software to design two target panels for P. falciparum. The panels generate 3-4 kbp reads for eight and sixteen targets respectively, covering key drug-resistance associated genes, diagnostic test antigens, polymorphic markers and the vaccine target csp. We validate our approach on mock and field samples, demonstrating robust sequencing coverage, accurate variant calls within coding sequences, the ability to explore P. falciparum within-sample diversity and to detect deletions underlying rapid diagnostic test failure.

In vitro delayed response to dihydroartemisinin of malaria parasites infecting sickle cell erythocytes

BACKGROUND

Decreased efficacy of artemisinin-based combination therapy (ACT) for Plasmodium falciparum malaria has been previously reported in patients with sickle cell disease (SCD). The main purpose of this study was to investigate the in vitro susceptibility of isolates to dihydro-artemisinin (DHA) to provide a hypothesis to explain this treatment failure.

METHODS

Isolates were collected from patients attending health centres in Abidjan with uncomplicated P. falciparum malaria. The haemoglobin type has been identified and in vitro drug sensitivity tests were conducted with the ring stage assay and maturation inhibition assay.

RESULTS

134 isolates were obtained. Parasitaemia and haemoglobin levels at inclusion were lower in patients with haemoglobin HbSS and HbSC than in patients with normal HbAA. After ex vivo RSA and drug inhibition assays, the lowest rate of parasitic growth was found with isolates from HbAS red cells. Conversely, a significantly higher survival rate of parasites ranging from 15 to 34% were observed in isolates from HbSS. Isolates with in vitro reduced DHA sensitivity correlate with lower RBC count and haematocrit and higher parasitaemia at inclusion compared to those with isolates with normal DHA sensitivity. However, this decrease of in vitro sensitivity to DHA was not associated with Kelch 13-Propeller gene polymorphism.

CONCLUSION

This study highlights an in vitro decreased sensitivity to DHA, for isolates collected from HbSS patients, not related to the Pfkelch13 gene mutations. These results are in line with recent studies pointing out the role of the redox context in the efficacy of the drug. Indeed, SCD red cells harbour a highly different ionic and redox context in comparison with normal red cells. This study offers new insights into the understanding of artemisinin selective pressure on the malaria parasite in the context of haemoglobinopathies in Africa.

Distribution of Plasmodium falciparum K13 gene polymorphisms across transmission settings in Ghana

Malaria is a significant global health concern, with a majority of cases in Sub-Saharan African nations. Numerous antimalarial drugs have been developed to counter the rampant prevalence of Plasmodium falciparum malaria. Artemisinin-based Combination Therapy (ACT) has served as the primary treatment of uncomplicated malaria in Ghana since 2005. However, a growing concern has emerged due to the escalating reports of ACT resistance, particularly in Southeast Asia, and its encroachment into Africa. Specifically, mutations in the Kelch propeller domain on chromosome 13 (Pfk13) have been linked to ACT resistance. Yet, our understanding of mutation prevalence in Africa remains largely uncharted. In this study, we compared Pfk13 sequences obtained from 172 P. falciparum samples across three ecological and transmission zones in Ghana. We identified 27 non-synonymous mutations among these sequences, of which two of the mutations, C580Y (found in two samples from the central region) and Y493H (found in one sample from the north), had previously been validated for their association with artemisinin resistance, a phenomenon widespread in Southeast Asia. The Pfk13 gene diversity was most pronounced in the northern savannah than the central forest and south coastal regions, where transmission rates are lower. The observed mutations were not significantly associated with geographical regions, suggesting a frequent spread of mutations across the country. The ongoing global surveillance of artemisinin resistance remains pivotal, and our findings provides insights into the potential spread of resistant parasites in West Africa. Furthermore, the identification of novel codon mutations in this study raises their potential association to ACT resistance, warranting further investigation through in vitro assays to ascertain their functional significance.

Mapping the genomic landscape of multidrug resistance in Plasmodium falciparum and its impact on parasite fitness

Drug-resistant Plasmodium falciparum parasites have swept across Southeast Asia and now threaten Africa. By implementing a P. falciparum genetic cross using humanized mice, we report the identification of key determinants of resistance to artemisinin (ART) and piperaquine (PPQ) in the dominant Asian KEL1/PLA1 lineage. We mapped k13 as the central mediator of ART resistance in vitro and identified secondary markers. Applying bulk segregant analysis, quantitative trait loci mapping using 34 recombinant haplotypes, and gene editing, our data reveal an epistatic interaction between mutant PfCRT and multicopy plasmepsins 2/3 in mediating high-grade PPQ resistance. Susceptibility and parasite fitness assays implicate PPQ as a driver of selection for KEL1/PLA1 parasites. Mutant PfCRT enhanced susceptibility to lumefantrine, the first-line partner drug in Africa, highlighting a potential benefit of opposing selective pressures with this drug and PPQ. We also identified that the ABCI3 transporter can operate in concert with PfCRT and plasmepsins 2/3 in mediating multigenic resistance to antimalarial agents.

Role of Seasonal Importation and Random Genetic Drift on Selection for Drug-Resistant Genotypes of Plasmodium falciparum in High Transmission Settings

Historically Plasmodium falciparum has followed a pattern of drug resistance first appearing in low transmission settings before spreading to high transmission settings. Several features of low-transmission regions are hypothesized as explanations: higher chance of symptoms and treatment seeking, better treatment access, less within-host competition among clones, and lower rates of recombination. Here, we test whether importation of drug-resistant parasites is more likely to lead to successful emergence and establishment in low-transmission or high-transmission periods of the same epidemiological setting, using a spatial, individual-based stochastic model of malaria and drug-resistance evolution calibrated for Burkina Faso. Upon controlling for the timing of importation of drug-resistant genotypes and examination of key model variables, we found that drug-resistant genotypes imported during the low transmission season were, (1) more susceptible to stochastic extinction due to the action of random genetic drift, and (2) more likely to lead to establishment of drug resistance when parasites are able to survive early stochastic loss due to drift. This implies that rare importation events are more likely to lead to establishment if they occur during a high-transmission season, but that constant importation (e.g., neighboring countries with high levels of resistance) may produce a greater risk during low-transmission periods.

A Plasmodium falciparum genetic cross reveals the contributions of pfcrt and plasmepsin II/III to piperaquine drug resistance

Piperaquine (PPQ) is widely used in combination with dihydroartemisinin (DHA) as a first-line treatment against malaria parasites. Multiple genetic drivers of PPQ resistance have been reported, including mutations in the Plasmodium falciparum chloroquine resistance transporter (pfcrt) and increased copies of plasmepsin II/III (pm2/3). We generated a cross between a Cambodia-derived multi-drug resistant KEL1/PLA1 lineage isolate (KH004) and a drug susceptible parasite isolated in Malawi (Mal31). Mal31 harbors a wild-type (3D7-like) pfcrt allele and a single copy of pm2/3, while KH004 has a chloroquine-resistant (Dd2-like) pfcrt allele with an additional G367C substitution and four copies of pm2/3. We recovered 104 unique recombinant progeny and examined a targeted set of progeny representing all possible combinations of variants at pfcrt and pm2/3 for detailed analysis of competitive fitness and a range of PPQ susceptibility phenotypes, including PPQ survival assay (PSA), area under the dose-response curve (AUC), and a limited point IC50 (LP-IC50). We find that inheritance of the KH004 pfcrt allele is required for PPQ resistance, whereas copy number variation in pm2/3 further enhances resistance but does not confer resistance in the absence of PPQ-R-associated mutations in pfcrt. Deeper investigation of genotype-phenotype relationships demonstrates that progeny clones from experimental crosses can be used to understand the relative contributions of pfcrt, pm2/3, and parasite genetic background, to a range of PPQ-related traits and confirm the critical role of the PfCRT G367C substitution in PPQ resistance.

How do the mutations in PfK13 protein promote anti-malarial drug resistance?

Plasmodium falciparum develops resistance to artemisinin upon exposure to the anti-malarial drug. Various mutations in the Plasmodium falciparum Kelch13 (PfK13) protein such as Y493H, R539T, I543T and C580Y have been associated with anti-malarial drug resistance. These mutations impede the regular ubiquitination process that eventually invokes drug resistance. However, the relationship between the mutation and the mechanism of drug resistance has not yet been fully elucidated. The comparative protein dynamics are studied by performing the classical molecular dynamics (MD) simulations and subsequent analysis of the trajectories adopting root-mean-square fluctuations, the secondary-structure predictions and the dynamical cross-correlation matrix analysis tools. Here, we observed that the mutations in the Kelch-domain do not have any structural impact on the mutated site; however, it significantly alters the overall dynamics of the protein. The loop-region of the BTB-domain especially for Y493H and C580Y mutants is found to have the enhanced dynamical fluctuations. The enhanced fluctuations in the BTB-domain could affect the protein-protein (PfK13-Cullin) binding interactions in the ubiquitination process and eventually lead to anti-malarial drug resistance.Communicated by Ramaswamy H. Sarma.

Molecular surveillance of Plasmodium falciparum drug-resistance markers in Vietnam using multiplex amplicon sequencing (2000-2016)

Emergence and spread of Plasmodium falciparum resistance to artemisinin-based combination therapies (ACT) is a major challenge for Greater Mekong Subregion countries in their goal to eliminate malaria by 2030. Tools to efficiently monitor drug resistance beyond resource-demanding therapeutic efficacy studies are necessary. A custom multiplex amplicon sequencing assay based on Illumina technology was designed to target the marker of partial resistance to artemisinin (K13), five candidate modulators of artemisinin resistance, the marker of resistance to chloroquine (crt), and four neutral microsatellite loci. The assay was used to genotype 635 P. falciparum-positive blood samples collected across seven provinces of Vietnam and one of Cambodia between 2000 and 2016. Markers of resistance to artemisinin partner-drugs piperaquine (copy number of plasmepsin-2) and mefloquine (copy number of multidrug-resistance 1) were determined by qPCR. Parasite population structure was further assessed using a 101-SNP barcode. Validated mutations of artemisinin partial resistance in K13 were found in 48.1% of samples, first detection was in 2000, and by 2015 prevalence overcame > 50% in Central Highlands and Binh Phuoc province. K13-C580Y variant became predominant country-wide, quickly replacing an outbreak of K13-I543T in Central Highlands. Mutations in candidate artemisinin resistance modulator genes paralleled the trends of K13 mutants, whereas resistance to piperaquine and mefloquine remained low (≈ 10%) by 2015-2016. Genomic tools applied to malaria surveillance generate comprehensive information on dynamics of drug resistance and population structure and reflect drug efficacy profiles from in vivo studies.

A human antibody epitope map of the malaria vaccine antigen Pfs25

Pfs25 is a leading antigen for a malaria transmission-blocking vaccine and shows moderate transmission-blocking activity and induction of rapidly decreasing antibody titers in clinical trials. A comprehensive definition of all transmission-reducing epitopes of Pfs25 will inform structure-guided design to enhance Pfs25-based vaccines, leading to potent transmission-blocking activity. Here, we compiled a detailed human antibody epitope map comprising epitope binning data and structures of multiple human monoclonal antibodies, including three new crystal structures of Pfs25 in complex with transmission-reducing antibodies from Malian volunteers immunized with Pfs25 conjugated to EPA and adjuvanted with AS01. These structures revealed additional epitopes in Pfs25 capable of reducing transmission and expanded this characterization to malaria-exposed humans. This work informs immunogen design to focus the antibody response to transmission-reducing epitopes of Pfs25, enabling development of more potent transmission-blocking vaccines for malaria.

Multiplexed ddPCR-amplicon sequencing reveals isolated Plasmodium falciparum populations amenable to local elimination in Zanzibar, Tanzania

Zanzibar has made significant progress toward malaria elimination, but recent stagnation requires novel approaches. We developed a highly multiplexed droplet digital PCR (ddPCR)-based amplicon sequencing method targeting 35 microhaplotypes and drug-resistance loci, and successfully sequenced 290 samples from five districts covering both main islands. Here, we elucidate fine-scale Plasmodium falciparum population structure and infer relatedness and connectivity of infections using an identity-by-descent (IBD) approach. Despite high genetic diversity, we observe pronounced fine-scale spatial and temporal parasite genetic structure. Clusters of near-clonal infections on Pemba indicate persistent local transmission with limited parasite importation, presenting an opportunity for local elimination efforts. Furthermore, we observe an admixed parasite population on Unguja and detect a substantial fraction (2.9%) of significantly related infection pairs between Zanzibar and the mainland, suggesting recent importation. Our study provides a high-resolution view of parasite genetic structure across the Zanzibar archipelago and provides actionable insights for prioritizing malaria elimination efforts.

Targeted and whole-genome sequencing reveal a north-south divide in P. falciparum drug resistance markers and genetic structure in Mozambique

Mozambique is one of the four African countries which account for over half of all malaria deaths worldwide, yet little is known about the parasite genetic structure in that country. We performed P. falciparum amplicon and whole genome sequencing on 2251 malaria-infected blood samples collected in 2015 and 2018 in seven provinces of Mozambique to genotype antimalarial resistance markers and interrogate parasite population structure using genome-wide microhaplotyes. Here we show that the only resistance-associated markers observed at frequencies above 5% were pfmdr1-184F (59%), pfdhfr-51I/59 R/108 N (99%) and pfdhps-437G/540E (89%). The frequency of pfdhfr/pfdhps quintuple mutants associated with sulfadoxine-pyrimethamine resistance increased from 80% in 2015 to 89% in 2018 (p < 0.001), with a lower expected heterozygosity and higher relatedness of microhaplotypes surrounding pfdhps mutants than wild-type parasites suggestive of recent selection. pfdhfr/pfdhps quintuple mutants also increased from 72% in the north to 95% in the south (2018; p < 0.001). This resistance gradient was accompanied by a concentration of mutations at pfdhps-436 (17%) in the north, a south-to-north increase in the genetic complexity of P. falciparum infections (p = 0.001) and a microhaplotype signature of regional differentiation. The parasite population structure identified here offers insights to guide antimalarial interventions and epidemiological surveys.

Mapping the genomic landscape of multidrug resistance in Plasmodium falciparum and its impact on parasite fitness

Drug-resistant Plasmodium falciparum parasites have swept across Southeast Asia and now threaten Africa. By implementing a P. falciparum genetic cross using humanized mice, we report the identification of key determinants of resistance to artemisinin (ART) and piperaquine (PPQ) in the dominant Asian KEL1/PLA1 lineage. We mapped k13 as the central mediator of ART resistance and identified secondary markers. Applying bulk segregant analysis, quantitative trait loci mapping and gene editing, our data reveal an epistatic interaction between mutant PfCRT and multicopy plasmepsins 2/3 in mediating high-grade PPQ resistance. Susceptibility and parasite fitness assays implicate PPQ as a driver of selection for KEL1/PLA1 parasites. Mutant PfCRT enhanced susceptibility to lumefantrine, the first-line partner drug in Africa, highlighting a potential benefit of opposing selective pressures with this drug and PPQ. We also identified that the ABCI3 transporter can operate in concert with PfCRT and plasmepsins 2/3 in mediating multigenic resistance to antimalarial agents.

Malaria outbreak in Laos driven by a selective sweep for Plasmodium falciparum kelch13 R539T mutants: a genetic epidemiology analysis

BACKGROUND

Malaria outbreaks are important public health concerns that can cause resurgence in endemic regions approaching elimination. We investigated a Plasmodium falciparum outbreak in Attapeu Province, Laos, during the 2020-21 malaria season, using genomic epidemiology methods to elucidate parasite population dynamics and identify its causes.

METHODS

In this genetic analysis, 2164 P falciparum dried blood spot samples were collected from southern Laos between Jan 1, 2017, and April 1, 2021, which included 249 collected during the Attapeu outbreak between April 1, 2020, and April 1, 2021, by routine surveillance. Genetic barcodes obtained from these samples were used to investigate epidemiological changes underpinning the outbreak, estimate population diversity, and analyse population structure. Whole-genome sequencing data from additional historical samples were used to reconstruct the ancestry of outbreak strains using identity-by-descent analyses.

FINDINGS

The outbreak parasite populations were characterised by unprecedented loss of genetic diversity, primarily caused by rapid clonal expansion of a multidrug-resistant strain (LAA1) carrying the kelch13 Arg539Thr (R539T) mutation. LAA1 replaced kelch13 Cys580Tyr (C580Y) mutants resistant to dihydroartemisinin-piperaquine (KEL1/PLA1) as the dominant strain. LAA1 inherited 58·8% of its genome from a strain circulating in Cambodia in 2008. A secondary outbreak strain (LAA2) carried the kelch13 C580Y allele, and a genome that is essentially identical to a Cambodian parasite from 2009. A third, low-frequency strain (LAA7) was a recombinant of KEL1/PLA1 with a kelch13 R539T mutant.

INTERPRETATION

These results strongly suggest that the outbreak was driven by a selective sweep, possibly associated with multidrug-resistant phenotypes of the outbreak strains. Established resistant populations can circulate at low frequencies for years before suddenly overwhelming dominant strains when the conditions for selection become favourable-eg, when front-line therapies change. Genetic surveillance can support elimination by characterising key properties of outbreaks such as population diversity, drug resistance marker prevalence, and the origins of outbreak strains.

FUNDING

Bill & Melinda Gates Foundation; The Global Fund to Fight AIDS, Tuberculosis and Malaria; Wellcome Trust.

TRANSLATION

For the Lao translation of the abstract see Supplementary Materials section.

Novel Plasmodium falciparum k13 gene polymorphisms from Kisii County, Kenya during an era of artemisinin-based combination therapy deployment

BACKGROUND

Currently, chemotherapy stands out as the major malaria intervention strategy, however, anti-malarial resistance may hamper global elimination programs. Artemisinin-based combination therapy (ACT) stands as the drug of choice for the treatment of Plasmodium falciparum malaria. Plasmodium falciparum kelch13 gene mutations are associated with artemisinin resistance. Thus, this study was aimed at evaluating the circulation of P. falciparum k13 gene polymorphisms from Kisii County, Kenya during an era of ACT deployment.

METHODS

Participants suspected to have malaria were recruited. Plasmodium falciparum was confirmed using the microscopy method. Malaria-positive patients were treated with artemether-lumefantrine (AL). Blood from participants who tested positive for parasites after day 3 was kept on filter papers. DNA was extracted using chelex-suspension method. A nested polymerase chain reaction (PCR) was conducted and the second-round products were sequenced using the Sanger method. Sequenced products were analysed using DNAsp 5.10.01 software and then blasted on the NCBI for k13 propeller gene sequence identity using the Basic Local Alignment Search Tool (BLAST). To assess the selection pressure in P. falciparum parasite population, Tajima' D statistic and Fu & Li's D test in DnaSP software 5.10.01 was used.

RESULTS

Out of 275 enrolled participants, 231 completed the follow-up schedule. 13 (5.6%) had parasites on day 28 hence characterized for recrudescence. Out of the 13 samples suspected of recrudescence, 5 (38%) samples were positively amplified as P. falciparum, with polymorphisms in the k13-propeller gene detected. Polymorphisms detected in this study includes R539T, N458T, R561H, N431S and A671V, respectively. The sequences have been deposited in NCBI with bio-project number PRJNA885380 and accession numbers SAMN31087434, SAMN31087433, SAMN31087432, SAMN31087431 and SAMN31087430 respectively.

CONCLUSIONS

WHO validated polymorphisms in the k13-propeller gene previously reported to be associated with ACT resistance were not detected in the P. falciparum isolates from Kisii County, Kenya. However, some previously reported un-validated k13 resistant single nucleotide polymorphisms were reported in this study but with limited occurrences. The study has also reported new SNPs. More studies need to be carried out in the entire country to understand the association of reported mutations if any, with ACT resistance.

Assessment of artemisinin tolerance in Plasmodium falciparum clinical isolates in children with uncomplicated malaria in Ghana

BACKGROUND

Artemisinin-based combination therapy (ACT) is the first-line treatment for uncomplicated malaria in Ghana. Artemisinin (ART) tolerance in Plasmodium falciparum has arisen in Southeast Asia and recently, in parts of East Africa. This is ascribed to the survival of ring-stage parasites post treatment. The present study sought to assess and characterize correlates of potential ART tolerance based on post-treatment parasite clearance, ex vivo and in vitro drug sensitivity, and molecular markers of drug resistance in P. falciparum isolates from children with uncomplicated malaria in Ghana.

METHODS

Six months to fourteen years old children presenting with acute uncomplicated malaria (n = 115) were enrolled in two hospitals and a Health Centre in Ghana's Greater Accra region and treated with artemether-lumefantrine (AL) according to body weight. Pre- and post-treatment parasitaemia (day 0 and day 3) was confirmed by microscopy. The ex vivo ring-stage survival assay (RSA) was used to detect percent ring survival while the 72 h SYBR Green I assay was used to measure the 50% inhibition concentration (IC_50s) of ART and its derivatives and partner drugs. Genetic markers of drug tolerance /resistance were evaluated using selective whole genome sequencing.

RESULTS

Of the total of 115 participants, 85 were successfully followed up on day 3 post-treatment and 2/85 (2.4%) had parasitaemia. The IC₅₀ values of ART, artesunate (AS), artemether (AM), dihydroartemisinin (DHA), amodiaquine (AQ), and lumefantrine (LUM) were not indicative of drug tolerance. However, 7/90 (7.8%) pre-treatment isolates had > 10% ring survival rates against DHA. Of the four isolates (2 RSA positive and 2 RSA negative) with high genomic coverage, P. falciparum (Pf) kelch 13 K188* and Pfcoronin V424I mutations were only present in the two RSA positive isolates with > 10% ring survival rates.

CONCLUSIONS

The observed low proportion of participants with day-3 post-treatment parasitaemia is consistent with rapid ART clearance. However, the increased rates of survival observed in the ex vivo RSA against DHA, maybe a pointer of an early start of ART tolerance. Furthermore, the role of two novel mutations in PfK13 and Pfcoronin genes, harboured by the two RSA positive isolates that had high ring survival in the present study, remains to be elucidated.

Characterization of a novel Plasmodium falciparum merozoite surface antigen and potential vaccine target

Introduction

Detailed analyses of genetic diversity, antigenic variability, protein localization and immunological responses are vital for the prioritization of novel malaria vaccine candidates. Comprehensive approaches to determine the most appropriate antigen variants needed to provide broad protection are challenging and consequently rarely undertaken.

Methods

Here, we characterized PF3D7_1136200, which we named Asparagine-Rich Merozoite Antigen (ARMA) based on the analysis of its sequence, localization and immunogenicity. We analyzed IgG and IgM responses against the common variants of ARMA in independent prospective cohort studies in Burkina Faso (N = 228), Kenya (N = 252) and Mali (N = 195) using a custom microarray, Div-KILCHIP.

Results

We found a marked population structure between parasites from Africa and Asia. African isolates shared 34 common haplotypes, including a dominant pair although the overall selection pressure was directional (Tajima's D = -2.57; Fu and Li's F = -9.69; P < 0.02). ARMA was localized to the merozoite surface, IgG antibodies induced Fc-mediated degranulation of natural killer cells and strongly inhibited parasite growth in vitro. We found profound serological diversity, but IgG and IgM responses were highly correlated and a hierarchical clustering analysis identified only three major serogroups. Protective IgG and IgM antibodies appeared to target both cross-reactive and distinct epitopes across variants. However, combinations of IgG and IgM antibodies against selected variants were associated with complete protection against clinical episodes of malaria.

Discussion

Our systematic strategy exploits genomic data to deduce the handful of antigen variants with the strongest potential to induce broad protection and may be broadly applicable to other complex pathogens for which effective vaccines remain elusive.

Ring-stage growth arrest: Metabolic basis of artemisinin tolerance in Plasmodium falciparum

The emergence and spread of artemisinin-tolerant malaria parasites threatens malaria control programmes worldwide. Mutations in the propeller domain of the Kelch13 protein confer Plasmodium falciparum artemisinin resistance (ART-R). ART-R is linked to the reduced susceptibility of temporary growth-arrested ring-stage parasites, but the metabolic mechanisms remain elusive. We generated two PfKelch13 mutant lines via CRISPR-Cas9 gene editing which displayed a reduced susceptibility accompanied by an extended ring stage. The metabolome of ART-induced ring-stage growth arrest parasites carrying PfKelch13 mutations showed significant alterations in the tricarboxylic acid (TCA) cycle, glycolysis, and amino acids metabolism, pointing to altered energy and porphyrin metabolism with metabolic plasticity. The critical role of these pathways was further confirmed by altering metabolic flow or through chemical inhibition. Our findings uncover that the growth arrestment associated with ART-R is potentially attributed to the adaptative metabolic plasticity, indicating that the defined metabolic remodeling turns out to be the trigger for ART-R.

A world free of malaria: It is time for Africa to actively champion and take leadership of elimination and eradication strategies

The global burden of malaria seems unabated. Africa carries the greatest burden accounting for over 95% of the annual cases of malaria. For the vision of a world free of malaria by Global Technical Strategy to be achieved, Africa must take up the stakeholder's role. It is therefore imperative that Africa rises up to the challenge of malaria and champion the fight against it. The fight against malaria may just be a futile or mere academic venture if Africans are not directly and fully involved. This work reviews the roles playable by Africans in order to curb the malaria in Africa and the world at large.

Gene Coverage Count and Classification (GC3): a locus sequence coverage assessment tool using short-read whole genome sequencing data, and its application to identify and classify histidine-rich protein 2 and 3 deletions in Plasmodium falciparum

BACKGROUND

The ability of malaria rapid diagnostic tests (RDTs) to effectively detect active infections is being compromised by the presence of malaria strains with genomic deletions at the hrp2 and hrp3 loci, encoding the antigens most commonly targeted in diagnostics for Plasmodium falciparum detection. The presence of such deletions can be determined in publically available P. falciparum whole genome sequencing (WGS) datasets. A computational approach was developed and validated, termed Gene Coverage Count and Classification (GC3), to analyse genome-wide sequence coverage data and provide informative outputs to assess presence and coverage profile of a target locus in WGS data. GC3 was applied to detect deletions at hrp2 and hrp3 (hrp2/3) and flanking genes in different geographic regions and across time points.

METHODS

GC3 uses Python and R scripts to extract locus read coverage metrics from mapped WGS data according to user-defined parameters and generates relevant tables and figures. GC3 was tested using WGS data for laboratory reference strains with known hrp2/3 genotypes, and its results compared to those of a hrp2/3-specific qPCR assay. Samples with at least 25% of coding region positions with zero coverage were classified as having a deletion. Publicly available sequence data was analysed and compared with published deletion frequency estimates.

RESULTS

GC3 results matched the expected coverage of known laboratory reference strains. Agreement between GC3 and a hrp2/3-specific qPCR assay reported for 19/19 (100%) hrp2 deletions and 18/19 (94.7%) hrp3 deletions. Among Cambodian (n = 127) and Brazilian (n = 20) WGS datasets, which had not been previously analysed for hrp2/3 deletions, GC3 identified hrp2 deletions in three and four samples, and hrp3 deletions in 10 and 15 samples, respectively. Plots of hrp2/3 coding regions, grouped by year of sample collection, showed a decrease in median standardized coverage among Malawian samples (n = 150) suggesting the importance of a careful, properly controlled follow up to determine if an increase in frequency of deletions has occurred between 2007-2008 and 2014-2015. Among Malian (n = 90) samples, median standardized coverage was lower in 2002 than 2010, indicating widespread deletions present at the gene locus in 2002.

CONCLUSIONS

The GC3 tool accurately classified hrp2/3 deletions and provided informative tables and figures to analyse targeted gene coverage. GC3 is an appropriate tool when performing preliminary and exploratory assessment of locus coverage data.

Structural analysis of Plasmodium falciparum ookinete surface antigen Pfs28 relevant for malaria vaccine design

Pfs28 is a Plasmodium falciparum malaria transmission-blocking vaccine candidate that is anchored to the parasite surface through a C-terminal glycosylphosphatidylinositol (GPI) moiety, and plays a role in parasite survival in the mosquito midgut. Pfs28 contains epidermal growth factor (EGF)-like domains and is part of a family of sexual stage malaria proteins that includes the related vaccine antigen Pfs25. The lack of structural definition of Pfs28 and the immune response to this candidate has limited further malaria vaccine development for this antigen. Here, we present the crystal structure of Pfs28, examine its conservation with P. vivax Pvs28, and evaluate the cross-reactivity of Pfs28 to antibodies that recognize Pfs25. Pfs28 is comprised of four EGF-like domains stabilized by ten disulfide bridges with an overall architecture that highly resembles Pfs25. Despite the high sequence and structural similarity between these antigens, no cross reactivity of Pfs28 to anti-Pfs25 monoclonal antibodies could be demonstrated.

Genomic approaches for monitoring transmission dynamics of malaria: A case for malaria molecular surveillance in Sub-Saharan Africa

Transmission dynamics is an important indicator for malaria control and elimination. As we move closer to eliminating malaria in Sub-Saharan Africa (sSA), transmission indices with higher resolution (genomic approaches) will complement our current measurements of transmission. Most of the present programmatic knowledge of malaria transmission patterns are derived from assessments of epidemiologic and clinical data, such as case counts, parasitological estimates of parasite prevalence, and Entomological Inoculation Rates (EIR). However, to eliminate malaria from endemic areas, we need to track changes in the parasite population and how they will impact transmission. This is made possible through the evolving field of genomics and genetics, as well as the development of tools for more in-depth studies on the diversity of parasites and the complexity of infections, among other topics. If malaria elimination is to be achieved globally, country-specific elimination activities should be supported by parasite genomic data from regularly collected blood samples for diagnosis, surveillance and possibly from other programmatic interventions. This presents a unique opportunity to track the spread of malaria parasites and shed additional light on intervention efficacy. In this review, various genetic techniques are highlighted along with their significance for an enhanced understanding of transmission patterns in distinct topological settings throughout Sub-Saharan Africa. The importance of these methods and their limitations in malaria surveillance to guide control and elimination strategies, are explored.

Multidisciplinary Investigations of Sustained Malaria Transmission in the Greater Mekong Subregion

In the course of malaria elimination in the Greater Mekong Subregion (GMS), malaria epidemiology has experienced drastic spatiotemporal changes with residual transmission concentrated along international borders and the rising predominance of Plasmodium vivax. The emergence of Plasmodium falciparum parasites resistant to artemisinin and partner drugs renders artemisinin-based combination therapies less effective while the potential spread of multidrug-resistant parasites elicits concern. Vector behavioral changes and insecticide resistance have reduced the effectiveness of core vector control measures. In recognition of these problems, the Southeast Asian International Center of Excellence for Malaria Research (ICEMR) has been conducting multidisciplinary research to determine how human migration, antimalarial drug resistance, vector behavior, and insecticide resistance sustain malaria transmission at international borders. These efforts allow us to comprehensively understand the ecology of border malaria transmission and develop population genomics tools to identify and track parasite introduction. In addition to employing in vivo, in vitro, and molecular approaches to monitor the emergence and spread of drug-resistant parasites, we also use genomic and genetic methods to reveal novel mechanisms of antimalarial drug resistance of parasites. We also use omics and population genetics approaches to study insecticide resistance in malaria vectors and identify changes in mosquito community structure, vectorial potential, and seasonal dynamics. Collectively, the scientific findings from the ICEMR research activities offer a systematic view of the factors sustaining residual malaria transmission and identify potential solutions to these problems to accelerate malaria elimination in the GMS.

Molecular Epidemiology of Drug Resistance Genes in Plasmodium falciparum Isolates Imported from Nigeria between 2016 and 2020: Continued Emergence of Fully Resistant Pfdhfr - Pfdhps Alleles

Malaria poses public health threats worldwide. Nigeria accounted for the highest numbers of cases (26.8%) and deaths (31.9%) among countries where malaria is endemic in 2020. Currently, monitoring molecular markers in imported malaria cases provides an efficient means to screen for emerging drug resistance in countries where malaria is endemic, particularly in those where field surveillance is challenging. Here, we investigated 165 Plasmodium falciparum infections imported from Nigeria to Zhejiang Province, China, between 2016 and 2020. Multiple molecular markers in k13, Pfcrt, Pfmdr1, Pfdhfr, and Pfdhps were detected. The prevalences and patterns of mutations were analyzed. Polymorphism of k13 was limited to 5 of 156 (3.21%) isolates. The wild-type CVMNK allele of Pfcrt became predominant (65.36%) compared with the triple mutation CVIET. A low frequency (4.73%) of double mutations (N86Y and Y184F) in Pfmdr1 was observed. The dominant haplotypes of Pfdhfr and Pfdhps were IRNDI (92.41%) and ISGKAA (36.84%), respectively. The newly discovered mutant I431V was identified in 21.71% of isolates. A “fully resistant” combination of Pfdhfr-Pfdhps, IRN-GE, was found in eight (5.67%) samples, which was hardly seen in Nigeria. The current study demonstrated a high frequency of wild-type Pfcrt. Limited polymorphism of Pfmdr1 but a high prevalence of Pfdhfr and Pfdhps mutations was illustrated. Our data so far serve as comprehensive surveillance of molecular markers of the k13, Pfcrt, Pfmdr1, Pfdhfr, and Pfdhps genes. Based on our findings, it has become crucial to evaluate the impact of the emerging fully resistant type of Pfdhfr-Pfdhps as well as its combination with I431V on the efficacy of sulfadoxine-pyrimethamine (SP) in Nigeria.

IMPORTANCE Monitoring the current resistance to antimalarial drugs is critical to enable timely action to prevent its spread and limit its impact. The high prevalence of wild-type Pfcrt found in our study is an optimistic signal to reevaluate chloroquine (CQ) sensitivity in Nigeria, which is cost-effective and once played a crucial role in the fight against malaria. Based on the continued emergence of fully resistant Pfdhfr-Pfdhps alleles illustrated in the current investigation, actions are needed in Nigeria, such as national systemic surveillance to monitor their updated epidemiology as well as assessments of their influence on SP efficacy to minimize any public health impact. These findings urge a response to the threat of drug resistance to facilitate appropriate drug policies in the study area.

Leveraging genome editing to functionally evaluate Plasmodium diversity

The ambitious goal of malaria elimination requires an in-depth understanding of the parasite's biology to counter the growing threat of antimalarial resistance and immune evasion. Timely assessment of the functional impact of antigenic diversity in the early stages of vaccine development will be critical for achieving the goal of malaria control, elimination, and ultimately eradication. Recent advances in targeted genome editing enabled the functional validation of resistance-associated markers in Plasmodium falciparum, the deadliest malaria-causing pathogen and strain-specific immune neutralization. This review explores recent advances made in leveraging genome editing to aid the functional evaluation of Plasmodium diversity and highlights how these techniques can assist in prioritizing both therapeutic and vaccine candidates.

Comparison of antibody responses and parasite clearance in artemisinin therapeutic efficacy studies in Democratic Republic of Congo and Asia

Background

Understanding the effect of immunity on Plasmodium falciparum clearance is essential for interpreting therapeutic efficacy studies designed to monitor emergence of artemisinin drug resistance. In low-transmission areas of Southeast Asia, where resistance has emerged, P. falciparum antibodies confound parasite clearance measures. However, variation in naturally acquired antibodies across Asian and sub-Saharan African epidemiological contexts and their impact on parasite clearance re yet to be quantified.

Methods

In an artemisinin therapeutic efficacy study, antibodies to 12 pre-erythrocytic and erythrocytic P. falciparum antigens were measured in 118 children with uncomplicated P. falciparum malaria in the Democratic Republic of Congo (DRC) and compared with responses in patients from Asian sites, described elsewhere.

Results

Parasite clearance half-life was shorter in DRC patients (median, 2 hours) compared with most Asian sites (median, 2–7 hours), but P. falciparum antibody levels and seroprevalences were similar. There was no evidence for an association between antibody seropositivity and parasite clearance half-life (mean difference between seronegative and seropositive, −0.14 to +0.40 hour) in DRC patients.

Conclusions

In DRC, where artemisinin remains highly effective, the substantially shorter parasite clearance time compared with Asia was not explained by differences in the P. falciparum antibody responses studied.

Novel pfk13 polymorphisms in Plasmodium falciparum population in Ghana

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.

Identification of Co-Existing Mutations and Gene Expression Trends Associated With K13-Mediated Artemisinin Resistance in Plasmodium falciparum

Plasmodium falciparum infects millions and kills thousands of people annually the world over. With the emergence of artemisinin and/or multidrug resistant strains of the pathogen, it has become even more challenging to control and eliminate the disease. Multiomics studies of the parasite have started to provide a glimpse into the confounding genetics and mechanisms of artemisinin resistance and identified mutations in Kelch13 (K13) as a molecular marker of resistance. Over the years, thousands of genomes and transcriptomes of artemisinin-resistant/sensitive isolates have been documented, supplementing the search for new genes/pathways to target artemisinin-resistant isolates. This meta-analysis seeks to recap the genetic landscape and the transcriptional deregulation that demarcate artemisinin resistance in the field. To explore the genetic territory of artemisinin resistance, we use genomic single-nucleotide polymorphism (SNP) datasets from 2,517 isolates from 15 countries from the MalariaGEN Network (The Pf3K project, pilot data release 4, 2015) to dissect the prevalence, geographical distribution, and co-existing patterns of genetic markers associated with/enabling artemisinin resistance. We have identified several mutations which co-exist with the established markers of artemisinin resistance. Interestingly, K13-resistant parasites harbor α-ß hydrolase and putative HECT domain-containing protein genes with the maximum number of SNPs. We have also explored the multiple, publicly available transcriptomic datasets to identify genes from key biological pathways whose consistent deregulation may be contributing to the biology of resistant parasites. Surprisingly, glycolytic and pentose phosphate pathways were consistently downregulated in artemisinin-resistant parasites. Thus, this meta-analysis highlights the genetic and transcriptomic features of resistant parasites to propel further exploratory studies in the community to tackle artemisinin resistance.

Therapeutic efficacy of artemether-lumefantrine, artesunate-amodiaquine and dihydroartemisinin-piperaquine in the treatment of uncomplicated Plasmodium falciparum malaria in Sub-Saharan Africa: A systematic review and meta-analysis

Background

Sub-Saharan Africa has the highest burden of malaria in the world. Artemisinin-based combination therapies (ACTs) have been the cornerstone in the efforts to reduce the global burden of malaria. In the effort to facilitate early detection of resistance for artemisinin derivatives and partner drugs, WHO recommends monitoring of ACT’s efficacy in the malaria endemic countries. The present systematic meta-analysis study summarises the evidence of therapeutic efficacy of the commonly used artemisinin-based combinations for the treatment of uncomplicated P. falciparum malaria in Sub-Saharan Africa after more than a decade since the introduction of the drugs.

Methods

Fifty two studies carried out from 2010 to 2020 on the efficacy of artemether-lumefantrine or dihydro-artemisinin piperaquine or artesunate amodiaquine in patients with uncomplicated P. falciparum malaria in Sub-Saharan Africa were searched for using the Google Scholar, Cochrane Central Register of controlled trials (CENTRAL), PubMed, Medline, LILACS, and EMBASE online data bases. Data was extracted by two independent reviewers. Random analysis effect was performed in STATA 13. Heterogeneity was established using I² statistics.

Results

Based on per protocol analysis, unadjusted cure rates in malaria infected patients treated with artemether-lumefantrine (ALU), artesunate-amodiaquine (ASAQ) and dihydroartemisinin-piperaquine (DHP) were 89%, 94% and 91% respectively. However, the cure rates after PCR correction were 98% for ALU, 99% for ASAQ and 99% for DHP.

Conclusion

The present meta-analysis reports the overall high malaria treatment success for artemether-lumefantrine, artesunate-amodiaquine and dihydroartemisinin-piperaquine above the WHO threshold value in Sub-Saharan Africa.

Design and implementation of multiplexed amplicon sequencing panels to serve genomic epidemiology of infectious disease: a malaria case study

Multiplexed PCR amplicon sequencing (AmpSeq) is an increasingly popular application for cost-effective monitoring of threatened species and managed wildlife populations, and shows strong potential for genomic epidemiology of infectious disease. AmpSeq data from infectious microbes can inform disease control in multiple ways, including measuring drug resistance marker prevalence, distinguishing imported from local cases, and determining the effectiveness of therapeutics. We describe the design and comparative evaluation of two new AmpSeq assays for Plasmodium falciparum malaria parasites: a four-locus panel ('4CAST') composed of highly diverse antigens, and a 129-locus panel ('AMPLseq') composed of drug resistance markers, highly diverse loci for inferring relatedness, and a locus to detect Plasmodium vivax co-infection. We explore the performance of each panel in various public health use cases with in silico simulations as well as empirical experiments. The 4CAST panel appears highly suitable for evaluating the number of distinct parasite strains within samples (complexity of infection), showing strong performance across a wide range of parasitemia levels without a DNA pre-amplification step. For relatedness inference, the larger AMPLseq panel performs similarly to two existing panels of comparable size, despite differences in the data and approach used for designing each panel. Finally, we describe an R package (paneljudge) that facilitates the design and comparative evaluation of genetic panels for relatedness estimation, and we provide general guidance on the design and implementation of AmpSeq panels for genomic epidemiology of infectious disease.

Optimizing bulk segregant analysis of drug resistance using Plasmodium falciparum genetic crosses conducted in humanized mice

Classical malaria parasite genetic crosses involve isolation, genotyping, and phenotyping of progeny parasites, which is time consuming and laborious. We tested a rapid alternative approach-bulk segregant analysis (BSA)-that utilizes sequencing of bulk progeny populations with and without drug selection for rapid identification of drug resistance loci. We used dihydroartemisinin (DHA) selection in two genetic crosses and investigated how synchronization, cryopreservation, and the drug selection regimen impacted BSA success. We detected a robust quantitative trait locus (QTL) at kelch13 in both crosses but did not detect QTLs at four other candidate loci. QTLs were detected using synchronized, but not unsynchronized progeny pools, consistent with the stage-specific action of DHA. We also successfully applied BSA to cryopreserved progeny pools, expanding the utility of this approach. We conclude that BSA provides a powerful approach for investigating the genetic architecture of drug resistance in Plasmodium falciparum.

An open dataset of Plasmodium vivax genome variation in 1,895 worldwide samples

This report describes the MalariaGEN Pv4 dataset, a new release of curated genome variation data on 1,895 samples of Plasmodium vivax collected at 88 worldwide locations between 2001 and 2017. It includes 1,370 new samples contributed by MalariaGEN and VivaxGEN partner studies in addition to previously published samples from these and other sources. We provide genotype calls at over 4.5 million variable positions including over 3 million single nucleotide polymorphisms (SNPs), as well as short indels and tandem duplications. This enlarged dataset highlights major compartments of parasite population structure, with clear differentiation between Africa, Latin America, Oceania, Western Asia and different parts of Southeast Asia. Each sample has been classified for drug resistance to sulfadoxine, pyrimethamine and mefloquine based on known markers at the dhfr, dhps and mdr1 loci. The prevalence of all of these resistance markers was much higher in Southeast Asia and Oceania than elsewhere. This open resource of analysis-ready genome variation data from the MalariaGEN and VivaxGEN networks is driven by our collective goal to advance research into the complex biology of P. vivax and to accelerate genomic surveillance for malaria control and elimination.

Evaluations of candidate markers of dihydroartemisinin-piperaquine resistance in Plasmodium falciparum isolates from the China-Myanmar, Thailand-Myanmar, and Thailand-Cambodia borders

Background

The fast-declining clinical efficacy of dihydroartemisinin-piperaquine (DHA-PPQ) in Cambodia is a warning of the underlying westward dissemination of piperaquine resistance in the Greater Mekong Subregion (GMS). Mutations in the Plasmodium falciparum Kelch 13-propeller (PfK13) and the P. falciparum chloroquine resistance transporter (PfCRT), as well as plasmepsin 2/3 gene amplification, have been discovered as molecular markers for predicting DHA-PPQ treatment failure. Determining whether these genetic variations of P. falciparum are linked to DHA-PPQ resistance is critical, especially along the China–Myanmar (CM) border, where PPQ has been utilized for decades.

Methods

A total of 173 P. falciparum samples of dried blood spots (DBS) were collected along the CM border between 2007 and 2010, the Thailand–Cambodia (TC) border between 2009 and 2013, and the Thailand–Myanmar (TM) border between 2012 and 2014. PCR and sequencing were used to identified PfCRT mutations, while qPCR was used to determine the copy number of plasmepsin 2/3. The prevalence of DHA-PPQ resistance in three locations was investigated using data paired with K13 mutations.

Results

Three fragments of the pfcrt gene were amplified for all 173 samples, and seven SNPs were identified (M74I, N75E/D, K76T, H97L, I218F, A220S, I356L). No new PfCRT mutations conferring resistance to PPQ (T93S, H97Y, F145I, M343L, and G353V) were discovered, except for one mutant I218F identified in the TM border (2.27%, 1/44). Additionally, mutant H97L was found in the TC, TM, and CM borders at 3.57% (1/28), 6.82% (3/44), and 1% (1/101), respectively. A substantial K13 C580Y variant prevalence was found in the TC and TM border, accounting for 64.29% (18/28) and 43.18% (19/44), respectively, while only 1% (1/101) was found in the CM border. The K13 F446I variant was only identified and found to reach a high level (28.71%, 29/101) in the CM border. Furthermore, 10.71% (3/28) of TC isolates and 2.27% (1/44) of TM isolates carried more than one copy of plasmepsin 2/3 and K13 C580Y variant, while no plasmepsin 2/3 amplification was identified in the CM isolates.

Conclusions

Compared with the P. falciparum samples collected from the TC and TM borders, fewer parasites carried plasmepsin 2/3 amplification and novel PfCRT variants, while more parasites carried predominant K13 mutations at position F446I, in the CM border. Clear evidence of DHA-PPQ resistance associated with candidate markers was not found in this border region suggesting a further evaluation of these markers and continuous surveillance is warranted.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-022-05239-1.

Tenfold difference in DNA recovery rate: systematic comparison of whole blood vs. dried blood spot sample collection for malaria molecular surveillance

BACKGROUND

Molecular and genomic surveillance is becoming increasingly used to track malaria control and elimination efforts. Blood samples can be collected as whole blood and stored at - 20 °C until DNA extraction, or as dried blood spots (DBS), circumventing the need for a cold chain. Despite the wide use of either method, systematic comparisons of how the method of blood sample preservation affects the limit of detection (LOD) of molecular diagnosis and the proportion of DNA recovered for downstream applications are lacking.

METHODS

Extractions based on spin columns, magnetic beads, Tween-Chelex, and direct PCR without prior extraction were compared for whole blood and dried blood spots (DBS) using dilution series of Plasmodium falciparum culture samples. Extracted DNA was quantified by qPCR and droplet digital PCR (ddPCR).

RESULTS

DNA recovery was 5- to 10-fold higher for whole blood compared to DBS, resulting in a 2- to 3-fold lower LOD for both extraction methods compared to DBS. For whole blood, a magnetic bead-based method resulted in a DNA recovery rate of 88-98% when extracting from whole blood compared to 17-33% for a spin-column based method. For extractions from DBS, the magnetic bead-based method resulted in 8-20% DNA recovery, while the spin-column based method resulted in only 2% DNA recovery. The Tween-Chelex method was superior to other methods with 15-21% DNA recovery, and even more sensitive than extractions from whole blood samples. The direct PCR method was found to have the lowest LOD overall for both, whole blood and DBS.

CONCLUSIONS

Pronounced differences in LOD and DNA yield need to be considered when comparing prevalence estimates based on molecular methods and when selecting sampling protocols for other molecular surveillance applications.

Reduction in Plasmodium falciparum Pfk13 and pfg377 allele diversity through time in southern Vietnam

Background

Plasmodium falciparum has acquired resistance to artemisinin in Southeast Asia, with mutations in the P. falciparum Kelch-13 (Pfk13) gene associated with the resistance phenotype. The widespread use of Artemisinin-based combination therapy (ACT)s in Southeast Asia has led to the selection and spread of parasites carrying mutations in Pfk13. We characterised the allele diversity of Pfk13 and pfg377, an artemisinin-resistance neutral polymorphic gene, in parasite DNA extracted human blood from in southern Vietnam in 2003, 2012, 2015 and 2018.

Method

This study was conducted in Bu Gia Map commune, Binh Phuoc province, Vietnam, from May 2018 to January 2019. Twenty-four samples from 2018 to 2019, 30 from 2003, 24 from 2012 and 32 from 2015 were analysed. Malaria-infected human blood was collected by finger-prick and used for molecular analysis. A nested-PCR targeting the small subunit ribosomal RNA gene was used for Plasmodium species identification, followed by amplification and nucleotide sequencing of Pfk13 and region 3 of pfg377. Archived blood samples collected in the same region in 2012 and 2015 were also analysed as above for comparison.

Results

The genetic diversity of Pfk13 and pfg377 was lower in 2018–2019 compared to 2012 and 2015. The number of distinct Pfk13 mutants decreased from three in 2012 and 2015, P553L, V568G and C580Y, to one, C580Y in 2018–2019. In 2018–2019, the frequency of C580Y mutant strains was 71% (17/24 isolates). All samples were wild type in 2003. In 2012 and 2015, there were single-strain infections as well as co-infections with two mutant strains or with mutant and wild strains, whereas there were no co-infections in 2018. pfg377 allele diversity decreased from five alleles in 2012 to two alleles in 2018–2019.

Conclusion

The genetic diversity of P. falciparum was reduced at the two genetic loci surveyed in this study, Pfk13 and pfg377. In the case of the former gene, we observed an increase in the prevalence of parasites carrying the C580Y gene, known to confer reduced susceptibility to ACTs. The reduction in the diversity of pfg377 may be linked to the clonal expansion of parasite strains carrying the C580Y mutation, leading to an overall reduction in parasite genetic diversity across the population.

Supplementary Information

The online version contains supplementary material available at 10.1186/s41182-022-00409-4.

5WBF: a low-cost and straightforward whole blood filtration method suitable for whole-genome sequencing of Plasmodium falciparum clinical isolates

Background

Whole-genome sequencing (WGS) is becoming increasingly helpful to assist malaria control programmes. A major drawback of this approach is the large amount of human DNA compared to parasite DNA extracted from unprocessed whole blood. As red blood cells (RBCs) have a diameter of about 7–8 µm and exhibit some deformability, it was hypothesized that cheap and commercially available 5 µm filters might retain leukocytes but much less of Plasmodium falciparum-infected RBCs. This study aimed to test the hypothesis that such a filtration method, named 5WBF (for 5 µm Whole Blood Filtration), may provide highly enriched parasite material suitable for P. falciparum WGS.

Methods

Whole blood was collected from five patients experiencing a P. falciparum malaria episode (ring-stage parasitaemia range: 0.04–5.5%) and from mock samples obtained by mixing synchronized, ring-stage cultured P. falciparum 3D7 parasites with uninfected human whole blood (final parasitaemia range: 0.02–1.1%). These whole blood samples (50 to 400 µL) were diluted in RPMI 1640 medium or PBS 1× buffer and filtered with a syringe connected to a 5 µm commercial filter. DNA was extracted from 5WBF-treated and unfiltered counterpart blood samples using a commercial kit. The 5WBF method was evaluated on the ratios of parasite:human DNA assessed by qPCR and by sequencing depth and percentages of coverage from WGS data (Illumina NextSeq 500). As a comparison, the popular selective whole-genome amplification (sWGA) method, which does not rely on blood filtration, was applied to the unfiltered counterpart blood samples.

Results

After applying 5WBF, qPCR indicated an average of twofold loss in the amount of parasite template DNA (Pf ARN18S gene) and from 4096- to 65,536-fold loss of human template DNA (human β actin gene). WGS analyses revealed that > 95% of the  parasite nuclear and organellar genomes were all covered at ≥ 10× depth for all samples tested. In sWGA counterparts, the organellar genomes were poorly covered and from 47.7 to 82.1% of the nuclear genome was covered at ≥ 10× depth depending on parasitaemia. Sequence reads were homogeneously distributed across gene sequences for 5WBF-treated samples (n = 5460 genes; mean coverage: 91×; median coverage: 93×; 5th percentile: 70×; 95th percentile: 103×), allowing the identification of gene copy number variations such as for gch1. This later analysis was not possible for sWGA-treated samples, as a much more heterogeneous distribution of reads across gene sequences was observed (mean coverage: 80×; median coverage: 51×; 5th percentile: 7×; 95th percentile: 245×).

Conclusions

The novel 5WBF leucodepletion method is simple to implement and based on commercially available, standardized 5 µm filters which cost from 1.0 to 1.7€ per unit depending on suppliers. 5WBF permits extensive genome-wide analysis of P. falciparum ring-stage isolates from minute amounts of whole blood even with parasitaemias as low as 0.02%.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12936-022-04073-1.

Temporal evolution of the resistance genotypes of Plasmodium falciparum in isolates from Equatorial Guinea during 20 years (1999 to 2019)

Background

Malaria is one of the deadliest diseases in the world, particularly in Africa. As such, resistance to anti-malarial drugs is one of the most important problems in terms of global malaria control. This study assesses the evolution of the different resistance markers over time and the possible influence of interventions and treatment changes that have been made in Equatorial Guinea.

Methods

A total of 1223 biological samples obtained in the period 1999 to 2019 were included in the study. Screening for mutations in the pfdhfr, pfdhps, pfmdr1, and pfcrt genes was carried out by nested PCR and restriction-fragment length polymorphisms (RFLPs), and the study of pfk13 genes was carried out by nested PCR, followed by sequencing to determine the presence of mutations.

Results

The partially and fully resistant haplotypes (pfdhfr + pfdhps) were found to increase over time. Moreover, in 2019, the fully resistant haplotype was found to be increasing, although its super-resistant counterpart remains much less prevalent. A continued decline in pfmdr1 and pfcrt gene mutations over time was also found. The number of mutations detected in pfk13 has increased since 2008, when artemisinin-based combination therapy (ACT) were first introduced, with more mutations being observed in 2019, with two synonymous and five non-synonymous mutations being detected, although these are not related to resistance to ACT. In addition, the non-synonymous A578S mutation, which is the most frequent on the African continent, was detected in 2013, although not in the following years.

Conclusions

Withdrawal of the use of chloroquine (CQ) as a treatment in Equatorial Guinea has been shown to be effective over time, as wild-type parasite populations outnumber mutant populations. The upward trend observed in sulfadoxine-pyrimethamine (SP) resistance markers suggest its misuse, either alone or in combination with artesunate (AS) or amodiaquine (AQ), in some areas of the country, as was found in a previous study conducted by this group, which allows selective pressure from SP to continue. Single nucleotide polymorphisms (SNPs) 540E and 581G do not exceed the limit of 50 and 10%, respectively, thus meaning that SP is still effective as an intermittent preventive treatment (IPT) in this country. As for the pfk13 gene, no mutations have been detected in relation to resistance to ACT. However, in 2019 there is a greater accumulation of non-synonymous mutations compared to years prior to 2008.

Graphical Abstract

Malaria research in Australia: looking through the lens of the past towards the future

Malaria remains a global health priority, with substantial resources devoted to control and intervention since the causative parasite was first identified in 1880. Major advances have been made in discovery and translational research activities aimed at prevention, treatment and control. Laboratory-based, clinical, and field-based studies have complemented public health approaches. Australian scientists have played important roles, developing and applying innovative approaches, novel research tools and cutting-edge technologies in animal and human models of disease, as well as in disease-endemic settings. This article will provide an insight into 50 years of Australian efforts to discover mechanisms and targets of immunity and pathogenesis; develop new diagnostics, drugs, vaccines, and therapeutics; and assess new public health interventions and control measures in malaria-endemic settings.

Plasmodium falciparum rosetting protects schizonts against artemisinin

Background

Artemisinin (ART) resistance in Plasmodium falciparum is thought to occur during the early stage of the parasite's erythrocytic cycle. Here, we identify a novel factor associated with the late stage parasite development that contributes to ART resistance.

Methods

Rosetting rates of clinical isolates pre- and post- brief (one hour) exposure to artesunate (AS, an ART derivative) were evaluated. The effects of AS-mediated rosetting on the post-AS-exposed parasite's replication and survival, as well as the extent of protection by AS-mediated rosetting on different parasite stages were investigated. The rosetting ligands, mechanisms, and gene mutations involved were studied.

Findings

Brief AS exposure stimulated rosetting, with AS-resistant isolates forming more rosettes in a more rapid manner. AS-mediated rosetting enabled infected erythrocytes (IRBC) to withstand AS exposure for several hours and protected the IRBC from phagocytosis. When their rosetting ability was blocked experimentally, the post-AS exposure survival advantage by the AS-resistant parasites was abrogated. Deletions in two genes coding for PfEMP1 exon 2 (PF3D7_0200300 and PF3D7_0223300) were found to be associated with AS-mediated rosetting, and these mutations were significantly selected through time in the parasite population under study, along with the K13 mutations, a molecular marker of ART-resistance.

Interpretation

Rapid ART parasite clearance is driven by the direct oxidative damages on IRBC by ART and the phagocytic destruction of the damaged IRBC. Rosetting serves as a rapid ‘buying time’ strategy that allows more parasites to complete schizont maturation, reinvasion and subsequent development into the intrinsically less ART-susceptible ring stage.

Funding

A*STAR, NMRC-OF-YIRG, HRC e-ASIA, Wellcome.

First Detection in West Africa of a Mutation That May Contribute to Artemisinin Resistance Plasmodium falciparum

Background: The spread of drug resistance has seriously impacted the effective treatment of infection with the malaria parasite, Plasmodium falciparum. Continuous monitoring of molecular marker polymorphisms associated with drug resistance in parasites is essential for malaria control and elimination efforts. Our study describes mutations observed in the resistance genes Pfkelch13, Pfcrt, and Pfmdr1 in imported malaria and identifies additional potential drug resistance-associated molecular markers.

Methods: Chinese patients infected in Africa with P. falciparum were treated with intravenous (IV) injections of artesunate 240–360 mg for 3–5 days while hospitalized and treated with oral dihydroartemisinin-piperaquine (DHP) for 3 days after hospital discharge. Blood samples were collected and PCR sequencing performed on genes Pfkelch13, Pfcrt, and Pfmdr1 from all isolates.

Results: We analyzed a total of 225 patients from Guangxi, China with P. falciparum malaria acquired in Africa between 2016 and 2018. All patients were cured completely after treatment. The F446I mutation of the Pfkelch13 gene was detected for the first time from samples of West African P. falciparum, with a frequency of 1.0%. Five haplotypes of Pfcrt that encode residues 72–76 were found, with the wild-type CVMNK sequence predominating (80.8% of samples), suggesting that the parasites might be chloroquine sensitive. For Pfmdr1, N86Y (13.1%) and Y184F (58.8%) were the most prevalent, suggesting that artemether-lumefantrine may not, in general, be a suitable treatment for the group.

Conclusions: For the first time, this study detected the F446I mutation of the Pfkelch13 gene from Africa parasites that lacked clinical evidence of resistance. This study provides the latest data for molecular marker surveillance related to antimalarial drug resistance genes Pfkelch13, Pfcrt, and Pfmdr1 imported from Africa, in Guangxi, China from Chinese migrate workers.

Clinical Trial Registration: ChiCTROPC17013106.

Efficacy of dihydroartemisinin/piperaquine in patients with non-complicated Plasmodium falciparum malaria in Yaoundé, Cameroon

BACKGROUND

Dihydroartemisinin/piperaquine is increasingly used for the treatment of uncomplicated Plasmodium falciparum malaria in Africa. The efficacy of this combination in Cameroon is poorly documented, while resistance to dihydroartemisinin/piperaquine readily spreads in Southeast Asia.

OBJECTIVES

This study evaluated the clinical efficacy of dihydroartemisinin/piperaquine in Cameroon, as well as the molecular profile and phenotypic susceptibility of collected isolates to dihydroartemisinin and piperaquine.

PATIENTS AND METHODS

Dihydroartemisinin/piperaquine efficacy in 42 days was followed-up for 138 patients presenting non-complicated falciparum malaria. Piperaquine concentration was determined at day 7 for 124 patients. kelch13 gene polymorphisms (n = 150) and plasmepsin2 gene amplification (n = 148) were determined as molecular markers of resistance to dihydroartemisinin and piperaquine, respectively. Parasite susceptibility to dihydroartemisinin and piperaquine was determined using validated in vitro survival assays.

RESULTS

The efficacy of dihydroartemisinin/piperaquine treatment was 100% after PCR correction. The reinfections were not associated with a variation of piperaquine concentration at day 7. Ninety-six percent (144/150) of the samples presented a WT allele of the kelch13 gene. Two percent (3/150) presented the non-synonymous mutation A578S, which is not associated with resistance to dihydroartemisinin. No duplication of the plasmepsin2 gene was observed (0/148). All the samples tested in vitro by survival assays (n = 87) were susceptible to dihydroartemisinin and piperaquine.

CONCLUSIONS

Dihydroartemisinin/piperaquine has demonstrated excellent therapeutic efficacy with no evidence of emerging artemisinin or piperaquine resistance in Yaoundé, Cameroon. This observation suggests that dihydroartemisinin/piperaquine could be a sustainable therapeutic solution for P. falciparum malaria if implemented in areas previously free of artemisinin- and piperaquine-resistant parasites, unlike Southeast Asia.

Mutation in Plasmodium falciparum BTB/POZ domain of K13 protein confers artemisinin resistance

Partial artemisinin resistance, defined in patients as a delayed parasite clearance following artemisinin-based treatment, is conferred by non-synonymous mutations in the Kelch beta-propeller domain of the Plasmodium falciparum k13 (pfk13) gene. Here, we carried out in vitro selection over a one-year period on a West African P. falciparum strain isolated from Kolle (Mali) under a dose-escalating artemisinin regimen. After 18 cycles of sequential drug pressure, the selected parasites exhibited enhanced survival to dihydroartemisinin in the ring-stage survival assay (RSA^0-3h = 9.2%). Sanger and whole-genome sequence analyses identified the PfK13 P413A mutation, localized in the BTB/POZ domain, upstream of the propeller domain. This mutation was sufficient to confer in vitro artemisinin resistance when introduced into the PfK13 coding sequence of the parasite strain Dd2 by CRISPR/Cas9 gene editing. These results together with structural studies of the protein demonstrate that the propeller domain is not the sole in vitro mediator of PfK13-mediated artemisinin resistance, and highlight the importance of monitoring for mutations throughout PfK13.

Imported Malaria in Portugal: Prevalence of Polymorphisms in the Anti-Malarial Drug Resistance Genes pfmdr1 and pfk13

Malaria is one of the ‘big three’ killer infectious diseases, alongside tuberculosis and HIV. In non-endemic areas, malaria may occur in travelers who have recently been to or visited endemic regions. The number of imported malaria cases in Portugal has increased in recent years, mostly due to the close relationship with the community of Portuguese language countries. Samples were collected from malaria-infected patients attending Centro Hospitalar Lisboa Ocidental (CHLO) or the outpatient clinic of Instituto de Higiene e Medicina Tropical (IHMT-NOVA) between March 2014 and May 2021. Molecular characterization of Plasmodium falciparum pfk13 and pfmdr1 genes was performed. We analyzed 232 imported malaria cases. The majority (68.53%) of the patients came from Angola and only three patients travelled to a non-African country; one to Brazil and two to Indonesia. P. falciparum was diagnosed in 81.47% of the cases, P. malariae in 7.33%, P. ovale 6.47% and 1.72% carried P. vivax. No mutations were detected in pfk13. Regarding pfmdr1, the wild-type haplotype (N86/Y184/D1246) was also the most prevalent (64.71%) and N86/184F/D1246 was detected in 26.47% of the cases. The typical imported malaria case was middle-aged male, traveling from Angola, infected with P. falciparum carrying wild type pfmdr1 and pfk13. Our study highlights the need for constant surveillance of malaria parasites imported into Portugal as an important pillar of public health.