Genomic epidemiology of artemisinin resistant malaria

Evidence of Artemisinin-Resistant Malaria in Africa

BACKGROUND

In the six Southeast Asian countries that make up the Greater Mekong Subregion, Plasmodium falciparum has developed resistance to derivatives of artemisinin, the main component of first-line treatments for malaria. Clinical resistance to artemisinin monotherapy in other global regions, including Africa, would be problematic.

METHODS

In this longitudinal study conducted in Northern Uganda, we treated patients who had P. falciparum infection with intravenous artesunate (a water-soluble artemisinin derivative) and estimated the parasite clearance half-life. We evaluated ex vivo susceptibility of the parasite using a ring-stage survival assay and genotyped resistance-related genes.

RESULTS

From 2017 through 2019, a total of 14 of 240 patients who received intravenous artesunate had evidence of in vivo artemisinin resistance (parasite clearance half-life, >5 hours). Of these 14 patients, 13 were infected with P. falciparum parasites with mutations in the A675V or C469Y allele in the kelch13 gene. Such mutations were associated with prolonged parasite clearance half-lives (geometric mean, 3.95 hours for A675V and 3.30 hours for C469Y, vs. 1.78 hours for wild-type allele; P<0.001 and P = 0.05, respectively). The ring-stage survival assay showed a higher frequency of parasite survival among organisms with the A675V allele than among those with the wild-type allele. The prevalence of parasites with kelch13 mutations increased significantly, from 3.9% in 2015 to 19.8% in 2019, due primarily to the increased frequency of the A675V and C469Y alleles (P<0.001 and P = 0.004, respectively). Single-nucleotide polymorphisms flanking the A675V mutation in Uganda were substantially different from those in Southeast Asia.

CONCLUSIONS

The independent emergence and local spread of clinically artemisinin-resistant P. falciparum has been identified in Africa. The two kelch13 mutations may be markers for detection of these resistant parasites. (Funded by the Japan Society for the Promotion of Science and others.).

New insights into the spread of resistance to artemisinin and its analogues

Plasmodium falciparum, the causative agent of malaria, has been developing resistance to several drugs worldwide for more than five decades. Initially, resistance was against drugs such as chloroquine, pyrimethamine, sulfadoxine, mefloquine and quinine. Research studies are now reporting parasites with resistance to the most effective and novel drug used against malaria infection worldwide, namely artemisinin. For this reason, the first-line treatment strategy of artemisinin-based combination therapy is becoming unsuccessful in areas where drug resistance is highly prevalent. The increase in artemisinin-resistant P. falciparum strains has threatened international efforts to eliminate malarial infections and to reduce the disease burden. Detection of several phenotypes that display artemisinin resistance, specification of basic genetic factors, the discovery of molecular pathways, and evaluation of its clinical outcome are possible by the current series of research on genomics and transcriptomic levels in Asia and Africa. In artemisinin resistance, slow parasite clearance among malaria-infected patients and enhanced in vitro survival of parasites occurs at the early ring stage. This resistance is due to single nucleotide polymorphisms within the Kelch 13 gene of the parasite and is related to significantly upregulated resistance signalling pathways; thus, the pro-oxidant action of artemisinins can be antagonised. New strategies are required to halt the spread of artemisinin-resistant malarial parasites.

Limited polymorphism in k13 gene of Plasmodium falciparum and k12 of Plasmodium vivax isolates imported from African and Asian countries between 2014 and 2019 in Hangzhou city, China

Background

Malaria causes major public health problems globally and drug resistance hinders its control and elimination. Molecular markers associated with drug resistance are considered as a beneficial tool to monitor the disease trends, evolution and distribution so as to help improve drug policy.

Methods

We collected 148 Plasmodium falciparum and 20 Plasmodium vivax isolates imported into Hangzhou city, China between 2014 and 2019. k13 gene of P. falciparum and k12 of P. vivax were sequenced. Polymorphisms and prevalence of k13 and k12 were analyzed.

Results

Most (98.65%, 146/148) P. falciparum infections were imported from Africa, and half P. vivax cases came from Africa and the other half from Asia. Nucleotide mutation prevalence was 2.03% (3/148) and the proportion of amino acid mutations was 0.68% (1/148). The amino acid mutation, A676S, was observed in an isolate from Nigeria. No mutation of k12 was observed from the parasites from African and Asian countries.

Conclusions

Limited polymorphism in k13 gene of P. falciparum isolates imported from African countries, but no evidence for the polymorphism of k12 in P. vivax samples from African and Asian countries was found. These results provide information for drug policy update in study region.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12879-021-06579-6.

Genetic surveillance in the Greater Mekong subregion and South Asia to support malaria control and elimination

Background

National Malaria Control Programmes (NMCPs) currently make limited use of parasite genetic data. We have developed GenRe-Mekong, a platform for genetic surveillance of malaria in the Greater Mekong Subregion (GMS) that enables NMCPs to implement large-scale surveillance projects by integrating simple sample collection procedures in routine public health procedures.

Methods

Samples from symptomatic patients are processed by SpotMalaria, a high-throughput system that produces a comprehensive set of genotypes comprising several drug resistance markers, species markers and a genomic barcode. GenRe-Mekong delivers Genetic Report Cards, a compendium of genotypes and phenotype predictions used to map prevalence of resistance to multiple drugs.

Results

GenRe-Mekong has worked with NMCPs and research projects in eight countries, processing 9623 samples from clinical cases. Monitoring resistance markers has been valuable for tracking the rapid spread of parasites resistant to the dihydroartemisinin-piperaquine combination therapy. In Vietnam and Laos, GenRe-Mekong data have provided novel knowledge about the spread of these resistant strains into previously unaffected provinces, informing decision-making by NMCPs.

Conclusions

GenRe-Mekong provides detailed knowledge about drug resistance at a local level, and facilitates data sharing at a regional level, enabling cross-border resistance monitoring and providing the public health community with valuable insights. The project provides a rich open data resource to benefit the entire malaria community.

Funding

The GenRe-Mekong project is funded by the Bill and Melinda Gates Foundation (OPP11188166, OPP1204268). Genotyping and sequencing were funded by the Wellcome Trust (098051, 206194, 203141, 090770, 204911, 106698/B/14/Z) and Medical Research Council (G0600718). A proportion of samples were collected with the support of the UK Department for International Development (201900, M006212), and Intramural Research Program of the National Institute of Allergy and Infectious Diseases.

The newly discovered role of endocytosis in artemisinin resistance

Artemisinin and its derivatives (ART) are the cornerstone of malaria treatment as part of artemisinin combination therapy (ACT). However, reduced susceptibility to artemisinin as well as its partner drugs threatens the usefulness of ACTs. Single point mutations in the parasite protein Kelch13 (K13) are necessary and sufficient for the reduced sensitivity of malaria parasites to ART but several alternative mechanisms for this resistance have been proposed. Recent work found that K13 is involved in the endocytosis of host cell cytosol and indicated that this is the process responsible for resistance in parasites with mutated K13. These studies also identified a series of further proteins that act together with K13 in the same pathway, including previously suspected resistance proteins such as UBP1 and AP-2μ. Here, we give a brief overview of artemisinin resistance, present the recent evidence of the role of endocytosis in ART resistance and discuss previous hypotheses in light of this new evidence. We also give an outlook on how the new insights might affect future research.

Genomic Approaches to Drug Resistance in Malaria

Although the last two decades have seen a substantial decline in malaria incidence and mortality due to the use of insecticide-treated bed nets and artemisinin combination therapy, the threat of drug resistance is a constant obstacle to sustainable malaria control. Given that patients can die quickly from this disease, public health officials and doctors need to understand whether drug resistance exists in the parasite population, as well as how prevalent it is so they can make informed decisions about treatment. As testing for drug efficacy before providing treatment to malaria patients is impractical, researchers need molecular markers of resistance that can be more readily tracked in parasite populations. To this end, much work has been done to unravel the genetic underpinnings of drug resistance in Plasmodium falciparum. The aim of this review is to provide a broad overview of common genomic approaches that have been used to discover the alleles that drive drug response phenotypes in the most lethal human malaria parasite.

Molecular Mechanisms of Drug Resistance in Plasmodium falciparum Malaria

Understanding and controlling the spread of antimalarial resistance, particularly to artemisinin and its partner drugs, is a top priority. Plasmodium falciparum parasites resistant to chloroquine, amodiaquine, or piperaquine harbor mutations in the P. falciparum chloroquine resistance transporter (PfCRT), a transporter resident on the digestive vacuole membrane that in its variant forms can transport these weak-base 4-aminoquinoline drugs out of this acidic organelle, thus preventing these drugs from binding heme and inhibiting its detoxification. The structure of PfCRT, solved by cryogenic electron microscopy, shows mutations surrounding an electronegative central drug-binding cavity where they presumably interact with drugs and natural substrates to control transport. P. falciparum susceptibility to heme-binding antimalarials is also modulated by overexpression or mutations in the digestive vacuole membrane-bound ABC transporter PfMDR1 (P. falciparum multidrug resistance 1 transporter). Artemisinin resistance is primarily mediated by mutations in P. falciparum Kelch13 protein (K13), a protein involved in multiple intracellular processes including endocytosis of hemoglobin, which is required for parasite growth and artemisinin activation. Combating drug-resistant malaria urgently requires the development of new antimalarial drugs with novel modes of action.

Plasmodium falciparum K13 mutations in Africa and Asia impact artemisinin resistance and parasite fitness

The emergence of mutant K13-mediated artemisinin (ART) resistance in Plasmodium falciparum malaria parasites has led to widespread treatment failures across Southeast Asia. In Africa, K13-propeller genotyping confirms the emergence of the R561H mutation in Rwanda and highlights the continuing dominance of wild-type K13 elsewhere. Using gene editing, we show that R561H, along with C580Y and M579I, confer elevated in vitro ART resistance in some African strains, contrasting with minimal changes in ART susceptibility in others. C580Y and M579I cause substantial fitness costs, which may slow their dissemination in high-transmission settings, in contrast with R561H that in African 3D7 parasites is fitness neutral. In Cambodia, K13 genotyping highlights the increasing spatio-temporal dominance of C580Y. Editing multiple K13 mutations into a panel of Southeast Asian strains reveals that only the R561H variant yields ART resistance comparable to C580Y. In Asian Dd2 parasites C580Y shows no fitness cost, in contrast with most other K13 mutations tested, including R561H. Editing of point mutations in ferredoxin or mdr2, earlier associated with resistance, has no impact on ART susceptibility or parasite fitness. These data underline the complex interplay between K13 mutations, parasite survival, growth and genetic background in contributing to the spread of ART resistance.

Toward New Transmission-Blocking Combination Therapies: Pharmacokinetics of 10-Amino-Artemisinins and 11-Aza-Artemisinin and Comparison with Dihydroartemisinin and Artemether

As artemisinin combination therapies (ACTs) are compromised by resistance, we are evaluating triple combination therapies (TACTs) comprising an amino-artemisinin, a redox drug, and a third drug with a different mode of action. Thus, here we briefly review efficacy data on artemisone, artemiside, other amino-artemisinins, and 11-aza-artemisinin and conduct absorption, distribution, and metabolism and excretion (ADME) profiling in vitro and pharmacokinetic (PK) profiling in vivo via intravenous (i.v.) and oral (p.o.) administration to mice. The sulfamide derivative has a notably long murine microsomal half-life (t1/2 > 150 min), low intrinsic liver clearance and total plasma clearance rates (CLint 189.4, CLtot 32.2 ml/min/kg), and high relative bioavailability (F = 59%). Kinetics are somewhat similar for 11-aza-artemisinin (t1/2 > 150 min, CLint = 576.9, CLtot = 75.0 ml/min/kg), although bioavailability is lower (F = 14%). In contrast, artemether is rapidly metabolized to dihydroartemisinin (DHA) (t1/2 = 17.4 min) and eliminated (CLint = 855.0, CLtot = 119.7 ml/min/kg) and has low oral bioavailability (F) of 2%. While artemisone displays low t1/2 of <10 min and high CLint of 302.1, it displays a low CLtot of 42.3 ml/min/kg and moderate bioavailability (F) of 32%. Its active metabolite M1 displays a much-improved t1/2 of >150 min and a reduced CLint of 37.4 ml/min/kg. Artemiside has t1/2 of 12.4 min, CLint of 673.9, and CLtot of 129.7 ml/kg/min, likely a reflection of its surprisingly rapid metabolism to artemisone, reported here for the first time. DHA is not formed from any amino-artemisinin. Overall, the efficacy and PK data strongly support the development of selected amino-artemisinins as components of new TACTs.

An open dataset of Plasmodium falciparum genome variation in 7,000 worldwide samples

MalariaGEN is a data-sharing network that enables groups around the world to work together on the genomic epidemiology of malaria. Here we describe a new release of curated genome variation data on 7,000 Plasmodium falciparum samples from MalariaGEN partner studies in 28 malaria-endemic countries. High-quality genotype calls on 3 million single nucleotide polymorphisms (SNPs) and short indels were produced using a standardised analysis pipeline. Copy number variants associated with drug resistance and structural variants that cause failure of rapid diagnostic tests were also analysed.  Almost all samples showed genetic evidence of resistance to at least one antimalarial drug, and some samples from Southeast Asia carried markers of resistance to six commonly-used drugs. Genes expressed during the mosquito stage of the parasite life-cycle are prominent among loci that show strong geographic differentiation. By continuing to enlarge this open data resource we aim to facilitate research into the evolutionary processes affecting malaria control and to accelerate development of the surveillance toolkit required for malaria elimination.

An open dataset of Plasmodium falciparum genome variation in 7,000 worldwide samples

MalariaGEN is a data-sharing network that enables groups around the world to work together on the genomic epidemiology of malaria. Here we describe a new release of curated genome variation data on 7,000 Plasmodium falciparum samples from MalariaGEN partner studies in 28 malaria-endemic countries. High-quality genotype calls on 3 million single nucleotide polymorphisms (SNPs) and short indels were produced using a standardised analysis pipeline. Copy number variants associated with drug resistance and structural variants that cause failure of rapid diagnostic tests were also analysed.  Almost all samples showed genetic evidence of resistance to at least one antimalarial drug, and some samples from Southeast Asia carried markers of resistance to six commonly-used drugs. Genes expressed during the mosquito stage of the parasite life-cycle are prominent among loci that show strong geographic differentiation. By continuing to enlarge this open data resource we aim to facilitate research into the evolutionary processes affecting malaria control and to accelerate development of the surveillance toolkit required for malaria elimination.

Artemisinin resistance surveillance in African Plasmodium falciparum isolates from imported malaria cases to Italy

BACKGROUND

Plasmodium falciparum (P. falciparum) malaria is a significant public health problem in returning travellers, and artemisinin combination therapy (ACT) remains the first choice for treatment. Several single nucleotide polymorphisms (SNPs) in the P. falciparum kelch 13 (Pfk13) gene have been associated with artemisinin (ART) resistance. Moreover, the increase in the P. falciparum plasmepsin 2 (Pfpm2) gene copy number was shown to be linked with reduced susceptibility of P. falciparum to piperaquine (PPQ), a partner drug in an ACT regimen. Active molecular surveillance for imported drug-resistant malaria parasites is a pivotal activity to provide adequate chemoprophylaxis and treatment guidelines.

METHODS

A retrospective study to review imported P. falciparum malaria in patients admitted to Spallanzani Institute between 2014 and 2015 was conducted. Information collected included clinic and epidemiological characteristics such as age, gender, country of origin, time since arrival to our country, travel history. All P.falciparum isolates were analysed for SNPs in the Pfk13 gene and for copy number variations in the Pfpm2 gene.

RESULTS

P. falciparum malaria was identified in 54 travellers. The mean age was 37 years, 44 were males. All cases were imported from non-EU countries. In the Pfk13 gene two mutations (R561R and F673L) were detected. Six P. falciparum isolates carried two copies of Pfpm2 gene, and one three copies, representing ≈16% of the analysed isolates.

CONCLUSIONS

None of the SNPs known to be associated with ART resistance were detected in the examined parasites. Our results provide evidence that Pfpm2 duplications (associated with piperaquine resistance) occur in Africa, emphasizing the necessity to better decode the genetic background associated with PPQ resistance. Further epidemiological investigations in Pfpm2 amplification along with mutations in the Pfk13 gene will be useful for developing and updating anti-malarial guidance in travellers.

Plasmodium falciparum kelch 13 Mutations, 9 Countries in Africa, 2014-2018

The spread of drug resistance to antimalarial treatments poses a serious public health risk globally. To combat this risk, molecular surveillance of drug resistance is imperative. We report the prevalence of mutations in the Plasmodium falciparum kelch 13 propeller domain associated with partial artemisinin resistance, which we determined by using Sanger sequencing samples from patients enrolled in therapeutic efficacy studies from 9 sub-Saharan countries during 2014-2018. Of the 2,865 samples successfully sequenced before treatment (day of enrollment) and on the day of treatment failure, 29 (1.0%) samples contained 11 unique nonsynonymous mutations and 83 (2.9%) samples contained 27 unique synonymous mutations. Two samples from Kenya contained the S522C mutation, which has been associated with delayed parasite clearance; however, no samples contained validated or candidate artemisinin-resistance mutations.

Protein Prenylation and Hsp40 in Thermotolerance of Plasmodium falciparum Malaria Parasites

During its complex life cycle, the malaria parasite survives dramatic environmental stresses, including large temperature shifts. Protein prenylation is required during asexual replication of Plasmodium falciparum, and the canonical heat shock protein 40 protein (HSP40; PF3D7_1437900) is posttranslationally modified with a 15-carbon farnesyl isoprenyl group. In other organisms, farnesylation of Hsp40 orthologs controls their localization and function in resisting environmental stress. In this work, we find that plastidial isopentenyl pyrophosphate (IPP) synthesis and protein farnesylation are required for malaria parasite survival after cold and heat shock. Furthermore, loss of HSP40 farnesylation alters its membrane attachment and interaction with proteins in essential pathways in the parasite. Together, this work reveals that farnesylation is essential for parasite survival during temperature stress. Farnesylation of HSP40 may promote thermotolerance by guiding distinct chaperone-client protein interactions.

Efficacy of dihydroartemisinin/piperaquine and artesunate monotherapy for the treatment of uncomplicated Plasmodium falciparum malaria in Central Vietnam

BACKGROUND

Artemisinin-based combination therapies (ACTs) have significantly contributed to reduce Plasmodium falciparum malaria burden in Vietnam, but their efficacy is challenged by treatment failure of dihydroartemisinin/piperaquine ACT in Southern provinces.

OBJECTIVES

To assess the efficacy of dihydroartemisinin/piperaquine for uncomplicated P. falciparum malaria in Gia Lai, Central Vietnam, and determine parasite resistance to artemisinin (ClinicalTrials.gov identifier NCT02604966).

METHODS

Sixty patients received either dihydroartemisinin/piperaquine (4 mg/kg/day, 3 days; n = 33) or artesunate monotherapy (4 mg/kg/day, 3 days; n = 27) followed by dihydroartemisinin/piperaquine (AS + DHA/PPQ). Clinical phenotypes were determined during a 42 day follow-up and analysed together with ex vivo susceptibility to antimalarials and molecular markers of drug resistance.

RESULTS

Day 3 positivity rate was significantly higher in the AS + DHA/PPQ arm compared with dihydroartemisinin/piperaquine (70.4% versus 39.4%, P = 0.016). Parasite clearance time was 95.2 h (AS + DHA/PPQ) versus 71.9 h (dihydroartemisinin/piperaquine, P = 0.063) and parasite clearance half-life was 7.4 h (AS + DHA/PPQ) versus 7.0 h (dihydroartemisinin/piperaquine, P = 0.140). Adequate clinical and parasitological response at Day 42 was 100% in both arms. By RT-qPCR, 36% (19/53) patients remained positive until Day 7. No recurrences were detected. kelch13 artemisinin resistance mutations were found in 87% (39/45) of isolates and 50% (20/40) were KEL1/C580Y. The piperaquine resistance marker plasmepsin-2 was duplicated in 10.4% (5/48). Isolates from Day 3-positive patients (n = 18) had higher ex vivo survival rates to artemisinin compounds (P < 0.048) and prevalence of kelch13 mutations (P = 0.005) than Day 3-negative patients (n = 5). The WHO definition of artemisinin resistance was fulfilled in 60% (24/40) of cases.

CONCLUSIONS

Although dihydroartemisinin/piperaquine remained effective to treat P. falciparum, the high Day 3 positivity rate and prevalence of KEL1 strains calls for continuous monitoring of dihydroartemisinin/piperaquine efficacy in Central Vietnam.

The Plasmodium falciparum Rh5 invasion protein complex reveals an excess of rare variant mutations

BACKGROUND

The invasion of the red blood cells by Plasmodium falciparum merozoites involves the interplay of several proteins that are also targets for vaccine development. The proteins PfRh5-PfRipr-PfCyRPA-Pfp113 assemble into a complex at the apical end of the merozoite and are together essential for erythrocyte invasion. They have also been shown to induce neutralizing antibodies and appear to be less polymorphic than other invasion-associated proteins, making them high priority blood-stage vaccine candidates. Using available whole genome sequencing data (WGS) and new capillary sequencing data (CS), this study describes the genetic polymorphism in the Rh5 complex in P. falciparum isolates obtained from Kilifi, Kenya.

METHODS

162 samples collected in 2013 and 2014 were genotyped by capillary sequencing (CS) and re-analysed WGS from 68 culture-adapted P. falciparum samples obtained from a drug trial conducted from 2005 to 2007. The frequency of polymorphisms in the merozoite invasion proteins, PfRh5, PfRipr, PfCyRPA and PfP113 were examined and where possible polymorphisms co-occurring in the same isolates.

RESULTS

From a total 70 variants, including 2 indels, 19 SNPs [27.1%] were identified by both CS and WGS, while an additional 15 [21.4%] and 36 [51.4%] SNPs were identified only by either CS or WGS, respectively. All the SNPs identified by CS were non-synonymous, whereas WGS identified 8 synonymous and 47 non-synonymous SNPs. CS identified indels in repeat regions in the p113 gene in codons 275 and 859 that were not identified in the WGS data. The minor allele frequencies of the SNPs ranged between 0.7 and 34.9% for WGS and 1.1-29.6% for CS. Collectively, 12 high frequency SNPs (> 5%) were identified: four in Rh5 codon 147, 148, 203 and 429, two in p113 at codons 7 and 267 and six in Ripr codons 190, 259, 524, 985, 1003 and 1039.

CONCLUSION

This study reveals that the majority of the polymorphisms are rare variants and confirms a low level of genetic polymorphisms in all proteins within the Rh5 complex.

Isolation of Mutants With Reduced Susceptibility to Piperaquine From a Mutator of the Rodent Malaria Parasite Plasmodium berghei

Elucidation of the mechanisms of drug resistance in malaria parasites is crucial for combatting the emergence and spread of resistant parasites, which can be achieved by tracing resistance-associated mutations and providing useful information for drug development. Previously, we produced a novel genetic tool, a Plasmodium berghei mutator (PbMut), whose base substitution rate is 36.5 times higher than that of wild-type parasites. Here, we report the isolation of a mutant with reduced susceptibility to piperaquine (PPQ) from PbMut under PPQ pressure by sequential nine-cycle screening and named it PbMut-PPQ-R-P9. The ED50 of PbMut-PPQ-R-P9 was 1.79 times higher than that of wild-type parasites, suggesting that its PPQ resistance is weak. In the 1st screen, recrudescence occurred in the mice infected with PbMut but not in those infected with wild-type parasites, suggesting earlier emergence of PPQ-resistant parasites from PbMut. Whole-genome sequence analysis of PbMut-PPQ-R-P9 clones revealed that eight nonsynonymous mutations were conserved in all clones, including N331I in PbCRT, the gene encoding chloroquine resistance transporter (CRT). The PbCRT(N331I) mutation already existed in the parasite population after the 2nd screen and was predominant in the population after the 8th screen. An artificially inserted PbCRT(N331I) mutation gave rise to reduced PPQ susceptibility in genome-edited parasites (PbCRT-N331I). The PPQ susceptibility and growth rates of PbCRT-N331I parasites were significantly lower than those of PbMut-PPQ-R-P9, implying that additional mutations in the PbMut-PPQ-R9 parasites could compensate for the fitness cost of the PbCRT(N331I) mutation and contribute to reduced PPQ susceptibility. In summary, PbMut could serve as a novel genetic tool for predicting gene mutations responsible for drug resistance. Further study on PbMut-PPQ-R-P9 could identify genetic changes that compensate for fitness costs owing to drug resistance acquisition.

Malaria in the 'Omics Era'

Genomics has revolutionised the study of the biology of parasitic diseases. The first Eukaryotic parasite to have its genome sequenced was the malaria parasite Plasmodium falciparum. Since then, Plasmodium genomics has continued to lead the way in the study of the genome biology of parasites, both in breadth-the number of Plasmodium species' genomes sequenced-and in depth-massive-scale genome re-sequencing of several key species. Here, we review some of the insights into the biology, evolution and population genetics of Plasmodium gained from genome sequencing, and look at potential new avenues in the future genome-scale study of its biology.

An integrated virtual screening and drug repurposing strategy for the discovery of new antimalarial drugs against Plasmodium falciparum phosphatidylinositol 3-kinase

The emergence and spread of drug resistance in Plasmodium falciparum, the parasite causing the most severe form of human malaria, is a major threat to malaria control and elimination programs around the globe. With P. falciparum having evolved widespread resistance against a number of previously widely used drugs, currently, artemisinin (ART) and its derivatives are the cornerstones of first-line treatments of uncomplicated malaria. However, growing incidences of ART failure reflect the spread of ART-resistant P. falciparum strains. Despite current efforts to understand the primary cause of ART resistance due to mutations in the Kelch 13 gene (PfK13), the mechanism underlying ART resistance is still not completely unclear and no feasible strategies to counteract the causes and thereby restoring the efficiency of ART have been developed. We use a polypharmacology approach to identify potential drugs that can be used for the novel purpose (target). Of note, we have designed a multimodal stratagem to identify approved drugs with a potential antimalarial activity using computational drug reprofiling. Our investigations suggest that oxetacaine, simvastatin, repaglinide, aclidinium, propafenone, and lovastatin could be repurposed for malaria control and prevention.

Orthosteric-allosteric dual inhibitors of PfHT1 as selective antimalarial agents

There is an urgent need for alternative antimalarials with the emergence of artemisinin-resistant malaria parasites. Blocking sugar uptake in Plasmodium falciparum by selectively inhibiting the hexose transporter P. falciparum hexose transporter 1 (PfHT1) kills the blood-stage parasites without affecting the host cells, making PfHT1 a promising therapeutic target. Here, we report the development of a series of small-molecule inhibitors that simultaneously target the orthosteric and the allosteric binding sites of PfHT1. These inhibitors all exhibit selective potency on the P. falciparum strains over human cell lines. Our findings establish the basis for the rational design of next-generation antimalarial drugs.

Artemisinin-resistant malaria parasites have emerged and have been spreading, posing a significant public health challenge. Antimalarial drugs with novel mechanisms of action are therefore urgently needed. In this report, we exploit a “selective starvation” strategy by inhibiting Plasmodium falciparum hexose transporter 1 (PfHT1), the sole hexose transporter in P. falciparum, over human glucose transporter 1 (hGLUT1), providing an alternative approach to fight against multidrug-resistant malaria parasites. The crystal structure of hGLUT3, which shares 80% sequence similarity with hGLUT1, was resolved in complex with C3361, a moderate PfHT1-specific inhibitor, at 2.3-Å resolution. Structural comparison between the present hGLUT3-C3361 and our previously reported PfHT1-C3361 confirmed the unique inhibitor binding-induced pocket in PfHT1. We then designed small molecules to simultaneously block the orthosteric and allosteric pockets of PfHT1. Through extensive structure–activity relationship studies, the TH-PF series was identified to selectively inhibit PfHT1 over hGLUT1 and potent against multiple strains of the blood-stage P. falciparum. Our findings shed light on the next-generation chemotherapeutics with a paradigm-shifting structure-based design strategy to simultaneously target the orthosteric and allosteric sites of a transporter.

Whole-genome analysis of Malawian Plasmodium falciparum isolates identifies possible targets of allele-specific immunity to clinical malaria

Individuals acquire immunity to clinical malaria after repeated Plasmodium falciparum infections. Immunity to disease is thought to reflect the acquisition of a repertoire of responses to multiple alleles in diverse parasite antigens. In previous studies, we identified polymorphic sites within individual antigens that are associated with parasite immune evasion by examining antigen allele dynamics in individuals followed longitudinally. Here we expand this approach by analyzing genome-wide polymorphisms using whole genome sequence data from 140 parasite isolates representing malaria cases from a longitudinal study in Malawi and identify 25 genes that encode possible targets of naturally acquired immunity that should be validated immunologically and further characterized for their potential as vaccine candidates.

Targeted deep amplicon sequencing of antimalarial resistance markers in Plasmodium falciparum isolates from Cameroon

BACKGROUND

Recent studies showed the first emergence of the R561H artemisinin-associated resistance marker in Africa, which highlights the importance of continued molecular surveillance to assess the selection and spread of this and other drug resistance markers in the region.

METHOD

In this study, we used targeted amplicon deep sequencing of 116 isolates collected in two areas of Cameroon to genotype the major drug resistance genes, k13, crt, mdr1, dhfr, and dhps, and the cytochrome b gene (cytb) in Plasmodium falciparum.

RESULTS

No confirmed or associated artemisinin resistance markers were observed in Pfk13. In comparison, both major and minor alleles associated with drug resistance were found in Pfcrt, Pfmdr1, Pfdhfr, and Pfdhps. Notably, a high frequency of other nonsynonymous mutations was observed across all the genes, except for Pfcytb, suggesting continued selection pressure.

CONCLUSIONS

The results from this study supported the continued use of artemisinin-based combination therapy and administration of sulfadoxine-pyrimethamine for intermittent preventive therapy in pregnant women, and for seasonal chemoprevention in these study sites in Cameroon.

Plasmodium falciparum Malaria Parasites in Ghana Show Signatures of Balancing Selection at Artemisinin Resistance Predisposing Background Genes

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

State of Artemisinin and Partner Drug Susceptibility in Plasmodium falciparum Clinical Isolates from Colombia

Delayed parasite clearance time observed in Southeast Asia provided the first evidence of Plasmodium falciparum resistance to artemisinins. The ex vivo ring-stage survival assay (RSA) mimics parasite exposure to pharmacologically relevant artemisinin concentrations. Mutations in the C-terminal propeller domain of the putative kelch protein Pf3D7_1343700 (K13) are associated with artemisinin resistance. Variations in the pfmdr1 gene are associated with reduced susceptibility to the artemisinin partner drugs mefloquine (MQ) and lumefantrine (LF). To clarify the unknown landscape of artemisinin resistance in Colombia, 71 patients with uncomplicated P. falciparum malaria were enrolled in a non-randomized observational study in three endemic localities in 2014-2015. Each patient's parasite isolate was assessed for ex vivo RSA, K13-propeller mutations, pfmdr1 copy number, and pfmdr1 mutations at codons 86, 184, 1034, 1042, and 1246, associated with reduced susceptibility, and 50% inhibitory concentration (IC₅₀) for other antimalarial drugs. Ex vivo RSAs were successful in 56% (40/71) of samples, and nine isolates showed survival rates > 1%. All isolates had wild-type K13-propeller sequences. All isolates harbored either of two pfmdr1 haplotypes, NFSDD (79.3%) and NFSDY (20.7%), and 7.1% of isolates had > 1 pfmdr1 gene. In vitro IC₅₀ assays showed that variable proportions of isolates had decreased susceptibility to chloroquine (52.4%, > 100 nM), amodiaquine (31.2%, > 30 nM), MQ (34.3%, > 30 nM), and LF (3.2%, > 10 nM). In this study, we report ex vivo RSA and K13 data on P. falciparum isolates from Colombia. The identification of isolates with increased ex vivo RSA rates in the absence of K13-propeller mutations and no positivity at day three requires further investigation.

Artemisinin-resistant K13 mutations rewire Plasmodium falciparum's intra-erythrocytic metabolic program to enhance survival

The emergence and spread of artemisinin resistance, driven by mutations in Plasmodium falciparum K13, has compromised antimalarial efficacy and threatens the global malaria elimination campaign. By applying systems-based quantitative transcriptomics, proteomics, and metabolomics to a panel of isogenic K13 mutant or wild-type P. falciparum lines, we provide evidence that K13 mutations alter multiple aspects of the parasite's intra-erythrocytic developmental program. These changes impact cell-cycle periodicity, the unfolded protein response, protein degradation, vesicular trafficking, and mitochondrial metabolism. K13-mediated artemisinin resistance in the Cambodian Cam3.II line was reversed by atovaquone, a mitochondrial electron transport chain inhibitor. These results suggest that mitochondrial processes including damage sensing and anti-oxidant properties might augment the ability of mutant K13 to protect P. falciparum against artemisinin action by helping these parasites undergo temporary quiescence and accelerated growth recovery post drug elimination.

Polymorphism Analysis of pfmdr1 and pfcrt from Plasmodium falciparum Isolates in Northwestern Nigeria Revealed the Major Markers Associated with Antimalarial Resistance

Suspicion of failure in the effectiveness of artemisinin-based combination therapies (currently the first-line treatment of malaria, worldwide) is leading to the unofficial use of alternative antimalarials, including chloroquine and sulfadoxine/pyrimethamine, across northern Nigeria. To facilitate evidence-based resistance management, antimalarial resistance mutations were investigated in Plasmodium falciparum multidrug resistance-1 (pfmdr1) and chloroquine resistance transporter (pfcrt), in isolates from Kano, northwestern Nigeria. Out of the 88 samples genotyped for pfmdr1N86Y mutation using PCR/restriction fragment length polymorphism, one sample contained the 86Y mutation (86Yfrequency = 1.14%). The analysis of 610 bp fragments of pfmdr1 from 16 isolates revealed two polymorphic sites and low haplotype diversity (Hd = 0.492), with only 86 Y mutations in one isolate, and 184 F replacements in five isolates (184Ffrequency = 31.25%). The analysis of 267 bp fragments of pfcrt isolates revealed high polymorphism (Hd = 0.719), with six haplotypes and seven non-synonymous polymorphic sites. Eleven isolates (61.11%) were chloroquine-resistant, CQR (C72V73I74E75T76 haplotype), two of which had an additional mutation, D57E. An additional sequence was CQR, but of the C72V73M74E75T76 haplotype, while the rest of the sequences (33.33%) were chloroquine susceptible (C72V73M74N75K76 haplotype). The findings of these well characterized resistance markers should be considered when designing resistance management strategies in the northwestern Nigeria.

Increase in Kelch 13 Polymorphisms in Plasmodium falciparum, Southern Rwanda

Artemisinin resistance in Plasmodium falciparum is associated with nonsynonymous mutations in the Kelch 13 (K13) propeller domain. We found that 12.1% (8/66) of clinical P. falciparum isolates from Huye district, Rwanda, exhibited K13 mutations, including R561H, a validated resistance marker. K13 mutations appear to be increasing in this region.

Emergence of artemisinin-resistant Plasmodium falciparum with kelch13 C580Y mutations on the island of New Guinea

The rapid and aggressive spread of artemisinin-resistant Plasmodium falciparum carrying the C580Y mutation in the kelch13 gene is a growing threat to malaria elimination in Southeast Asia, but there is no evidence of their spread to other regions. We conducted cross-sectional surveys in 2016 and 2017 at two clinics in Wewak, Papua New Guinea (PNG) where we identified three infections caused by C580Y mutants among 239 genotyped clinical samples. One of these mutants exhibited the highest survival rate (6.8%) among all parasites surveyed in ring-stage survival assays (RSA) for artemisinin. Analyses of kelch13 flanking regions, and comparisons of deep sequencing data from 389 clinical samples from PNG, Indonesian Papua and Western Cambodia, suggested an independent origin of the Wewak C580Y mutation, showing that the mutants possess several distinctive genetic features. Identity by descent (IBD) showed that multiple portions of the mutants' genomes share a common origin with parasites found in Indonesian Papua, comprising several mutations within genes previously associated with drug resistance, such as mdr1, ferredoxin, atg18 and pnp. These findings suggest that a P. falciparum lineage circulating on the island of New Guinea has gradually acquired a complex ensemble of variants, including kelch13 C580Y, which have affected the parasites' drug sensitivity. This worrying development reinforces the need for increased surveillance of the evolving parasite populations on the island, to contain the spread of resistance.

Genetic background and PfKelch13 affect artemisinin susceptibility of PfCoronin mutants in Plasmodium falciparum

Malaria continues to impose a significant health burden in the continent of Africa with 213 million cases in 2018 alone, representing 93% of cases worldwide. Because of high transmission of malaria within the continent, the selection pressures to develop drug resistance in African parasites are distinct compared to the rest of the world. In light of the spread of resistance to artemisinin conferred by the C580Y mutation in the PfKelch13 propeller domain in Southeast Asia, and its independent emergence in South America, it is important to study genetic determinants of resistance in the African context using African parasites. Through in vitro evolution of Senegalese parasites, we had previously generated the artemisinin-resistant parasites Pikine_R and Thiès_R and established pfcoronin mutations to be sufficient to confer artemisinin resistance in the standard ring-stage survival assay (RSA). In the current study, we used genetic analysis of revertants to demonstrate pfcoronin to be the major driver of elevated RSA in the artemisinin-resistant parasites Pikine_R and Thiès_R evolved in vitro. We interrogated the role of a second gene PF3D7_1433800, which also had mutations in both the Pikine_R and Thiès_R selected lines, but found no evidence of a contribution to reduced susceptibility in the RSA survival assay. Nevertheless, our genetic analysis demonstrates that parasite genetic background is important in the level of pfcoronin mediated RSA survival, and therefore we cannot rule out a role for PF3D7_1433800 in other genetic backgrounds. Finally, we tested the potential synergy between the mutations of pfcoronin and pfkelch13 through the generation of single and double mutants in the Pikine genetic background and found that the contribution of pfcoronin to reduced susceptibility is masked by the presence of pfkelch13. This phenomenon was also observed in the 3D7 background, suggesting that pfcoronin may mediate its effects via the same pathway as pfkelch13. Investigating the biology of proteins containing the beta-propeller domain could further elucidate the different pathways that the parasite could use to attain resistance.

Efficacy of Artemether-Lumefantrine on various Plasmodium falciparum Kelch 13 and Pfmdr1 genes isolated in Ghana

Introduction

Artemether-Lumefantrine (A-L) remains the drug of choice for the treatment of uncomplicated malaria in Ghana. However, the pharmaco-activity of A-L has not been assessed on various Plasmodium falciparum Kelch 13 and Pfmdr1 genes. Therefore, this study sought to determine the therapeutic efficacy of A-L on P. falciparum parasites isolated from Ghana.

Methods

The clinical study was done in Ga West Municipality, Ghana, where 78 uncomplicated malaria patients were recruited with prior consent. The patients were treated orally with A-L according to national treatment guidelines. Baseline parasitaemia was determined before treatment and 8-hourly parasitaemia posttreatment were determined till initial clearance of parasitaemia and at days 7, 14, 21, and 28. Kelch 13 and Pfmdr1 genes were genotyped by sequencing using baseline samples. Parasite clearance characteristics were determined using Parasite Clearance Estimator beta 0.9 application.

Results

Five Kelch 13 (F446I, S466N, R539I, A578S, and A676S) and three Pfmdr1 mutations (N86Y, Y184F and D1246Y) were identified in 78 infected samples. About 8% of the samples contained two Pfmdr1 double mutations (N86Y & D1246Y and Y184F & N86Y). Additionally, three samples (3.8%) were found to contain both Kelch 13 mutations and Pfmdr1 wild type genes. In all patients, parasitaemia persisted within the first 24 h of A-L therapy. However, at hour 40, only two patients were parasitaemic while all patients were aparasitaemic at hour 48. The genotypic profiles of the two persistent parasites at hour 40 were F446I and D1246Y, and R539I, Y184F, and N86Y. The slope half-life of the former was 6.4 h while the latter was 6.9 h and their respective PCT99 were 47.9 h and 49.2 h as well as a clearance rate constants of 0.109 and 0.092 respectively.

Conclusion

This study reports the effectiveness of A-L on various P. falciparum mutant alleles. However, continuous surveillance of Kelch 13 mutations and Pfmdr1 gene in Ghana and regular assessment of the therapeutic efficacy of A-L and other artemisinin derivatives is recommended.

Kelch 13-propeller polymorphisms in Plasmodium falciparum from Jazan region, southwest Saudi Arabia

BACKGROUND

Artemisinin-based combination therapy (ACT) is recommended at the initial phase for treatment of Plasmodium falciparum, to reduce morbidity and mortality in all countries where malaria is endemic. Polymorphism in portions of P. falciparum gene encoding kelch (K13)-propeller domains is associated with delayed parasite clearance after ACT. Of about 124 different non-synonymous mutations, 46 have been identified in Southeast Asia (SEA), 62 in sub-Saharan Africa (SSA) and 16 in both the regions. This is the first study designed to analyse the prevalence of polymorphism in the P. falciparum k13-propeller domain in the Jazan region of southwest Saudi Arabia, where malaria is endemic.

METHODS

One-hundred and forty P. falciparum samples were collected from Jazan region of southwest Saudi Arabia at three different times: 20 samples in 2011, 40 samples in 2016 and 80 samples in 2020 after the implementation of ACT. Plasmodium falciparum kelch13 (k13) gene DNA was extracted, amplified, sequenced, and analysed using a basic local alignment search tool (BLAST).

RESULTS

This study obtained 51 non-synonymous (NS) mutations in three time groups, divided as follows: 6 single nucleotide polymorphisms (SNPs) '11.8%' in samples collected in 2011 only, 3 (5.9%) in 2011and 2016, 5 (9.8%) in 2011 and 2020, 5 (9.8%) in 2016 only, 8 (15.7%) in 2016 and 2020, 14 (27.5%) in 2020 and 10 (19.6%) in all the groups. The BLAST revealed that the 2011 isolates were genetically closer to African isolates (53.3%) than Asian ones (46.7%). Interestingly, this proportion changed completely in 2020, to become closer to Asian isolates (81.6%) than to African ones (18.4%).

CONCLUSIONS

Despite the diversity of the identified mutations in the k13-propeller gene, these data did not report widespread artemisinin-resistant polymorphisms in the Jazan region where these samples were collected. Such a process would be expected to increase frequencies of mutations associated with the resistance of ACT.

Plasmodium berghei K13 Mutations Mediate In Vivo Artemisinin Resistance That Is Reversed by Proteasome Inhibition

Recent successes in malaria control have been seriously threatened by the emergence of Plasmodium falciparum parasite resistance to the frontline artemisinin drugs in Southeast Asia. P. falciparum artemisinin resistance is associated with mutations in the parasite K13 protein, which associates with a delay in the time required to clear the parasites upon drug treatment. Gene editing technologies have been used to validate the role of several candidate K13 mutations in mediating P. falciparum artemisinin resistance in vitro under laboratory conditions. Nonetheless, the causal role of these mutations under in vivo conditions has been a matter of debate. Here, we have used CRISPR/Cas9 gene editing to introduce K13 mutations associated with artemisinin resistance into the related rodent-infecting parasite, Plasmodium berghei. Phenotyping of these P. berghei K13 mutant parasites provides evidence of their role in mediating artemisinin resistance in vivo, which supports in vitro artemisinin resistance observations. However, we were unable to introduce some of the P. falciparum K13 mutations (C580Y and I543T) into the corresponding amino acid residues, while other introduced mutations (M476I and R539T equivalents) carried pronounced fitness costs. Our study provides evidence of a clear causal role of K13 mutations in modulating susceptibility to artemisinins in vitro and in vivo using the well-characterized P. berghei model. We also show that inhibition of the P. berghei proteasome offsets parasite resistance to artemisinins in these mutant lines.

The recent emergence of Plasmodium falciparum parasite resistance to the first line antimalarial drug artemisinin is of particular concern. Artemisinin resistance is primarily driven by mutations in the P. falciparum K13 protein, which enhance survival of early ring-stage parasites treated with the artemisinin active metabolite dihydroartemisinin in vitro and associate with delayed parasite clearance in vivo. However, association of K13 mutations with in vivo artemisinin resistance has been problematic due to the absence of a tractable model. Herein, we have employed CRISPR/Cas9 genome editing to engineer selected orthologous P. falciparum K13 mutations into the K13 gene of an artemisinin-sensitive Plasmodium berghei rodent model of malaria. Introduction of the orthologous P. falciparum K13 F446I, M476I, Y493H, and R539T mutations into P. berghei K13 yielded gene-edited parasites with reduced susceptibility to dihydroartemisinin in the standard 24-h in vitro assay and increased survival in an adapted in vitro ring-stage survival assay. Mutant P. berghei K13 parasites also displayed delayed clearance in vivo upon treatment with artesunate and achieved faster recrudescence upon treatment with artemisinin. Orthologous C580Y and I543T mutations could not be introduced into P. berghei, while the equivalents of the M476I and R539T mutations resulted in significant growth defects. Furthermore, a Plasmodium-selective proteasome inhibitor strongly synergized dihydroartemisinin action in these P. berghei K13 mutant lines, providing further evidence that the proteasome can be targeted to overcome artemisinin resistance. Taken together, our findings provide clear experimental evidence for the involvement of K13 polymorphisms in mediating susceptibility to artemisinins in vitro and, most importantly, under in vivo conditions.

Genetic Markers of Adaptation of Plasmodium falciparum to Transmission by American Vectors Identified in the Genomes of Parasites from Haiti and South America

The malaria parasite, Plasmodium falciparum, was introduced into Hispaniola and other regions of the Americas through the slave trade spanning the 16th through the 19th centuries. During this period, more than 12 million Africans were brought across the Atlantic to the Caribbean and other regions of the Americas. Since malaria is holoendemic in West Africa, a substantial percentage of these individuals carried the parasite. St. Domingue on Hispaniola, now modern-day Haiti, was a major port of disembarkation, and malaria is still actively transmitted there. We undertook a detailed study of the phylogenetics of the Haitian parasites and those from Colombia and Peru utilizing whole-genome sequencing. Principal-component and phylogenetic analyses, based upon single nucleotide polymorphisms (SNPs) in protein coding regions, indicate that, despite the potential for millions of introductions from Africa, the Haitian parasites share an ancestral relationship within a well-supported monophyletic clade with parasites from South America, while belonging to a distinct lineage. This result, in stark contrast to the historical record of parasite introductions, is best explained by a severe population bottleneck experienced by the parasites introduced into the Americas. Here, evidence is presented for targeted selection of rare African alleles in genes which are expressed in the mosquito stages of the parasite's life cycle. These genetic markers support the hypothesis that the severe population bottleneck was caused by the required adaptation of the parasite to transmission by new definitive hosts among the Anopheles (Nyssorhynchus) spp. found in the Caribbean and South America.IMPORTANCE Historical data suggest that millions of P. falciparum parasite lineages were introduced into the Americas during the trans-Atlantic slave trade, which would suggest a paraphyletic origin of the extant isolates in the Western Hemisphere. Our analyses of whole-genome variants show that the American parasites belong to a well-supported monophyletic clade. We hypothesize that the required adaptation to American vectors created a severe bottleneck, reducing the effective introduction to a few lineages. In support of this hypothesis, we discovered genes expressed in the mosquito stages of the life cycle that have alleles with multiple, high-frequency or fixed, nonsynonymous mutations in the American populations which are rarely found in African isolates. These alleles appear to be in gene products critical for transmission through the anopheline vector. Thus, these results may inform efforts to develop novel transmission-blocking vaccines by identifying parasite proteins functionally interacting with the vector that are important for successful transmission. Further, to the best of our knowledge, these are the first whole-genome data available from Haitian P. falciparum isolates. Defining the genome of these parasites provides genetic markers useful for mapping parasite populations and monitoring parasite movements/introductions.

Microbe hunting in the modern era: reflecting on a decade of microbial genomic epidemiology

Since the first recognition that infectious microbes serve as the causes of many human diseases, physicians and scientists have sought to understand and control their spread. For the past 150+ years, these ‘microbe hunters’ have learned to combine epidemiological information with knowledge of the infectious agent(s). In this essay, I reflect on the evolution of microbe hunting, beginning with the history of pre-germ theory epidemiological studies, through the microbiological and molecular eras. Now in the genomic age, modern-day microbe hunters are combining pathogen whole-genome sequencing with epidemiological data to enhance epidemiological investigations, advance our understanding of the natural history of pathogens and drivers of disease, and ultimately reshape our plans and priorities for global disease control and eradication. Indeed, as we have seen during the ongoing Covid-19 pandemic, the role of microbe hunters is now more important than ever. Despite the advances already made by microbial genomic epidemiology, the field is still maturing, with many more exciting developments on the horizon.

A cautionary note on the use of unsupervised machine learning algorithms to characterise malaria parasite population structure from genetic distance matrices

Genetic surveillance of malaria parasites supports malaria control programmes, treatment guidelines and elimination strategies. Surveillance studies often pose questions about malaria parasite ancestry (e.g. how antimalarial resistance has spread) and employ statistical methods that characterise parasite population structure. Many of the methods used to characterise structure are unsupervised machine learning algorithms which depend on a genetic distance matrix, notably principal coordinates analysis (PCoA) and hierarchical agglomerative clustering (HAC). PCoA and HAC are sensitive to both the definition of genetic distance and algorithmic specification. Importantly, neither algorithm infers malaria parasite ancestry. As such, PCoA and HAC can inform (e.g. via exploratory data visualisation and hypothesis generation), but not answer comprehensively, key questions about malaria parasite ancestry. We illustrate the sensitivity of PCoA and HAC using 393 Plasmodium falciparum whole genome sequences collected from Cambodia and neighbouring regions (where antimalarial resistance has emerged and spread recently) and we provide tentative guidance for the use and interpretation of PCoA and HAC in malaria parasite genetic epidemiology. This guidance includes a call for fully transparent and reproducible analysis pipelines that feature (i) a clearly outlined scientific question; (ii) a clear justification of analytical methods used to answer the scientific question along with discussion of any inferential limitations; (iii) publicly available genetic distance matrices when downstream analyses depend on them; and (iv) sensitivity analyses. To bridge the inferential disconnect between the output of non-inferential unsupervised learning algorithms and the scientific questions of interest, tailor-made statistical models are needed to infer malaria parasite ancestry. In the absence of such models speculative reasoning should feature only as discussion but not as results.

Plasmodium falciparum Knockout for the GPCR-Like PfSR25 Receptor Displays Greater Susceptibility to 1,2,3-Triazole Compounds That Block Malaria Parasite Development

The search for new compounds with antimalarial activity is urgent, as resistance to ones in the classical drug, has already been described in more than one continent. Compounds derived from 1,2,3-triazoles are effective against parasites and bacteria. Here, we evaluated the potential antimalarial activity against the human malaria parasite Plasmodium falciparum in a culture of fifty-four triazole compounds derived from 1H-and 2H-1,2,3-triazole. We identified thirty-one compounds with potential antimalarial activity at concentrations in the micromolar order (µM) and IC₅₀ values ranging from 2.80 µM (9) to 29.27 µM (21). Then, we selected some of these compounds to perform the same tests on the PfSR25- strain (knockout for P. falciparum G-protein coupled receptor-like, SR25). Our experiences with the PfSR25- strain showed that both compounds with higher antimalarial activity for the 3D7 strain and those with less activity resulted in lower IC₅₀ values for the knockout strain. The cytotoxicity of the compounds was evaluated in human renal embryonic cells (HEK 293), using MTT assays. This demonstrated that the compounds with the highest activity (9, 13, 19, 22, 24, 29), showed no toxicity at the tested concentrations.

Emergence and clonal expansion of in vitro artemisinin-resistant Plasmodium falciparum kelch13 R561H mutant parasites in Rwanda

Artemisinin resistance (delayed P. falciparum clearance following artemisinin-based combination therapy), is widespread across Southeast Asia but to date has not been reported in Africa^1–4. Here we genotyped the P. falciparum K13 (Pfkelch13) propeller domain, mutations in which can mediate artemisinin resistance^5,6, in pretreatment samples collected from recent dihydroarteminisin-piperaquine and artemether-lumefantrine efficacy trials in Rwanda⁷. While cure rates were >95% in both treatment arms, the Pfkelch13 R561H mutation was identified in 19 of 257 (7.4%) patients at Masaka. Phylogenetic analysis revealed the expansion of an indigenous R561H lineage. Gene editing confirmed that this mutation can drive artemisinin resistance in vitro. This study provides evidence for the de novo emergence of Pfkelch13-mediated artemisinin resistance in Rwanda, potentially compromising the continued success of antimalarial chemotherapy in Africa.

Identification in Rwanda of mutations in Plasmodium falciparum capable of conferring in vitro resistance to artemisinin, an essential medicine for the treatment of malaria, underscore the crucial need for surveillance in Africa to safeguard efficacy of life-saving therapies.

Characterization of putative drug resistant biomarkers in Plasmodium falciparum isolated from Ghanaian blood donors

BACKGROUND

Plasmodium falciparum parasites, which could harbour anti-malaria drug resistance genes, are commonly detected in blood donors in malaria-endemic areas. Notwithstanding, anti-malaria drug resistant biomarkers have not been characterized in blood donors with asymptomatic P. falciparum infection.

METHODS

A total of 771 blood donors were selected from five districts in the Greater Accra Region, Ghana. Each donor sample was screened with malaria rapid diagnostic test (RDT) kit and parasitaemia quantified microscopically. Dried blood spots from malaria positive samples were genotyped for P. falciparum chloroquine resistance transporter (Pfcrt), P. falciparum multi-drug resistance (Pfmdr1), P. falciparum dihydropteroate-synthetase (Pfdhps), P. falciparum dihydrofolate-reductase (Pfdhfr) and Kelch 13 propeller domain on chromosome 13 (Kelch 13) genes.

RESULTS

Of the 771 blood donors, 91 (11.8%) were positive by RDT. Analysis of sequence reads indicated successful genotyping of Pfcrt, Pfmdr1, Pfdhfr, Pfdhps and Kelch 13 genes in 84.6, 81.3, 86.8, 86.9 and 92.3% of the isolates respectively. Overall, 21 different mutant haplotypes were identified in 69 isolates (75.8%). In Pfcrt, CVIET haplotype was observed in 11.6% samples while in Pfmdr1, triple mutation (resulting in YFN haplotype) was detected in 8.1% of isolates. In Pfdhfr gene, triple mutation resulting in IRNI haplotype and in Pfdhps gene, quintuple mutation resulting in AGESS haplotype was identified in 17.7% parasite isolates. Finally, five non-synonymous Kelch 13 alleles were detected; C580Y (3.6%), P615L (4.8%), A578S (4.8%), I543V (2.4%) and A676S (1.2%) were detected.

CONCLUSION

Results obtained in this study indicated various frequencies of mutant alleles in Pfcrt, Pfmdr1, Pfdhfr, Pfdhps and Kelch 13 genes from P. falciparum infected blood donors. These alleles could reduce the efficacy of standard malaria treatment in transfusion-transmitted malaria cases. Incorporating malaria screening into donor screening protocol to defer infected donors is therefore recommended.

Identification of Mutations in Antimalarial Resistance Gene Kelch13 from plasmodium falciparum Isolates in Kano, Nigeria

Malaria control relies on first-line treatments that use artemisinin-combination therapies (ACT). Unfortunately, mutations in the plasmodium falciparum kelch13 gene result in delayed parasite clearance. Research on what is causing ACT failure is non-existent in northwestern Nigeria. Thus, the presence of mutations in kelch13 in P. falciparum isolates from Kano, Nigeria was investigated in this study. Microscopic examination of 154 blood samples obtained from patients revealed a high prevalence of P. falciparum infection (114 positive individuals, slide positivity rate = 74.03%). The 114 patients were administered Cartef® (ACT) and out of the 50 patients that returned for the 14-day follow up, 11 were positive for P. falciparum (slide positivity rate = 22%). On day 0, 80 samples out of 114 and 11 samples on day 14 (91 out of 125 microscopy-positive samples) were positive with Plasmodium according to the PCR of cytochrome oxidase I, which corresponds to 72.8%. A fragment of the kelch13 gene encompassing the propeller domains was sequenced in 49 samples, alongside samples of the susceptible strain pf_3D7. Low polymorphism was observed, suggesting a lack of selection on this gene, and only six mutations (Glu433Gly, Phe434Ile, Phe434Ser, Ile684Asn, Ile684Thr and Glu688Lys) were found. The epidemiologic impact of these mutations and their potential role in ACT resistance needs to be investigated further.

Different Plasmodium falciparum clearance times in two Malian villages following artesunate monotherapy

•

High prevalence of residual parasitemia at day 3 post-artesunate monotherapy treatment using qPCR while no parasites were detected by microscopy at the same timepoint.

•

A longer parasite clearance time observed in a Malian village.

•

Artesunate treatment is still efficacious on Plasmodium falciparum in Mali.

Background

Artemisinin resistance described as increased parasite clearance time (PCT) is rare in Africa. More sensitive methods such as qPCR might better characterize the clearance phenotype in sub-Saharan Africa.

Methods

PCT is explored in Mali using light microscopy and qPCR after artesunate for uncomplicated malaria. In two villages, patients were followed for 28 days. Blood smears and spots were collected respectively for microscopy and qPCR. Parasitemia slope half-life was calculated after microscopy. Patient residual parasitemia were measured by qPCR.

Results

Uncorrected adequate clinical and parasitological responses (ACPR) observed in Faladje and Bougoula-Hameau were 78% and 92%, respectively (p = 0.01). This reached 100% for both after molecular correction.

Proportions of 24H microscopy positive patients in Faladje and Bougoula-Hameau were 97.2% and 72%, respectively (p < 0.0001).

Slope half-life was 2.8 h in Faladje vs 2H in Bougoula-Hameau (p < 0.001) and

Proportions of 72H patients with residual parasitemia were 68.5% and 40% in Faladje and Bougoula-Hameau, respectively (p = 0.003). The mean residual parasitemia was 2.9 in Faladje vs. 0.008 in Bougoula-Hameau (p = 0.002).

Although artesunate is efficacious in Mali, the longer parasite clearance time with submicroscopic parasitemia observed may represent early signs of developing P. falciparum resistance to artemisinins.

Detecting geospatial patterns of Plasmodium falciparum parasite migration in Cambodia using optimized estimated effective migration surfaces

BACKGROUND

Understanding the genetic structure of natural populations provides insight into the demographic and adaptive processes that have affected those populations. Such information, particularly when integrated with geospatial data, can have translational applications for a variety of fields, including public health. Estimated effective migration surfaces (EEMS) is an approach that allows visualization of the spatial patterns in genomic data to understand population structure and migration. In this study, we developed a workflow to optimize the resolution of spatial grids used to generate EEMS migration maps and applied this optimized workflow to estimate migration of Plasmodium falciparum in Cambodia and bordering regions of Thailand and Vietnam.

METHODS

The optimal density of EEMS grids was determined based on a new workflow created using density clustering to define genomic clusters and the spatial distance between genomic clusters. Topological skeletons were used to capture the spatial distribution for each genomic cluster and to determine the EEMS grid density; i.e., both genomic and spatial clustering were used to guide the optimization of EEMS grids. Model accuracy for migration estimates using the optimized workflow was tested and compared to grid resolutions selected without the optimized workflow. As a test case, the optimized workflow was applied to genomic data generated from P. falciparum sampled in Cambodia and bordering regions, and migration maps were compared to estimates of malaria endemicity, as well as geographic properties of the study area, as a means of validating observed migration patterns.

RESULTS

Optimized grids displayed both high model accuracy and reduced computing time compared to grid densities selected in an unguided manner. In addition, EEMS migration maps generated for P. falciparum using the optimized grid corresponded to estimates of malaria endemicity and geographic properties of the study region that might be expected to impact malaria parasite migration, supporting the validity of the observed migration patterns.

CONCLUSIONS

Optimized grids reduce spatial uncertainty in the EEMS contours that can result from user-defined parameters, such as the resolution of the spatial grid used in the model. This workflow will be useful to a broad range of EEMS users as it can be applied to analyses involving other organisms of interest and geographic areas.

Insights into the intracellular localization, protein associations and artemisinin resistance properties of Plasmodium falciparum K13

The emergence of artemisinin (ART) resistance in Plasmodium falciparum intra-erythrocytic parasites has led to increasing treatment failure rates with first-line ART-based combination therapies in Southeast Asia. Decreased parasite susceptibility is caused by K13 mutations, which are associated clinically with delayed parasite clearance in patients and in vitro with an enhanced ability of ring-stage parasites to survive brief exposure to the active ART metabolite dihydroartemisinin. Herein, we describe a panel of K13-specific monoclonal antibodies and gene-edited parasite lines co-expressing epitope-tagged versions of K13 in trans. By applying an analytical quantitative imaging pipeline, we localize K13 to the parasite endoplasmic reticulum, Rab-positive vesicles, and sites adjacent to cytostomes. These latter structures form at the parasite plasma membrane and traffic hemoglobin to the digestive vacuole wherein artemisinin-activating heme moieties are released. We also provide evidence of K13 partially localizing near the parasite mitochondria upon treatment with dihydroartemisinin. Immunoprecipitation data generated with K13-specific monoclonal antibodies identify multiple putative K13-associated proteins, including endoplasmic reticulum-resident molecules, mitochondrial proteins, and Rab GTPases, in both K13 mutant and wild-type isogenic lines. We also find that mutant K13-mediated resistance is reversed upon co-expression of wild-type or mutant K13. These data help define the biological properties of K13 and its role in mediating P. falciparum resistance to ART treatment.

Amplicon deep sequencing of kelch13 in Plasmodium falciparum isolates from Senegal

BACKGROUND

In 2006, the Senegalese National Malaria Control Programme recommended artemisinin-based combination therapy (ACT) with artemether-lumefantrine as the first-line treatment for uncomplicated Plasmodium falciparum malaria. To date, multiple mutations associated with artemisinin delayed parasite clearance have been described in Southeast Asia in the Pfk13 gene, such as Y493H, R539T, I543T and C580Y. Even though ACT remains clinically and parasitologically efficacious in Senegal, the spread of resistance is possible as shown by the earlier emergence of resistance to chloroquine in Southeast Asia that subsequently spread to Africa. Therefore, surveillance of artemisinin resistance in malaria endemic regions is crucial and requires the implementation of sensitive tools, such as next-generation sequencing (NGS) which can detect novel mutations at low frequency.

METHODS

Here, an amplicon sequencing approach was used to identify mutations in the Pfk13 gene in eighty-one P. falciparum isolates collected from three different regions of Senegal.

RESULTS

In total, 10 SNPs around the propeller domain were identified; one synonymous SNP and nine non-synonymous SNPs, and two insertions. Three of these SNPs (T478T, A578S and V637I) were located in the propeller domain. A578S, is the most frequent mutation observed in Africa, but has not previously been reported in Senegal. A previous study has suggested that A578S could disrupt the function of the Pfk13 propeller region.

CONCLUSION

As the genetic basis of possible artemisinin resistance may be distinct in Africa and Southeast Asia, further studies are necessary to assess the new SNPs reported in this study.

Mapping the travel patterns of people with malaria in Bangladesh

BACKGROUND

Spread of malaria and antimalarial resistance through human movement present major threats to current goals to eliminate the disease. Bordering the Greater Mekong Subregion, southeast Bangladesh is a potentially important route of spread to India and beyond, but information on travel patterns in this area are lacking.

METHODS

Using a standardised short survey tool, 2090 patients with malaria were interviewed at 57 study sites in 2015-2016 about their demographics and travel patterns in the preceding 2 months.

RESULTS

Most travel was in the south of the study region between Cox's Bazar district (coastal region) to forested areas in Bandarban (31% by days and 45% by nights), forming a source-sink route. Less than 1% of travel reported was between the north and south forested areas of the study area. Farmers (21%) and students (19%) were the top two occupations recorded, with 67 and 47% reporting travel to the forest respectively. Males aged 25-49 years accounted for 43% of cases visiting forests but only 24% of the study population. Children did not travel. Women, forest dwellers and farmers did not travel beyond union boundaries. Military personnel travelled the furthest especially to remote forested areas.

CONCLUSIONS

The approach demonstrated here provides a framework for identifying key traveller groups and their origins and destinations of travel in combination with knowledge of local epidemiology to inform malaria control and elimination efforts. Working with the NMEP, the findings were used to derive a set of policy recommendations to guide targeting of interventions for elimination.

Uptake of Plasmodium falciparum Gametocytes During Mosquito Bloodmeal by Direct and Membrane Feeding

Plasmodium falciparum remains one of the leading causes of child mortality, and nearly half of the world's population is at risk of contracting malaria. While pathogenesis results from replication of asexual forms in human red blood cells, it is the sexually differentiated forms, gametocytes, which are responsible for the spread of the disease. For transmission to succeed, both mature male and female gametocytes must be taken up by a female Anopheles mosquito during its blood meal for subsequent differentiation into gametes and mating inside the mosquito gut. Observed circulating numbers of gametocytes in the human host are often surprisingly low. A pre-fertilization behavior, such as skin sequestration, has been hypothesized to explain the efficiency of human-to-mosquito transmission but has not been sufficiently tested due to a lack of appropriate tools. In this study, we describe the optimization of a qPCR tool that enables the relative quantification of gametocytes within very small input samples. Such a tool allows for the quantification of gametocytes in different compartments of the host and the vector that could potentially unravel mechanisms that enable highly efficient malaria transmission. We demonstrate the use of our gametocyte quantification method in mosquito blood meals from both direct skin feeding on Plasmodium gametocyte carriers and standard membrane feeding assay. Relative gametocyte abundance was not different between mosquitoes fed through a membrane or directly on the skin suggesting that there is no systematic enrichment of gametocytes picked up in the skin.

Role of Plasmodium falciparum Kelch 13 Protein Mutations in P. falciparum Populations from Northeastern Myanmar in Mediating Artemisinin Resistance

Artemisinin resistance has emerged in Southeast Asia, endangering the substantial progress in malaria elimination worldwide. It is associated with mutations in the PfK13 protein, but how PfK13 mediates artemisinin resistance is not completely understood. Here we used a new antibody against PfK13 to show that the PfK13 protein is expressed in all stages of the asexual intraerythrocytic cycle as well as in gametocytes and is partially localized in the endoplasmic reticulum. By introducing four PfK13 mutations into the 3D7 strain and reverting these mutations in field parasite isolates, we determined the impacts of these mutations identified in the parasite populations from northern Myanmar on the ring stage using the in vitro ring survival assay. The introduction of the N458Y mutation into the 3D7 background significantly increased the survival rates of the ring-stage parasites but at the cost of the reduced fitness of the parasites. Introduction of the F446I mutation, the most prevalent PfK13 mutation in northern Myanmar, did not result in a significant increase in ring-stage survival after exposure to dihydroartemisinin (DHA), but these parasites showed extended ring-stage development. Further, parasites with the F446I mutation showed only a marginal loss of fitness, partially explaining its high frequency in northern Myanmar. Conversely, reverting all these mutations, except for the C469Y mutation, back to their respective wild types reduced the ring-stage survival of these isolates in response to in vitro DHA treatment.

Mutations in the Plasmodium falciparum Kelch 13 (PfK13) protein are associated with artemisinin resistance. PfK13 is essential for asexual erythrocytic development, but its function is not known. We tagged the PfK13 protein with green fluorescent protein in P. falciparum to study its expression and localization in asexual and sexual stages. We used a new antibody against PfK13 to show that the PfK13 protein is expressed ubiquitously in both asexual erythrocytic stages and gametocytes and is localized in punctate structures, partially overlapping an endoplasmic reticulum marker. We introduced into the 3D7 strain four PfK13 mutations (F446I, N458Y, C469Y, and F495L) identified in parasites from the China-Myanmar border area and characterized the in vitro artemisinin response phenotypes of the mutants. We found that all the parasites with the introduced PfK13 mutations showed higher survival rates in the ring-stage survival assay (RSA) than the wild-type (WT) control, but only parasites with N458Y displayed a significantly higher RSA value (26.3%) than the WT control. After these PfK13 mutations were reverted back to the WT in field parasite isolates, all revertant parasites except those with the C469Y mutation showed significantly lower RSA values than their respective parental isolates. Although the 3D7 parasites with introduced F446I, the predominant PfK13 mutation in northern Myanmar, did not show significantly higher RSA values than the WT, they had prolonged ring-stage development and showed very little fitness cost in in vitro culture competition assays. In comparison, parasites with the N458Y mutations also had a prolonged ring stage and showed upregulated resistance pathways in response to artemisinin, but this mutation produced a significant fitness cost, potentially leading to their lower prevalence in the Greater Mekong subregion.

Falciparum malaria from coastal Tanzania and Zanzibar remains highly connected despite effective control efforts on the archipelago

BACKGROUND

Tanzania's Zanzibar archipelago has made significant gains in malaria control over the last decade and is a target for malaria elimination. Despite consistent implementation of effective tools since 2002, elimination has not been achieved. Importation of parasites from outside of the archipelago is thought to be an important cause of malaria's persistence, but this paradigm has not been studied using modern genetic tools.

METHODS

Whole-genome sequencing (WGS) was used to investigate the impact of importation, employing population genetic analyses of Plasmodium falciparum isolates from both the archipelago and mainland Tanzania. Ancestry, levels of genetic diversity and differentiation, patterns of relatedness, and patterns of selection between these two populations were assessed by leveraging recent advances in deconvolution of genomes from polyclonal malaria infections.

RESULTS

Significant decreases in the effective population sizes were inferred in both populations that coincide with a period of decreasing malaria transmission in Tanzania. Identity by descent analysis showed that parasites in the two populations shared long segments of their genomes, on the order of 5 cM, suggesting shared ancestry within the last 10 generations. Even with limited sampling, two of isolates between the mainland and Zanzibar were identified that are related at the expected level of half-siblings, consistent with recent importation.

CONCLUSIONS

These findings suggest that importation plays an important role for malaria incidence on Zanzibar and demonstrate the value of genomic approaches for identifying corridors of parasite movement to the island.

Strains used in whole organism Plasmodium falciparum vaccine trials differ in genome structure, sequence, and immunogenic potential

BACKGROUND

Plasmodium falciparum (Pf) whole-organism sporozoite vaccines have been shown to provide significant protection against controlled human malaria infection (CHMI) in clinical trials. Initial CHMI studies showed significantly higher durable protection against homologous than heterologous strains, suggesting the presence of strain-specific vaccine-induced protection. However, interpretation of these results and understanding of their relevance to vaccine efficacy have been hampered by the lack of knowledge on genetic differences between vaccine and CHMI strains, and how these strains are related to parasites in malaria endemic regions.

METHODS

Whole genome sequencing using long-read (Pacific Biosciences) and short-read (Illumina) sequencing platforms was conducted to generate de novo genome assemblies for the vaccine strain, NF54, and for strains used in heterologous CHMI (7G8 from Brazil, NF166.C8 from Guinea, and NF135.C10 from Cambodia). The assemblies were used to characterize sequences in each strain relative to the reference 3D7 (a clone of NF54) genome. Strains were compared to each other and to a collection of clinical isolates (sequenced as part of this study or from public repositories) from South America, sub-Saharan Africa, and Southeast Asia.

RESULTS

While few variants were detected between 3D7 and NF54, we identified tens of thousands of variants between NF54 and the three heterologous strains. These variants include SNPs, indels, and small structural variants that fall in regulatory and immunologically important regions, including transcription factors (such as PfAP2-L and PfAP2-G) and pre-erythrocytic antigens that may be key for sporozoite vaccine-induced protection. Additionally, these variants directly contributed to diversity in immunologically important regions of the genomes as detected through in silico CD8+ T cell epitope predictions. Of all heterologous strains, NF135.C10 had the highest number of unique predicted epitope sequences when compared to NF54. Comparison to global clinical isolates revealed that these four strains are representative of their geographic origin despite long-term culture adaptation; of note, NF135.C10 is from an admixed population, and not part of recently formed subpopulations resistant to artemisinin-based therapies present in the Greater Mekong Sub-region.

CONCLUSIONS

These results will assist in the interpretation of vaccine efficacy of whole-organism vaccines against homologous and heterologous CHMI.