Molecular markers of artemisinin resistance during falciparum malaria elimination in Eastern Myanmar.

BACKGROUND: Artemisinin resistance in Plasmodium falciparum threatens global malaria elimination efforts. To contain and then eliminate artemisinin resistance in Eastern Myanmar a network of community-based malaria posts was instituted and targeted mass drug administration (MDA) with dihydroartemisinin-piperaquine (three rounds at monthly intervals) was conducted. The prevalence of artemisinin resistance during the elimination campaign (2013-2019) was characterized. METHODS: Throughout the six-year campaign Plasmodium falciparum positive blood samples from symptomatic patients and from cross-sectional surveys were genotyped for mutations in kelch-13-a molecular marker of artemisinin resistance. RESULT: The program resulted in near elimination of falciparum malaria. Of 5162 P. falciparum positive blood samples genotyped, 3281 (63.6%) had K13 mutations. The prevalence of K13 mutations was 73.9% in 2013 and 64.4% in 2019. Overall, there was a small but significant decline in the proportion of K13 mutants (p < 0.001). In the MDA villages there was no significant change in the K13 proportions before and after MDA. The distribution of different K13 mutations changed substantially; F446I and P441L mutations increased in both MDA and non-MDA villages, while most other K13 mutations decreased. The proportion of C580Y mutations fell from 9.2% (43/467) before MDA to 2.3% (19/813) after MDA (p < 0.001). Similar changes occurred in the 487 villages where MDA was not conducted. CONCLUSION: The malaria elimination program in Kayin state, eastern Myanmar, led to a substantial reduction in falciparum malaria. Despite the intense use of artemisinin-based combination therapies, both in treatment and MDA, this did not select for artemisinin resistance.

Short tandem repeat polymorphism in the promoter region of cyclophilin 19B drives its transcriptional upregulation and contributes to drug resistance in the malaria parasite Plasmodium falciparum.

Resistance of the human malaria parasites, Plasmodium falciparum, to artemisinins is now fully established in Southeast Asia and is gradually emerging in Sub-Saharan Africa. Although nonsynonymous SNPs in the pfk13 Kelch-repeat propeller (KREP) domain are clearly associated with artemisinin resistance, their functional relevance requires cooperation with other genetic factors/alterations of the P. falciparum genome, collectively referred to as genetic background. Here we provide experimental evidence that P. falciparum cyclophilin 19B (PfCYP19B) may represent one putative factor in this genetic background, contributing to artemisinin resistance via its increased expression. We show that overexpression of PfCYP19B in vitro drives limited but significant resistance to not only artemisinin but also piperaquine, an important partner drug in artemisinin-based combination therapies. We showed that PfCYP19B acts as a negative regulator of the integrated stress response (ISR) pathway by modulating levels of phosphorylated eIF2α (eIF2α-P). Curiously, artemisinin and piperaquine affect eIF2α-P in an inverse direction that in both cases can be modulated by PfCYP19B towards resistance. Here we also provide evidence that the upregulation of PfCYP19B in the drug-resistant parasites appears to be maintained by a short tandem repeat (SRT) sequence polymorphism in the gene's promoter region. These results support a model that artemisinin (and other drugs) resistance mechanisms are complex genetic traits being contributed to by altered expression of multiple genes driven by genetic polymorphism at their promoter regions.

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.

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

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

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.

Artemisinin resistance in the malaria parasite, Plasmodium falciparum, originates from its initial transcriptional response.

The emergence and spread of artemisinin-resistant Plasmodium falciparum, first in the Greater Mekong Subregion (GMS), and now in East Africa, is a major threat to global malaria elimination ambitions. To investigate the artemisinin resistance mechanism, transcriptome analysis was conducted of 577 P. falciparum isolates collected in the GMS between 2016-2018. A specific artemisinin resistance-associated transcriptional profile was identified that involves a broad but discrete set of biological functions related to proteotoxic stress, host cytoplasm remodelling, and REDOX metabolism. The artemisinin resistance-associated transcriptional profile evolved from initial transcriptional responses of susceptible parasites to artemisinin. The genetic basis for this adapted response is likely to be complex.

Triple therapy with artemether-lumefantrine plus amodiaquine versus artemether-lumefantrine alone for artemisinin-resistant, uncomplicated falciparum malaria: an open-label, randomised, multicentre trial.

BACKGROUND: Late treatment failures after artemisinin-based combination therapies (ACTs) for falciparum malaria have increased in the Greater Mekong subregion in southeast Asia. Addition of amodiaquine to artemether-lumefantrine could provide an efficacious treatment for multidrug-resistant infections. METHODS: We conducted an open-label, randomised trial at five hospitals or health centres in three locations (western Cambodia, eastern Cambodia, and Vietnam). Eligible participants were male and female patients aged 2-65 years with uncomplicated Plasmodium falciparum malaria. Patients were randomly allocated (1:1 in blocks of eight to 12) to either artemether-lumefantrine alone (dosed according to WHO guidelines) or artemether-lumefantrine plus amodiaquine (10 mg base per kg/day), both given orally as six doses over 3 days. All received a single dose of primaquine (0·25 mg/kg) 24 h after the start of study treatment to limit transmission of the parasite. Parasites were genotyped, identifying artemisinin resistance. The primary outcome was Kaplan-Meier 42-day PCR-corrected efficacy against recrudescence of the original parasite, assessed by intent-to-treat. Safety was a secondary outcome. This completed trial is registered at ClinicalTrials.gov (NCT03355664). FINDINGS: Between March 18, 2018, and Jan 30, 2020, 310 patients received randomly allocated treatment; 154 received artemether-lumefantrine alone and 156 received artemether-lumefantrine plus amodiaquine. Parasites from 305 of these patients were genotyped. 42-day PCR-corrected treatment efficacy was noted in 151 (97%, 95% CI 92-99) of 156 patients with artemether-lumefantrine plus amodiaquine versus 146 (95%, 89-97) of 154 patients with artemether-lumefantrine alone; hazard ratio (HR) for recrudescence 0·6 (95% CI 0·2-1·9, p=0·38). Of the 13 recrudescences, 12 were in 174 (57%) of 305 infections with pfkelch13 mutations indicating artemisinin resistance, for which 42-day efficacy was noted in 89 (96%) of 93 infections with artemether-lumefantrine plus amodiaquine versus 73 (90%) of 81 infections with artemether-lumefantrine alone; HR for recrudescence 0·44 (95% CI 0·14-1·40, p=0·17). Artemether-lumefantrine plus amodiaquine was generally well tolerated, but the number of mild (grade 1-2) adverse events, mainly gastrointestinal, was greater in this group compared with artemether-lumefantrine alone (vomiting, 12 [8%] with artemether-lumefantrine plus amodiaquine vs three [2%] with artemether-lumefantrine alone, p=0·03; and nausea, 11 [7%] with artemether-lumefantrine plus amodiaquine vs three [2%] with artemether-lumefantrine alone, p=0·05). Early vomiting within 1 h of treatment, requiring retreatment, occurred in no patients of 154 with artemether-lumefantrine alone versus five (3%) of 156 with artemether-lumefantrine plus amodiaquine, p=0·06. Bradycardia (≤54 beats/min) of any grade was noted in 59 (38%) of 154 patients with artemether-lumefantrine alone and 95 (61%) of 156 with artemether-lumefantrine plus amodiaquine, p=0·0001. INTERPRETATION: Artemether-lumefantrine plus amodiaquine provides an alternative to artemether-lumefantrine alone as first-line treatment for multidrug-resistant P falciparum malaria in the Greater Mekong subregion, and could prolong the therapeutic lifetime of artemether-lumefantrine in malaria-endemic populations. FUNDING: Bill & Melinda Gates Foundation, Wellcome Trust.

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.

Arterolane-piperaquine-mefloquine versus arterolane-piperaquine and artemether-lumefantrine in the treatment of uncomplicated Plasmodium falciparum malaria in Kenyan children: a single-centre, open-label, randomised, non-inferiority trial.

BACKGROUND: Triple antimalarial combination therapies combine potent and rapidly cleared artemisinins or related synthetic ozonides, such as arterolane, with two, more slowly eliminated partner drugs to reduce the risk of resistance. We aimed to assess the safety, tolerability, and efficacy of arterolane-piperaquine-mefloquine versus arterolane-piperaquine and artemether-lumefantrine for the treatment of uncomplicated falciparum malaria in Kenyan children. METHODS: In this single-centre, open-label, randomised, non-inferiority trial done in Kilifi County Hospital, Kilifi, coastal Kenya, children with uncomplicated Plasmodium falciparum malaria were recruited. Eligible patients were aged 2-12 years and had an asexual parasitaemia of 5000-250 000 parasites per µL. The exclusion criteria included the presence of an acute illness other than malaria, the inability to tolerate oral medications, treatment with an artemisinin derivative in the previous 7 days, a known hypersensitivity or contraindication to any of the study drugs, and a QT interval corrected for heart rate (QTc interval) longer than 450 ms. Patients were randomly assigned (1:1:1), by use of blocks of six, nine, and 12, and opaque, sealed, and sequentially numbered envelopes, to receive either arterolane-piperaquine, arterolane-piperaquine-mefloquine, or artemether-lumefantrine. Laboratory staff, but not the patients, the patients' parents or caregivers, clinical or medical officers, nurses, or trial statistician, were masked to the intervention groups. For 3 days, oral artemether-lumefantrine was administered twice daily (target dose 5-24 mg/kg of bodyweight of artemether and 29-144 mg/kg of bodyweight of lumefantrine), and oral arterolane-piperaquine (arterolane dose 4 mg/kg of bodyweight; piperaquine dose 20 mg/kg of bodyweight) and oral arterolane-piperaquine-mefloquine (mefloquine dose 8 mg/kg of bodyweight) were administered once daily. All patients received 0·25 mg/kg of bodyweight of oral primaquine at hour 24. All patients were admitted to Kilifi County Hospital for at least 3 consecutive days and followed up at day 7 and, thereafter, weekly for up to 42 days. The primary endpoint was 42-day PCR-corrected efficacy, defined as the absence of treatment failure in the first 42 days post-treatment, of arterolane-piperaquine-mefloquine versus artemether-lumefantrine, and, along with safety, was analysed in the intention-to-treat population, which comprised all patients who received at least one dose of a study drug. The 42-day PCR-corrected efficacy of arterolane-piperaquine-mefloquine versus arterolane-piperaquine was an important secondary endpoint and was also analysed in the intention-to-treat population. The non-inferiority margin for the risk difference between treatments was -7%. The study is registered in ClinicalTrials.gov, NCT03452475, and is completed. FINDINGS: Between March 7, 2018, and May 2, 2019, 533 children with P falciparum were screened, of whom 217 were randomly assigned to receive either arterolane-piperaquine (n=73), arterolane-piperaquine-mefloquine (n=72), or artemether-lumefantrine (n=72) and comprised the intention-to-treat population. The 42-day PCR-corrected efficacy after treatment with arterolane-piperaquine-mefloquine (100%, 95% CI 95-100; 72/72) was non-inferior to that after treatment with artemether-lumefantrine (96%, 95% CI 88-99; 69/72; risk difference 4%, 95% CI 0-9; p=0·25). The 42-day PCR-corrected efficacy of arterolane-piperaquine-mefloquine was non-inferior to that of arterolane-piperaquine (100%, 95% CI 95-100; 73/73; risk difference 0%). Vomiting rates in the first hour post-drug administration were significantly higher in patients treated with arterolane-piperaquine (5%, 95% CI 2-9; ten of 203 drug administrations; p=0·0013) or arterolane-piperaquine-mefloquine (5%, 3-9; 11 of 209 drug administrations; p=0·0006) than in patients treated with artemether-lumefantrine (1%, 0-2; three of 415 drug administrations). Upper respiratory tract complaints (n=26 for artemether-lumefantrine; n=19 for arterolane-piperaquine-mefloquine; n=23 for arterolane-piperaquine), headache (n=13; n=4; n=5), and abdominal pain (n=7; n=5; n=5) were the most frequently reported adverse events. There were no deaths. INTERPRETATION: This study shows that arterolane-piperaquine-mefloquine is an efficacious and safe treatment for uncomplicated falciparum malaria in children and could potentially be used to prevent or delay the emergence of antimalarial resistance. FUNDING: UK Department for International Development, The Wellcome Trust, The Bill & Melinda Gates Foundation, Sun Pharmaceutical Industries.

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.

Emergence of artemisinin-resistant Plasmodium falciparum with kelch13 C580Y mutations on the island of New Guinea.

The rapid and aggressive spread of artemisinin-resistant Plasmodium falciparum carrying the C580Y mutation in the kelch13 gene is a growing threat to malaria elimination in Southeast Asia, but there is no evidence of their spread to other regions. We conducted cross-sectional surveys in 2016 and 2017 at two clinics in Wewak, Papua New Guinea (PNG) where we identified three infections caused by C580Y mutants among 239 genotyped clinical samples. One of these mutants exhibited the highest survival rate (6.8%) among all parasites surveyed in ring-stage survival assays (RSA) for artemisinin. Analyses of kelch13 flanking regions, and comparisons of deep sequencing data from 389 clinical samples from PNG, Indonesian Papua and Western Cambodia, suggested an independent origin of the Wewak C580Y mutation, showing that the mutants possess several distinctive genetic features. Identity by descent (IBD) showed that multiple portions of the mutants' genomes share a common origin with parasites found in Indonesian Papua, comprising several mutations within genes previously associated with drug resistance, such as mdr1, ferredoxin, atg18 and pnp. These findings suggest that a P. falciparum lineage circulating on the island of New Guinea has gradually acquired a complex ensemble of variants, including kelch13 C580Y, which have affected the parasites' drug sensitivity. This worrying development reinforces the need for increased surveillance of the evolving parasite populations on the island, to contain the spread of resistance.

Molecular epidemiology of resistance to antimalarial drugs in the Greater Mekong subregion: an observational study.

BACKGROUND: The Greater Mekong subregion is a recurrent source of antimalarial drug resistance in Plasmodium falciparum malaria. This study aimed to characterise the extent and spread of resistance across this entire region between 2007 and 2018. METHODS: P falciparum isolates from Myanmar, Thailand, Laos, and Cambodia were obtained from clinical trials and epidemiological studies done between Jan 1, 2007, and Dec 31, 2018, and were genotyped for molecular markers (pfkelch, pfcrt, pfplasmepsin2, and pfmdr1) of antimalarial drug resistance. Genetic relatedness was assessed using microsatellite and single nucleotide polymorphism typing of flanking sequences around target genes. FINDINGS: 10 632 isolates were genotyped. A single long pfkelch Cys580Tyr haplotype (from -50 kb to +31·5 kb) conferring artemisinin resistance (PfPailin) now dominates across the eastern Greater Mekong subregion. Piperaquine resistance associated with pfplasmepsin2 gene amplification and mutations in pfcrt downstream of the Lys76Thr chloroquine resistance locus has also developed. On the Thailand-Myanmar border a different pfkelch Cys580Tyr lineage rose to high frequencies before it was eliminated. Elsewhere in Myanmar the Cys580Tyr allele remains widespread at low allele frequencies. Meanwhile a single artemisinin-resistant pfkelch Phe446Ile haplotype has spread across Myanmar. Despite intense use of dihydroartemisinin-piperaquine in Kayin state, eastern Myanmar, both in treatment and mass drug administrations, no selection of piperaquine resistance markers was observed. pfmdr1 amplification, a marker of resistance to mefloquine, remains at low prevalence across the entire region. INTERPRETATION: Artemisinin resistance in P falciparum is now prevalent across the Greater Mekong subregion. In the eastern Greater Mekong subregion a multidrug resistant P falciparum lineage (PfPailin) dominates. In Myanmar a long pfkelch Phe446Ile haplotype has spread widely but, by contrast with the eastern Greater Mekong subregion, there is no indication of artemisinin combination therapy (ACT) partner drug resistance from genotyping known markers, and no evidence of spread of ACT resistant P falciparum from the east to the west. There is still a window of opportunity to prevent global spread of ACT resistance. FUNDING: Thailand Science Research and Innovation, Initiative 5%, Expertise France, Wellcome Trust.

Implementing parasite genotyping into national surveillance frameworks: feedback from control programmes and researchers in the Asia-Pacific region.

The Asia-Pacific region faces formidable challenges in achieving malaria elimination by the proposed target in 2030. Molecular surveillance of Plasmodium parasites can provide important information on malaria transmission and adaptation, which can inform national malaria control programmes (NMCPs) in decision-making processes. In November 2019 a parasite genotyping workshop was held in Jakarta, Indonesia, to review molecular approaches for parasite surveillance and explore ways in which these tools can be integrated into public health systems and inform policy. The meeting was attended by 70 participants from 8 malaria-endemic countries and partners of the Asia Pacific Malaria Elimination Network. The participants acknowledged the utility of multiple use cases for parasite genotyping including: quantifying the prevalence of drug resistant parasites, predicting risks of treatment failure, identifying major routes and reservoirs of infection, monitoring imported malaria and its contribution to local transmission, characterizing the origins and dynamics of malaria outbreaks, and estimating the frequency of Plasmodium vivax relapses. However, the priority of each use case varies with different endemic settings. Although a one-size-fits-all approach to molecular surveillance is unlikely to be applicable across the Asia-Pacific region, consensus on the spectrum of added-value activities will help support data sharing across national boundaries. Knowledge exchange is needed to establish local expertise in different laboratory-based methodologies and bioinformatics processes. Collaborative research involving local and international teams will help maximize the impact of analytical outputs on the operational needs of NMCPs. Research is also needed to explore the cost-effectiveness of genetic epidemiology for different use cases to help to leverage funding for wide-scale implementation. Engagement between NMCPs and local researchers will be critical throughout this process.

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.

Mapping the travel patterns of people with malaria in Bangladesh.

BACKGROUND: Spread of malaria and antimalarial resistance through human movement present major threats to current goals to eliminate the disease. Bordering the Greater Mekong Subregion, southeast Bangladesh is a potentially important route of spread to India and beyond, but information on travel patterns in this area are lacking. METHODS: Using a standardised short survey tool, 2090 patients with malaria were interviewed at 57 study sites in 2015-2016 about their demographics and travel patterns in the preceding 2 months. RESULTS: Most travel was in the south of the study region between Cox's Bazar district (coastal region) to forested areas in Bandarban (31% by days and 45% by nights), forming a source-sink route. Less than 1% of travel reported was between the north and south forested areas of the study area. Farmers (21%) and students (19%) were the top two occupations recorded, with 67 and 47% reporting travel to the forest respectively. Males aged 25-49 years accounted for 43% of cases visiting forests but only 24% of the study population. Children did not travel. Women, forest dwellers and farmers did not travel beyond union boundaries. Military personnel travelled the furthest especially to remote forested areas. CONCLUSIONS: The approach demonstrated here provides a framework for identifying key traveller groups and their origins and destinations of travel in combination with knowledge of local epidemiology to inform malaria control and elimination efforts. Working with the NMEP, the findings were used to derive a set of policy recommendations to guide targeting of interventions for elimination.

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.

Mapping imported malaria in Bangladesh using parasite genetic and human mobility data.

For countries aiming for malaria elimination, travel of infected individuals between endemic areas undermines local interventions. Quantifying parasite importation has therefore become a priority for national control programs. We analyzed epidemiological surveillance data, travel surveys, parasite genetic data, and anonymized mobile phone data to measure the spatial spread of malaria parasites in southeast Bangladesh. We developed a genetic mixing index to estimate the likelihood of samples being local or imported from parasite genetic data and inferred the direction and intensity of parasite flow between locations using an epidemiological model integrating the travel survey and mobile phone calling data. Our approach indicates that, contrary to dogma, frequent mixing occurs in low transmission regions in the southwest, and elimination will require interventions in addition to reducing imported infections from forested regions. Unlike risk maps generated from clinical case counts alone, therefore, our approach distinguishes areas of frequent importation as well as high transmission.

The impact of targeted malaria elimination with mass drug administrations on falciparum malaria in Southeast Asia: A cluster randomised trial.

BACKGROUND: The emergence and spread of multidrug-resistant Plasmodium falciparum in the Greater Mekong Subregion (GMS) threatens global malaria elimination efforts. Mass drug administration (MDA), the presumptive antimalarial treatment of an entire population to clear the subclinical parasite reservoir, is a strategy to accelerate malaria elimination. We report a cluster randomised trial to assess the effectiveness of dihydroartemisinin-piperaquine (DP) MDA in reducing falciparum malaria incidence and prevalence in 16 remote village populations in Myanmar, Vietnam, Cambodia, and the Lao People's Democratic Republic, where artemisinin resistance is prevalent. METHODS AND FINDINGS: After establishing vector control and community-based case management and following intensive community engagement, we used restricted randomisation within village pairs to select 8 villages to receive early DP MDA and 8 villages as controls for 12 months, after which the control villages received deferred DP MDA. The MDA comprised 3 monthly rounds of 3 daily doses of DP and, except in Cambodia, a single low dose of primaquine. We conducted exhaustive cross-sectional surveys of the entire population of each village at quarterly intervals using ultrasensitive quantitative PCR to detect Plasmodium infections. The study was conducted between May 2013 and July 2017. The investigators randomised 16 villages that had a total of 8,445 residents at the start of the study. Of these 8,445 residents, 4,135 (49%) residents living in 8 villages, plus an additional 288 newcomers to the villages, were randomised to receive early MDA; 3,790 out of the 4,423 (86%) participated in at least 1 MDA round, and 2,520 out of the 4,423 (57%) participated in all 3 rounds. The primary outcome, P. falciparum prevalence by month 3 (M3), fell by 92% (from 5.1% [171/3,340] to 0.4% [12/2,828]) in early MDA villages and by 29% (from 7.2% [246/3,405] to 5.1% [155/3,057]) in control villages. Over the following 9 months, the P. falciparum prevalence increased to 3.3% (96/2,881) in early MDA villages and to 6.1% (128/2,101) in control villages (adjusted incidence rate ratio 0.41 [95% CI 0.20 to 0.84]; p = 0.015). Individual protection was proportional to the number of completed MDA rounds. Of 221 participants with subclinical P. falciparum infections who participated in MDA and could be followed up, 207 (94%) cleared their infections, including 9 of 10 with artemisinin- and piperaquine-resistant infections. The DP MDAs were well tolerated; 6 severe adverse events were detected during the follow-up period, but none was attributable to the intervention. CONCLUSIONS: Added to community-based basic malaria control measures, 3 monthly rounds of DP MDA reduced the incidence and prevalence of falciparum malaria over a 1-year period in areas affected by artemisinin resistance. P. falciparum infections returned during the follow-up period as the remaining infections spread and malaria was reintroduced from surrounding areas. Limitations of this study include a relatively small sample of villages, heterogeneity between villages, and mobility of villagers that may have limited the impact of the intervention. These results suggest that, if used as part of a comprehensive, well-organised, and well-resourced elimination programme, DP MDA can be a useful additional tool to accelerate malaria elimination. TRIAL REGISTRATION: ClinicalTrials.gov NCT01872702.

Genomic Analysis of Plasmodium vivax in Southern Ethiopia Reveals Selective Pressures in Multiple Parasite Mechanisms.

The Horn of Africa harbors the largest reservoir of Plasmodium vivax in the continent. Most of sub-Saharan Africa has remained relatively vivax-free due to a high prevalence of the human Duffy-negative trait, but the emergence of strains able to invade Duffy-negative reticulocytes poses a major public health threat. We undertook the first population genomic investigation of P. vivax from the region, comparing the genomes of 24 Ethiopian isolates against data from Southeast Asia to identify important local adaptions. The prevalence of the Duffy binding protein amplification in Ethiopia was 79%, potentially reflecting adaptation to Duffy negativity. There was also evidence of selection in a region upstream of the chloroquine resistance transporter, a putative chloroquine-resistance determinant. Strong signals of selection were observed in genes involved in immune evasion and regulation of gene expression, highlighting the need for a multifaceted intervention approach to combat P. vivax in the region.