Multigenic architecture of piperaquine resistance trait in Plasmodium falciparum

Artemisinin resistance-associated gene mutations in Plasmodium falciparum: A case study of severe malaria from Mozambique

BACKGROUND

The effectiveness of artemisinin-based combination therapies (ACT) in treating Plasmodium falciparum, is vital for global malaria control efforts, particularly in sub-Saharan Africa. The examination of imported cases from endemic areas holds implications for malaria chemotherapy on a global scale.

METHOD

A 45-year-old male presented with high fever, dry cough, diarrhoea and generalized muscle pain, following a two-week trip to Mozambique. P. falciparum infection with hiperparasitemia was confirmed and the patient was treated initially with quinine and doxycycline, then intravenous artesunate. To assess drug susceptibility, ex vivo half-maximal inhibitory concentration assays were conducted, and the isolated P. falciparum genome was deep sequenced.

RESULTS

The clinical isolate exhibited elevated ex vivo half-maximal inhibitory concentration values to dihydroartemisinin, lumefantrine, mefloquine and piperaquine. Genomic analysis identified a I416V mutation in the P. falciparum Kelch13 (PF3D7_1343700) gene, and several mutations at the Kelch13 interaction candidate genes, pfkics (PF3D7_0813000, PF3D7_1138700, PF3D7_1246300), including the ubiquitin carboxyl-terminal hydrolase 1, pfubp1 (PF3D7_0104300). Mutations at the drug transporters and genes linked to next-generation antimalarial drug resistance were also present.

CONCLUSIONS

This case highlights the emergence of P. falciparum strains carrying mutations in artemisinin resistance-associated genes in Mozambique, couple with a reduction in ex vivo susceptibility to ACT drugs. Continuous surveillance of mutations linked to drug resistance and regular monitoring of drug susceptibility are imperative to anticipate the spread of potential resistant strains emerging in Mozambique and to maintain effective malaria control strategies.

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.

Plasmodium falciparum resistance to artemisinin-based combination therapies

Multidrug-resistant Plasmodium falciparum parasites are a major threat to public health in intertropical regions. Understanding the mechanistic basis, origins, and spread of resistance can inform strategies to mitigate its impact and reduce the global burden of malaria. The recent emergence in Africa of partial resistance to artemisinins, the core component of first-line combination therapies, is particularly concerning. Here, we review recent advances in elucidating the mechanistic basis of artemisinin resistance, driven primarily by point mutations in P. falciparum Kelch13, a key regulator of hemoglobin endocytosis and parasite response to artemisinin-induced stress. We also review resistance to partner drugs, including piperaquine and mefloquine, highlighting a key role for plasmepsins 2/3 and the drug and solute transporters P. falciparum chloroquine-resistance transporter and P. falciparum multidrug-resistance protein-1.

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.

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.

The newly discovered role of endocytosis in artemisinin resistance

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

Genomic Approaches to Drug Resistance in Malaria

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

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.

Expansion of a Specific Plasmodium falciparum PfMDR1 Haplotype in Southeast Asia with Increased Substrate Transport

Global efforts to eliminate malaria depend on the continued success of artemisinin-based combination therapies (ACTs) that target Plasmodium asexual blood-stage parasites. Resistance to ACTs, however, has emerged, creating the need to define the underlying mechanisms. Mutations in the P. falciparum multidrug resistance protein 1 (PfMDR1) transporter constitute an important determinant of resistance. Applying gene editing tools combined with an analysis of a public database containing thousands of parasite genomes, we show geographic selection and expansion of a pfmdr1 gene amplification encoding the N86/184F haplotype in Southeast Asia. Parasites expressing this PfMDR1 variant possess a higher transport capacity that modulates their responses to antimalarials. These data could help tailor and optimize antimalarial drug usage in different regions where malaria is endemic by taking into account the regional prevalence of pfmdr1 polymorphisms.

Artemisinin-based combination therapies (ACTs) have been vital in reducing malaria mortality rates since the 2000s. Their efficacy, however, is threatened by the emergence and spread of artemisinin resistance in Southeast Asia. The Plasmodium falciparum multidrug resistance protein 1 (PfMDR1) transporter plays a central role in parasite resistance to ACT partner drugs through gene copy number variations (CNV) and/or single nucleotide polymorphisms (SNPs). Using genomic epidemiology, we show that multiple pfmdr1 copies encoding the N86 and 184F haplotype are prevalent across Southeast Asia. Applying genome editing tools on the Southeast Asian Dd2 strain and using a surrogate assay to measure transporter activity in infected red blood cells, we demonstrate that parasites harboring multicopy N86/184F PfMDR1 have a higher Fluo-4 transport capacity compared with those expressing the wild-type N86/Y184 haplotype. Multicopy N86/184F PfMDR1 is also associated with decreased parasite susceptibility to lumefantrine. These findings provide evidence of the geographic selection and expansion of specific multicopy PfMDR1 haplotypes associated with multidrug resistance in Southeast Asia.

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.

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.

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.