Skip to main content

Limited polymorphism in k13 gene of Plasmodium falciparum and k12 of Plasmodium vivax isolates imported from African and Asian countries between 2014 and 2019 in Hangzhou city, China

Abstract

Background

Malaria causes major public health problems globally and drug resistance hinders its control and elimination. Molecular markers associated with drug resistance are considered as a beneficial tool to monitor the disease trends, evolution and distribution so as to help improve drug policy.

Methods

We collected 148 Plasmodium falciparum and 20 Plasmodium vivax isolates imported into Hangzhou city, China between 2014 and 2019. k13 gene of P. falciparum and k12 of P. vivax were sequenced. Polymorphisms and prevalence of k13 and k12 were analyzed.

Results

Most (98.65%, 146/148) P. falciparum infections were imported from Africa, and half P. vivax cases came from Africa and the other half from Asia. Nucleotide mutation prevalence was 2.03% (3/148) and the proportion of amino acid mutations was 0.68% (1/148). The amino acid mutation, A676S, was observed in an isolate from Nigeria. No mutation of k12 was observed from the parasites from African and Asian countries.

Conclusions

Limited polymorphism in k13 gene of P. falciparum isolates imported from African countries, but no evidence for the polymorphism of k12 in P. vivax samples from African and Asian countries was found. These results provide information for drug policy update in study region.

Peer Review reports

Introduction

Malaria is one of the most important public health problems worldwide. According to the 2020 World Malaria Report, there were an estimated 229 million cases in 2019. Of the five malaria species, Plasmodium falciparum caused the largest number of cases worldwide, and Plasmodium vivax was the most widely distributed [1].

Since the use of antimalarial drugs began in the 1940s, drugs such as chloroquine, mefloquine, quinine, pyrimethamine and sulfadoxine have played an important role in malaria prevention and control. However, with the long-term and widespread use of antimalarial drugs, drug resistance has developed [2, 3]. To date, P. falciparum has developed varying degrees of resistance to all commonly used antimalarial drugs, such as chloroquine, artemisinin and sulfadoxine–pyrimethamine [4,5,6]. P. vivax has also shown resistance to antimalarial drugs [7]. Artemisinin-based combination therapy (ACT) has been suggested as the first-line treatment by the World Health Organization (WHO) since 2006, and its widespread use has led to concerns about artemisinin resistance being raised by researchers and local health authorities. The emergence and development of antimalarial drug resistance not only threatens malaria patients, but also seriously hinders the process of malaria elimination. Therefore, it is important to monitor antimalarial drug resistance, thereby gaining experience on which to base further improvements of treatment protocols.

There are four methods used to monitor drug resistance [8]. Compared to therapeutic efficacy studies, in vitro tests, and measurement of drug concentration, use of molecular markers has the following advantages: resistance-related mutations easily obtained by molecular experiments, avoiding host confounding factors, and blood filter paper samples are more easily transported and stored [8]. Furthermore, multiple samples can be detected in a batch. So, it was preferred as a tool to monitor drug resistance. For P. falciparum, Kelch 13-propeller(k13) mutations have been identified as molecular makers of partial artemisinin resistance [1]. To date, 10 markers have been validated as being associated with resistance to artemisinin: F446I, N458Y, M476I, Y493H, R539T, I543T, P553L, R561H, P574L and C580Y, whereas 11 markers are considered as candidates [1, 9,10,11,12]. Regarding P. vivax, it is assumed that the widespread use of ACT also exerts similar selective pressure on P. vivax parasites [13]. k12, k13 orthologue in P. vivax (here named k12 gene because of its location on chromosome 12), has been investigated by recent studies. Studies from Cambodia and the China–Myanmar border have found k12 mutations (V552I in Cambodia and M124I on China–Myanmar border) [13, 14]. However, another study from Cambodia by longitudinal pooled deep sequencing of k12 gene suggested a lack of selection by artemisinin [15]. Therefore, more evidence is required to clarify its polymorphism in parasite population as well as its potential association with artemisinin resistance.

Hangzhou city, located in the east of China. About thirty malaria cases were imported to Hangzhou each year from Africa and other countries of Asia. To reduce the risk of reemergence due to malaria importation, a highly effective surveillance system was implemented to find each imported case and stop malaria spreading, following national “1-3-7” strategy. And each infection was laboratory-confirmed and epidemiologically investigated. In this study, P. falciparum and P. vivax samples were obtained from imported malaria cases in Hangzhou from 2014 to 2019. k13 gene of P. falciparum and k12 of P. vivax were sequenced to identify and track the prevalence of drug resistance in the countries of origin, so as to provide evidence for refinement of drug policy.

Materials and methods

Study site

Hangzhou city (longitude 118°20′–120°37′, latitude 29°11′–30°34′) is located in the eastern Zhejiang province, China. The dominant mosquito species in Hangzhou city has been Anopheles sinensis over the past 40 years. Since 2010, all the malaria cases in Hangzhou were imported, most of which returned back from Africa, followed by southeast Asia.

Sample collection and DNA extraction

We investigated 148 cases of P. falciparum and 20 of P. vivax malaria that were imported into Hangzhou city between 2014 and 2019. According to the national guideline on malaria diagnosis, individuals with malaria-related symptoms, and positive by microscopy or rapid diagnostic tests (RDTs) would be reported by hospitals or clinics in Hangzhou. Local Center for Disease Control and Prevention carried out epidemiological investigation for each patient, including detailed travelling history to track parasite origins. Infections were double-checked for species by PCR in Zhejiang Provincial Center for Disease Control and Prevention. According to national regulation of antimalarial drugs, for uncomplicated P. falciparum infections, Dihydroartemisinin and Piperaquine phosphate was the first-choice. And for P. vivax patients, Chloroquine Phosphate was the first line drug, followed by Dihydroartemisinin and Piperaquine phosphate. Approximately 1 ml of venous blood was obtained from each patient before treatment. All the blood samples were stored at − 80 °C until use. Genomic DNA was extracted by QIAamp DNA Mini kit (Qiagen Inc., Hilden, Germany).

DNA amplification and sequencing

The propeller region of k13 gene of P. falciparum isolates was amplified by nested PCR as previously described [16], whereas k12 gene of P. vivax isolates was amplified by regular PCR using proof-reading polymerase (see Additional file 1: Table S1). Because of its long sequence, k12 gene was divided into two segments for amplification. Amplification products were then sequenced by Sangon Biotech Co. Ltd. (Shanghai, China). Oligonucleotide primers and cycling conditions are listed in Additional file 1: Table S1). A final 25-μl reaction volume was used.

Data analysis

Nucleotide and amino acid sequences of k13 gene were aligned and compared with reference sequences from NCBI database by Mega version7.0.26 (https://www.megasoftware.net/). The k13 sequence of PF3D7_1343700(MF285413.1) and k12 sequence of Salvador I (Sal-I) genome (PVX_083080) were used as references. A database was constructed using Microsoft Excel 2017. The mutant and wild-type alleles of the collected samples were used to generate the prevalence of the alleles.

Results

General information

We collected 168 samples of malaria infections imported from Africa and other countries of Asia to Hangzhou city, China, including 148 P. falciparum and 20 P. vivax cases (Tables 1 and 2). Specifically, 98.65% (146/148) P. falciparum infections were imported from Africa, while two cases returned back from Asia. Nigeria (31.08%, 46/148), Cameroon (11.49%,17/148), Democratic Republic of the Congo (8.11%, 12/148) and Ghana (6.76%, 10/148) were the main importing countries (Fig. 1). For imported cases of vivax malaria, half came from Africa and half from Asia. Most P. vivax infections (30%,6/20) were imported from Pakistan, followed by Ethiopia (25%, 5/20) (Fig. 2).

Table 1 Distribution of imported P. falciparum cases of Hangzhou city between 2014 and 2019
Table 2 Distribution of imported P. vivax cases of Hangzhou city between 2014 and 2019
Fig. 1
figure 1

Geographical distribution of imported P. falciparum cases from Africa and Asia

Fig. 2
figure 2

Geographical distribution of imported P. vivax cases from Africa and Asia

k13 propeller polymorphisms

The sequences of 148 k13 nested PCR products were successfully obtained from 148 isolates. The k13 SNPs were analyzed by comparing with the reference. Sequence alignment revealed that three samples harbored a nucleotide mutation, G to T at position 2026, in one isolate from Nigeria, C to G at position 1767 in one isolate from Nigeria, and C to T at position 1479 in one isolate from Côte d'Ivoire. The nucleotide mutation prevalence was 2.03% (3/148). Nucleotide mutation at position 2026 resulted in non-synonymous mutation of amino acid A676S, while nucleotide mutations at positions 1767 and 1479 did not bring in non-synonymous amino acid mutations. The prevalence of amino acid mutation was 0.68% (1/148) (Table 3).

Table 3 Prevalence of mutations in the k13 propeller domain of imported P. falciparum infections of Hangzhou city between 2014 and 2019

k12 propeller polymorphisms

k12 of all P. vivax isolates was successfully amplified and sequenced, and compared with the reference. All the k12 of 20 imported P. vivax cases were identical to the reference sequence. Namely, no mutation of k12 was observed from the isolates from Africa and Asia.

Discussion

The current study aimed at providing surveillance information for artemisinin resistance by analyzing the k13 propeller domain in P. falciparum and k12 in P. vivax isolates imported from African and Asian countries between 2014 and 2019 in Hangzhou city, China. Our results showed that both nucleotide and amino acid mutations in P. falciparum isolates were at very low frequencies. As for P. vivax, we demonstrated that k12 in P. vivax isolates was highly conserved.

Mutations in the k13 propeller domain are important determinants of artemisinin resistance [9]. Therefore, k13 marker studies help to identify and track the prevalence of molecular mutations associated with artemisinin partial resistance, as WHO has established a validated list of mutation markers. Since the first discovery of artemisinin resistance in western Cambodia and the border between Cambodia and Thailand in 2008 [5, 6], k13 mutations have become widespread across southeast Asia [17]. In Africa, even slow-clearing infections after ACT treatment were observed at a low frequency of 1%. A578S was the most frequently reported mutation of k13 propeller domain in Africa, and it did not confer artemisinin resistance in vivo or in vitro [18], so the situation concerning artemisinin resistance was thought to be optimistic [19]. However, a recent study reported that an indigenous k13 R561H mutation was identified in 19 samples between 2013 and 2015 from Rwanda and gene editing confirmed that this mutation can drive artemisinin resistance in vitro [10]. It raised international concern for the epidemiology of the k13 marker in Africa. So, the survey primarily explored a diversity of k13 mutations from imported cases imported from Africa between 2014 and 2019. Limited polymorphism of k13 was observed in our study, which indicated that k13 propeller domain of P. falciparum imported from African countries was still relatively conserved. However, the sample size of the current study was small, and more samples from different regions are required for confirmation. Notably, A676S substitution was verified in a 2106 isolate from Nigeria in the current survey. It was previously discovered in Ghana and southeast Asia [20,21,22], but it might be the first time in Nigeria. As for the association between A676S substitution and artemisinin resistance, therapeutic efficacy data from Ghana support that P. falciparum with this mutation is still susceptible to artemether–lumefantrine [20], however, there is still no evidence from in vitro tests. This warrants continuous surveillance to monitor this molecular marker, so as to evaluate its potential impact and association with drug resistance.

Although several publications are showing that short term culture of P. vivax is possible [23, 24], it is still difficult to reveal the underlying mechanisms of antimalarial drug resistance of P. vivax. However, evidence from several studies indicates that multiple molecular markers of P. vivax are involved in drug resistance, which are conferred from homologous genes in P. falciparum, such as Pvdhps and Pvdhfr [25, 26]. As k13 orthologue in P. vivax, k12 is a concern for researchers, although its role in artemisinin resistance is not demonstrated. There are several reports documenting the prevalence of mutant k12 in Cambodia, Thailand, Lao People’s Democratic Republic and the China–Myanmar border [13, 14, 27]. Generally speaking, previous studies have demonstrated that k12 mutation is rare, and all mutations are present at very low prevalence [13, 14, 27]. In our study, polymorphism of k12 was not observed in parasites from Africa or Asia countries, indicating that k12 is conserved in this region.

The current study was limited by a lack of clinical and in vitro data for the association between k13 mutation and artemisinin resistance. Further, there was a small sample size of P. vivax. and the contribution of several sites is limited to one or two isolates, which limited interpretation.

In conclusion, the current study demonstrated limited polymorphism in k13 gene of imported P. falciparum isolates from African countries, and our data presented, for the first time, A676S mutation of k13 in Nigeria. No evidence for the polymorphism of k12 in imported P. vivax samples from African and Asian countries was found. These data might be helpful for updating the evidence-based drug guidelines in study region.

Availability of data and materials

The datasets analyzed in this study are available from the corresponding author on reasonable request. DNA sequences of k13 with the point mutations were deposited in the NCBI database under GenBank accession number MZ170790 to MZ170792.

References

  1. WHO.: World Malaria Report. Geneva: WHO. 2020.

  2. Andrea E, Lehane AM, Jérôme C, Fidock DA. PfCRT and its role in antimalarial drug resistance. Trends Parasitol. 2012;28(11):504–14.

    Article  Google Scholar 

  3. Lakshmanan V, Bray PG, Verdier-Pinard D, Johnson DJ, Horrocks P, Muhle RA, Alakpa GE, Hughes RH, Ward SA, Krogstad DJ, et al. A critical role for PfCRT K76T in Plasmodium falciparum verapamil-reversible chloroquine resistance. EMBO J. 2005;24(13):2294–305.

    Article  CAS  Google Scholar 

  4. Picot S, Olliaro P, Monbrison FD, Bienvenu AL, Price RN, Ringwald P. A systematic review and meta-analysis of evidence for correlation between molecular markers of parasite resistance and treatment outcome in falciparum malaria. Malar J. 2009;8(1):89.

    Article  Google Scholar 

  5. Harald N, Youry S, Kurt S, Smith BL, Duong S, Fukuda MM. Evidence of artemisinin-resistant malaria in western Cambodia. N Engl J Med. 2008;359(24):2619–20.

    Article  Google Scholar 

  6. Dondorp AM, Nosten F, Yi P, Das D, Phyo AP, Tarning J, Lwin KM, Ariey F, Hanpithakpong W, Lee SJ, et al. Artemisinin resistance in Plasmodium falciparum malaria. N Engl J Med. 2009;361(5):455–67.

    Article  CAS  Google Scholar 

  7. Barnadas C, Ratsimbasoa A, Tichit M, Bouchier C, Jahevitra M, Picot S, Menard D. Plasmodium vivax resistance to chloroquine in Madagascar: clinical efficacy and polymorphisms in pvmdr1 and pvcrt-o genes. Antimicrob Agents Chemother. 2008;52(12):4233–40.

    Article  CAS  Google Scholar 

  8. WHO: Global report on antimalaria drug efficacy and drug resistance: 2000–2010.

  9. Ariey F, Witkowski B, Amaratunga C, Beghain J, Langlois AC, Khim N, Kim S, Duru V, Bouchier C, Ma L, et al. A molecular marker of artemisinin-resistant Plasmodium falciparum malaria. Nature. 2014;505(7481):50–5.

    Article  Google Scholar 

  10. Uwimana A, Legrand E. Emergence and clonal expansion of in vitro artemisinin-resistant Plasmodium falciparum kelch13 R561H mutant parasites in Rwanda. Nat Med. 2020;26(10):1602–8.

    Article  CAS  Google Scholar 

  11. WHO: Artemisinin resistance and artemisinin-based combination therapy efficacy. Geneva: WHO 2018.

  12. Straimer J, Gnädig NF, Witkowski B, Amaratunga C, Duru V, Ramadani AP, Dacheux M, Khim N, Zhang L, Lam S, et al. Drug resistance. K13-propeller mutations confer artemisinin resistance in Plasmodium falciparum clinical isolates. Science (80-). 2015;347(6220):428–31.

    Article  CAS  Google Scholar 

  13. Deng S, Ruan Y, Bai Y, Hu Y, Deng Z, He Y, Ruan R, Wu Y, Yang Z, Cui L. Genetic diversity of the Pvk12 gene in Plasmodium vivax from the China–Myanmar border area. Malar J. 2016;15(1):528.

    Article  Google Scholar 

  14. Popovici J, Kao S, Eal L, Bin S, Kim S, Ménard D. Reduced polymorphism in the kelch propeller domain in Plasmodium vivax isolates from Cambodia. Antimicrob Agents Chemother. 2015;59(1):730–3.

    Article  Google Scholar 

  15. Brazeau NF, Hathaway N, Parobek CM, Lin JT, Bailey JA, Lon C, Saunders DL, Juliano JJ. Longitudinal pooled deep sequencing of the Plasmodium vivax K12 Kelch gene in Cambodia reveals a lack of selection by Artemisinin. Am J Trop Med Hyg. 2016;95(6):1409–12.

    Article  CAS  Google Scholar 

  16. Voumbo-Matoumona DF, Kouna LC, Madamet M, Maghendji-Nzondo S, Pradines B, Lekana-Douki JB. Prevalence of Plasmodium falciparum antimalarial drug resistance genes in Southeastern Gabon from 2011 to 2014. Infect Drug Resist. 2018;11:1329–38.

    Article  CAS  Google Scholar 

  17. Ashley EA, Dhorda M, Fairhurst RM, Amaratunga C, Lim P, Suon S, Sreng S, Anderson JM, Mao S, Sam B, et al. Spread of artemisinin resistance in Plasmodium falciparum malaria. N Engl J Med. 2014;371(5):411–23.

    Article  Google Scholar 

  18. Menard D, Khim N, Beghain J, Adegnika AA, Shafiul-Alam M, Amodu O, Rahim-Awab G, Barnadas C, Berry A, Boum Y, et al. A Worldwide map of Plasmodium falciparum K13-propeller polymorphisms. N Engl J Med. 2016;374(25):2453–64.

    Article  Google Scholar 

  19. Dahal P, d’Alessandro U, Dorsey G, Guerin PJ, Nsanzabana C, Price RN, Sibley CH, Stepniewska K, Talisuna AO. Clinical determinants of early parasitological response to ACTs in African patients with uncomplicated falciparum malaria: a literature review and meta-analysis of individual patient data. BMC Med. 2015;13:212.

    Article  Google Scholar 

  20. Aninagyei E, Tetteh CD, Oppong M, Boye A, Acheampong DO. Efficacy of Artemether-Lumefantrine on various Plasmodium falciparum Kelch 13 and Pfmdr1 genes isolated in Ghana. Parasite Epidemiol Control. 2020;11: e00190.

    Article  Google Scholar 

  21. Aninagyei E, Duedu KO, Rufai T, Tetteh CD, Chandi MG, Ampomah P, Acheampong DO. Characterization of putative drug resistant biomarkers in Plasmodium falciparum isolated from Ghanaian blood donors. BMC Infect Dis. 2020;20(1):533.

    Article  CAS  Google Scholar 

  22. Project MPfC. Genomic epidemiology of artemisinin resistant malaria. Elife. 2016. https://0-doi-org.brum.beds.ac.uk/10.7554/eLife.08714

  23. Bermúdez M, Moreno-Pérez DA, Arévalo-Pinzón G, Curtidor H, Patarroyo MA. Plasmodium vivax in vitro continuous culture: the spoke in the wheel. Malar J. 2018;17(1):301.

    Article  Google Scholar 

  24. Gural N, Mancio-Silva L, Miller AB, Galstian A, Butty VL, Levine SS, Patrapuvich R, Desai SP, Mikolajczak SA, Kappe SHI, et al. In vitro culture, drug sensitivity, and transcriptome of Plasmodium Vivax hypnozoites. Cell Host Microbe. 2018;23(3):395-406.e394.

    Article  CAS  Google Scholar 

  25. Nyunt MH, Han JH, Wang B, Aye KM, Aye KH, Lee SK, Htut Y, Kyaw MP, Han KT, Han ET. Clinical and molecular surveillance of drug resistant vivax malaria in Myanmar (2009–2016). Malar J. 2017;16(1):117.

    Article  Google Scholar 

  26. Htun MW, Mon NCN, Aye KM, Hlaing CM, Kyaw MP, Handayuni I, Trimarsanto H, Bustos D, Ringwald P, Price RN, et al. Chloroquine efficacy for Plasmodium vivax in Myanmar in populations with high genetic diversity and moderate parasite gene flow. Malar J. 2017;16(1):281.

    Article  Google Scholar 

  27. Duanguppama J, Mathema VB, Tripura R, Day NPJ, Maxay M, Nguon C, von Seidlein L, Dhorda M, Peto TJ, Nosten F, et al. Polymorphisms in Pvkelch12 and gene amplification of Pvplasmepsin4 in Plasmodium vivax from Thailand, Lao PDR and Cambodia. Malar J. 2019;18(1):114.

    Article  Google Scholar 

Download references

Acknowledgements

We are grateful for the physicians of local center for disease control and prevention for their invaluable assistance in sample collection and epidemiological investigation.

Funding

This work was funded by the Hangzhou agricultural and social development research initiative design project (No. 20190101A12) and Hangzhou health and family planning science and technology project (No. 2018A79).

Author information

Authors and Affiliations

Authors

Contributions

X.Y.J, S.J.Z and W.M.X carried out the molecular studies. J.F.C conceived the study. X.X.W and W.R. analyzed the data and drafted the manuscript. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Wei Ruan or Xiaoxiao Wang.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the Ethical Review Committee of Zhejiang Center for Disease Control and Prevention. All methods were carried out in accordance with relevant guideline and regulation. An informed consent was waived by Ethical Review Committee of Zhejiang Center for Disease Control and Prevention because malaria infection investigation was conducted as part of control and management for national notifiable diseases according to local regulations.

Consent for publication

Not applicable.

Competing interest

The authors have no conflicts of interest to report.

Additional information

Publisher's Note

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

Supplementary Information

Additional file 1: Table S1.

Primers and cycling conditions for K13 and K12 genotyping assay.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Jin, X., Zhu, S., Xu, W. et al. Limited polymorphism in k13 gene of Plasmodium falciparum and k12 of Plasmodium vivax isolates imported from African and Asian countries between 2014 and 2019 in Hangzhou city, China. BMC Infect Dis 21, 853 (2021). https://0-doi-org.brum.beds.ac.uk/10.1186/s12879-021-06579-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s12879-021-06579-6

Keywords