key: cord-0049935-a6za55dq
authors: Huang, Fang; Shrestha, Biraj; Liu, Hui; Tang, Lin-Hua; Zhou, Shui-Sen; Zhou, Xiao-Nong; Takala-Harrison, Shannon; Ringwald, Pascal; Nyunt, Myaing M.; Plowe, Christopher V.
title: No evidence of amplified Plasmodium falciparum plasmepsin II gene copy number in an area with artemisinin-resistant malaria along the China–Myanmar border
date: 2020-09-14
journal: Malar J
DOI: 10.1186/s12936-020-03410-6
sha: 3eb1448a31a19cbc03b332559b61cece8735f0e2
doc_id: 49935
cord_uid: a6za55dq

BACKGROUND: The emergence and spread of artemisinin resistance in Plasmodium falciparum poses a threat to malaria eradication, including China’s plan to eliminate malaria by 2020. Piperaquine (PPQ) resistance has emerged in Cambodia, compromising an important partner drug that is widely used in China in the form of dihydroartemisinin (DHA)-PPQ. Several mutations in a P. falciparum gene encoding a kelch protein on chromosome 13 (k13) are associated with artemisinin resistance and have arisen spread in the Great Mekong subregion, including the China–Myanmar border. Multiple copies of the plasmepsin II/III (pm2/3) genes, located on chromosome 14, have been shown to be associated with PPQ resistance. METHODS: The therapeutic efficacy of DHA-PPQ for the treatment of uncomplicated P. falciparum was evaluated along the China–Myanmar border from 2010 to 2014. The dry blood spots samples collected in the efficacy study prior DHA-PPQ treatment and from the local hospital by passive detection were used to amplify k13 and pm2. Polymorphisms within k13 were genotyped by capillary sequencing and pm2 copy number was quantified by relative-quantitative real-time polymerase chain reaction. Treatment outcome was evaluated with the World Health Organization protocol. A linear regression model was used to estimate the association between the day 3 positive rate and k13 mutation and the relationship of the pm2 copy number variants and k13 mutations. RESULTS: DHA-PPQ was effective for uncomplicated P. falciparum infection in Yunnan Province with cure rates > 95%. Twelve non synonymous mutations in the k13 domain were observed among the 268 samples with the prevalence of 44.0% and the predominant mutation was F446I with a prevalence of 32.8%. Only one sample was observed with multi-copies of pm2, including parasites with and without k13 mutations. The therapeutic efficacy of DHA-PPQ was > 95% along the China–Myanmar border, consistent with the lack of amplification of pm2. CONCLUSION: DHA-PPQ for uncomplicated P. falciparum infection still showed efficacy in an area with artemisinin-resistant malaria along the China–Myanmar border. There was no evidence to show PPQ resistance by clinical study and molecular markers survey. Continued monitoring of the parasite population using molecular markers will be important to track emergence and spread of resistance in this region.

Malaria remains one of the global major public health problems. According to the latest world malaria report, there were an estimated 228 million cases and 405,000 deaths from malaria globally in 2018, compared with 416,000 estimated deaths in 2017 and 585,000 in 2010 [1] .

Artemisinin-based combination therapy (ACT), which combines a fast-acting, rapidly eliminated artemisinin derivative with another slower-acting partner drug with a longer half-life, have played an indispensable role in reducing global malaria-associated mortality and morbidity [2] . Currently, five artemisinin-based combinations are recommended by the World Health Organization (WHO): artesunate-amodiaquine, artemether-lumefantrine, artesunate-mefloquine, artesunate-sulfadoxinepyrimethamine, and dihydroartemisinin-piperaquine (DHA-PPQ) [3] . ACT has been integral to the recent success of global malaria control and protecting its efficacy for the treatment of malaria has been a global health priority. However, the emergence and spread of artemisinin resistance in Plasmodium falciparum poses a threat to malaria control and eradication goals in the Greater Mekong subregion (GMS), where resistance has emerged independently and spread [4] [5] [6] [7] . The WHO has implemented a strategy to eliminate P. falciparum from the six countries located in the GMS by 2025 to respond the threat of an untreatable multi-drug resistant parasite [8] .

Several mutations in a P. falciparum gene encoding a kelch protein on chromosome 13 (k13) are associated with artemisinin resistance and have arisen multiple times and spread in the GMS, including the China-Myanmar border [9] . Over 200 nonsynonymous k13 mutations have been reported to date, of which nine variants (F446I, N458Y, M476I, Y493H, R539T, I543T, P553L, R561H, and C580Y) have been validated, and over 20 k13 mutations are considered candidates or associated markers [9] .

From 2010, PPQ became a partner drug with DHA to treat falciparum malaria recommended by the WHO [10] . DHA-PPQ, as the first-line drug for uncomplicated P. falciparum treatment, is widely used in China [11] and high efficacy has been reported [12] [13] [14] [15] . Recently, the emergence of DHA-PPQ resistance was observed in Cambodia with high treatment failure rate and recrudescent infections. Resistance then rapidly spread to other countries in Southeast Asia [16] [17] [18] [19] [20] [21] [22] , and subsequent investigations confirmed the presence of high treatment failure and PPQ resistance in Cambodia and Vietnam [16, 21, 23] . The evidence prompted a treatment policy change in Cambodia and Vietnam from DHA-PPQ to other artemisinin-based combinations in areas where DHA-PPQ was failing. Multiple copies of the P. falciparum plasmepsin II (pm2) genes (PF3D7_1408000), located on chromosome 14, encoding a protease involved in haemoglobin degradation, have recently been shown to be associated with PPQ resistance in Cambodia [19, 24] . This situation raises concern about resistance to the partner drugs of ACT.

There had been zero indigenous malaria infections in China from 2017 and China has vowed to interrupt local malaria transmission by 2020 [25] . Yunnan Province located in Southern China, bordering Myanmar, Laos and Vietnam is the key focus of the national malaria elimination program. In this study, we investigated the therapeutic efficacy of DHA-PPQ for treatment of uncomplicated P. falciparum along the China-Myanmar borders from 2010 to 2014. Sequencing of the propeller domain of the k13 gene was performed to identify mutations and pm2 copy number was quantified to determine the copy variants from both parasites with and without k13 mutations.

The clinical studies were one-arm prospective evaluations to observe treatment of uncomplicated malaria using a standard WHO therapeutic efficacy study (TES) protocol. The studies were conducted in four counties (Yingjiang, Tengchong, Menglian and Ruili) in Southwest Yunnan along the China-Myanmar border between 2010 and 2014 ( Fig. 1 ).

Patients aged > 6 months with fever (axillary temperature ≥ 37.5 ℃) or a history of fever in the previous 48 h were screened for inclusion. Inclusion criteria were: mono-infection with P. falciparum; parasitaemia density between 500 and 100,000 asexual parasites/μl; no history of anti-malarial use in the past 14 days and no signs of severe malaria or danger signs. After written informed consent was provided, a detailed medical history, clinical examination and blood smears were performed for each participant. DHA-PPQ (Zhejiang Holley Nanhu Pharmaceutical Co. Ltd, China) was administered at a total adult dose of 2.5 mg/kg dihydroartemisnin and 20 mg/kg PPQ for 3 days. All the anti-malarial drugs were provided by the WHO.

After the first day of treatment (day 0), clinical and laboratory tests, including axillary temperature Keywords: Plasmodium falciparum, Artemisinin resistance, Piperaquine, Plasmepsin II, China-Myanmar border measurement and thick/thin blood smear preparation and examination were performed at daily for the first 3 days with DHA-PPQ treatment to evaluate the day 3 positive rate. Post-discharge parasite examinations were performed on days 7, 14, 21, 28, 35 and 42. The parasite species and parasitaemia were identified by microscopists certified as Level 1 by the WHO. If the discrepancies of the two values were > 30%, they were re-evaluated by a third microscopist to obtain a final diagnostic consensus (the average of the two close values of the parasitaemia).

Treatment efficacy was evaluated based on clinical and parasitological outcomes and study end points in accordance with WHO guidelines for TES studies (Additional file 1). The outcomes were classified as early treatment failure (ETF), late clinical failure (LCF), late parasitological failure (LPF) and adequate clinical and parasitological response (ACPR). Primary endpoint was ACPR on day 42. PCR genotyping of msp1, msp2 and glurp, comparing day 0 and day of failure samples, was used to differentiate recrudescence (same parasite strain) from new infection with another strain parasite [26] .

Dry blood spots on filter paper (Whatman ™ 903, GE Healthcare, USA) were collected from participants from TES studies and clinical patients by passive surveillance from 2010 to 2014. All the filter paper samples were used for k13 sequencing and pm2 amplification. DNA was extracted from dried blood spots using QIAamp 96 DNA Blood Kit (Valencia, CA, USA).

A nested PCR amplification method (Takara PCR kit) to amplify the k13 gene from codon 433 to 702 (850 bp) as described previously with some minor modifications [27] . The primers for nested PCR, sequencing PCR and cycling conditions for k13 was shown in Additional file 2. PCR products were purified using filter plates (Edge Biosystems, Gaithersburg, MD, USA) and directly sequenced and analyzed on an ABI 3730XL automatic sequencer. The amplification products were analyzed by 1.5% agarose gel electrophoresis before sequencing. Purified products were sequenced by an ABI 3730XL automatic sequencer. Bi-directional sequencing was used and all the products were sequenced twice using independently amplified PCR products. pm2 amplification ASYBR-green based quantitative PCR (ThermoFisher Scientific, Waltham, MA, USA) was used to determine the copy number of pm2 (PF3D7_1408000) as previously described with modifications [24] . The copy numbers were quantified with quantitative PCR thermal cycle at 98 °C for 3 min, followed by 45 cycles at 98 °C for 15 s, 63 °C for 20 s, and 72 °C for 20 s on a C1000 Thermal Cycler (Bio-Rad, Marnes-la-Coquette, France) with the CFX96 Real-Time System (Bio-Rad). All the procedures were executed in triplicate. Plasmodium falciparum β-tubulin (PF3D7_1008700) gene with single copy was used as the internal non-duplicated standard and the P. falciparum 3D7 clone as a parallel 1 copy control. The primers sequences for pm2 and internal standard of P. falciparum β-tubulin were shown in Table 1 . The size of PCR products was 79 bp. The P. falciparum NF54 strain was also included in each amplification as the control. The copy number was calculated by the 2-ΔCt method (ΔCt = Ct PfPM2 −Ct Pf β-tubulin , where Ct is the threshold cycle). Copy number > 1.6 was defined as amplification of the gene. The quantitative PCR primers are listed in Table 1 and protocols are shown in Additional file 3.

The outcome of the therapeutic efficacy study was evaluated with the WHO protocol [28] . The output sequence data of k13 were assembled, edited, and aligned using Sequencher (version 5.1) software. PlasmoDB accession number PF3D7_1343700 (http://www.plasm odb.org) was used as a reference sequence. Fisher's exact test was used for statistical analysis. A regression test was used to estimate the association between the day 3 positive rate and k13 mutation and the relationship of the pm2 copy number variants (CNV) and k13 mutations. Variables with P < 0.05 were considered statistically significant.

This study was approved by the Institutional Review Board of the National Institute of Parasitic Diseases, Chinese Center for Disease Control and Prevention and of the WHO Western Pacific Regional Office. Written informed consent was obtained from patients or guardians. The studies were also registered as clinical trials at https ://www.anzct r.org.au under the numbers ACTRN12610001008011 and ACTRN12610001028099.

A total of 174 participants received directly observed anti-malarial treatment with DHA-PPQ. The treatment outcome of DHA-PPQ was adjusted and corrected by PCR [29] . Table 2 ). One patient from Menglian County in 2014 showed treatment failure after DHA-PPQ treatment on day 35 was confirmed using PCR amplification as P. vivax infection. Another participant from Yingjiang County in 2012 experienced early treatment failure with parasitaemia and axillary temperature ≥ 37.5 °C on day 3 following treatment with DHA-PPQ, but this patient cleared parasitaemia without rescue treatment and was finally classified to be ACPR [26] . Seven of 174 participants (4.0%) were positive for parasitaemia on day 3 after DHA-PPQ treatment.

The k13 gene was successfully sequenced in 268 samples from therapeutic efficacy studies (39.6%, 99/268) and the hospitals by passive surveillance (63.1%, 169/268) ( Table 3) . Twelve non synonymous mutations (F446I, N458Y, Q467H, C469Y, F483S, F495L, P553L, E556D, R561H, A578S, V589G, and Q661R) in the k13 propeller domain were observed in 118 samples, and the prevalence of k13 mutations was 44.0% (118/268). The prevalence of k13 mutation in TES samples and passive surveillance was 46.5% (46/99) and 46.2% (72/169), respectively. F446I was the predominant k13 mutation and accounted for 74.6% (88/118) of all k13 mutations. The k13 mutation from samples collected from Yingjiang and Tengchong had a higher prevalent than from Ruili and Menglian. In addition, five of the seven patients who were positive for parasitaemia on day 3 after DHA-PPQ treatment were observed with k13 mutation (N458Y, C469Y, and P553L). The k13 mutation was not significantly associated with the day 3 positive parasitaemia after DHA-PPQ treatment (P = 0.245).

The 224 samples comprising 91 with k13 mutant and 133 with k13 wild-type were tested on pm2 copy number amplification ( Table 3) . Proportion of different k13 genotype of the samples for pm2 amplification was shown in Fig. 2 . Only one sample was observed with multi-copies of pm2 from k13 mutant samples and other k13 mutant or wild-type samples showed no copy variants of pm2 (99.6%, 223/224). These data were consistent with the therapeutic efficacy of DHA-PPQ. Since there was only one sample with pm2 multiple copies, the data did not meet the conditions for statistical analysis of the association between the pm2 copy number variants and k13 mutations. Nevertheless, it was necessary to note that this sample with multi-copies of pm2 was also carrying k13 mutation.

The emergence and rapid spread of multidrug-resistant P. falciparum across the GMS is a threat to the therapeutic efficacy of DHA-PPQ for treatment of uncomplicated P. falciparum. According to the WHO definition, the partial resistance to artemisinin is characterized as much slower clearance of parasitaemia in the first 3 days of treatment following artemisinin mono-therapy or ACT [30] . Such resistance does not usually lead to treatment failure; however, if the artemisinin component is less effective, the partner drug has to clear a greater parasitic mass, jeopardizing the future efficacy of the partner drug. PPQ resistance was first reported in Cambodia [20] . Recently, alarmingly high rates of treatment failure have occurred in Cambodia, Vietnam, and Thailand, which will have contributed to increased transmission of P. falciparum [21] . In Africa, artemisinin resistance emergence has been reported [31] , but artemisinin (partial) resistance has not been confirmed to date and first-line ACT remains efficacious in all malaria-endemic settings [32] [33] [34] . DHA-PPQ was increasingly deployed as an anti-malarial drugs in Africa. The pm2 duplications was observed in recurrent infections in Mali within 2 months after DHA-PPQ treatment [35] . Another study reported that only a single copy of pm2 was detected in two isolates from Ethiopia and Cameroon after DHA-PPQ failures [36, 37] . These findings indicated the DHA-PPQ failures or pm2 multi-copies in very few cases in Africa, but still raised concerns about the long-term efficacy of DHA-PPQ treatment.

In a previous study, P. falciparum showed markedly delayed clearance following 7 days' treatment with artesunate along the China-Myanmar border, and k13 F446I was the predominant mutation and was associated with delayed parasite clearance [7] . The present study aimed to evaluate the therapeutic efficacy of DHA-PPQ, which is the most widely used artemisininbased combinations in China, as the first-line drug for uncomplicated P. falciparum infection. DHA-PPQ was still efficacy in Southern Yunnan Province with > 95% ACPR outcome. In the 1970s to 1980s, PPQ alone was used for prophylaxis and treatment in the malaria endemic areas in Southern China [15, 38] . In 2009, PPQ was used as a partner drug of artemisinin for first-line treatment of P. falciparum infection in China and now Table 2 Outcome of the therapeutic efficacy study of DHA-PPQ treatment of uncomplicated P. falciparum from 2010 to 2014 ACPR, adequate clinical parasitological response; LCF, late clinical failure; LPF, late parasitological failure; ETF, early treatment failure; WTH, withdrawal; LFU, loss to follow-up a One of 43 participants had parasite and fever on day 3, but this case was continually followed up and cleared parasitaemia without rescue treatment and was finally classified to be ACPR b P. vivax infection was positive on day 35 c Part data of the therapeutic efficacy study of DHA-PPQ has been presented in another study [7] Year DHA-PPQ is the most widely used artemisinin-based combination. Although no resistance to artemisinin and its partner drug PPQ was confirmed in falciparum malaria in Yunnan Province, the efficacy of ACT was decreasing. The prevalence of day 3 positive parasitaemia was > 3% in most of the sentinel sites and artemisinin resistance has been confirmed along the China-Myanmar border [7, 13] . Only one sample was observed with pm2 multi-copies with k13 mutation out of 91 samples and no copy variants of k13 wild-type samples. The association between k13 polymorphisms and pm2 amplification was strong in the Cambodian parasites, which reflected the history of drug selection in Cambodia [24] . The isolates from the China-Myanmar border in our study showed high prevalence of single copy of pm2 compared with the high multi-copies in Cambodia, which also indicated the low connectivity of genetic population between parasites from the China-Myanmar border and Cambodia [39] . Furthermore, the data for pm2 copy number variants provide evidence that no PPQ resistance has emerged along the China-Myanmar border, which is consistent with the therapeutic efficacy of DHA-PPQ.

China is approaching malaria elimination and zero indigenous cases were reported from 2017 [25, 40] . There are still several challenges in the post malaria elimination stage. One big challenge is how to maintain surveillance and response capacity. Monitoring the efficacy of antimalarial drugs for the treatment of imported cases was recommended by the WHO in "A framework for malaria elimination" [41] , which was called integrating drug efficacy surveillance (iDES). China has set up an antimalarial drug surveillance network that is responsible for implementing iDES of DHA-PPQ and anti-malarial molecular markers surveillance throughout the country.

This study only analyzed the pm2 copy number variants of P. falciparum isolates along the China-Myanmar border. Recent studies have identified that mutations in P. falciparum chloroquine resistance transporter (pfcrt) could contribute to a multi-factorial basis of PPQ resistance [42] [43] [44] , which indicates that there are multi-genes or multi-mutations involved in PPQ resistance. Therefore, further molecular surveillance of PPQ resistance including multi-genes will be necessary.

DHA-PPQ for uncomplicated P. falciparum infection was still efficacy in the area with artemisinin-resistant malaria along the China-Myanmar border. There was no evidence to show PPQ resistance in the clinical study and molecular markers survey of pm2. Continued monitoring of the parasite population using molecular markers will be important to track emergence and spread of resistance in this geographic region.

Supplementary information accompanies this paper at https ://doi. org/10.1186/s1293 6-020-03410 -6.

Additional file 1. Classification of treatment outcomes in the therapeutic study.

Additional file 2: Table S1 . Primers for nested PCR and sequencing PCR and cycling conditions for k13.

Additional file 3. The protocol for pm2 copy number amplification.

Geneva: World Health Organization

Artemisinin-based combination treatment of falciparum malaria

Guidelines for the treatment of malaria

A worldwide map of Plasmodium falciparum K13-propeller polymorphisms

Independent emergence of artemisinin resistance mutations among Plasmodium falciparum in Southeast Asia

Spread of artemisinin resistance in Plasmodium falciparum malaria

A single mutation in K13 predominates in Southern China and is associated with delayed clearance of Plasmodium falciparum following artemisinin treatment

Strategy for malaria elimination in the GMS (2015-2030). Geneva: World Health Organization

Status report: artemisinin resistance and artemisinin-based combination therapy efficacy. Geneva: World Health Organization

Guidelines for the treatment of malaria

Technical regulation for application antimalarials of China. Beijing: National Health and Family Planning Commission of China

High antimalarial efficacy of dihydroartemisinin-piperaquine on the China-Myanmar border: the calm before the storm

Clinical efficacy of dihydroartemisinin-piperaquine for the treatment of uncomplicated Plasmodium falciparum malaria at the China-Myanmar border

In vivo monitoring of dihydroartemisinin-piperaquine sensitivity in Plasmodium falciparum along the China-Myanmar border of Yunnan Province

Efficacy of dihydroartemisinin-piperaquine and artemether-lumefantrine in the treatment of uncomplicated falciparum malaria in Hainan

Treatment failure of dihydroartemisinin/piperaquine for Plasmodium falciparum malaria

Ex vivo piperaquine resistance developed rapidly in Plasmodium falciparum isolates in northern Cambodia compared to Thailand

Evidence of Plasmodium falciparum malaria multidrug resistance to artemisinin and piperaquine in Western Cambodia: dihydroartemisinin-piperaquine open-label multicenter clinical assessment

Genetic markers associated with dihydroartemisinin-piperaquine failure in Plasmodium falciparum malaria in Cambodia: a genotype-phenotype association study

Dihydroartemisinin-piperaquine failure in Cambodia

Determinants of dihydroartemisinin-piperaquine treatment failure in Plasmodium falciparum malaria in Cambodia, Thailand, and Vietnam: a prospective clinical, pharmacological, and genetic study

Evolution and expansion of multidrug-resistant malaria in southeast Asia: a genomic epidemiology study

Efficacy of dihydroartemisinin/piperaquine and artesunate monotherapy for the treatment of uncomplicated Plasmodium falciparum malaria in Central Vietnam

A surrogate marker of piperaquine-resistant Plasmodium falciparum malaria: a phenotype-genotype association study

Ready for malaria elimination: zero indigenous case reported in the People's Republic of China

Methods and techniques for clinical trials on antimalarial efficacy: genotyping to identify parasite populations. Geneva: World Health Organization

A molecular mechanism of artemisinin resistance in Plasmodium falciparum malaria

Methods for the surveillance of antimalarial drug efficacy. Geneva: World Health Organization

High sensitivity of detection of human malaria parasites by the use of nested polymerase chain reaction

Antimalarial drug resistance: a threat to malaria elimination. Cold Spring Harb Perspect Med

Emergence of indigenous artemisinin-resistant Plasmodium falciparum in Africa

Absence of putative artemisinin resistance mutations among Plasmodium falciparum in Sub-Saharan Africa: a molecular epidemiologic study

Molecular surveillance for artemisinin resistance in Africa

Scale-up of malaria rapid diagnostic tests and artemisinin-based combination therapy: challenges and perspectives in sub-Saharan Africa

Plasmodium falciparum plasmepsin 2 duplications

Plasmodium falciparum recrudescence two years after treatment of an uncomplicated infection without return to an area where malaria is endemic

Dihydroartemisinin-piperaquine treatment failure in uncomplicated Plasmodium falciparum malaria case imported from Ethiopia

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

In vitro sensitivity of Plasmodium falciparum isolates from China-Myanmar border region to chloroquine, piperaquine and pyronaridine

Genomic structure and diversity of Plasmodium falciparum in Southeast Asia reveal recent parasite migration patterns

Protecting the gains of malaria elimination in China

A framework for malaria elimination. Geneva: World Health Organization

Multigenic architecture of piperaquine resistance trait in Plasmodium falciparum

Association of a novel mutation in the Plasmodium falciparum chloroquine resistance transporter with decreased piperaquine sensitivity

A variant PfCRT isoform can contribute to Plasmodium falciparum resistance to the first-line partner drug piperaquine

Publisher's Note

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

We thank all the participants involved in the clinical study for their contribution work for the study, and local staff from Yingjiang Center for Disease Control and Prevention (CDC), Tengchong CDC, Ruili CDC and Menglian CDC.

This study was presented at the annual meeting of the American Society of Tropical Medicine and Hygiene in Baltimore, MD, November 2017.

Authors' contributions FH, ST and CVP conceived and designed the study; HL, LHT, SZ carried out the TES study; HL collected the passive detection samples; FH, BS, XZ carried out the data analysis; BS and FH conducted the laboratory work. PR, MMN, and CVP provided technique support for the data analysis. FH drafted the manuscript. All authors read and approved the final manuscript. 

The authors declare no competing interests. PR is a staff member of the World Health Organization. The author alone is responsible for the views expressed in this publication and does not necessarily represent the decisions, policy or views of the World Health Organization.