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

Temporal and spatial dynamics of Plasmodium falciparum clonal lineages in Guyana

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

Lack of selection of antimalarial drug resistance markers after intermittent preventive treatment of schoolchildren (IPTsc) against malaria in northeastern Tanzania

OBJECTIVE

Intermittent Preventive Treatment of schoolchildren (IPTsc) is recommended by WHO as a strategy to protect against malaria; To explore whether IPTsc with dihydroartemisinin-piperaquine (DP) or artesunate-amodiaquine (ASAQ) cause a selection of molecular markers in Plasmodium falciparum genes associated with resistance in children in seven schools in Tanga region, Tanzania.

METHODS

SNPs in P. falciparum genes Pfmdr1, Pfexo, Pfkelch13 and Pfcrt and copy number variations in Pfplasmepsin-2 and Pfmdr1 were assessed in samples collected at 12 months (visit 4, n=74) and 20 months (visit 6, n=364) after initiation of IPTsc and compared with the baseline prevalence (n=379).

RESULTS

The prevalence of Pfmdr1 N86 and Pfexo 415G was >99% and 0%, respectively without any temporal differences observed. The prevalence of Pfmdr1 184F changed significantly from baseline (52.2%) to visit 6 (64.6%) (χ2=6.11, p=0.013), but no differences were observed between the treatment arms (χ2=0.05, p=0.98). Finally, only minor differences in the amplification of Pfmdr1 were observed; from 10.2% at baseline to 16.7% at visit 6 (χ2=0.98, p=0.32).

CONCLUSIONS

The IPTsc strategy does not seem to pose a risk for the selection of markers associated with DP or ASAQ resistance. Continuously and timely surveillance of markers of antimalarial drug resistance is recommended.

Antimalarial resistance risk in Mozambique detected by a novel quadruplex droplet digital PCR assay

While the Plasmodium falciparum malaria parasite continues to cause severe disease globally, Mozambique is disproportionally represented in malaria case totals. Acquisition of copy number variations (CNVs) in the parasite genome contributes to antimalarial drug resistance through overexpression of drug targets. Of interest, piperaquine resistance is associated with plasmepsin 2 and 3 CNVs (pfpmp2 and pfpmp3, respectively), while CNVs in the multidrug efflux pump, multidrug resistance-1 (pfmdr1), increase resistance to amodiaquine and lumefantrine. These antimalarials are partner drugs in artemisinin combination therapies (ACTs) and therefore, CNV detection with accurate and efficient tools is necessary to track ACT resistance risk. Here, we evaluated ~300 clinically derived samples collected from three sites in Mozambique for resistance-associated CNVs. We developed a novel, medium-throughput, quadruplex droplet digital PCR (ddPCR) assay to simultaneously quantify the copy number of pfpmp3, pfpmp2, and pfmdr1 loci in these clinical samples. By using DNA from laboratory parasite lines, we show that this nanodroplet-based method is capable of detecting picogram levels of parasite DNA, which facilitates its application for low yield and human host-contaminated clinical surveillance samples. Following ddPCR and the application of quality control standards, we detected CNVs in 13 of 229 high-quality samples (prevalence of 5.7%). Overall, our study revealed a low number of resistance CNVs present in the parasite population across all three collection sites, including various combinations of pfmdr1, pfpmp2, and pfpmp3 CNVs. The potential for future ACT resistance across Mozambique emphasizes the need for continued molecular surveillance across the region.

PHARE: a bioinformatics pipeline for compositional profiling of multiclonal Plasmodium falciparum infections from long-read Nanopore sequencing data

BACKGROUND

The emergence of drug-resistant clones of Plasmodium falciparum is a major public health concern, and the ability to detect and track the spread of these clones is crucial for effective malaria control and treatment. However, in endemic settings, malaria infected people often carry multiple P. falciparum clones simultaneously making it likely to miss drug-resistant clones using traditional molecular typing methods.

OBJECTIVES

Our goal was to develop a bioinformatics pipeline for compositional profiling in multiclonal P. falciparum samples, sequenced using the Oxford Nanopore Technologies MinION platform.

METHODS

We developed the 'Finding P. falciparum haplotypes with resistance mutations in polyclonal infections' (PHARE) pipeline using existing bioinformatics tools and custom scripts written in python. PHARE was validated on three control datasets containing P. falciparum DNA of four laboratory strains at varying mixing ratios. Additionally, the pipeline was tested on clinical samples from children admitted to a paediatric hospital in the Central African Republic.

RESULTS

The PHARE pipeline achieved high recall and accuracy rates in all control datasets. The pipeline can be used on any gene and was tested with amplicons of the P. falciparum drug resistance marker genes pfdhps, pfdhfr and pfK13.

CONCLUSIONS

The PHARE pipeline helps to provide a more complete picture of drug resistance in the circulating P. falciparum population and can help to guide treatment recommendations. PHARE is freely available under the GNU Lesser General Public License v.3.0 on GitHub: https://github.com/Fippu/PHARE.

Temporal and spatial dynamics of Plasmodium falciparum clonal lineages in Guyana

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

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.

Targeted and whole-genome sequencing reveal a north-south divide in P. falciparum drug resistance markers and genetic structure in Mozambique

Mozambique is one of the four African countries which account for over half of all malaria deaths worldwide, yet little is known about the parasite genetic structure in that country. We performed P. falciparum amplicon and whole genome sequencing on 2251 malaria-infected blood samples collected in 2015 and 2018 in seven provinces of Mozambique to genotype antimalarial resistance markers and interrogate parasite population structure using genome-wide microhaplotyes. Here we show that the only resistance-associated markers observed at frequencies above 5% were pfmdr1-184F (59%), pfdhfr-51I/59 R/108 N (99%) and pfdhps-437G/540E (89%). The frequency of pfdhfr/pfdhps quintuple mutants associated with sulfadoxine-pyrimethamine resistance increased from 80% in 2015 to 89% in 2018 (p < 0.001), with a lower expected heterozygosity and higher relatedness of microhaplotypes surrounding pfdhps mutants than wild-type parasites suggestive of recent selection. pfdhfr/pfdhps quintuple mutants also increased from 72% in the north to 95% in the south (2018; p < 0.001). This resistance gradient was accompanied by a concentration of mutations at pfdhps-436 (17%) in the north, a south-to-north increase in the genetic complexity of P. falciparum infections (p = 0.001) and a microhaplotype signature of regional differentiation. The parasite population structure identified here offers insights to guide antimalarial interventions and epidemiological surveys.

Malaria outbreak in Laos driven by a selective sweep for Plasmodium falciparum kelch13 R539T mutants: a genetic epidemiology analysis

BACKGROUND

Malaria outbreaks are important public health concerns that can cause resurgence in endemic regions approaching elimination. We investigated a Plasmodium falciparum outbreak in Attapeu Province, Laos, during the 2020-21 malaria season, using genomic epidemiology methods to elucidate parasite population dynamics and identify its causes.

METHODS

In this genetic analysis, 2164 P falciparum dried blood spot samples were collected from southern Laos between Jan 1, 2017, and April 1, 2021, which included 249 collected during the Attapeu outbreak between April 1, 2020, and April 1, 2021, by routine surveillance. Genetic barcodes obtained from these samples were used to investigate epidemiological changes underpinning the outbreak, estimate population diversity, and analyse population structure. Whole-genome sequencing data from additional historical samples were used to reconstruct the ancestry of outbreak strains using identity-by-descent analyses.

FINDINGS

The outbreak parasite populations were characterised by unprecedented loss of genetic diversity, primarily caused by rapid clonal expansion of a multidrug-resistant strain (LAA1) carrying the kelch13 Arg539Thr (R539T) mutation. LAA1 replaced kelch13 Cys580Tyr (C580Y) mutants resistant to dihydroartemisinin-piperaquine (KEL1/PLA1) as the dominant strain. LAA1 inherited 58·8% of its genome from a strain circulating in Cambodia in 2008. A secondary outbreak strain (LAA2) carried the kelch13 C580Y allele, and a genome that is essentially identical to a Cambodian parasite from 2009. A third, low-frequency strain (LAA7) was a recombinant of KEL1/PLA1 with a kelch13 R539T mutant.

INTERPRETATION

These results strongly suggest that the outbreak was driven by a selective sweep, possibly associated with multidrug-resistant phenotypes of the outbreak strains. Established resistant populations can circulate at low frequencies for years before suddenly overwhelming dominant strains when the conditions for selection become favourable-eg, when front-line therapies change. Genetic surveillance can support elimination by characterising key properties of outbreaks such as population diversity, drug resistance marker prevalence, and the origins of outbreak strains.

FUNDING

Bill & Melinda Gates Foundation; The Global Fund to Fight AIDS, Tuberculosis and Malaria; Wellcome Trust.

TRANSLATION

For the Lao translation of the abstract see Supplementary Materials section.

Impact of parasite genomic dynamics on the sensitivity of Plasmodium falciparum isolates to piperaquine and other antimalarial drugs

BACKGROUND

Dihydroartemisinin-piperaquine (DHA-PPQ) is an alternative first-line antimalarial to artemether-lumefantrine in Kenya. However, recent reports on the emergence of PPQ resistance in Southeast Asia threaten its continued use in Kenya and Africa. In line with the policy on continued deployment of DHA-PPQ, it is imperative to monitor the susceptibility of Kenyan parasites to PPQ and other antimalarials.

METHODS

Parasite isolates collected between 2008 and 2021 from individuals with naturally acquired P. falciparum infections presenting with uncomplicated malaria were tested for in vitro susceptibility to piperaquine, dihydroartemisinin, lumefantrine, artemether, and chloroquine using the malaria SYBR Green I method. A subset of the 2019-2021 samples was further tested for ex vivo susceptibility to PPQ using piperaquine survival assay (PSA). Each isolate was also characterized for mutations associated with antimalarial resistance in Pfcrt, Pfmdr1, Pfpm2/3, Pfdhfr, and Pfdhps genes using real-time PCR and Agena MassARRAY platform. Associations between phenotype and genotype were also determined.

RESULTS

The PPQ median IC₅₀ interquartile range (IQR) remained stable during the study period, 32.70 nM (IQR 20.2-45.6) in 2008 and 27.30 nM (IQR 6.9-52.8) in 2021 (P=0.1615). The median ex vivo piperaquine survival rate (IQR) was 0% (0-5.27) at 95% CI. Five isolates had a PSA survival rate of ≥10%, consistent with the range of PPQ-resistant parasites, though they lacked polymorphisms in Pfmdr1 and Plasmepsin genes. Lumefantrine and artemether median IC₅₀s rose significantly to 62.40 nM (IQR 26.9-100.8) (P = 0.0201); 7.00 nM (IQR 2.4-13.4) (P = 0.0021) in 2021 from 26.30 nM (IQR 5.1-64.3); and 2.70 nM (IQR 1.3-10.4) in 2008, respectively. Conversely, chloroquine median IC₅₀s decreased significantly to 10.30 nM (IQR 7.2-20.9) in 2021 from 15.30 nM (IQR 7.6-30.4) in 2008, coinciding with a decline in the prevalence of Pfcrt 76T allele over time (P = 0.0357). The proportions of piperaquine-resistant markers including Pfpm2/3 and Pfmdr1 did not vary significantly. A significant association was observed between PPQ IC₅₀ and Pfcrt K76T allele (P=0.0026).

CONCLUSIONS

Circulating Kenyan parasites have remained sensitive to PPQ and other antimalarials, though the response to artemether (ART) and lumefantrine (LM) is declining. This study forms a baseline for continued surveillance of current antimalarials for timely detection of resistance.

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.

Artesunate and dihydroartemisinin-piperaquine treatment failure in a severe Plasmodium falciparum malaria case imported from Republic of Côte d'Ivoire

A 68-year-old man returning from Republic of Côte d'Ivoire (Ivory Coast) was diagnosed with severe Plasmodium falciparum malaria and treated with intravenous artesunate followed by oral dihydroartemisinin-piperaquine (DHA-PPQ). A month later the patient experienced a new P. falciparum episode; analysis of pfmsp-1 and pfmsp-2 revealed that the infection was caused by a genetic strain identical to the strain that caused the initial episode, indicating resurgence of the previous infection. No mutations in genes associated with resistance to artemisinin derivatives (pfk13) or piperaquine (pfexonuclease, pfplasmepsin 2/3) were detected, suggesting that treatment failure could have been caused by drug malabsorption or poor drug manufacturing practices. A second treatment with atovaquone-proguanil was successful in eliminating the infection, with no further relapses. To our knowledge, this is the first description of a treatment failure with both artesunate and DHA-PPQ in a traveler returning from a malaria-endemic region. Analysis of molecular markers of resistance to antimalarial drugs revealed mutations associated with resistance to sulfadoxine (pfdhps) and pyrimethamine (pfdhfr), highlighting the important contribution of surveillance of imported malaria cases to the monitoring of drug resistance globally.

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.

Molecular surveillance for operationally relevant genetic polymorphisms in Plasmodium falciparum in Southern Chad, 2016–2017

Background

Resistance to anti-malarials is a serious threat to the efforts to control and eliminate malaria. Surveillance based on simple field protocols with centralized testing to detect molecular markers associated with anti-malarial drug resistance can be used to identify locations where further investigations are needed.

Methods

Dried blood spots were collected from 398 patients (age range 5–59 years, 99% male) with Plasmodium falciparum infections detected using rapid diagnostic tests over two rounds of sample collection conducted in 2016 and 2017 in Komé, South-West Chad. Specimens were genotyped using amplicon sequencing or qPCR for validated markers of anti-malarial resistance including partner drugs used in artemisinin-based combination therapy (ACT).

Results

No mutations in the pfk13 gene known to be associated with artemisinin resistance were found but a high proportion of parasites carried other mutations, specifically K189T (190/349, 54.4%, 95%CI 49.0–59.8%). Of 331 specimens successfully genotyped for pfmdr1 and pfcrt, 52% (95%CI 46.4–57.5%) carried the NFD-K haplotype, known to be associated with reduced susceptibility to lumefantrine. Only 20 of 336 (6.0%, 95%CI 3.7–9.0%) had parasites with the pfmdr1-N86Y polymorphism associated with increased treatment failures with amodiaquine. Nearly all parasites carried at least one mutation in pfdhfr and/or pfdhps genes but ‘sextuple’ mutations in pfdhfr—pfdhps including pfdhps -A581G were rare (8/336 overall, 2.4%, 95%CI 1.2–4.6%). Only one specimen containing parasites with pfmdr1 gene amplification was detected.

Conclusions

These results provide information on the likely high efficacy of artemisinin-based combinations commonly used in Chad, but suggest decreasing levels of sensitivity to lumefantrine and high levels of resistance to sulfadoxine-pyrimethamine used for seasonal malaria chemoprevention and intermittent preventive therapy in pregnancy. A majority of parasites had mutations in the pfk13 gene, none of which are known to be associated with artemisinin resistance. A therapeutic efficacy study needs to be conducted to confirm the efficacy of artemether-lumefantrine.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12936-022-04095-9.

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.

Boosting the impact of seasonal malaria chemoprevention (SMC) through simultaneous screening and treatment of household members of children receiving SMC in Burkina Faso: a protocol for a randomized open label trial

BACKGROUND

Plasmodium falciparum malaria remains a major public health concern in sub-Sahara Africa. Seasonal malaria chemoprevention (SMC) with amodiaquine + sulfadoxine-pyrimethamine is one of the most important preventive interventions. Despite its implementation, the burden of malaria is still very high in children under five years old in Burkina Faso, suggesting that the expected impact of this promising strategy might not be attained. Development of innovative strategies to improve the efficacy of these existing malaria control measures is essential. In such context, we postulate that screening and treatment of malaria in household members of children receiving SMC could greatly improve the impact of SMC intervention and reduce malaria transmission in endemic settings.

METHODS

This randomized superiority trial will be carried out in the Nanoro health district, Burkina Faso. The unit of randomisation will be the household and all eligible children from a household will be allocated to the same study group. Households with 3-59 months old children will be assigned to either (i) control group (SMC alone) or (ii) intervention (SMC+ screening of household members with standard Histidin Rich Protein Rapid Diagnostic Test (HRP2-RDT) and treatment if positive). The sample size will be 526 isolated households per arm, i.e., around 1052 children under SMC coverage and an expected 1315 household members. Included children will be followed-up for 24 months to fully cover two consecutive malaria transmission seasons and two SMC cycles. Children will be actively followed-up during the malaria transmission seasons while in the dry seasons the follow-up will be passive.

CONCLUSION

The study will respond to a major public health concern by providing evidence of the efficacy of an innovative strategy to boost the impact of SMC intervention.

Evolution of multidrug resistance in Plasmodium falciparum: a longitudinal study of genetic resistance markers in the Greater Mekong Subregion

Increasing resistance in Plasmodium falciparum to artemisinins and their artemisinin combination therapy (ACT) partner drugs jeopardizes effective antimalarial treatment. Resistance is worst in the Greater Mekong subregion. Monitoring genetic markers of resistance can help to guide antimalarial therapy. Markers of resistance to artemisinins (PfKelch mutations), mefloquine (amplification of P. falciparum multidrug resistance-1 [PfMDR1]), and piperaquine (PfPlasmepsin2/3 amplification and specific P. falciparum chloroquine resistance transporter [PfCRT] mutations) were assessed in 6,722 P. falciparum samples from Vietnam, Lao People’s Democratic Republic (PDR), Cambodia, Thailand, and Myanmar between 2007 and 2019. Against a high background prevalence of PfKelch mutations, PfMDR1 and PfPlasmepsin2/3 amplification closely followed regional drug pressures over time. PfPlasmepsin2/3 amplification preceded piperaquine resistance-associated PfCRT mutations in Cambodia and reached a peak prevalence of 23/28 (82%) in 2015. This declined to 57/156 (38%) after first-line treatment was changed from dihydroartemisinin-piperaquine to artesunate-mefloquine (ASMQ) between 2014 and 2017. The frequency of PfMDR1 amplification increased from 0/293 (0%) between 2012 and 2017 to 12/156 (8%) in 2019. Amplification of PfMDR1 and PfPlasmepsin2/3 in the same parasites was extremely rare (4/6,722 [0.06%]) and was dispersed over time. The mechanisms conferring mefloquine and piperaquine resistance may be counterbalancing. This supports the development of ASMQ plus piperaquine as a triple artemisinin combination therapy.

Advances and opportunities in malaria population genomics

Almost 20 years have passed since the first reference genome assemblies were published for Plasmodium falciparum, the deadliest malaria parasite, and Anopheles gambiae, the most important mosquito vector of malaria in sub-Saharan Africa. Reference genomes now exist for all human malaria parasites and nearly half of the ~40 important vectors around the world. As a foundation for genetic diversity studies, these reference genomes have helped advance our understanding of basic disease biology and drug and insecticide resistance, and have informed vaccine development efforts. Population genomic data are increasingly being used to guide our understanding of malaria epidemiology, for example by assessing connectivity between populations and the efficacy of parasite and vector interventions. The potential value of these applications to malaria control strategies, together with the increasing diversity of genomic data types and contexts in which data are being generated, raise both opportunities and challenges in the field. This Review discusses advances in malaria genomics and explores how population genomic data could be harnessed to further support global disease control efforts.

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.