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.

An open dataset of Plasmodium vivax genome variation in 1,895 worldwide samples.

This report describes the MalariaGEN Pv4 dataset, a new release of curated genome variation data on 1,895 samples of Plasmodium vivax collected at 88 worldwide locations between 2001 and 2017. It includes 1,370 new samples contributed by MalariaGEN and VivaxGEN partner studies in addition to previously published samples from these and other sources. We provide genotype calls at over 4.5 million variable positions including over 3 million single nucleotide polymorphisms (SNPs), as well as short indels and tandem duplications. This enlarged dataset highlights major compartments of parasite population structure, with clear differentiation between Africa, Latin America, Oceania, Western Asia and different parts of Southeast Asia. Each sample has been classified for drug resistance to sulfadoxine, pyrimethamine and mefloquine based on known markers at the dhfr, dhps and mdr1 loci. The prevalence of all of these resistance markers was much higher in Southeast Asia and Oceania than elsewhere. This open resource of analysis-ready genome variation data from the MalariaGEN and VivaxGEN networks is driven by our collective goal to advance research into the complex biology of P. vivax and to accelerate genomic surveillance for malaria control and elimination.

Comparison of Antibody Responses and Parasite Clearance in Artemisinin Therapeutic Efficacy Studies in the Democratic Republic of Congo and Asia.

BACKGROUND: Understanding the effect of immunity on Plasmodium falciparum clearance is essential for interpreting therapeutic efficacy studies designed to monitor emergence of artemisinin drug resistance. In low-transmission areas of Southeast Asia, where resistance has emerged, P. falciparum antibodies confound parasite clearance measures. However, variation in naturally acquired antibodies across Asian and sub-Saharan African epidemiological contexts and their impact on parasite clearance re yet to be quantified. METHODS: In an artemisinin therapeutic efficacy study, antibodies to 12 pre-erythrocytic and erythrocytic P. falciparum antigens were measured in 118 children with uncomplicated P. falciparum malaria in the Democratic Republic of Congo (DRC) and compared with responses in patients from Asian sites, described elsewhere. RESULTS: Parasite clearance half-life was shorter in DRC patients (median, 2 hours) compared with most Asian sites (median, 2-7 hours), but P. falciparum antibody levels and seroprevalences were similar. There was no evidence for an association between antibody seropositivity and parasite clearance half-life (mean difference between seronegative and seropositive, -0.14 to +0.40 hour) in DRC patients. CONCLUSIONS: In DRC, where artemisinin remains highly effective, the substantially shorter parasite clearance time compared with Asia was not explained by differences in the P. falciparum antibody responses studied.

Anti-Gametocyte Antigen Humoral Immunity and Gametocytemia During Treatment of Uncomplicated Falciparum Malaria: A Multi-National Study.

Introduction: Understanding the human immune response to Plasmodium falciparum gametocytes and its association with gametocytemia is essential for understanding the transmission of malaria as well as progressing transmission blocking vaccine candidates. Methods: In a multi-national clinical efficacy trial of artemisinin therapies (13 sites of varying transmission over South-East Asia and Africa), we measured Immunoglobulin G (IgG) responses to recombinant P. falciparum gametocyte antigens expressed on the gametocyte plasma membrane and leading transmission blocking vaccine candidates Pfs230 (Pfs230c and Pfs230D1M) and Pfs48/45 at enrolment in 1,114 participants with clinical falciparum malaria. Mixed effects linear and logistic regression were used to determine the association between gametocyte measures (gametocytemia and gametocyte density) and antibody outcomes at enrolment. Results: Microscopy detectable gametocytemia was observed in 11% (127/1,114) of participants at enrolment, and an additional 9% (95/1,114) over the follow-up period (up to day 42) (total 20% of participants [222/1,114]). IgG levels in response to Pfs230c, Pfs48/45 and Pfs230D1M varied across study sites at enrolment (p < 0.001), as did IgG seroprevalence for anti-Pfs230c and D1M IgG (p < 0.001), but not for anti-Pfs48/45 IgG (p = 0.159). In adjusted analyses, microscopy detectable gametocytemia at enrolment was associated with an increase in the odds of IgG seropositivity to the three gametocyte antigens (Pfs230c OR [95% CI], p: 1.70 [1.10, 2.62], 0.017; Pfs48/45: 1.45 [0.85, 2.46], 0.174; Pfs230D1M: 1.70 [1.03, 2.80], 0.037), as was higher gametocyte density at enrolment (per two-fold change in gametocyte density Pfs230c OR [95% CI], p: 1.09 [1.02, 1.17], 0.008; Pfs48/45: 1.05 [0.98, 1.13], 0.185; Pfs230D1M: 1.07 [0.99, 1.14], 0.071). Conclusion: Pfs230 and Pfs48/45 antibodies are naturally immunogenic targets associated with patent gametocytemia and increasing gametocyte density across multiple malaria endemic settings, including regions with emerging artemisinin-resistant P. falciparum.

Impact of Sickle Cell Trait Hemoglobin on the Intraerythrocytic Transcriptional Program of Plasmodium falciparum.

Sickle-trait hemoglobin (HbAS) confers nearly complete protection from severe, life-threatening falciparum malaria in African children. Despite this clear protection, the molecular mechanisms by which HbAS confers these protective phenotypes remain incompletely understood. As a forward genetic screen for aberrant parasite transcriptional responses associated with parasite neutralization in HbAS red blood cells (RBCs), we performed comparative transcriptomic analyses of Plasmodium falciparum in normal (HbAA) and HbAS erythrocytes during both in vitro cultivation of reference parasite strains and naturally occurring P. falciparum infections in Malian children with HbAA or HbAS. During in vitro cultivation, parasites matured normally in HbAS RBCs, and the temporal expression was largely unperturbed of the highly ordered transcriptional program that underlies the parasite's maturation throughout the intraerythrocytic development cycle (IDC). However, differential expression analysis identified hundreds of transcripts aberrantly expressed in HbAS, largely occurring late in the IDC. Surprisingly, transcripts encoding members of the Maurer's clefts were overexpressed in HbAS despite impaired parasite protein export in these RBCs, while parasites in HbAS RBCs underexpressed transcripts associated with the endoplasmic reticulum and those encoding serine repeat antigen proteases that promote parasite egress. Analyses of P. falciparum transcriptomes from 32 children with uncomplicated malaria identified stage-specific differential expression: among infections composed of ring-stage parasites, only cyclophilin 19B was underexpressed in children with HbAS, while trophozoite-stage infections identified a range of differentially expressed transcripts, including downregulation in HbAS of several transcripts associated with severe malaria in collateral studies. Collectively, our comparative transcriptomic screen in vitro and in vivo indicates that P. falciparum adapts to HbAS by altering its protein chaperone and folding machinery, oxidative stress response, and protein export machinery. Because HbAS consistently protects from severe P. falciparum, modulation of these responses may offer avenues by which to neutralize P. falciparum parasites. IMPORTANCE Sickle-trait hemoglobin (HbAS) confers nearly complete protection from severe, life-threatening malaria, yet the molecular mechanisms that underlie HbAS protection from severe malaria remain incompletely understood. Here, we used transcriptome sequencing (RNA-seq) to measure the impact of HbAS on the blood-stage transcriptome of Plasmodium falciparum in in vitro time series experiments and in vivo samples from natural infections. Our in vitro time series data reveal that, during its blood stage, P. falciparum's gene expression in HbAS is impacted primarily through alterations in the abundance of gene products as opposed to variations in the timing of gene expression. Collectively, our in vitro and in vivo data indicate that P. falciparum adapts to HbAS by altering its protein chaperone and folding machinery, oxidative stress response, and protein export machinery. Due to the persistent association of HbAS and protection from severe disease, these processes that are modified in HbAS may offer strategies to neutralize P. falciparum.

Development of weight and age-based dosing of daily primaquine for radical cure of vivax malaria.

BACKGROUND: In many endemic areas, Plasmodium vivax malaria is predominantly a disease of young adults and children. International recommendations for radical cure recommend fixed target doses of 0.25 or 0.5 mg/kg/day of primaquine for 14 days in glucose-6-phosphate dehydrogenase normal patients of all ages. However, for many anti-malarial drugs, including primaquine, there is evidence that children have lower exposures than adults for the same weight-adjusted dose. The aim of the study was to develop 14-day weight-based and age-based primaquine regimens against high-frequency relapsing tropical P. vivax. METHODS: The recommended adult target dose of 0.5 mg/kg/day (30 mg in a 60 kg patient) is highly efficacious against tropical P. vivax and was assumed to produce optimal drug exposure. Primaquine doses were calculated using allometric scaling to derive a weight-based primaquine regimen over a weight range from 5 to 100 kg. Growth curves were constructed from an anthropometric database of 53,467 individuals from the Greater Mekong Subregion (GMS) to define weight-for-age relationships. The median age associated with each weight was used to derive an age-based dosing regimen from the weight-based regimen. RESULTS: The proposed weight-based regimen has 5 dosing bands: (i) 5-7 kg, 5 mg, resulting in 0.71-1.0 mg/kg/day; (ii) 8-16 kg, 7.5 mg, 0.47-0.94 mg/kg/day; (iii) 17-40 kg, 15 mg, 0.38-0.88 mg/kg/day; (iv) 41-80 kg, 30 mg, 0.37-0.73 mg/kg/day; and (v) 81-100 kg, 45 mg, 0.45-0.56 mg/kg/day. The corresponding age-based regimen had 4 dosing bands: 6-11 months, 5 mg, 0.43-1.0 mg/kg/day; (ii) 1-5 years, 7.5 mg, 0.35-1.25 mg/kg/day; (iii) 6-14 years, 15 mg, 0.30-1.36 mg/kg/day; and (iv) ≥ 15 years, 30 mg, 0.35-1.07 mg/kg/day. CONCLUSION: The proposed weight-based regimen showed less variability around the primaquine dose within each dosing band compared to the age-based regimen and is preferred. Increased dose accuracy could be achieved by additional dosing bands for both regimens. The age-based regimen might not be applicable to regions outside the GMS, which must be based on local anthropometric data. Pharmacokinetic data in small children are needed urgently to inform the proposed regimens.

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.

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.

Sickle-trait hemoglobin reduces adhesion to both CD36 and EPCR by Plasmodium falciparum-infected erythrocytes.

Sickle-trait hemoglobin protects against severe Plasmodium falciparum malaria. Severe malaria is governed in part by the expression of the Plasmodium falciparum Erythrocyte Membrane Protein 1 (PfEMP1) that are encoded by var genes, specifically those variants that bind Endothelial Protein C Receptor (EPCR). In this study, we investigate the effect of sickle-trait on parasite var gene expression and function in vitro and in field-collected parasites. We mapped var gene reads generated from RNA sequencing in parasite cultures in normal and sickle-cell trait blood throughout the asexual lifecycle. We investigated sickle-trait effect on PfEMP1 interactions with host receptors CD36 and EPCR using static adhesion assays and flow cytometry. Var expression in vivo was compared by assembling var domains sequenced from total RNA in parasites infecting Malian children with HbAA and HbAS. Sickle-trait did not alter the abundance or type of var gene transcripts in vitro, nor the abundance of overall transcripts or of var functional domains in vivo. In adhesion assays using recombinant host receptors, sickle-trait reduced adhesion by 73-86% to CD36 and 83% to EPCR. Similarly, sickle-trait reduced the surface expression of EPCR-binding PfEMP1. In conclusion, Sickle-cell trait does not directly affect var gene transcription but does reduce the surface expression and function of PfEMP1. This provides a direct mechanism for protection against severe malaria conferred by sickle-trait hemoglobin. Trial Registration: ClinicalTrials.gov Identifier: NCT02645604.

Artemisinin-resistant K13 mutations rewire Plasmodium falciparum's intra-erythrocytic metabolic program to enhance survival.

The emergence and spread of artemisinin resistance, driven by mutations in Plasmodium falciparum K13, has compromised antimalarial efficacy and threatens the global malaria elimination campaign. By applying systems-based quantitative transcriptomics, proteomics, and metabolomics to a panel of isogenic K13 mutant or wild-type P. falciparum lines, we provide evidence that K13 mutations alter multiple aspects of the parasite's intra-erythrocytic developmental program. These changes impact cell-cycle periodicity, the unfolded protein response, protein degradation, vesicular trafficking, and mitochondrial metabolism. K13-mediated artemisinin resistance in the Cambodian Cam3.II line was reversed by atovaquone, a mitochondrial electron transport chain inhibitor. These results suggest that mitochondrial processes including damage sensing and anti-oxidant properties might augment the ability of mutant K13 to protect P. falciparum against artemisinin action by helping these parasites undergo temporary quiescence and accelerated growth recovery post drug elimination.

State of Artemisinin and Partner Drug Susceptibility in Plasmodium falciparum Clinical Isolates from Colombia.

Delayed parasite clearance time observed in Southeast Asia provided the first evidence of Plasmodium falciparum resistance to artemisinins. The ex vivo ring-stage survival assay (RSA) mimics parasite exposure to pharmacologically relevant artemisinin concentrations. Mutations in the C-terminal propeller domain of the putative kelch protein Pf3D7_1343700 (K13) are associated with artemisinin resistance. Variations in the pfmdr1 gene are associated with reduced susceptibility to the artemisinin partner drugs mefloquine (MQ) and lumefantrine (LF). To clarify the unknown landscape of artemisinin resistance in Colombia, 71 patients with uncomplicated P. falciparum malaria were enrolled in a non-randomized observational study in three endemic localities in 2014-2015. Each patient's parasite isolate was assessed for ex vivo RSA, K13-propeller mutations, pfmdr1 copy number, and pfmdr1 mutations at codons 86, 184, 1034, 1042, and 1246, associated with reduced susceptibility, and 50% inhibitory concentration (IC50) for other antimalarial drugs. Ex vivo RSAs were successful in 56% (40/71) of samples, and nine isolates showed survival rates > 1%. All isolates had wild-type K13-propeller sequences. All isolates harbored either of two pfmdr1 haplotypes, NFSDD (79.3%) and NFSDY (20.7%), and 7.1% of isolates had > 1 pfmdr1 gene. In vitro IC50 assays showed that variable proportions of isolates had decreased susceptibility to chloroquine (52.4%, > 100 nM), amodiaquine (31.2%, > 30 nM), MQ (34.3%, > 30 nM), and LF (3.2%, > 10 nM). In this study, we report ex vivo RSA and K13 data on P. falciparum isolates from Colombia. The identification of isolates with increased ex vivo RSA rates in the absence of K13-propeller mutations and no positivity at day three requires further investigation.

The risk of Plasmodium vivax parasitaemia after P. falciparum malaria: An individual patient data meta-analysis from the WorldWide Antimalarial Resistance Network.

BACKGROUND: There is a high risk of Plasmodium vivax parasitaemia following treatment of falciparum malaria. Our study aimed to quantify this risk and the associated determinants using an individual patient data meta-analysis in order to identify populations in which a policy of universal radical cure, combining artemisinin-based combination therapy (ACT) with a hypnozoitocidal antimalarial drug, would be beneficial. METHODS AND FINDINGS: A systematic review of Medline, Embase, Web of Science, and the Cochrane Database of Systematic Reviews identified efficacy studies of uncomplicated falciparum malaria treated with ACT that were undertaken in regions coendemic for P. vivax between 1 January 1960 and 5 January 2018. Data from eligible studies were pooled using standardised methodology. The risk of P. vivax parasitaemia at days 42 and 63 and associated risk factors were investigated by multivariable Cox regression analyses. Study quality was assessed using a tool developed by the Joanna Briggs Institute. The study was registered in the International Prospective Register of Systematic Reviews (PROSPERO: CRD42018097400). In total, 42 studies enrolling 15,341 patients were included in the analysis, including 30 randomised controlled trials and 12 cohort studies. Overall, 14,146 (92.2%) patients had P. falciparum monoinfection and 1,195 (7.8%) mixed infection with P. falciparum and P. vivax. The median age was 17.0 years (interquartile range [IQR] = 9.0-29.0 years; range = 0-80 years), with 1,584 (10.3%) patients younger than 5 years. 2,711 (17.7%) patients were treated with artemether-lumefantrine (AL, 13 studies), 651 (4.2%) with artesunate-amodiaquine (AA, 6 studies), 7,340 (47.8%) with artesunate-mefloquine (AM, 25 studies), and 4,639 (30.2%) with dihydroartemisinin-piperaquine (DP, 16 studies). 14,537 patients (94.8%) were enrolled from the Asia-Pacific region, 684 (4.5%) from the Americas, and 120 (0.8%) from Africa. At day 42, the cumulative risk of vivax parasitaemia following treatment of P. falciparum was 31.1% (95% CI 28.9-33.4) after AL, 14.1% (95% CI 10.8-18.3) after AA, 7.4% (95% CI 6.7-8.1) after AM, and 4.5% (95% CI 3.9-5.3) after DP. By day 63, the risks had risen to 39.9% (95% CI 36.6-43.3), 42.4% (95% CI 34.7-51.2), 22.8% (95% CI 21.2-24.4), and 12.8% (95% CI 11.4-14.5), respectively. In multivariable analyses, the highest rate of P. vivax parasitaemia over 42 days of follow-up was in patients residing in areas of short relapse periodicity (adjusted hazard ratio [AHR] = 6.2, 95% CI 2.0-19.5; p = 0.002); patients treated with AL (AHR = 6.2, 95% CI 4.6-8.5; p < 0.001), AA (AHR = 2.3, 95% CI 1.4-3.7; p = 0.001), or AM (AHR = 1.4, 95% CI 1.0-1.9; p = 0.028) compared with DP; and patients who did not clear their initial parasitaemia within 2 days (AHR = 1.8, 95% CI 1.4-2.3; p < 0.001). The analysis was limited by heterogeneity between study populations and lack of data from very low transmission settings. Study quality was high. CONCLUSIONS: In this meta-analysis, we found a high risk of P. vivax parasitaemia after treatment of P. falciparum malaria that varied significantly between studies. These P. vivax infections are likely attributable to relapses that could be prevented with radical cure including a hypnozoitocidal agent; however, the benefits of such a novel strategy will vary considerably between geographical areas.

Analysis of Plasmodium vivax schizont transcriptomes from field isolates reveals heterogeneity of expression of genes involved in host-parasite interactions.

Plasmodium vivax gene regulation remains difficult to study due to the lack of a robust in vitro culture method, low parasite densities in peripheral circulation and asynchronous parasite development. We adapted an RNA-seq protocol "DAFT-seq" to sequence the transcriptome of four P. vivax field isolates that were cultured for a short period ex vivo before using a density gradient for schizont enrichment. Transcription was detected from 78% of the PvP01 reference genome, despite being schizont-enriched samples. This extensive data was used to define thousands of 5' and 3' untranslated regions, some of which overlapped with neighbouring transcripts, and to improve the gene models of 352 genes, including identifying 20 novel gene transcripts. This dataset has also significantly increased the known amount of heterogeneity between P. vivax schizont transcriptomes from individual patients. The majority of genes found to be differentially expressed between the isolates lack Plasmodium falciparum homologs and are predicted to be involved in host-parasite interactions, with an enrichment in reticulocyte binding proteins, merozoite surface proteins and exported proteins with unknown function. An improved understanding of the diversity within P. vivax transcriptomes will be essential for the prioritisation of novel vaccine targets.

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.

Detecting geospatial patterns of Plasmodium falciparum parasite migration in Cambodia using optimized estimated effective migration surfaces.

BACKGROUND: Understanding the genetic structure of natural populations provides insight into the demographic and adaptive processes that have affected those populations. Such information, particularly when integrated with geospatial data, can have translational applications for a variety of fields, including public health. Estimated effective migration surfaces (EEMS) is an approach that allows visualization of the spatial patterns in genomic data to understand population structure and migration. In this study, we developed a workflow to optimize the resolution of spatial grids used to generate EEMS migration maps and applied this optimized workflow to estimate migration of Plasmodium falciparum in Cambodia and bordering regions of Thailand and Vietnam. METHODS: The optimal density of EEMS grids was determined based on a new workflow created using density clustering to define genomic clusters and the spatial distance between genomic clusters. Topological skeletons were used to capture the spatial distribution for each genomic cluster and to determine the EEMS grid density; i.e., both genomic and spatial clustering were used to guide the optimization of EEMS grids. Model accuracy for migration estimates using the optimized workflow was tested and compared to grid resolutions selected without the optimized workflow. As a test case, the optimized workflow was applied to genomic data generated from P. falciparum sampled in Cambodia and bordering regions, and migration maps were compared to estimates of malaria endemicity, as well as geographic properties of the study area, as a means of validating observed migration patterns. RESULTS: Optimized grids displayed both high model accuracy and reduced computing time compared to grid densities selected in an unguided manner. In addition, EEMS migration maps generated for P. falciparum using the optimized grid corresponded to estimates of malaria endemicity and geographic properties of the study region that might be expected to impact malaria parasite migration, supporting the validity of the observed migration patterns. CONCLUSIONS: Optimized grids reduce spatial uncertainty in the EEMS contours that can result from user-defined parameters, such as the resolution of the spatial grid used in the model. This workflow will be useful to a broad range of EEMS users as it can be applied to analyses involving other organisms of interest and geographic areas.

Triple artemisinin-based combination therapies versus artemisinin-based combination therapies for uncomplicated Plasmodium falciparum malaria: a multicentre, open-label, randomised clinical trial.

BACKGROUND: Artemisinin and partner-drug resistance in Plasmodium falciparum are major threats to malaria control and elimination. Triple artemisinin-based combination therapies (TACTs), which combine existing co-formulated ACTs with a second partner drug that is slowly eliminated, might provide effective treatment and delay emergence of antimalarial drug resistance. METHODS: In this multicentre, open-label, randomised trial, we recruited patients with uncomplicated P falciparum malaria at 18 hospitals and health clinics in eight countries. Eligible patients were aged 2-65 years, with acute, uncomplicated P falciparum malaria alone or mixed with non-falciparum species, and a temperature of 37·5°C or higher, or a history of fever in the past 24 h. Patients were randomly assigned (1:1) to one of two treatments using block randomisation, depending on their location: in Thailand, Cambodia, Vietnam, and Myanmar patients were assigned to either dihydroartemisinin-piperaquine or dihydroartemisinin-piperaquine plus mefloquine; at three sites in Cambodia they were assigned to either artesunate-mefloquine or dihydroartemisinin-piperaquine plus mefloquine; and in Laos, Myanmar, Bangladesh, India, and the Democratic Republic of the Congo they were assigned to either artemether-lumefantrine or artemether-lumefantrine plus amodiaquine. All drugs were administered orally and doses varied by drug combination and site. Patients were followed-up weekly for 42 days. The primary endpoint was efficacy, defined by 42-day PCR-corrected adequate clinical and parasitological response. Primary analysis was by intention to treat. A detailed assessment of safety and tolerability of the study drugs was done in all patients randomly assigned to treatment. This study is registered at ClinicalTrials.gov, NCT02453308, and is complete. FINDINGS: Between Aug 7, 2015, and Feb 8, 2018, 1100 patients were given either dihydroartemisinin-piperaquine (183 [17%]), dihydroartemisinin-piperaquine plus mefloquine (269 [24%]), artesunate-mefloquine (73 [7%]), artemether-lumefantrine (289 [26%]), or artemether-lumefantrine plus amodiaquine (286 [26%]). The median age was 23 years (IQR 13 to 34) and 854 (78%) of 1100 patients were male. In Cambodia, Thailand, and Vietnam the 42-day PCR-corrected efficacy after dihydroartemisinin-piperaquine plus mefloquine was 98% (149 of 152; 95% CI 94 to 100) and after dihydroartemisinin-piperaquine was 48% (67 of 141; 95% CI 39 to 56; risk difference 51%, 95% CI 42 to 59; p<0·0001). Efficacy of dihydroartemisinin-piperaquine plus mefloquine in the three sites in Myanmar was 91% (42 of 46; 95% CI 79 to 98) versus 100% (42 of 42; 95% CI 92 to 100) after dihydroartemisinin-piperaquine (risk difference 9%, 95% CI 1 to 17; p=0·12). The 42-day PCR corrected efficacy of dihydroartemisinin-piperaquine plus mefloquine (96% [68 of 71; 95% CI 88 to 99]) was non-inferior to that of artesunate-mefloquine (95% [69 of 73; 95% CI 87 to 99]) in three sites in Cambodia (risk difference 1%; 95% CI -6 to 8; p=1·00). The overall 42-day PCR-corrected efficacy of artemether-lumefantrine plus amodiaquine (98% [281 of 286; 95% CI 97 to 99]) was similar to that of artemether-lumefantrine (97% [279 of 289; 95% CI 94 to 98]; risk difference 2%, 95% CI -1 to 4; p=0·30). Both TACTs were well tolerated, although early vomiting (within 1 h) was more frequent after dihydroartemisinin-piperaquine plus mefloquine (30 [3·8%] of 794) than after dihydroartemisinin-piperaquine (eight [1·5%] of 543; p=0·012). Vomiting after artemether-lumefantrine plus amodiaquine (22 [1·3%] of 1703) and artemether-lumefantrine (11 [0·6%] of 1721) was infrequent. Adding amodiaquine to artemether-lumefantrine extended the electrocardiogram corrected QT interval (mean increase at 52 h compared with baseline of 8·8 ms [SD 18·6] vs 0·9 ms [16·1]; p<0·01) but adding mefloquine to dihydroartemisinin-piperaquine did not (mean increase of 22·1 ms [SD 19·2] for dihydroartemisinin-piperaquine vs 20·8 ms [SD 17·8] for dihydroartemisinin-piperaquine plus mefloquine; p=0·50). INTERPRETATION: Dihydroartemisinin-piperaquine plus mefloquine and artemether-lumefantrine plus amodiaquine TACTs are efficacious, well tolerated, and safe treatments of uncomplicated P falciparum malaria, including in areas with artemisinin and ACT partner-drug resistance. FUNDING: UK Department for International Development, Wellcome Trust, Bill & Melinda Gates Foundation, UK Medical Research Council, and US National Institutes of Health.

Evolution and expansion of multidrug-resistant malaria in southeast Asia: a genomic epidemiology study.

BACKGROUND: A multidrug-resistant co-lineage of Plasmodium falciparum malaria, named KEL1/PLA1, spread across Cambodia in 2008-13, causing high rates of treatment failure with the frontline combination therapy dihydroartemisinin-piperaquine. Here, we report on the evolution and spread of KEL1/PLA1 in subsequent years. METHODS: For this genomic epidemiology study, we analysed whole genome sequencing data from P falciparum clinical samples collected from patients with malaria between 2007 and 2018 from Cambodia, Laos, northeastern Thailand, and Vietnam, through the MalariaGEN P falciparum Community Project. Previously unpublished samples were provided by two large-scale multisite projects: the Tracking Artemisinin Resistance Collaboration II (TRAC2) and the Genetic Reconnaissance in the Greater Mekong Subregion (GenRe-Mekong) project. By investigating genome-wide relatedness between parasites, we inferred patterns of shared ancestry in the KEL1/PLA1 population. FINDINGS: We analysed 1673 whole genome sequences that passed quality filters, and determined KEL1/PLA1 status in 1615. Before 2009, KEL1/PLA1 was only found in western Cambodia; by 2016-17 its prevalence had risen to higher than 50% in all of the surveyed countries except for Laos. In northeastern Thailand and Vietnam, KEL1/PLA1 exceeded 80% of the most recent P falciparum parasites. KEL1/PLA1 parasites maintained high genetic relatedness and low diversity, reflecting a recent common origin. Several subgroups of highly related parasites have recently emerged within this co-lineage, with diverse geographical distributions. The three largest of these subgroups (n=84, n=79, and n=47) mostly emerged since 2016 and were all present in Cambodia, Laos, and Vietnam. These expanding subgroups carried new mutations in the crt gene, which arose on a specific genetic background comprising multiple genomic regions. Four newly emerging crt mutations were rare in the early period and became more prevalent by 2016-17 (Thr93Ser, rising to 19·8%; His97Tyr to 11·2%; Phe145Ile to 5·5%; and Ile218Phe to 11·1%). INTERPRETATION: After emerging and circulating for several years within Cambodia, the P falciparum KEL1/PLA1 co-lineage diversified into multiple subgroups and acquired new genetic features, including novel crt mutations. These subgroups have rapidly spread into neighbouring countries, suggesting enhanced fitness. These findings highlight the urgent need for elimination of this increasingly drug-resistant parasite co-lineage, and the importance of genetic surveillance in accelerating malaria elimination efforts. FUNDING: Wellcome Trust, Bill & Melinda Gates Foundation, UK Medical Research Council, and UK Department for International Development.

Determinants of dihydroartemisinin-piperaquine treatment failure in Plasmodium falciparum malaria in Cambodia, Thailand, and Vietnam: a prospective clinical, pharmacological, and genetic study.

BACKGROUND: The emergence and spread of resistance in Plasmodium falciparum malaria to artemisinin combination therapies in the Greater Mekong subregion poses a major threat to malaria control and elimination. The current study is part of a multi-country, open-label, randomised clinical trial (TRACII, 2015-18) evaluating the efficacy, safety, and tolerability of triple artemisinin combination therapies. A very high rate of treatment failure after treatment with dihydroartemisinin-piperaquine was observed in Thailand, Cambodia, and Vietnam. The immediate public health importance of our findings prompted us to report the efficacy data on dihydroartemisinin-piperaquine and its determinants ahead of the results of the overall trial, which will be published later this year. METHODS: Patients aged between 2 and 65 years presenting with uncomplicated P falciparum or mixed species malaria at seven sites in Thailand, Cambodia, and Vietnam were randomly assigned to receive dihydroartemisinin-piperaquine with or without mefloquine, as part of the TRACII trial. The primary outcome was the PCR-corrected efficacy at day 42. Next-generation sequencing was used to assess the prevalence of molecular markers associated with artemisinin resistance (kelch13 mutations, in particular Cys580Tyr) and piperaquine resistance (plasmepsin-2 and plasmepsin-3 amplifications and crt mutations). This study is registered with ClinicalTrials.gov, number NCT02453308. FINDINGS: Between Sept 28, 2015, and Jan 18, 2018, 539 patients with acute P falciparum malaria were screened for eligibility, 292 were enrolled, and 140 received dihydroartemisinin-piperaquine. The overall Kaplan-Meier estimate of PCR-corrected efficacy of dihydroartemisinin-piperaquine at day 42 was 50·0% (95% CI 41·1-58·3). PCR-corrected efficacies for individual sites were 12·7% (2·2-33·0) in northeastern Thailand, 38·2% (15·9-60·5) in western Cambodia, 73·4% (57·0-84·3) in Ratanakiri (northeastern Cambodia), and 47·1% (33·5-59·6) in Binh Phuoc (southwestern Vietnam). Treatment failure was associated independently with plasmepsin2/3 amplification status and four mutations in the crt gene (Thr93Ser, His97Tyr, Phe145Ile, and Ile218Phe). Compared with the results of our previous TRACI trial in 2011-13, the prevalence of molecular markers of artemisinin resistance (kelch13 Cys580Tyr mutations) and piperaquine resistance (plasmepsin2/3 amplifications and crt mutations) has increased substantially in the Greater Mekong subregion in the past decade. INTERPRETATION: Dihydroartemisinin-piperaquine is not treating malaria effectively across the eastern Greater Mekong subregion. A highly drug-resistant P falciparum co-lineage is evolving, acquiring new resistance mechanisms, and spreading. Accelerated elimination of P falciparum malaria in this region is needed urgently, to prevent further spread and avoid a potential global health emergency. FUNDING: UK Department for International Development, Wellcome Trust, Bill & Melinda Gates Foundation, Medical Research Council, and National Institutes of Health.

Contribution of Functional Antimalarial Immunity to Measures of Parasite Clearance in Therapeutic Efficacy Studies of Artemisinin Derivatives.

BACKGROUND: Antibodies to the blood stages of malaria parasites enhance parasite clearance and antimalarial efficacy. The antibody subclass and functions that contribute to parasite clearance during antimalarial treatment and their relationship to malaria transmission intensity have not been characterized. METHODS: Levels of immunoglobulin G (IgG) subclasses and C1q fixation in response to Plasmodium falciparum merozoite antigens (erythrocyte-binding antigen [EBA] 175RIII-V, merozoite surface protein 2 [MSP-2], and MSP-142) and opsonic phagocytosis of merozoites were measured in a multinational trial assessing the efficacy of artesunate therapy across 11 Southeast Asian sites. Regression analyses assessed the effects of antibody seropositivity on the parasite clearance half-life (PC½), having a PC½ of ≥5 hours, and having parasitemia 3 days after treatment. RESULTS: IgG3, followed by IgG1, was the predominant IgG subclass detected (seroprevalence range, 5%-35% for IgG1 and 27%-41% for IgG3), varied across study sites, and was lowest in study sites with the lowest transmission intensity and slowest mean PC½. IgG3, C1q fixation, and opsonic-phagocytosis seropositivity were associated with a faster PC½ (range of the mean reduction in PC½, 0.47-1.16 hours; P range, .001-.03) and a reduced odds of having a PC½ of ≥5 hours and having parasitemia 3 days after treatment. CONCLUSIONS: The prevalence of IgG3, complement-fixing antibodies, and merozoite phagocytosis vary according to transmission intensity, are associated with faster parasite clearance, and may be sensitive surrogates of an augmented clearance capacity of infected erythrocytes. Determining the functional immune mechanisms associated with parasite clearance will improve characterization of artemisinin resistance.

Identification and Characterization of Functional Human Monoclonal Antibodies to Plasmodium vivax Duffy-Binding Protein.

Plasmodium vivax invasion of reticulocytes relies on distinct receptor-ligand interactions between the parasite and host erythrocytes. Engagement of the highly polymorphic domain II of the P. vivax Duffy-binding protein (DBPII) with the erythrocyte's Duffy Ag receptor for chemokines (DARC) is essential. Some P. vivax-exposed individuals acquired Abs to DBPII that block DBPII-DARC interaction and inhibit P. vivax reticulocyte invasion, and Ab levels correlate with protection against P. vivax malaria. To better understand the functional characteristics and fine specificity of protective human Abs to DBPII, we sorted single DBPII-specific IgG+ memory B cells from three individuals with high blocking activity to DBPII. We identified 12 DBPII-specific human mAbs from distinct lineages that blocked DBPII-DARC binding. All mAbs were P. vivax strain transcending and targeted known binding motifs of DBPII with DARC. Eleven mAbs competed with each other for binding, indicating recognition of the same or overlapping epitopes. Naturally acquired blocking Abs to DBPII from individuals with high levels residing in different P. vivax-endemic areas worldwide competed with mAbs, suggesting broadly shared recognition sites. We also found that mAbs inhibited P. vivax entry into reticulocytes in vitro. These findings suggest that IgG+ memory B cell activity in individuals with P. vivax strain-transcending Abs to DBPII display a limited clonal response with inhibitory blocking directed against a distinct region of the molecule.