Combinatorial Genetic Modeling of pfcrt-Mediated Drug Resistance Evolution in Plasmodium falciparum

Additional PfCRT mutations driven by selective pressure for improved fitness can result in the loss of piperaquine resistance and altered Plasmodium falciparum physiology

IMPORTANCE

Our study leverages gene editing techniques in Plasmodium falciparum asexual blood stage parasites to profile novel mutations in mutant PfCRT, an important mediator of piperaquine resistance, which developed in Southeast Asian field isolates or in parasites cultured for long periods of time. We provide evidence that increased parasite fitness of these lines is the primary driver for the emergence of these PfCRT variants. These mutations differentially impact parasite susceptibility to piperaquine and chloroquine, highlighting the multifaceted effects of single point mutations in this transporter. Molecular features of drug resistance and parasite physiology were examined in depth using proteoliposome-based drug uptake studies and peptidomics, respectively. Energy minimization calculations, showing how these novel mutations might impact the PfCRT structure, suggested a small but significant effect on drug interactions. This study reveals the subtle interplay between antimalarial resistance, parasite fitness, PfCRT structure, and intracellular peptide availability in PfCRT-mediated parasite responses to changing drug selective pressures.

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

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

Evolution of genetic markers for drug resistance after the introduction of dihydroartemisinin-piperaquine as first-line anti-malarial treatment for uncomplicated falciparum malaria in Indonesia

BACKGROUND

Dihydroartemisinin-piperaquine has been Indonesia's first-line anti-malarial treatment since 2008. Annual therapeutic efficacy studies (TES) done in the last 12 years showed continued high treatment efficacy in uncomplicated Plasmodium falciparum malaria. Although these studies did not show evidence for artemisinin resistance, a slight increase in Late Treatment Failure was observed over time. It is highlight to explore the evolution of genetic markers for ACT partner drug resistance since adopting DHA-PPQ.

METHODS

Dry blood spots were identified from a mass blood survey of uncomplicated falciparum malaria patients (N = 50) in Sumba from 2010 to 2018. Analysis of genotypic profile (N = 51) and a Therapeutic Efficacy Study (TES) from Papua (N = 142) from 2020 to 2021, 42-day follow-up. PCR correction using msp1, msp2, and glurp was used to distinguish recrudescence and reinfection. Parasite DNA from DBSs was used for genotyping molecular markers for antimalaria drug resistance, including in Pfk13, pfcrt, and pfmdr1, as well as gene copy number variation in pfpm2/3 and pfmdr1.

RESULTS

The study revealed the absence of SNPs associated with ART resistance and several novel SNPs such as L396F, I526V, M579I and N537S (4.25%). In Sumba, the mutant haplotype SDD of pfmdr1 was found in one-third of the isolates, while only 8.9% in Papua. None of the pfcrt mutations linked to piperaquine resistance were observed, but 71% of isolates had pfcrt I356L. Amplification of the pfpm2/3 genes was in Sumba (17.02%) and Papua (13.7%), while pfmdr1 copy number prevalence was low (3.8%) in both areas. For the TES study, ten recurrences of infection were observed on days 28, 35, and 42. Late parasitological failure (LPF) was observed in 10/117 (8.5%) subjects by microscopy. PCR correction revealed that all nine cases were re-infections and one was confirmed as recrudescence.

CONCLUSION

This study revealed that DHA-PPQ is still highly effective against P. falciparum. The genetic architecture of the parasite P. falciparum isolates during 2010-2021 revealed single copy of Pfpm2 and pfmdr1 were highly prevalent. The slight increase in DHA-PPQ LTF alerts researchers to start testing other ACTs as alternatives to DHA-PPQ for baseline data in order to get a chance of achieving malaria elimination wants by 2030.

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

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

An Alternative Autophagy-Related Mechanism of Chloroquine Drug Resistance in the Malaria Parasite

We generated highly chloroquine (CQ)-resistant (ResCQ) Plasmodium yoelii parasites by stepwise exposure to increasing concentrations of CQ and CQ-sensitive parasites (SenCQ) by parallel mock treatments. No mutations in genes that are associated with drug resistance were detected in ResCQ clones. Autophagy-related genes were highly upregulated in SenCQ compared to ResCQ parasites during CQ treatment. This indicates that CQ resistance can be developed in the malaria parasite by the inhibition of autophagy as an alternative drug resistance mechanism.

Mutant PfCRT Can Mediate Piperaquine Resistance in African Plasmodium falciparum With Reduced Fitness and Increased Susceptibility to Other Antimalarials

BACKGROUND

Additional therapeutic strategies could benefit efforts to reverse the recent increase in malaria cases in sub-Saharan Africa, which mostly affects young children. A primary candidate is dihydroartemisinin + piperaquine (DHA + PPQ), which is effective for uncomplicated malaria treatment, seasonal malaria chemoprevention, and intermittent preventive treatment. In Southeast Asia, Plasmodium falciparum parasites acquired PPQ resistance, mediated primarily by mutations in the P falciparum chloroquine resistance transporter PfCRT. The recent emergence in Africa of DHA-resistant parasites creates an imperative to assess whether PPQ resistance could emerge in African parasites with distinct PfCRT isoforms.

METHODS

We edited 2 PfCRT mutations known to mediate high-grade PPQ resistance in Southeast Asia into GB4 parasites from Gabon. Gene-edited clones were profiled in antimalarial concentration-response and fitness assays.

RESULTS

The PfCRT F145I mutation mediated moderate PPQ resistance in GB4 parasites but with a substantial fitness cost. No resistance was observed with the PfCRT G353V mutant. Both edited clones became significantly more susceptible to amodiaquine, chloroquine, and quinine.

CONCLUSIONS

A single PfCRT mutation can mediate PPQ resistance in GB4 parasites, but with a growth defect that may preclude its spread without further genetic adaptations. Our findings support regional use of drug combinations that exert opposing selective pressures on PfCRT.

Barcoding Genetically Distinct Plasmodium falciparum Strains for Comparative Assessment of Fitness and Antimalarial Drug Resistance

The repeated emergence of antimalarial drug resistance in Plasmodium falciparum, including to the current frontline antimalarial artemisinin, is a perennial problem for malaria control. Next-generation sequencing has greatly accelerated the identification of polymorphisms in resistance-associated genes but has also highlighted the need for more sensitive and accurate laboratory tools to profile current and future antimalarials and to quantify the impact of drug resistance acquisition on parasite fitness. The interplay of fitness and drug response is of fundamental importance in understanding why particular genetic backgrounds are better at driving the evolution of drug resistance in natural populations, but the impact of parasite fitness landscapes on the epidemiology of drug resistance has typically been laborious to accurately quantify in the lab, with assays being limited in accuracy and throughput. Here we present a scalable method to profile fitness and drug response of genetically distinct P. falciparum strains with well-described sensitivities to several antimalarials. We leverage CRISPR/Cas9 genome-editing and barcode sequencing to track unique barcodes integrated into a nonessential gene (pfrh3). We validate this approach in multiplex competitive growth assays of three strains with distinct geographical origins. Furthermore, we demonstrate that this method can be a powerful approach for tracking artemisinin response as it can identify an artemisinin resistant strain within a mix of multiple parasite lines, suggesting an approach for scaling the laborious ring-stage survival assay across libraries of barcoded parasite lines. Overall, we present a novel high-throughput method for multiplexed competitive growth assays to evaluate parasite fitness and drug response. IMPORTANCE The complex interplay between antimalarial resistance and parasite fitness has important implications for understanding the development and spread of drug resistance alleles and the impact of genetic background on transmission. One limitation with current methodologies to measure parasite fitness is the ability to scale this beyond simple head-to-head competition experiments between a wildtype control line and test line, with a need for a scalable approach that allows tracking of parasite growth in complex mixtures. In our study, we have used CRISPR editing to insert unique DNA barcodes into a safe-harbor genomic locus to tag multiple parasite strains and use next-generation sequencing to read out strain dynamics. We observe inherent fitness differences between the strains, as well as sensitive modulation of responses to challenge with clinically relevant antimalarials, including artemisinin.

Genome-wide SNP analysis of Plasmodium falciparum shows differentiation at drug-resistance-associated loci among malaria transmission settings in southern Mali

Plasmodium falciparum malaria cases in Africa represent over 90% of the global burden with Mali being amongst the 11 highest burden countries that account for 70% of this annual incidence. The persistence of P. falciparum despite massive global interventions is because of its genetic diversity that drives its ability to adapt to environmental changes, develop resistance to drugs, and evade the host immune system. Knowledge on P. falciparum genetic diversity across populations and intervention landscape is thus critical for the implementation of new strategies to eliminate malaria. This study assessed genetic variation with 12,177 high-quality SNPs from 830 Malian P. falciparum isolates collected between 2007 and 2017 from seven locations. The complexity of infections remained high, varied between sites, and showed a trend toward overall decreasing complexity over the decade. Though there was no significant substructure, allele frequencies varied geographically, partly driven by temporal variance in sampling, particularly for drug resistance and antigen loci. Thirty-two mutations in known drug resistance markers (pfcrt, pfdhps, pfdhfr, pfmdr1, pfmdr2, and pfk13) attained a frequency of at least 2% in the populations. SNPs within and around the major markers of resistance to quinolines (pfmdr1 and pfcrt) and antifolates (pfdhfr and pfdhps) varied temporally and geographically, with strong linkage disequilibrium and signatures of directional selection in the genome. These geo-temporal populations also differentiated at alleles in immune-related loci, including, protein E140, pfsurfin8, pfclag8, and pfceltos, as well as pftrap, which showed signatures of haplotype differentiation between populations. Several regions across the genomes, including five known drug resistance loci, showed signatures of differential positive selection. These results suggest that drugs and immune pressure are dominant selective forces against P. falciparum in Mali, but their effect on the parasite genome varies temporally and spatially. Interventions interacting with these genomic variants need to be routinely evaluated as malaria elimination strategies are implemented.

Efficient generation of mNeonGreen Plasmodium falciparum reporter lines enables quantitative fitness analysis

CRISPR editing has enabled the rapid creation of fluorescent Plasmodium transgenic lines, facilitating a deeper understanding of parasite biology. The impact of genetic perturbations such as gene disruption or the introduction of drug resistance alleles on parasite fitness is typically quantified in competitive growth assays between the query line and a wild type reference. Although fluorescent reporter lines offer a facile and frequently used method to measure relative growth, this approach is limited by the strain background of the existing reporter, which may not match the growth characteristics of the query strains, particularly if these are slower-growing field isolates. Here, we demonstrate an efficient CRISPR-based approach to generate fluorescently labelled parasite lines using mNeonGreen derived from the LanYFP protein in Branchiostoma lanceolatum, which is one of the brightest monomeric green fluorescent proteins identified. Using a positive-selection approach by insertion of an in-frame blasticidin S deaminase marker, we generated a Dd2 reporter line expressing mNeonGreen under the control of the pfpare (P. falciparum Prodrug Activation and Resistance Esterase) locus. We selected the pfpare locus as an integration site because it is highly conserved across P. falciparum strains, expressed throughout the intraerythrocytic cycle, not essential, and offers the potential for negative selection to further enrich for integrants. The mNeonGreen@pare line demonstrates strong fluorescence with a negligible fitness defect. In addition, the construct developed can serve as a tool to fluorescently tag other P. falciparum strains for in vitro experimentation.

Evidence for the early emergence of piperaquine-resistant Plasmodium falciparum malaria and modeling strategies to mitigate resistance

Multidrug-resistant Plasmodium falciparum parasites have emerged in Cambodia and neighboring countries in Southeast Asia, compromising the efficacy of first-line antimalarial combinations. Dihydroartemisinin + piperaquine (PPQ) treatment failure rates have risen to as high as 50% in some areas in this region. For PPQ, resistance is driven primarily by a series of mutant alleles of the P. falciparum chloroquine resistance transporter (PfCRT). PPQ resistance was reported in China three decades earlier, but the molecular driver remained unknown. Herein, we identify a PPQ-resistant pfcrt allele (China C) from Yunnan Province, China, whose genotypic lineage is distinct from the PPQ-resistant pfcrt alleles currently observed in Cambodia. Combining gene editing and competitive growth assays, we report that PfCRT China C confers moderate PPQ resistance while re-sensitizing parasites to chloroquine (CQ) and incurring a fitness cost that manifests as a reduced rate of parasite growth. PPQ transport assays using purified PfCRT isoforms, combined with molecular dynamics simulations, highlight differences in drug transport kinetics and in this transporter’s central cavity conformation between China C and the current Southeast Asian PPQ-resistant isoforms. We also report a novel computational model that incorporates empirically determined fitness landscapes at varying drug concentrations, combined with antimalarial susceptibility profiles, mutation rates, and drug pharmacokinetics. Our simulations with PPQ-resistant or -sensitive parasite lines predict that a three-day regimen of PPQ combined with CQ can effectively clear infections and prevent the evolution of PfCRT variants. This work suggests that including CQ in combination therapies could be effective in suppressing the evolution of PfCRT-mediated multidrug resistance in regions where PPQ has lost efficacy.

The recent emergence of Plasmodium falciparum parasite resistance to the antimalarial drug piperaquine (PPQ) has contributed to frequent treatment failures across Southeast Asia, originating in Cambodia. Here, we show that earlier reports of PPQ resistance in Yunnan Province, China could be explained by the unique China C variant of the P. falciparum chloroquine resistance transporter PfCRT. Gene-edited parasites show a loss of fitness and parasite resensitization to the chemically related former first-line antimalarial chloroquine, while acquiring PPQ resistance via drug efflux. Molecular features of drug resistance were examined using biochemical assays to measure mutant PfCRT-mediated drug transport and molecular dynamics simulations with the recently solved PfCRT structure to assess changes in the central drug-binding cavity. We also describe a new computational model that incorporates parasite mutation rates, fitness costs, antimalarial susceptibilities, and drug pharmacological profiles to predict how infections with parasite strains expressing distinct PfCRT variants can evolve and be selected in response to different drug pressures and regimens. Simulations predict that a three-day regimen of PPQ plus chloroquine would be fully effective at preventing recrudescence of drug-resistant infections.

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.

The power and promise of genetic mapping from Plasmodium falciparum crosses utilizing human liver-chimeric mice

Genetic crosses are most powerful for linkage analysis when progeny numbers are high, parental alleles segregate evenly and numbers of inbred progeny are minimized. We previously developed a novel genetic crossing platform for the human malaria parasite Plasmodium falciparum, an obligately sexual, hermaphroditic protozoan, using mice carrying human hepatocytes (the human liver-chimeric FRG NOD huHep mouse) as the vertebrate host. We report on two genetic crosses-(1) an allopatric cross between a laboratory-adapted parasite (NF54) of African origin and a recently patient-derived Asian parasite, and (2) a sympatric cross between two recently patient-derived Asian parasites. We generated 144 unique recombinant clones from the two crosses, doubling the number of unique recombinant progeny generated in the previous 30 years. The allopatric African/Asian cross has minimal levels of inbreeding and extreme segregation distortion, while in the sympatric Asian cross, inbred progeny predominate and parental alleles segregate evenly. Using simulations, we demonstrate that these progeny provide the power to map small-effect mutations and epistatic interactions. The segregation distortion in the allopatric cross slightly erodes power to detect linkage in several genome regions. We greatly increase the power and the precision to map biomedically important traits with these new large progeny panels.

Structural Insights into Transporter-Mediated Drug Resistance in Infectious Diseases

Infectious diseases present a major threat to public health globally. Pathogens can acquire resistance to anti-infectious agents via several means including transporter-mediated efflux. Typically, multidrug transporters feature spacious, dynamic, and chemically malleable binding sites to aid in the recognition and transport of chemically diverse substrates across cell membranes. Here, we discuss recent structural investigations of multidrug transporters involved in resistance to infectious diseases that belong to the ATP-binding cassette (ABC) superfamily, the major facilitator superfamily (MFS), the drug/metabolite transporter (DMT) superfamily, the multidrug and toxic compound extrusion (MATE) family, the small multidrug resistance (SMR) family, and the resistance-nodulation-division (RND) superfamily. These structural insights provide invaluable information for understanding and combatting multidrug resistance.

Historical trends and new surveillance of Plasmodium falciparum drug resistance markers in Angola

BACKGROUND

Plasmodium falciparum resistance to chloroquine (CQ) and sulfadoxine-pyrimethamine (SP) has historically posed a major threat to malaria control throughout the world. The country of Angola officially replaced CQ with artemisinin-based combination therapy (ACT) as a first-line treatment in 2006, but malaria cases and deaths have recently been rising. Many classic resistance mutations are relevant for the efficacy of currently available drugs, making it important to continue monitoring their frequency in Angola.

METHODS

Plasmodium falciparum DNA was sampled from the blood of 50 hospital patients in Cabinda, Angola from October-December of 2018. Each infection was genotyped for 13 alleles in the genes crt, mdr1, dhps, dhfr, and kelch13, which are collectively involved in resistance to six common anti-malarials. To compare frequency patterns over time, P. falciparum genotype data were also collated from studies published from across Angola in the last two decades.

RESULTS

The two most important alleles for CQ resistance, crt 76T and mdr1 86Y, were found at respective frequencies of 71.4% and 6.5%. Historical data suggest that mdr1 N86 has been steadily replacing 86Y throughout Angola in the last decade, while the frequency of crt 76T has been more variable across studies. Over a third of new samples from Cabinda were 'quintuple mutants' for SP resistance in dhfr/dhps, with a sixth mutation at dhps A581G present at 9.6% frequency. The markers dhfr 51I, dhfr 108N, and dhps 437G have been nearly fixed in Angola since the early 2000s, whereas dhfr 59R may have risen to high frequency more recently. Finally, no non-synonymous polymorphisms were detected in kelch13, which is involved in artemisinin resistance in Southeast Asia.

CONCLUSIONS

Genetic markers of P. falciparum resistance to CQ are likely declining in frequency in Angola, consistent with the official discontinuation of CQ in 2006. The high frequency of multiple genetic markers of SP resistance is consistent with the continued public and private use of SP. In the future, more complete haplotype data from mdr1, dhfr, and dhps will be critical for understanding the changing efficacy of multiple anti-malarial drugs. These data can be used to support effective drug policy decisions in Angola.

In vitro growth competition experiments that suggest consequences of the substandard artemisinin epidemic that may be accelerating drug resistance in P. falciparum malaria

Over the past decade, artemisinin (ART)-combination therapies (ACTs) have shown declining efficacy within Southeast Asia (SEA). These resistance-like phenomena manifest as a delayed clearance phenotype (DCP) in some patients treated with ACTs. ACTs are currently the recommended treatment for P. falciparum infections by the World Health Organization (WHO), and they are our last line of defense to effectively treat all strains of malaria. Acceleration of antimicrobial resistance (AMR) is often theorized to be exacerbated by the use of subtherapeutic dosages of drugs ("substandard" drug), which for ACTs has been well documented over the last decade. Troublingly, in 2017, the WHO estimated that nearly 1 in 10 medical products tested in low- and middle-income countries failed to meet quality standards. We have developed a tissue culture-based approach for testing possible connections between substandard treatment and the spread of ACT resistant blood stage forms of P. falciparum. Via sequencing of pfk13, a molecular marker that is predictive for ART resistance (ARTR), we monitor competition of sensitive vs resistant strains over time and under various conditions and define conditions that favor emergence of ARTR parasites. Our findings help to define the conditions under which substandard drug treatments might favor the proliferation of mutant PfK13-mediated drug resistant strains over drug sensitive.

Novel Antimalarial Tetrazoles and Amides Active against the Hemoglobin Degradation Pathway in Plasmodium falciparum

Malaria control programs continue to be threatened by drug resistance. To identify new antimalarials, we conducted a phenotypic screen and identified a novel tetrazole-based series that shows fast-kill kinetics and a relatively low propensity to develop high-level resistance. Preliminary structure-activity relationships were established including identification of a subseries of related amides with antiplasmodial activity. Assaying parasites with resistance to antimalarials led us to test whether the series had a similar mechanism of action to chloroquine (CQ). Treatment of synchronized Plasmodium falciparum parasites with active analogues revealed a pattern of intracellular inhibition of hemozoin (Hz) formation reminiscent of CQ's action. Drug selections yielded only modest resistance that was associated with amplification of the multidrug resistance gene 1 (pfmdr1). Thus, we have identified a novel chemical series that targets the historically druggable heme polymerization pathway and that can form the basis of future optimization efforts to develop a new malaria treatment.

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.

Inferring a complete genotype-phenotype map from a small number of measured phenotypes

Understanding evolution requires detailed knowledge of genotype-phenotype maps; however, it can be a herculean task to measure every phenotype in a combinatorial map. We have developed a computational strategy to predict the missing phenotypes from an incomplete, combinatorial genotype-phenotype map. As a test case, we used an incomplete genotype-phenotype dataset previously generated for the malaria parasite’s ‘chloroquine resistance transporter’ (PfCRT). Wild-type PfCRT (PfCRT^3D7) lacks significant chloroquine (CQ) transport activity, but the introduction of the eight mutations present in the ‘Dd2’ isoform of PfCRT (PfCRT^Dd2) enables the protein to transport CQ away from its site of antimalarial action. This gain of a transport function imparts CQ resistance to the parasite. A combinatorial map between PfCRT^3D7 and PfCRT^Dd2 consists of 256 genotypes, of which only 52 have had their CQ transport activities measured through expression in the Xenopus laevis oocyte. We trained a statistical model with these 52 measurements to infer the CQ transport activity for the remaining 204 combinatorial genotypes between PfCRT^3D7 and PfCRT^Dd2. Our best-performing model incorporated a binary classifier, a nonlinear scale, and additive effects for each mutation. The addition of specific pairwise- and high-order-epistatic coefficients decreased the predictive power of the model. We evaluated our predictions by experimentally measuring the CQ transport activities of 24 additional PfCRT genotypes. The R² value between our predicted and newly-measured phenotypes was 0.90. We then used the model to probe the accessibility of evolutionary trajectories through the map. Approximately 1% of the possible trajectories between PfCRT^3D7 and PfCRT^Dd2 are accessible; however, none of the trajectories entailed eight successive increases in CQ transport activity. These results demonstrate that phenotypes can be inferred with known uncertainty from a partial genotype-phenotype dataset. We also validated our approach against a collection of previously published genotype-phenotype maps. The model therefore appears general and should be applicable to a large number of genotype-phenotype maps.

Biological macromolecules are built from chains of building blocks. The function of a macromolecule depends on the specific chemical properties of the building blocks that make it up. Macromolecules evolve through mutations that swap one building block for another. Understanding how biomolecules work and evolve therefore requires knowledge of the effects of mutations. The effects of mutations can be measured experimentally; however, because there are a vast number of possible combinations of mutations, it is often difficult to make enough measurements to understand biomolecular function and evolution. In this paper, we describe a simple method to predict the effects of mutations on biomolecules from a small number of measurements. This method works by appropriately averaging the effects of mutations seen in different contexts. We test the method by predicting the effects of mutations on a PfCRT—a macromolecule from the malarial parasite that confers drug resistance. We find that our method is fast and effective. Using a small number of measurements, we were able to gain insight into the evolutionary steps by which this macromolecule conferred drug resistance. To make this method accessible to other researchers, we have released it as an open-source software package: https://gpseer.readthedocs.io.

Role of a Concentration Gradient in Malaria Drug Resistance Evolution: A Combined within- and between-Hosts Modelling Approach

Resistance to antimalarial drugs is currently a growing public health problem, resulting in more cases with treatment failure. Although previous studies suggested that a concentration gradient facilitates the antibiotic resistance evolution in bacteria, no attempt has been made to investigate the roles of a concentration gradient in malaria drug resistance. Unlike the person-to-person mode of transmission of bacteria, the malaria parasites need to switch back and forth between the human and mosquito hosts to complete the life cycle and to spread the resistant alleles. Here we developed a stochastic combined within- and between-hosts evolutionary dynamics model specific to malaria parasites in order to investigate the influence of an antimalarial concentration gradient on the evolutionary dynamics of malaria drug resistance. Every stage of malaria development in both human and mosquito hosts are individually modelled using the tau-leaping algorithm. We found that the concentration gradient can accelerate antimalarial resistance evolution. The gain in resistance evolution was improved by the increase in the parasite mutation rate and the mosquito biting rate. In addition, even though the rate of resistance evolution is not sensitive to the changes in parasite reduction ratios (PRRs) of antimalarial drugs, the probability of finding the antimalarial drug resistant parasites decreases when the PRR increases.

Elucidating Mechanisms of Drug-Resistant Plasmodium falciparum

Intensified treatment and control efforts since the early 2000s have dramatically reduced the burden of Plasmodium falciparum malaria. However, drug resistance threatens to derail this progress. In this review, we present four antimalarial resistance case studies that differ in timeline, technical approaches, mechanisms of action, and categories of resistance: chloroquine, sulfadoxine-pyrimethamine, artemisinin, and piperaquine. Lessons learned from prior losses of treatment efficacy, drug combinations, and control strategies will help advance mechanistic research into how P. falciparum parasites acquire resistance to current first-line artemisinin-based combination therapies. Understanding resistance in the clinic and laboratory is essential to prolong the effectiveness of current antimalarial drugs and to optimize the pipeline of future medicines.

Pairwise growth competitions identify relative fitness relationships among artemisinin resistant Plasmodium falciparum field isolates

Background

Competitive outcomes between co-infecting malaria parasite lines can reveal fitness disparities in blood stage growth. Blood stage fitness costs often accompany the evolution of drug resistance, with the expectation that relatively fitter parasites will be more likely to spread in populations. With the recent emergence of artemisinin resistance, it is important to understand the relative competitive fitness of the metabolically active asexual blood stage parasites. Genetically distinct drug resistant parasite clones with independently evolved sets of mutations are likely to vary in asexual proliferation rate, contributing to their chance of transmission to the mosquito vector.

Methods

An optimized in vitro 96-well plate-based protocol was used to quantitatively measure-head-to-head competitive fitness during blood stage development between seven genetically distinct field isolates from a hotspot of emerging artemisinin resistance and the laboratory strain, NF54. These field isolates were isolated from patients in Southeast Asia carrying different alleles of kelch13 and included both artemisinin-sensitive and artemisinin-resistant isolates. Fluorescent labeled microsatellite markers were used to track the relative densities of each parasite throughout the co-growth period of 14–60 days. All-on-all competitions were conducted for the panel of eight parasite lines (28 pairwise competitions) to determine their quantitative competitive fitness relationships.

Results

Twenty-eight pairwise competitive growth outcomes allowed for an unambiguous ranking among a set of seven genetically distinct parasite lines isolated from patients in Southeast Asia displaying a range of both kelch13 alleles and clinical clearance times and a laboratory strain, NF54. This comprehensive series of assays established the growth relationships among the eight parasite lines. Interestingly, a clinically artemisinin resistant parasite line that carries the wild-type form of kelch13 outcompeted all other parasites in this study. Furthermore, a kelch13 mutant line (E252Q) was competitively more fit without drug than lines with other resistance-associated kelch13 alleles, including the C580Y allele that has expanded to high frequencies under drug pressure in Southeast Asian resistant populations.

Conclusions

This optimized competitive growth assay can be employed for assessment of relative growth as an index of fitness during the asexual blood stage growth between natural lines carrying different genetic variants associated with artemisinin resistance. Improved understanding of the fitness costs of different parasites proliferating in human blood and the role different resistance mutations play in the context of specific genetic backgrounds will contribute to an understanding of the potential for specific mutations to spread in populations, with the potential to inform targeted strategies for malaria therapy.

Phosphomimetic substitution at Ser-33 of the chloroquine resistance transporter PfCRT reconstitutes drug responses in Plasmodium falciparum

The chloroquine resistance transporter PfCRT of the human malaria parasite Plasmodium falciparum confers resistance to the former first-line antimalarial drug chloroquine, and it modulates the responsiveness to a wide range of quinoline and quinoline-like compounds. PfCRT is post-translationally modified by phosphorylation, palmitoylation, and, possibly, ubiquitination. However, the impact of these post-translational modifications on P. falciparum biology and, in particular, the drug resistance-conferring activity of PfCRT has remained elusive. Here, we confirm phosphorylation at Ser-33 and Ser-411 of PfCRT of the chloroquine-resistant P. falciparum strain Dd2 and show that kinase inhibitors can sensitize drug responsiveness. Using CRISPR/Cas9 genome editing to generate genetically engineered PfCRT variants in the parasite, we further show that substituting Ser-33 with alanine reduced chloroquine and quinine resistance by ∼50% compared with the parental P. falciparum strain Dd2, whereas the phosphomimetic amino acid aspartic acid could fully and glutamic acid could partially reconstitute the level of chloroquine/quinine resistance. Transport studies conducted in the parasite and in PfCRT-expressing Xenopus laevis oocytes linked phosphomimetic substitution at Ser-33 to increased transport velocity. Our data are consistent with phosphorylation of Ser-33 relieving an autoinhibitory intramolecular interaction within PfCRT, leading to a stimulated drug transport activity. Our findings shed additional light on the function of PfCRT and suggest that chloroquine could be reevaluated as an antimalarial drug by targeting the kinase in P. falciparum that phosphorylates Ser-33 of PfCRT.

Genetic diversity and drug resistance surveillance of Plasmodium falciparum for malaria elimination: is there an ideal tool for resource-limited sub-Saharan Africa?

The intensification of malaria control interventions has resulted in its global decline, but it remains a significant public health burden especially in sub-Saharan Africa (sSA). Knowledge on the parasite diversity, its transmission dynamics, mechanisms of adaptation to environmental and interventional pressures could help refine or develop new control and elimination strategies. Critical to this is the accurate assessment of the parasite’s genetic diversity and monitoring of genetic markers of anti-malarial resistance across all susceptible populations. Such wide molecular surveillance will require selected tools and approaches from a variety of ever evolving advancements in technology and the changing epidemiology of malaria. The choice of an effective approach for specific endemic settings remains challenging, particularly for countries in sSA with limited access to advanced technologies. This article examines the current strategies and tools for Plasmodium falciparum genetic diversity typing and resistance monitoring and proposes how the different tools could be employed in resource-poor settings. Advanced approaches enabling targeted deep sequencing is valued as a sensitive method for assessing drug resistance and parasite diversity but remains out of the reach of most laboratories in sSA due to the high cost of development and maintenance. It is, however, feasible to equip a limited number of laboratories as Centres of Excellence in Africa (CEA), which will receive and process samples from a network of peripheral laboratories in the continent. Cheaper, sensitive and portable real-time PCR methods can be used in peripheral laboratories to pre-screen and select samples for targeted deep sequence or genome wide analyses at these CEAs.

Global Spread of Mutant PfCRT and Its Pleiotropic Impact on Plasmodium falciparum Multidrug Resistance and Fitness

The global spread of Plasmodium falciparum chloroquine resistance transporter (PfCRT) variant haplotypes earlier caused the widespread loss of chloroquine (CQ) efficacy. In Asia, novel PfCRT mutations that emerged on the Dd2 allelic background have recently been implicated in high-level resistance to piperaquine, and N326S and I356T have been associated with genetic backgrounds in which resistance emerged to artemisinin derivatives. By analyzing large-scale genome sequencing data, we report that the predominant Asian CQ-resistant Dd2 haplotype is undetectable in Africa. Instead, the GB4 and previously unexplored Cam783 haplotypes predominate, along with wild-type, drug-sensitive PfCRT that has reemerged as the major haplotype. To interrogate how these alleles impact drug susceptibility, we generated pfcrt-modified isogenic parasite lines spanning the mutational interval between GB4 and Dd2, which includes Cam783 and involves amino acid substitutions at residues 326 and 356. Relative to Dd2, the GB4 and Cam783 alleles were observed to mediate lower degrees of resistance to CQ and the first-line drug amodiaquine, while resulting in higher growth rates. These findings suggest that differences in growth rates, a surrogate of parasite fitness, influence selection in the context of African infections that are frequently characterized by high transmission rates, mixed infections, increased immunity, and less recourse to treatment. We also observe that the Asian Dd2 allele affords partial protection against piperaquine yet does not directly impact artemisinin efficacy. Our results can help inform the regional recommendations of antimalarials, whose activity is influenced by and, in certain cases, enhanced against select PfCRT variant haplotypes.IMPORTANCE Our study defines the allelic distribution of pfcrt, an important mediator of multidrug resistance in Plasmodium falciparum, in Africa and Asia. We leveraged whole-genome sequence analysis and gene editing to demonstrate how current drug combinations can select different allelic variants of this gene and shape region-specific parasite population structures. We document the ability of PfCRT mutations to modulate parasite susceptibility to current antimalarials in dissimilar, pfcrt allele-specific ways. This study underscores the importance of actively monitoring pfcrt genotypes to identify emerging patterns of multidrug resistance and help guide region-specific treatment options.

Plasmodium falciparum In Vitro Drug Resistance Selections and Gene Editing

Malaria continues to be a global health burden, threatening over 40% of the world's population. Drug resistance in Plasmodium falciparum, the etiological agent of the majority of human malaria cases, is compromising elimination efforts. New approaches to treating drug-resistant malaria benefit from defining resistance liabilities of known antimalarial agents and compounds in development and defining genetic changes that mediate loss of parasite susceptibility. Here, we present protocols for in vitro selection of drug-resistant parasites and for site-directed gene editing of candidate resistance mediators to test for causality.

Stable Integration and Comparison of hGrx1-roGFP2 and sfroGFP2 Redox Probes in the Malaria Parasite Plasmodium falciparum

Studying redox metabolism in malaria parasites is of great interest for understanding parasite biology, parasite-host interactions, and mechanisms of drug action. Genetically encoded fluorescent redox sensors have recently been described as powerful tools for determining the glutathione-dependent redox potential in living parasites. In the present study, we genomically integrated and expressed the ratiometric redox sensors hGrx1-roGFP2 (human glutaredoxin 1 fused to reduction-oxidation sensitive green fluorescent protein) and sfroGFP2 (superfolder roGFP2) in the cytosol of NF54- attB blood-stage Plasmodium falciparum parasites. Both sensors were evaluated in vitro and in cell culture with regard to their fluorescence properties and reactivity. As genomic integration allows for the stable expression of redox sensors in parasites, we systematically compared single live-cell imaging with plate reader detection. For these comparisons, short-term effects of redox-active compounds were analyzed along with mid- and long-term effects of selected antimalarial agents. Of note, the single components of the redox probes themselves did not influence the redox balance of the parasites. Our analyses revealed comparable results for both the hGrx1-roGFP2 and sfroGFP2 probes, with sfroGFP2 exhibiting a more pronounced fluorescence intensity in cellulo. Accordingly, the sfroGFP2 probe was employed to monitor the fluorescence signals throughout the parasites' asexual life cycle. Through the use of stable genomic integration, we demonstrate a means of overcoming the limitations of transient transfection, allowing more detailed in-cell studies as well as high-throughput analyses using plate reader-based approaches.

Evidence for Regulation of Hemoglobin Metabolism and Intracellular Ionic Flux by the Plasmodium falciparum Chloroquine Resistance Transporter

Plasmodium falciparum multidrug resistance constitutes a major obstacle to the global malaria elimination campaign. Specific mutations in the Plasmodium falciparum chloroquine resistance transporter (PfCRT) mediate resistance to the 4-aminoquinoline drug chloroquine and impact parasite susceptibility to several partner agents used in current artemisinin-based combination therapies, including amodiaquine. By examining gene-edited parasites, we report that the ability of the wide-spread Dd2 PfCRT isoform to mediate chloroquine and amodiaquine resistance is substantially reduced by the addition of the PfCRT L272F mutation, which arose under blasticidin selection. We also provide evidence that L272F confers a significant fitness cost to asexual blood stage parasites. Studies with amino acid-restricted media identify this mutant as a methionine auxotroph. Metabolomic analysis also reveals an accumulation of short, hemoglobin-derived peptides in the Dd2 + L272F and Dd2 isoforms, compared with parasites expressing wild-type PfCRT. Physiologic studies with the ionophores monensin and nigericin support an impact of PfCRT isoforms on Ca²⁺ release, with substantially reduced Ca²⁺ levels observed in Dd2 + L272F parasites. Our data reveal a central role for PfCRT in regulating hemoglobin catabolism, amino acid availability, and ionic balance in P. falciparum, in addition to its role in determining parasite susceptibility to heme-binding 4-aminoquinoline drugs.

Fitness Loss under Amino Acid Starvation in Artemisinin-Resistant Plasmodium falciparum Isolates from Cambodia

Artemisinin is the most rapidly effective drug for Plasmodium falciparum malaria treatment currently in clinical use. Emerging artemisinin-resistant parasites pose a great global health risk. At present, the level of artemisinin resistance is still relatively low with evidence pointing towards a trade-off between artemisinin resistance and fitness loss. Here we show that artemisinin-resistant P. falciparum isolates from Cambodia manifested fitness loss, showing fewer progenies during the intra-erythrocytic developmental cycle. The loss in fitness was exacerbated under the condition of low exogenous amino acid supply. The resistant parasites failed to undergo maturation, whereas their drug-sensitive counterparts were able to complete the erythrocytic cycle under conditions of amino acid deprivation. The artemisinin-resistant phenotype was not stable, and loss of the phenotype was associated with changes in the expression of a putative target, Exp1, a membrane glutathione transferase. Analysis of SNPs in haemoglobin processing genes revealed associations with parasite clearance times, suggesting changes in haemoglobin catabolism may contribute to artemisinin resistance. These findings on fitness and protein homeostasis could provide clues on how to contain emerging artemisinin-resistant parasites.

Emerging Southeast Asian PfCRT mutations confer Plasmodium falciparum resistance to the first-line antimalarial piperaquine

The widely used antimalarial combination therapy dihydroartemisinin + piperaquine (DHA + PPQ) has failed in Cambodia. Here, we perform a genomic analysis that reveals a rapid increase in the prevalence of novel mutations in the Plasmodium falciparum chloroquine resistance transporter PfCRT following DHA + PPQ implementation. These mutations occur in parasites harboring the K13 C580Y artemisinin resistance marker. By introducing PfCRT mutations into sensitive Dd2 parasites or removing them from resistant Cambodian isolates, we show that the H97Y, F145I, M343L, or G353V mutations each confer resistance to PPQ, albeit with fitness costs for all but M343L. These mutations sensitize Dd2 parasites to chloroquine, amodiaquine, and quinine. In Dd2 parasites, multicopy plasmepsin 2, a candidate molecular marker, is not necessary for PPQ resistance. Distended digestive vacuoles were observed in pfcrt-edited Dd2 parasites but not in Cambodian isolates. Our findings provide compelling evidence that emerging mutations in PfCRT can serve as a molecular marker and mediator of PPQ resistance.

Increasing resistance of Plasmodium falciparum strains to piperaquine (PPQ) in Southeast Asia is of concern and resistance mechanisms are incompletely understood. Here, Ross et al. show that mutations in the P. falciparum chloroquine resistance transporter are rapidly increasing in prevalence in Cambodia and confer resistance to PPQ.

Chloroquine exposure triggers distinct cellular responses in sensitive versus resistant Plasmodium falciparum parasites

Chloroquine (CQ) treatment failure in Plasmodium falciparum parasites has been documented for decades, but the pharmacological explanation of this phenotype is not fully understood. Current concepts attribute CQ resistance to reduced accumulation of the drug at a given external CQ concentration ([CQ]_ex) in resistant compared to sensitive parasites. The implication of this explanation is that the mechanisms of CQ-induced toxicity in resistant and sensitive strains are similar once lethal internal concentrations have been reached. To test this hypothesis, we investigated the mechanism of CQ-induced toxicity in CQ-sensitive (CQS) versus CQ-resistant (CQR) parasites by analyzing the time-course of cellular responses in these strains after exposure to varying [CQ]_ex as determined in 72 h toxicity assays. Parasite killing was delayed in CQR parasites for up to 10 h compared to CQS parasites when exposed to equipotent [CQ]_ex. In striking contrast, brief exposure (1 h) to lethal [CQ]_ex in CQS but not CQR parasites caused the appearance of hitherto undescribed hemozoin (Hz)-containing compartments in the parasite cytosol. Hz-containing compartments were very rarely observed in CQR parasites even after CQ exposures sufficient to cause irreversible cell death. These findings challenge current concepts that CQ killing of malaria parasites is solely concentration-dependent, and instead suggest that CQS and CQR strains fundamentally differ in the consequences of CQ exposure.

Variation in intronic microsatellites and exon 2 of the Plasmodium falciparum chloroquine resistance transporter gene during modification of artemisinin combination therapy in Thailand

The amino acid substitution at residue 76 of the food vacuolar transmembrane protein encoded by the chloroquine resistance transporter gene of Plasmodium falciparum (Pfcrt) is an important, albeit imperfect, determinant of chloroquine susceptibility status of the parasite. Other mutations in Pfcrt can modulate susceptibility of P. falciparum to other antimalarials capable of interfering with heme detoxification process, and may exert compensatory effect on parasite growth rate. To address whether nationwide implementation of artemisinin combination therapy (ACT) in Thailand could affect sequence variation in exon 2 and introns of Pfcrt, we analyzed 136 P. falciparum isolates collected during 1997 and 2016 from endemic areas bordering Myanmar, Cambodia and Malaysia. Results revealed 6 haplotypes in exon 2 of Pfcrt with 2 novel substitutions at c.243A > G (p.R81) and c.251A > T (p.N84I). Positive selection was observed at amino acid residues 75, 76 and 97. Four, 3, and 2 alleles of microsatellite (AT/TA) repeats occurred in introns 1, 2 and 4, respectively, resulting in 7 different 3-locus haplotypes. The number of haplotypes and haplotype diversity of exon 2, and introns 1, 2 and 4 were significantly greater among isolates collected during 2009 and 2016 than those collected during 1997 and 2008 when 3-day ACT and 2-day ACT regimens were implemented nationwide, respectively (p < 0.05). By contrast, the number of haplotypes and haplotype diversity of the merozoite surface proteins 1 and 2 of these parasite populations did not differ significantly between these periods. Therefore, the Pfcrt locus of P. falciparum in Thailand continues to evolve and could have been affected by selective pressure from modification of ACT regimen.

Analysis of Plasmodium vivax Chloroquine Resistance Transporter Mutant Isoforms

Chloroquine (CQ) resistance (CQR) in Plasmodium falciparum malaria is widespread and has limited the use of CQ in many regions of the globe. Malaria caused by the related human parasite P. vivax is as widespread as is P. falciparum malaria and has been treated with CQ as extensively as has P. falciparum, suggesting that P. vivax parasites have been selected with CQ as profoundly as have P. falciparum parasites. Indeed, a growing number of clinical reports have presented data suggesting increased P. vivax CQR. Cytostatic (growth inhibitory) CQR for P. falciparum is caused by Plasmodium falciparum chloroquine resistance transporter (PfCRT) mutations, and it has been proposed that mutations in the PvCRT orthologue may simliarly cause P. vivax CQR via increasing CQ transport from the P. vivax digestive vacuole. Here we report the first quantitative analysis of drug transport mediated by all known mutant isoforms of Plasmodium vivax chloroquine resistance transporter (PvCRT) in order to test the protein's potential link to growing P. vivax CQR phenomena. Small, but statistically significant, differences in the transport of CQ and other quinoline antimalarial drugs were found for multiple PvCRT isoforms, relative to wild type PvCRT, suggesting that mutations in PvCRT can contribute to P. vivax CQR and other examples of quinoline antimalarial drug resistance.

Antimalarial drug resistance: linking Plasmodium falciparum parasite biology to the clinic

The global adoption of artemisinin-based combination therapies (ACTs) in the early 2000s heralded a new era in effectively treating drug-resistant Plasmodium falciparum malaria. However, several Southeast Asian countries have now reported the emergence of parasites that have decreased susceptibility to artemisinin (ART) derivatives and ACT partner drugs, resulting in increasing rates of treatment failures. Here we review recent advances in understanding how antimalarials act and how resistance develops, and discuss new strategies for effectively combatting resistance, optimizing treatment and advancing the global campaign to eliminate malaria.

Plasmid-free CRISPR/Cas9 genome editing in Plasmodium falciparum confirms mutations conferring resistance to the dihydroisoquinolone clinical candidate SJ733

Genetic manipulation of the deadly malaria parasite Plasmodium falciparum remains challenging, but the rise of CRISPR/Cas9-based genome editing tools is increasing the feasibility of altering this parasite’s genome in order to study its biology. Of particular interest is the investigation of drug targets and drug resistance mechanisms, which have major implications for fighting malaria. We present a new method for introducing drug resistance mutations in P. falciparum without the use of plasmids or the need for cloning homologous recombination templates. We demonstrate this method by introducing edits into the sodium efflux channel PfATP4 by transfection of a purified CRISPR/Cas9-guide RNA ribonucleoprotein complex and a 200-nucleotide single-stranded oligodeoxynucleotide (ssODN) repair template. Analysis of whole genome sequencing data with the variant-finding program MinorityReport confirmed that only the intended edits were made, and growth inhibition assays confirmed that these mutations confer resistance to the antimalarial SJ733. The method described here is ideally suited for the introduction of mutations that confer a fitness advantage under selection conditions, and the novel finding that an ssODN can function as a repair template in P. falciparum could greatly simplify future editing attempts regardless of the nuclease used or the delivery method.

A Variant PfCRT Isoform Can Contribute to Plasmodium falciparum Resistance to the First-Line Partner Drug Piperaquine

Current efforts to reduce the global burden of malaria are threatened by the rapid spread throughout Asia of Plasmodium falciparum resistance to artemisinin-based combination therapies, which includes increasing rates of clinical failure with dihydroartemisinin plus piperaquine (PPQ) in Cambodia. Using zinc finger nuclease-based gene editing, we report that addition of the C101F mutation to the chloroquine (CQ) resistance-conferring PfCRT Dd2 isoform common to Asia can confer PPQ resistance to cultured parasites. Resistance was demonstrated as significantly higher PPQ concentrations causing 90% inhibition of parasite growth (IC₉₀) or 50% parasite killing (50% lethal dose [LD₅₀]). This mutation also reversed Dd2-mediated CQ resistance, sensitized parasites to amodiaquine, quinine, and artemisinin, and conferred amantadine and blasticidin resistance. Using heme fractionation assays, we demonstrate that PPQ causes a buildup of reactive free heme and inhibits the formation of chemically inert hemozoin crystals. Our data evoke inhibition of heme detoxification in the parasite’s acidic digestive vacuole as the primary mode of both the bis-aminoquinoline PPQ and the related 4-aminoquinoline CQ. Both drugs also inhibit hemoglobin proteolysis at elevated concentrations, suggesting an additional mode of action. Isogenic lines differing in their pfmdr1 copy number showed equivalent PPQ susceptibilities. We propose that mutations in PfCRT could contribute to a multifactorial basis of PPQ resistance in field isolates.

The global agenda to eliminate malaria depends on the continued success of artemisinin-based combination therapies (ACTs), which target the asexual blood stages of the intracellular parasite Plasmodium. Partial resistance to artemisinin, however, is now established in Southeast Asia, exposing the partner drugs to increased selective pressure. Plasmodium falciparum resistance to the first-line partner piperaquine (PPQ) is now spreading rapidly in Cambodia, resulting in clinical treatment failures. Here, we report that a variant form of the Plasmodium falciparum chloroquine resistance transporter, harboring a C101F mutation edited into the chloroquine (CQ)-resistant Dd2 isoform prevalent in Asia, can confer PPQ resistance in cultured parasites. This was accompanied by a loss of CQ resistance. Biochemical assays showed that PPQ, like CQ, inhibits the detoxification of reactive heme that is formed by parasite-mediated catabolism of host hemoglobin. We propose that novel PfCRT variants emerging in the field could contribute to a multigenic basis of PPQ resistance.

Plasmodium falciparum K13 Mutations Differentially Impact Ozonide Susceptibility and Parasite Fitness In Vitro

The emergence and spread in Southeast Asia of Plasmodium falciparum resistance to artemisinin (ART) derivatives, the cornerstone of first-line artemisinin-based combination therapies (ACTs), underscore the urgent need to identify suitable replacement drugs. Discovery and development efforts have identified a series of ozonides with attractive chemical and pharmacological properties that are being touted as suitable replacements. Partial resistance to ART, defined as delayed parasite clearance in malaria patients treated with an ART derivative or an ACT, has been associated with mutations in the P. falciparum K13 gene. In light of reports showing that ART derivatives and ozonides share similar modes of action, we have investigated whether parasites expressing mutant K13 are cross-resistant to the ozonides OZ439 (artefenomel) and OZ227 (arterolane). This work used a panel of culture-adapted clinical isolates from Cambodia that were genetically edited to express variant forms of K13. Phenotypic analyses employed ring-stage survival assays (ring-stage survival assay from 0 to 3 h [RSA_0–3h]), whose results have earlier been shown to correlate with parasite clearance rates in patients. Our results document cross-resistance between OZ277 and dihydroartemisinin (DHA), a semisynthetic derivative of ART, in parasites carrying the K13 mutations C580Y, R539T, and I543T. For OZ439, we observed cross-resistance only for parasites that carried the rare K13 I543T mutation, with no evidence of cross-resistance afforded by the prevalent C580Y mutation. Mixed-culture competition experiments with isogenic lines carrying modified K13 revealed variable growth deficits depending on the K13 mutation and parasite strain and provide a rationale for the broad dissemination of the fitness-neutral K13 C580Y mutation throughout strains currently circulating in Southeast Asia.

ACTs have helped halve the malaria disease burden in recent years; however, emerging resistance to ART derivatives threatens to reverse this substantial progress. Resistance is driven primarily by mutations in the P. falciparum K13 gene. These mutations pose a threat to ozonides, touted as promising alternatives to ARTs that share a similar mode of action. We report that DHA was considerably more potent than OZ439 and OZ277 against ART-sensitive asexual blood-stage parasites cultured in vitro. We also document that mutant K13 significantly compromised the activity of the registered drug OZ277. In contrast, OZ439 remained effective against most parasite lines expressing mutant K13, with the exception of I543T that merits further monitoring in field-based OZ439 efficacy studies. K13 mutations differed considerably in their impact on parasite growth rates, in a strain-dependent context, with the most prevalent C580Y mutation being fitness neutral in recently culture-adapted strains from Cambodia, the epicenter of emerging ART resistance.

Toward a quantitative understanding of antibiotic resistance evolution

The rising prevalence of antibiotic resistant bacteria is an increasingly serious public health challenge. To address this problem, recent work ranging from clinical studies to theoretical modeling has provided valuable insights into the mechanisms of resistance, its emergence and spread, and ways to counteract it. A deeper understanding of the underlying dynamics of resistance evolution will require a combination of experimental and theoretical expertise from different disciplines and new technology for studying evolution in the laboratory. Here, we review recent advances in the quantitative understanding of the mechanisms and evolution of antibiotic resistance. We focus on key theoretical concepts and new technology that enables well-controlled experiments. We further highlight key challenges that can be met in the near future to ultimately develop effective strategies for combating resistance.

Evolution of Fitness Cost-Neutral Mutant PfCRT Conferring P. falciparum 4-Aminoquinoline Drug Resistance Is Accompanied by Altered Parasite Metabolism and Digestive Vacuole Physiology

Southeast Asia is an epicenter of multidrug-resistant Plasmodium falciparum strains. Selective pressures on the subcontinent have recurrently produced several allelic variants of parasite drug resistance genes, including the P. falciparum chloroquine resistance transporter (pfcrt). Despite significant reductions in the deployment of the 4-aminoquinoline drug chloroquine (CQ), which selected for the mutant pfcrt alleles that halted CQ efficacy decades ago, the parasite pfcrt locus is continuously evolving. This is highlighted by the presence of a highly mutated allele, Cam734 pfcrt, which has acquired the singular ability to confer parasite CQ resistance without an associated fitness cost. Here, we used pfcrt-specific zinc-finger nucleases to genetically dissect this allele in the pathogenic setting of asexual blood-stage infection. Comparative analysis of drug resistance and growth profiles of recombinant parasites that express Cam734 or variants thereof, Dd2 (the most common Southeast Asian variant), or wild-type pfcrt, revealed previously unknown roles for PfCRT mutations in modulating parasite susceptibility to multiple antimalarial agents. These results were generated in the GC03 strain, used in multiple earlier pfcrt studies, and might differ in natural isolates harboring this allele. Results presented herein show that Cam734-mediated CQ resistance is dependent on the rare A144F mutation that has not been observed beyond Southeast Asia, and reveal distinct impacts of this and other Cam734-specific mutations on CQ resistance and parasite growth rates. Biochemical assays revealed a broad impact of mutant PfCRT isoforms on parasite metabolism, including nucleoside triphosphate levels, hemoglobin catabolism and disposition of heme, as well as digestive vacuole volume and pH. Results from our study provide new insights into the complex molecular basis and physiological impact of PfCRT-mediated antimalarial drug resistance, and inform ongoing efforts to characterize novel pfcrt alleles that can undermine the efficacy of first-line antimalarial drug regimens.

Point mutations in the Plasmodium falciparum chloroquine resistance transporter (PfCRT) earlier thwarted the clinical efficacy of chloroquine, the former gold standard, and constitute a major determinant of parasite susceptibility to antimalarial drugs. Recently, we reported that the highly mutated Cambodian PfCRT isoform Cam734 is fitness-neutral in terms of parasite growth, unlike other less fit isoforms such as Dd2 that are outcompeted by wild-type parasites in the absence of CQ pressure. Using pfcrt-specific zinc-finger nucleases to genetically dissect the Cam734 allele, we report that its unique constituent mutations directly contribute to CQ resistance and collectively offset fitness costs associated with intermediate mutational steps. We also report that these mutations can contribute to resistance or increased sensitivity to multiple first-line partner drugs. Using isogenic parasite lines, we provide evidence of changes in parasite metabolism associated with the Cam734 allele compared to Dd2. We also observe a close correlation between CQ inhibition of hemozoin formation and parasite growth, and provide evidence that Cam734 PfCRT can modulate drug potency depending on its membrane electrochemical gradient. Our data highlight the capacity of PfCRT to evolve new states of antimalarial drug resistance and to offset associated fitness costs through its impact on parasite physiology and hemoglobin catabolism.

Genome-Wide Association Studies of Drug-Resistance Determinants

Population genetic strategies that leverage association, selection, and linkage have identified drug-resistant loci. However, challenges and limitations persist in identifying drug-resistance loci in malaria. In this review we discuss the genetic basis of drug resistance and the use of genome-wide association studies, complemented by selection and linkage studies, to identify and understand mechanisms of drug resistance and response. We also discuss the implications of nongenetic mechanisms of drug resistance recently reported in the literature, and present models of the interplay between nongenetic and genetic processes that contribute to the emergence of drug resistance. Throughout, we examine artemisinin resistance as an example to emphasize challenges in identifying phenotypes suitable for population genetic studies as well as complications due to multiple-factor drug resistance.

Drug resistance in eukaryotic microorganisms

Eukaryotic microbial pathogens are major contributors to illness and death globally. Although much of their impact can be controlled by drug therapy as with prokaryotic microorganisms, the emergence of drug resistance has threatened these treatment efforts. Here, we discuss the challenges posed by eukaryotic microbial pathogens and how these are similar to, or differ from, the challenges of prokaryotic antibiotic resistance. The therapies used for several major eukaryotic microorganisms are then detailed, and the mechanisms that they have evolved to overcome these therapies are described. The rapid emergence of resistance and the restricted pipeline of new drug therapies pose considerable risks to global health and are particularly acute in the developing world. Nonetheless, we detail how the integration of new technology, biological understanding, epidemiology and evolutionary analysis can help sustain existing therapies, anticipate the emergence of resistance or optimize the deployment of new therapies.

CRISPR-Cas9-modified pfmdr1 protects Plasmodium falciparum asexual blood stages and gametocytes against a class of piperazine-containing compounds but potentiates artemisinin-based combination therapy partner drugs

Emerging resistance to first-line antimalarial combination therapies threatens malaria treatment and the global elimination campaign. Improved therapeutic strategies are required to protect existing drugs and enhance treatment efficacy. We report that the piperazine-containing compound ACT-451840 exhibits single-digit nanomolar inhibition of the Plasmodium falciparum asexual blood stages and transmissible gametocyte forms. Genome sequence analyses of in vitro-derived ACT-451840-resistant parasites revealed single nucleotide polymorphisms in pfmdr1, which encodes a digestive vacuole membrane-bound ATP-binding cassette transporter known to alter P. falciparum susceptibility to multiple first-line antimalarials. CRISPR-Cas9 based gene editing confirmed that PfMDR1 point mutations mediated ACT-451840 resistance. Resistant parasites demonstrated increased susceptibility to the clinical drugs lumefantrine, mefloquine, quinine and amodiaquine. Stage V gametocytes harboring Cas9-introduced pfmdr1 mutations also acquired ACT-451840 resistance. These findings reveal that PfMDR1 mutations can impart resistance to compounds active against asexual blood stages and mature gametocytes. Exploiting PfMDR1 resistance mechanisms provides new opportunities for developing disease-relieving and transmission-blocking antimalarials.

Evidence of a Mild Mutator Phenotype in Cambodian Plasmodium falciparum Malaria Parasites

Malaria control efforts have been continuously stymied by drug-resistant strains of Plasmodium falciparum, which typically originate in Southeast Asia prior to spreading into high-transmission settings in Africa. One earlier proposed explanation for Southeast Asia being a hotbed of resistance has been the hypermutability or "Accelerated Resistance to Multiple Drugs" (ARMD) phenotype, whereby multidrug-resistant Southeast Asian parasites were reported to exhibit 1,000-fold higher rates of resistance to unrelated antimalarial agents when compared to drug-sensitive parasites. However, three recent studies do not recapitulate this hypermutability phenotype. Intriguingly, genome sequencing of recently derived multidrug-resistant Cambodian isolates has identified a high proportion of DNA repair gene mutations in multidrug-resistant parasites, suggesting their potential role in shaping local parasite evolution. By adapting fluctuation assays for use in P. falciparum, we have examined the in vitro mutation rates of five recent Cambodian isolates and three reference laboratory strains. For these studies we also generated a knockout parasite line lacking the DNA repair factor Exonuclease I. In these assays, parasites were typed for their ability to acquire resistance to KAE609, currently in advanced clinical trials, yielding 13 novel mutations in the Na+/H+-ATPase PfATP4, the primary resistance determinant. We observed no evidence of hypermutability. Instead, we found evidence of a mild mutator (up to a 3.4-fold increase in mutation rate) phenotype in two artemisinin-resistant Cambodian isolates, which carry DNA repair gene mutations. We observed that one such mutation in the Mismatch Repair protein Mlh1 contributes to the mild mutator phenotype when modeled in yeast (scmlh1-P157S). Compared to basal rates of mutation, a mild mutator phenotype may provide a greater overall benefit for parasites in Southeast Asia in terms of generating drug resistance without incurring detrimental fitness costs.