# A novel field-based molecular assay to detect validated artemisinin-resistant k13 mutants

Laurence Vachot-Ganée, Nimol Khim, Alexandra Iannello, Eric Legrand, Saorin Kim, Rotha Eam, Chanra Khean, Malen Ken, Elizabeth Ashley, Kyaw Myo Tun, et al.

► To cite this version:

Laurence Vachot-Ganée, Nimol Khim, Alexandra Iannello, Eric Legrand, Saorin Kim, et al. A novel field-based molecular assay to detect validated artemisinin-resistant k13 mutants. Malaria Journal, 2018, 17 (1), pp.175. 10.1186/s12936-018-2329-y. pasteur-01813852

HAL Id: pasteur-01813852

https://pasteur.hal.science/pasteur-01813852

Submitted on 12 Jun 2018

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>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

# A novel field-based molecular assay to detect validated artemisinin-resistant *k13* mutants

Laurence Vachot-Ganée¹, Nimol Khim², Alexandra Iannello¹, Eric Legrand³,⁴, Saorin Kim², Rotha Eam², Chanra Khean², Malen Ken², Elizabeth Ashley⁵, Kyaw Myo Tun⁶,⁷, Mehul Dhorda⁵,⁸, François Nosten⁹,¹⁰, Issa Mahamat Souleymane¹¹, Sophie Blein¹, Alexandre Pachot¹, Frédéric Ariey¹²,¹³, Karine Kaiser¹ and Didier Ménard²,⁴*

**Abstract**

**Background:** Given the risk of artemisinin resistance spreading from the Greater Mekong sub-region, prospective monitoring in sub-Saharan Africa should be expedited. Molecular biology techniques used for monitoring rely on the detection of *k13* validated mutants by using PCR and Sanger sequencing approach, usually not available in malaria endemic areas.

**Methods:** A semi-automated workflow based on the easyMAG® platform and the Argene Solution® (bioMérieux, Marcy l'Etoile, France) as a field-based surveillance tool operable at national level was developed in four steps. Clinical and analytical performances of this tool detecting five of the most frequent and validated *k13* mutants (Y493H, I543T, R539T, F446I and C580Y) from dried blood spots (DBS) were compared to the gold standard approach (PCR and Sanger sequencing).

**Results:** By using the ARMS (amplification-refractory mutation system) strategy, the best multiplexing options were found in 3 separate real-time PCR duplexes (IC as internal control/I543T, C580Y/Y493H and F446I/R539T) with limits of detection ranging from 50 (C580Y) to 6.25 parasites/µL (Y493H). In field conditions, using 642 clinical DBS (from symptomatic patients and asymptomatic individuals) collected from Cambodia, Myanmar and Africa (Chad), the overall sensitivity and specificity of the K13 bMx prototype assay developed by bioMérieux were ≥ 90%. Areas under the ROC curves were estimated to be > 0.90 for all *k13* mutants in samples from symptomatic patients.

**Conclusion:** The K13 ready-to-use bMx prototype assay, considered by the end-users as a user-friendly assay to perform (in shorter time than the K13 reference assay) and easy to interpret, was found to require less budget planning and had fewer logistical constraints. Its excellent performance qualifies the prototype as a reliable screening tool usable in malaria endemic countries recognized to be at risk of emergence or spread of validated *k13* mutants. Additional multi-site studies are needed to evaluate the performances of the K13 bMx prototype assay in different epidemiological contexts such as Africa, India, or South America.

**Keywords:** Malaria, *Plasmodium falciparum*, Artemisinin resistance, *k13* mutation detection, Surveillance

**Background**

In 2016, an estimated 216 million (196–263 million) cases of malaria and 445,000 deaths were recorded worldwide,

mainly in sub-Saharan Africa (90%), followed by Southeast Asia (3%) and Eastern Mediterranean region (2%) [1]. Between 2010 and 2016, the incidence rate of malaria decreased by 18% globally [1]. This trend is commonly attributed to the widespread implementation of artemisinin-based combination therapy (an estimated 409 million treatment courses of ACT were procured in 2016)

<footnote>*Correspondence: dmenard@pasteur.fr
² Malaria Molecular Epidemiology Unit, Institute Pasteur in Cambodia, Phnom Penh, Cambodia
Full list of author information is available at the end of the article</footnote>
and, therefore, to the improvement of the management of falciparum uncomplicated malaria cases [2]. ACT combined a potent and fast-acting artemisinin derivative that reduces rapidly the parasite biomass with a long half-life partner drug that eliminates any remaining parasites [3].

Unfortunately, in 2007–2008, partial resistance to artemisinin derivatives (ART-R) was detected in Southeast Asia (first in Cambodia and later in Myanmar, Thailand, Vietnam and Lao PDR) followed more recently by the emergence of resistance to piperaquine (PPQ-R) [4–6]. At present, the reduced efficacy of ART and PPQ translates into late treatment failures and prolonged parasite carriage, thereby increasing the transmission potential of multidrug resistant infections [7–11]. Dihydroartemisinin–piperaquine (DHA–PPQ) treatment failures are estimated to reach 60% in some Cambodian provinces, indicating the dramatic spread of both ART-R and PPQ-R across Cambodia and beyond. Indeed, recent genomic investigations have demonstrated that *Plasmodium falciparum* ART-R and PPQ-R have spread from western Cambodia to neighbouring provinces in Cambodia and then to north-eastern Thailand, southern Laos and southern Vietnam, leading to a substantial increase in DHA–PPQ failure rates. This rapid expansion was clearly shown to be associated to the hard-selective sweep of a single fit ART-R *P. falciparum* parasite lineage, outcompeting the other resistant malaria parasites, and subsequently acquiring PPQ-R [9, 12].

Since 2014, non-synonymous mutations in the propeller domain of a *kelch* gene located on the chromosome 13 (*k13*, PF3D7_1343700) have been shown to be a major determinant of ART-R [13–15]. Currently, the molecular surveillance of the global extent of ART-R malaria parasites relies on these markers [16]. To date, almost 200 *k13* mutations have been described worldwide. In Southeast Asia, *k13* mutants are frequent and distributed between Cambodia–Vietnam–Lao PDR (C580Y, R539T, Y493H and I543T) and western Thailand–Myanmar–China (C580Y, F446L, N458Y, P574L and R561H) with almost no overlap (except C580Y and P553L which has been found in the two areas) [17]. The C580Y mutant is largely predominant in Southeast Asia. In Africa, a broad array of rare non-synonymous mutations is observed at very low frequency. Clinical and biological investigations have shown that not all non-synonymous *k13* mutations confer ART-R [17]. To be validated as ART-R molecular marker, a *k13* mutant has to be correlated with delayed parasite clearance in clinical studies and reduced drug in vitro susceptibility (survival rate ≥ 1% expressed by the ring-stage survival assay, RSA⁰⁻³h) in fresh isolates (ex vivo assays) or culture-adapted field parasites or K13 genome-edited parasites (in vitro assays) [18]. According to the latest WHO update on artemisinin resistance

(released in April 2017), only six *k13* mutants are validated (C580Y, Y493H, R539T, I543T, N458Y, R561H) [16]. The F446I mutant, which is highly prevalent in Myanmar, is strongly suspected of being associated to ART-R [17, 19–22].

In such an epidemiological context, the risk of ART-R parasites spreading from the Greater Mekong sub-region (GMS) to Africa, as happened previously with chloroquine and sulfadoxine/pyrimethamine-resistant parasites, is a major concern threatening the world's malaria elimination efforts [23]. This alarming situation mandates close surveillance and containment of the dissemination of ART-R parasites in Southeast Asia, but above all implementation of proactive monitoring of anti-malarial drug resistance in sub-Saharan Africa. The emergence and the spread of ART-R outside Southeast Asia can now be easily tracked in real-time by the detection of *k13* validated mutants in dried blood spots (DBS) collected from finger prick [24]. Molecular biology techniques used for ART-R monitoring rely mainly on polymerase chain reaction (PCR) and Sanger sequencing. However, sequencing platforms are not available in most malaria endemic areas, thus obligating the shipment of the samples or the PCR products to appropriate facilities. Substantial work remains to be done to harmonize and validate collected data. This includes the development of field-based molecular biology techniques, the provision of control samples, or the implementation of a quality assurance system [25].

To overcome this issue, a semi-automated workflow based on the easyMAG® platform and the Argene Solution® (bioMérieux, Marcy l'Etoile, France) as a field-based surveillance tool operable at national level was developed and evaluated. This tool (referred as K13 bMx prototype) allows the detection of the five most frequent and validated ART-R *k13* mutants from DBS and can be deployed in most laboratories in malaria endemic areas easily (Additional file 1). The clinical and analytical performances of the K13 bMx prototype were compared to the gold standard approach (QiaAmp DNA blood mini kit extraction followed by PCR Sanger sequencing referred to as K13 reference assay) [15] and are presented in this manuscript.

## Methods

### K13 bMx prototype: development and validation

The prototype was developed in four sequential steps.

**Step 1** The easyMAG® DNA extraction protocol from DBS was optimized. First, a new format of 3 MM filter paper was designed to facilitate the semi-automated extraction protocol and avoid inter-sample DNA contamination (Additional file 2). Second, the limit of detection (LoD), the repeatability and the reproducibility of the easyMAG® DNA extraction protocol were
assessed by using artificial DBS (prepared from the 3D7 culture-adapted line at parasite densities ranging from 5 to 0.6 parasite/µL). Real-time PCR assay targeting the *cytochrome c oxidase subunit 1* gene (PlasmoDB gene ID: mal_mito_2 cytochrome c oxidase subunit 1) was used to evaluate the extraction efficiency. This PCR assay was thereafter integrated as internal control (IC) in the K13 bMx prototype. Third, results of real-time PCR assay targeting the *Plasmodium cytochrome b* gene (PlasmoDB gene ID: mal_mito_3 cytochrome b) using DNA extracted with the QiaAmp DNA blood mini kit (Qiagen, Hilden, Germany) and real-time PCR assay targeting the *cytochrome c oxidase subunit 1* gene using DNA extracted with the easyMAG® DNA extraction protocol were compared. DNA were obtained from *P. falciparum* blood samples (clinical *falciparum* malaria cases with parasitaemia ranging from 0.1 to 2% and asymptomatic individuals harbouring submicroscopic *P. falciparum* PCR positive infection with a parasite density<0.01%) as well as negative blood samples spotted onto filter paper (50 µL of blood).

**Step 2** The analytical performances of single and multiplex real-time PCR assays designed to detect C580Y, Y493H, R539T, I543T, F446I *k13* mutants (PlasmoDB gene ID: PF3D7_1343700) and the *cytochrome c oxidase subunit 1* gene (used as internal control, IC) were assessed. First, in silico strategies were used to design specific PCR primers. Second, single and then multiplex real-time PCRs for each primer set were optimized to reach the best analytical performances (LoD) by using gDNA extracted from *P. falciparum* culture-adapted isolates harbouring different *k13* alleles (wild-type, C580Y, Y493H, R539T, I543T, F446I, parasite densities ranging from 4000 to 0 parasites/µL) or from DBS containing 50 µL blood of *P. falciparum* culture-adapted isolates containing decreasing parasite densities of each *k13* mutant (from 100 to 6.25 parasites/µL) extracted with the easyMAG®.

**Step 3** The performances of the K13 bMx prototype (from DBS to results) were assessed at bioMerieux (Grenoble, France). The LoD, repeatability and reproducibility were then evaluated by using artificial DBS (obtained from wild-type and *k13* mutants culture-adapted lines at parasite densities ranging from 500 to 6.25 parasites/µL). The specificity (absence of cross reactivity with other *Plasmodium* species) was assessed by using pure DNA from *P. vivax*, *Plasmodium ovale* and *Plasmodium malariae* (from DBS at 4000, 100 and 10 parasites/µL per species) or DNA mixed at the same parasite densities with *P. falciparum* 3D7 DNA (4000 parasites/µL). Finally, the performances were evaluated with DBS collected in Cambodia and Myanmar from malaria cases (symptomatic patients and asymptomatic individuals)

or malaria-free individuals. K13 reference assay used as gold standard method was performed at Institute Pasteur Cambodia (IP Cambodia, Phnom Penh, Cambodia) on the same specimens. The technicians were blinded to results from the other method.

**Step 4** The final performances of the K13 bMx prototype used in field conditions were evaluated at IP Cambodia. Tested DBS samples were obtained in Cambodia, Myanmar and Africa from symptomatic patients, asymptomatic and malaria-free individuals. Molecular tests (K13 bMx prototype and K13 reference assays) were performed independently by technicians with no knowledge of each other's results.

## Samples

For steps 1, 2 and 3, a set of *P. falciparum* isolates harbouring different *k13* alleles (C580Y, Y493H, R539T, I543T, F446I) and the 3D7 reference line (*k13* wild-type) were adapted for long-term in vitro culture. Culture adaptation of the field Cambodian isolates was performed as described previously [26]. Artificial DBS specimens, mimicking a whole-blood sample obtained from patients, were prepared from culture-adapted parasites by spotting 50 µL blood (50% haematocrit) with decreasing parasite densities of each *k13* wild-type and mutants onto the 3 MM filter paper. Parasite-free blood samples were obtained from the blood bank (Phnom Penh, Cambodia) and were used to prepare artificial malaria-free DBS.

For steps 1, 3 and 4, blood samples (284 and 642 clinical DBS for steps 3 and 4, respectively) were obtained from febrile patients suspected of having malaria and seeking anti-malarial treatment in health facilities or in villages (through the village malaria workers network) across Cambodia, Myanmar and Chad. Malaria diagnosis was achieved in the field either by microscopy of Giemsa-stained malaria blood films or by malaria rapid diagnostic test (CareStart Malaria HRP2/pLDH Pf/PAN Combo, Access Bio, USA) detecting both *P. falciparum* and non-*P. falciparum* infections. After obtained written informed consent, DBS were prepared from finger prick (using the new format of 3 MM filter paper or the Whatman FTA card). Fresh capillary blood samples spotted onto the 3 MM filter paper were obtained from asymptomatic individuals living in villages in Rattanakiri and Preah Vihear province (Cambodia) during cross sectional surveys conducted in 2016–2017 as previously described [27, 28].

## Parasite DNA extraction

Parasite DNA was extracted at IP Cambodia from red blood cell pellets obtained from culture-adapted isolates using the QiaAmp DNA blood mini kit (Qiagen) following manufacturer's recommendations. Parasite DNA
was extracted from clinical or artificial DBS using the QiaAmp DNA blood mini kit (Qiagen) at IP Cambodia or the easyMAG® system in both sites (IP Cambodia and bioMérieux), according to the supplier's recommendations. DNA obtained from the easyMAG® Generic 2.0.1 protocol was extracted using 2 mL lysis buffer with 50 µL of silica and eluted in a volume of 50 µL.

## Detection of Plasmodium parasites and k13 non-synonymous mutants

At IP Cambodia, *Plasmodium* DNA was detected using a two-step real-time PCR assay targeting the *Plasmodium cytochrome b* gene, as previously described [15]. Briefly, each blood sample was screened using genus-specific primers for malaria infection. Only positive samples were tested for *Plasmodium* species identification using species-specific primers for the four main human malaria species (*Plasmodium vivax*, *P. falciparum*, *P. ovale*, *P. malariae*). In the K13 reference assay, amplification by nested-PCR of the K13-propeller domain (codons 443–666, i.e. 720 bp) was used to detect *k13* mutations by sequencing, according to the operating procedures previously developed [15, 17]. PCR products were sequenced by Macrogen (Seoul, South Korea). Electrophoregrams were analysed on both strands with CEQ 2000 genetic analysis system software (Beckman Coulter, Villepinte, France), using PF3D7_1343700 as reference sequence. Isolates with mixed alleles were considered as mutants.

In the K13 bMx prototype assay, *P. falciparum* DNA (IC, targeting the *cytochrome c oxidase subunit 1* gene) and *k13* mutants (C580Y, Y493H, R539T, I543T, F446I) were detected in 3 separate real-time PCR duplexes (IC/ I543T, C580Y/Y493H and F446I/R539T) using FAM and HEX fluorophores. Real-time PCR duplexes reactions were done with 15 µL of ready-to-use reaction mix and 10 µL of DNA extract and amplified on the CFX96 PCR instrument (Bio-Rad, Marnes-la-Coquette, France). PCR amplifications were performed under the following conditions: heating at 95 °C for 15 min, followed by 45 cycles of heating at 95 °C for 10 s, and of annealing/ extension at 60 °C for 40 s. For each run, positive (gDNA from culture-adapted Cambodian *P. falciparum* parasites harbouring each of the *k13* mutant alleles and from the 3D7 line) and no template (elution buffer) controls were included.

The “automatic threshold mode” was applied to define Cq and fluorescence values (usually set up at ~10% of the average fluorescence). A specimen was considered positive for the tested mutation if the delta Cq (Cq mutant PCR–Cq IC PCR) was below a specific cut-off. This cut-off was set at 6 for I543T mutant, 9 for R539T and F446 mutants and 11 for Y493H and C580Y mutants.

## Statistical analysis

Data were recorded and analysed using Excel software and MedCalc (MedCalc Software, Belgium). The F-test was used to compare variances of Cq values.

To assess the analytical performance of the K13 bMx prototype assay (vs. the K13 reference assay), standard diagnostic test measures were determined. Sensitivity (Se) was the proportion of DBS that the K13 bMx prototype assay classified as mutant among DBS classified as *k13* mutant (C580Y, Y493H, R539T, I543T, or F446I) with the K13 reference assay (true positive rate). Specificity (Sp) was the proportion of DBS that the K13 bMx prototype assay classified as *k13* non-mutant (C580Y, Y493H, R539T, I543T, or F446I) among DBS classified as *k13* non-mutant (included in the panel assay) with the K13 reference assay (true negative rate). Positive predictive value (PPV) was the probability that a *k13* mutant included in the panel assay was present in DBS when the K13 bMx prototype assay classified DBS as *k13* mutant and negative predictive value (NPV), the probability that a *k13* mutant included in the panel assay was not present in DBS when the K13 bMx prototype assay classified DBS as no mutant.

PPV and NPV were calculated based on the proportion of *k13* mutants (C580Y, Y493H, R539T, I543T, or F446I, referred as PP) found in the set of tested DBS as following:

$$ 
\begin{aligned}
\text{PPV} &= \text{Se} \times \text{PP}/\text{Se} \times \text{PP} + (1 - \text{Sp}) \times (1 - \text{PP}) \\
\text{and NPV} &= \text{Sp} \times (1 - \text{PP})/\text{Sp} \times (1 - \text{PP}) \\
&\quad + (1 - \text{Se}) \times \text{PP}.
\end{aligned}
$$

The AURoC values (Area under a ROC curve) were estimated to assess the accuracy of the K13 bMx prototype assay to detect the five *k13* mutants included in the panel, based on the results of the K13 reference assay. An AUC value of 1 represents a perfect test (i.e. a test that would classify all screened mutants exactly as the reference test) and an AUC value of 0.5 represents a worthless test. AUC value > 0.90 was considered to represent a diagnostic test with an excellent accuracy.

The 95% confidence intervals (CI 95%) for sensitivity and specificity are “exact” Clopper–Pearson confidence intervals. CI 95% for the PPV and NPV are the standard logit confidence intervals given by Mercaldo et al. [29]. The CI 95% for AUC values were used to test the hypothesis that the theoretical area is 0.5. All reported P-values are two-sided and were considered statistically significant if < 0.05.

## Results

### Performances of the easyMAG® extraction protocol (step 1)
The Cq values obtained from artificial DBS (3D7 line at 5–0.6 parasites/µL) are presented in Fig. 1. The LoD of the real-time PCR assay targeting the *cytochrome c*
<figure><img src="image_3.png" /><figcaption>**Fig. 1** Performances of the easyMAG® extraction protocol: Cq values from real-time PCR assay targeting the *cytochrome c oxidase subunit 1* gene with DNA extracted from artificial DBS (3D7 line at 5–0.6 parasites/µL)</figcaption></figure>

*cytochrome c oxidase subunit 1* gene using DNA extracted with the easyMAG® DNA extraction protocol was < 0.6 parasite/ µL. Runs assessing the reproducibility and the repeatability showed no significant difference of the median Cq values according to the operators (N=2, P=0.47, F-test), batches of reagents (N=2, P=0.93, F-test) or the easyMAG® system (N=2, P=0.73, F-test).

Seventy clinical DBS were tested both at IP Cambodia (using QiaAmp DNA blood mini kit Qiagen/two-step real-time PCR assays targeting the *Plasmodium cytochrome b* gene) as previously described [30, 31] and at bioMérieux (using the easyMAG® system and the real-time PCR targeting the *cytochrome c oxidase subunit 1* gene). The proportion of concordant results between both methods was very high (64/70 DBS, 91%), more frequent with DBS collected from symptomatic patients (29/30 DBS, 97%) than DBS collected from asymptomatic individuals (35/40 DBS, 88%) (Additional file 3). Six samples were found to have discordant results: one asymptomatic case was detected as *P. falciparum* infection at IP Cambodia and negative at bioMérieux; three asymptomatic cases were detected as negative at IP Cambodia and as *P. falciparum* infection at bioMérieux and two cases (one symptomatic and one asymptomatic) were detected as *P. vivax* infections at IP Cambodia and *P. falciparum* infections at bioMérieux.

## Development of multiplex real-time PCR assays to detect ART-R validated *k13* mutations (step 2)

Two molecular strategies were explored to detect the presence of SNPs (single nucleotide polymorphisms): TaqMan using mutant specific probes and Amplification-refractory mutation system (ARMS). The ARMS technique, based on allele-specific primers [32] (Additional file 4) was found to be the best option (high multiplexing capacity). First, allele-specific ARMS primers were designed and tested to detect each *k13* mutations in simplex reactions. Second, multiplex strategies by mixing simplex reactions were tested and optimized (testing different concentrations of MgCl₂, dNTP and primers/probe) using gDNA extracted from *k13* wild-type and mutant culture-adapted

<table><caption>**Table 1** Summary of the results obtained for the multiplex strategies by mixing simplex reactions</caption><thead><tr><th>Duplexes</th><th colspan="6">Target 2</th></tr><tr><th>Target 1</th><th>IC</th><th>F446I</th><th>Y493H</th><th>R539T</th><th>I543T</th><th>C580Y</th></tr></thead><tbody><tr><th>IC</th><td></td><td>Failed<sup>a</sup></td><td>Failed<sup>a</sup></td><td>OK*</td><td>OK</td><td>OK</td></tr><tr><th>F446I</th><td></td><td></td><td>Failed<sup>b</sup></td><td>OK</td><td>Failed<sup>b</sup></td><td>OK*</td></tr><tr><th>Y493H</th><td></td><td></td><td></td><td>Not tested</td><td>Not tested</td><td>OK</td></tr><tr><th>R539T</th><td></td><td></td><td></td><td></td><td>Failed<sup>d</sup></td><td>Failed<sup>c</sup></td></tr><tr><th>I543T</th><td></td><td></td><td></td><td></td><td></td><td>Failed<sup>b</sup></td></tr></tbody></table>

<footnote>* Of note, similar performances with those selected were obtained with 446/580 and IC/539 duplexes, demonstrating the multiplexing flexibility of the ARMS strategy
Failed means: <sup>a</sup> decrease of sensitivity to detect both SNPs, <sup>b</sup> presence of non-specific amplification curves, <sup>c</sup> decrease of sensitivity to detect one SNP (R539T), <sup>d</sup> overlap of the primers set</footnote>
*P. falciparum* isolates. The best analytical performances were observed in 3 separate real-time PCR duplexes (IC/ I543T, C580Y/Y493H and F446I/R539T) using FAM and HEX fluorophores, as presented in Table 1.

## Analytical performances of the K13 bMx prototype assay (from DBS to results) performed in controlled condition (step 3)

No signal of amplification was obtained with pure *P. vivax*, *P. ovale*, *P. malariae* DNA or mixed DNA with *P. falciparum* 3D7 (K13 wild-type). The LoDs of each duplex was investigated by using DNA extracted by the easyMAG® system from artificial DBS containing decreasing parasite densities of each *k13* mutant (from

100 to 6.25 parasites/µL). The LoDs (>95% detection in replicates) ranged from 50 (C580Y) to 25 (I543T, R539T, F446I) or 6.25 parasites/µL (Y493H). Data are shown in the Fig. 2a. The ranges of the mean Cq values according to the parasite density of *k13* mutants are displayed in the Fig. 2b (from 33.8 to 37.2 at 100 parasites/µL to 38.1–39.8 at 6.25 parasites/µL for the K13 mutant gene detection).

The reproducibility and the repeatability of the C580Y/ Y493H and F446I/R539T duplexes were investigated. No significant differences in mean Cq values were observed according to the operators (N=2, P=0.82, F-test), days when the runs were performed (N=2, P=0.17, F-test), the easyMAG® system used (N=2, P=0.83, F-test), or the batches of reagents tested (N=2, P=0.82, F-test).

<figure><img src="image_4.png" /><figcaption>**Fig. 2** a Limits of detection (% of positive replicates, 20 replicates per parasite density for K13 mutant) of each duplex obtained from DNA extracted by the easyMAG® system (Generic 2.0.1 protocol) from artificial DBS containing decreasing parasite densities of each K13 mutant (from 100 to 6.25 parasites/µL). b Mean Cq values determined according to the parasite density of *k13* mutant (from 100 to 6.25 parasites/µL)</figcaption></figure>
<table><caption>**Table 2** Evaluation of the analytical performances of the whole process (from DBS to results) of the K13 bMx prototype using 284 DBS collected in Cambodia and Myanmar and the K13 reference assay as gold standard method</caption><thead><tr><th rowspan="2">Assay</th><th rowspan="2">K13 allele</th><th colspan="8">K13 reference assay</th><th rowspan="2">Negative</th><th rowspan="2">Total</th></tr><tr><th>Wild-type</th><th>Other K13 alleles *</th><th>F446 I</th><th>Y493H</th><th>R539T</th><th>I543T</th><th>C580Y</th></tr></thead><tbody><tr><td rowspan="5">K13 bMx prototype</td><td>No mutant</td><td>100</td><td>7</td><td></td><td></td><td></td><td></td><td>1 + 1</td><td></td><td>109</td></tr><tr><td>F446I</td><td></td><td></td><td>4</td><td></td><td></td><td></td><td></td><td>1</td><td>5</td></tr><tr><td>Y493H</td><td></td><td></td><td></td><td>6</td><td></td><td></td><td></td><td></td><td>6</td></tr><tr><td>R539T</td><td></td><td></td><td></td><td></td><td>4</td><td></td><td></td><td></td><td>4</td></tr><tr><td>I543T</td><td></td><td></td><td></td><td></td><td></td><td>4</td><td></td><td></td><td>4</td></tr><tr><td rowspan="2">Total</td><td></td><td>100</td><td>7</td><td>5</td><td>6</td><td>4</td><td>4</td><td>125</td><td>33</td><td>284</td></tr></tbody></table>

Cells coloured in green displayed concordant results and those coloured in grey discordant results. Numbers in black, green, red and blue fonts are results obtained from DBS collected from *P. falciparum* symptomatic malaria cases, *P. falciparum* asymptomatic carriers, *P. vivax* asymptomatic carriers, malaria-free individuals, respectively

<table><caption>**Table 3** Origins of the DBS used to define the clinical performances of the K13 bMx prototype (from DBS to results) in field conditions</caption><thead><tr><th rowspan="2">Plasmodium species</th><th rowspan="2">Clinical status</th><th colspan="2">Asia</th><th colspan="2">Africa</th><th rowspan="2">Total</th></tr><tr><th>Cambodia</th><th>Myanmar</th><th></th><th></th></tr></thead><tbody><tr><td>–</td><td>Malaria-free</td><td>10</td><td></td><td></td><td></td><td>10</td></tr><tr><td rowspan="2">P. falciparum</td><td>Symptomatic</td><td>296</td><td>221</td><td>31</td><td></td><td>548</td></tr><tr><td>Asymptomatic</td><td>63</td><td></td><td></td><td></td><td>63</td></tr><tr><td rowspan="2">P. falciparum/P. vivax</td><td>Symptomatic</td><td>17</td><td></td><td></td><td></td><td>17</td></tr><tr><td>Asymptomatic</td><td>2</td><td></td><td></td><td></td><td>2</td></tr><tr><td rowspan="2">P. vivax</td><td>Symptomatic</td><td>1</td><td></td><td></td><td></td><td>1</td></tr><tr><td>Asymptomatic</td><td>1</td><td></td><td></td><td></td><td>1</td></tr><tr><td>Total</td><td></td><td>390</td><td>221</td><td>31</td><td></td><td>642</td></tr></tbody></table>

A total of 284 clinical DBS were tested both at IP Cambodia (K13 reference assay) and at bioMérieux (K13 bMx prototype assay). The overall concordance between assays was excellent (97.5%, CI 95% 96.9–99.8%). Negative (no mutant) and positive (k13 mutants included in the panel assay) concordances were 100% (CI 95% 96.5–100%) and 98.4% (CI 95% 94.3–99.6%), respectively. Only results from 7 DBS were found to be discordant (2.5%). The K13 bMx prototype misclassified two C580Y mutant samples as wild-type, four *P. falciparum* negative samples as F446I mutant (N=1) or C580Y mutants (N=3) and one F446I mutant sample as negative. Of note, seven *P. falciparum* samples harbouring k13 mutations not included in the panel (A626E, C469F, G449A, P553L, P574L and Y511H) were correctly classified as no mutant by the K13 bMx prototype. All data are presented in Table 2.

Clinical and analytical performances of the K13 bMx prototype (from DBS to results) performed in field conditions (step 4)

To evaluate the performances of the K13 bMx prototype in field conditions, 642 clinical DBS were tested with both assays at IP Cambodia. DBS included 390 specimens from Cambodia, 221 specimens from Myanmar and 31 specimens from Africa (Table 3).

Three technicians at IP Cambodia conducted the K13 bMx prototype assay and were blinded to the results of the K13 reference assay. Sequencing data were available for 606 samples, as presented in Table 4. The mutant I543T was not found in the set of clinical DBS. Among the 36 samples with no K13 sequencing valid data (10 malaria-free, 14 *P. falciparum* and 1 *P. vivax* infections from symptomatic patients and 10 *P. falciparum* and 1 *P. vivax* infections from asymptomatic individuals), the K13 bMx prototype assay detected 8 C580Y and 2 F446I mutants from *P. falciparum* symptomatic patients and 11 no mutant (4 from *P. falciparum* symptomatic patients, 5 from *P. falciparum* asymptomatic individuals, 1 from *P. vivax* symptomatic patient, 1 from *P. vivax* asymptomatic individual). The final analysis included only the 606 samples with available sequencing data.

The proportion of k13 wild-type and mutants (PPs) among the 606 samples with available sequencing data were: wild-type or k13 mutations not included in the panel assay (39.3%, 238/606), C580Y (42.2%, 256/606), F446I (11.7%, 71/606), R539T (3.5%, 21/606) and Y493H (3.3%, 20/606).

The overall concordance between results of both assays was 90.1% (546/606, CI 95% 82.7–98.0%), ranging from 90.0% (18/20, CI 95% 53.3–100%) (Y493H) to 95.8%

<footnote>* Other k13 alleles include A626E, C469F, G449A, P553L, P574L and Y511H mutants</footnote>
<table><caption>Table 4 K13 sequencing data of 642 DBS tested with the K13 reference assay in field conditions</caption><thead><tr><th rowspan="2">Sites</th><th rowspan="2">K13 reference assay</th><th colspan="4">Plasmodium species</th><th rowspan="2">Total</th></tr><tr><th>P. falciparum</th><th>P. falciparum/P. vivax</th><th>P. vivax</th><th>Negative</th></tr></thead><tbody><tr><td rowspan="12">ASIA</td><td>No data</td><td>21</td><td>3</td><td>2</td><td>10</td><td>36</td></tr><tr><td>C580Y</td><td>244</td><td>12</td><td></td><td></td><td>256</td></tr><tr><td>F446I</td><td>71</td><td></td><td></td><td></td><td>71</td></tr><tr><td>R539T</td><td>21</td><td></td><td></td><td></td><td>21</td></tr><tr><td>Y493H</td><td>20</td><td></td><td></td><td></td><td>20</td></tr><tr><td>G538 V</td><td>1</td><td></td><td></td><td></td><td>1</td></tr><tr><td>N458Y</td><td>8</td><td></td><td></td><td></td><td>8</td></tr><tr><td>P553L</td><td>2</td><td></td><td></td><td></td><td>2</td></tr><tr><td>P574L</td><td>1</td><td></td><td></td><td></td><td>1</td></tr><tr><td>R561H</td><td>13</td><td></td><td></td><td></td><td>13</td></tr><tr><td>S459L</td><td>1</td><td></td><td></td><td></td><td>1</td></tr><tr><td>WT</td><td>177</td><td>4</td><td></td><td></td><td>181</td></tr><tr><td>Total</td><td>580</td><td>19</td><td>2</td><td>10</td><td>611</td></tr><tr><td rowspan="3">AFRICA</td><td>V510V</td><td>1</td><td></td><td></td><td></td><td>1</td></tr><tr><td>WT</td><td>30</td><td></td><td></td><td></td><td>30</td></tr><tr><td>Total</td><td>31</td><td>0</td><td>0</td><td>0</td><td>31</td></tr></tbody></table>

(68/71, CI 95% 74.4–100%) (F446I) (Additional file 5). The concordance was higher using DBS from *P. falciparum* symptomatic patients (93.6%, 516/551, CI 95% 85.7–100.0%) compared to that obtained with DBS from *P. falciparum* asymptomatic individuals (54.5%, 30/55, CI 95% 36.8–77.9%, P = 0.02).

By excluding the 16 DBS with low DNA yield (DBS classified as no valid data with the K13 bMx prototype assay), the final performances of the K13 bMx prototype assay were excellent, as shown in the Table 5. The sensitivity and specificity of the K13 bMx prototype assay were ≥ 90% (except for the detection of the C580Y in samples from asymptomatic individuals, sensitivity = 62.8%). Similarly, PPV and NPV were > 86% (except for the detection of wild-type in samples from asymptomatic individuals, PPV = 16.7%). As expected, due to the higher amount of DNA, the performances of K13 bMx prototype assay with DBS from *P. falciparum* symptomatic patients were better compared to those obtained with DBS collected from the *P. falciparum* asymptomatic individuals (Table 5). The areas under the ROC curves were estimated to be > 0.93 for all *k13* mutants in samples from symptomatic patients (Fig. 3). Raw data are presented in Additional files 5, 6 and 7.

## Discussion

Over a period of 18 months, the K13 bMx prototype was designed, optimized, developed and finally evaluated using the classical PCR and Sanger sequencing approach (K13 reference assay) as gold standard [15, 17]. The main

objective of this project was to propose a robust and accurate field-based molecular assay to detect the five most prevalent validated *k13* mutations in an adapted throughput format. Results for 24 samples could be obtained in less than 4 h (including DNA extraction and real-time PCR runs) (Additional file 1). Quality assurance of the whole process was ensured by the semi-automated nucleic acid extraction platform from a new format of DBS (the most common matrix for blood sample collection used in surveillance surveys) and the inclusion of an internal control in the K13 bMx prototype assay.

Overall, the analytical and clinical performances of the K13 bMx prototype obtained from field clinical samples collected from symptomatic patients were excellent both in controlled (carried out at bioMerieux, Grenoble, France) and field conditions (carried out at IP Cambodia, Phnom Penh, Cambodia). However, those obtained from samples collected from asymptomatic individuals were largely lower. For instance, the K13 bMx prototype gave a result for almost all DBS from symptomatic patients tested in field conditions (544/551, 98.7% compared to 30/55, 54.5%, for DBS from asymptomatic individuals). The *k13* mutants were correctly detected: for C580Y in 95.8% (204/213 symptomatic patients) and 62.8% (27/43 asymptomatic individuals), for F446I in 95.8% (68/71), for the R539T in 95.2% (20/21) and for the Y493H in 90.0% (18/20). Wild-type allele or *k13* mutants not included in the panel assay were correctly classified as ‘no mutant’ in 91.1% (206/226) and 25.0% (3/12) in DBS collected from symptomatic patients and asymptomatic individuals, respectively.
<table><caption>Table 5 Summary of the performances of K13 bMx prototype according to the origin of the samples (all samples, N = 590; samples collected from *P. falciparum* symptomatic patients, N = 544; and samples collected from the *P. falciparum* asymptomatic individuals, N = 46)</caption><thead><tr><th rowspan="2">Samples from</th><th rowspan="2">Performances</th><th colspan="4">K13 mutants</th><th rowspan="2">K13 wild-type</th></tr><tr><th>C580Y</th><th>F446I</th><th>R539T</th><th>Y493H</th></tr></thead><tbody><tr><td rowspan="6">All samples</td><td>Sensitivity (CI 95%)</td><td>90.2%<br>(85.9–93.6)</td><td>95.8%<br>(88.1–99.1)</td><td>95.2%<br>(76.1–99.9)</td><td>90.0%<br>(68.3–98.8)</td><td>94.1%<br>(90.2–96.8)</td></tr><tr><td>Specificity (CI 95%)</td><td>99.4%<br>(97.6–99.9)</td><td>98.0%<br>(96.2–99.0)</td><td>100%<br>(99.4–100)</td><td>99.8%<br>(99.0–99.9)</td><td>92.4%<br>(89.2–94.9)</td></tr><tr><td>AUC<br>(CI 95%)</td><td>0.948<br>(0.927–0.965)</td><td>0.968<br>(0.951–0.981)</td><td>0.976<br>(0.961–0.987)</td><td>0.949<br>(0.928–0.965)</td><td>0.933<br>(0.909–0.952)</td></tr><tr><td>PPV<br>(CI 95%)</td><td>99.1%<br>(96.6–99.8)</td><td>86.1%<br>(77.5–91.7)</td><td>100%<br>(71.6–99.2)</td><td>94.7%<br>(83.9–91.4)</td><td>88.2%<br>(83.9–91.4)</td></tr><tr><td>NPV<br>(CI 95%)</td><td>93.0%<br>(90.1–95.1)</td><td>99.4%<br>(98.2–99.8)</td><td>99.8%<br>(98.9–99.9)</td><td>99.7%<br>(98.7–99.9)</td><td>96.3%<br>(93.9–97.8)</td></tr><tr><td>Sensitivity (CI 95%)</td><td>95.8%<br>(92.1–98.0)</td><td>95.8%<br>(88.1–99.1)</td><td>95.2%<br>(76.2–99.9)</td><td>90.0%<br>(68.3–98.8)</td><td>94.0%<br>(90.1–96.8)</td></tr><tr><td rowspan="6">P. falciparum symptomatic patients</td><td>Specificity (CI 95%)</td><td>99.4%<br>(97.8–99.9)</td><td>97.6%<br>(95.7–98.8)</td><td>100%<br>(99.3–100)</td><td>99.8%<br>(98.9–99.9)</td><td>96.0%<br>(93.3–97.9)</td></tr><tr><td>AUC<br>(CI 95%)</td><td>0.976<br>(0.958–0.987)</td><td>0.967<br>(0.949–0.981)</td><td>0.976<br>(0.960–0.987)</td><td>0.949<br>(0.927–0.966)</td><td>0.950<br>(0.929–0.967)</td></tr><tr><td>PPV<br>(CI 95%)</td><td>99.0%<br>(96.2–99.8)</td><td>86.1%<br>(77.5–91.7)</td><td>100%<br>(71.6–99.2)</td><td>94.7%<br>(90.3–96.4)</td><td>94.1%<br>(90.3–96.4)</td></tr><tr><td>NPV<br>(CI 95%)</td><td>97.3%<br>(95.1–98.6)</td><td>99.3%<br>(98.1–99.8)</td><td>99.8%<br>(98.7–99.9)</td><td>99.6%<br>(98.5–99.9)</td><td>96.0%<br>(93.4–97.6)</td></tr><tr><td>Sensitivity (CI 95%)</td><td>62.8%<br>(46.7–77.0)</td><td></td><td></td><td></td><td>100%<br>(29.2–100)</td></tr><tr><td>Specificity (CI 95%)</td><td>–</td><td></td><td></td><td></td><td>65.1%<br>(49.1–79.0)</td></tr><tr><td rowspan="6">P. falciparum asymptomatic individuals</td><td>AUC<br>(CI 95%)</td><td>0.314<br>(0.185–0.468)</td><td></td><td></td><td></td><td>0.826<br>(0.699–0.915)</td></tr><tr><td>PPV<br>(CI 95%)</td><td>90.0%<br>(87.7–91.9)</td><td></td><td></td><td></td><td>16.7%<br>(11.7–23.1)</td></tr><tr><td>NPV<br>(CI 95%)</td><td>–</td><td></td><td></td><td></td><td>100%<br>(93.4–97.6)</td></tr><tr><td>Sensitivity (CI 95%)</td><td>62.8%<br>(46.7–77.0)</td><td></td><td></td><td></td><td>100%<br>(29.2–100)</td></tr><tr><td>Specificity (CI 95%)</td><td>–</td><td></td><td></td><td></td><td>65.1%<br>(49.1–79.0)</td></tr><tr><td>AUC<br>(CI 95%)</td><td>0.314<br>(0.185–0.468)</td><td></td><td></td><td></td><td>0.826<br>(0.699–0.915)</td></tr></tbody></table>

AUC areas under the ROC curves, PPV positive predictive value, NPV negative predictive value, CI 95% confidence intervals 95%

The sensitivity and specificity of the K13 bMx prototype assay were ≥ 90%, except for samples collected from asymptomatic individuals. The sensitivity of the K13 bMx prototype to detect the C580Y mutant (the predominant allele in Southeast Asia) was significantly lower in assays using DBS collected from asymptomatic individuals (62.81%) compared to those performed with DBS from symptomatic patients (95.8%). This decrease in sensitivity was mainly due to the low amount of parasite DNA in DBS collected from asymptomatic individuals who were likely carrying parasites at densities lower than the LoD of the assay (ranging from 6.25 to 50 parasites/µL for the 5 mutants in the panel).

In field conditions, at IP Cambodia, the K13 ready-to-use bMx prototype assay was considered by the end-users as a user-friendly assay to perform. The interpretation of the results was facilitated by the design of an Excel file template. The end-users appreciated the shorter time to result than the K13 reference assay. The

major advantages of the K13 bMx prototype over the K13 reference assay raised by the end-users were less requirement of budget planning (budget can be easily estimated per test) and fewer logistical issues linked to reagents purchase and the basic laboratory equipment needed to perform this assay. The findings seem to qualify the K13 bMx prototype assay as a screening tool usable in malaria endemic countries recognized to be at risk of spreading of ART-R parasites from Southeast Asia in their central laboratory dedicated to the monitoring of the anti-malarial drug resistance. However, rare occurrences of ART-R parasites that lack *k13* mutations with an increased RSA<sup>0–3h</sup> survival rates or delayed parasite clearance times have been described [19, 33, 34] and suggest loci other than K13 may modulate ART-R. It seems, therefore, important to continue phenotypic assessment of ART-R, in addition to monitoring emergence of validated K13 mutants.
<figure><img src="image_5.png" /><figcaption>**Fig. 3** Area under the ROC curve values obtained from the K13 bMx prototype assay according to the *k13* mutant detected and the origin of the samples collected in Cambodia, Myanmar and Africa</figcaption></figure>

However, the K13 bMx prototype assay still has some limitations. First, the K13 bMx prototype assay has been developed to detect only 5 *k13* mutants. Nevertheless, the ARMS strategy and the flexible format allows us to update the K13 bMx prototype assay according to the epidemiological change (i.e. emergence of a novel and validated *k13* mutant) and to add new reactions detecting new *k13* mutants. Second, the K13 bMx prototype assay is restricted to the detection of ART-R parasites. However, based on the easyMAG® platform and the Argene Solution®, the development of new field-based molecular assays detecting molecular signatures associated to anti-malarial drug resistances is feasible. New prototypes could include molecular markers associated with other anti-malarials used as ACT partner drugs (mefloquine, piperaquine) as well as molecular markers associated with anti-malarial drugs used in prophylaxis (sulfadoxine–pyrimethamine, atovaquone). To date molecular markers associated to amodiaquine or lumefantrine resistance are not well established, but marker panels could be adapted to detect mutations in *Pfcrt* and *Pfmdr1* genes.

Third, additional and multi-site studies are further needed to evaluate the performances of the K13 bMx prototype assay in field conditions, in different epidemiological contexts, especially settings recognized to be at risk of spreading of ART-R parasites from Southeast Asia (e.g. in Africa, India, South America). The ARTEMIS (for Artemisinin ResisTance fiEld MultIsite Study) project is planned to be conducted in 2018 to extend this first evaluation and to, generate additional evidences to consider the K13 bMx prototype assay a reference epidemiological field-based molecular tool for the global surveillance of ART-R.

## Additional files

**Additional file 1.** K13 bMx prototype assay 5-step workflow (24 DBS in less than 4 h).

**Additional file 2.** New 3 MM filter paper format design to facilitate the semi-automated extraction protocol and avoid inter-samples DNA contamination.

**Additional file 3.** Raw data showing the results between DBS tested at IP Cambodia using an in-house protocol (QiaAmp DNA blood mini kit Qiagen and two-step real-time PCR assays targeting the *Plasmodium cytochrome b* gene) and at bioMérieux using the easyMAG® system (Generic 2.0.1 protocol) and the real-time PCR targeting the *cytochrome c oxidase subunit 1* gene (IC) according to the origin of the DBS.

**Additional file 4.** Principe of the Amplification-refractory mutation system (ARMS) strategy, based on allele-specific primers. In red is the mismatch position at the 3' end of the primer. In blue an additional mismatch at the previous nucleotide.

**Additional file 5.** Overall clinical performances of the K13 reference assay obtained from 642 DBS collected from *P. falciparum* symptomatic patients in Cambodia, Myanmar and Africa and *P. falciparum* asymptomatic individuals in Cambodia. Cells coloured in green present concordant data obtained between the K13 bMx prototype assay and the K13 reference assay; Cells coloured in pale red present discordant data obtained between the K13 bMx prototype assay and the K13 reference assay.

**Additional file 6.** Clinical performances of the K13 reference assay obtained from 565 DBS collected from *P. falciparum* symptomatic patients in Cambodia, Myanmar and Africa. Cells coloured in green present concordant data obtained between the K13 bMx prototype assay and the K13 reference assay; Cells coloured in pale red present discordant data obtained between the K13 bMx prototype assay and the K13 reference assay.

**Additional file 7.** Clinical performances of the K13 reference assay obtained from 65 DBS collected from *P. falciparum* asymptomatic individuals in Cambodia. Cells coloured in green present concordant data obtained between the K13 bMx prototype assay and the K13 reference assay; Cells coloured in pale red present discordant data obtained between the K13 bMx prototype assay and the K13 reference assay.
## Authors' contributions

Sample collections in Cambodia, Myanmar and Africa were coordinated by IP Cambodia (SK), Mahidol-Oxford Tropical Medicine Research Unit and the Shoklo Malaria Research Unit in Thailand (EA, MD, FN), Department of Preventive & Social Medicine in Myanmar (KMT) and National Malaria Control Programs in Chad (IMS). The K13 bMx prototype assay was developed and performed at bioMérieux by AI and LVG. NK, RE, CK, and MK conducted the K13 reference and the bMx prototype assays at IP Cambodia. The study was designed by LVG, EL, FA, AP, KK and DM. Data entry and quality control were performed at IP Cambodia by NK and at bioMérieux by AI. Statistical analyses were performed by SB, LVG and DM. The first draft of the manuscript was written by DM. All authors read and approved the final manuscript.

## Author details

¹ Medical Diagnostic Discovery Department (MD3), bioMérieux, Grenoble and Marcy l'Etoile, France. ² Malaria Molecular Epidemiology Unit, Institute Pasteur in Cambodia, Phnom Penh, Cambodia. ³ Genetics and Genomics of Insect Vectors Unit, Institute Pasteur, Paris, France. ⁴ Malaria Genetics and Resistance Group, Biology of Host-Parasite Interactions Unit, Institute Pasteur, Paris, France. ⁵ Mahidol-Oxford Tropical Medicine Research Unit, Bangkok, Thailand. ⁶ Department of Preventive & Social Medicine, Defence Services Medical Academy, Yangon, Myanmar. ⁷ Myanmar Oxford Clinical Research Unit, Yangon, Myanmar. ⁸ Worldwide Antimalarial Resistance Network, Asia Regional Centre, Bangkok, Thailand. ⁹ Centre for Tropical Medicine and Global Health, Nuffield Department of Medicine Research Building, University of Oxford, Old Road Campus, Oxford, UK. ¹⁰ Shoklo Malaria Research Unit, Mahidol-Oxford Tropical Medicine Research Unit, Faculty of Tropical Medicine, Mahidol University, Mae Sot, Thailand. ¹¹ Programme National de Lutte Contre le Paludisme au Tchad, Ndjamina, Chad. ¹² Institute Cochin Inserm U1016, Université Paris-Descartes, Sorbonne Paris Cité, Paris, France. ¹³ Laboratoire de Parasitologie-Mycologie, Hôpital Cochin, Paris, France.

## Acknowledgements

We are grateful to all patients and field staff without whom this work could not have been done We would like to thank colleagues from the Direction des Applications de la Recherche et des Relations Industrielles (DARRI) at IP Paris (Chloé Pellicier, Karine Mignon-Godefroy and Isabelle Buckle) for their support. We would like to thank Virginie Moucadel (bioMérieux) for her support at the beginning of the project. SMRU is part of the Mahidol Oxford University Research Unit, supported by the Wellcome Trust of Great Britain,

## Competing interests

LVG, AI, SB, AP and KK are bioMérieux employees. The authors declare that they have no competing interests.

## Availability of data and materials

The datasets supporting the conclusions of this paper are included within the paper. Raw data may be obtained from the corresponding author on request.

## Consent for publication

Not applicable.

## Ethics approval and consent to participate

The study protocol was reviewed and approved by the Ethics Committee of the Cambodian Ministry of Health (No 397 NECHR, November 6, 2015) and the Clinical Research Committee of Institut Pasteur, Paris, France (No. 2015-082, October 8, 2015).

DBS from Myanmar and Africa were provided by Pr. Francois Nosten (Shoklo Malaria Research Unit, Mae Sot, Thailand), Dr. Mehul Dhorda (on the behalf of the Tracking Resistance to Artemisinin Collaboration and Mahidol-Oxford Tropical Research Unit, Bangkok, Thailand) and Dr. Issa Mahamat Souleymane (National Malaria Control Programme, Ndjamena, Chad).

## Funding

This study was supported by bioMérieux and Institute Pasteur Paris. The funders had a role in study design but no role in data collection and interpretation, or the decision to submit the work for publication.

## Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Received: 22 January 2018 Accepted: 18 April 2018

Published online: 24 April 2018

## References

1. WHO. World Malaria Report. Geneva: World Health Organization; 2017.

2. Bhatt S, Weiss DJ, Cameron E, Bisanzio D, Mappin B, Dalrymple U, et al. The effect of malaria control on *Plasmodium falciparum* in Africa between 2000 and 2015. Nature. 2015;526:207–11.

3. Nosten F, White NJ. Artemisinin-based combination treatment of falciparum malaria. Am J Trop Med Hyg. 2007;77:181–92.

4. Dondorp AM, Nosten F, Yi P, Das D, Phyó AP, Tarning J, et al. Artemisinin resistance in *Plasmodium falciparum* malaria. N Engl J Med. 2009;361:455–67.

5. Noedl H, Se Y, Schaecher K, Smith BL, Socheat D, Fukuda MM. Evidence of artemisinin-resistant malaria in western Cambodia. N Engl J Med. 2008;359:2619–20.

6. Saunders DL, Vanachayangkul P, Lon C, US Army Military Malaria Research Program, National Center for Parasitology Entomology and Malaria Control, Royal Cambodian Armed Forces. Dihydroartemisinin–piperaquine failure in Cambodia. N Engl J Med. 2014;371:484–5.

7. Witkowski B, Duru V, Khim N, Ross LS, Saintpierre B, Beghain J, et al. A surrogate marker of piperaquine-resistant *Plasmodium falciparum* malaria: a phenotype–genotype association study. Lancet Infect Dis. 2017;17:174–83.

8. Phuc BQ, Rasmussen C, Duong TT, Dong LT, Loi MA, Menard D, et al. Treatment failure of dihydroartemisinin/piperaquine for *Plasmodium falciparum* malaria, Vietnam. Emerg Infect Dis. 2017;23:715–7.

9. Amato R, Lim P, Miotto O, Amaratunga C, Dek D, Pearson RD, et al. Genetic markers associated with dihydroartemisinin–piperaquine failure in *Plasmodium falciparum* malaria in Cambodia: a genotype–phenotype association study. Lancet Infect Dis. 2017;17:164–73.

10. Leang R, Taylor WR, Bouth DM, Song L, Tarning J, Char MC, et al. Evidence of *Plasmodium falciparum* malaria multidrug resistance to artemisinin and piperaquine in Western Cambodia: dihydroartemisinin–piperaquine open-label multicenter clinical assessment. Antimicrob Agents Chemother. 2015;59:4719–26.

11. Duru V, Khim N, Leang R, Kim S, Domergue A, Kloeung N, et al. *Plasmodium falciparum* dihydroartemisinin-piperaquine failures in Cambodia are associated with mutant K13 parasites presenting high survival rates in novel piperaquine in vitro assays: retrospective and prospective investigations. BMC Med. 2015;13:305.

12. Imwong M, Suwannasin K, Kunasol C, Sutawong K, Mayxay M, Rekol H, et al. The spread of artemisinin-resistant *Plasmodium falciparum* in the Greater Mekong subregion: a molecular epidemiology observational study. Lancet Infect Dis. 2017;17:491–7.

13. Straimer J, Gnadig NF, Stokes BH, Ehrenberger M, Crane AA, Fidock DA. *Plasmodium falciparum* K13 mutations differentially impact ozonide susceptibility and parasite fitness *in vitro*. MBio. 2017;8:e00172–817.

14. Ghorbal M, Gorman M, Macpherson CR, Martins RM, Scherf A, Lopez-Rubio JJ. Genome editing in the human malaria parasite *Plasmodium falciparum* using the CRISPR–Cas9 system. Nat Biotechnol. 2014;32:819–21.

15. Ariey F, Witkowski B, Amaratunga C, Beghain J, Langlois AC, Khim N, et al. A molecular marker of artemisinin-resistant *Plasmodium falciparum* malaria. Nature. 2014;505:50–5.

16. WHO. Status report on artemisinin and ACT resistance. Geneva: World Health Organization; 2017.

17. Menard D, Khim N, Beghain J, Adegnika AA, Shafiul-Alam M, Amodu O, et al. A worldwide map of *Plasmodium falciparum* K13-propeller polymorphisms. N Engl J Med. 2016;374:2453–64.

18. Witkowski B, Amaratunga C, Khim N, Sreng S, Chim P, Kim S, et al. Novel phenotypic assays for the detection of artemisinin-resistant *Plasmodium falciparum* malaria in Cambodia: in vitro and ex vivo drug-response studies. Lancet Infect Dis. 2013;13:1043–9.
19. Boulle M, Witkowski B, Duru V, Sriprawat K, Nair SK, McDew-White M. et al. Artemisinin-resistant *Plasmodium falciparum* K13 mutant alleles, Thailand–Myanmar border. Emerg Infect Dis. 2016;22:1503–5.

20. Wang Z, Wang Y, Cabrera M, Zhang Y, Gupta B, Wu Y, et al. Artemisinin resistance at the China–Myanmar border and association with mutations in the K13 propeller gene. Antimicrob Agents Chemother. 2015;59:6952–9.

21. Tun KM, Imwong M, Lwin KM, Win AA, Hlaing TM, Hlaing T, et al. Spread of artemisinin-resistant *Plasmodium falciparum* in Myanmar: a cross-sectional survey of the K13 molecular marker. Lancet Infect Dis. 2015;15:415–21.

22. Nyunt MH, Hlaing T, Oo HW, Tin-Oo LL, Phway HP, Wang B, et al. Molecular assessment of artemisinin resistance markers, polymorphisms in the k13 propeller, and a multidrug-resistance gene in the eastern and western border areas of Myanmar. Clin Infect Dis. 2015;60:1208–15.

23. Menard D, Dondorp A. Antimalarial drug resistance: a threat to malaria elimination. Cold Spring Harb Perspect Med. 2017;7:a025619.

24. Menard D, Ariey F. Towards real-time monitoring of artemisinin resistance. Lancet Infect Dis. 2015;15:367–8.

25. Nsanzabana C, Djalle D, Guerin PJ, Menard D, Gonzalez IJ. Tools for surveillance of anti-malarial drug resistance: an assessment of the current landscape. Malar J. 2018;17:75.

26. Witkowski B, Khim N, Chim P, Kim S, Ke S, Kloeung N, et al. Reduced artemisinin susceptibility of *Plasmodium falciparum* ring stages in western Cambodia. Antimicrob Agents Chemother. 2013;57:914–23.

27. Rossi G, Van den Bergh R, Nguon C, Debackere M, Vernaeve L, Khim N, et al. Adapting reactive case detection strategies for falciparum malaria in a low-transmission area in Cambodia. Clin Infect Dis. 2017;66:296–8.

28. Vantaux A, Samreth R, Piv EP, Khim N, Kim S, Berne L, et al. Contribution to malaria transmission of symptomatic and asymptomatic parasite carriers in Cambodia. J Infect Dis. 2018. https://doi.org/10.1093/infdis/jiy060 (Epub ahead of print)

29. Mercaldo ND, Lau KF, Zhou XH. Confidence intervals for predictive values with an emphasis to case–control studies. Stat Med. 2007;26:2170–83.

30. Canier L, Khim N, Kim S, Eam R, Khean C, Loch K, et al. Malaria PCR detection in Cambodian low-transmission settings: dried blood spots versus venous blood samples. Am J Trop Med Hyg. 2015;92:573–7.

31. Canier L, Khim N, Kim S, Sluydts V, Heng S, Dourng D, et al. An innovative tool for moving malaria PCR detection of parasite reservoir into the field. Malar J. 2013;12:405.

32. Newton CR, Graham A, Heptinstall LE, Powell SJ, Summers C, et al. Analysis of any point mutation in DNA. The amplification refractory mutation system (ARMS). Nucleic Acids Res. 1989;17:2503–16.

33. Miotto O, Amato R, Ashley EA, MacInnis B, Almagro-Garcia J, Amaratunga C, et al. Genetic architecture of artemisinin-resistant *Plasmodium falciparum*. Nat Genet. 2015;47:226–34.

34. Mukherjee A, Bopp S, Magistrado P, Wong W, Daniels R, Demas A, et al. Artemisinin resistance without pfkelch13 mutations in *Plasmodium falciparum* isolates from Cambodia. Malar J. 2017;16:195.

Ready to submit your research? Choose BMC and benefit from:

* fast, convenient online submission

* thorough peer review by experienced researchers in your field

* rapid publication on acceptance

* support for research data, including large and complex data types

* gold Open Access which fosters wider collaboration and increased citations

* maximum visibility for your research: over 100M website views per year

At BMC, research is always in progress.

Learn more biomedcentral.com/submissions

<figure><img src="image_6.png" /></figure>