- Open Access
Mutations in Plasmodium falciparum K13 propeller gene from Bangladesh (2009–2013)
Malaria Journal volume 13, Article number: 431 (2014)
Bangladesh is a malaria hypo-endemic country sharing borders with India and Myanmar. Artemisinin combination therapy (ACT) remains successful in Bangladesh. An increase of artemisinin-resistant malaria parasites on the Thai-Cambodia and Thai-Myanmar borders is worrisome. K13 propeller gene (PF3D7_1343700 or PF13_0238) mutations have been linked to both in vitro artemisinin resistance and in vivo slow parasite clearance rates. This group undertook to evaluate if mutations seen in Cambodia have emerged in Bangladesh where ACT use is now standard for a decade.
Samples were obtained from Plasmodium falciparum-infected malaria patients from Upazila health complexes (UHC) between 2009 and 2013 in seven endemic districts of Bangladesh. These districts included Khagrachari (Matiranga UHC), Rangamati (Rajasthali UHC), Cox’s Bazar (Ramu and Ukhia UHC), Bandarban (Lama UHC), Mymensingh (Haluaghat UHC), Netrokona (Durgapur and Kalmakanda UHC), and Moulvibazar (Sreemangal and Kamalganj UHC).
Out of 296 microscopically positive P. falciparum samples, 271 (91.6%) were confirmed as mono-infections by both real-time PCR and nested PCR. The K13 propeller gene from 253 (93.4%) samples was sequenced bi-directionally. One non-synonymous mutation (A578S) was found in Bangladeshi clinical isolates. The A578S mutation was confirmed and lies adjacent to the C580Y mutation, the major mutation causing delayed parasite clearance in Cambodia. Based on computational modeling A578S should have a significant effect on tertiary structure of the protein.
The data suggest that P. falciparum in Bangladesh remains free of the C580Y mutation linked to delayed parasite clearance. However, the mutation A578S is present and based on structural analysis could affect K13 gene function. Further in vivo clinical studies are required to validate the effect of this mutation.
The term ‘artemisinin resistance’ in Plasmodium falciparum has not been precisely defined. The clinical interpretation is a relatively slow parasite clearance rate in patients receiving artemisinin or artemisinin combination therapy (ACT) . In low malaria transmission areas, parasite clearance studies require screening of thousands of febrile individuals to enroll a few patients. Therefore, these studies can be logistically and financially difficult and inconvenient for patients . From a biological perspective, in vitro studies have sought a suitable molecular marker to identify artemisinin resistance in P. falciparum parasites. Mutation and variable expression of several genes such as PfMDR and PfATPase6 have been suggested but not proven . Recently, the K13 propeller gene PF3D7_1343700 (PF13_0238) has been linked to in vitro artemisinin resistance and in vivo slow parasite clearance rate and therefore proposed as marker of artemisinin resistance. K13 propeller has been predicted to consist of three domains 225 amino acid in length: i) Plasmodium-specific and well conserved N-terminal domain; ii) a BTB/POZ domain; and, iii) a six-blade C-terminal propeller domain formed of canonical kelch motifs. Kelch motif containing proteins are evolutionary conserved across different species and grouped into KLHL type proteins. The human host contains 42 of these KLHL-type proteins . KLHL19, also known as KEAP1 protein, has maximum homology with Plasmodium K13 protein . Human KEAP1 protein is a negative regulator of the inducible nuclear erythroid 2-related factor 2 (Nrf2)-dependent cytoprotective response, sequestering Nrf2 in the cytoplasm under steady state [6, 7]. The transcription factor Nrf2 binds to the antioxidant response element (ARE) present in promoters of genes involved in phase II detoxification and oxidative stress responses. Nrf2 forms a heterodimer to activate Maf transcription factor protein that binds to the ARE and activates transcription through the Maf recognition element (MARE). KEAP1 is involved Nrf-2 degradation by targeting it for ubiquitination through the cullin 3 ligase complex. Therefore, it is assumed that the K13 propeller performs a similar function in the P. falciparum and mutation in the gene impairs anti-oxidant/cytoprotective function. However, no orthologue of Nrf2 has been identified in Plasmodium parasite genome .
Bangladesh is a malaria hypo-endemic country, sharing its border with India and Myanmar. Bangladesh has made great strides in malaria control with a 65% decrease in total malaria cases and 91% reduction in malaria-associated mortality from 2008 to 2012. The number of national cases ranged between 26, 891 and 63, 873 between 2009 and 2013 according to the National Malaria Control Programme. High coverage and increased use of insecticide-treated nets, use of rapid diagnostic tests for case detection at the community level, anti-malarial treatments with ACT, and a high number of community health workers and health facilities have contributed significantly to this achievement . Ministry of Health and Family Welfare (MoHFW) has introduced artemether-lumefantrine (Coartem©) as ACT to treat uncomplicated falciparum malaria since 2004  and the process has accelerated with impetus from the Global Fund in 2007 . Artemisinins are considered as the last line of defense against P. falciparum-associated malaria. Since countrywide implementation from 2004, ACT remains quite successful in Bangladesh [10–12]. No day 3 positivity data is available at this time to objectively gauge ACT efficacy. However, the increase of heritable artemisinin-resistant malaria parasite in the Thai-Cambodia border , and subsequently in the Thai-Myanmar  border to southern Myanmar , is worrisome for the Bangladesh-India-Myanmar tri-country border area. Exact population movement in this border area is not defined because huge illegal movements are involved. Rohinga communities of Myanmar have already formed several refugee camps in Bangladesh and still their number is increasing due to civil war in Myanmar. Both transport of resistance via migration and de novo emergence of artemisinin resistance are possible. Previously, chloroquine and sulphadoxine/pyrimethamine-resistant parasites emerged in Thai-Cambodia border, consequently spreading to sub-Saharan Africa across the Asia [16, 17]. Hence, monitoring of artemisinin resistance in Bangladesh is an important global health issue. The current study aimed to identify the presence or absence of mutations in the K13 propeller gene of Bangladeshi P. falciparum clinical isolates in conjunction with structural analysis to assess the functional implications of mutations present.
Plasmodium falciparum- infected malaria patients from corresponding Upazila health complex (UHC) of seven endemic districts of Bangladesh were included in this study. These districts include Khagrachari (Matiranga UHC), Rangamati (Rajasthali UHC), Cox’s Bazar (Ramu and Ukhia UHC), Bandarban (Lama UHC), Mymensingh (Haluaghat UHC), Netrokona (Durgapur and Kalmakanda UHC), and Moulvibazar (Sreemangal and Kamalgonj UHC). Informed consents were obtained from the patients or their legal guardians in the case of children and assent from children aged 11–17 years. A total of 296 P. falciparum positive blood samples confirmed by field and laboratory microscopy in between May 2009 and 2013 was collected for this study. Consent was also taken for future use of blood samples. The study was approved by the institutional Ethics Review Committee (ERC) of International Centre for Diarrhoeal Disease Research, Bangladesh ( icddr,b).
Real time PCR and nested PCR
DNA was extracted by Qiagen Blood mini kit (Carlsbad, CA) from the collected blood specimens. Microscopy results were confirmed by both real time PCR and nested PCR [18, 19]. If any discrepancies were found by different diagnostic methods, those samples were excluded from the study.
K13 propeller gene amplification and sequencing
The K13 propeller gene was amplified by the nested PCR method described by Ariey et al. with some modifications. The first step of the nested PCR remained unaltered while for the second round of PCR, annealing temperature was adjusted to 65.5°C and annealing time was reduced to 45 sec. Platinum Taq Polymerase (Life Technologies, Ontario, Canada) was used for all of the PCR reactions according to its recommended settings. Amplified products were purified by Qiaquick PCR purification kit (Carlsbad, CA) and then sequenced by Applied Biosystems 3730XL 96 capillary sequencer. ClustalW software  was used to align all of the sequences under Unix interface. The sequences in the multiple sequence alignment was edited manually by Jalview software  to remove the gaps. ClustalX colour scheme in the Jalview helped to identify the mismatched bases. The sequence of K13 propeller gene was obtained from Genbank (ID: AL844509.2) and used it as reference sequence to locate the point mutations in the clinical isolates.
Protein structure prediction
Amino acid sequences of the wild type and mutated proteins were used to predict the structure of the protein and effect of mutations found in Bangladeshi clinical isolates. Phyre2 online protein structure prediction tool  was used for modelling. The chosen template (PDB ID: 2WOZ) resulted in a model with a confidence score of 100% based on homology assessment and model prediction quality.
Out of 296 microscopically confirmed P. falciparum samples, 271 (91.56%) were confirmed by both real time PCR and nested PCR to have pure P. falciparum infection. Other (n = 25) specimens were excluded from the study. With nested PCR, 253 (93.35%) amplicons (848 bp) were obtained from the DNA of 271 samples as visualized after SYBR®Safe DNA Gel Stain (Invitrogen) staining of the gel. After ClustalW alignment and Jalview analysis of the retrieved sequences, this study found one synonymous and one non-synonymous mutation in Bangladeshi clinical isolates (Table 1). Geospatial mapping of the districts from which samples were obtained are depicted (Figure 1). The mapping demonstrates that the non-synonymous mutations are in the K13 protein locus of P. falciparum samples (n = 2) from the Matiranga UHC of Khagrachari district. Two samples from unrelated individuals living 15 km apart contained the same A578S mutation in the K13 propeller gene. Genotyping of the isolates using msp1 and msp2 demonstrated the isolates were genetically distinct (data not shown). Matiranga UHC is close to the border with India.
Computational modeling was performed using the beta-propeller domain of the btb-kelch protein krp1 (PDB ID: 2WOZ). The resulting model (Figure 2a) resembled closely the topology and statistics of the predicted model by Ariey  as assessed by sequence alignment to the coordinates of the template and prediction quality. The model predicts a 6-bladed propeller structure with homology to previously characterized kelch domain proteins [23, 24]. Mutational sensitivity prediction by the PHYRE2 Investigator suggests that the A578S mutation observed in this study is likely to disrupt the P. falciparum K13 propeller protein function. K13-propeller protein mutations are reported to alter the biological function of the protein by modification of surface charges that disrupt interactions with other proteins. Superimposition of the homology models of the wild type and mutated proteins show qualitatively that A578S mutant (in addition to the known mutant C580Y seen in South East Asia) is predicted to cause structural changes in the overall topology of the protein (Figure 2b). Due to these changes, the function of the protein with these mutations is quite likely disrupted. Indeed, these modeling predictions are consistent with the findings of the study by Ashley  who report that mutations of the K13 propeller protein after amino acid position 440 were associated with a mean increase in the parasite clearance half-life of 116% (95% CI, 103 to 131; P < 0.001) .
Although the history of ACT use in Bangladesh dates to only 2004, the country lies adjacent to countries where resistance to ACT is emerging. Parasites from South East Asian countries appear to have a propensity to become multi-drug resistant . Researchers have observed that mutations such as M476I in the K13 propeller gene are linked to artemisinin resistance in vitro using the ring stage assay (RSA) developed specifically to test artemisinin activity in the test tube . Clinical studies have also confirmed 17 mutations in the K13 propeller gene were associated with delayed parasite clearance in Cambodia , contrasting data showing high efficacy of ACT in Bangladesh [10–12] in the Bandarban district. This study has collected blood specimens from different parts of the country and evaluated the presence of K13 propeller mutations. Two (0.78%) specimens from unrelated individuals appear to contain structure-altering mutations (A578S) in the K13 propeller gene from genetically distinct isolates. These data suggest that mutations in K13 arise de novo rather through clonal expansion. The Matiranga UHC in the district of Khagrachari is in the tri-border area with Myanmar and India. There is no report of artemisinin resistance from India so far. While ACT remains efficacious and these mutations are rare based on this surveillance, ACT usage will impose positive selection on the resistant parasite. It has been suggested that the high diversity of vector species in Bangladesh [27, 28] can contribute to the spread of these mutated strains of P. falciparum. Moreover, a very low rate of transmission can result in inadequate immunity to clear the parasite that might have survived after ACT treatment, thereby increasing the selection pressure on resistant clones . Monitoring of possible resistant parasite clones is of paramount importance. A limitation of this study is the lack of clinical efficacy data. Also no other genetic markers of resistance were evaluated.
Protein modelling data suggest that A578S can alter the function of the K13 propeller protein. A578S is very close to the C580Y mutation, the major mutation observed in the K13 propeller protein of the Cambodian isolates . The tracking resistance to artemisinin collaboration (TRAC) group also found the A578S mutation in one specimen within their study areas but this mutation was not associated with increased parasite clearance half-life . Importantly, mutations such as A578S lie after the position 400 and alter the structure of K13 significantly and thus may have an effect on artemisinin activity. Crude measures like parasite clearance time are multi-factorial phenotypes influenced by pharmacokinetics, immunity, and patient adherence, and thus may not immediately correlate with important mutations without a large sample size. Ramu UHC in Cox’s Bazar district is one of the TRAC study sites where slow clearing parasites were observed without any mutation on K13 propeller gene which again speaks to the multi-factorial nature of this phenotype. Similarly, clinical specimens collected from Ramu in the study presented here did not demonstrate mutations in the K13 gene. Takala-Harrison  recently demonstrated that the vast majority of isolates in their collection were wild type at the K13 locus with A578S not present.
Up-to-date information on whether artemisinin resistance has already disseminated or independently emerged beyond the Cambodia-Thailand region is a significant issue. The WHO, with the help numerous collaborators, has declared the global plan for artemisinin resistance containment . If resistance were limited to a small, well-defined area, then elimination of drug-resistant malaria from these regions would be possible . Unfortunately, resistance is not limited to the Thai-Cambodian border, it has already spread to the Vietnam, Laos and Myanmar border . Our data demonstrate that mutations strongly linked to artemisinin resistance are not present at this time in Bangladesh.
Dondorp AM, Nosten F, Yi P, Das D, Phyo AP, Tarning J, Lwin KM, Ariey F, Hanpithakpong W, Lee SJ: Artemisinin resistance in Plasmodium falciparum malaria. New Engl J Med. 2009, 361: 455-467. 10.1056/NEJMoa0808859.
Witkowski B, Amaratunga C, Khim N, Sreng S, Chim P, Kim S, Lim P, Mao S, Sopha C, Sam B: Novel phenotypic assays for the detection of artemisinin-resistant Plasmodium falciparum malaria in Cambodia: in-vitro and ex-vivo drug-response studies. Lancet Infect Dis. 2013, 13: 1043-1049. 10.1016/S1473-3099(13)70252-4.
O’Brien C, Henrich PP, Passi N, Fidock DA: Recent clinical and molecular insights into emerging artemisinin resistance in Plasmodium falciparum. Curr Opin Infect Dis. 2011, 24: 570-10.1097/QCO.0b013e32834cd3ed.
Dhanoa BS, Cogliati T, Satish AG, Bruford EA, Friedman JS: Update on the Kelch-like (KLHL) gene family. Hum Genomics. 2013, 7: 13-10.1186/1479-7364-7-13.
Ariey F, Witkowski B, Amaratunga C, Beghain J, Langlois AC, Khim N, Kim S, Duru V, Bouchier C, Ma L, Lim P, Leang R, Duong S, Sreng S, Suon S, Chuor CM, Bout DM, Menard S, Rogers WO, Genton B, Fandeur T, Miotto O, Ringwald P, Le Bras J, Berry A, Barale JC, Fairhurst RM, Benoit-Vical F, Mercereau-Puijalon O, Menard D: A molecular marker of artemisinin-resistant Plasmodium falciparum malaria. Nature. 2014, 505: 50-55.
Itoh K, Wakabayashi N, Katoh Y, Ishii T, Igarashi K, Engel JD, Yamamoto M: Keap1 represses nuclear activation of antioxidant responsive elements by Nrf2 through binding to the amino-terminal Neh2 domain. Gene Dev. 1999, 13: 76-86. 10.1101/gad.13.1.76.
Zhang DD, Hannink M: Distinct cysteine residues in Keap1 are required for Keap1-dependent ubiquitination of Nrf2 and for stabilization of Nrf2 by chemopreventive agents and oxidative stress. Mol Cell Biol. 2003, 23: 8137-8151. 10.1128/MCB.23.22.8137-8151.2003.
Haque U, Overgaard HJ, Clements ACA, Norris DE, Islam N, Karim J, Roy S, Haque W, Kabir M, Smith DL: Malaria burden and control in Bangladesh and prospects for elimination: an epidemiological and economic assessment. Lancet Glob Health. 2014, 2: e98-e105. 10.1016/S2214-109X(13)70176-1.
JICPD: Joint International Continued Professional Development Meeting. Bangladesh College of Physicians and Surgeons: 22–23 February, 2004. 2004, Dhaka: BCPS
Haque R, Thriemer K, Wang Z, Sato K, Wagatsuma Y, Salam MA, Akther S, Akter J, Fukuda M, Miller RS: Therapeutic efficacy of artemether-lumefantrine for the treatment of uncomplicated Plasmodium falciparum malaria in Bangladesh. Am J Trop Med Hyg. 2007, 76: 39-41.
Starzengruber P, Swoboda P, Fuehrer H-P, Khan WA, Hofecker V, Siedl A, Fally M, Graf O, Teja-Isavadharm P, Haque R: Current status of artemisinin-resistant falciparum malaria in South Asia: a randomized controlled artesunate monotherapy trial in Bangladesh. PLoS One. 2012, 7: e52236-10.1371/journal.pone.0052236.
Takala-Harrison S, Clark TG, Jacob CG, Cummings MP, Miotto O, Dondorp AM, Fukuda MM, Nosten F, Noedl H, Imwong M: Genetic loci associated with delayed clearance of Plasmodium falciparum following artemisinin treatment in Southeast Asia. Proc Natl Acad Sci U S A. 2012, 110: 240-245.
Phyo AP, Nkhoma S, Stepniewska K, Ashley EA, Nair S, McGready R, Al-Saai S, Dondorp AM, Lwin KM, Singhasivanon P: Emergence of artemisinin-resistant malaria on the western border of Thailand: a longitudinal study. Lancet. 2012, 379: 1960-1966. 10.1016/S0140-6736(12)60484-X.
Phompradit P, Muhamad P, Wisedpanichkij R, Chaijaroenkul W, Na-Bangchang K: Four years’ monitoring of in vitro sensitivity and candidate molecular markers of resistance of Plasmodium falciparum to artesunate-mefloquine combination in the Thai-Myanmar border. Malar J. 2014, 13: 23-10.1186/1475-2875-13-23.
Kyaw MP, Nyunt MH, Chit K, Aye MM, Aye KH, Lindegardh N, Tarning J, Imwong M, Jacob CG, Rasmussen C: Reduced susceptibility of Plasmodium falciparum to artesunate in southern Myanmar. PLoS One. 2013, 8: e57689-10.1371/journal.pone.0057689.
Roper C, Pearce R, Nair S, Sharp B, Nosten F, Anderson T: Intercontinental spread of pyrimethamine-resistant malaria. Science. 2004, 305 (5687): 1124-10.1126/science.1098876.
Verdrager J: Epidemiology of the emergence and spread of drug-resistant falciparum malaria in South-East Asia and Australasia. J Trop Med Hyg. 1986, 89 (6): 277-289.
Alam MS, Mohon AN, Mustafa S, Khan WA, Islam N, Karim MJ, Khanum H, Sullivan DJ, Haque R: Real-time PCR assay and rapid diagnostic tests for the diagnosis of clinically suspected malaria patients in Bangladesh. Malar J. 2011, 10: 175-10.1186/1475-2875-10-175.
Snounou G, Viriyakosol S, Zhu XP, Jarra W, Pinheiro L, Do Rosario VE, Thaithong S, Brown KN: High sensitivity of detection of human malaria parasites by the use of nested polymerase chain reaction. Mol Biochem Parasit. 1993, 61: 315-320. 10.1016/0166-6851(93)90077-B.
Thompson JD, Higgins DG, Gibson TJ, Clustal W: Improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22: 4673-4680. 10.1093/nar/22.22.4673.
Clamp M, Cuff J, Searle SM, Barton GJ: The jalview java alignment editor. Bioinformatics. 2004, 20: 426-427. 10.1093/bioinformatics/btg430.
Kelley LA, Sternberg MJE: Protein structure prediction on the Web: a case study using the Phyre server. Nat Prot. 2009, 4: 363-371. 10.1038/nprot.2009.2.
Boyden LM, Choi M, Choate KA, Nelson-Williams CJ, Farhi A, Toka HR, Tikhonova IR, Bjornson R, Mane SM, Colussi G: Mutations in kelch-like 3 and cullin 3 cause hypertension and electrolyte abnormalities. Nature. 2012, 482: 98-102. 10.1038/nature10814.
Padmanabhan B, Tong KI, Ohta T, Nakamura Y, Scharlock M, Ohtsuji M, Kang M-I, Kobayashi A, Yokoyama S, Yamamoto M: Structural basis for defects of Keap1 activity provoked by its point mutations in lung cancer. Mol Cell. 2006, 21: 689-700. 10.1016/j.molcel.2006.01.013.
Ashley EA, Dhorda M, Fairhurst RM, Amaratunga C, Lim P, Suon S, Sreng S, Anderson JM, Mao S, Sam B: Spread of artemisinin resistance in Plasmodium falciparum malaria. N Engl J Med. 2014, 371: 411-423. 10.1056/NEJMoa1314981.
White NJ: Antimalarial drug resistance. J Clin Invest. 2004, 113: 1084-1092. 10.1172/JCI21682.
Alam MS, Chakma S, Khan WA, Glass GE, Mohon AN, Elahi R, Norris LC, Podder MP, Ahmed S, Haque R, Sack DA, Sullivan DJ, Norris DE: Diversity of anopheline species and their Plasmodium infection status in rural Bandarban, Bangladesh. Parasit Vectors. 2012, 5: 150-10.1186/1756-3305-5-150.
Alam MS, Khan MG, Chaudhury N, Deloer S, Nazib F, Bangali AM, Haque R: Prevalence of anopheline species and their Plasmodium infection status in epidemic-prone border areas of Bangladesh. Malar J. 2010, 9: 15-10.1186/1475-2875-9-15.
Takala-Harrison S, Jacob CG, Arze C, Cummings MP, Silva JC, Dondorp AM, Fukuda MM, Hien TT, Mayxay M, Noedl H: Independent emergence of artemisinin resistance mutations among Plasmodium falciparum in Southeast Asia. J Infect Dis. 2014, jiu491-Epub ahead of print
World Health Organization: Global Plan for Artemisinin Resistance Containment (GPARC). 2011, Geneva: WHO
Maude RJ, Pontavornpinyo W, Saralamba S, Aguas R, Yeung S, Dondorp AM, Day NPJ, White NJ, White LJ: The last man standing is the most resistant: eliminating artemisinin-resistant malaria in Cambodia. Malar J. 2009, 8: 31-10.1186/1475-2875-8-31.
The authors are indebted to the people of Bangladesh who consented to participate in the study and the doctors and staff of different UHCs for their extended support. The authors appreciate the valuable contribution to the study of Md Khaja Mohiuddin, Shariar Mustafa, Rubayet Elahi and H M Al-Amin.”.
The authors declare that they have no competing interests.
ANM, MSA, RH, and DRP conceptualized and designed the study, collected and identified samples, analysed the data and drafted the manuscript. DS did the protein modelling. ANM, AF and AGB performed sample analyses and made critical revisions of the manuscript. All the authors read the final version of the manuscript and approved.
Abu Naser Mohon, Mohammad Shafiul Alam contributed equally to this work.