The newly discovered role of endocytosis in artemisinin resistance

Temporal and spatial dynamics of Plasmodium falciparum clonal lineages in Guyana

Plasmodium parasites, the causal agents of malaria, are eukaryotic organisms that obligately undergo sexual recombination within mosquitoes. In low transmission settings, parasites recombine with themselves, and the clonal lineage is propagated rather than broken up by outcrossing. We investigated whether stochastic/neutral factors drive the persistence and abundance of Plasmodium falciparum clonal lineages in Guyana, a country with relatively low malaria transmission, but the only setting in the Americas in which an important artemisinin resistance mutation (pfk13 C580Y) has been observed. We performed whole genome sequencing on 1,727 Plasmodium falciparum samples collected from infected patients across a five-year period (2016-2021). We characterized the relatedness between each pair of monoclonal infections (n = 1,409) through estimation of identity-by-descent (IBD) and also typed each sample for known or candidate drug resistance mutations. A total of 160 multi-isolate clones (mean IBD ≥ 0.90) were circulating in Guyana during the study period, comprising 13 highly related clusters (mean IBD ≥ 0.40). In the five-year study period, we observed a decrease in frequency of a mutation associated with artemisinin partner drug (piperaquine) resistance (pfcrt C350R) and limited co-occurence of pfcrt C350R with duplications of plasmepsin 2/3, an epistatic interaction associated with piperaquine resistance. We additionally observed 61 nonsynonymous substitutions that increased markedly in frequency over the study period as well as a novel pfk13 mutation (G718S). However, P. falciparum clonal dynamics in Guyana appear to be largely driven by stochastic factors, in contrast to other geographic regions, given that clones carrying drug resistance polymorphisms do not demonstrate enhanced persistence or higher abundance than clones carrying polymorphisms of comparable frequency that are unrelated to resistance. The use of multiple artemisinin combination therapies in Guyana may have contributed to the disappearance of the pfk13 C580Y mutation.

Artemisinin resistance mutations in Pfcoronin impede hemoglobin uptake

Artemisinin (ART) combination therapies have been critical in reducing malaria morbidity and mortality, but these important drugs are threatened by growing resistance associated with mutations in Pfcoronin and Pfkelch13 . Here, we describe the mechanism of Pfcoronin -mediated ART resistance. Pf Coronin interacts with Pf Actin and localizes to the parasite plasma membrane (PPM), the digestive vacuole (DV) membrane, and membrane of a newly identified preDV compartment-all structures involved in the trafficking of hemoglobin from the RBC for degradation in the DV. Pfcoronin mutations alter Pf Actin homeostasis and impair the development and morphology of the preDV. Ultimately, these changes are associated with decreased uptake of red blood cell cytosolic contents by ring-stage Plasmodium falciparum . Previous work has identified decreased hemoglobin uptake as the mechanism of Pfkelch 13-mediated ART resistance. This work demonstrates that Pf Coronin appears to act via a parallel pathway. For both Pfkelch13 -mediated and Pfcoronin -mediated ART resistance, we hypothesize that the decreased hemoglobin uptake in ring stage parasites results in less heme-based activation of the artemisinin endoperoxide ring and reduced cytocidal activity. This study deepens our understanding of ART resistance, as well as hemoglobin uptake and development of the DV in early-stage parasites.

Whole genome sequencing identifies novel mutations in malaria parasites resistant to artesunate (ATN) and to ATN + mefloquine combination

Introduction

The global evolution of resistance to Artemisinin-based Combination Therapies (ACTs) by malaria parasites, will severely undermine our ability to control this devastating disease.

Methods

Here, we have used whole genome sequencing to characterize the genetic variation in the experimentally evolved Plasmodium chabaudi parasite clone AS-ATNMF1, which is resistant to artesunate + mefloquine.

Results and discussion

Five novel single nucleotide polymorphisms (SNPs) were identified, one of which was a previously undescribed E738K mutation in a 26S proteasome subunit that was selected for under artesunate pressure (in AS-ATN) and retained in AS-ATNMF1. The wild type and mutated three-dimensional (3D) structure models and molecular dynamics simulations of the P. falciparum 26S proteasome subunit Rpn2 suggested that the E738K mutation could change the toroidal proteasome/cyclosome domain organization and change the recognition of ubiquitinated proteins. The mutation in the 26S proteasome subunit may therefore contribute to altering oxidation-dependent ubiquitination of the MDR-1 and/or K13 proteins and/or other targets, resulting in changes in protein turnover. In light of the alarming increase in resistance to artemisin derivatives and ACT partner drugs in natural parasite populations, our results shed new light on the biology of resistance and provide information on novel molecular markers of resistance that may be tested (and potentially validated) in the field.

Temporal and spatial dynamics of Plasmodium falciparum clonal lineages in Guyana

Plasmodium parasites, the causal agents of malaria, are eukaryotic organisms that obligately undergo sexual recombination within mosquitoes. However, in low transmission settings where most mosquitoes become infected with only a single parasite clone, parasites recombine with themselves, and the clonal lineage is propagated rather than broken up by outcrossing. We investigated whether stochastic/neutral factors drive the persistence and abundance of Plasmodium falciparum clonal lineages in Guyana, a country with relatively low malaria transmission, but the only setting in the Americas in which an important artemisinin resistance mutation (pfk13 C580Y) has been observed. To investigate whether this clonality was potentially associated with the persistence and spatial spread of the mutation, we performed whole genome sequencing on 1,727 Plasmodium falciparum samples collected from infected patients across a five-year period (2016-2021). We characterized the relatedness between each pair of monoclonal infections (n=1,409) through estimation of identity by descent (IBD) and also typed each sample for known or candidate drug resistance mutations. A total of 160 clones (mean IBD ≥ 0.90) were circulating in Guyana during the study period, comprising 13 highly related clusters (mean IBD ≥ 0.40). In the five-year study period, we observed a decrease in frequency of a mutation associated with artemisinin partner drug (piperaquine) resistance (pfcrt C350R) and limited co-occurence of pfcrt C350R with duplications of plasmepsin 2/3, an epistatic interaction associated with piperaquine resistance. We additionally report polymorphisms exhibiting evidence of selection for drug resistance or other phenotypes and reported a novel pfk13 mutation (G718S) as well as 61 nonsynonymous substitutions that increased markedly in frequency. However, P. falciparum clonal dynamics in Guyana appear to be largely driven by stochastic factors, in contrast to other geographic regions. The use of multiple artemisinin combination therapies in Guyana may have contributed to the disappearance of the pfk13 C580Y mutation.

Efficacies and ADME properties of redox active methylene blue and phenoxazine analogues for use in new antimalarial triple drug combinations with amino-artemisinins

Efforts to develop new artemisinin triple combination therapies effective against artemisinin-tolerant strains of Plasmodium falciparum based on rational combinations comprising artemisone or other amino-artemisinins, a redox active drug and a third drug with a different mode of action have now been extended to evaluation of three potential redox partners. These are the diethyl analogue AD01 of methylene blue (MB), the benzo [α]phenoxazine PhX6, and the thiosemicarbazone DpNEt. IC50 values in vitro against CQ-sensitive and resistant P. falciparum strains ranged from 11.9 nM for AD01-41.8 nM for PhX6. PhX6 possessed the most favourable pharmacokinetic (PK) profile: intrinsic clearance rate CLint was 21.47 ± 1.76 mL/min/kg, bioavailability was 60% and half-life was 7.96 h. AD01 presented weaker, but manageable pharmacokinetic properties with a rapid CLint of 74.41 ± 6.68 mL/min/kg leading to a half-life of 2.51 ± 0.07 h and bioavailability of 15%. DpNEt exhibited a half-life of 1.12 h and bioavailability of 8%, data which discourage its further examination, despite a low CLint of 10.20 mL/min/kg and a high Cmax of 6.32 µM. Efficacies of AD01 and PhX6 were enhanced synergistically when each was paired with artemisone against asexual blood stages of P. falciparum NF54 in vitro. The favourable pharmacokinetics of PhX6 indicate this is the best partner among the compounds examined thus far for artemisone. Future work will focus on extending the drug combination studies to artemiside in vitro, and conducting efficacy studies in vivo for artemisone with each of PhX6 and the related benzo[α]phenoxazine SSJ-183.

The many paths to artemisinin resistance in Plasmodium falciparum

Emerging resistance against artemisinin (ART) poses a major challenge in controlling malaria. Parasites with mutations in PfKelch13, the major marker for ART resistance, are known to reduce hemoglobin endocytosis, induce unfolded protein response (UPR), elevate phosphatidylinositol-3-phosphate (PI3P) levels, and stimulate autophagy. Nonetheless, PfKelch13-independent resistance is also reported, indicating extensive complementation by reconfiguration in the parasite metabolome and transcriptome. These findings implicate that there may not be a single 'universal identifier' of ART resistance. This review sheds light on the molecular, transcriptional, and metabolic pathways associated with ART resistance, while also highlighting the interplay between cellular heterogeneity, environmental stress, and ART sensitivity.

Malaria Molecular Surveillance in the Peruvian Amazon with a Novel Highly Multiplexed Plasmodium falciparum AmpliSeq Assay

Molecular surveillance for malaria has great potential to support national malaria control programs (NMCPs). To bridge the gap between research and implementation, several applications (use cases) have been identified to align research, technology development, and public health efforts. For implementation at NMCPs, there is an urgent need for feasible and cost-effective tools. We designed a new highly multiplexed deep sequencing assay (Pf AmpliSeq), which is compatible with benchtop sequencers, that allows high-accuracy sequencing with higher coverage and lower cost than whole-genome sequencing (WGS), targeting genomic regions of interest. The novelty of the assay is its high number of targets multiplexed into one easy workflow, combining population genetic markers with 13 nearly full-length resistance genes, which is applicable for many different use cases. We provide the first proof of principle for hrp2 and hrp3 deletion detection using amplicon sequencing. Initial sequence data processing can be performed automatically, and subsequent variant analysis requires minimal bioinformatic skills using any tabulated data analysis program. The assay was validated using a retrospective sample collection (n = 254) from the Peruvian Amazon between 2003 and 2018. By combining phenotypic markers and a within-country 28-single-nucleotide-polymorphism (SNP) barcode, we were able to distinguish different lineages with multiple resistance haplotypes (in dhfr, dhps, crt and mdr1) and hrp2 and hrp3 deletions, which have been increasing in recent years. We found no evidence to suggest the emergence of artemisinin (ART) resistance in Peru. These findings indicate a parasite population that is under drug pressure but is susceptible to current antimalarials and demonstrate the added value of a highly multiplexed molecular tool to inform malaria strategies and surveillance systems. IMPORTANCE While the power of next-generation sequencing technologies to inform and guide malaria control programs has become broadly recognized, the integration of genomic data for operational incorporation into malaria surveillance remains a challenge in most countries where malaria is endemic. The main obstacles include limited infrastructure, limited access to high-throughput sequencing facilities, and the need for local capacity to run an in-country analysis of genomes at a large-enough scale to be informative for surveillance. In addition, there is a lack of standardized laboratory protocols and automated analysis pipelines to generate reproducible and timely results useful for relevant stakeholders. With our standardized laboratory and bioinformatic workflow, malaria genetic surveillance data can be readily generated by surveillance researchers and malaria control programs in countries of endemicity, increasing ownership and ensuring timely results for informed decision- and policy-making.

A Plasmodium falciparum RING Finger E3 Ubiquitin Ligase Modifies the Roles of PfMDR1 and PfCRT in Parasite Drug Responses

Protein ubiquitination is an important posttranslational regulation mechanism that mediates Plasmodium development and modifies parasite responses to antimalarial drugs. Although mutations in several parasite ubiquitination enzymes have been linked to increased drug tolerance, the molecular mechanisms by which ubiquitination pathways mediate these parasite responses remain largely unknown. Here, we investigate the roles of a Plasmodium falciparum ring finger ubiquitin ligase (PfRFUL) in parasite development and in responses to antimalarial drugs. We engineered a transgenic parasite having the Pfrful gene tagged with an HA-2A-NeoR-glmS sequence to knockdown (KD) Pfrful expression using glucosamine (GlcN). A Western blot analysis of the proteins from GlcN-treated pSLI-HA-NeoR-glmS-tagged (PfRFULg) parasites, relative to their wild-type (Dd2) controls, showed changes in the ubiquitination of numerous proteins. PfRFUL KD rendered the parasites more sensitive to multiple antimalarial drugs, including mefloquine, piperaquine, amodiaquine, and dihydroartemisinin. PfRFUL KD also decreased the protein level of the P. falciparum multiple drug resistance 1 protein (PfMDR1) and altered the ratio of two bands of the P. falciparum chloroquine resistance transporter (PfCRT), suggesting contributions to the changed drug responses by the altered ubiquitination of these two molecules. The inhibition of proteasomal protein degradation by epoxomicin increased the PfRFUL level, suggesting the degradation of PfRFUL by the proteasome pathways, whereas the inhibition of E3 ubiquitin ligase activities by JNJ26854165 reduced the PfRFUL level. This study reveals the potential mechanisms of PfRFUL in modifying the expression of drug transporters and their roles in parasite drug responses. PfRFUL could be a potential target for antimalarial drug development.

Ring stage dormancy of Plasmodium falciparum tolerant to artemisinin and its analogues - A genetically regulated "Sleeping Beauty"

The appearance in 2008 in western Cambodia of Plasmodium falciparum tolerant to artemisinin, defined by longer parasite clearance time following drug administration and in vitro by a slightly higher survival rate of the ring stage after a 3-h treatment with 700 nM artemisinin (or analogues, collectively termed ART), has raised concerns of the possible loss of this frontline antimalarial [used in the form of an artemisinin combination therapy (ACT)], with its low IC₅₀ value against the ring stage and pleiotropic pro-drug/poison property. The key genetic marker of ART tolerance phenotype is a number of non-synonymous mutations in Pfkelch13 propeller domain. This results in defective assembly at the ring stage of a cytostome structure located at cytoplasmic side of the parasite membrane required for invagination of a double-membrane endosome carrying host cytosol haemoglobin to the digestive vacuole. The consequential deprivation of amino acids initiates ring stage parasites bearing the causal mutations in PfK13 (or other key cytostome components) entry into a dormant state ("Sleeping Beauty"), which, after a duration longer than that the short-lived ART, "Sleeping Beauty" ring parasite resumes its normal, but accelerated, development to maintain the 48-h intra-erythrocytic life-cycle. We posit that when ART-tolerant P. falciparum has acquired under ART stress the causative PfK13 mutation (not obligatory if mutations occur in other critical cytostome components), together with other necessary mutations to adjust to the new normalcy and to provide survival competitiveness, ART-tolerant parasite has now evolved into a genetically programmed "Sleeping Beauty". The onus of preventing the spread of ART-tolerant P. falciparum lies with the efficacy of ACT partner drug, hence the recommendation of a triple ACT (TACT). Nevertheless, attention should also be focussed on understanding the mechanisms of dormancy, such as induction, maintenance and recovery, to enable discovery and development of novel antimalarials targeting this unique parasite stage.

Naturally Acquired Kelch13 Mutations in Plasmodium falciparum Strains Modulate In Vitro Ring-Stage Artemisinin-Based Drug Tolerance and Parasite Survival in Response to Hyperoxia

The ring-stage survival assay was utilized to assess the impact of physiological hyperoxic stress on dihydroartemisinin (DHA) tolerance for a panel of Plasmodium falciparum strains with and without Kelch13 mutations. Strains without naturally acquired Kelch13 mutations or the postulated genetic background associated with delayed parasite clearance time demonstrated reduced proliferation under hyperoxic conditions in the subsequent proliferation cycle. Dihydroartemisinin tolerance in three isolates with naturally acquired Kelch13 mutations but not two genetically manipulated laboratory strains was modulated by in vitro hyperoxic stress exposure of early-ring-stage parasites in the cycle before drug exposure. Reduced parasite tolerance to additional derivatives, including artemisinin, artesunate, and OZ277, was observed within the second proliferation cycle. OZ439 and epoxomicin completely prevented parasite survival under both hyperoxia and normoxic in vitro culture conditions, highlighting the unique relationship between DHA tolerance and Kelch13 mutation-associated genetic background. IMPORTANCE Artemisinin-based combination therapy (ACT) for treating malaria is under intense scrutiny following treatment failures in the Greater Mekong subregion of Asia. This is further compounded by the potential for extensive loss of life if treatment failures extend to the African continent. Although Plasmodium falciparum has become resistant to all antimalarial drugs, artemisinin "resistance" does not present in the same way as resistance to other antimalarial drugs. Instead, a partial resistance or tolerance is demonstrated, associated with the parasite's genetic profile and linked to a molecular marker referred to as K13. It is suggested that parasites may have adapted to drug treatment, as well as the presence of underlying population health issues such as hemoglobinopathies, and/or environmental pressures, resulting in parasite tolerance to ACT. Understanding parasite evolution and control of artemisinin tolerance will provide innovative approaches to mitigate the development of artemisinin tolerance and thereby artemisinin-based drug treatment failure and loss of life globally to malaria infections.

On the contribution of the rodent model Plasmodium chabaudi for understanding the genetics of drug resistance in malaria

Malaria is a devastating disease that still claims over half a million lives every year, mostly in sub-Saharan Africa. One of the main barriers to malaria control is the evolution and propagation of drug-resistant mutant parasites. Knowing the genes and respective mutations responsible for drug resistance facilitates the design of drugs with novel modes of action and allows predicting and monitoring drug resistance in natural parasite populations in real-time. The best way to identify these mutations is to experimentally evolve resistance to the drug in question and then comparing the genomes of the drug-resistant mutants to that of the sensitive progenitor parasites. This simple evolutive concept was the starting point for the development of a paradigm over the years, based on the use of the rodent malaria parasite Plasmodium chabaudi to unravel the genetics of drug resistance in malaria. It involves the use of a cloned parasite isolate (P. chabaudi AS) whose genome is well characterized, to artificially select resistance to given drugs through serial passages in mice under slowly increasing drug pressure. The end resulting parasites are cloned and the genetic mutations are then discovered through Linkage Group Selection, a technique conceived by Prof. Richard Carter and his group, and/or Whole Genome Sequencing. The precise role of these mutations can then be interrogated in malaria parasites of humans through allelic replacement experiments and/or genotype-phenotype association studies in natural parasite populations. Using this paradigm, all the mutations underlying resistance to the most important antimalarial drugs were identified, most of which were pioneering and later shown to also play a role in drug resistance in natural infections of human malaria parasites. This supports the use of P. chabaudi a fast-track predictive model to identify candidate genetic markers of resistance to present and future antimalarial drugs and improving our understanding of the biology of resistance.