# Genomic epidemiology of artemisinin resistant malaria

Roberto Amato, Olivo Miotto, Charles J Woodrow, Jacob Almagro-Garcia, Ipsita Sinha, Susana Campino, Daniel Mead, Eleanor Drury, Mihir Kekre, Mandy Sanders, et al.

► To cite this version:

Roberto Amato, Olivo Miotto, Charles J Woodrow, Jacob Almagro-Garcia, Ipsita Sinha, et al. Genomic epidemiology of artemisinin resistant malaria. eLife, 2016, 48 (8), pp.959-964. 10.7554/eLife.08714. pasteur-01971950

HAL Id: pasteur-01971950

https://riip.hal.science/pasteur-01971950

Submitted on 7 Jan 2019

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
<figure><img src="image_2.png" /></figure>

# Genomic epidemiology of artemisinin resistant malaria

MalariaGEN Plasmodium falciparum Community Project*

**Abstract** 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: 10.7554/eLife.08714.001

*For correspondence: dominic@sanger.ac.uk (DPK); ra4@sanger.ac.uk (RA); olivo@tropmedres.ac (OM)

**Group author details:**

MalariaGEN Plasmodium falciparum Community Project See page 21

**Competing interest:** See page 23

**Funding:** See page 23

**Received:** 14 May 2015

**Accepted:** 23 December 2015

**Published:** 04 March 2016

**Reviewing editor:** Richard A Neher, Max Planck Institute for Developmental Biology, Germany

© This is an open-access article, free of all copyright, and may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose. The work is made available under the [Creative Commons CC0 public domain dedication](https://creativecommons.org/licenses/by-nc-nd/1.0/).

## Introduction

Artemisinin combination therapy (ACT), the frontline treatment for *P. falciparum* infection, has played a major part in reducing the number of deaths due to malaria over the past decade (*World Health Organization, 2014b*). However the increasing prevalence of artemisinin resistant *P. falciparum* across large parts of Southeast Asia threatens to destabilise malaria control worldwide (*Dondorp et al., 2009; Hien et al., 2012; Phyo et al., 2012; Kyaw et al., 2013; Ashley et al., 2014; World Health Organization, 2014a*). One of the main contemporary challenges in global health is to prevent artemisinin resistance from becoming established in Africa, where the consequences for childhood mortality could be disastrous (*Dondorp and Ringwald, 2013*).

Understanding the epidemiological and evolutionary driving forces behind artemisinin resistance is essential to develop effective strategies to stop it spreading. At the molecular level, artemisinin resistance is caused by mutations in a kelch protein encoded by PF3D7_1343700 on *P. falciparum* chromosome 13, referred to here as *kelch13*. More specifically, non-synonymous mutations in the *kelch13* propeller and BTB-POZ domains (KPBD) result in reduced sensitivity of *P. falciparum* to artemisinin, as demonstrated by multiple lines of evidence including laboratory studies of artificially acquired resistance, genetic association studies of natural resistance and allelic replacement experiments (*Ariey et al., 2014; Ghorbal et al., 2014; Miotto et al., 2015; Straimer et al., 2015; Takala-Harrison et al., 2015). Parasites with KBPD mutations tend to grow more slowly in the early part of the erythrocytic cycle, and have an enhanced unfolded protein response, both of which might act to protect against oxidative damage caused by artemisinin (*Dogovski et al., 2015; Mok et al., 2015*). It has also recently been shown that the *PI3K* protein is the target of artemisinin action, and that *kelch13* binds to *PI3K* to mark it for degradation; by affecting this binding, KBPD mutations result in high levels of *PI3K* which counteract the effect of artemisinin (*Mbengue et al., 2015*).

A striking characteristic of the current distribution of artemisinin resistance is that it is caused by multiple independent KPBD mutations emerging in different locations, i.e. it does not originate from a single mutational event. More than 20 KPBD SNPs have been associated with delayed parasite clearance during artemisinin treatment and there are several documented instances of the same
**eLife digest** Malaria is an infectious disease caused by a microscopic parasite called *Plasmodium*, which is transferred between humans by mosquitos. One species of malaria parasite called *Plasmodium falciparum* can cause particularly severe and life-threatening forms of the disease. Currently, the most widely used treatment for *P. falciparum* infections is artemisinin combination therapy, a treatment that combines the drug artemisinin (or a closely related molecule) with another antimalarial drug. However, resistance to artemisinin has started to spread throughout Southeast Asia.

Artemisinin resistance is caused by mutations in a parasite gene called *kelch13*, and researchers have identified over 20 different mutations in *P. falciparum* that confer artemisinin resistance. The diversity of mutations involved, and the fact that the same mutation can arise independently in different locations, make it difficult to track the spread of resistance using conventional molecular marker approaches.

Here, Amato, Miotto et al. sequenced the entire genomes of more than 3,000 clinical samples of *P. falciparum* from Southeast Asia and Africa, collected as part of a global network of research groups called the MalariaGEN *Plasmodium falciparum* Community Project. Amato, Miotto et al. found that African parasites had independently acquired many of the same *kelch13* mutations that are known to cause resistance to artemisinin in Southeast Asia. However the *kelch13* mutations seen in Africa remained at low levels in the parasite population, and appeared to be under much less pressure for evolutionary selection than those found in Southeast Asia.

These findings demonstrate that the emergence and spread of resistance to antimalarial drugs does not depend solely on the mutational process, but also on other factors that influence whether the mutations will spread in the population. Understanding how this is affected by different patterns of drug treatments and other environmental conditions will be important in developing more effective strategies for combating malaria.

DOI: 10.7554/eLife.08714.002

allele arising independently in different locations (*Ashley et al., 2014; Miotto et al., 2015; Takala-Harrison et al., 2015*). These are classic features of a soft selective sweep which, according to evolutionary theory, is most likely to arise in populations where the selected alleles are already present as standing genetic variations or have been repeatedly introduced by *de novo* mutations (*Hermisson and Pennings, 2005; Pennings and Hermisson, 2006; Messer and Petrov, 2013*). There is ongoing debate among evolutionary biologists about how commonly soft selective sweeps occur in nature (*Jensen, 2014; Garud et al., 2015*) but they have clearly played a role in previous forms of antimalarial drug resistance (*Nair et al., 2007; Salgueiro et al., 2010*) and the current epidemic of artemisinin resistance is the most extreme example of a soft selective sweep thus far observed in eukaryotes.

This creates a practical problem in monitoring the global spread of resistance. Artemisinin resistance can be measured directly, by following the rate of parasite clearance in patients (*Flegg et al., 2011*) or by testing parasite isolates in vitro (*Witkowski et al., 2013*), but these phenotypic assays are resource intensive and impractical for large-scale screening in resource-poor settings. Genetic approaches are therefore preferable for practical implementation of large-scale surveillance, but the soft selective sweep of artemisinin resistance produces much more heterogeneous genetic signatures than previous global waves of chloroquine and pyrimethamine resistance, where hard selective sweeps were the dominant mode of spread. Thus there is considerable uncertainty about the epidemiological significance of the growing number of non-synonymous KPBD mutations reported in Africa (*Ashley et al., 2014; Hopkins Sibley, 2015; Kamau et al., 2015; Taylor et al., 2015*). Previous studies in Africa have not identified variants known to cause resistance in Southeast Asia (*Kamau et al., 2015; Taylor et al., 2015*). At the same time, there are documented instances of African parasites that can be rapidly cleared by artemisinin treatment, although they carry KPBD mutations (*Ashley et al., 2014*). In the absence of comprehensive phenotypic data, it is not known which of these mutations, if any, are markers of resistance. This is a limitation of conventional molecular epidemiology, which tracks specific mutations and haplotypes and is poorly equipped to monitor
soft selective sweeps where new mutations are continually arising on different haplotypic backgrounds, making it difficult to keep track of their phenotypic effects and evolutionary trajectories.

Here we explore how genomic data might help overcome these practical obstacles to monitoring the current epidemic of artemisinin resistance. This analysis includes genome sequencing data for 3,411 clinical samples of *P. falciparum* obtained from 43 locations in 23 countries. This large dataset was generated by the MalariaGEN Plasmodium falciparum Community Project, (www.malariagen.net/projects/parasite/pf), a collaborative study in which multiple research groups working on different scientific questions are sharing genome variation data to generate an integrated view of polymorphism in the global parasite population. Data resources arising from the present analysis, including genotype calls and sample metadata, can be obtained at www.malariagen.net/resource/16. The MalariaGEN Plasmodium falciparum Community Project also developed a user-friendly web application, with interactive tools for exploring and querying the latest version of the data: www.malariagen.net/apps/pf.

## Results

### Africa and Southeast Asia both have many *kelch13* polymorphisms

Paired-end sequence reads were generated using the Illumina platform and aligned to the *P. falciparum* 3D7 reference genome, applying a series of quality control filters as previously described (see Materials and methods) (Manske et al., 2012; Miotto et al., 2013). The initial alignments identified 4,305,957 potential SNPs, which after quality control filtering produced a set of 935,601 exonic SNPs that could be genotyped with high confidence in the majority of samples, and that form the basis for the current analysis.

As summarised in Table 1, the dataset comprised 1,648 samples from Africa and 1,599 samples from Southeast Asia, allowing us to compare these two groups directly without the need for sample size corrections. We identified a total of 155 SNPs in the *kelch13* gene, of which 128 were seen in Africa and 62 in Southeast Asia (Table 2). Studies in Southeast Asia have found that artemisinin resistance is associated with non-synonymous polymorphisms in the KPBD. Out of a total of 46 non-synonymous SNPs in these domains we found a similar number in Africa (n = 26) and Southeast Asia (n = 34), with 14 seen in both places (Table 3). Seven of those observed in Africa have previously been associated with artemisinin resistance in Southeast Asia, including C580Y, the most common allele in resistant parasites (Miotto et al., 2015).

### *kelch13* polymorphisms in Africa appear to be indigenous

We asked whether *kelch13* polymorphisms seen in Africa had emerged independently, as opposed to migrating from Southeast Asia. We started by grouping samples according to genome-wide genetic similarity, based on a neighbour-joining (NJ) tree (Figure 1A). African and Southeast Asian samples formed two well-separated and distinct clusters, suggesting that gene flow between the two regions is very modest or negligible, a view supported by an alternative analysis using principal coordinate analysis (Figure 1B). None of the African parasites carrying *kelch13* mutations grouped with the Southeast Asian population, supporting the idea that these mutations are indigenous.

We considered that recombination between African and imported parasites could result in the transfer of DNA segments of Asian origin into genomes that otherwise appear to be local, and therefore not be easily detectable through genome-wide analysis. Thus, we analyzed the regions immediately flanking *kelch13* in both directions, using a probabilistic method to reconstruct the most likely geographical origin of haplotypes observed in African *kelch13* mutants (Figure 2). The vast majority of *kelch13* mutants showed no evidence that flanking regions may have been imported, which indicates that most mutations observed in Africa do not have a common origin with those in Asia, and are likely to have emerged independently on a variety of haplotypes (Figure 2—figure supplement 1).

Only 5 samples out of 56 were assigned a significant probability of Southeast Asian origin in at least one flanking sequence, and only two of these samples carried haplotypes similar to those in SEA mutants (Figure 2—figure supplement 2). Both in the NJ tree and in the PCA, these two African *kelch13* mutants do not cluster with the bulk of African samples, but occupy a somewhat isolated position (Figure 1). We note that they exhibit unusually high levels of heterozygosity (F<sub>WS</sub> < 0.4), with mixed calls randomly distributed across the genome, consistent with a mix of African and
<table><caption>**Table 1.** Count of samples analysed in this study.</caption><thead><tr><th>Region</th><th>Code</th><th>Sample counts</th><th>Country</th><th>Code</th><th>Sample counts</th></tr></thead><tbody><tr><td rowspan="7">West Africa</td><td rowspan="7">WAF</td><td rowspan="7">957</td><td>Burkina Faso</td><td>BF</td><td>56</td></tr><tr><td>Cameroon</td><td>CM</td><td>134</td></tr><tr><td>Ghana</td><td>GH</td><td>478</td></tr><tr><td>Gambia, The</td><td>GM</td><td>73</td></tr><tr><td>Guinea</td><td>GN</td><td>124</td></tr><tr><td>Mali</td><td>ML</td><td>87</td></tr><tr><td>Nigeria</td><td>NG</td><td>5</td></tr><tr><td rowspan="5">East Africa</td><td rowspan="5">EAF</td><td rowspan="5">412</td><td>Kenya</td><td>KE</td><td>52</td></tr><tr><td>Madagascar</td><td>MG</td><td>18</td></tr><tr><td>Malawi</td><td>MW</td><td>262</td></tr><tr><td>Tanzania</td><td>TZ</td><td>68</td></tr><tr><td>Uganda</td><td>UG</td><td>12</td></tr><tr><td>Central Africa</td><td>CAF</td><td>279</td><td>D. R. Congo</td><td>CD</td><td>279</td></tr><tr><td rowspan="2">South America</td><td rowspan="2">SAM</td><td rowspan="2">27</td><td>Colombia</td><td>CO</td><td>16</td></tr><tr><td>Peru</td><td>PE</td><td>11</td></tr><tr><td>South Asia</td><td>SAS</td><td>75</td><td>Bangladesh</td><td>BD</td><td>75</td></tr><tr><td rowspan="2">West Southeast Asia</td><td rowspan="2">WSEA</td><td rowspan="2">497</td><td>Myanmar</td><td>MM</td><td>111</td></tr><tr><td>Thailand</td><td>TH</td><td>386</td></tr><tr><td rowspan="3">East Southeast Asia</td><td rowspan="3">ESEA</td><td rowspan="3">1102</td><td>Cambodia</td><td>KH</td><td>762</td></tr><tr><td>Laos</td><td>LA</td><td>120</td></tr><tr><td>Vietnam</td><td>VN</td><td>220</td></tr><tr><td rowspan="2">Oceania</td><td rowspan="2">OCE</td><td rowspan="2">62</td><td>Indonesia (Papua)</td><td>ID</td><td>17</td></tr><tr><td>Papua New Guinea</td><td>PG</td><td>45</td></tr></tbody></table>

DOI: 10.7554/eLife.08714.003

SEA parasites. Continuous genetic monitoring of the parasite population will determine whether these are indeed just isolated cases, or they constitute very early evidence of gene flow between the two regions. Since these samples passed through a number of different laboratories, we cannot absolutely rule out the possibility that these mixtures could be the result of biological contamination during sample preparation and processing.

## Across the genome there are many more rare variants in Africa

Historical demographic changes such as population expansions and bottlenecks (Tanabe et al., 2010), and epidemiological and environmental factors (Prugnolle et al., 2010) are highly influential

<table><caption>**Table 2.** Worldwide distribution of *kelch13* mutations. Number of *kelch13* propeller and BTB-POZ domain (KPBD) mutations present (not necessarily exclusively) in 5 populations (AFR = Africa, SEA = Southeast Asia, SAS = South Asia, OCE = Oceania, SAM = South America). Sample size is reported for each population.</caption><thead><tr><th></th><th></th><th>AFR<br>(N = 1,648)</th><th>SEA<br>(N = 1,599)</th><th>SAS<br>(N = 75)</th><th>OCE<br>(N = 62)</th><th>SAM<br>(N = 27)</th></tr></thead><tbody><tr><td rowspan="2">Non-synonymous</td><td>KPBD</td><td>26</td><td>34</td><td>1</td><td>1</td><td>0</td></tr><tr><td>Upstream region</td><td>42</td><td>16</td><td>2</td><td>1</td><td>1</td></tr><tr><td rowspan="2">Synonymous</td><td>KPBD</td><td>38</td><td>9</td><td>1</td><td>0</td><td>0</td></tr><tr><td>Upstream region</td><td>22</td><td>3</td><td>1</td><td>0</td><td>0</td></tr></tbody></table>

DOI: 10.7554/eLife.08714.004
<table><caption>**Table 3.** List of non-synonymous KPBD mutations. Non-synonymous mutations found in the *kelch13* propeller and BTB-POZ domains (KPBD) and their position on chromosome Pf3D7_13_v3. For each mutation is reported where it has been observed and in how many samples. Where known, we reported if the mutation has been previously observed in patients with a prolonged parasite clearance half-life (>5 hr) by Miotto et al. 2014 and/or Ashley et al 2014. Sample size is reported for each population.</caption><thead><tr><th>Mutation</th><th>Genomic coordinates</th><th>AFR (N = 1,648)</th><th>SEA (N = 1,599)</th><th>SAS (N = 75)</th><th>OCE (N = 62)</th><th>SAM (N = 27)</th><th>Observed in ART-R samples?</th></tr></thead><tbody><tr><td>D353Y</td><td>1725941</td><td>-</td><td>4</td><td>-</td><td>-</td><td>-</td><td>Yes</td></tr><tr><td>F395Y</td><td>1725814</td><td>-</td><td>1</td><td>-</td><td>-</td><td>-</td><td>No</td></tr><tr><td>I416V</td><td>1725752</td><td>1</td><td>1</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>I416M</td><td>1725750</td><td>1</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>K438N</td><td>1725684</td><td>-</td><td>1</td><td>-</td><td>-</td><td>-</td><td>No</td></tr><tr><td>P441L</td><td>1725676</td><td>-</td><td>27</td><td>-</td><td>-</td><td>-</td><td>Yes</td></tr><tr><td>P443S</td><td>1725671</td><td>-</td><td>1</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>F446I</td><td>1725662</td><td>-</td><td>7</td><td>-</td><td>-</td><td>-</td><td>Yes</td></tr><tr><td>G449A</td><td>1725652</td><td>-</td><td>7</td><td>-</td><td>-</td><td>-</td><td>Yes</td></tr><tr><td>S459L</td><td>1725622</td><td>2</td><td>2</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>A481V</td><td>1725556</td><td>-</td><td>4</td><td>-</td><td>-</td><td>-</td><td>Yes</td></tr><tr><td>S485N</td><td>1725544</td><td>-</td><td>1</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Y493H</td><td>1725521</td><td>1</td><td>76</td><td>-</td><td>-</td><td>-</td><td>Yes</td></tr><tr><td>V520I</td><td>1725440</td><td>1</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>S522C</td><td>1725434</td><td>2</td><td>1</td><td>-</td><td>-</td><td>-</td><td>Yes</td></tr><tr><td>P527H</td><td>1725418</td><td>1</td><td>5</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>C532S</td><td>1725404</td><td>1</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>V534L</td><td>1725398</td><td>2</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>N537I</td><td>1725388</td><td>1</td><td>1</td><td>-</td><td>-</td><td>-</td><td>No</td></tr><tr><td>G538V</td><td>1725385</td><td>-</td><td>19</td><td>-</td><td>-</td><td>-</td><td>Yes</td></tr><tr><td>R539T</td><td>1725382</td><td>-</td><td>63</td><td>-</td><td>-</td><td>-</td><td>Yes</td></tr><tr><td>I543T</td><td>1725370</td><td>-</td><td>34</td><td>-</td><td>-</td><td>-</td><td>Yes</td></tr><tr><td>P553L</td><td>1725340</td><td>2</td><td>24</td><td>-</td><td>-</td><td>-</td><td>Yes</td></tr><tr><td>A557S</td><td>1725329</td><td>1</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>R561H</td><td>1725316</td><td>1</td><td>24</td><td>-</td><td>-</td><td>-</td><td>Yes</td></tr><tr><td>V568G</td><td>1725295</td><td>-</td><td>6</td><td>-</td><td>-</td><td>-</td><td>Yes</td></tr><tr><td>T573S</td><td>1725280</td><td>2</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>P574L</td><td>1725277</td><td>-</td><td>12</td><td>-</td><td>-</td><td>-</td><td>Yes</td></tr><tr><td>R575K</td><td>1725274</td><td>-</td><td>3</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>A578S</td><td>1725266</td><td>18</td><td>-</td><td>-</td><td>-</td><td>-</td><td>No</td></tr><tr><td>C580Y</td><td>1725259</td><td>2</td><td>423</td><td>-</td><td>1</td><td>-</td><td>Yes</td></tr><tr><td>D584V</td><td>1725247</td><td>-</td><td>3</td><td>-</td><td>-</td><td>-</td><td>Yes</td></tr><tr><td>V589I</td><td>1725233</td><td>2</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>T593S</td><td>1725221</td><td>1</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>E612D</td><td>1725162</td><td>1</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Q613E</td><td>1725161</td><td>5</td><td>1</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Q613L</td><td>1725160</td><td>1</td><td>-</td><td>-</td><td>-</td><td>-</td><td>No</td></tr><tr><td>F614L</td><td>1725158</td><td>-</td><td>1</td><td>-</td><td>-</td><td>-</td><td>No</td></tr><tr><td>Y630F</td><td>1725109</td><td>2</td><td>1</td><td>-</td><td>-</td><td>-</td><td>-</td></tr></tbody></table>

Table 3 continued on next page
<table><caption>Table 3 continued</caption><thead><tr><th>Mutation</th><th>Genomic coordinates</th><th>AFR (N = 1,648)</th><th>SEA (N = 1,599)</th><th>SAS (N = 75)</th><th>OCE (N = 62)</th><th>SAM (N = 27)</th><th>Observed in ART-R samples?</th></tr></thead><tbody><tr><td>V637I</td><td>1725089</td><td>2</td><td>-</td><td>-</td><td>-</td><td>-</td><td></td></tr><tr><td>P667A</td><td>1724999</td><td>-</td><td>2</td><td>-</td><td>-</td><td>-</td><td></td></tr><tr><td>P667L</td><td>1724998</td><td>-</td><td>2</td><td>1</td><td>-</td><td>-</td><td></td></tr><tr><td>F673I</td><td>1724981</td><td>-</td><td>3</td><td>-</td><td>-</td><td>-</td><td>Yes</td></tr><tr><td>A675V</td><td>1724974</td><td>1</td><td>18</td><td>-</td><td>-</td><td>-</td><td>Yes</td></tr><tr><td>A676S</td><td>1724972</td><td>2</td><td>3</td><td>-</td><td>-</td><td>-</td><td></td></tr><tr><td>H719N</td><td>1724843</td><td>1</td><td>8</td><td>-</td><td>-</td><td>-</td><td>Yes</td></tr></tbody></table>

DOI: 10.7554/eLife.08714.005

forces that shaped the allele frequency spectra of *P. falciparum* populations across the globe (Nielsen et al., 2009; Manske et al., 2012). In order to properly contextualize the numbers and frequencies of *kelch13* mutations, it is therefore important to characterize genomic variation patterns in different geographical regions.

One of the most striking features of this dataset is the high number of rare variations in the high-quality SNP list. At more than half of all polymorphic sites, the minor allele was only carried by a single sample (referred to as *singletons*, n = 330,783 or 35%) or by two samples (*doubletons*, n = 214,179 or 23%), often in heterozygous calls. By contrast, only 13,383 polymorphisms (1.4%) had a minor allele in ≥5% of samples. Rare alleles, however, are not evenly distributed geographically. There is a large excess of polymorphisms with minor allele frequency (MAF) below 0.1% in Africa (71% of all SNPs, vs. 17% in SEA), while numbers in the two regions are similar for SNPs with MAF>1% (2% of all SNPs, Figure 3a). Rare variations in Africa are not confined to a limited set of highly variable genes, but evenly distributed across the genome, as attested by the distribution of variants across all genes: SNP density in Africa (median = 67 SNPs/kbp, interquartile range = 51–84) is approximately 3.9 times higher than in SEA (median = 17, IQR = 13–22, p < 10⁻¹⁶, Figure 3b). Very similar ratios are estimated in both non-synonymous (Africa median SNP density = 43, IQR = 27–58; SEA median = 11, IQR = 7–15) and synonymous variants (Africa median SNP density=25, IQR = 20–30, SEA median = 6, IQR = 4–8). Accordingly, we found virtually identical distributions of the ratio of non-synonymous to synonymous mutations (N/S ratio) in the two regions (Figure 3c). This suggests that the huge disparity in SNP density between the two regions is more likely to be the result of different demographic histories and epidemiological characteristics, such as changes in effective population size (Joy et al., 2003), rather than the product of different selective constraints.

In summary, we observe many more rare variants in Africa than in SEA; however we expect N/S ratios to be similar in these two regions in genes that are not subjected to selective pressures.

## Comparing *kelch13* with other parts of the genome

The density of *kelch13* synonymous variations in the two continents is roughly consistent with that observed in the rest of the genome (Africa: 28 SNPs/kbp, SEA: 6; Figure 4a), which is expected since synonymous changes are less likely to be affected by selection. The excess of African non-synonymous mutations in the upstream region is also consistent with expectations (Africa: 44 SNPs/kbp, SEA: 16). In contrast, non-synonymous polymorphisms in the KPBD show a reversal of this relationship: SEA parasites possess about 30% more polymorphisms than African ones (34 in SEA vs. 26 in AFR). In addition, all but two non-synonymous changes in the KPBD are observed in Africa at very low frequency (singletons and doubletons), while in SEA more than half of the changes are observed in >2 samples (Table 4 and Figure 3—figure supplement 1).

At a first approximation, these observations are consistent with a high number of non-synonymous changes that have risen in frequency in SEA parasites because of their association with artemisinin resistance, and a low number of mostly rare alleles in Africa, where artemisinin has been introduced more recently and resistance is yet to be reported.
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Differentiation between African and Asian genomes. (A) Neighbour-joining tree showing population structure across all sampling sites, with sample branches coloured according to country groupings (Table 1). Branches with large open tip circle indicate the sample is a *kelch13* mutant, while those with a small black symbol are mixed infections (i.e. mixture of wild-type and mutant parasites or two mutant parasites with different mutations). Branches without tip symbols are *kelch13* wild type. African *kelch13* mutants are, at a genomic level, similar to other African samples. (B) Plot of second principal component (PC2) against the first (PC1), computed from a principal coordinate analysis (PCoA) of all samples in the present study, based on the same pairwise genetic distance matrix used for the tree of Figure 1A. PC1 clearly separates African samples from those collected in SEA, while PC2 is mainly driven by extreme population structure in ESEA.

DOI: 10.7554/eLife.08714.006</figcaption></figure>
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Local origin of African *kelch13* mutations. (A) Chromosome painting (see Materials and methods) of the 52 African *kelch13* mutants across the two 250 kbp flanking regions on each side of the *kelch13* gene. Each genome chunk is coloured according to the aggregated posterior probabilities that it originated in the African (red) or SEA (blue) population, according to the scale shown. (B) Detail of the flanking regions over a span of approximately 15 kbp, using the same colour scheme. The country of origin is indicated on the left, followed by the proportion of African chunks identified; the *kelch13* mutation carried by the sample is shown on the right. Samples are sorted by the proportion of Asian chunks in this window, and the same order was applied to Panel A. Only five samples (lower region of the panel) show strong probability of Asian origin of the chunks closest to *kelch13*.</figcaption></figure>

DOI: 10.7554/eLife.08714.007

The following figure supplements are available for figure 2:

**Figure supplement 1.** Diversity within *kelch13*-flanking haplotypes in African mutants.

DOI: 10.7554/eLife.08714.008

**Figure supplement 2.** Origin of flanking haplotypes in selected African *kelch13* mutants.

DOI: 10.7554/eLife.08714.009
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Number and density of variants in Africa and Southeast Asia. (a) Allele frequency spectrum for Africa (red) and SEA (blue). Polymorphisms were binned by their minor allele frequency (MAF), and the counts in each bin were plotted against frequency, shown on a logarithmic scale. Although the number of high-frequency variations is consistent between the two regions, samples from Africa possess an excess of low-frequency variations. (b) SNP density per gene: for each gene, the number of variants found in the two regions is normalized by the length of the coding region (in kbp). African samples have on average 3.9 times more mutations than parasites from SEA. (c) Non-synonymous/synonymous ratio per gene: ratios of non-synonymous to synonymous mutations found per gene are similar in the two regions. To reduce artifacts due to small numbers, only genes with at least 10 SNP were considered in both analyses.</figcaption></figure>

DOI: 10.7554/eLife.08714.010

The following figure supplement is available for figure 3:

**Figure supplement 1.** Site frequency spectrum for non-synonymous mutations in the KPBD.

DOI: 10.7554/eLife.08714.011
Comparing *kelch13* to other highly conserved genes

Although the function of *kelch13* is as yet unclear, an alignment of its homologous gene sequences in eight *Plasmodium* species shows that the propeller and BTB-POZ domains are part of a highly conserved region (Figure 4b), suggesting a crucial role in parasite fitness. A reconstruction of ancestral alleles from this alignment suggests that *P. falciparum* accumulated only five conservative amino acid changes in the *kelch13 propeller* domain since diverging from other species 55 Myr ago (Escalante and Ayala, 1995) (Table 5). Given this extreme level of conservation, non-synonymous polymorphisms may appear surprisingly numerous in the present dataset, both in SEA (n = 34) and in Africa (n = 26). Such elevated numbers may be produced by selection processes; alternatively, they may be present in a large neutrally-evolving population, in which low-frequency variations continually emerge, but can only be detected for a brief span of time before they are removed by genetic drift and/or purifying selection. The question is, then, whether neutral evolution can account for the pattern of *kelch13* mutations observed here.

To answer this question, we compared patterns of *kelch13* mutations to those in the rest of the *P. falciparum* genome Supplementary file 1. Since fewer non-synonymous mutations are expected in more conserved genes, we applied *genomic calibration*, i.e. we stratified these analyses by evolutionary conservation. Each gene was assigned a *conservation score* determined from a sequence alignment of the *P. falciparum* gene with its *P. chabaudi* homologue, using a substitution matrix corrected for the AT bias in the *Pf* genome (Brick and Pizzi, 2008). *P. chabaudi* was chosen since it was the member of the group most differentiated from *P. falciparum* (rodent plasmodia) with the most complete reference sequence. A genome-wide non-linear negative correlation between gene conservation and

<figure><img src="image_6.png" /><figcaption>**Figure 4.** Structure of *kelch13*, positioning of mutations in Africa and Southeast Asia, and sequence conservation. (a) The amino acid positions of *kelch13* polymorphisms observed in Africa (red) and SEA (blue) are shown. Coloured rectangles describe the extents of the resistance domains (BTB/POZ: aa. 349–448; kelch propeller: aa. 443–721) and upstream region, with the locations of non-synonymous changes indicated above, and that of synonymous changes below. Short lines represent singleton/doubleton mutations, while longer lines represent more frequent mutations. (b) Conservation score across amino acid residues of *kelch13*, derived by applying a CCF53P62 matrix on alignments of the *P. falciparum* gene coding sequence with its homologues in six other *Plasmodium* species for which high-quality sequence data were available: *P. reichenowi*, *P. vivax*, *P. knowlesi*, *P. yoelii*, *P. berghei*, and *P. chabaudi* (see Materials and methods). Although the region below position 200 is less conserved, particularly in rodent species (*P. yoelii*, *berghei*, and *chabaudi*), there is remarkably high conservation across all species over the rest of the gene, which includes the KPBD.

DOI: 10.7554/eLife.08714.012</figcaption></figure>
<table><caption>**Table 4.** Frequency of the non-synonymous KPBD mutations. Counts of non-synonymous mutations in the conserved propeller and BTB-POZ domains of *kelch13* are shown for each geographical region, stratified by the number of samples in which they are observed. Sample size for each population is reported.</caption><thead><tr><th></th><th>AFR<br>(N = 1,648)</th><th>SEA<br>(N = 1,599)</th><th>SAS<br>(N = 75)</th><th>OCE<br>(N = 62)</th><th>SAM<br>(N = 27)</th></tr></thead><tbody><tr><td>1–2 samples</td><td>24</td><td>13</td><td>1</td><td>1</td><td>0</td></tr><tr><td>3–5 samples</td><td>1</td><td>7</td><td>0</td><td>0</td><td>0</td></tr><tr><td>&gt;5 samples</td><td>1</td><td>14</td><td>0</td><td>0</td><td>0</td></tr></tbody></table>

DOI: 10.7554/eLife.08714.013

N/S ratio is clearly observable; this trend is almost identical in the two populations (Figure 5a). Although *kelch13* did not diverge significantly from this relationship in Africa (P = 0.2), its N/S ratio in SEA was the highest observed at its level of conservation, far exceeding the expected ratio (P<0.001). Accordingly, *kelch13* showed the most significant difference in N/S ratios between Africa and SEA genome-wide (3.7-fold, P = 2 × 10⁻⁴ by Fisher's exact test, Figure 5b and Supplementary file 1), even compared to other well-known drug resistance genes (Figure 5—figure supplement 1). Such unusually high N/S ratio in SEA parasites is mainly due to an excess of high frequency non-synonymous variations (Figure 5—figure supplement 2), suggesting that multiple independent origins of artemisinin resistance (Miotto et al., 2015; Takala-Harrison et al., 2015) have produced an unusually large number of common non-synonymous mutations.

From this analysis we conclude that the high prevalence of *kelch13* non-synonymous variants in SEA is not explainable by neutral evolution, but is consistent with selection of artemisinin-resistance alleles. In Africa, on the other hand, the observed non-synonymous changes appear to constitute a "physiological" level of variation consistent with a population rich in low-frequency alleles.

## In Southeast Asia there are more radical substitutions in *kelch13*

The different *kelch13* mutation repertoires in Africa and in SEA raise the question of whether these sets of mutations have different structural and functional properties. While there is high conservation across the whole of the propeller domain, it is unlikely that all possible amino acid changes have the same functional relevance or that they all carry the same fitness cost for parasites. Although direct measures of functional relevance are not yet available, and the exact function of *kelch13* is hitherto unknown, we can make statistical comparisons of some properties of the observed changes, in at least two respects. First, assuming that *kelch13* function is conserved across *Plasmodium* species, we can assess the strength of evolutionary constraints at any given position by examining whether amino acid substitutions between species are conservative or radical. Second, given that *kelch* proteins have been shown in other organisms to play an adapter role, with key binding sites defined by the arrangement of hydrophobic β-strands in the propeller domain (Adams et al., 2000), we can assess changes in hydrophobicity caused by the observed mutations, which may be informative of their functional importance.

Detailed mapping against the secondary structure of the propeller domain suggests that the polymorphisms found in SEA parasites occur in different blades, preferentially at positions proximal to the first and second β-strand of the propeller's blades (Figure 6). This may indicate that these two strands play a role in defining the binding interface to the PI3K protein involved in artemisinin resistance (Mbengue et al., 2015), but this needs to be confirmed by in-depth structural analyses.

We characterized changes in the propeller domain by a conservation score derived from a substitution matrix specific to AT-rich genomes (Brick and Pizzi, 2008), and assigned a hydrophobicity score to each site, estimated from the Kyte-Doolittle (KD) hydropathicity score (Kyte and Doolittle, 1982). The five putative derived alleles that have become established in the *P. falciparum* *kelch13* propeller domain since its divergence from other *Plasmodia* are all conservative changes at hydrophilic sites (Figure 6—figure supplement 1a). Common mutations in Africa have characteristics broadly consistent with this conservative history of change (Figure 6—figure supplement 1b and Table 6). Polymorphisms in SEA parasites, on the other hand, show a pattern of changes that are
<table><caption>**Table 5.** *Kelch13 propeller* domain mutations in different *Plasmodium* species. Here we report amino acid allele differences in a multiple sequence alignments of *kelch13* homologues for seven species of *Plasmodium* parasites for which high-quality sequence data were available: *P. falciparum* (Pf), *P. reichenowi* (Pr), *P. vivax* (Pv), *P. knowlesi* (Pk), *P. yoelii* (Py), *P. berghei* (Pb), and *P. chabaudi* (Pc). The species formed three groups by similarity: Laverania (Pf, Pr), primate Plasmodia (Pv, Pk) and rodent Plasmodia (Py, Pb and Pc). An allele shared by all members of two different groups was identified as a putative ancestral allele. The table shows, for each position where at least one species exhibits a difference from the others: the amino acid position in the *Pf kelch13* sequence; the putative ancestral amino acid allele; the alleles in the various species (columns with heading listing multiple species show mutations common to those species); and a substitution score of the mutation, based on a CCF53P62 substitution matrix (see Materials and methods). All substitution scores are ≥0, denoting conservative substitutions.</caption><thead><tr><th>Pf Position</th><th>Ancestral Allele</th><th>Pf,Pr</th><th>Pv,Pk</th><th>Pk</th><th>Py,Pb,Pc</th><th>Py,Pb</th><th>Pb</th><th>Pc</th><th>Substitution Score</th></tr></thead><tbody><tr><td>434</td><td>F</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>Y</td><td>2</td></tr><tr><td>447</td><td>C</td><td>-</td><td>-</td><td>-</td><td>-</td><td>S</td><td>-</td><td>-</td><td>1</td></tr><tr><td>448</td><td>I</td><td>-</td><td>M</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>2</td></tr><tr><td>517</td><td>T</td><td>V</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>0</td></tr><tr><td>519</td><td>F</td><td>Y</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>2</td></tr><tr><td>520</td><td>V</td><td>-</td><td>-</td><td>-</td><td>L</td><td>-</td><td>-</td><td>-</td><td>0</td></tr><tr><td>534</td><td>V</td><td>-</td><td>-</td><td>-</td><td>I</td><td>-</td><td>-</td><td>-</td><td>2</td></tr><tr><td>550</td><td>S</td><td>-</td><td>-</td><td>C</td><td>-</td><td>-</td><td>-</td><td>-</td><td>1</td></tr><tr><td>566</td><td>V</td><td>-</td><td>-</td><td>-</td><td>I</td><td>-</td><td>-</td><td>-</td><td>2</td></tr><tr><td>568</td><td>V</td><td>-</td><td>I</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>2</td></tr><tr><td>578</td><td>A</td><td>-</td><td>S</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>0</td></tr><tr><td>584</td><td>D</td><td>-</td><td>-</td><td>E</td><td>-</td><td>-</td><td>-</td><td>-</td><td>1</td></tr><tr><td>590</td><td>I</td><td>-</td><td>-</td><td>-</td><td>-</td><td>V</td><td>-</td><td>-</td><td>2</td></tr><tr><td>593</td><td>T</td><td>-</td><td>-</td><td>-</td><td>A</td><td>-</td><td>-</td><td>-</td><td>0</td></tr><tr><td>605</td><td>D</td><td>E</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>1</td></tr><tr><td>613</td><td>Q</td><td>-</td><td>-</td><td>-</td><td>-</td><td>N</td><td>-</td><td>K</td><td>0</td></tr><tr><td>648</td><td>D</td><td>-</td><td>-</td><td>-</td><td>E</td><td>-</td><td>-</td><td>-</td><td>1</td></tr><tr><td>666</td><td>V</td><td>-</td><td>-</td><td>-</td><td>I</td><td>-</td><td>-</td><td>-</td><td>2</td></tr><tr><td>676</td><td>A</td><td>-</td><td>-</td><td>-</td><td>T</td><td>-</td><td>-</td><td>-</td><td>0</td></tr><tr><td>691</td><td>D</td><td>E</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>1</td></tr><tr><td>708</td><td>I</td><td>L</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>0</td></tr><tr><td>711</td><td>S</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>P</td><td>-</td><td>0</td></tr><tr><td>723</td><td>I</td><td>-</td><td>-</td><td>-</td><td>V</td><td>-</td><td>-</td><td>-</td><td>2</td></tr></tbody></table>

DOI: 10.7554/eLife.08714.014

more radical than those in Africa (P = 10⁻³) and more commonly found at hydrophobic sites (Figure 6—figure supplement 1c and Table 6).

Taken together, the above results suggest that the propeller domain has long been under very strong evolutionary constraints, and that the number and nature of the changes observed in African parasites is consistent with these constraints, once we discard the abundant rare alleles expected in such a large population. In contrast, mutations in SEA parasites are not only far more numerous than expected, but they produce radical changes that are likely to be important determinants of the binding properties of the *kelch13* protein, consistent with recent findings that binding of *kelch13* to the *PI3K* protein is a critical factor in *P. falciparum* response to artemisinin (Mbengue et al., 2015).

## Genetic background

A recent study has shown that resistance-causing KPBD mutations are significantly more likely to arise in parasites with a particular genetic background (Miotto et al., 2015). This predisposing genetic background is marked by specific SNP alleles of the genes encoding ferredoxin (fd), apicoplast ribosomal protein S10 (arps10), multidrug resistance protein 2 (mdr2) and chloroquine
<figure><img src="image_7.png" /><figcaption>**Figure 5.** Genome-wide analysis of N/S ratio. (a) For each gene with more than 2 synonymous or non-synonymous SNPs, the N/S ratio in Africa (red points) and in SEA (blue points) are plotted against the conservation score of the gene coding sequence. The *kelch13* gene values are represented by larger circles. For each region, a solid line shows median values, while dotted lines delimit 95% of the genes at varying levels of conservation. This plot is truncated on the y-axis to show more clearly the bulk of the distribution; the full range is shown in **Figure 5—figure supplement 2**. (b) Histogram showing the distribution of the ratio of N/S ratios in SEA and Africa, for all genes with ≥5 synonymous and ≥5 non-synonymous SNPs on each region. An arrow shows the placement of *kelch13*.</figcaption></figure>

DOI: 10.7554/eLife.08714.015

The following figure supplements are available for figure 5:

**Figure supplement 1.** N/S ratio of well-known drug resistance genes.

DOI: 10.7554/eLife.08714.016

**Figure supplement 2.** Genome-wide analysis of high-frequency SNP density.

DOI: 10.7554/eLife.08714.017
resistance transporter (crt). Here we extend this analysis, confirming that this particular combination of variants is extremely common in the parts of Southeast Asia where artemisinin resistance is known to be established, and is absent from Africa and other regions sampled here (Table 7).

Discussion

This study demonstrates the value of genomic data in characterising the current epidemic of artemisinin resistance, which is problematic for conventional molecular epidemiology since new resistance-causing mutations are continually emerging on different haplotypic backgrounds. A key problem is to define the geographical origin of KBPD mutations, and we show that this can be solved by using genomic epidemiological data to analyse ancestral relationships between samples, thereby demonstrating that the KPBD mutations observed in Africa are of local origin.

<figure><img src="image_8.png" /><figcaption>**Figure 6.** Structure of the *kelch13* propeller domain, showing the position of mutations in Southeast Asia and Africa. The sequence of the *kelch13* propeller domain (amino acids 443–726) is organized according to its 6-blade tertiary structure, with the four β-strands characterizing each blade highlighted in colour. Polymorphisms observed in SEA (top panel) and Africa (bottom panel) are shown by circles above the mutated position. Small circles indicate very rare mutations (singletons and doubletons), while larger circles are used for more frequent mutations.</figcaption></figure>

DOI: 10.7554/eLife.08714.018

The following figure supplement is available for figure 6:

<caption>**Figure supplement 1.** Characterization of *kelch13* mutations in Africa and Southeast Asia.</caption>

DOI: 10.7554/eLife.08714.019
<table><caption>**Table 6.** Conservation score of KPBD mutations. The table shows, for each non-synonymous KPBD mutation observed in the dataset, the number of samples carrying the mutation in Africa (AFR), in Southeast Asia (SEA), and a substitution score of the mutation, based on a CCF53P62 substitution matrix; lower values indicate more radical substitutions. Mutations observed in Africa tend to have higher conservation score, whereas in SEA mutations tend to be more radical.</caption><thead><tr><th>Mutation</th><th>AFR</th><th>SEA</th><th>CFF53P62</th></tr></thead><tbody><tr><td>Q613E</td><td>5</td><td>1</td><td>2</td></tr><tr><td>Y630F</td><td>2</td><td>1</td><td>2</td></tr><tr><td>V637I</td><td>2</td><td>0</td><td>2</td></tr><tr><td>V589I</td><td>2</td><td>0</td><td>2</td></tr><tr><td>T573S</td><td>2</td><td>0</td><td>2</td></tr><tr><td>Y493H</td><td>1</td><td>76</td><td>2</td></tr><tr><td>I416V</td><td>1</td><td>1</td><td>2</td></tr><tr><td>T593S</td><td>1</td><td>0</td><td>2</td></tr><tr><td>V520I</td><td>1</td><td>0</td><td>2</td></tr><tr><td>I416M</td><td>1</td><td>0</td><td>2</td></tr><tr><td>F395Y</td><td>0</td><td>1</td><td>2</td></tr><tr><td>S522C</td><td>2</td><td>1</td><td>1</td></tr><tr><td>E612D</td><td>1</td><td>0</td><td>1</td></tr><tr><td>C532S</td><td>1</td><td>0</td><td>1</td></tr><tr><td>R575K</td><td>0</td><td>3</td><td>1</td></tr><tr><td>A578S</td><td>18</td><td>0</td><td>0</td></tr><tr><td>A676S</td><td>2</td><td>3</td><td>0</td></tr><tr><td>V534L</td><td>2</td><td>0</td><td>0</td></tr><tr><td>R561H</td><td>1</td><td>24</td><td>0</td></tr><tr><td>A675V</td><td>1</td><td>18</td><td>0</td></tr><tr><td>H719N</td><td>1</td><td>8</td><td>0</td></tr><tr><td>P527H</td><td>1</td><td>5</td><td>0</td></tr><tr><td>A557S</td><td>1</td><td>0</td><td>0</td></tr><tr><td>R539T</td><td>0</td><td>63</td><td>0</td></tr><tr><td>I543T</td><td>0</td><td>34</td><td>0</td></tr><tr><td>G449A</td><td>0</td><td>7</td><td>0</td></tr><tr><td>F446I</td><td>0</td><td>7</td><td>0</td></tr><tr><td>A481V</td><td>0</td><td>4</td><td>0</td></tr><tr><td>F673I</td><td>0</td><td>3</td><td>0</td></tr><tr><td>P667A</td><td>0</td><td>2</td><td>0</td></tr><tr><td>F614L</td><td>0</td><td>1</td><td>0</td></tr><tr><td>S485N</td><td>0</td><td>1</td><td>0</td></tr><tr><td>P443S</td><td>0</td><td>1</td><td>0</td></tr><tr><td>C580Y</td><td>2</td><td>423</td><td>-2</td></tr><tr><td>D353Y</td><td>0</td><td>4</td><td>-2</td></tr><tr><td>K438N</td><td>0</td><td>1</td><td>-2</td></tr><tr><td>P553L</td><td>2</td><td>24</td><td>-3</td></tr><tr><td>P441L</td><td>0</td><td>27</td><td>-3</td></tr><tr><td>G538V</td><td>0</td><td>19</td><td>-3</td></tr><tr><td>P574L</td><td>0</td><td>12</td><td>-3</td></tr><tr><td>V568G</td><td>0</td><td>6</td><td>-3</td></tr></tbody></table>

Table 6 continued on next page
Table 6 continued

<table><thead><tr><th>Mutation</th><th>AFR</th><th>SEA</th><th>CFF53P62</th></tr></thead><tbody><tr><td>P667L</td><td>0</td><td>2</td><td>-3</td></tr><tr><td>S459L</td><td>2</td><td>2</td><td>-4</td></tr><tr><td>N537I</td><td>1</td><td>1</td><td>-4</td></tr><tr><td>Q613L</td><td>1</td><td>0</td><td>-4</td></tr><tr><td>D584V</td><td>0</td><td>3</td><td>-5</td></tr></tbody></table>

DOI: 10.7554/eLife.08714.020

Another important question is whether KPBD mutations are under positive selection in Africa, which is difficult to determine by standard haplotype-based methods because so many independent mutations are involved. This question is further complicated by marked geographical variation in normal levels of genetic diversity, i.e. there are many more rare variants in Africa than Southeast Asia, most likely due to the larger population size and other demographic factors (Manske et al., 2012). Here we address this question by comparing *kelch13* against other genes in the same samples, a process that we refer to as genomic calibration. We show that for most genes the ratio of non-synonymous to synonymous mutations is relatively constant across geographical regions, despite geographic differences in genetic diversity, and that this ratio is correlated with the level of sequence conservation across different *Plasmodium* species. When calibrated against other genes with the same level of cross-species conservation, allowing for geographical differences in the overall level of genetic diversity, *kelch13* shows a marked excess of non-synonymous substitutions in Southeast Asia, but appears normal in Africa. Moreover, KPBD mutations causing radical amino acid changes at highly conserved positions are found at relatively high frequency in Southeast Asia but remain at very low frequency in Africa. Taken together, these findings indicate that non-synonymous KBPD mutations are undergoing strong evolutionary selection in Southeast Asia, whereas those seen in Africa have originated locally and most likely reflect normal variation.

These findings have practical implications for the prevention of artemisinin resistance in Africa, where there is evidently a deep reservoir of low frequency genetic variations that could potentially allow resistance to emerge rapidly as the levels of selective pressure increase. In most parts of Africa, the selective pressure of artemisinin is probably relatively low at present, for several reasons. Artemisinin has been widely used in Southeast Asia for over two decades, whereas its usage in Africa is more recent, and it is estimated that only 20% of infected African children currently have access to frontline treatment ACT medication (World Health Organization, 2014b). Another factor is that people living in regions of high malaria endemicity acquire partial immunity, resulting in asymptomatic infection, so that there is a large reservoir of parasites in Africa that are not exposed to antimalarial drugs because asymptomatic individuals do not seek treatment (Hastings, 2003). The situation could change dramatically as malaria control efforts are intensified, and it will be vital to monitor the effects of major interventions on the emergence of resistance, particularly in African countries that have already achieved relatively low levels of malaria transmission. A key question for the future is whether parasite populations in certain locations possess genetic features that predispose to the emergence of artemisinin resistance, as suggested by the strong association of certain *fd*, *arps10*, *mdr2* and *crt* alleles with resistance-causing KPBD mutations in Southeast Asia (Miotto et al., 2015).

<table><caption>**Table 7.** Frequency of the genetic background alleles across the world. Frequency of the four genetic background alleles identified in *Miotto et al.* (2015) for each geographical region. For each SNP, we show mutation name; chromosome number; nucleotide position; and frequencies of the mutant allele in the various populations.</caption><thead><tr><th>Mutation</th><th>Chr</th><th>Pos</th><th>AFR</th><th>SAS</th><th>SEA</th><th>PNG</th><th>SAM</th></tr></thead><tbody><tr><td>arps10-V127M</td><td>14</td><td>2481070</td><td>0.0%</td><td>0.0%</td><td>59.4%</td><td>0.0%</td><td>0.0%</td></tr><tr><td>fd-D193Y</td><td>13</td><td>748395</td><td>0.1%</td><td>2.2%</td><td>62.8%</td><td>23.9%</td><td>0.0%</td></tr><tr><td>mdr2-T484I</td><td>14</td><td>1956225</td><td>0.1%</td><td>5.7%</td><td>64.2%</td><td>0.4%</td><td>0.0%</td></tr><tr><td>crt-N326S</td><td>7</td><td>405362</td><td>0.8%</td><td>28.2%</td><td>68.6%</td><td>0.1%</td><td>0.0%</td></tr></tbody></table>

DOI: 10.7554/eLife.08714.021
Another concern is that growing resistance to ACT partner drugs, now emerging in Southeast Asia (Saunders et al., 2014) may spread to Africa, or evolve independently there, and lead to increased selective pressure for artemisinin resistance there.

The global spread of resistance to chloroquine and sulfadoxine-pyrimethamine was dominated by hard sweeps of specific haplotypes originating in Southeast Asia, although it is clear that there were also localised emergences (Mita et al., 2009). Although we still know relatively little about the functional properties of different KPBD mutations, it is clear that some are more successful than others, e.g. the C580Y allele has emerged at multiple locations in Southeast Asia and Africa, and a specific C580Y haplotype is approaching fixation in large parts of Western Cambodia (Miotto et al., 2015). The high level of sequence conservation of KPBD across Plasmodium species indicates that mutations in these domains incur fitness costs, and this is supported by the observation that although there have been multiple independent origins of resistance-causing *kelch13* mutations, individual *kelch13* mutations appear to have limited spread, implying that there is a substantial fitness cost in the absence of sustained drug pressure. The danger is that these fitness costs may be compensated by other genetic variants, either in *kelch13* or elsewhere in the genome, and that, as a result of this continuing evolutionary process, parasites in Southeast Asia will progressively acquire higher levels of artemisinin resistance (World Health Organization, 2014a) coupled with strong biological fitness and the ability to propagate across a wide range of vector species. Under these circumstances, the current soft sweep of artemisinin resistance could give way to a pervasive hard sweep with potentially disastrous consequences.

These findings demonstrate the utility of applying genomic epidemiology to identify features of parasite demography and evolution that affect how drug resistance spreads. Future strategies to combat resistance will require better understanding of the evolutionary consequences of malaria control interventions, e.g. how the selective advantage of a resistance allele is counterbalanced by its fitness cost under different control regimes and in different geographical settings. It is now possible to approach this problem prospectively, by conducting systematic spatiotemporal sampling and genome sequencing of the parasite population as an integral part of public health interventions to prevent resistance spreading.

# Materials and methods

## Ethical approval

All samples in this study were derived from blood samples obtained from patients with P. falciparum malaria, collected with informed consent from the patient or a parent or guardian. At each location, sample collection was approved by the appropriate local and institutional ethics committees. The following local and institutional committees that gave ethical approval for the partner studies: Comité d'Éthique, Ministère de la Santé, Bobo-Dioulasso, Burkina Faso; Navrongo Health Research Centre Institutional Review Board, Navrongo, Ghana; Kintampo Health Research Centre Institutional Ethics Committee, Kintampo, Ghana; Noguchi Memorial Institute for Medical Research Institutional Review Committee, University of Ghana, Legon, Ghana; Ghana Health Service Ethical Review Committee, Accra, Ghana; Gambia Government/MRC Joint Ethics Committee, Banjul, The Gambia; Comité d'Ethique National Pour la Recherche en Santé, Guinea; Ethics Committee of Faculté de Médecine, de Pharmacie et d'Odonto-Stomatologie, University of Bamako, Bamako, Mali; Ethical Review Committee, University of Ilorin Teaching Hospital, Ilorin, Nigeria; Institutional Review Board, Faculty of Health Sciences, University of Buea, Cameroon; Comité d'Ethique, Ecole de Santé Publique, Université de Kinshasa, Ministère de l'Enseignement Superieur, Universitaire et Recherche Scientifique, D. R. Congo; Comité National d'Ethique auprès du Ministère de la Santé Publique, Madagascar; Institutional Review Committee, Med Biotech Laboratories, Kampala, Uganda and Uganda National Council for Sciences and Technology (UNCST); College of Medicine Research Ethics Committee, University of Malawi, Blantyre, Malawi; KEMRI National Ethical Review Committee, Kenya; Medical Research Coordinating Committee of the National Institute for Medical Research, Tanzania; Ethical Review Committee, Bangladesh Medical Research Council, Bangladesh; Ethics Committee of the International Centre for Diarrheal Disease Research, Bangladesh; Institutional Ethical Review Committee, Department of Medical Research (Lower Myanmar); Ministry of Health, Government of The Republic of the Union of Myanmar; National Ethics Committee for
Health Research, Ministry of Health, Phnom Penh, Cambodia; Ministry of Health National Ethics Committee For Health Research, Laos; Ethics Committee, Faculty of Tropical Medicine, Mahidol University, Bangkok, Thailand; Ethical Committee, Hospital for Tropical Diseases, Ho Chi Minh City, Vietnam; Eijkman Institute Research Ethics Commission, Jakarta, Indonesia; Institutional Review Board, Papua New Guinea Institute of Medical Research, Goroka, Papua New Guinea; Institutional Review Board, International Center for Medical Research and Training, Cali, Colombia; Institutional Review Board, Universidad Nacional de la Amazonia Peruana, Iquitos, Peru; Human Research Ethics Committee of NT Department of Health and Families and Menzies School of Health Research, Darwin, Australia; Institutional Review Board, New York University Medical School, NY, USA; Institutional Review Board, National Institute of Allergy and Infectious Diseases, Bethesda, MD, USA; Institutional Review Board, Walter Reed Army Institute of Research, Washington DC, USA; Ethics Review Committee, World Health Organization, Geneva, Switzerland; Ethics Committee of the Faculty of Medicine, Heidelberg, Germany; Ethics Committee of the Medical University of Vienna; Ethics Committee, London School of Hygiene and Tropical Medicine, London, UK; Oxford Tropical Research Ethics Committee, Oxford, UK.

## Sample preparation, sequencing and genotyping

DNA was extracted directly from blood samples taken from patients at admission time, after leukocyte depletion to minimize human DNA. Leukocyte depletion was achieved by CF11 filtration in most samples (Venkatesan et al., 2012), or alternatively by Lymphoprep density gradient centrifugation (Axis-Shield, Dundee, UK) followed by Plasmodipur filtration (Euro-Diagnostica, Malmö, Sweden) (Auburn et al., 2011) or by Plasmodipur filtration alone. Genomic DNA was extracted using the QIAamp DNA Blood Midi or Maxi Kit (Qiagen, Hilden, Germany), and quantities of human and Plasmodium DNA were determined by fluorescence analysis using a Qubit instrument (Invitrogen, Carlsbad, California) and multi-species quantitative PCR (Q-PCR) using the Roche Lightcycler 480 II system, as described previously (Manske et al., 2012). Samples with >50 ng DNA and <80% human DNA contamination were selected for sequencing on the Illumina HiSeq platform following the manufacturer's standard protocols (Bentley et al., 2008). Paired-end sequencing reads of length 200-300 bp were obtained, generating approximately 1 Gbp of read data per sample. All short read sequence data have been deposited in the European Nucleotide Archive (http://www.ebi.ac.uk/ena/data/search/?query=plasmodium), and metadata will be released at the time of publication.

Polymorphism discovery, quality control and sample genotyping followed a process described elsewhere (Manske et al., 2012). Short sequence reads from 3,411 *P. falciparum* samples included in the MalariaGEN Plasmodium falciparum Community Project (http://www.malariagen.net/projects/parasite/pf) were aligned against the *P. falciparum* 3D7 reference sequence V3 (ftp://ftp.sanger.ac.uk/pub/pathogens/Plasmodium/falciparum/3D7/3D7.latest_version/version3/), using the bwa program (Li and Durbin, 2009) (http://bio-bwa.sourceforge.net/) as previously described (Manske et al., 2012), to identify an initial global set of 4,305,957 potential SNPs. This list was then used to guide stringent re-alignment using the SNP-o-matic algorithm (Manske and Kwiatkowski, 2009), to reduce misalignment errors. The stringent alignments were then examined by a series of quality filters, with the aim of removing alignment artefacts and their sources. In particular, the following were removed: a) non-coding SNPs; b) SNPs where polymorphisms have extremely low support (<10 reads in one sample); c) SNPs with more than two alleles, with the exception of loci known to be important for drug resistance, which were manually verified for artifacts; d) SNPs where coverage across samples is lower than the 25th percentile and higher than the 95th percentile of coverage in coding SNPs (these thresholds were determined from artifact analysis); e) SNPs located in regions of relatively low uniqueness (Manske et al., 2012); f) SNPs where heterozygosity levels were found to be inconsistent with the heterozygosity distribution at the SNP's allele frequency; and g) SNPs where genotype could not be established in at least 70% of the samples. These analyses produced a final list of 935,601 high-quality SNPs in the 14 chromosomes of the nuclear genome, whose genotypes were used for analysis in this study.

All samples were genotyped at each high-quality SNP by a single allele, based on the number of reads observed for the two alleles at that position in the sample. At positions with fewer than 5 reads, the genotype was undetermined (no call was made). At all other positions, the sample was determined to be *heterozygous* if both alleles were each observed in more than 2 reads; otherwise,
the sample was called as homozygous for the allele observed in the majority of reads. For the purposes of estimating allele frequencies and genetic distances, a within-sample allele frequency ($f_w$) was also assigned to each valid call. For heterozygous calls, $f_w$ was estimated as ratio of non-reference read count to reference read count; homozygous calls were assigned $f_w = 0$ when called with the reference allele, and $f_w = 1$ when called with the non-reference allele.

The genotype of *kelch13* was derived from read counts at non-synonymous SNP in the KPBD using a procedure described previously (Miotto et al., 2015)

## Frequency estimation and clustering

For a given population *P*, we estimated the non-reference allele frequency (NRAF) at a given SNP as the mean of the within-sample allele frequency ($f_w$) for all samples in *P* which have a valid genotype at that SNP. The minor allele frequency (MAF) at is the computed as $min(NRAF, (1 - NRAF))$.

To investigate the global population structure, we started by computing an NxN pairwise distance matrix, where *N* is the number of samples. Each cell of the matrix contained an estimate of genetic distance between the relevant pair of samples, obtained by summing the pairwise distance, estimated from within-sample allele frequency ($f_w$), at each SNP in the 100 kbp window considered. When comparing a pair of samples $s_A$ and $s_B$ at a single SNP *i* where a genotype could be called in each sample, with within-sample allele frequencies $f_A$ and $f_B$ respectively, the distance $d_{AB}$ was estimated as $d_{AB} = f_A (1 - f_B) + f_B (1 - f_A)$. The genome-wide distance $D_{AB}$ between the two samples is then calculated as

$$D_{AB} = \frac{\alpha}{n_{AB}} \sum_{i} w_{i} d_{AB}$$

where $n_{AB}$ is the number of SNPs where both samples could be genotyped, $w_i$ is an LD weighting factor (see below) and $\alpha$ is a scaling constant, equal to 70% of the number of coding positions in the genome (since our genotyping covers approximately 70% of the coding genome). The exact value of $\alpha$ is uninfluential towards the analyses conducted in this study. The LD weighting factor, which corrects for the cumulative contribution of physically linked polymorphisms, was computed at each SNP *i* with $MAF \ge 0.1$ in our sample set, by considering a window of *m* SNPs ($j = 0.. m$) centred at *i*. For each *j*, we computed the squared correlation coefficient $r_{ij}^2$ between SNPs *i* and *j*. Ignoring positions *j* where $r_{ij}^2 < 0.1$, the weighting $w_i$ was computed by

$$w_i = \frac{1}{1 + \sum_j r_{ij}^2}$$

A neighbour-joining tree was then produced using the nj implementation in the R ape package. Principal coordinate analysis (PCoA) was performed using the same pairwise distance matrices using the Classical Multidimensional Scaling (MDS) method (Gower, 1966) PCoA is a computationally efficient variant of principal component analysis (PCA) in which a pairwise distance matrix is used as input, rather than a table of genotypes. The matrix was supplied as input to the MDS algorithm, using the R language cmdscale implementation.

## Conservation scoring and ancestral allele analysis

We analyzed homologous protein sequences of *kelch13* genes for seven *Plasmodium* species for which high-quality sequence data were available: *P. falciparum*, *P. reichenowi*, *P. vivax*, *P. knowlesi*, *P. yoelii*, *P. berghei*, and *P. chabaudi*. The sequences were retrieved from the OrthoMCL cluster ORTHOMCL894 in GeneDB (http://www.genedb.org/), and a multiple alignment was obtained using ClustalW (Larkin et al., 2007) at default settings. In turn, we considered each pair alignment of *P. falciparum* with one of the remaining species, assigning a substitution score to each *kelch13* amino acid position, derived from the CCF53P62 substitution matrix. Although this matrix was chosen due to its suitability for AT-rich codon biases (Brick and Pizzi, 2008), we found that use of the more commonly used BLOSUM62 matrix (Henikoff and Henikoff, 1992) did not have a significant effect on the results. Finally, each amino acid position was assigned a conservation score for the pair alignment, equal to the mean of the substitution scores in a 9-residue window centered at that position.

To reconstruct putative *ancestral* and *derived* alleles in the propeller domain, we catalogued all polymorphic positions in the multiple sequence alignment. We organized the seven species into
three groups by similarity: Laverania (P. falciparum, P. reichenowi), primate Plasmodia (P. vivax, P. knowlesi) and rodent Plasmodia (P. yoelii, P. berghei, and P. chabaudi), and observed that at each position, only one of the groups presented an allele different from that in the remaining groups. This group-specific allele was labelled as a putative derived allele, and the alternative allele as ancestral. (see Table 5). Rodent species were found to carry the highest number of derived alleles, and therefore deemed to be good comparators for genome-wide conservation scoring. P. chabaudi was selected as a representative species in this group, and used for subsequent comparative analyses.

We estimated a gene conservation score for every P. falciparum gene for which a P. chabaudi orthologue sequence could be obtained from PlasmoDB (http://www.plasmodb.org/). The details of the method are described elsewhere (Gardner et al., 2011). Briefly, alignments of orthologous protein sequences were performed using ClustalW (Larkin et al., 2007) at default settings, and each amino acid position was assigned a CCF53P62 substitution score (see above). The gene conservation score assigned was equal to the mean substitution score for all amino acid positions in the gene.

## Hydrophobicity analysis

Each amino acid position in the *kelch13* was assigned a hydrophobicity score, estimated by computing the mean of the Kyte-Doolittle index in a 14-residue window centered at the position, using the protein sequence translated from the 3D7 reference sequence.

## Chromosome painting

To reconstruct the probable origin of *kelch13* flanking haplotypes in African mutants, we used chromosome painting (Lawson et al., 2012), a method that compares haplotypes in a sample to those in the remaining samples, and estimates probabilities that a genome fragment originates each population, by identifying individuals that share the same haplotype. For all African and SEA samples, we applied chromosome painting across the 250 kbp flanking regions on each side of the *kelch13* gene. For each sample, haplotypes surrounding genome loci (chunks) were assigned posterior copying probabilities with respect to all remaining samples (unrestricted painting). We aggregated these probabilities according to the geographical origin of the donor samples, assigning to each chunk a probability that it originates from each of the populations. For each sample, we then estimated the expected fraction of chunks copied from each population.

In this analysis, we assumed a mutation rate per base per generation of 3.9⁻¹⁰ (Claessens et al., 2014) and a uniform recombination map. Since the two populations differ substantially in effective recombination rates (Mu et al., 2005), we assumed a conservative recombination rate of 30 kbp/cM. Effective population size was initially set to 10,000 and optimized by 10 iterations of the expectation-maximization procedure (Lawson et al., 2012). To account for the presence of heterozygous genotypes due to mixed infections, we modified the matrix of emission probabilities by introducing a novel parameter (ε, set to 10⁻⁸) to represent the probability of emitting a mixed call. We repeated the analysis varying this parameter set, to assess the effects of misspecification, and results were found to be very similar qualitatively (data not shown).

## Data access and URLs

Illumina sequence reads have been submitted to the European Nucleotide Archive with study accessions ERP000190 (www.ebi.ac.uk/ena/data/view/ERP000190) and ERP000197 (www.ebi.ac.uk/ena/data/view/ERP000197). ENA accession numbers and metadata for the samples used in this paper can be obtained via the MalariaGEN website, which also provides access to genotype calls on individual samples (www.malariagen.net/resource/16). Further details of all SNPs reported in this dataset including their genome coverage, mapping quality and allele frequencies in different populations, together with tools for querying the data, can be explored at www.malariagen.net/apps/pf.

## Acknowledgements

This study was conducted by the MalariaGEN Plasmodium falciparum Community Project, and was made possible by clinical parasite samples contributed by partner studies, whose investigators are represented in the author list. RA and OM contributed equally. In addition, the authors would like to thank the following individuals, who contributed to partner studies or to the MalariaGEN Resource Centre, making this study possible: James Abugri, Nicholas Amoako, Steven M Kiara, John Okombo,
Rogelin Raherinjafy, Seheno Razanatsiorimalala, Hongying Jiang, Xin-zhuan Su. The sequencing, analysis, informatics and management of the Community Project are supported by the Wellcome Trust through Sanger Institute core funding (098051) and a Strategic Award (090770/Z/09/Z) and by the MRC Centre for Genomics and Global Health which is jointly funded by the Medical Research Council and the Department for International Development (DFID) (G0600718; M006212). AEB and IM acknowledge the Victorian State Government Operational Infrastructure Support and Australian Government NHMRC IRIISS.

## Additional information

### Group author details

MalariaGEN Plasmodium falciparum Community Project

Roberto Amato: Wellcome Trust Sanger Institute, Cambridge, United Kingdom; Centre for Genomics and Global Health, Wellcome Trust Centre for Human Genetics, Oxford, United Kingdom; Olivo Miotto: Mahidol-Oxford Tropical Medicine Research Unit, Mahidol University, Bangkok, Thailand; Centre for Genomics and Global Health, Wellcome Trust Centre for Human Genetics, Oxford, United Kingdom; Wellcome Trust Sanger Institute, Cambridge, United Kingdom; Charles J Woodrow: Mahidol-Oxford Tropical Medicine Research Unit, Mahidol University, Bangkok, Thailand; Centre for Tropical Medicine and Global Health, Nuffield Department of Medicine, University of Oxford, Oxford, United Kingdom; Jacob Almagro-Garcia: Centre for Genomics and Global Health, Wellcome Trust Centre for Human Genetics, Oxford, United Kingdom; Wellcome Trust Sanger Institute, Cambridge, United Kingdom; Ipsita Sinha: Mahidol-Oxford Tropical Medicine Research Unit, Mahidol University, Bangkok, Thailand; Susana Campino: Wellcome Trust Sanger Institute, Cambridge, United Kingdom; Daniel Mead: Wellcome Trust Sanger Institute, Cambridge, United Kingdom; Eleanor Drury: Wellcome Trust Sanger Institute, Cambridge, United Kingdom; Mihir Kekre: Wellcome Trust Sanger Institute, Cambridge, United Kingdom; Mandy Sanders: Wellcome Trust Sanger Institute, Cambridge, United Kingdom; Alfred Amambua-Ngwa: Medical Research Council Unit, Banjul, The Gambia; Chanaki Amaratunga: National Institute of Allergy and Infectious Diseases, National Institutes of Health, Rockville, United States; Lucas Amenga-Etego: Navrongo Health Research Centre, Navrongo, Ghana; Voahangy Andrianaranjaka: Malaria Research Unit, Institut Pasteur de Madagascar, Antananarivo, Madagascar; Tobias Apinjoh: University of Buea, Buea, Cameroon; Elizabeth Ashley: Mahidol-Oxford Tropical Medicine Research Unit, Mahidol University, Bangkok, Thailand; Sarah Auburn: Global and Tropical Health Division, Menzies School of Health Research and Charles Darwin University, Darwin, Australia; Gordon A Awandare: West African Centre for Cell Biology of Infectious Pathogens, College of Basic and Applied Sciences, University of Ghana, Legon, Ghana; Vito Baraka: National Institute for Medical Research, Tanga - Centre, Tanga, Tanzania; International Health Unit, Department of Epidemiology, University of Antwerp, Antwerp, Belgium; Alyssa Barry: Division of Population Health and Immunity, Walter and Eliza Hall Institute of Medical Research, Parkville, Australia; Department of Medical Biology, University of Melbourne, Carlton, Australia; Maciej F Boni: University of Oxford, Ho Chi Minh City, Vietnam; Centre for Tropical Medicine and Global Health, Nuffield Department of Medicine, University of Oxford, Oxford, United Kingdom; Steffen Borrmann: Kenya Medical Research Institute/Wellcome Trust Collaborative Research Program, Kilifi, Kenya; Institute for Tropical Medicine, University of Tübingen, Tübingen, Germany; Teun Bousema: London School of Hygiene and Tropical Medicine, London, United Kingdom; Oralee Branch: New York University School of Medicine, New York, United States; Peter C Bull: Kenya Medical Research Institute/Wellcome Trust Collaborative Research Program, Kilifi, Kenya; Centre for Tropical Medicine and Global Health, Nuffield Department of Medicine, University of Oxford, Oxford, United Kingdom; Kesinee Chotivanich: Faculty of Tropical Medicine, Mahidol University, Bangkok, Thailand; David J Conway: Medical Research Council Unit, Banjul, The Gambia; London School of Hygiene and Tropical Medicine, London, United Kingdom; Alister Craig: Liverpool School of Tropical Medicine, Liverpool, United Kingdom; Nicholas P Day: Mahidol-Oxford Tropical Medicine Research Unit, Mahidol University, Bangkok, Thailand; Centre for Tropical Medicine and Global Health, Nuffield Department of Medicine, University of Oxford, Oxford, United Kingdom; Abdoulaye Djimdé: Malaria Research and
Training Centre, Faculty of Medicine, University of Bamako, Bamako, Mali; Christiane Dolecek: Oxford University Clinical Research Unit, Ho Chi Minh City, Vietnam; Centre for Tropical Medicine and Global Health, Nuffield Department of Medicine, University of Oxford, Oxford, United Kingdom; Arjen M Dondorp: Mahidol-Oxford Tropical Medicine Research Unit, Mahidol University, Bangkok, Thailand; Centre for Tropical Medicine and Global Health, Nuffield Department of Medicine, University of Oxford, Oxford, United Kingdom; Chris Drakeley: London School of Hygiene and Tropical Medicine, London, United Kingdom; Patrick Duffy: National Institute of Allergy and Infectious Diseases, National Institutes of Health, Rockville, United States; Diego F Echeverry: Department of Entomology, Purdue University, West Lafayette, United States; International Center for Medical Research and Training, CIDEIM, Cali, Colombia; Thomas G Egwang: Med Biotech Laboratories, Kampala, Uganda; Rick M Fairhurst: National Institute of Allergy and Infectious Diseases, National Institutes of Health, Rockville, United States; Md. Abul Faiz: Malaria Research Group and Dev Care Foundation, Dhaka, Bangladesh; Caterina I Fanello: Mahidol-Oxford Tropical Medicine Research Unit, Mahidol University, Bangkok, Thailand; Tran Tinh Hien: Oxford University Clinical Research Unit, Ho Chi Minh City, Vietnam; Centre for Tropical Medicine and Global Health, Nuffield Department of Medicine, University of Oxford, Oxford, United Kingdom; Abraham Hodgson: Navrongo Health Research Centre, Navrongo, Ghana; Mallika Imwong: Mahidol-Oxford Tropical Medicine Research Unit, Mahidol University, Bangkok, Thailand; Deus Ishengoma: National Institute for Medical Research, Tanga - Centre, Tanga, Tanzania; Pharath Lim: National Institute of Allergy and Infectious Diseases, National Institutes of Health, Rockville, United States; National Centre for Parasitology, Entomology and Malaria Control, Phnom Penh, Cambodia; Chanthap Lon: Department of Immunology and Medicine, US Army Medical Component, Armed Forces Research Institute of Medical Sciences, Bangkok, Thailand; Jutta Marfurt: Global and Tropical Health Division, Menzies School of Health Research and Charles Darwin University, Darwin, Australia; Kevin Marsh: Kenya Medical Research Institute/Wellcome Trust Collaborative Research Program, Kilifi, Kenya; Mayfong Mayxay: Lao-Oxford-Mahosot Hospital-Wellcome Trust Research Unit, Vientiane, Lao PDR; Faculty of Postgraduate Studies, University of Health Sciences, Vientiane, Lao PDR; Pascal Michon: Faculty of Medicine and Health Sciences, Divine Word University, Madang, Papua New Guinea; Victor Mobegi: Medical Research Council Unit, Banjul, The Gambia; London School of Hygiene and Tropical Medicine, London, United Kingdom; Olugbenga A Mokuolu: Department of Paediatrics and Child Health, University of Ilorin, Ilorin, Nigeria; Jacqui Montgomery: Department of Biology and Entomology, Pennsylvania State University, University Park, United States; Ivo Mueller: Department of Medical Biology, University of Melbourne, Carlton, Australia; Myat Phone Kyaw: Department of Medical Research Lower Myanmar, Ministry of Health, Yangon, Myanmar; Paul N Newton: Lao-Oxford-Mahosot Hospital-Wellcome Trust Research Unit, Vientiane, Lao PDR; Francois Nosten: Mahidol-Oxford Tropical Medicine Research Unit, Mahidol University, Mae Sot, Thailand; Shoklo Malaria Research Unit, Mae Sot, Thailand; Rintis Noviyanti: Eijkman Institute for Molecular Biology, Jakarta, Indonesia; Alexis Nzila: Department of Biology, King Fahd University of Petroleum and Minerals, Dhahran, Saudi Arabia; Harold Ocholla: Malawi Liverpool Wellcome Trust Clinical Research Programme, Malawi; Abraham Oduro: Navrongo Health Research Centre, Navrongo, Ghana; Marie Onyamboko: University of Kinshasa, Kinshasa, Democratic Republic of Congo; Jean-Bosco Ouedraogo: Institut de Recherche en Sciences de la Sante, Direction Régionale de l'Ouést, Bobo-Dioulasso, Burkina Faso; Aung Pyae P Phyo: Shoklo Malaria Research Unit, Bangkok, Thailand; Shoklo Malaria Research Unit, Mae Sot, Thailand; Christopher Plowe: Institute for Global Health, University of Maryland School of Medicine, Baltimore, United States; Ric N Price: Centre for Tropical Medicine and Global Health, Nuffield Department of Medicine, University of Oxford, Oxford, United Kingdom; Global and Tropical Health Division, Menzies School of Health Research and Charles Darwin University, Darwin, Australia; Sasithon Pukrittayakamee: Faculty of Tropical Medicine, Mahidol University, Bangkok, Thailand; Milijaona Randrianarivelojosia: Malaria Research Unit, Institut Pasteur de Madagascar, Antananarivo, Madagascar; Pascal Ringwald: Global Malaria Programme, World Health Organization, Geneva, Switzerland; Lastenia Ruiz: Universidad Nacional de la Amazonia Peruana, Iquitos, Peru; David Saunders: Department of Immunology and Medicine, US Army Medical Component, Armed Forces Research Institute of Medical Sciences, Bangkok, Thailand; Alex Shayo: Department of Biotechnology and Bioinformatics, The University of Dodoma, Dodoma, Tanzania; Peter Siba: Papua New Guinea Institute of Medical Research, Madang, Papua New Guinea; Shannon Takala-Harrison: Institute for Global Health, University of Maryland School of
Medicine, Baltimore, United States; Thuy-Nhien N Thanh: Centre for Tropical Medicine and Global Health, Nuffield Department of Medicine, University of Oxford, Oxford, United Kingdom; Oxford University Clinical Research Unit, Ho Chi Minh City, Vietnam; Vandana Thathy: Kenya Medical Research Institute/Wellcome Trust Collaborative Research Program, Kilifi, Kenya; Federica Verra: London School of Hygiene and Tropical Medicine, London, United Kingdom; Jason Wendler: Centre for Genomics and Global Health, Wellcome Trust Centre for Human Genetics, Oxford, United Kingdom; Nicholas J White: Mahidol-Oxford Tropical Medicine Research Unit, Mahidol University, Bangkok, Thailand; Centre for Tropical Medicine and Global Health, Nuffield Department of Medicine, University of Oxford, Oxford, United Kingdom; Htut Ye: Department of Medical Research Lower Myanmar, Ministry of Health, Yangon, Myanmar; Victoria J Cornelius: Centre for Genomics and Global Health, Wellcome Trust Centre for Human Genetics, Oxford, United Kingdom; Rachel Giacomantonio: Centre for Genomics and Global Health, Wellcome Trust Centre for Human Genetics, Oxford, United Kingdom; Wellcome Trust Sanger Institute, Cambridge, United Kingdom; Dawn Muddyman: Wellcome Trust Sanger Institute, Cambridge, United Kingdom; Christa Henrichs: Centre for Genomics and Global Health, Wellcome Trust Centre for Human Genetics, Oxford, United Kingdom; Cinzia Malangone: Wellcome Trust Sanger Institute, Cambridge, United Kingdom; Dushyanth Jyothi: Wellcome Trust Sanger Institute, Cambridge, United Kingdom; Richard D Pearson: Centre for Genomics and Global Health, Wellcome Trust Centre for Human Genetics, Oxford, United Kingdom; Wellcome Trust Sanger Institute, Cambridge, United Kingdom; Julian C Rayner: Wellcome Trust Sanger Institute, Cambridge, United Kingdom; Gilean McVean: Centre for Genomics and Global Health, Wellcome Trust Centre for Human Genetics, Oxford, United Kingdom; Kirk A Rockett: Centre for Genomics and Global Health, Wellcome Trust Centre for Human Genetics, Oxford, United Kingdom; Wellcome Trust Sanger Institute, Cambridge, United Kingdom; Alistair Miles: Centre for Genomics and Global Health, Wellcome Trust Centre for Human Genetics, Oxford, United Kingdom; Wellcome Trust Sanger Institute, Cambridge, United Kingdom; Paul Vauterin: Centre for Genomics and Global Health, Wellcome Trust Centre for Human Genetics, Oxford, United Kingdom; Ben Jeffery: Centre for Genomics and Global Health, Wellcome Trust Centre for Human Genetics, Oxford, United Kingdom; Magnus Manske: Wellcome Trust Sanger Institute, Cambridge, United Kingdom; Jim Stalker: Wellcome Trust Sanger Institute, Cambridge, United Kingdom; Bronwyn Maclnnis: Wellcome Trust Sanger Institute, Cambridge, United Kingdom; Dominic P Kwiatkowski: Wellcome Trust Sanger Institute, Cambridge, United Kingdom; Centre for Genomics and Global Health, Wellcome Trust Centre for Human Genetics, Oxford, United Kingdom

## Competing interests

GM: Reviewing editor, eLife. The other authors declare that no competing interests exist.

## Funding

<table><thead><tr><td>Funder</td><td>Grant reference number</td><td>Author</td></tr></thead><tbody><tr><td>Wellcome Trust</td><td>WT Sanger Institute Core Funding, 98051</td><td>Roberto Amato<br>Susana Campino<br>Daniel Mead<br>Eleanor Drury<br>Mihir Kekre<br>Mandy Sanders<br>Victoria J Cornelius<br>Rachel Giacomantonio<br>Dawn Muddyman<br>Cinzia Malangone<br>Dushyanth Jyothi<br>Richard D Pearson<br>Julian C Rayner<br>Kirk A Rockett<br>Alistair Miles<br>Heinrich M Manske<br>Jim Stalker<br>Bronwyn Maclnnis<br>Dominic P Kwiatkowski</td></tr></tbody></table>
<table><tr><td>Wellcome Trust</td><td>WT Centre for Human Genetics Core Funding, 090532/Z/09/Z</td><td>Roberto Amato<br>Jacob Almagro-Garcia<br>Victoria J Cornelius<br>Rachel Giacomantonio<br>Christa Henrichs<br>Richard D Pearson<br>Gilean McVean<br>Kirk A Rockett<br>Alistair Miles<br>Paul Vauterin<br>Ben Jeffery<br>Dominic P Kwiatkowski</td></tr><tr><td>Wellcome Trust</td><td>WT Strategic Award, 90770</td><td>Roberto Amato<br>Olivo Miotto<br>Jacob Almagro-Garcia<br>Susana Campino<br>Daniel Mead<br>Eleanor Drury<br>Mihir Kekre<br>Lucas Amenga-Etego<br>Tobias Apinjoh<br>Abdoulaye Djimdé<br>Deus Ishengoma<br>Milijaona Randrianarivelojosia<br>Victoria J Cornelius<br>Rachel Giacomantonio<br>Dawn Muddyman<br>Christa Henrichs<br>Cinzia Malangone<br>Dushyanth Jyothi<br>Richard D Pearson<br>Kirk A Rockett<br>Alistair Miles<br>Paul Vauterin<br>Ben Jeffery<br>Heinrich M Manske<br>Jim Stalker<br>Bronwyn MacInnis<br>Dominic P Kwiatkowski</td></tr><tr><td>Medical Research Council</td><td>Centre for Genomics and Global Health, G0600718</td><td>Roberto Amato<br>Olivo Miotto<br>Jacob Almagro-Garcia<br>Lucas Amenga-Etego<br>Tobias Apinjoh<br>Abdoulaye Djimdé<br>Deus Ishengoma<br>Victoria J Cornelius<br>Rachel Giacomantonio<br>Dawn Muddyman<br>Gilean McVean<br>Kirk A Rockett<br>Alistair Miles<br>Paul Vauterin<br>Ben Jeffery<br>Heinrich M Manske<br>Jim Stalker<br>Bronwyn MacInnis<br>Dominic P Kwiatkowski</td></tr><tr><td>Bill and Melinda Gates Foundation</td><td>Grant OPP1040463</td><td>Roberto Amato<br>Olivo Miotto<br>Charles J Woodrow<br>Ipsita Sinha<br>Elizabeth Ashley<br>Nicholas P Day<br>Arjen M Dondorp<br>Mallika Imwong<br>Nicholas J White<br>Dominic P Kwiatkowski</td></tr></table>
Wellcome Trust

WT Mahidol University
Oxford Tropical Medicine
Research Programme,
0892375/H/09/Z

Olivo Miotto
Charles J Woodrow
Ipsita Sinha
Elizabeth Ashley
Nicholas P Day
Arjen M Dondorp
Caterina I Fanello
Mayfong Mayxay
Paul N Newton
François Nosten
Nicholas J White

National Institute of Allergy and Infectious Diseases

Intramural Research Program

Chanaki Amaratunga
Rick M Fairhurst
Pharath Lim

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

## Author contributions

RA, OM, DPK, Conception and design, Analysis and interpretation of data, Drafting or revising the article, Approved the submitted manuscript; DMe, ED, SC, MK, MS, Developed and implemented sample processing and sequencing library preparation methods, Drafting or revising the article, Approved the submitted manuscript; AAN, CA, LAE, VA, TA, EA, SA, VB, AB, MFB, TB, PCB, DJC, AC, NPD, AD, CDo, AMD, CDr, PD, DFE, TGE, RMF, MAF, CIF, TTH, AH, MI, DI, PL, CL, JMa, KM, MMay, PM, VM, OAM, JMo, IM, MPK, PNN, FN, RN, AN, HO, AO, MO, J-BO, APPP, CP, RNP, SP, MR, PR, LR, DS, AS, PS, STH, T-NNT, VT, FV, JW, NJW, HY, GAA, SB, OB, KC, Carried out clinical studies and/or laboratory work to obtain parasite samples, Drafting or revising the article, Approved the submitted manuscript; CM, DJ, MMan, JS, Development and management of data production pipelines, Drafting or revising the article, Approved the submitted manuscript; CJW, JAG, IS, AM, RDP, PV, BJ, Developed analysis, software tools and methods, Drafting or revising the article, Approved the submitted manuscript; VJC, RG, DMu, CH, GM, KAR, BM, JCR, Contributed to study design and management, Drafting or revising the article, Approved the submitted manuscript.

## Additional files

### Supplementary files

* Supplementary file 1. Table of N/S ratio and conservation per gene. For each *P. falciparum* gene, this table lists: the systematic ID, position and description of the gene; the ortholog gene in *P. chabaudi*; the conservation score; the count of all (N)on-synonymous and (S)ynonymous mutation in Africa (AFR) and Southeast Asia (SEA); the count of the rare (i.e. present in only 1 or 2 samples) mutations; the N/S log fold-change in SEA vs AFR; the Fisher's test p-value of N/S in AFR vs SEA.

DOI: 10.7554/eLife.08714.022

### Major datasets

The following datasets were generated:

<table><thead><tr><td>Author(s)</td><td>Year</td><td>Dataset title</td><td>Dataset URL</td><td>Database, license, and accessibility information</td></tr></thead><tbody><tr><td>MalariaGEN Plasmodium falciparum Community Project</td><td>2016</td><td>Sample data and catalogue of genetic variations</td><td>http://www.malariagen.net/resource/16</td><td>Access information stated in Web page</td></tr><tr><td>MalariaGEN Plasmodium falciparum Community Project</td><td>2016</td><td>Plasmodium falciparum natural genome variation</td><td>http://www.ebi.ac.uk/ena/data/view/ERP000190</td><td>Publicly available at European Nucleotide Archive (accession no. ERP000190)</td></tr><tr><td>MalariaGEN Plasmodium falciparum Community Project</td><td>2016</td><td>Plasmodium falciparum Illumina sequencing R&amp;D</td><td>http://www.ebi.ac.uk/ena/data/view/ERP000197</td><td>Publicly available at the European Nucleotide Archive (accession no. ERP000197)</td></tr></tbody></table>
## References

Adams J, Kelso R, Cooley L. 2000. The kelch repeat superfamily of proteins: propellers of cell function. *Trends in Cell Biology* **10**:17-24. doi: 10.1016/S0962-8924(99)01673-6

Ariey F, Witkowski B, Amaratunga C, Beghain J, Langlois A-C, Khim N, Kim S, Duru V, Bouchier C, Ma L, Lim P, Leang R, Duong S, Sreng S, Suon S, Chuor CM, Bout DM, Ménard S, Rogers WO, Genton B, Fandeur T, Miotto O, Ringwald P, Le Bras J, Berry A, Barale J-C, Fairhurst RM, Benoit-Vical F, Mercereau-Puijalon O, Ménard D. 2014. A molecular marker of artemisinin-resistant plasmodium falciparum malaria. *Nature* **505**:50-55. doi: 10.1038/nature12876

Ashley EA, Dhorda M, Fairhurst RM, Amaratunga C, Lim P, Suon S, Sreng S, Anderson JM, Mao S, Sam B, Sopha C, Chuor CM, Nguon C, Sovannaroth S, Pukrittayakamee S, Jittamala P, Chotivanich K, Chutasmit K, Suchatsoonthorn C, Runcharoen R, Hien TT, Thuy-Nhien NT, Thanh NV, Phu NH, Htut Y, Han KT, Aye KH, Mokuolu OA, Olaosebikan RR, Folaranmi OO, Mayxay M, Khanthavong M, Hongvanthong B, Newton PN, Onyamboko MA, Fanello CI, Tshefu AK, Mishra N, Valecha N, Phyo AP, Nosten F, Yi P, Tripura R, Borrmann S, Bashraheil M, Peshu J, Faiz MA, Ghose A, Hossain MA, Samad R, Rahman MR, Hasan MM, Islam A, Miotto O, Amato R, MacInnis B, Stalker J, Kwiatkowski DP, Bozdech Z, Jeeyapant A, Cheah PY, Sakulthaew T, Chalk J, Intharabut B, Silamut K, Lee SJ, Vihokhern B, Kunasol C, Imwong M, Tarning J, Taylor WJ, Yeung S, Woodrow CJ, Flegg JA, Das D, Smith J, Venkatesan M, Plowe CV, Stepniewska K, Guerin PJ, Dondorp AM, Day NP, White NJ. Tracking Resistance to Artemisinin Collaboration (TRAC) Tracking Resistance to Artemisinin Collaboration (TRAC) Tracking Resistance to Artemisinin Collaboration (TRAC). 2014. Spread of artemisinin resistance in plasmodium falciparum malaria. *New England Journal of Medicine* **371**:411-423. doi: 10.1056/NEJMoa1314981

Auburn S, Campino S, Clark TG, Djimde AA, Zongo I, Pinches R, Manske M, Mangano V, Alcock D, Anastasi E, Maslen G, MacInnis B, Rockett K, Modiano D, Doumbo OK, Ouédraogo JB, Kwiatkowski DP, Dimopoulos G. 2011. An effective method to purify plasmodium falciparum DNA directly from clinical blood samples for whole genome high-throughput sequencing. *PLoS ONE* **6**:e22213. doi: 10.1371/journal.pone.0022213

Bentley DR, Balasubramanian S, Swerdlow HP, Smith GP, Milton J, Brown CG, Hall KP, Evers DJ, Barnes CL, Bignell HR, Boutell JM, Bryant J, Carter RJ, Keira Cheetham R, Cox AJ, Ellis DJ, Flatbush MR, Gormley NA, Humphray SJ, Irving LJ, Karbelashvili MS, Kirk SM, Li H, Liu X, Maisinger KS, Murray LJ, Obradovic B, Ost T, Parkinson ML, Pratt MR, Rasolonjatovo IMJ, Reed MT, Rigatti R, Rodighiero C, Ross MT, Sabot A, Sankar SV, Scally A, Schroth GP, Smith ME, Smith VP, Spiridou A, Torrance PE, Tzonev SS, Vermaas EH, Walter K, Wu X, Zhang L, Alam MD, Anastasi C, Aniebo IC, Bailey DMD, Bancarz IR, Banerjee S, Barbour SG, Baybayan PA, Benoit VA, Benson KF, Bevis C, Black PJ, Boodhun A, Brennan JS, Bridgham JA, Brown RC, Brown AA, Buermann DH, Bundu AA, Burrows JC, Carter NP, Castillo N, Chiara E. Catenazzi M, Chang S, Neil Cooley R, Crake NR, Dada OO, Diakoumakos KD, Dominguez-Fernandez B, Earnshaw DJ, Egbujor UC, Elmore DW, Etchin SS, Ewan MR, Fedurco M, Fraser LJ, Fuentes Fajardo KV, Scott Furey W, George D, Gietzen KJ, Goddard CP, Golda GS, Granieri PA, Green DE, Gustafson DL, Hansen NF, Harnish K, Haudenschild CD, Heyer NI, Hims MM, Ho JT, Horgan AM, Hoschler K, Hurwitz S, Ivanov DV, Johnson MQ, James T, Huw Jones TA, Kang G-D, Kerelska TH, Kersey AD, Khrebtukova I, Kindwall AP, Kingsbury Z, Kokko-Gonzales PI, Kumar A, Laurent MA, Lawley CT, Lee SE, Lee X, Liao AK, Loch JA, Lok M, Luo S, Mammen RM, Martin JW, McCauley PG, McNitt P, Mehta P, Moon KW, Mullens JW, Newington T, Ning Z, Ling Ng B, Novo SM, O'Neill MJ, Osborne MA, Osnowski A, Ostadan O, Paraschos LL, Pickering L, Pike AC, Pike AC, Chris Pinkard D, Pliskin DP, Podhasky J, Quijano VJ, Raczy C, Rae VH, Rawlings SR, Chiva Rodriguez A, Roe PM, Rogers J, Rogert Bacigalupo MC, Romanov N, Romieu A, Roth RK, Rourke NJ, Ruediger ST, Rusman E, Sanches-Kuiper RM, Schenker MR, Seoane JM, Shaw RJ, Shiver MK, Short SW, Sizto NL, Sluis JP, Smith MA, Ernest Sohna Sohna J, Spence EJ, Stevens K, Sutton N, Szajkowski L, Tregidgo CL, Turcatti G, vandeVondele S, Verhovsky Y, Virk SM, Wakelin S, Walcott GC, Wang J, Worsley GJ, Yan J, Yau L, Zuerlein M, Rogers J, Mullikin JC, Hurles ME, McCooke NJ, West JS, Oaks FL, Lundberg PL, Klenerman D, Durbin R, Smith AJ. 2008. Accurate whole human genome sequencing using reversible terminator chemistry. *Nature* **456**:53-59. doi: 10.1038/nature07517

Brick K, Pizzi E. 2008. A novel series of compositionally biased substitution matrices for comparing plasmodium proteins. *BMC Bioinformatics* **9**:236. doi: 10.1186/1471-2105-9-236

Claessens A, Hamilton WL, Kekre M, Otto TD, Faizullabhoy A, Rayner JC, Kwiatkowski D, Deitsch K. 2014. Generation of antigenic diversity in plasmodium falciparum by structured rearrangement of var genes during mitosis. *PLoS Genetics* **10**:e1004812. doi: 10.1371/journal.pgen.1004812

Dogovski C, Xie SC, Burgio G, Bridgford J, Mok S, McCaw JM, Chotivanich K, Kenny S, Gnädig N, Straimer J, Bozdech Z, Fidock DA, Simpson JA, Dondorp AM, Foote S, Klonis N, Tilley L, Schneider DS. 2015. Targeting the cell stress response of plasmodium falciparum to overcome artemisinin resistance. *PLOS Biology* **13**: e1002132. doi: 10.1371/journal.pbio.1002132

Dondorp AM, Nosten F, Yi P, Das D, Phyo AP, Tarning J, Lwin KM, Ariey F, Hanpithakpong W, Lee SJ, Ringwald P, Silamut K, Imwong M, Chotivanich K, Lim P, Herdman T, An SS, Yeung S, Singhasivanon P, Day NPJ, Lindegardh N, Socheat D, White NJ. 2009. Artemisinin resistance in plasmodium falciparum malaria. *New England Journal of Medicine* **361**:455-467. doi: 10.1056/NEJMoa0808859

Dondorp AM, Ringwald P. 2013. Artemisinin resistance is a clear and present danger. *Trends in Parasitology* **29**: 359-360. doi: 10.1016/j.pt.2013.05.005
Escalante AA, Ayala FJ. 1995. Evolutionary origin of plasmodium and other apicomplexa based on rRNA genes. *Proceedings of the National Academy of Sciences of the United States of America* **92**:5793-5797. doi: 10.1073/pnas.92.13.5793

Flegg JA, Guerin PJ, White NJ, Stepniewska K. 2011. Standardizing the measurement of parasite clearance in falciparum malaria: the parasite clearance estimator. *Malaria Journal* **10**:339. doi: 10.1186/1475-2875-10-339

Gardner KB, Sinha I, Bustamante LY, Day NPJ, White NJ, Woodrow CJ. 2011. Protein-based signatures of functional evolution in plasmodium falciparum. *BMC Evolutionary Biology* **11**:257. doi: 10.1186/1471-2148-11-257

Garud NR, Messer PW, Buzbas EO, Petrov DA, Copenhaver GP. 2015. Recent selective sweeps in north american drosophila melanogaster show signatures of soft sweeps. *PLOS Genetics* **11**:e1005004. doi: 10.1371/journal.pgen.1005004

Ghorbal M, Gorman M, Macpherson CR, Martins RM, Scherf A, Lopez-Rubio J-J. 2014. Genome editing in the human malaria parasite plasmodium falciparum using the CRISPR-Cas9 system. *Nature Biotechnology* **32**:819-821. doi: 10.1038/nbt.2925

Gower JC. 1966. Some distance properties of latent root and vector methods used in multivariate analysis. *Biometrika* **53**:325-338. doi: 10.1093/biomet/53.3.4-325

Hastings IM. 2003. Malaria control and the evolution of drug resistance: an intriguing link. *Trends in Parasitology* **19**:70-73. doi: 10.1016/S1471-4922(02)00017-X

Henikoff S, Henikoff JG. 1992. Amino acid substitution matrices from protein blocks. *Proceedings of the National Academy of Sciences of the United States of America* **89**:10915-10919. doi: 10.1073/pnas.89.22.10915

Hermisson J, Pennings PS. 2005. Soft sweeps: molecular population genetics of adaptation from standing genetic variation. *Genetics* **169**:2335-2352. doi: 10.1534/genetics.104.036947

Hien T, Thuy-Nhien N, Phu N, Boni MF, Thanh N, Nha-Ca N, Thai L, Thai C, Van Toi P, Thuan P, Long L, Dong L, Merson L, Dolecek C, Stepniewska K, Ringwald P, White NJ, Farrar J, Wolbers M. 2012. In vivo susceptibility of plasmodium falciparum to artesunate in binh phuoc province, vietnam. *Malaria Journal* **11**:355. doi: 10.1186/1475-2875-11-355

Hopkins Sibley C. 2015. Artemisinin resistance: the more we know, the more complicated it appears. *Journal of Infectious Diseases* **211**:667-669. doi: 10.1093/infdis/jiu469

Jensen JD. 2014. On the unfounded enthusiasm for soft selective sweeps. *Nature Communications* **5**:5281. doi: 10.1038/ncomms6281

Joy DA, Feng X, Mu J, Furuya T, Chotivanich K, Krettli AU, Ho M, Wang A, White NJ, Suh E, Beerli P, Su XZ. 2003. Early origin and recent expansion of plasmodium falciparum. *Science* **300**:318-321. doi: 10.1126/science.1081449

Kamau E, Campino S, Amenga-Etego L, Drury E, Ishengoma D, Johnson K, Mumba D, Kekre M, Yavo W, Mead D, Bouyou-Akotet M, Apinjoh T, Golassa L, Randrianarivelojosia M, Andagalu B, Maiga-Ascofare O, Amambua-Ngwa A, Tindana P, Ghansah A, MacInnis B, Kwiatkowski D, Djimde AA. 2015. K13-propeller polymorphisms in plasmodium falciparum parasites from sub-saharan africa. *The Journal of Infectious Diseases* **211**:1352-1355. doi: 10.1093/infdis/jiu608

Kyaw MP, Nyunt MH, Chit K, Aye MM, Aye KH, Aye MM, Lindegardh N, Tarning J, Imwong M, Jacob CG, Rasmussen C, Perin J, Ringwald P, Nyunt MM, Borrmann S. 2013. Reduced susceptibility of plasmodium falciparum to artesunate in southern myanmar. *PLoS ONE* **8**:e57689. doi: 10.1371/journal.pone.0057689

Kyte J, Doolittle RF. 1982. A simple method for displaying the hydropathic character of a protein. *Journal of Molecular Biology* **157**:105-132. doi: 10.1016/0022-2836(82)90515-0

Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, Valentin F, Wallace IM, Wilm A, Lopez R, Thompson JD, Gibson TJ, Higgins DG. 2007. Clustal w and clustal x version 2.0. *Bioinformatics* **23**:2947-2948. doi: 10.1093/bioinformatics/btm404

Lawson DJ, Hellenthal G, Myers S, Falush D, Copenhaver GP. 2012. Inference of population structure using dense haplotype data. *PLoS Genetics* **8**:e1002453. doi: 10.1371/journal.pgen.1002453

Li H, Durbin R. 2009. Fast and accurate short read alignment with burrows-wheeler transform. *Bioinformatics* **25**:1754-1760. doi: 10.1093/bioinformatics/btp324

Manske HM, Kwiatkowski DP. 2009. SNP-o-matic. *Bioinformatics* **25**:2434-2435. doi: 10.1093/bioinformatics/btp403

Manske M, Miotto O, Campino S, Auburn S, Almagro-Garcia J, Maslen G, O'Brien J, Djimde A, Doumbo O, Zongo I, Ouedraogo J-B, Michon P, Mueller I, Siba P, Nzila A, Borrmann S, Kiara SM, Marsh K, Jiang H, Su X-Z, Amaratunga C, Fairhurst R, Socheat D, Nosten F, Imwong M, White NJ, Sanders M, Anastasi E, Alcock D, Drury E, Oyola S, Quail MA, Turner DJ, Ruano-Rubio V, Jyothi D, Amenga-Etego L, Hubbart C, Jeffreys A, Rowlands K, Sutherland C, Roper C, Mangano V, Modiano D, Tan JC, Ferdig MT, Amambua-Ngwa A, Conway DJ, Takala-Harrison S, Plowe CV, Rayner JC, Rockett KA, Clark TG, Newbold CI, Berriman M, MacInnis B, Kwiatkowski DP. 2012. Analysis of plasmodium falciparum diversity in natural infections by deep sequencing. *Nature* **487**:375-379. doi: 10.1038/nature11174

Mbengue A, Bhattacharjee S, Pandharkar T, Liu H, Estiu G, Stahelin RV, Rizk SS, Njimoh DL, Ryan Y, Chotivanich K, Nguon C, Ghorbal M, Lopez-Rubio J-J, Pfrender M, Emrich S, Mohandas N, Dondorp AM, Wiest O, Haldar K. 2015. A molecular mechanism of artemisinin resistance in plasmodium falciparum malaria. *Nature* **520**:683-687. doi: 10.1038/nature14412

Messer PW, Petrov DA. 2013. Population genomics of rapid adaptation by soft selective sweeps. *Trends in Ecology & Evolution* **28**:659-669. doi: 10.1016/j.tree.2013.08.003
Miotto O, Almagro-Garcia J, Manske M, MacInnis B, Campino S, Rockett KA, Amaratunga C, Lim P, Suon S, Sreng S, Anderson JM, Duong S, Nguon C, Chuer CM, Saunders D, Se Y, Lon C, Fukuda MM, Amenga-Etego L, Hodgson AVO, Asoala V, Imwong M, Takala-Harrison S, Nosten F, Su Xin-zhuan, Ringwald P, Ariey F, Dolecek C, Hien TT, Boni MF, Thai CQ, Amambua-Ngwa A, Conway DJ, Djimdé AA, Doumbo OK, Zongo I, Ouedraogo J-B, Alcock D, Drury E, Auburn S, Koch O, Sanders M, Hubbart C, Maslen G, Ruano-Rubio V, Jyothi D, Miles A, O'Brien J, Gamble C, Oyola SO, Rayner JC, Newbold CI, Berriman M, Spencer CCA, McVean G, Day NP, White NJ, Bethell D, Dondorp AM, Plowe CV, Fairhurst RM, Kwiatkowski DP. 2013. Multiple populations of artemisinin-resistant plasmodium falciparum in cambodia. *Nature Genetics* **45**:648–655. doi: **10.1038/ng.2624**

Miotto O, Amato R, Ashley EA, MacInnis B, Almagro-Garcia J, Amaratunga C, Lim P, Mead D, Oyola SO, Dhorda M, Imwong M, Woodrow C, Manske M, Stalker J, Drury E, Campino S, Amenga-Etego L, Thanh T-NN, Tran HT, Ringwald P, Bethell D, Nosten F, Phyto AP, Pukrittayakamee S, Chotivanich K, Chuer CM, Nguon C, Suon S, Sreng S, Newton PN, Mayxay M, Khanthavong M, Hongvanthong B, Htut Y, Han KT, Kyaw MP, Faiz MA, Fanello CI, Onyamboko M, Mokuolu OA, Jacob CG, Takala-Harrison S, Plowe CV, Day NP, Dondorp AM, Spencer CCA, McVean G, Fairhurst RM, White NJ, Kwiatkowski DP. 2015. Genetic architecture of artemisinin-resistant plasmodium falciparum. *Nature Genetics* **47**:226–234. doi: **10.1038/ng.3189**

Mita T, Tanabe K, Kita K. 2009. Spread and evolution of plasmodium falciparum drug resistance. *Parasitology International* **58**:201–209. doi: **10.1016/j.parint.2009.04.004**

Mok S, Ashley EA, Ferreira PE, Zhu L, Lin Z, Yeo T, Chotivanich K, Imwong M, Pukrittayakamee S, Dhorda M, Nguon C, Lim P, Amaratunga C, Suon S, Hien TT, Htut Y, Faiz MA, Onyamboko MA, Mayxay M, Newton PN, Tripura R, Woodrow CJ, Miotto O, Kwiatkowski DP, Nosten F, Day NPJ, Preiser PR, White NJ, Dondorp AM, Fairhurst RM, Bozdech Z. 2015. Population transcriptomics of human malaria parasites reveals the mechanism of artemisinin resistance. *Science* **347**:431–435. doi: **10.1126/science.1260403**

Mu J, Awadalla P, Duan J, McGee KM, Joy DA, McVean GAT, Su Xin-zhuan, Hey J. 2005. Recombination hotspots and population structure in plasmodium falciparum. *PLoS Biology* **3**:e335. doi: **10.1371/journal.pbio.0030335**

Nair S, Nash D, Sudimack D, Jaidee A, Barends M, Uhlemann A-C, Krishna S, Nosten F, Anderson TJC. 2007. Recurrent gene amplification and soft selective sweeps during evolution of multidrug resistance in malaria parasites. *Molecular Biology and Evolution* **24**:562–573. doi: **10.1093/molbev/msl185**

Nielsen R, Hubisz MJ, Hellmann I, Torgerson D, Andres AM, Albrechtsen A, Gutenkunst R, Adams MD, Cargill M, Boyko A, Indap A, Bustamante CD, Clark AG. 2009. Darwinian and demographic forces affecting human protein coding genes. *Genome Research* **19**:838–849. doi: **10.1101/gr.088336.108**

Pennings PS, Hermisson J. 2006. Soft sweeps II–molecular population genetics of adaptation from recurrent mutation or migration. *Molecular Biology and Evolution* **23**:1076–1084. doi: **10.1093/molbev/msj117**

Phyo AP, Nkhoma S, Stepniewska K, Ashley EA, Nair S, McGready R, ler Moo C, Al-Saai S, Dondorp AM, Lwin KM, Singhasivanon P, Day NPJ, White NJ, Anderson TJC, Nosten F. 2012. Emergence of artemisinin-resistant malaria on the western border of thailand: a longitudinal study. *The Lancet* **379**:1960–1966. doi: **10.1016/S0140-6736(12)60484-X**

Prugnolle F, Durand P, Renaud F, Rousset F. 2010. Effective size of the hierarchically structured populations of the agent of malaria: a coalescent-based model. *Heredity* **104**:371–377. doi: **10.1038/hdy.2009.127**

Salgueiro P, Vicente JL, Ferreira C, Teófilo V, Galvão A, do Rosário VE, Cravo P, Pinto J. 2010. Tracing the origins and signatures of selection of antifolate resistance in island populations of plasmodium falciparum. *BMC Infectious Diseases* **10**:163. doi: **10.1186/1471-2334-10-163**

Saunders DL, Vanachayangkul P, Lon C, U.S. Army Military Malaria Research Program, National Center for Parasitology, Entomology, and Malaria Control (CNM), Royal Cambodian Armed Forces. 2014. Dihydroartemisinin–piperaquine failure in cambodia. *New England Journal of Medicine* **371**:484–485. doi: **10.1056/NEJMc1403007**

Straimer J, Gnadig NF, Witkowski B, Amaratunga C, Duru V, Ramadani AP, Dacheux M, Khim N, Zhang L, Lam S, Gregory PD, Urnov FD, Mercereau-Puijalon O, Benoit-Vical F, Fairhurst RM, Menard D, Fidock DA. 2015. K13-propeller mutations confer artemisinin resistance in plasmodium falciparum clinical isolates. *Science* **347**:428–431. doi: **10.1126/science.1260867**

Takala-Harrison S, Jacob CG, Arze C, Cummings MP, Silva JC, Dondorp AM, Fukuda MM, Hien TT, Mayxay M, Noedl H, Nosten F, Kyaw MP, Nhien NTT, Imwong M, Bethell D, Se Y, Lon C, Tyner SD, Saunders DL, Ariey F, Mercereau-Puijalon O, Menard D, Newton PN, Khanthavong M, Hongvanthong B, Starzengruber P, Fuehrer H-P, Swoboda P, Khan WA, Phyto AP, Nyunt MM, Nyunt MH, Brown TS, Adams M, Pepin CS, Bailey J, Tan JC, Ferdig MT, Clark TG, Miotto O, MacInnis B, Kwiatkowski DP, White NJ, Ringwald P, Plowe CV. 2015. Independent emergence of artemisinin resistance mutations among plasmodium falciparum in southeast asia. *Journal of Infectious Diseases* **211**:670–679. doi: **10.1093/infdis/jiu491**

Tanabe K, Mita T, Jombart T, Eriksson A, Horibe S, Palacpac N, Ranford-Cartwright L, Sawai H, Sakihama N, Ohmae H, Nakamura M, Ferreira MU, Escalante AA, Prugnolle F, Björkman A, Färnert A, Kaneko A, Horii T, Manica A, Kishino H, Balloux F. 2010. Plasmodium falciparum accompanied the human expansion out of Africa. *Current Biology* **20**:1283–1289. doi: **10.1016/j.cub.2010.05.053**

Taylor SM, Parobek CM, DeConti DK, Kayentao K, Coulibaly SO, Greenwood BM, Tagbor H, Williams J, Bojang K, Njie F, Desai M, Kariuki S, Gutman J, Mathanga DP, Martensson A, Ngasala B, Conrad MD, Rosenthal PJ, Tshefu AK, Moormann AM, Vulule JM, Doumbo OK, ter Kuile FO, Meshnick SR, Bailey JA, Juliano JJ. 2015. Absence of putative artemisinin resistance mutations among plasmodium falciparum in sub-saharan Africa: a molecular epidemiologic study. *Journal of Infectious Diseases* **211**:680–688. doi: **10.1093/infdis/jiu467**
Venkatesan M, Amaratunga C, Campino S, Auburn S, Koch O, Lim P, Uk S, Socheat D, Kwiatkowski DP, Fairhurst RM, Plowe CV. 2012. Using CF11 cellulose columns to inexpensively and effectively remove human DNA from plasmodium falciparum-infected whole blood samples. *Malaria Journal* **11**:41. doi: 10.1186/1475-2875-11-41

Witkowski B, Amaratunga C, Khim N, Sreng S, Chim P, Kim S, Lim P, Mao S, Sopha C, Sam B, Anderson JM, Duong S, Chuer CM, Taylor WRJ, Suon S, Mercereau-Puijalon O, Fairhurst RM, Menard D. 2013. Novel phenotypic assays for the detection of artemisinin-resistant plasmodium falciparum malaria in cambodia: in-vitro and ex-vivo drug-response studies. *The Lancet Infectious Diseases* **13**:1043–1049. doi: 10.1016/S1473-3099(13)70252-4

World Health Organization. 2014a. Update on Artemisinin Resistance - September 2014. Geneva: World Health Organization. http://www.who.int/entity/malaria/publications/atoz/status-rep-artemisinin-resistance-sep2014.pdf.

World Health Organization. 2014b. World Malaria Report 2014. Geneva: World Health Organization. http://www.who.int/malaria/publications/world_malaria_report_2014/en/.