We report a large multicenter genome-wide association study of Plasmodium falciparum resistance to artemisinin, the frontline antimalarial drug. Across 15 locations in Southeast Asia, we identified at least 20 mutations in kelch13 (PF3D7_1343700) affecting the encoded propeller and BTB/POZ domains, which were associated with a slow parasite clearance rate after treatment with artemisinin derivatives. Nonsynonymous polymorphisms in fd (ferredoxin), arps10 (apicoplast ribosomal protein S10), mdr2 (multidrug resistance protein 2) and crt (chloroquine resistance transporter) also showed strong associations with artemisinin resistance. Analysis of the fine structure of the parasite population showed that the fd, arps10, mdr2 and crt polymorphisms are markers of a genetic background on which kelch13 mutations are particularly likely to arise and that they correlate with the contemporary geographical boundaries and population frequencies of artemisinin resistance. These findings indicate that the risk of new resistance-causing mutations emerging is determined by specific predisposing genetic factors in the underlying parasite population.

Malaria elimination strategies require surveillance of the parasite population for genetic changes that demand a public health response, such as new forms of drug resistance. 1,2 Here we describe methods for large-scale analysis of genetic variation in Plasmodium falciparum by deep sequencing of parasite DNA obtained from the blood of patients with malaria, either directly or after short term culture. Analysis of 86,158 exonic SNPs that passed genotyping quality control in 227 samples from Africa, Asia and Oceania provides genome-wide estimates of allele frequency distribution, population structure and linkage disequilibrium. By comparing the genetic diversity of individual infections with that of the local parasite population, we derive a metric of within-host diversity that is related to the level of inbreeding in the population. An open-access web application has been established for exploration of regional differences in allele frequency and of highly differentiated loci in the P. falciparum genome.

We describe an analysis of genome variation in 825 Plasmodium falciparum samples from Asia and Africa that reveals an unusual pattern of parasite population structure at the epicentre of artemisinin resistance in western Cambodia. Within this relatively small geographical area we have discovered several distinct but apparently sympatric parasite subpopulations with extremely high levels of genetic differentiation. Of particular interest are three subpopulations, all associated with clinical resistance to artemisinin, which have skewed allele frequency spectra and remarkably high levels of haplotype homozygosity, indicative of founder effects and recent population expansion. We provide a catalogue of SNPs that show high levels of differentiation in the artemisinin-resistant subpopulations, including codon variants in various transporter proteins and DNA mismatch repair proteins. These data provide a population genetic framework for investigating the biological origins of artemisinin resistance and for defining molecular markers to assist its elimination.

Summary Background As the prevalence of artemisinin-resistant Plasmodium falciparum malaria increases in the Greater Mekong Subregion (GMS), emerging resistance to partner drugs in artemisinin combination therapies (ACTs) seriously threatens global efforts to treat and eliminate this disease. Molecular markers for ACT failure are urgently needed to monitor the spread of partner drug resistance, and to recommend alternative treatments in Southeast Asia and beyond. Methods We performed a genome-wide association study (GWAS) of 297 P. falciparum isolates from Cambodia to investigate the relationship of 11,630 exonic single-nucleotide polymorphisms (SNPs) and 43 copy number variations (CNVs) with in-vitro piperaquine 50% inhibitory concentrations (IC50s), and tested whether these genetic variants are markers of dihydroartemisinin-piperaquine failures. We then performed a survival analysis of 133 patients to determine whether candidate molecular markers predicted parasite recrudescence following dihydroartemisinin-piperaquine treatment. Findings Piperaquine IC50s increased significantly from 2011 to 2013 in 3 Cambodian provinces. Genome-wide analysis of SNPs identified a chromosome 13 region that associates with elevated piperaquine IC50s. A nonsynonymous SNP (encoding a Glu415Gly substitution) in this region, within a gene encoding an exonuclease, associates with parasite recrudescence following dihydroartemisinin-piperaquine treatment. Genome-wide analysis of CNVs revealed that a single copy of the mdr1 gene on chromosome 5 and a novel amplification of the plasmepsin II and plasmepsin III genes on chromosome 14 also associate with elevated piperaquine IC50s. After adjusting for covariates, both exo-E415G and plasmepsin II-III markers significantly associate with decreased treatment efficacy (0.38 and 0.41 survival rates, respectively). Interpretation The exo-E415G SNP and plasmepsin II-III amplification are markers of piperaquine resistance and dihydroartemisinin-piperaquine failures in Cambodia, and can help monitor the spread of these phenotypes into GMS countries, and elucidate the mechanism of piperaquine resistance. Since plasmepsins are involved in the parasite’s haemoglobin-to-haemozoin conversion pathway, targeted by related antimalarials, plasmepsin II-III amplification likely mediates piperaquine resistance. Funding Intramural Research Program of the US National Institute of Allergy and Infectious Diseases, National Institutes of Health; Wellcome Trust; Bill and Melinda Gates Foundation; Medical Research Council; and UK Department for International Development.

The widespread distribution and relapsing nature of Plasmodium vivax infection present major challenges for malaria elimination. To characterise the genetic diversity of this parasite within individual infections and across the population, we performed deep genome sequencing of >200 clinical samples collected across the Asia-Pacific region, and analysed data on >300,000 SNPs and 9 regions of the genome with large copy number variations. Individual infections showed complex patterns of genetic structure, with variation not only in the number of dominant clones but also in their level of relatedness and inbreeding. At the population level, we observed strong signals of recent evolutionary selection both in known drug resistance genes and at novel loci, and these varied markedly between geographical locations. These findings reveal a dynamic landscape of local evolutionary adaptation in P. vivax populations, and provide a foundation for genomic surveillance to guide effective strategies for control and elimination.

The current epidemic of artemisinin resistant Plasmodium falciparum in Southeast Asia is the result of a soft selective sweep involving at least 20 independent kelch13 mutations. In a large global survey, we find that kelch13 mutations which cause resistance in Southeast Asia are present at low frequency in Africa. We show that African kelch13 mutations have originated locally, and that kelch13 shows a normal variation pattern relative to other genes in Africa, whereas in Southeast Asia there is a great excess of non-synonymous mutations, many of which cause radical amino-acid changes. Thus, kelch13 is not currently undergoing strong selection in Africa, despite a deep reservoir of variations that could potentially allow resistance to emerge rapidly. The practical implications are that public health surveillance for artemisinin resistance should not rely on kelch13 data alone, and interventions to prevent resistance must account for local evolutionary conditions, shown by genomic epidemiology to differ greatly between geographical regions.DOI: http://dx.doi.org/10.7554/eLife.08714.001

SummaryBackgroundAntimalarial resistance is rapidly spreading across parts of southeast Asia where dihydroartemisinin–piperaquine is used as first-line treatment for Plasmodium falciparum malaria. The first published reports about resistance to antimalarial drugs came from western Cambodia in 2013. Here, we analyse genetic changes in the P falciparum population of western Cambodia in the 6 years before those reports.MethodsWe analysed genome sequence data on 1492 P falciparum samples from 11 locations across southeast Asia, including 464 samples collected in western Cambodia between 2007 and 2013. Different epidemiological origins of resistance were identified by haplotypic analysis of the kelch13 artemisinin resistance locus and the plasmepsin 2–3 piperaquine resistance locus.FindingsWe identified more than 30 independent origins of artemisinin resistance, of which the KEL1 lineage accounted for 140 (91%) of 154 parasites resistant to dihydroartemisinin–piperaquine. In 2008, KEL1 combined with PLA1, the major lineage associated with piperaquine resistance. By 2013, the KEL1/PLA1 co-lineage had reached a frequency of 63% (24/38) in western Cambodia and had spread to northern Cambodia.InterpretationThe KEL1/PLA1 co-lineage emerged in the same year that dihydroartemisinin–piperaquine became the first-line antimalarial drug in western Cambodia and spread rapidly thereafter, displacing other artemisinin-resistant parasite lineages. These findings have important implications for management of the global health risk associated with the current outbreak of multidrug-resistant malaria in southeast Asia.FundingWellcome Trust, Bill & Melinda Gates Foundation, Medical Research Council, UK Department for International Development, and the Intramural Research Program of the National Institute of Allergy and Infectious Diseases.

The increasing availability of large genomic data sets as well as the advent of Bayesian phylogenetics facilitates the investigation of phylogenetic incongruence, which can result in the impossibility of representing phylogenetic relationships using a single tree. While sometimes considered as a nuisance, phylogenetic incongruence can also reflect meaningful biological processes as well as relevant statistical uncertainty, both of which can yield valuable insights in evolutionary studies. We introduce a new tool for investigating phylogenetic incongruence through the exploration of phylogenetic tree landscapes. Our approach, implemented in the R package treespace, combines tree metrics and multivariate analysis to provide low‐dimensional representations of the topological variability in a set of trees, which can be used for identifying clusters of similar trees and group‐specific consensus phylogenies. treespace also provides a user‐friendly web interface for interactive data analysis and is integrated alongside existing standards for phylogenetics. It fills a gap in the current phylogenetics toolbox in R and will facilitate the investigation of phylogenetic results.