Genomics of Plasmodium vivax in Colombia reveals evidence of local bottle-necking and inter-country connectivity in the Americas.

Colombia aims to eliminate malaria by 2030 but remains one of the highest burden countries in the Americas. Plasmodium vivax contributes half of all malaria cases, with its control challenged by relapsing parasitaemia, drug resistance and cross-border spread. Using 64 Colombian P. vivax genomes collected between 2013 and 2017, we explored diversity and selection in two major foci of transmission: Chocó and Córdoba. Open-access data from other countries were used for comparative assessment of drug resistance candidates and to assess cross-border spread. Across Colombia, polyclonal infections were infrequent (12%), and infection connectivity was relatively high (median IBD = 5%), consistent with low endemicity. Chocó exhibited a higher frequency of polyclonal infections (23%) than Córdoba (7%), although the difference was not significant (P = 0.300). Most Colombian infections carried double pvdhfr (95%) and single pvdhps (71%) mutants, but other drug resistance mutations were less prevalent (< 10%). There was no evidence of selection at the pvaat1 gene, whose P. falciparum orthologue has recently been implicated in chloroquine resistance. Global population comparisons identified other putative adaptations. Within the Americas, low-level connectivity was observed between Colombia and Peru, highlighting potential for cross-border spread. Our findings demonstrate the potential of molecular data to inform on infection spread and adaptation.

Genomic analysis of Plasmodium vivax describes patterns of connectivity and putative drivers of adaptation in Ethiopia.

Ethiopia has the greatest burden of Plasmodium vivax in Africa, but little is known about the epidemiological landscape of parasites across the country. We analysed the genomic diversity of 137 P. vivax isolates collected nine Ethiopian districts from 2012 to 2016. Signatures of selection were detected by cross-country comparisons with isolates from Thailand (n = 104) and Indonesia (n = 111), representing regions with low and high chloroquine resistance respectively. 26% (35/137) of Ethiopian infections were polyclonal, and 48.5% (17/35) of these comprised highly related clones (within-host identity-by-descent > 25%), indicating frequent co-transmission and superinfection. Parasite gene flow between districts could not be explained entirely by geographic distance, with economic and cultural factors hypothesised to have an impact on connectivity. Amplification of the duffy binding protein gene (pvdbp1) was prevalent across all districts (16-75%). Cross-population haplotype homozygosity revealed positive selection in a region proximal to the putative chloroquine resistance transporter gene (pvcrt-o). An S25P variant in amino acid transporter 1 (pvaat1), whose homologue has recently been implicated in P. falciparum chloroquine resistance evolution, was prevalent in Ethiopia (96%) but not Thailand or Indonesia (35-53%). The genomic architecture in Ethiopia highlights circulating variants of potential public health concern in an endemic setting with evidence of stable transmission.

Drug resistance and vaccine target surveillance of Plasmodium falciparum using nanopore sequencing in Ghana.

Malaria results in over 600,000 deaths annually, with the highest burden of deaths in young children living in sub-Saharan Africa. Molecular surveillance can provide important information for malaria control policies, including detection of antimalarial drug resistance. However, genome sequencing capacity in malaria-endemic countries is limited. We designed and implemented an end-to-end workflow to detect Plasmodium falciparum antimalarial resistance markers and diversity in the vaccine target circumsporozoite protein (csp) using nanopore sequencing in Ghana. We analysed 196 clinical samples and showed that our method is rapid, robust, accurate and straightforward to implement. Importantly, our method could be applied to dried blood spot samples, which are readily collected in endemic settings. We report that P. falciparum parasites in Ghana are mostly susceptible to chloroquine, with persistent sulfadoxine-pyrimethamine resistance and no evidence of artemisinin resistance. Multiple single nucleotide polymorphisms were identified in csp, but their significance is uncertain. Our study demonstrates the feasibility of nanopore sequencing for malaria genomic surveillance in endemic countries.

Chloroquine resistance evolution in Plasmodium falciparum is mediated by the putative amino acid transporter AAT1.

Malaria parasites break down host haemoglobin into peptides and amino acids in the digestive vacuole for export to the parasite cytoplasm for growth: interrupting this process is central to the mode of action of several antimalarial drugs. Mutations in the chloroquine (CQ) resistance transporter, pfcrt, located in the digestive vacuole membrane, confer CQ resistance in Plasmodium falciparum, and typically also affect parasite fitness. However, the role of other parasite loci in the evolution of CQ resistance is unclear. Here we use a combination of population genomics, genetic crosses and gene editing to demonstrate that a second vacuolar transporter plays a key role in both resistance and compensatory evolution. Longitudinal genomic analyses of the Gambian parasites revealed temporal signatures of selection on a putative amino acid transporter (pfaat1) variant S258L, which increased from 0% to 97% in frequency between 1984 and 2014 in parallel with the pfcrt1 K76T variant. Parasite genetic crosses then identified a chromosome 6 quantitative trait locus containing pfaat1 that is selected by CQ treatment. Gene editing demonstrated that pfaat1 S258L potentiates CQ resistance but at a cost of reduced fitness, while pfaat1 F313S, a common southeast Asian polymorphism, reduces CQ resistance while restoring fitness. Our analyses reveal hidden complexity in CQ resistance evolution, suggesting that pfaat1 may underlie regional differences in the dynamics of resistance evolution, and modulate parasite resistance or fitness by manipulating the balance between both amino acid and drug transport.

Generation of a mutator parasite to drive resistome discovery in Plasmodium falciparum.

In vitro evolution of drug resistance is a powerful approach for identifying antimalarial targets, however, key obstacles to eliciting resistance are the parasite inoculum size and mutation rate. Here we sought to increase parasite genetic diversity to potentiate resistance selections by editing catalytic residues of Plasmodium falciparum DNA polymerase δ. Mutation accumulation assays reveal a ~5-8 fold elevation in the mutation rate, with an increase of 13-28 fold in drug-pressured lines. Upon challenge with the spiroindolone PfATP4-inhibitor KAE609, high-level resistance is obtained more rapidly and at lower inocula than wild-type parasites. Selections also yield mutants with resistance to an "irresistible" compound, MMV665794 that failed to yield resistance with other strains. We validate mutations in a previously uncharacterised gene, PF3D7_1359900, which we term quinoxaline resistance protein (QRP1), as causal for resistance to MMV665794 and a panel of quinoxaline analogues. The increased genetic repertoire available to this "mutator" parasite can be leveraged to drive P. falciparum resistome discovery.

T-cell activation, senescence, and exhaustion in asymptomatic HIV/Leishmania infantum co-infection.

Background: Leishmania infantum is an opportunistic parasitic infection. An immunocompromised state increases the risk of converting asymptomatic infection to symptomatic visceral leishmaniasis (VL), which has a ~5% fatality rate even with treatment. HIV coinfection increases the risk of death from VL. Methods: A cross-sectional study was performed between 2014 and 2016 to determine the prevalence of L. infantum infection in HIV positive subjects residing in the state of Rio Grande do Norte, Brazil (n=1,372) and of these a subgroup of subjects were followed longitudinally. Subsequent incident cases of VL were ascertained from a public health database through 2018. A subgroup (n=69) of the cross-sectional study subjects was chosen to assess immune status (T cell activation, senescence, exhaustion) and outcome. The data were compared between asymptomatic HIV+/L. infantum+ (HIV/Leish), symptomatic visceral leishmaniasis (VL), recovered VL, DTH+ (Delayed-Type Hypersensitivity response - Leishmanin skin test), AIDS/VL, HIV+ only (HIV+), and Non-HIV/Non L. infantum infection (control subjects). Results: The cross-sectional study showed 24.2% of HIV+ subjects had positive anti-IgG Leishmania antibodies. After 3 years, 2.4% (8 of 333) of these HIV/Leish coinfected subjects developed AIDS/VL, whereas 1.05% (11 of 1,039) of HIV subjects with negative leishmania serology developed AIDS/VL. Poor adherence to antiretroviral therapy (p=0.0008) or prior opportunistic infections (p=0.0007) was associated with development of AIDS/VL. CD4+ (p=0.29) and CD8+ (p=0.38) T cells counts or viral load (p=0.34) were similar between asymptomatic HIV/Leish and HIV subjects. However, activated CD8+CD38+HLA-DR+ T cells were higher in asymptomatic HIV/Leish than HIV group. Likewise, senescent (CD57+) or exhausted (PD1+) CD8+ T cells were higher in asymptomatic HIV/Leish than in AIDS/VL or HIV groups. Conclusion: Although asymptomatic HIV/Leish subjects had normal and similar CD4+ and CD8+ T cells counts, their CD8+T cells had increased activation, senescence, and exhaustion, which could contribute to risk of developing VL.

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.

A molecular barcode and web-based data analysis tool to identify imported Plasmodium vivax malaria.

Traditionally, patient travel history has been used to distinguish imported from autochthonous malaria cases, but the dormant liver stages of Plasmodium vivax confound this approach. Molecular tools offer an alternative method to identify, and map imported cases. Using machine learning approaches incorporating hierarchical fixation index and decision tree analyses applied to 799 P. vivax genomes from 21 countries, we identified 33-SNP, 50-SNP and 55-SNP barcodes (GEO33, GEO50 and GEO55), with high capacity to predict the infection's country of origin. The Matthews correlation coefficient (MCC) for an existing, commonly applied 38-SNP barcode (BR38) exceeded 0.80 in 62% countries. The GEO panels outperformed BR38, with median MCCs > 0.80 in 90% countries at GEO33, and 95% at GEO50 and GEO55. An online, open-access, likelihood-based classifier framework was established to support data analysis (vivaxGEN-geo). The SNP selection and classifier methods can be readily amended for other use cases to support malaria control programs.

Geographical distribution and genetic diversity of Plasmodium vivax reticulocyte binding protein 1a correlates with patient antigenicity.

Plasmodium vivax is the most widespread cause of human malaria. Recent reports of drug resistant vivax malaria and the challenge of eradicating the dormant liver forms increase the importance of vaccine development against this relapsing disease. P. vivax reticulocyte binding protein 1a (PvRBP1a) is a potential vaccine candidate, which is involved in red cell tropism, a crucial step in the merozoite invasion of host reticulocytes. As part of the initial evaluation of the PvRBP1a vaccine candidate, we investigated its genetic diversity and antigenicity using geographically diverse clinical isolates. We analysed pvrbp1a genetic polymorphisms using 202 vivax clinical isolates from six countries. Pvrbp1a was separated into six regions based on specific domain features, sequence conserved/polymorphic regions, and the reticulocyte binding like (RBL) domains. In the fragmented gene sequence analysis, PvRBP1a region II (RII) and RIII (head and tail structure homolog, 152-625 aa.) showed extensive polymorphism caused by random point mutations. The haplotype network of these polymorphic regions was classified into three clusters that converged to independent populations. Antigenicity screening was performed using recombinant proteins PvRBP1a-N (157-560 aa.) and PvRBP1a-C (606-962 aa.), which contained head and tail structure region and sequence conserved region, respectively. Sensitivity against PvRBP1a-N (46.7%) was higher than PvRBP1a-C (17.8%). PvRBP1a-N was reported as a reticulocyte binding domain and this study identified a linear epitope with moderate antigenicity, thus an attractive domain for merozoite invasion-blocking vaccine development. However, our study highlights that a global PvRBP1a-based vaccine design needs to overcome several difficulties due to three distinct genotypes and low antigenicity levels.

Malaria protection due to sickle haemoglobin depends on parasite genotype.

Host genetic factors can confer resistance against malaria1, raising the question of whether this has led to evolutionary adaptation of parasite populations. Here we searched for association between candidate host and parasite genetic variants in 3,346 Gambian and Kenyan children with severe malaria caused by Plasmodium falciparum. We identified a strong association between sickle haemoglobin (HbS) in the host and three regions of the parasite genome, which is not explained by population structure or other covariates, and which is replicated in additional samples. The HbS-associated alleles include nonsynonymous variants in the gene for the acyl-CoA synthetase family member2-4 PfACS8 on chromosome 2, in a second region of chromosome 2, and in a region containing structural variation on chromosome 11. The alleles are in strong linkage disequilibrium and have frequencies that covary with the frequency of HbS across populations, in particular being much more common in Africa than other parts of the world. The estimated protective effect of HbS against severe malaria, as determined by comparison of cases with population controls, varies greatly according to the parasite genotype at these three loci. These findings open up a new avenue of enquiry into the biological and epidemiological significance of the HbS-associated polymorphisms in the parasite genome and the evolutionary forces that have led to their high frequency and strong linkage disequilibrium in African P. falciparum populations.

Disseminated Cutaneous Leishmaniasis and Alcohol Misuse, Northeast Brazil, 2015-2018.

Disseminated cutaneous leishmaniasis (DCL) is an uncommon form of Leishmania braziliensis infection. It remains unknown why some people develop this clinical condition. We describe 14 DCL patients in Northeast Brazil during 2015-2018. These patients regularly drank large amounts of alcohol, possibly increasing their risk for DCL.

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.

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

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

Intrinsic multiplication rate variation and plasticity of human blood stage malaria parasites.

Pathogen multiplication rate is theoretically an important determinant of virulence, although often poorly understood and difficult to measure accurately. We show intrinsic asexual blood stage multiplication rate variation of the major human malaria parasite Plasmodium falciparum to be associated with blood-stage infection intensity in patients. A panel of clinical isolates from a highly endemic West African population was analysed repeatedly during five months of continuous laboratory culture, showing a range of exponential multiplication rates at all timepoints tested, mean rates increasing over time. All isolates had different genome sequences, many containing within-isolate diversity that decreased over time in culture, but increases in multiplication rates were not primarily attributable to genomic selection. New mutants, including premature stop codons emerging in a few isolates, did not attain sufficiently high frequencies to substantially affect overall multiplication rates. Significantly, multiplication rate variation among the isolates at each of the assayed culture timepoints robustly correlated with parasite levels seen in patients at clinical presentation, indicating innate parasite control of multiplication rate that contributes to virulence.

Genetic diversity and neutral selection in Plasmodium vivax erythrocyte binding protein correlates with patient antigenicity.

Plasmodium vivax is the most widespread and difficult to treat cause of human malaria. The development of vaccines against the blood stages of P. vivax remains a key objective for the control and elimination of vivax malaria. Erythrocyte binding-like (EBL) protein family members such as Duffy binding protein (PvDBP) are of critical importance to erythrocyte invasion and have been the major target for vivax malaria vaccine development. In this study, we focus on another member of EBL protein family, P. vivax erythrocyte binding protein (PvEBP). PvEBP was first identified in Cambodian (C127) field isolates and has subsequently been showed its preferences for binding reticulocytes which is directly inhibited by antibodies. We analysed PvEBP sequence from 316 vivax clinical isolates from eight countries including China (n = 4), Ethiopia (n = 24), Malaysia (n = 53), Myanmar (n = 10), Papua New Guinea (n = 16), Republic of Korea (n = 10), Thailand (n = 174), and Vietnam (n = 25). PvEBP gene exhibited four different phenotypic clusters based on the insertion/deletion (indels) variation. PvEBP-RII (179-479 aa.) showed highest polymorphism similar to other EBL family proteins in various Plasmodium species. Whereas even though PvEBP-RIII-V (480-690 aa.) was the most conserved domain, that showed strong neutral selection pressure for gene purifying with significant population expansion. Antigenicity of both of PvEBP-RII (16.1%) and PvEBP-RIII-V (21.5%) domains were comparatively lower than other P. vivax antigen which expected antigens associated with merozoite invasion. Total IgG recognition level of PvEBP-RII was stronger than PvEBP-RIII-V domain, whereas total IgG inducing level was stronger in PvEBP-RIII-V domain. These results suggest that PvEBP-RII is mainly recognized by natural IgG for innate protection, whereas PvEBP-RIII-V stimulates IgG production activity by B-cell for acquired immunity. Overall, the low antigenicity of both regions in patients with vivax malaria likely reflects genetic polymorphism for strong positive selection in PvEBP-RII and purifying selection in PvEBP-RIII-V domain. These observations pose challenging questions to the selection of EBP and point out the importance of immune pressure and polymorphism required for inclusion of PvEBP as a vaccine candidate.

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.

Frequent expansion of Plasmodium vivax Duffy Binding Protein in Ethiopia and its epidemiological significance.

Plasmodium vivax invasion of human erythrocytes depends on the Duffy Binding Protein (PvDBP) which interacts with the Duffy antigen. PvDBP copy number has been recently shown to vary between P. vivax isolates in Sub-Saharan Africa. However, the extent of PvDBP copy number variation, the type of PvDBP multiplications, as well as its significance across broad samples are still unclear. We determined the prevalence and type of PvDBP duplications, as well as PvDBP copy number variation among 178 Ethiopian P. vivax isolates using a PCR-based diagnostic method, a novel quantitative real-time PCR assay and whole genome sequencing. For the 145 symptomatic samples, PvDBP duplications were detected in 95 isolates, of which 81 had the Cambodian and 14 Malagasy-type PvDBP duplications. PvDBP varied from 1 to >4 copies. Isolates with multiple PvDBP copies were found to be higher in symptomatic than asymptomatic infections. For the 33 asymptomatic samples, PvDBP was detected with two copies in two of the isolates, and both were the Cambodian-type PvDBP duplication. PvDBP copy number in Duffy-negative heterozygotes was not significantly different from that in Duffy-positives, providing no support for the hypothesis that increased copy number is a specific association with Duffy-negativity, although the number of Duffy-negatives was small and further sampling is required to test this association thoroughly.

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.

The origins and relatedness structure of mixed infections vary with local prevalence of P. falciparum malaria.

Individual malaria infections can carry multiple strains of Plasmodium falciparum with varying levels of relatedness. Yet, how local epidemiology affects the properties of such mixed infections remains unclear. Here, we develop an enhanced method for strain deconvolution from genome sequencing data, which estimates the number of strains, their proportions, identity-by-descent (IBD) profiles and individual haplotypes. Applying it to the Pf3k data set, we find that the rate of mixed infection varies from 29% to 63% across countries and that 51% of mixed infections involve more than two strains. Furthermore, we estimate that 47% of symptomatic dual infections contain sibling strains likely to have been co-transmitted from a single mosquito, and find evidence of mixed infections propagated over successive infection cycles. Finally, leveraging data from the Malaria Atlas Project, we find that prevalence correlates within Africa, but not Asia, with both the rate of mixed infection and the level of IBD.