# Effect of mass dihydroartemisinin-piperaquine administration in southern Mozambique on the carriage of molecular markers of antimalarial resistance

Himanshu Gupta, Beatriz Galatas, Arlindo Chidimatembue, Silvie Huijben, Pau Cisteró, Gloria Matambisso, Lidia Nhamussua, Wilson Simone, Quique Bassat, Didier Menard, et al.

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

## To cite this version:

Himanshu Gupta, Beatriz Galatas, Arlindo Chidimatembue, Silvie Huijben, Pau Cisteró, et al. Effect of mass dihydroartemisinin-piperaquine administration in southern Mozambique on the carriage of molecular markers of antimalarial resistance. PLoS ONE, 2020, 15 (10), pp.e0240174. 10.1371/journal.pone.0240174. pasteur-03274274

HAL Id: pasteur-03274274

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

Submitted on 29 Jun 2021

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published 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_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH ARTICLE

# Effect of mass dihydroartemisinin-piperaquine administration in southern Mozambique on the carriage of molecular markers of antimalarial resistance

Himanshu Gupta¹,²ᵃ, Beatriz Galatas¹,², Arlindo Chidimatembue², Silvie Huijben¹,²ᵇ, Pau Cisteró¹, Gloria Matambisso², Lidia Nhamussua¹,², Wilson Simone², Quique Bassat¹,²,³,⁴,⁵, Didier Ménard⁶,⁷,⁸, Pascal Ringwald⁹, N. Regina Rabinovich¹,¹⁰, Pedro L. Alonso¹,²ᶜ, Francisco Saúte², Pedro Aide²,¹¹, Alfredo Mayor¹,²,³*

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

1 ISGlobal, Barcelona Institute for Global Health, Hospital Clínic—Universitat de Barcelona, Barcelona, Spain, 2 Centro de Investigação em Saúde de Manhiça, Manhica, Mozambique, 3 ICREA, Pg. Lluís Companys, Barcelona, Spain, 4 Spanish Consortium for Research in Epidemiology and Public Health (CIBERESP), Madrid, Spain, 5 Pediatric Infectious Diseases Unit, Pediatrics Department, Hospital Sant Joan de Déu (University of Barcelona), Barcelona, Spain, 6 Institut Pasteur, Paris, France, 7 INSERM U1201, Paris, France, 8 CNRS ERL9195, Paris, France, 9 World Health Organization (WHO), Global Malaria Programme, Geneva, Switzerland, 10 Harvard T.H. Chan School of Public Health, Boston, Massachusetts, United States of America, 11 National Institute of Health, Ministry of Health, Manhica, Mozambique

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

OPEN ACCESS

**Citation:** Gupta H, Galatas B, Chidimatembue A, Huijben S, Cisteró P, Matambisso G, et al. (2020) Effect of mass dihydroartemisinin-piperaquine administration in southern Mozambique on the carriage of molecular markers of antimalarial resistance. PLoS ONE 15(10): e0240174. https://doi.org/10.1371/journal.pone.0240174

**Editor:** Luzia Helena Carvalho, Instituto Rene Rachou, BRAZIL

**Received:** June 27, 2020

**Accepted:** September 21, 2020

**Published:** October 19, 2020

**Copyright:** © 2020 Gupta et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All relevant aggregated data are within the manuscript and its Supporting Information files. The individual level-data are available upon request from CISM's institutional ethics committee (ten.acihnam@oicaca.ohnizos), Hospital Clínic of Barcelona's Ethics Committee (tac.cinilc@CIEC), and the Mozambican Ministry of Health National Bioethics Committee (moc.liamg@hcablawhcslfj) for researchers who meet the criteria for access to confidential data.

## Abstract

### Background

Mass drug administration (MDA) can rapidly reduce the burden of *Plasmodium falciparum* (*Pf*). However, concerns remain about its contribution to select for antimalarial drug resistance.

### Methods

We used Sanger sequencing and real-time PCR to determine the proportion of molecular markers associated with antimalarial resistance (*k13*, *pfpm2*, *pfmdr1* and *pfcr1*) in *Pf* isolates collected before (n = 99) and after (n = 112) the implementation of two monthly MDA rounds with dihydroartemisinin-piperaquine (DHAp) for two consecutive years in Magude district of Southern Mozambique.

### Results

None of the *k13* polymorphisms associated with artemisinin resistance were observed in the *Pf* isolates analyzed. The proportion of *Pf* isolates with multiple copies of *pfpm2*, an amplification associated with piperaquine resistance, was similar in pre- (4.9%) and post-MDA

<footnote>ᵃa Current address: Department of Infection Biology, Faculty of Infectious and Tropical Diseases, London School of Hygiene and Tropical Medicine, London, United Kingdom</footnote>

<footnote>ᵇb Current address: Center for Evolution and Medicine, School of Life Sciences, Arizona State University, Tempe, AZ, United States of America</footnote>

<footnote>ᶜc Current address: Global Malaria Program, World Health Organization, Geneva, Switzerland</footnote>

<footnote>* alfredo.mayor@isglobal.org</footnote>
**Funding:** We would like to thank La Caixa and Bill and Melinda Gates Foundations for providing the funds for this study (grant number INV-008483 / OPP1115265). A.M. is supported by the Departament d'Universitats i Recerca de la Generalitat de Catalunya, Agència de Gestió d'Ajuts Universitaris i de Recerca (2017SGR664). HG was supported (January 2017-January 2019) by the Science and Engineering Research Board (SERB), Department of Science & Technology, Government of India (Overseas Postdoctoral Fellowship, SB/OS/ PDF043/201516). We acknowledge support from the Spanish Ministry of Science and Innovation through the "Centro de Excelencia Severo Ochoa 2019-2023" Program (CEX2018-000806-S), and support from the Generalitat de Catalunya through the CERCA Program. The Centro de Investigação em Saúde de Manhica (CISM) is supported by the Government of Mozambique and the Spanish Agency for International Development Cooperation (AECID). This research is part of ISGlobal's Program on the Molecular Mechanisms of Malaria, which is partially supported by the Fundación Ramón Areces. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**Competing interests:** The authors have declared that no competing interests exist.

groups (3.4%; p = 1.000). No statistically significant differences were observed between pre- and post-MDA groups in the proportion of *Pf* isolates neither with mutations in *pfcrt* and *pfmdr1* genes, nor with the carriage of *pfmdr1* multiple copies (p>0.05).

## Conclusions

This study does not show any evidence of increased frequency of molecular makers of antimalarial resistance after MDA with DHAp in southern Mozambique where markers of antimalarial resistance were absent or low at the beginning of the intervention.

## Introduction

The administration of drugs to whole populations irrespective of disease status aims to prevent and reduce morbidity on the one hand and reduce transmission on the other, altogether improving global health [1]. This strategy, known as mass drug administration (MDA), is recommended by the World Health Organization (WHO) to control or eliminate several neglected tropical pathogens, including bacteria and helminths. Recent studies suggest that MDA, when used as part of a comprehensive and well-organized elimination programme, can be a useful additional tool to accelerate the path towards malaria elimination [2–6].

The selection and subsequent spread of drug resistance is a major concern when administering any antimicrobial agent on a mass scale, especially if the pathogen is being targeted with only a single drug [1]. However, evidence for the establishment of drug resistance at large scale driven by MDA is limited, in part due to the scarcity of programmes which have monitored changes in drug efficacy or potential drug resistance [1]. Few reports have shown emergence of azithromycin-resistant *Treponema pallidum* [7] and temporary increases in carriage of macrolide-resistance following azithromycin MDAs [8,9], although no evidence of drug resistance has been documented after long-term MDA in other situations [1,10]. With regards to anti-malarial MDAs, circumstantial evidence has linked indirect MDA using medicated salts to the emergence of chloroquine (CQ) resistance in the 1980s [5].

Concerns remain as to whether current MDA strategies based on the use of artemisinin-based combination therapies (ACT) might contribute to the emergence and spread of antimalarial drug resistance [4,5]. MDA leads to opposing forces on the selection for resistance and it is not at all that obvious that MDA would always lead to resistance evolution. On the one hand, MDA drastically decreases total number of malaria parasites, which reduces to probability of *de novo* resistance mutations. However, it would also maximize selective pressure for any resistant mutant that exists or does arise [11]. How both these forces play out in different epidemiological contexts is not yet known. Mathematical models have estimated that the implementation of MDA with atovaquone-proguanil would lead to rapid selection for high-level resistance, even after a single round of MDA [12]. Modeling predictions suggested reduction in the effectiveness of subsequent rounds of treatment, with total loss of efficacy within 4–5 years, although mutations in *cytochrome b* gene were later found to be lethal in the mosquito host [13] and thus to impede transmission of resistant parasites. The same models suggested a lower risk in the selection of drug resistance by ACTs such as dihydroartemisinin-piperaquine (DHAp), due to the weaker resistance phenotype resulting from artemisinin resistant mutations known to date [14–16].

Several polymorphisms in the *kelch13* (*k13*) propeller gene of *P. falciparum* have been associated with artemisinin resistance [14–17]. The reduced copy number of *P. falciparum* multidrug resistance 1 (*pfmdr1*) gene has also been linked to increased sensitivity to artemisinin of
trophozoite stages [18], on the contrary, multiple *pfmdr1* copies associated with declining effi-
cacy of mefloquine–artesunate in the Thai–Myanmar border [19]. Recently, *plasmepsin-2*
(*pfpm2*) copy number and polymorphisms in *P. falciparum* chloroquine resistance transporter
(*pfcrt*) gene (H97Y, C101F, F145I, M343L or G353V) have been associated with decreased
piperaquine sensitivity and high DHAp treatment failure rates in settings where artemisinin
resistance is common [20–25]. In addition, evidence suggests that the presence of 86N and
184F alleles in *pfmdr1* gene may reduce the susceptibility to piperaquine in parasites expressing
the CVIET haplotype in the *pfcrt* gene [26].

A before-after study was conducted in southern Mozambique to evaluate the impact of a
package of interventions with the aim to interrupt *Plasmodium falciparum* (*Pf*) malaria trans-
mission [27]. Two rounds of MDA with DHAp per year over 2015 and 2017, together with
annual indoor residual spraying (IRS), programmatically distributed long-lasting insecticide
treated nets (LLINs) and standard case management, lead to a 71.3% reduction of all-age para-
site prevalence by rapid diagnostic test (from 9.1% to 2.6%), and a 61.5% reduction in case
incidence at health facility level (from 195 to 75 cases per 1000). Here we aimed to determine
the effect of MDA with DHAp on the selection of dihydroartemisinin and piperaquine resis-
tance at the population level. To achieve this, we assessed molecular markers of resistance
(*k13*, *pfmdr1*, *pfcrt* and *pfpm2*) among the circulating parasite population before and after the
implementation of MDA with DHAp in Magude district in southern Mozambique.

**Material and methods**

Study site

Two MDA rounds separated by a period of 4–6 weeks were deployed at the beginning of the rainy seasons of 2015–16 (November and January-February) and 2016–17 (December and January-February) in Magude District (Southern Mozambique). The MDAs were targeted to the entire population of Magude (48,488 residents, 10,965 households in 2015) and reached 72.3%, 58%, 66.6% and 64.8% of the population in each of the four rounds, respectively [27].

**Sample collection and study design**

Blood samples were collected on Whatman 3mm filter papers from a random subsample of
individuals during the first MDA round in November 2015 prior to drug administration
(n = 6858) as well as from 3752 of the 3865 individuals participating in an age-stratified cross-
sectional survey conducted during May 2017 (after the four rounds of MDA) among all-age
individuals randomly selected from the Magude population census. qPCR was conducted on a
random selection of 1271 samples collected in November 2015 and in all samples collected
during May 2017.

DNA extraction and *P. falciparum* detection by real time quantitative PCR
(qPCR)

DNA was extracted from a half cut of the dried blood spot (corresponding to approximately
25µL blood) by using a QIAamp DNA Mini kit (Qiagen), as per the manufacturer's instruc-
tions, with a final elution in 100µL of elution buffer. The ABI PRISM 7500 HT Real-Time Sys-
tem (Applied Biosystems) was used to amplify purified parasite DNA templates, using a
previously described method that targets the 18S rRNA gene [28,29]. A standard curve was
prepared from a synchronized *in vitro* culture of 3D7 strain containing known numbers of
ring-infected erythrocytes and performed in triplicate for each test with five serially diluted
points. *P. falciparum* parasitemia was quantified in the samples by interpolating the Ct values
against the standard curve. DNAs extracted from post-MDA samples with low parasitemia (less than 5 parasites/µl) were pre-amplified using PicoPLEX™ WGA Kit as per the manufacturer's instructions (Rubicon Genomics). Amplified samples were not used for copy number estimations using qPCR. PCR amplification and Sanger sequencing were performed on 21 randomly selected samples before and after pre-amplification, to confirm that the pre-amplification step was not affecting allele frequencies obtained by Sanger sequencing.

## Analysis of molecular markers of antimalarial resistance

To assess polymorphisms in *k13*, *pfcrt* and *pfmdr1* genes, purified DNA templates were amplified using 2720 Thermal Cycler (Applied Biosystems) following protocols described elsewhere for *k13* (aminoacids 427–709 of 3D7 [30,31]), *pfcrt* (aminoacids 35–120 [30]), and *pfmdr1* (two separate fragments covering aminoacids 45–209 and 984–1277 [30]) and sent to Genewiz for bi-directional sequencing. Six positive controls with known *k13* alleles and four parasite lines (3D7, 7G8, Dd2 and V1/S) with known *pfcrt* and *pfmdr1* alleles [30], as well as a negative control (water instead of template DNA), were also included. The variations in the test sequences of *k13*, *pfcrt* and *pfmdr1* were identified by sequence alignment against PF3D7_13 43700, PF3D7_0709000 and PF3D7_0523000 reference sequence of 3D7 respectively, retrieved from PlasmoDB. Polymorphisms were confirmed in both sequences obtained using forward and reverse primers. Isolates with mixed alleles were considered as mutant-type for the purposes of polymorphism frequency estimation.

## pfpm2 and pfmdr1 copy number

A qPCR was used to assess variations in the copy number of *pfpm2* and *pfmdr1* genes as described elsewhere [30]. For each run, the *pfpm2* and *pfmdr1* copy numbers of each sample were measured in triplicate and the *pfβ-tubulin* gene was used as an endogenous control. The PCR efficiencies of the *pfpm2*, *pfmdr1* and *pfβ-tubulin* genes were measured using ten-fold dilutions of 3D7 DNA. The specificity of three primer pairs against human gDNA was also determined. Along with no template control, positive controls with the known copies (3–4) of *pfpm2* [30] and *pfmdr1* (Dd2 parasite line) were also included. All samples with $C_t > 35$ for *pfpm2*, *pfmdr1* and *pfβ-tubulin* were not considered for the copy number analysis. The copy number of *pfmdr1* and *pfpm2* genes was estimated as described previously [23]. Samples with estimated copy number values above 1.5 were defined as having multiple copies and confirmed in independent qPCRs [30].

## Data analysis

Sex and residence area of pre- and post-MDA study participants were compared by Fisher's exact test, while age and log-transformed qPCR parasite densities were compared using Student t test. The proportion of mutant alleles for each specific gene was calculated based on the frequencies of samples with wild-type and mutant alleles. The percentage of isolates with multiple copy numbers was also determined. Fisher's exact test was used to compare the proportion of *Pf* isolates with resistant genotypes, as well as with multiple gene copy number, before and after MDA. The statistical significance was defined as a p-value<0.05.

## Ethics

Study ethical approval was obtained from the National Mozambican Ethical Review Committee (Mozambique), pharmaceutical department of the MoH (IRB), and Hospital Clínic
(Barcelona, Spain) ethics review committees, and signed written informed consent was obtained from all participants or from guardians/parents in the case of minors.

## Results

### Study participants and samples analyzed

Dried blood spots (DBS) were collected from a random selection of consenting MDA1 participants and from individuals included in an age-stratified community-based cross-sectional survey, conducted in May 2017 (three months after MDA4 round). A further random selection of 1271 (November 2015) and 3752 (May 2017) DBS analyzed by qPCR targeting *Pf* 18S rRNA identified 168 and 139 *Pf* infections, respectively. Among these, samples which successfully amplified *k13* gene (99 and 112 pre- and post-MDA, respectively) were further selected for the analysis of molecular markers of resistance.

DNAs extracted from samples with low parasitemia (less than 5 parasites/µl, n = 29) in post-MDA isolates, were pre-amplified using PicoPLEX™ WGA Kit before PCR amplification of targeted genes. Allele frequencies were similar in the 21 *Pf* isolates that were tested with and without pre-amplification, meaning that pre-amplification was not affecting allele frequencies (S1 Fig). Study participants pre- and post-MDA were similar in area of residence, age, sex and qPCR-determined *Pf* densities (S1 Table).

### *k13* polymorphisms

The sequences of *k13* were successfully determined in 99 (100%) and 107 (96%) pre- and post-MDA isolates, respectively. As expected, all positive controls sequencing analysis revealed the existence of wild-type and mutant alleles of *k13* polymorphisms. *k13* polymorphisms reported in Cambodian isolates [15] were absent in the *Pf* isolates analyzed in this study. However, 2 novel synonymous polymorphisms were observed at amino acid positions 477 (0.9% [1/107]) and 548 (1.9% [2/107]), as well as a non-synonymous mutation from “aspartic acid” to “asparagine” at codon 641 (2% [2/99]). The polymorphism G690G previously described in *Pf* field isolates collected in 2015 from Mozambique [30] was also observed in the studied isolates (1.9% [2/107]). No statistically significant difference was noticed when the frequency of polymorphisms was compared between pre- and post-MDA groups (p>0.05; Table 1; Fig 1).

### *pfmdr1* polymorphisms

The sequences of *pfmdr1* were successfully determined in 93 (94%) and 105 (94%) pre- and post-MDA isolates, respectively. All positive controls sequencing analysis revealed the existence of *pfmdr1* wild- and mutant-type alleles. Nine polymorphisms (2 [22%] non-synonymous) in pre-MDA and thirteen polymorphisms (8 [61.5%] non-synonymous) in post-MDA groups were identified in the *pfmdr1* gene (Table 2). Y184F was the most common

<table><caption>Table 1. Proportion of *P. falciparum* isolates with *k13* gene polymorphisms in pre- and post-MDA groups.</caption><thead><tr><th rowspan="2"></th><th>Pre-MDA (n = 99)</th><th>Post-MDA (n = 107)</th><th rowspan="2">p*</th></tr><tr><th>n(%)</th><th>n(%)</th></tr></thead><tbody><tr><td>S477S</td><td>0 (0)</td><td>1 (0.9)</td><td>1.000</td></tr><tr><td>G548G</td><td>0 (0)</td><td>2 (1.9)</td><td>0.498</td></tr><tr><td>D641N</td><td>2 (2)</td><td>0 (0)</td><td>0.230</td></tr><tr><td>G690G</td><td>0 (0)</td><td>2 (1.9)</td><td>0.498</td></tr></tbody></table>

<footnote>* Fisher's exact test.</footnote>

<footnote>https://doi.org/10.1371/journal.pone.0240174.t001</footnote>
<figure><img src="image_5.png" /><figcaption>Fig 1. Distribution of *Plasmodium falciparum* k13, *pfmdr1* and *pfcrt* polymorphism frequencies among pre- and post-MDA isolates.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0240174.g001

<table><caption>Table 2. Proportion of *P. falciparum* isolates with *pfmdr1* gene polymorphisms in pre- and post-MDA groups.</caption><thead><tr><th rowspan="2"></th><th>Pre-MDA (n = 93)</th><th>Post-MDA (n = 105)</th><th rowspan="2">p*</th></tr><tr><th>n(%)</th><th>n(%)</th></tr></thead><tbody><tr><td>T70T</td><td>1 (1.1)</td><td>0 (0)</td><td>0.470</td></tr><tr><td>N86Y</td><td>0 (0)</td><td>1 (1.0)</td><td>1.000</td></tr><tr><td>G102G</td><td>3 (3.2)</td><td>3 (2.9)</td><td>1.000</td></tr><tr><td>E160K</td><td>0 (0)</td><td>4 (3.8)</td><td>1.000</td></tr><tr><td>G182G</td><td>5 (5.4)</td><td>6 (5.7)</td><td>1.000</td></tr><tr><td>Y184F</td><td>57 (61.3)</td><td>53 (50.5)</td><td>0.152</td></tr><tr><td>L1030L</td><td>1 (1.1)</td><td>0 (0)</td><td>0.470</td></tr><tr><td>S1034C</td><td>0 (0)</td><td>3 (2.9)</td><td>0.249</td></tr><tr><td>N1042D</td><td>0 (0)</td><td>2 (1.9)</td><td>0.499</td></tr><tr><td>T1069T</td><td>5 (5.4)</td><td>1 (1.0)</td><td>0.101</td></tr><tr><td>S1137S</td><td>1 (1.1)</td><td>0 (0)</td><td>0.470</td></tr><tr><td>D1179D</td><td>2 (2.2)</td><td>2 (1.9)</td><td>1.000</td></tr><tr><td>Q1195Q</td><td>0 (0)</td><td>2 (1.9)</td><td>0.499</td></tr><tr><td>Y1197N</td><td>1 (1.1)</td><td>2 (1.9)</td><td>1.000</td></tr><tr><td>S1214L</td><td>0 (0)</td><td>2 (1.9)</td><td>0.499</td></tr><tr><td>D1246Y</td><td>0 (0)</td><td>2 (1.9)</td><td>0.499</td></tr></tbody></table>

* Fisher's exact test.

https://doi.org/10.1371/journal.pone.0240174.t002
polymorphism, observed at a frequency of 61.3% (57/93) and 50.5% (53/105) in pre- and post-MDA samples, respectively (p = 0.152). The rest of non-synonymous mutations were at frequencies less than 5%, with no evidence of statistically significant differences between pre- and post-MDA groups (Table 2; Fig 1).

## pfcrt polymorphisms

The sequences of *pfcrt* were successfully determined in 88 (89%) and 104 (93%) pre- and post-MDA isolates, respectively. Wild-type and mutant *pfcrt* polymorphisms present in positive controls were successfully detected. All *Pf* isolates carried wild-type allele at codon 72 (C), 74 (M), 75 (N), 76 (K), 97 (H) and 101 (C) in pre-MDA isolates. *P. falciparum* isolates collected post-MDA presented M74I, N75E and K76T polymorphisms along with newly identified I66T and K120E polymorphisms with allele frequencies <2%. No statistically significant differences were observed in polymorphism frequencies between pre- and post-MDA groups (Table 3; Fig 1).

## P. falciparum gene copy number

Sixty-one (62%) pre-MDA and 59 (53%) post-MDA isolates were successfully analyzed for copy number variation of *pfpm2* and *pfmdr1* genes. All positive control copy numbers for *pfpm2* and *pfmdr1* genes were estimated between 3–4 copies. The lowest to highest range of estimated copy numbers were 0.61 to 2.4 (*pfpm2*) and 0.63 to 2.1 (*pfmdr1*) for pre-MDA isolates, and 0.59 to 1.56 (*pfpm2*) and 0.60 to 1.89 (*pfmdr1*) for post-MDA isolates. Using a copy number threshold of 1.5 to define multiple gene copies isolates [30], 3/61 (4.9%) and 2/59 (3.4%) of *Pf* isolates had multiple copies of *pfpm2* in pre- and post-MDA isolates, respectively (Fig 2). Similarly, 3/61 (4.9%) and 4/59 (5.1%) of *Pf* isolates had multiple copies of *pfmdr1* in pre- and post-MDA isolates, respectively (S2 Table). No statistically significant differences were noticed in the proportion of *Pf* isolates with multiple copies of *pfpm2* (p = 1.000) and *pfmdr1* (p = 1.000) between pre- and post-MDA groups.

## Discussion

Although the use of MDA as a tool to rapidly reduce malaria transmission has become increasingly more popular [4,32–36], the question remains as to whether the distribution of antimalarial drugs at a population level could lead to an increase in drug resistance [37]. This study revealed that there was no emergence of drug resistance after the implementation of four MDA rounds with moderate coverage levels (58–72%). Therefore, there is no evidence that

<table><caption>Table 3. Proportion of *P. falciparum* isolates with *pfcrt* gene polymorphisms in pre- and post-MDA groups.</caption><thead><tr><th rowspan="2"></th><th>Pre-MDA (n = 88)</th><th>Post-MDA (n = 104)</th><th rowspan="2">p*</th></tr><tr><th>n(%)</th><th>n(%)</th></tr></thead><tbody><tr><td>I66T</td><td>0 (0)</td><td>1 (1.0)</td><td>1.000</td></tr><tr><td>M74I</td><td>0 (0)</td><td>2 (1.9)</td><td>0.501</td></tr><tr><td>N75E</td><td>0 (0)</td><td>2 (1.9)</td><td>0.501</td></tr><tr><td>K76T</td><td>0 (0)</td><td>2 (1.9)</td><td>0.501</td></tr><tr><td>H97Y</td><td>0 (0)</td><td>0 (0)</td><td>NA</td></tr><tr><td>C101F</td><td>0 (0)</td><td>0 (0)</td><td>NA</td></tr><tr><td>K120E</td><td>0 (0)</td><td>1 (1.0)</td><td>1.000</td></tr></tbody></table>

<footnote>* Fisher's exact test; NA-not applicable.</footnote>

<footnote>https://doi.org/10.1371/journal.pone.0240174.t003</footnote>
<figure><img src="image_6.png" /><figcaption>Fig 2. Copy number of *pfpm2* and *pfmdr1* genes in *P. falciparum* isolates from pre- and post-MDA groups.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0240174.g002

effectively implemented MDA with an ACT in low- to moderate-endemic settings will lead to drug resistance, and therefore supports the use of DHAp as a drug for future MDAs in areas with no detectable markers of resistance at the start of the intervention.

In line with previous studies conducted in Mozambique [30], polymorphisms in *k13* gene associated with artemisinin resistance were not found in the *Pf* isolates collected in Magude district, while the proportion of multiple *pfpm2* copies was below 5% (4.9% before MDA and 4.3% post-MDA). Other markers of piperaquine resistance in *pfcrt* gene (present in the 35–120 aminoacid region of the protein) were completely absent in these isolates. The absence or low levels of molecular markers of parasite resistance against artemisinin and piperaquine suggests good efficacy of these drugs against *Pf* isolates circulating in the study area. This is aligned with the substantial reduction in malaria incidence achieved with the MDA rounds conducted in Magude, together with annual indoor residual spraying, programmatically distributed long-lasting insecticide treated nets and standard case management [27]. Although comforting at this stage, piperaquine resistance may spread in Mozambique if piperaquine drug pressure increases and artemisinin resistance emerges, subsequently facilitating selection of resistance to ACT partner drugs, as observed in Southeast Asia [23,38,39]. However, no statistically significant differences in the proportion of *k13*, *pfcrt* and *pfmdr1* polymorphisms and multiple copies of *pfpm2* and *pfmdr1* genes were observed between pre- and post-MDA *Pf* isolates.

These results are in line with reports from other settings where artemisinin and piperaquine resistance is higher compared to Mozambique. Data from Myanmar show a stable proportion of *k13* wild-type parasites and no piperaquine resistance during a three-year period in which intense DHAp deployment for targeted MDA substantially decreased malaria incidence [40]. Similarly, no statistical differences were observed in the proportion of *k13* polymorphisms in *Pf* isolates collected before and after 3 monthly MDA rounds in Comoros [32]. In Cambodia, transmission of multidrug-resistant *Pf* parasites was interrupted using MDA with high coverage [36], and no reports of clinical cases were reported for at least 1 year. Similarly, other community-based trials of MDA with DHAp have demonstrated this drug to be efficacious, effective and safe in killing malaria parasites in pre-elimination settings [2, 33, 35, 41, 42], without any suggestion of decreased effectiveness due to expansion of antimalarial resistance. A recent study in South-East Asia has also shown that DHAp mass treatments have not selected resistance further in areas with high frequencies of parasites carrying mutations associated with DHAp resistance [43], thus supporting the use of targeted mass treatment. To our knowledge, recent studies describing MDA with DHAp implemented in Africa have not assessed the
molecular markers of DHAp resistance [44]. In contrast to the lack of studies reporting the molecular markers of DHAp resistance in Africa, seasonal malaria chemoprevention trial has reported the frequency of molecular markers of sulfadoxine-pyrimethamine (SP) plus amodiaquine (AQ) resistance, which did not increase significantly over the study period [45].

The absence of a signature of resistance evolution following MDA with DHAp in southern Mozambique and in other endemic settings could be explained by several factors. First, initial signals of artemisinin resistance characterized by slow parasite clearance [46] before the appearance of *in vivo* resistance [47] may lower the probability of successful transmission of a resistant parasite strains present in a population receiving MDA. This transmission may be further reduced by strong vector control approaches. Second, the administration of drugs with a fairly short half-life, such as derivatives of artemisinin, have lower risk compared to longer acting drugs of selecting resistant parasites, as subtherapeutic concentrations are expected to be available for the parasites for only a limited period of time. Additionally, although MDA will expose many asymptomatic infections to the antimalarial used, the low parasite burdens and the effective host defense mechanisms in these subclinical carriers will reduce the risk of emergence of resistance and the transmission potential of any recurrent infection [11]. Fourth, the homogeneity in the population's drug concentration profiles after MDA, which is given to everyone at the same time, reduces the opportunity for selection of a higher level of resistance as parasites are less likely to encounter high drug concentrations in the next host. Lastly, MDA decreases the incidence of symptomatic malaria in the overall population, which results in less use of the antimalarial, lower parasite biomass exposed to the antimalarials and thus fewer opportunities for the selection of drug resistance.

Mathematical modeling suggests that the spread of drug resistance is strongly dependent on treatment coverage, whereby initial allele frequencies and frequency of treatment play key roles as well [48]. Poor coverage, poor adherence or substantial migration, and therefore reduced effectiveness of MDA, will probably increase the probability of selecting resistance. The risk of selecting resistance may be likely reduced if MDAs are deployed when the parasite biomass is at the lowest levels (i.e., during the dry season in areas of low seasonal transmission) and when there is good adherence to the treatment provided. Therefore, areas with malaria elimination strategies in place that reduces malaria transmission, and hence immunity [49], may be more susceptible to a rapid emergence and spread of resistant parasites [50]. In this situation, MDA could select for resistant mutants that are introduced in the population (such as migration). Further studies would be needed to understand if parasite adaptations and compensatory responses to stress situations driven by steep reductions in malaria transmission, such as increased parasite investment on gametocyte production to maximize transmission [51,52], or reductions in the intensity of between-genotype competition within hosts [53], may have a long-term impact on the emergence and spread of antimalarial resistance.

The results of this study are subject to several limitations. First, there is a possibility the intervention did lead to some evolution of resistance, but this was too rare to be detected with the relatively low number of samples analyzed in this study. Furthermore, the presence of multiple clones in each infection may limit the chances to detect mutant alleles is these constitute a minor fraction in the infection. However, the description of DHAp molecular markers will provide the baseline information to identify the potential expansion of resistant parasites if malaria resurges in the area. Second, this study may have missed the emergence of mutations at densities below the sensitivity limit of the PCRs and Sanger sequencing assays. Third, DNA degradation during sample storage, processing, freezing and thawing could potentially have biased the assays towards the detection of wild type polymorphisms if parasites carrying mutations associated with resistance are more frequently being degraded due to their lower densities. Fourth, our sequencing approach [30] did not include polymorphisms in codons 145, 343
or 353 of *pfcrt*, which have been associated with decreased piperaquine sensitivity and high DHAp treatment failure rates in settings where artemisinin resistance is common [20–25]. Lastly, we cannot discard the possibility of emerging resistance following MDA if conducted on a much larger scale, although this study demonstrates the likelihood of such an event happening is relatively low.

Overall, this study shows that in sub-Saharan setting with absent or low background levels of artemisinin and piperaquine resistance, two monthly MDA rounds with DHAp moderate coverage levels for two consecutive years would not increase the frequency of molecular markers of antimalarial resistance in the general *Pf* parasite population. Based on the findings of the present study and other previous studies, there is no suggestion that a well-implemented MDA with an ACT has led to drug resistance [3,5]. Therefore, DHAp proves to be safe and effective tool in reducing the malaria burden [2,35,36,40,54] and could be useful for future MDAs following WHO recommendations for MDA implementation. Enabling good molecular surveillance systems should be a prerequisite for the use of community-wide distribution of antimalarials aiming the interruption of malaria transmission.

**Supporting information**

S1 Fig. Sequencing results of *P. falciparum* samples before (A) and after (B) pre-amplification.
(PDF)

S1 Table. Main characteristics of study participants.
(PDF)

S2 Table. Copy number of *pfpm2* and *pfmdr1* in *P. falciparum* isolates from pre- and post-MDA groups.
(PDF)

**Acknowledgments**

We would like to express our gratitude to the community of Magude and the district authorities, for allowing this project to take place and for greatly contributing to its successful implementation and evaluation. We would also like to thank all the field workers, the staff of the Hospitals, clinical officers, field supervisors, data manager and those who worked in the lab for their contribution to the collection and organization of samples shipment. We thank the Mozambican National Malaria Control Program and the Ministry of Health for supporting and providing guidance throughout the project. P.R. and P.L.A are staff members of the World Health Organization. Authors alone are responsible for the views expressed in this publication, and these do not necessarily represent the decisions, policy, or views of the World Health Organization.

**Author Contributions**

**Conceptualization:** Himanshu Gupta, Beatriz Galatas, Silvie Huijben, Pascal Ringwald, N. Regina Rabinovich, Alfredo Mayor.

**Data curation:** Himanshu Gupta, Beatriz Galatas, Alfredo Mayor.

**Formal analysis:** Himanshu Gupta, Arlindo Chidimatembue, Pau Cisteró, Didier Ménard, Alfredo Mayor.

**Funding acquisition:** Himanshu Gupta, Pedro L. Alonso, Francisco Saúte, Pedro Aide, Alfredo Mayor.
**Investigation:** Himanshu Gupta, Beatriz Galatas, Gloria Matambisso, Lidia Nhamussua, Wilson Simone, Quique Bassat.

**Methodology:** Himanshu Gupta.

**Project administration:** Beatriz Galatas, Francisco Saúte, Pedro Aide, Alfredo Mayor.

**Resources:** Didier Ménard.

**Supervision:** Alfredo Mayor.

**Writing – original draft:** Himanshu Gupta, Alfredo Mayor.

**Writing – review & editing:** Himanshu Gupta, Beatriz Galatas, Arlindo Chidimatembue, Silvie Huijben, Pau Cisteró, Gloria Matambisso, Lidia Nhamussua, Wilson Simone, Quique Bassat, Didier Ménard, Pascal Ringwald, N. Regina Rabinovich, Pedro L. Alonso, Francisco Saúte, Pedro Aide, Alfredo Mayor.

## References

1. Webster JP, Molyneux DH, Hotez PJ, Fenwick A. The contribution of mass drug administration to global health: past, present and future. Philos Trans R Soc Lond B Biol Sci. 2014; 369(1645):20130434. Epub 2014/05/14. https://doi.org/10.1098/rstb.2013.0434 PMID: 24821920; PubMed Central PMCID: PMC4024227.

2. von Seidlein L, Peto TJ, Landier J, Nguyen TN, Tripura R, Phommasone K, et al. The impact of targeted malaria elimination with mass drug administrations on falciparum malaria in Southeast Asia: A cluster randomised trial. PLoS Med. 2019; 16(2):e1002745. Epub 2019/02/16. https://doi.org/10.1371/journal.pmed.1002745 PMID: 30768615; PubMed Central PMCID: PMC6377128 stipend as a Specialty Consulting Editor for PLOS Medicine and serves on the journal's Editorial Board. NJW is a member of the Editorial Board of PLOS Medicine.

3. Newby G, Hwang J, Koita K, Chen I, Greenwood B, von Seidlein L, et al. Review of mass drug administration for malaria and its operational challenges. Am J Trop Med Hyg. 2015; 93(1):125–34. Epub 2015/05/28. https://doi.org/10.4269/ajtmh.14-0254 PMID: 26013371; PubMed Central PMCID: PMC4497884.

4. Poirot E, Skarbinski J, Sinclair D, Kachur SP, Slutsker L, Hwang J. Mass drug administration for malaria. Cochrane Database Syst Rev. 2013;(12):CD008846. Epub 2013/12/10. https://doi.org/10.1002/14651858.CD008846.pub2 PMID: 24318836; PubMed Central PMCID: PMC4468927.

5. von Seidlein L, Greenwood BM. Mass administrations of antimalarial drugs. Trends Parasitol. 2003; 19(10):452–60. Epub 2003/10/02. https://doi.org/10.1016/j.pt.2003.08.003 PMID: 14519583.

6. Kaneko A, Taleo G, Kalkoa M, Yamar S, Kobayakawa T, Bjorkman A. Malaria eradication on islands. Lancet. 2000; 356(9241):1560–4. Epub 2000/11/15. https://doi.org/10.1016/S0140-6736(00)03127-5 PMID: 11075770.

7. Mitja O, Godornes C, Houinei W, Kapa A, Paru R, Abel H, et al. Re-emergence of yaws after single mass azithromycin treatment followed by targeted treatment: a longitudinal study. Lancet. 2018; 391(10130):1599–607. Epub 2018/02/13. https://doi.org/10.1016/S0140-6736(18)30204-6 PMID: 29428183; PubMed Central PMCID: PMC5920722.

8. O'Brien KS, Emerson P, Hooper PJ, Reingold AL, Dennis EG, Keenan JD, et al. Antimicrobial resistance following mass azithromycin distribution for trachoma: a systematic review. Lancet Infect Dis. 2019; 19(1):e14–e25. Epub 2018/10/08. https://doi.org/10.1016/S1473-3099(18)30444-4 PMID: 30292480.

9. Ho DK, Sawicki C, Grassly N. Antibiotic Resistance in Streptococcus pneumoniae after Azithromycin Distribution for Trachoma. J Trop Med. 2015; 2015:917370. Epub 2015/11/12. https://doi.org/10.1155/2015/917370 PMID: 26557143; PubMed Central PMCID: PMC4628654.

10. West SK, Moncada J, Munoz B, Mkocha H, Storey P, Hardick J, et al. Is there evidence for resistance of ocular Chlamydia trachomatis to azithromycin after mass treatment for trachoma control? J Infect Dis. 2014; 210(1):65–71. Epub 2014/01/22. https://doi.org/10.1093/infdis/jiu046 PMID: 24446528.

11. White NJ. Does antimalarial mass drug administration increase or decrease the risk of resistance? Lancet Infect Dis. 2017; 17(1):e15–e20. Epub 2016/11/15. https://doi.org/10.1016/S1473-3099(16)30269-9 PMID: 27839929.

12. Maude RJ, Nguon C, Dondorp AM, White LJ, White NJ. The diminishing returns of atovaquone-proguanil for elimination of Plasmodium falciparum malaria: modelling mass drug administration and treatment.
Malar J. 2014; 13:380. https://doi.org/10.1186/1475-2875-13-380 PMID: 25249272; PubMed Central PMCID: PMC4192368.

13. Goodman CD, Siregar JE, Mollard V, Vega-Rodriguez J, Syafruddin D, Matsuoka H, et al. Parasites resistant to the antimalarial atovaquone fail to transmit by mosquitoes. Science. 2016; 352(6283):349–53. Epub 2016/04/16. https://doi.org/10.1126/science.aad9279 PMID: 27081071; PubMed Central PMCID: PMC5149070.

14. 2017. WHO. Artemisinin and artemisinin-based combination therapy resistance. Status report. Geneva: World Health Organization.

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(7481):50–5. https://doi.org/10.1038/nature12876 PMID: 24352242; PubMed Central PMCID: PMC5007947.

16. Ashley EA, Dhorda M, Fairhurst RM, Amaratunga C, Lim P, Suon S, et al. Spread of artemisinin resistance in Plasmodium falciparum malaria. N Engl J Med. 2014; 371(5):411–23. https://doi.org/10.1056/NEJMoa1314981 PMID: 25075834; PubMed Central PMCID: PMC4143591.

17. Group WKG-PS. Association of mutations in the Plasmodium falciparum Kelch13 gene (Pf3D7_1343700) with parasite clearance rates after artemisinin-based treatments-a WWARN individual patient data meta-analysis. BMC Med. 2019; 17(1):1. Epub 2019/01/18. https://doi.org/10.1186/s12916-018-1207-3 PMID: 30651111; PubMed Central PMCID: PMC6335805.

18. Sidhu AB, Uhlemann AC, Valderramos SG, Valderramos JC, Krishna S, Fidock DA. Decreasing pfmdr1 copy number in plasmodium falciparum malaria heightens susceptibility to mefloquine, lumefantrine, halofantrine, quinine, and artemisinin. J Infect Dis. 2006; 194(4):528–35. https://doi.org/10.1086/507115 PMID: 16845638; PubMed Central PMCID: PMC2978021.

19. Phyo AP, Ashley EA, Anderson TJC, Bozdech Z, Carrara VI, Sriprawat K, et al. Declining Efficacy of Artemisinin Combination Therapy Against P. Falciparum Malaria on the Thai-Myanmar Border (2003–2013): The Role of Parasite Genetic Factors. Clin Infect Dis. 2016; 63(6):784–91. Epub 2016/06/18. https://doi.org/10.1093/cid/ciw388 PMID: 27313266; PubMed Central PMCID: PMC4996140.

20. Agrawal S, Moser KA, Morton L, Cummings MP, Parihar A, Dwivedi A, et al. Association of a Novel Mutation in the Plasmodium falciparum Chloroquine Resistance Transporter With Decreased Piperaquine Sensitivity. J Infect Dis. 2017; 216(4):468–76. https://doi.org/10.1093/infdis/jix334 PMID: 28931241; PubMed Central PMCID: PMC5853219.

21. 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-phe-notype association study. Lancet Infect Dis. 2017; 17(2):164–73. https://doi.org/10.1016/S1473-3099(16)30409-1 PMID: 27818095; PubMed Central PMCID: PMC5564489.

22. Dhingra SK, Redhi D, Combrinck JM, Yeo T, Okombo J, Henrich PP, et al. A Variant PfCRT Isoform Can Contribute to Plasmodium falciparum Resistance to the First-Line Partner Drug Piperaquine. MBio. 2017; 8(3). https://doi.org/10.1128/mBio.00303-17 PMID: 28487425; PubMed Central PMCID: PMC5424201.

23. 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(2):174–83. https://doi.org/10.1016/S1473-3099(16)30415-7 PMID: 27818097; PubMed Central PMCID: PMC5266792.

24. Ross LS, Dhingra SK, Mok S, Yeo T, Wicht KJ, Kumpornsin K, et al. Emerging Southeast Asian PfCRT mutations confer Plasmodium falciparum resistance to the first-line antimalarial piperaquine. Nature communications. 2018; 9(1):3314. https://doi.org/10.1038/s41467-018-05652-0 PMID: 30115924; PubMed Central PMCID: PMC6095916.

25. Kim J, Tan YZ, Wicht KJ, Erramilli SK, Dhingra SK, Okombo J, et al. Structure and drug resistance of the Plasmodium falciparum transporter PfCRT. Nature. 2019; 576(7786):315–20. Epub 2019/11/30. https://doi.org/10.1038/s41586-019-1795-x PMID: 31776516; PubMed Central PMCID: PMC6911266.

26. Veiga MI, Dhingra SK, Henrich PP, Straimer J, Gnadig N, Uhlemann AC, et al. Globally prevalent PfMDR1 mutations modulate Plasmodium falciparum susceptibility to artemisinin-based combination therapies. Nat Commun. 2016; 7:11553. https://doi.org/10.1038/ncomms11553 PMID: 27189525; PubMed Central PMCID: PMC4873939.

27. Galatas B, Saute F, Marti-Soler H, Guinovart C, Nhamussua L, Simone W, et al. A multiphase program for malaria elimination in southern Mozambique (the Magude project): A before-after study. PLoS Med. 2020; 17(8):e1003227. Epub 2020/08/17. https://doi.org/10.1371/journal.pmed.1003227 PMID: 32797101; PubMed Central PMCID: PMC7428052 the Editorial Board of PLOS Medicine.

28. Mayor A, Serra-Casas E, Bardaji A, Sanz S, Puyol L, Cistero P, et al. Sub-microscopic infections and long-term recrudescence of Plasmodium falciparum in Mozambican pregnant women. Malar J. 2009; 8:9. https://doi.org/10.1186/1475-2875-8-9 PMID: 19134201; PubMed Central PMCID: PMC2633011.
29. Taylor SM, Mayor A, Mombo-Ngoma G, Kenguele HM, Ouedraogo S, Nam NT, et al. A quality control program within a clinical trial Consortium for PCR protocols to detect Plasmodium species. J Clin Microbiol. 2014; 52(6):2144-9. https://doi.org/10.1128/JCM.00565-14 PMID: 24740073; PubMed Central PMCID: PMC4042732.

30. Gupta H, Macete E, Bulo H, Salvador C, Warsame M, Carvalho E, et al. Drug-Resistant Polymorphisms and Copy Numbers in Plasmodium falciparum, Mozambique, 2015. Emerg Infect Dis. 2018; 24(1):40-8. https://doi.org/10.3201/eid2401.170864 PMID: 29260689; PubMed Central PMCID: PMC5749434.

31. 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(25):2453-64. https://doi.org/10.1056/NEJMoa1513137 PMID: 27332904; PubMed Central PMCID: PMC4955562.

32. Deng C, Huang B, Wang Q, Wu W, Zheng S, Zhang H, et al. Large-scale Artemisinin-Piperaquine Mass Drug Administration With or Without Primaquine Dramatically Reduces Malaria in a Highly Endemic Region of Africa. Clinical infectious diseases: an official publication of the Infectious Diseases Society of America. 2018. https://doi.org/10.1093/cid/ciy364 PMID: 29846536.

33. Eisele TP, Bennett A, Silumbe K, Finn TP, Chalwe V, Kamuliwo M, et al. Short-term Impact of Mass Drug Administration With Dihydroartemisinin Plus Piperaquine on Malaria in Southern Province Zambia: A Cluster-Randomized Controlled Trial. J Infect Dis. 2016; 214(12):1831-9. Epub 2016/12/08. https://doi.org/10.1093/infdis/jiw416 PMID: 27923947; PubMed Central PMCID: PMC5142084.

34. Pongvongsa T, Phommasone K, Adhikari B, Henriques G, Chotivanich K, Hanboonkunupakarn B, et al. The dynamic of asymptomatic Plasmodium falciparum infections following mass drug administrations with dihydroartemisinin-piperaquine plus a single low dose of primaquine in Savannakhet Province, Laos. Malar J. 2018; 17(1):405. Epub 2018/11/06. https://doi.org/10.1186/s12936-018-2541-9 PMID: 30390647; PubMed Central PMCID: PMC6215638.

35. Landier J, Kajeechiwa L, Thwin MM, Parker DM, Chaumeau V, Wiladphaingern J, et al. Safety and effectiveness of mass drug administration to accelerate elimination of artemisinin-resistant falciparum malaria: A pilot trial in four villages of Eastern Myanmar. Wellcome Open Res. 2017; 2:81. Epub 2017/10/25. https://doi.org/10.12688/wellcomeopenres.12240.1 PMID: 29062913; PubMed Central PMCID: PMC5635445.

36. Tripura R, Peto TJ, Chea N, Chan D, Mukaka M, Sirithiranont P, et al. A Controlled Trial of Mass Drug Administration to Interrupt Transmission of Multidrug-Resistant Falciparum Malaria in Cambodian Villages. Clin Infect Dis. 2018; 67(6):817-26. Epub 2018/03/10. https://doi.org/10.1093/cid/ciy196 PMID: 29522113; PubMed Central PMCID: PMC6117448.

37. Zuber JA, Takala-Harrison S. Multidrug-resistant malaria and the impact of mass drug administration. Infection and drug resistance. 2018; 11:299-306. https://doi.org/10.2147/IDR.S123887 PMID: 29535546; PubMed Central PMCID: PMC5840189.

38. 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. Antimicrobial agents and chemotherapy. 2015; 59(8):4719-26. https://doi.org/10.1128/AAC.00835-15 PMID: 26014949; PubMed Central PMCID: PMC4505193.

39. Phuc BQ, Rasmussen C, Duong TT, Dong LT, Loi MA, Menard D, et al. Treatment Failure of Dihydroartemisinin/Piperaquine for Plasmodium falciparum Malaria, Vietnam. Emerging infectious diseases. 2017; 23(4):715-7. https://doi.org/10.3201/eid2304.161872 PMID: 28322709; PubMed Central PMCID: PMC5367417.

40. Landier J, Parker DM, Thu AM, Lwin KM, Delmas G, Nosten FH, et al. Effect of generalised access to early diagnosis and treatment and targeted mass drug administration on Plasmodium falciparum malaria in Eastern Myanmar: an observational study of a regional elimination programme. Lancet. 2018; 391(10133):1916-26. Epub 2018/04/29. https://doi.org/10.1016/S0140-6736(18)30792-X PMID: 29703425; PubMed Central PMCID: PMC5946089.

41. Zongo I, Milligan P, Compaore YD, Some AF, Greenwood B, Tarning J, et al. Randomized Noninferiority Trial of Dihydroartemisinin-Piperaquine Compared with Sulfadoxine-Pyrimethamine plus Amodiaquine for Seasonal Malaria Chemoprevention in Burkina Faso. Antimicrob Agents Chemother. 2015; 59(8):4387-96. Epub 2015/04/29. https://doi.org/10.1128/AAC.04923-14 PMID: 25918149; PubMed Central PMCID: PMC4505196.

42. Some AF, Zongo I, Compaore YD, Sakande S, Nosten F, Ouedraogo JB, et al. Selection of drug resistance-mediating Plasmodium falciparum genetic polymorphisms by seasonal malaria chemoprevention in Burkina Faso. Antimicrob Agents Chemother. 2014; 58(7):3660-5. Epub 2014/04/16. https://doi.org/10.1128/AAC.02406-14 PMID: 24733476; PubMed Central PMCID: PMC4068591.

43. Imwong M, Dhorda M, Myo Tun K, Thu AM, Phyo AP, Proux S, et al. Molecular epidemiology of resistance to antimalarial drugs in the Greater Mekong subregion: an observational study. Lancet Infect Dis. 2020. Epub 2020/07/18. https://doi.org/10.1016/S1473-3099(20)30228-0 PMID: 32679084.
44. Eisele TP, Bennett A, Silumbe K, Finn TP, Porter TR, Chalwe V, et al. Impact of Four Rounds of Mass Drug Administration with Dihydroartemisinin-Piperaquine Implemented in Southern Province, Zambia. Am J Trop Med Hyg. 2020 Aug; 103(2_Suppl):7-18. https://doi.org/10.4269/ajtmh.19-0659 PMID: 32618247; PubMed Central PMCID: PMC7416977.

45. Cairns ME, Sagara I, Zongo I, Kuepfer I, Thera I, Nikiema F, et al. Evaluation of seasonal malaria chemoprevention in two areas of intense seasonal malaria transmission: Secondary analysis of a household-randomised, placebo-controlled trial in Hounde District, Burkina Faso and Bougouni District, Mali. PLoS Med. 2020; 17(8):e1003214. Epub 2020/08/22. https://doi.org/10.1371/journal.pmed.1003214 PMID: 32822362.

46. Dondorp AM, Nosten F, Yi P, Das D, Phyo AP, Tarning J, et al. Artemisinin resistance in Plasmodium falciparum malaria. N Engl J Med. 2009; 361(5):455-67. Epub 2009/07/31. https://doi.org/10.1056/NEJMoa0808859 PMID: 19641202; PubMed Central PMCID: PMC3495232.

47. Tun KM, Jeeyapant A, Imwong M, Thein M, Aung SS, Hlaing TM, et al. Parasite clearance rates in Upper Myanmar indicate a distinctive artemisinin resistance phenotype: a therapeutic efficacy study. Malar J. 2016; 15:185. Epub 2016/04/03. https://doi.org/10.1186/s12936-016-1240-7 PMID: 27036739; PubMed Central PMCID: PMC4815199.

48. Schwab AE, Churcher TS, Schwab AJ, Basanez MG, Prichard RK. Population genetics of concurrent selection with albendazole and ivermectin or diethylcarbamazine on the possible spread of albendazole resistance in Wuchereria bancrofti. Parasitology. 2006; 133(Pt 5):589-601. Epub 2006/07/13. https://doi.org/10.1017/S003118200600076X PMID: 16834821.

49. Mayor A, Bardaji A, Macete E, Nhampossa T, Fonseca AM, Gonzalez R, et al. Changing Trends in P. falciparum Burden, Immunity, and Disease in Pregnancy. N Engl J Med. 2015; 373(17):1607-17. Epub 2015/10/22. https://doi.org/10.1056/NEJMoa1406459 PMID: 26488692.

50. Scott N, Ataide R, Wilson DP, Hellard M, Price RN, Simpson JA, et al. Implications of population-level immunity for the emergence of artemisinin-resistant malaria: a mathematical model. Malar J. 2018; 17 (1):279. Epub 2018/08/04. https://doi.org/10.1186/s12936-018-2418-y PMID: 30071877; PubMed Central PMCID: PMC6071336.

51. Gadalla AA, Schneider P, Churcher TS, Nassir E, Abdel-Muhsin AA, Ranford-Cartwright LC, et al. Associations between Season and Gametocyte Dynamics in Chronic Plasmodium falciparum Infections. PLoS One. 2016; 11(11):e0166699. Epub 2016/11/22. https://doi.org/10.1371/journal.pone.0166699 PMID: 27870874; PubMed Central PMCID: PMC5117706.

52. Rono MK, Nyonda MA, Simam JJ, Ngoi JM, Mok S, Kortok MM, et al. Adaptation of Plasmodium falciparum to its transmission environment. Nat Ecol Evol. 2018; 2(2):377-87. Epub 2017/12/20. https://doi.org/10.1038/s41559-017-0419-9 PMID: 29255304.

53. Anderson TJ, Haubold B, Williams JT, Estrada-Franco JG, Richardson L, Mollinedo R, et al. Microsatellite markers reveal a spectrum of population structures in the malaria parasite Plasmodium falciparum. Mol Biol Evol. 2000; 17(10):1467-82. Epub 2000/10/06. https://doi.org/10.1093/oxfordjournals.molbev.a026247 PMID: 11018154.

54. Deng C, Huang B, Wang Q, Wu W, Zheng S, Zhang H, et al. Large-scale Artemisinin-Piperaquine Mass Drug Administration With or Without Primaquine Dramatically Reduces Malaria in a Highly Endemic Region of Africa. Clin Infect Dis. 2018; 67(11):1670-6. Epub 2018/05/31. https://doi.org/10.1093/cid/ciy364 PMID: 29846536; PubMed Central PMCID: PMC6455902.