Association of a Novel Mutation in the Plasmodium falciparum Chloroquine Resistance Transporter With Decreased Piperaquine Sensitivity

Strong positive selection biases identity-by-descent-based inferences of recent demography and population structure in Plasmodium falciparum

Malaria genomic surveillance often estimates parasite genetic relatedness using metrics such as Identity-By-Decent (IBD), yet strong positive selection stemming from antimalarial drug resistance or other interventions may bias IBD-based estimates. In this study, we use simulations, a true IBD inference algorithm, and empirical data sets from different malaria transmission settings to investigate the extent of this bias and explore potential correction strategies. We analyze whole genome sequence data generated from 640 new and 3089 publicly available Plasmodium falciparum clinical isolates. We demonstrate that positive selection distorts IBD distributions, leading to underestimated effective population size and blurred population structure. Additionally, we discover that the removal of IBD peak regions partially restores the accuracy of IBD-based inferences, with this effect contingent on the population's background genetic relatedness and extent of inbreeding. Consequently, we advocate for selection correction for parasite populations undergoing strong, recent positive selection, particularly in high malaria transmission settings.

Assessing emergence risk of double-resistant and triple-resistant genotypes of Plasmodium falciparum

Delaying and slowing antimalarial drug resistance evolution is a priority for malaria-endemic countries. Until novel therapies become available, the mainstay of antimalarial treatment will continue to be artemisinin-based combination therapy (ACT). Deployment of different ACTs can be optimized to minimize evolutionary pressure for drug resistance by deploying them as a set of co-equal multiple first-line therapies (MFT) rather than rotating therapies in and out of use. Here, we consider one potential detriment of MFT policies, namely, that the simultaneous deployment of multiple ACTs could drive the evolution of different resistance alleles concurrently and that these resistance alleles could then be brought together by recombination into double-resistant or triple-resistant parasites. Using an individual-based model, we compare MFT and cycling policies in malaria transmission settings ranging from 0.1% to 50% prevalence. We define a total risk measure for multi-drug resistance (MDR) by summing the area under the genotype-frequency curves (AUC) of double- and triple-resistant genotypes. When prevalence ≥ 1%, total MDR risk ranges from statistically similar to 80% lower under MFT policies than under cycling policies, irrespective of whether resistance is imported or emerges de novo. At 0.1% prevalence, there is little statistical difference in MDR risk between MFT and cycling.

The problem of antimalarial resistance and its implications for drug discovery

INTRODUCTION

Malaria remains a devastating infectious disease with hundreds of thousands of casualties each year. Antimalarial drug resistance has been a threat to malaria control and elimination for many decades and is still of concern today. Despite the continued effectiveness of current first-line treatments, namely artemisinin-based combination therapies, the emergence of drug-resistant parasites in Southeast Asia and even more alarmingly the occurrence of resistance mutations in Africa is of great concern and requires immediate attention.

AREAS COVERED

A comprehensive overview of the mechanisms underlying the acquisition of drug resistance in Plasmodium falciparum is given. Understanding these processes provides valuable insights that can be harnessed for the development and selection of novel antimalarials with reduced resistance potential. Additionally, strategies to mitigate resistance to antimalarial compounds on the short term by using approved drugs are discussed.

EXPERT OPINION

While employing strategies that utilize already approved drugs may offer a prompt and cost-effective approach to counter antimalarial drug resistance, it is crucial to recognize that only continuous efforts into the development of novel antimalarial drugs can ensure the successful treatment of malaria in the future. Incorporating resistance propensity assessment during this developmental process will increase the likelihood of effective and enduring malaria treatments.

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.

Antimalarial drug discovery: progress and approaches

Recent antimalarial drug discovery has been a race to produce new medicines that overcome emerging drug resistance, whilst considering safety and improving dosing convenience. Discovery efforts have yielded a variety of new molecules, many with novel modes of action, and the most advanced are in late-stage clinical development. These discoveries have led to a deeper understanding of how antimalarial drugs act, the identification of a new generation of drug targets, and multiple structure-based chemistry initiatives. The limited pool of funding means it is vital to prioritize new drug candidates. They should exhibit high potency, a low propensity for resistance, a pharmacokinetic profile that favours infrequent dosing, low cost, preclinical results that demonstrate safety and tolerability in women and infants, and preferably the ability to block Plasmodium transmission to Anopheles mosquito vectors. In this Review, we describe the approaches that have been successful, progress in preclinical and clinical development, and existing challenges. We illustrate how antimalarial drug discovery can serve as a model for drug discovery in diseases of poverty.

Antimalarial mass drug administration in large populations and the evolution of drug resistance

Mass drug administration (MDA) with antimalarials has been shown to reduce prevalence and interrupt transmission in small populations, in populations with reliable access to antimalarial drugs, and in populations where sustained improvements in diagnosis and treatment are possible. In addition, when MDA is effective it eliminates both drug-resistant parasites and drug-sensitive parasites, which has the long-term benefit of extending the useful therapeutic life of first-line therapies for all populations, not just the focal population where MDA was carried out. However, in order to plan elimination measures effectively, it is necessary to characterize the conditions under which failed MDA could exacerbate resistance. We use an individual-based stochastic model of Plasmodium falciparum transmission to evaluate this risk for MDA using dihydroartemisinin-piperaquine (DHA-PPQ), in populations where access to antimalarial treatments may not be uniformly high and where re-importation of drug-resistant parasites may be common. We find that artemisinin-resistance evolution at the kelch13 locus can be accelerated by MDA when all three of the following conditions are met: (1) strong genetic bottlenecking that falls short of elimination, (2) re-importation of artemisinin-resistant genotypes, and (3) continued selection pressure during routine case management post-MDA. Accelerated resistance levels are not immediate but follow the rebound of malaria cases post-MDA, if this is allowed to occur. Crucially, resistance is driven by the selection pressure during routine case management post-MDA and not the selection pressure exerted during the MDA itself. Second, we find that increasing treatment coverage post-MDA increases the probability of local elimination in low-transmission regions (prevalence < 2%) in scenarios with both low and high levels of drug-resistance importation. This emphasizes the importance of planning for and supporting high coverage of diagnosis and treatment post-MDA.

pH-dependence of the Plasmodium falciparum chloroquine resistance transporter is linked to the transport cycle

The chloroquine resistance transporter, PfCRT, of the human malaria parasite Plasmodium falciparum is sensitive to acidic pH. Consequently, PfCRT operates at 60% of its maximal drug transport activity at the pH of 5.2 of the digestive vacuole, a proteolytic organelle from which PfCRT expels drugs interfering with heme detoxification. Here we show by alanine-scanning mutagenesis that E207 is critical for pH sensing. The E207A mutation abrogates pH-sensitivity, while preserving drug substrate specificity. Substituting E207 with Asp or His, but not other amino acids, restores pH-sensitivity. Molecular dynamics simulations and kinetics analyses suggest an allosteric binding model in which PfCRT can accept both protons and chloroquine in a partial noncompetitive manner, with increased proton concentrations decreasing drug transport. Further simulations reveal that E207 relocates from a peripheral to an engaged location during the transport cycle, forming a salt bridge with residue K80. We propose that the ionized carboxyl group of E207 acts as a hydrogen acceptor, facilitating transport cycle progression, with pH sensing as a by-product.

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.

In vitro activity of rhinacanthin analogues against drug resistant Plasmodium falciparum isolates from Northeast Thailand

BACKGROUND

New anti-malarial drugs are needed urgently to address the increasing challenges of drug-resistant falciparum malaria. Two rhinacanthin analogues containing a naphthoquinone moiety resembling atovaquone showed promising in-vitro activity against a P. falciparum laboratory reference strain (K1). The anti-malarial activity of these 2 compounds was further evaluated for P. falciparum field isolates from an area of multi-drug resistance in Northeast Thailand.

METHODS

Using a pLDH enzyme-linked immunosorbent assay, four P. falciparum isolates from Northeast Thailand in 2018 were tested for in vitro sensitivity to the two synthetic rhinacanthin analogues 1 and 2 as well as established anti-malarials. Mutations in the P. falciparum cytochrome b gene, a marker for atovaquone (ATQ) resistance, were genotyped in all four field isolates as well as 100 other clinical isolates from the same area using PCR-artificial Restriction Fragment Length Polymorphisms. Pfkelch13 mutations, a marker for artemisinin (ART) resistance, were also examined in all isolates.

RESULTS

The 50% inhibitory concentrations (IC₅₀) of P. falciparum field isolates for rhinacanthin analogue 1 was 321.9-791.1 nM (median = 403.1 nM). Parasites were more sensitive to analogue 2: IC₅₀ 48.6-63.3 nM (median = 52.2 nM). Similar results were obtained against P. falciparum reference laboratory strains 3D7 and W2. The ART-resistant IPC-5202 laboratory strain was more sensitive to these compounds with a median IC₅₀ 45.9 and 3.3 nM for rhinacanthin analogues 1 and 2, respectively. The ATQ-resistant C2B laboratory strain showed high-grade resistance towards both compounds (IC₅₀ > 15,000 nM), and there was a strong positive correlation between the IC₅₀ values for these compounds and ATQ (r = 0.83-0.97, P < 0.001). There were no P. falciparum cytochrome b mutations observed in the field isolates, indicating that P. falciparum isolates from this area remained ATQ-sensitive. Pfkelch13 mutations and the ring-stage survival assay confirmed that most isolates were resistant to ART.

CONCLUSIONS

Two rhinacanthin analogues showed parasiticidal activity against multi-drug resistant P. falciparum isolates, although less potent than ATQ. Rhinacanthin analogue 2 was more potent than analogue 1, and can be a lead compound for further optimization as an anti-malarial in areas with multidrug resistance.

Molecular surveillance of artemisinin resistance-related Pfk13 and pfcrt polymorphisms in imported Plasmodium falciparum isolates reported in eastern China from 2015 to 2019

BACKGROUND

Artemisinin-based combination therapy (ACT) has been recommended as the first-line treatment by the World Health Organization to treat uncomplicated Plasmodium falciparum malaria. However, the emergence and spread of P. falciparum resistant to artemisinins and their partner drugs is a significant risk for the global effort to reduce disease burden facing the world. Currently, dihydroartemisinin-piperaquine (DHA-PPQ) is the most common drug used to treat P. falciparum, but little evidence about the resistance status targeting DHA (ACT drug) and its partner drug (PPQ) has been reported in Shandong Province, China.

METHODS

A retrospective study was conducted to explore the prevalence and spatial distribution of Pfk13 and Pfcrt polymorphisms (sites of 72-76, and 93-356) among imported P. falciparum isolates between years 2015-2019 in Shandong Province in eastern China. Individual epidemiological information was collected from a web-based reporting system were reviewed and analysed.

RESULTS

A total of 425 P. falciparum blood samples in 2015-2019 were included and 7.3% (31/425) carried Pfk13 mutations. Out of the isolates that carried Pfk13 mutations, 54.8% (17/31) were nonsynonymous polymorphisms. The mutant alleles A578S, Q613H, C469C, and S549S in Pfk13 were the more frequently detected allele, the mutation rate was the same as 9.7% (3/31). Another allele Pfk13 C580Y, closely associated with artemisinin (ART) resistance, was found as 3.2% (2/31), which was found in Cambodia. A total of 14 mutant isolates were identified in Western Africa countries (45.2%, 14/31). For the Pfcrt gene, the mutation rate was 18.1% (77/425). T₇₆T₃₅₆ and T₇₆ were more frequent in all 13 different haplotypes with 26.0% (20/77) and 23.4% (18/77). The CVIET and CVIKT mutant at loci 72-76 have exhibited a prevalence of 19.5% (15/77) and 3.9% (3/77), respectively. The CVIET was mainly observed in samples from Congo (26.7%, 4/15) and Mozambique (26.7%, 4/15). No mutations were found at loci 97, 101 and 145. For polymorphisms at locus 356, a total of 24 isolates were identified and mainly from Congo (29.2%, 7/24).

CONCLUSION

These findings indicate a low prevalence of Pfk13 in the African isolates. However, the emergence and increase in the new alleles Pfcrt I356T, reveals a potential risk of drug pressure in PPQ among migrant workers returned from Africa. Therefore, continuous molecular surveillance of Pfcrt mutations and in vitro susceptibility tests related to PPQ are necessary.

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

Background

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

Methods

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

Results

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

Conclusions

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

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.

Molecular markers of dihydroartemisinin-piperaquine resistance in northwestern Thailand

BACKGROUND

Dihydroartemisinin-piperaquine (DHA-PPQ) combination therapy is the current first-line treatment for Plasmodium falciparum malaria in Thailand. Since its introduction in 2015, resistance to this drug combination has emerged in the eastern part of the Greater Mekong Subregion including the eastern part of Thailand near Cambodia. This study aimed to assess whether the resistance genotypes have arisen the western part of country.

METHODS

Fifty-seven P. falciparum-infected blood samples were collected in Tak province of northwestern Thailand between 2013 and 2019. Resistance to DHA was examined through the single nucleotide polymorphisms (SNPs) of kelch13. PPQ resistance was examined through the copy number plasmepsin-2 and the SNPs of Pfcrt.

RESULTS

Among the samples whose kelch13 were successfully sequenced, approximately half (31/55; 56%) had mutation associated with artemisinin resistance, including G533S (23/55; 42%), C580Y (6/55; 11%), and G538V (2/55; 4%). During the study period, G533S mutation appeared and increased from 20% (4/20) in 2014 to 100% (9/9) in 2019. No plasmepsin-2 gene amplification was observed, but one sample (1/54) had the Pfcrt F145I mutation previously implicated in PPQ resistance.

CONCLUSIONS

Kelch13 mutation was common in Tak Province in 2013-2019. A new mutation G533S emerged in 2014 and rose to dominance in 2019. PPQ resistance marker Pfcrt F145I was also detected in 2019. Continued surveillance of treatment efficacy and drug resistance markers is warranted.

Piperaquine-resistant PfCRT mutations differentially impact drug transport, hemoglobin catabolism and parasite physiology in Plasmodium falciparum asexual blood stages

The emergence of Plasmodium falciparum parasite resistance to dihydroartemisinin + piperaquine (PPQ) in Southeast Asia threatens plans to increase the global use of this first-line antimalarial combination. High-level PPQ resistance appears to be mediated primarily by novel mutations in the P. falciparum chloroquine resistance transporter (PfCRT), which enhance parasite survival at high PPQ concentrations in vitro and increase the risk of dihydroartemisinin + PPQ treatment failure in patients. Using isogenic Dd2 parasites expressing contemporary pfcrt alleles with differential in vitro PPQ susceptibilities, we herein characterize the molecular and physiological adaptations that define PPQ resistance in vitro. Using drug uptake and cellular heme fractionation assays we report that the F145I, M343L, and G353V PfCRT mutations differentially impact PPQ and chloroquine efflux. These mutations also modulate proteolytic degradation of host hemoglobin and the chemical inactivation of reactive heme species. Peptidomic analyses reveal significantly higher accumulation of putative hemoglobin-derived peptides in the PPQ-resistant mutant PfCRT isoforms compared to parental PPQ-sensitive Dd2. Joint transcriptomic and metabolomic profiling of late trophozoites from PPQ-resistant or -sensitive isogenic lines reveals differential expression of genes involved in protein translation and cellular metabolism. PPQ-resistant parasites also show increased susceptibility to an inhibitor of the P. falciparum M17 aminopeptidase that operates on short globin-derived peptides. These results reveal unique physiological changes caused by the gain of PPQ resistance and highlight the potential therapeutic value of targeting peptide metabolism in P. falciparum.

In vitro susceptibility profile of Plasmodium falciparum clinical isolates from Ghana to antimalarial drugs and polymorphisms in resistance markers

Drug resistance in Plasmodium falciparum compromises the effectiveness of antimalarial therapy. This study aimed to evaluate the extent of drug resistance in parasites obtained from international travelers returning from Ghana to guide the management of malaria cases. Eighty-two clinical parasite isolates were obtained from patients returning from Ghana in 2016–2018, of which 29 were adapted to continuous in vitro culture. Their geometric mean IC₅₀ values to a panel of 11 antimalarial drugs, assessed using the standard SYBR Green-I drug sensitivity assay, were 2.1, 3.8, 1.0, 2.7, 17.2, 4.6, 8.3, 8.3, 19.6, 55.1, and 11,555 nM for artemether, artesunate, dihydroartemisinin, lumefantrine, mefloquine, piperaquine, naphthoquine, pyronaridine, chloroquine, quinine, and pyrimethamine, respectively. Except for chloroquine and pyrimethamine, the IC₅₀ values for other tested drugs were below the resistance threshold. The mean ring-stage survival assay value was 0.8%, with four isolates exceeding 1%. The mean piperaquine survival assay value was 2.1%, all below 10%. Mutations associated with chloroquine resistance (pfcrt K76T and pfmdr1 N86Y) were scarce, consistent with the discontinuation of chloroquine a decade ago. Instead, the pfmdr1 86N-184F-1246D haplotype was predominant, suggesting selection by the extensive use of artemether-lumefantrine. No mutations in the pfk13 propeller domain were detected. The pfdhfr/pfdhps quadruple mutant IRNGK associated with resistance to sulfadoxine-pyrimethamine reached an 82% prevalence. In addition, five isolates had pfgch1 gene amplification but, intriguingly, increased susceptibilities to pyrimethamine. This study showed that parasites originating from Ghana were susceptible to artemisinins and the partner drugs of artemisinin-based combination therapies. Genotyping drug resistance genes identified the signature of selection by artemether-lumefantrine. Parasites showed substantial levels of resistance to the antifolate drugs. Continuous resistance surveillance is necessary to guide timely changes in drug policy.

Multidisciplinary Investigations of Sustained Malaria Transmission in the Greater Mekong Subregion

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

Evaluation of the Effect of Gene Duplication by Genome Editing on Drug Resistance in Plasmodium falciparum

The emergence and spread of drug-resistant Plasmodium falciparum have compromised antimalarial efficacy and threatened the global malaria elimination campaign using artemisinin combination therapies. The impacts of amino acid substitutions in antimalarial drug resistance-associated genes on drug susceptibility have been investigated; however, the effects of amplification of those genes remain unexplored due to the lack of robust genetic approaches. Here, we generated transgenic P. falciparum parasites with an additional copy of a drug resistance-associated gene using the highly efficient CRISPR/Cas9 system and investigated their drug response. Insertion of a drug resistance-associated gene expression cassette in the genome resulted in a roughly twofold increase of mRNA levels of the target gene mdr1, which encodes multidrug resistance protein 1. The gene duplication event contributed to resistance to mefloquine, lumefantrine, and dihydroartemisinin, while the duplication of a genomic region encoding plasmepsin 2 and plasmepsin 3 did not affect resistance to antimalarial drugs, including piperaquine. We further demonstrated that mdr1 mRNA expression levels are strongly associated with mefloquine resistance in several field-derived P. falciparum lines with various genetic backgrounds. This study provides compelling evidence that mdr1 could serve as a molecular marker for the surveillance of mefloquine-resistant parasites. Long DNA integration into parasite genomes using the CRISPR/Cas9 system provides a useful tool for the evaluation of the effect of copy number variation on drug response.

How has mass drug administration with dihydroartemisinin-piperaquine impacted molecular markers of drug resistance? A systematic review

The World Health Organization (WHO) recommends surveillance of molecular markers of resistance to anti-malarial drugs. This is particularly important in the case of mass drug administration (MDA), which is endorsed by the WHO in some settings to combat malaria. Dihydroartemisinin-piperaquine (DHA-PPQ) is an artemisinin-based combination therapy which has been used in MDA. This review analyses the impact of MDA with DHA-PPQ on the evolution of molecular markers of drug resistance. The review is split into two parts. Section I reviews the current evidence for different molecular markers of resistance to DHA-PPQ. This includes an overview of the prevalence of these molecular markers in Plasmodium falciparum Whole Genome Sequence data from the MalariaGEN Pf3k project. Section II is a systematic literature review of the impact that MDA with DHA-PPQ has had on the evolution of molecular markers of resistance. This systematic review followed PRISMA guidelines. This review found that despite being a recognised surveillance tool by the WHO, the surveillance of molecular markers of resistance following MDA with DHA-PPQ was not commonly performed. Of the total 96 papers screened for eligibility in this review, only 20 analysed molecular markers of drug resistance. The molecular markers published were also not standardized. Overall, this warrants greater reporting of molecular marker prevalence following MDA implementation. This should include putative pfcrt mutations which have been found to convey resistance to DHA-PPQ in vitro.

Plasmodium falciparum Drug Resistance Genes pfmdr1 and pfcrt In Vivo Co-Expression During Artemether-Lumefantrine Therapy

Background: Artemisinin-based combination therapies (ACTs) are the global mainstay treatment of uncomplicated Plasmodium falciparum infections. PfMDR1 and PfCRT are two transmembrane transporters, associated with sensitivity to several antimalarials, found in the parasite food vacuole. Herein, we explore if their relatedness extends to overlapping patterns of gene transcriptional activity before and during ACT administration.

Methods: In a clinical trial performed in Tanzania, we explored the pfmdr1 and pfcrt transcription levels from 48 patients with uncomplicated P. falciparum malaria infections who underwent treatment with artemether-lumefantrine (AL). Samples analyzed were collected before treatment initiation and during the first 24 h of treatment. The frequency of PfMDR1 N86Y and PfCRT K76T was determined through PCR-RFLP or direct amplicon sequencing. Gene expression was analyzed by real-time quantitative PCR.

Results: A wide range of pre-treatment expression levels was observed for both genes, approximately 10-fold for pfcrt and 50-fold for pfmdr1. In addition, a significant positive correlation demonstrates pfmdr1 and pfcrt co-expression. After AL treatment initiation, pfmdr1 and pfcrt maintained the positive co-expression correlation, with mild downregulation throughout the 24 h post-treatment. Additionally, a trend was observed for PfMDR1 N86 alleles and higher expression before treatment initiation.

Conclusion:pfmdr1 and pfcrt showed significant co-expression patterns in vivo, which were generally maintained during ACT treatment. This observation points to relevant related roles in the normal parasite physiology, which seem essential to be maintained when the parasite is exposed to drug stress. In addition, keeping the simultaneous expression of both transporters might be advantageous for responding to the drug action.

A single point mutation in the Plasmodium falciparum 3'-5' exonuclease does not alter piperaquine susceptibility

Background

The rise in Plasmodium falciparum resistance to dihydroartemisinin–piperaquine (DHA–PPQ) treatment has been documented in the Greater Mekong Subregion with associations with mutations in the P. falciparum chloroquine resistance transporter (pfcrt) and plasmepsin 2 (pfpm2) genes. However, it is unclear whether other genes also play a role with PPQ resistance, such as the E415G mutation in the exonuclease (pfexo) gene. The aim of this study was to investigate the role of this mutation in PPQ resistance by generating transgenic parasites expressing the pfexo-E415G mutant allele.

Methods

Transgenic parasite clones carrying the E415G mutation in PfEXO of the B5 isolate were derived by CRISPR-Cas9 gene editing and verified using PCR and gene sequencing. Polymorphisms of pfkelch-13, pfcrt, and pfexo were examined by PCR while the copy number variations of pfpm2 were examined by both relative quantitative real-time PCR and the duplication breakpoint assay. Drug sensitivity against a panel of antimalarials, the ring-stage survival assay (RSA), the PPQ survival assay (PSA), and bimodal dose-response curves were used to evaluate antimalarial susceptibility.

Results

The transgenic line, B5-rexo-E415G-B8, was successfully generated. The PPQ-IC₉₀, %PPQ survival, and the bimodal dose-response clearly showed that E415G mutation in PfEXO of B5 isolate remained fully susceptible to PPQ. Furthermore, growth assays demonstrated that the engineered parasites grew slightly faster than the unmodified parental isolates whereas P. falciparum isolates harbouring pfkelch-13, pfcrt, and pfexo mutations with multiple copies of pfpm2 grew much more slowly.

Conclusions

Insertion of the E415G mutation in PfEXO did not lead to increased PPQ-IC₉₀ and %PPQ survival, suggesting that this mutation alone may not be associated with PPQ resistance, but could still be an important marker if used in conjunction with other markers for monitoring PPQ-resistant parasites. The results also highlight the importance of monitoring and evaluating suspected genetic mutations with regard to parasite fitness and resistance.

Supplementary information

The online version contains supplementary material available at 10.1186/s12936-022-04148-z.

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

Background

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

Methods

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

Results

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

Conclusions

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

Graphical Abstract

Supplementary Information

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

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.

Longitudinal ex vivo and molecular trends of chloroquine and piperaquine activity against Plasmodium falciparum and P. vivax before and after introduction of artemisinin-based combination therapy in Papua, Indonesia

Drug resistant Plasmodium parasites are a major threat to malaria control and elimination. After reports of high levels of multidrug resistant P. falciparum and P. vivax in Indonesia, in 2005, the national first-line treatment policy for uncomplicated malaria was changed in March 2006, to dihydroartemisinin-piperaquine against all species. This study assessed the temporal trends in ex vivo drug susceptibility to chloroquine (CQ) and piperaquine (PIP) for both P. falciparum and P. vivax clinical isolates collected between 2004 and 2018, by using schizont maturation assays, and genotyped a subset of isolates for known and putative molecular markers of CQ and PIP resistance by using Sanger and next generation whole genome sequencing. The median CQ IC50 values varied significantly between years in both Plasmodium species, but there was no significant trend over time. In contrast, there was a significant trend for increasing PIP IC50s in both Plasmodium species from 2010 onwards. Whereas the South American CQ resistant 7G8 pfcrt SVMNT isoform has been fixed since 2005 in the study area, the pfmdr1 86Y allele frequencies decreased and became fixed at the wild-type allele in 2015. In P. vivax isolates, putative markers of CQ resistance (no pvcrt-o AAG (K10) insertion and pvmdr1 Y967F and F1076L) were fixed at the mutant alleles since 2005. None of the putative PIP resistance markers were detected in P. falciparum. The ex vivo drug susceptibility and molecular analysis of CQ and PIP efficacy for P. falciparum and P. vivax after 12 years of intense drug pressure with DHP suggests that whilst the degree of CQ resistance appears to have been sustained, there has been a slight decline in PIP susceptibility, although this does not appear to have reached clinically significant levels. The observed decreasing trend in ex vivo PIP susceptibility highlights the importance of ongoing surveillance.

Protecting future antimalarials from the trap of resistance: Lessons from artemisinin-based combination therapy (ACT) failures

Having faced increased clinical treatment failures with dihydroartemisinin-piperaquine (DHA-PPQ), Cambodia swapped the first line artemisinin-based combination therapy (ACT) from DHA-PPQ to artesunate-mefloquine given that parasites resistant to piperaquine are susceptible to mefloquine. However, triple mutants have now emerged, suggesting that drug rotations may not be adequate to keep resistance at bay. There is, therefore, an urgent need for alternative treatment strategies to tackle resistance and prevent its spread. A proper understanding of all contributors to artemisinin resistance may help us identify novel strategies to keep artemisinins effective until new drugs become available for their replacement. This review highlights the role of the key players in artemisinin resistance, the current strategies to deal with it and suggests ways of protecting future antimalarial drugs from bowing to resistance as their predecessors did.

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

Background

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

Methods

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

Results

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

Conclusions

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

Funding

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

Plasmodium falciparum resistance to ACTs: Emergence, mechanisms, and outlook

Emergence and spread of resistance in Plasmodium falciparum to the frontline treatment artemisinin-based combination therapies (ACTs) in the epicenter of multidrug resistance of Southeast Asia threaten global malaria control and elimination. Artemisinin (ART) resistance (or tolerance) is defined clinically as delayed parasite clearance after treatment with an ART drug. The resistance phenotype is restricted to the early ring stage and can be measured in vitro using a ring-stage survival assay. ART resistance is associated with mutations in the propeller domain of the Kelch family protein K13. As a pro-drug, ART is activated primarily by heme, which is mainly derived from hemoglobin digestion in the food vacuole. Activated ARTs can react promiscuously with a wide range of cellular targets, disrupting cellular protein homeostasis. Consistent with this mode of action for ARTs, the molecular mechanisms of K13-mediated ART resistance involve reduced hemoglobin uptake/digestion and increased cellular stress response. Mutations in other genes such as AP-2μ (adaptor protein-2 μ subunit), UBP-1 (ubiquitin-binding protein-1), and Falcipain 2a that interfere with hemoglobin uptake and digestion also increase resistance to ARTs. ART resistance has facilitated the development of resistance to the partner drugs, resulting in rapidly declining ACT efficacies. The molecular markers for resistance to the partner drugs are mostly associated with point mutations in the two food vacuole membrane transporters PfCRT and PfMDR1, and amplification of pfmdr1 and the two aspartic protease genes plasmepsin 2 and 3. It has been observed that mutations in these genes can have opposing effects on sensitivities to different partner drugs, which serve as the principle for designing triple ACTs and drug rotation. Although clinical ACT resistance is restricted to Southeast Asia, surveillance for drug resistance using in vivo clinical efficacy, in vitro assays, and molecular approaches is required to prevent or slow down the spread of resistant parasites.

Integration of population and functional genomics to understand mechanisms of artemisinin resistance in Plasmodium falciparum

Resistance to antimalarial drugs, and in particular to the artemisinin derivatives and their partner drugs, threatens recent progress toward regional malaria elimination and eventual global malaria eradication. Population-level studies utilizing whole-genome sequencing approaches have facilitated the identification of regions of the parasite genome associated with both clinical and in vitro drug-resistance phenotypes. However, the biological relevance of genes identified in these analyses and the establishment of a causal relationship between genotype and phenotype requires functional characterization. Here we examined data from population genomic and transcriptomic studies in the context of data generated from recent functional studies, using a new population genetic approach designed to identify potential favored mutations within the region of a selective sweep (iSAFE). We identified several genes functioning in pathways now known to be associated with artemisinin resistance that were supported in early population genomic studies, as well as potential new drug targets/pathways for further validation and consideration for treatment of artemisinin-resistant Plasmodium falciparum. In addition, we establish the utility of iSAFE in identifying positively-selected mutations in population genomic studies, potentially accelerating the time to functional validation of candidate genes.

Genomic Approaches to Drug Resistance in Malaria

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

Molecular Mechanisms of Drug Resistance in Plasmodium falciparum Malaria

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

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

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

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

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

Plasmodium falciparum phenotypic and genotypic resistance profile during the emergence of Piperaquine resistance in Northeastern Thailand

Malaria remains a public health problem in Thailand, especially along its borders where highly mobile populations can contribute to persistent transmission. This study aimed to determine resistant genotypes and phenotypes of 112 Plasmodium falciparum isolates from patients along the Thai-Cambodia border during 2013-2015. The majority of parasites harbored a pfmdr1-Y184F mutation. A single pfmdr1 copy number had CVIET haplotype of amino acids 72-76 of pfcrt and no pfcytb mutations. All isolates had a single pfk13 point mutation (R539T, R539I, or C580Y), and increased % survival in the ring-stage survival assay (except for R539I). Multiple copies of pfpm2 and pfcrt-F145I were detected in 2014 (12.8%) and increased to 30.4% in 2015. Parasites containing either multiple pfpm2 copies with and without pfcrt-F145I or a single pfpm2 copy with pfcrt-F145I exhibited elevated IC₉₀ values of piperaquine. Collectively, the emergence of these resistance patterns in Thailand near Cambodia border mirrored the reports of dihydroartemisinin-piperaquine treatment failures in the adjacent province of Cambodia, Oddar Meanchey, suggesting a migration of parasites across the border. As malaria elimination efforts ramp up in Southeast Asia, host nations militaries and other groups in border regions need to coordinate the proposed interventions.

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

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

Plasmodium falciparum Atg18 localizes to the food vacuole via interaction with the multi-drug resistance protein 1 and phosphatidylinositol 3-phosphate

Autophagy, a lysosome-dependent degradative process, does not appear to be a major degradative process in malaria parasites and has a limited repertoire of genes. To better understand the autophagy process, we investigated Plasmodium falciparum Atg18 (PfAtg18), a PROPPIN family protein, whose members like S. cerevisiae Atg18 (ScAtg18) and human WIPI2 bind PI3P and play an essential role in autophagosome formation. Wild type and mutant PfAtg18 were expressed in P. falciparum and assessed for localization, the effect of various inhibitors and antimalarials on PfAtg18 localization, and identification of PfAtg18-interacting proteins. PfAtg18 is expressed in asexual erythrocytic stages and localized to the food vacuole, which was also observed with other Plasmodium Atg18 proteins, indicating that food vacuole localization is likely a shared feature. Interaction of PfAtg18 with the food vacuole-associated PI3P is essential for localization, as PfAtg18 mutants of PI3P-binding motifs neither bound PI3P nor localized to the food vacuole. Interestingly, wild type ScAtg18 interacted with PI3P, but its expression in P. falciparum showed complete cytoplasmic localization, indicating additional requirement for food vacuole localization. The food vacuole multi-drug resistance protein 1 (MDR1) was consistently identified in the immunoprecipitates of PfAtg18 and P. berghei Atg18, and also interacted with PfAtg18. In contrast with PfAtg18, ScAtg18 did not interact with MDR1, which, in addition to PI3P, could play a critical role in localization of PfAtg18. Chloroquine and amodiaquine caused cytoplasmic localization of PfAtg18, suggesting that these target PfAtg18 transport pathway. Thus, PI3P and MDR1 are critical mediators of PfAtg18 localization.

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.

In Vitro Susceptibility of Plasmodium falciparum Isolates from the China-Myanmar Border Area to Piperaquine and Association with Candidate Markers

Plasmodium falciparum from the Greater Mekong subregion has evolved resistance to the artemisinin-based combination therapy dihydroartemisinin and the partner drug piperaquine. To monitor the potential westward spread or independent evolution of piperaquine resistance, we evaluated the in vitro susceptibility of 120 P. falciparum isolates collected at the China-Myanmar border during 2007-2016. The parasite isolates displayed a relatively wide range of piperaquine susceptibility estimates. While 56.7% of the parasites showed bimodal drug response curves, all but five generated area-under-the-curve (AUC) estimates consistent with a susceptible phenotype. Using the piperaquine survival assay (PSA), 5.6% parasites showed reduced susceptibility. Of note, parasites from 2014-2016 showed the highest AUC value and the highest proportion with a bimodal curve, suggesting falling effectiveness in these later years. Unsupervised K-mean analysis of the combined data assigned parasites into three clusters and identified significant correlations between IC₅₀, IC_90, and AUC values. No parasites carried the E415G mutation in a putative exo-nuclease, new mutations in PfCRT, or amplification of the plasmepsin 2/3 genes, suggesting mechanisms of reduced piperaquine susceptibility that differ from those described in other countries of the region. The association of increased AUC, IC₅₀, and IC₉₀ values with major PfK13 mutations (F446I and G533S) suggests that piperaquine resistance may evolve in these PfK13 genetic backgrounds. Additionally, the Pfmdr1 F1226Y mutation was associated with significantly higher PSA values. Further elucidation of piperaquine resistance mechanisms and continuous surveillance are warranted.

Multidrug-Resistant Plasmodium falciparum Parasites in the Central Highlands of Vietnam Jeopardize Malaria Control and Elimination Strategies

Plasmodium falciparum resistance to dihydroartemisinin-piperaquine has spread through the Greater Mekong Subregion to southwestern Vietnam. In 2018 to 2019, we collected 127 P. falciparum isolates from Dak Nong (36), Dak Lak (55), Gia Lai (13), and Kon Tum (23) provinces in Vietnam's Central Highlands and found parasites bearing the Pfkelch13 C580Y mutation and multiple plasmepsin 2/3 genes (mean prevalence, 17.9%; range, 4.3% to 27.8%), conferring resistance to dihydroartemisinin-piperaquine. This information is important for drug policy decisions in Vietnam.

Antimalarial Drug Resistance and Implications for the WHO Global Technical Strategy

PURPOSE OF REVIEW

Five years have passed since the World Health Organization released its Global Technical Strategy for Malaria (GTS). In that time, progress against malaria has plateaued. This review focuses on the implications of antimalarial drug resistance for the GTS and how interim progress in parasite genomics and antimalarial pharmacology offer a bulwark against it.

RECENT FINDINGS

For the first time, drug resistance-conferring genes have been identified and validated before their global expansion in malaria parasite populations. More efficient methods for their detection and elaboration have been developed, although low-density infections and polyclonality remain a nuisance to be solved. Clinical trials of alternative regimens for multidrug-resistant malaria have delivered promising results. New agents continue down the development pipeline, while a nascent infrastructure in sub-Saharan Africa for conducting phase I trials and trials of transmission-blocking agents has come to fruition after years of preparation.

SUMMARY

These and other developments can help inform the GTS as the world looks ahead to the next two decades of its implementation. To remain ahead of the threat that drug resistance poses, wider application of genomic-based surveillance and optimization of existing and forthcoming antimalarial drugs are essential.

Distribution and temporal dynamics of P. falciparum chloroquine resistance transporter mutations associated with piperaquine resistance in Northern Cambodia

BACKGROUND

Newly emerged mutations within the Plasmodium falciparum chloroquine resistance transporter (PfCRT) can confer piperaquine resistance in the absence of amplified plasmepsin II (pfpm2). In this study, we estimated the prevalence of co-circulating piperaquine resistance mutations in P. falciparum isolates collected in northern Cambodia from 2009-2017.

METHODS

The sequence of pfcrt was determined for 410 P. falciparum isolates using PacBio amplicon sequencing or whole genome sequencing. Quantitative PCR was used to estimate pfpm2 and pfmdr1 copy number.

RESULTS

Newly emerged PfCRT mutations increased in prevalence after the change to dihydroartemisinin-piperaquine in 2010, with >98% of parasites harboring these mutations by 2017. After 2014, the prevalence of PfCRT F145I declined, being out-competed by parasites with less resistant, but more fit PfCRT alleles. After the change to artesunate-mefloquine, the prevalence of parasites with amplified pfpm2 decreased, with nearly half of piperaquine-resistant PfCRT mutants having single copy pfpm2.

CONCLUSIONS

The large proportion of PfCRT mutants that lack pfpm2 amplification emphasizes the importance of including PfCRT mutations as part of molecular surveillance for piperaquine resistance in this region. Likewise, it is critical to monitor for amplified pfmdr1 in these PfCRT mutants, as increased mefloquine pressure could lead to mutants resistant to both drugs.

Transmission of Artemisinin-Resistant Malaria Parasites to Mosquitoes under Antimalarial Drug Pressure

Resistance to artemisinin-based combination therapy (ACT) in the Plasmodium falciparum parasite is threatening to reverse recent gains in reducing global deaths from malaria. While resistance manifests as delayed parasite clearance in patients, the phenotype can only spread geographically via the sexual stages and mosquito transmission. In addition to their asexual killing properties, artemisinin and its derivatives sterilize sexual male gametocytes. Whether resistant parasites overcome this sterilizing effect has not, however, been fully tested.

Resistance to artemisinin-based combination therapy (ACT) in the Plasmodium falciparum parasite is threatening to reverse recent gains in reducing global deaths from malaria. While resistance manifests as delayed parasite clearance in patients, the phenotype can only spread geographically via the sexual stages and mosquito transmission. In addition to their asexual killing properties, artemisinin and its derivatives sterilize sexual male gametocytes. Whether resistant parasites overcome this sterilizing effect has not, however, been fully tested. Here, we analyzed P. falciparum clinical isolates from the Greater Mekong Subregion, each demonstrating delayed clinical clearance and known resistance-associated polymorphisms in the Kelch13 (PfK13^var) gene. As well as demonstrating reduced asexual sensitivity to drug, certain PfK13^var isolates demonstrated a marked reduction in sensitivity to artemisinin in an in vitro male gamete formation assay. Importantly, this same reduction in sensitivity was observed when the most resistant isolate was tested directly in mosquito feeds. These results indicate that, under artemisinin drug pressure, while sensitive parasites are blocked, resistant parasites continue transmission. This selective advantage for resistance transmission could favor acquisition of additional host-specificity or polymorphisms affecting partner drug sensitivity in mixed infections. Favored resistance transmission under ACT coverage could have profound implications for the spread of multidrug-resistant malaria beyond Southeast Asia.

Modulation of Triple Artemisinin-Based Combination Therapy Pharmacodynamics by Plasmodium falciparum Genotype

The first-line treatments for uncomplicated Plasmodium falciparum malaria are artemisinin-based combination therapies (ACTs), consisting of an artemisinin derivative combined with a longer acting partner drug. However, the spread of P. falciparum with decreased susceptibility to artemisinin and partner drugs presents a significant challenge to malaria control efforts. To stem the spread of drug resistant parasites, novel chemotherapeutic strategies are being evaluated, including the implementation of triple artemisinin-based combination therapies (TACTs). Currently, there is limited knowledge on the pharmacodynamic and pharmacogenetic interactions of proposed TACT drug combinations. To evaluate these interactions, we established an in vitro high-throughput process for measuring the drug concentration-response to three distinct antimalarial drugs present in a TACT. Sixteen different TACT combinations were screened against 15 parasite lines from Cambodia, with a focus on parasites with differential susceptibilities to piperaquine and artemisinins. Analysis revealed drug-drug interactions unique to specific genetic backgrounds, including antagonism between piperaquine and pyronaridine associated with gene amplification of plasmepsin II/III, two aspartic proteases that localize to the parasite digestive vacuole. From this initial study, we identified parasite genotypes with decreased susceptibility to specific TACTs, as well as potential TACTs that display antagonism in a genotype-dependent manner. Our assay and analysis platform can be further leveraged to inform drug implementation decisions and evaluate next-generation TACTs.

Effect of mass dihydroartemisinin-piperaquine administration in southern Mozambique on the carriage of molecular markers of antimalarial resistance

BACKGROUND

Mass drug administration (MDA) can rapidly reduce the burden of Plasmodium falciparum (Pf). However, concerns remain about its contribution to select for antimalarial drug resistance.

METHODS

We used Sanger sequencing and real-time PCR to determine the proportion of molecular markers associated with antimalarial resistance (k13, pfpm2, pfmdr1 and pfcrt) in Pf isolates collected before (n = 99) and after (n = 112) the implementation of two monthly MDA rounds with dihydroartemisinin-piperaquine (DHAp) for two consecutive years in Magude district of Southern Mozambique.

RESULTS

None of the k13 polymorphisms associated with artemisinin resistance were observed in the Pf isolates analyzed. The proportion of Pf isolates with multiple copies of pfpm2, an amplification associated with piperaquine resistance, was similar in pre- (4.9%) and post-MDA groups (3.4%; p = 1.000). No statistically significant differences were observed between pre- and post-MDA groups in the proportion of Pf isolates neither with mutations in pfcrt and pfmdr1 genes, nor with the carriage of pfmdr1 multiple copies (p>0.05).

CONCLUSIONS

This study does not show any evidence of increased frequency of molecular makers of antimalarial resistance after MDA with DHAp in southern Mozambique where markers of antimalarial resistance were absent or low at the beginning of the intervention.

No evidence of amplified Plasmodium falciparum plasmepsin II gene copy number in an area with artemisinin-resistant malaria along the China-Myanmar border

Background

The emergence and spread of artemisinin resistance in Plasmodium falciparum poses a threat to malaria eradication, including China’s plan to eliminate malaria by 2020. Piperaquine (PPQ) resistance has emerged in Cambodia, compromising an important partner drug that is widely used in China in the form of dihydroartemisinin (DHA)-PPQ. Several mutations in a P. falciparum gene encoding a kelch protein on chromosome 13 (k13) are associated with artemisinin resistance and have arisen spread in the Great Mekong subregion, including the China–Myanmar border. Multiple copies of the plasmepsin II/III (pm2/3) genes, located on chromosome 14, have been shown to be associated with PPQ resistance.

Methods

The therapeutic efficacy of DHA-PPQ for the treatment of uncomplicated P. falciparum was evaluated along the China–Myanmar border from 2010 to 2014. The dry blood spots samples collected in the efficacy study prior DHA-PPQ treatment and from the local hospital by passive detection were used to amplify k13 and pm2. Polymorphisms within k13 were genotyped by capillary sequencing and pm2 copy number was quantified by relative-quantitative real-time polymerase chain reaction. Treatment outcome was evaluated with the World Health Organization protocol. A linear regression model was used to estimate the association between the day 3 positive rate and k13 mutation and the relationship of the pm2 copy number variants and k13 mutations.

Results

DHA-PPQ was effective for uncomplicated P. falciparum infection in Yunnan Province with cure rates > 95%. Twelve non synonymous mutations in the k13 domain were observed among the 268 samples with the prevalence of 44.0% and the predominant mutation was F446I with a prevalence of 32.8%. Only one sample was observed with multi-copies of pm2, including parasites with and without k13 mutations. The therapeutic efficacy of DHA-PPQ was > 95% along the China–Myanmar border, consistent with the lack of amplification of pm2.

Conclusion

DHA-PPQ for uncomplicated P. falciparum infection still showed efficacy in an area with artemisinin-resistant malaria along the China–Myanmar border. There was no evidence to show PPQ resistance by clinical study and molecular markers survey. Continued monitoring of the parasite population using molecular markers will be important to track emergence and spread of resistance in this region.

The natural function of the malaria parasite's chloroquine resistance transporter

The Plasmodium falciparum chloroquine resistance transporter (PfCRT) is a key contributor to multidrug resistance and is also essential for the survival of the malaria parasite, yet its natural function remains unresolved. We identify host-derived peptides of 4-11 residues, varying in both charge and composition, as the substrates of PfCRT in vitro and in situ, and show that PfCRT does not mediate the non-specific transport of other metabolites and/or ions. We find that drug-resistance-conferring mutations reduce both the peptide transport capacity and substrate range of PfCRT, explaining the impaired fitness of drug-resistant parasites. Our results indicate that PfCRT transports peptides from the lumen of the parasite's digestive vacuole to the cytosol, thereby providing a source of amino acids for parasite metabolism and preventing osmotic stress of this organelle. The resolution of PfCRT's native substrates will aid the development of drugs that target PfCRT and/or restore the efficacy of existing antimalarials.

Piperaquine resistant Cambodian Plasmodium falciparum clinical isolates: in vitro genotypic and phenotypic characterization

BACKGROUND

High rates of dihydroartemisinin-piperaquine (DHA-PPQ) treatment failures have been documented for uncomplicated Plasmodium falciparum in Cambodia. The genetic markers plasmepsin 2 (pfpm2), exonuclease (pfexo) and chloroquine resistance transporter (pfcrt) genes are associated with PPQ resistance and are used for monitoring the prevalence of drug resistance and guiding malaria drug treatment policy.

METHODS

To examine the relative contribution of each marker to PPQ resistance, in vitro culture and the PPQ survival assay were performed on seventeen P. falciparum isolates from northern Cambodia, and the presence of E415G-Exo and pfcrt mutations (T93S, H97Y, F145I, I218F, M343L, C350R, and G353V) as well as pfpm2 copy number polymorphisms were determined. Parasites were then cloned by limiting dilution and the cloned parasites were tested for drug susceptibility. Isobolographic analysis of several drug combinations for standard clones and newly cloned P. falciparum Cambodian isolates was also determined.

RESULTS

The characterization of culture-adapted isolates revealed that the presence of novel pfcrt mutations (T93S, H97Y, F145I, and I218F) with E415G-Exo mutation can confer PPQ-resistance, in the absence of pfpm2 amplification. In vitro testing of PPQ resistant parasites demonstrated a bimodal dose-response, the existence of a swollen digestive vacuole phenotype, and an increased susceptibility to quinine, chloroquine, mefloquine and lumefantrine. To further characterize drug sensitivity, parental parasites were cloned in which a clonal line, 14-B5, was identified as sensitive to artemisinin and piperaquine, but resistant to chloroquine. Assessment of the clone against a panel of drug combinations revealed antagonistic activity for six different drug combinations. However, mefloquine-proguanil and atovaquone-proguanil combinations revealed synergistic antimalarial activity.

CONCLUSIONS

Surveillance for PPQ resistance in regions relying on DHA-PPQ as the first-line treatment is dependent on the monitoring of molecular markers of drug resistance. P. falciparum harbouring novel pfcrt mutations with E415G-exo mutations displayed PPQ resistant phenotype. The presence of pfpm2 amplification was not required to render parasites PPQ resistant suggesting that the increase in pfpm2 copy number alone is not the sole modulator of PPQ resistance. Genetic background of circulating field isolates appear to play a role in drug susceptibility and biological responses induced by drug combinations. The use of latest field isolates may be necessary for assessment of relevant drug combinations against P. falciparum strains and when down-selecting novel drug candidates.

Prevalence of mutations in the Plasmodium falciparum chloroquine resistance transporter, PfCRT, and association with ex vivo susceptibility to common anti-malarial drugs against African Plasmodium falciparum isolates

Background

The Plasmodium falciparum chloroquine transporter gene (pfcrt) is known to be involved in chloroquine and amodiaquine resistance, and more particularly the mutations on the loci 72 to 76 localized within the second exon. Additionally, new mutations (T93S, H97Y, C101F, F145I, M343L, C350R and G353V) were recently shown to be associated with in vitro reduced susceptibility to piperaquine in Asian or South American P. falciparum strains. However, very few data are available on the prevalence of these mutations and their effect on parasite susceptibility to anti-malarial drugs, and more particularly piperaquine in Africa.

Methods

A molecular investigation of these mutations was performed in 602 African P. falciparum parasites collected between 2017 and 2018 on malaria patients hospitalized in France after a travel in African countries. Associations between genotypes and in vitro susceptibilities to piperaquine and standard antimalarial drugs were assessed.

Results

None of the mutations, previously described as associated with piperaquine resistance, was found in the 602 P. falciparum African isolates. The K76T mutation is associated with resistance to chloroquine (p < 0.0002) and desethylamodiaquine (p < 0.002) in Africa. The K76T mutation is not associated with in vitro reduced susceptibility to piperaquine. The mutation I356T, identified in 54.7% (n = 326) of the African isolates, was significantly associated with reduced susceptibility to quinine (p < 0.02) and increased susceptibility to mefloquine (p < 0.04). The K76T and I356T mutations were significantly associated in West African isolates (p = 0.008).

Conclusion

None of the mutations in pfcrt found to be associated with piperaquine reduced susceptibility in Asia or South America (T93S, H97Y, C101F, F145I, M343L C350R and G353V) were found in the 602 African isolates including the three isolates with reduced susceptibility to piperaquine. The K76T mutation, involved in resistance to chloroquine and amodiaquine, and the I356T mutation were not associated with in vitro reduced susceptibility to piperaquine. Differences in mefloquine susceptibility between I356 and 356T isolates were, while statistically different, minimal. Further analyses are needed with a more important sample size from the same geographic area to confirm the role of the I356T mutation on quinine susceptibility.

Development of copy number assays for detection and surveillance of piperaquine resistance associated plasmepsin 2/3 copy number variation in Plasmodium falciparum

BACKGROUND

Long regarded as an epicenter of drug-resistant malaria, Southeast Asia continues to provide new challenges to the control of Plasmodium falciparum malaria. Recently, resistance to the artemisinin combination therapy partner drug piperaquine has been observed in multiple locations across Southeast Asia. Genetic studies have identified single nucleotide polymorphisms as well as copy number variations in the plasmepsin 2 and plasmepsin 3 genes, which encode haemoglobin-degrading proteases that associate with clinical and in vitro piperaquine resistance.

RESULTS

To accurately and quickly determine the presence of copy number variations in the plasmepsin 2/3 genes in field isolates, this study developed a quantitative PCR assay using TaqMan probes. Copy number estimates were validated using a separate SYBR green-based quantitative PCR assay as well as a novel PCR-based breakpoint assay to detect the hybrid gene product. Field samples from 2012 to 2015 across three sites in Cambodia were tested using DNA extracted from dried blood spots and whole blood to monitor the extent of plasmepsin 2/3 gene amplifications, as well as amplifications in the multidrug resistance transporter 1 gene (pfmdr1), a marker of mefloquine resistance. This study found high concordance across all methods of copy number detection. For samples derived from dried blood spots, a success rate greater than 80% was found in each assay, with more recent samples performing better. Evidence of extensive plasmepsin 2/3 copy number amplifications was observed in Pursat (94%, 2015) (Western Cambodia) and Preah Vihear (87%, 2014) (Northern Cambodia), and lower levels in Ratanakiri (16%, 2014) (Eastern Cambodia). A shift was observed from two copies of plasmepsin 2 in Pursat in 2013 to three copies in 2014-2015 (25% to 64%). Pfmdr1 amplifications were absent in all samples from Preah Vihear and Ratanakiri in 2014 and absent in Pursat in 2015.

CONCLUSIONS

The multiplex TaqMan assay is a robust tool for monitoring both plasmepsin 2/3 and pfmdr1 copy number variations in field isolates, and the SYBR-green and breakpoint assays are useful for monitoring plasmepsin 2/3 amplifications. This study shows increasing levels of plasmepsin 2 copy numbers across Cambodia from 2012 to 2015 and a complete reversion of multicopy pfmdr1 parasites to single copy parasites in all study locations.