ABSTRACT
The continuous surveillance of polymorphisms in the kelch propeller domain of Plasmodium falciparum from Africa is important for the discovery of the actual markers of artemisinin resistance in the region. The information on the markers is crucial for control strategies involving chemotherapy and chemoprophylaxis for residents and nonimmune travelers to the country. Polymorphisms in the kelch propeller domain of Ghanaian malaria parasites from three different ecological zones at several time periods were assessed. A total of 854 archived samples (2007 to 2016) collected from uncomplicated malaria patients aged ≤9 years old from 10 sentinel sites were used. Eighty-four percent had wild-type sequences (PF3D7_1343700), while many of the mutants had mostly nonsynonymous mutations clustered around codons 404 to 650. Variants with different amino acid changes of the codons associated with artemisinin (ART) resistance validated markers were observed in Ghanaian isolates: frequencies for I543I, I543S, I543V, R561P, R561R, and C580V were 0.12% each and 0.6% for R539I. Mutations reported from African parasites, A578S (0.23%) and Q613L (0.23%), were also observed. Three persisting nonsynonymous (NS) mutations, N599Y (0.005%), K607E (0.004%), and V637G (0.004%), were observed in 3 of the 5 time periods nationally. The presence of variants of the validated markers of artemisinin resistance as well as persisting polymorphisms after 14 years of artemisinin-based combination therapy use argues for continuous surveillance of the markers. The molecular markers of artemisinin resistance and the observed variants will be monitored subsequently as part of ongoing surveillance of antimalarial drug efficacy/resistance studies in the country.
INTRODUCTION
Malaria parasite resistance to antimalarial drugs is of paramount concern in control, elimination, and eradication programs throughout the world. Plasmodium falciparum resistance to artemisinin (ART) and its derivatives is well documented in Southeast Asia (SEA), but comparatively low-level ART resistance has been identified in Africa (1). Ongoing worldwide surveillance is necessary considering the potential tremendous public health impact resistance could have, especially on children under 5 years and pregnant women from Africa as well as nonimmune travelers (2). Artemisinin resistance as defined by the WHO is the slow parasite clearance rate with the use of artemisinin-based combination therapy (ACT) or microscopically detected persistent parasites after 3 days of ACT administration (2). Although the exact mechanism of parasite ART resistance is not clearly understood, genetic factors have been implicated for other antimalarials such as chloroquine (CQ; P. falciparum chloroquine resistance transporter gene pfcrt) and sulfadoxine pyrimethamine (SP; P. falciparum dihydrofolate reductase gene pfdhfr and P. falciparum dihydropteroate synthase gene pfdhps) (3, 4).
ART resistance in the SEA region has been linked to the kelch propeller domain on chromosome 13 (pfk13) (3–6). Single nucleotide polymorphisms (SNPs) and consequent amino acid changes in the gene that occurred as a result of drug pressure are used as molecular markers (7, 8). The pfk13 markers of resistance are nonsynonymous and they include N458Y, R539T, E556D, P574L, R575K, C580Y, S621F, Y493H, R539T, and I543T (3, 9). However, the principal mutations are C580Y, Y493H, R539T, I543T, and N458Y, which were observed in all isolates with the slow clearance phenotypic trait. To date, the principal mutations observed in African isolates are C580Y seen in two samples from Cameroon and Y493H seen in one sample from Ghana out of 1,648 samples from Africa (10). Other mutations that have been observed in ART-resistant parasites from the SEA region and have also been prevalent in African isolates include S522C, P553L, R561H, A675V, and H719N (10, 11). The predominant pfk13 mutation found in African isolates is A578S, which has been observed in all the African countries where the pfk13 mutations have been typed (10, 12–18). Since it is close to the principal mutation C580Y, it is alleged that A578S may be the mutant for ART resistance in African isolates (12); however, there is no phenotypic association with ART resistance (8). Of all the pfk13 mutations observed in African isolates, none has yet been directly linked to ART resistance; thus, further investigations are needed.
In Ghana, the malaria treatment policy was changed in 2005 due to P. falciparum resistance to CQ and SP (19). The ACTs used for uncomplicated malaria treatment are artesunate-amodiaquine (AA), which was started in 2005, and artemether-lumefantrine (AL) as well as dihydroartemisinin piperaquine (DHAP), which were subsequently added in 2008 (20). Recently, reports indicate delayed parasite clearance in isolates from the northern parts of Ghana, and the nationwide cure rates are 96.0% for AL and 99.2% for AA (21). Studies on the molecular markers of ART resistance showed increased P. falciparum multidrug resistance gene (pfmdr1) copy number in 18% of isolates and an increasing trend in the prevalence of haplotype pfmdr1 N86-F184-D1246 (linked to AL resistance) from 2003 to 2010 (22). In addition, the P. falciparum adaptor protein complex 2 gene (pfap2mu) S106N and ubiquitin specific protease 1 gene (pfubp1) E1528D and D1525E mutations, which have been linked to the delayed clearance of parasites to ACT in Kenya and recurrent imported malaria in Britain, were observed in 7.4%, 7.4%, and 4.9%, respectively, of Ghanaian parasites (23–25). This is indicative of the possibility of subtle levels of ART resistance in circulating parasites. As part of ongoing surveillance for molecular markers of ART resistance in Ghana, this study reports the diversity and the prevalence of both novel and known mutations in the pfk13 gene in Ghanaian P. falciparum isolates.
RESULTS
The information on the molecular markers of antimalarial drug resistance (a tool for monitoring drug resistance, which gives early warning of the development/emergence of drug resistance) is crucial for control strategies involving chemotherapy and chemoprophylaxis for residents and nonimmune travelers to countries where malaria is endemic. Therefore, polymorphisms in pfk13 of Ghanaian malaria parasites from three different ecological zones at several time periods were assessed. A total of 1,100 samples were sequenced for the pfk13 gene. Of these, 854 isolates with good sequences were used for the analysis. The majority of samples were from the 2015 to 2016 collection, because stored samples from prior years were limited due to usage in other studies. The samples contributed per year from sites in the three ecological zones are shown in Table 1.
Ecological zones and year of collection of archived samples at five time periods
Distribution of pfk13 wild type and mutants in the ecological zones.The majority of isolates 84% (716/854) had wild-type sequences like the 3D7 strain (PF3D7_1343700, PLASMODB). The mutants (16%) had SNPs clustered around codons 404 to 650 of the gene that were mostly nonsynonymous (NS) mutations (77.0%). Most of the mutations occurred only once in the samples. Multiple SNPs were observed in some of the samples, which were due to the multiple clones in Ghanaian malaria infections and was obvious from peaks shown in the DNA sequence chromatograms. Seventeen (12.4%) of the 138 samples with mutations had at least two mutations. For the ecological zones, the distributions of wild-type and mutant alleles in the isolates are shown in Fig. 1. Using the pooled data from Ghana, the proportions of wild-type alleles for the 5 time periods showed a significant declining trend (χ2 = 33.9, P = 0.001) over time. Figure 2 shows the distributions of wild-type and mutant alleles for Ghana over the five time periods.
The distribution of pfk13 wild-type and mutant alleles in parasite populations from three ecological zones in Ghana from 2007 to 2016. The orange shade is the proportion for wild-type alleles (wt) and the yellow shade for mutant alleles (mt). (Ghana map adapted with permission from www.netmaps.net.)
The distribution of wild-type and mutant alleles for the study time periods. The gray line depicts a decreasing trend of pfk13 wild-type alleles with time in Ghanaian parasite populations.
Novel pfk13 mutations in Ghanaian isolates.The mutations observed were diverse, and no particular mutation was exceptionally high in frequency in the population over the years. Likewise, no haplotypes were observed, and the majority of the mutations were observed once. Unique mutations with variants of both synonymous (SYN) and NS mutations were also observed. At codon 607 with amino acid lysine (K), variants of NS mutations, K607E (4), K607H (2), K607N (2), K607I (1), K607Q (1), and K607* (* indicates a stop codon) (1) were observed in 11 samples from 2010, 2012, 2013 to 2014, and 2015 to 2016. Another codon with two variants, G496G (13) and G496W (1), was seen in 14 samples from the forest zone in 2015 to 2016. In addition, K610 had three mutant variants, K610R (5), K610E (4), and K610K (1), in 10 samples from 2010, 2012, and 2016. Both SYN and NS mutations which occurred in more than two samples are shown in Fig. 3.
The pfk13 synonymous and nonsynonymous mutations that occurred in more than two samples; n is the number of samples with the mutations.
Other mutations were observed in specific ecological zones. The mutations seen from all zones were L422F, N408S, N599Y, S466N, and Y511*. Those seen in samples from forest and Guinea savannah zones were K607E, R539I, C542W, and R515I, while A627S seen in Guinea and coastal savannah samples. Those seen only in the forest zone were E606D, G496G, and L407R. The frequencies of the novel mutations and the time periods for the ecological zones are shown in Table 2. There were unique mutations to specific sites: Wa, I461I (2 samples); Hohoe, R528T (3 samples) and D421N (2 samples); Cape Coast, Q613P (2 samples); Begoro, A578S (2 samples); Navrongo, I590T (2 samples); and Bekwai, D584D (2 samples). The persisting NS mutations observed at three time periods include N599Y (2010, 2012, and 2013 to 2014), K607E (2010, 2013 to 2014, and 2015 to 2016), and V637G (2010, 2013 to 2014, and 2015 to 2016), and the zones in which they were observed are shown in Table 2.
Frequencies of pfk13 mutations unique to or shared by the ecological zones at different time periods
ART resistance validated and candidate pfk13 mutations in Ghanaian isolates.The codons bearing the mutations linked to artemisinin resistance in the SEA region were observed with different amino acid substitutions in P. falciparum strains from Ghana. Some of the mutations were found at specific sentinel sites and time periods. They include C580V, detected in one sample from Cape Coast, R539I, detected in 5 isolates from Hohoe (1), Begoro (1), Sunyani (1), and Wa (2), variants of codon 543, I543I, I543S, and I543V, in single isolates from Navrongo, Sunyani, and Hohoe, respectively, R561P from Hohoe, and R561R from Sunyani. The M476I mutation (found in the Tanzanian ART resistant isolate, African descent) was seen in 1 sample from Wa, and a variant mutant, M476V, was detected in Begoro.
The effect of amino acid changes in proteins using the method by Brick and Pizzi (26), CCF53_62 matrices, indicates that the smaller the stabilization score, the more negative and drastic the change is on the proteins. Referring to the CCF53_62 matrices, the change from C to Y in C580Y gives a stabilization score of −2, and the change from C to V in C580V gives a stabilization score of −1. This implies that both amino acid changes have effects on the protein structure and function. However, C580Y will have a more drastic change due to its more negative stabilization score. The change from R to T in R539T gives a stabilization sore of −1, and the change from R to I in R539I observed in Ghanaian isolates gives a stabilization score of −3, which is much lower, implying that the change from R to I may have a graver effect on the protein structure and function. Similarly, the change from I to T in I543T gives a score of −1, whereas the change from I to S at the same codon is −2.
Previously reported African pfK13 mutations observed in Ghanaian isolates.This study detected mutations that were previously described from African malaria parasites which are likely candidates for markers of drug resistance. The A578S seen in both SEA and African isolates was seen in 2 isolates from Begoro, while the additional variants A578P and A578A were seen in isolates from Wa and Sunyani, respectively. Variants of most of the reported mutations from African countries were observed except for C469C in 3 isolates from Begoro, Wa, and Navrongo, M608L and N629Y in samples from Cape Coast, and C532C from Begoro. The validated mutations as well as other candidate mutations for ART resistance in Africa reported in published articles and the variants observed in Ghanaian isolates are shown in Table 3.
Variants of ART resistance validated and candidate mutations in Ghanaian isolates
DISCUSSION
The pfk13 mutations linked to ART resistance in the SEA region are being used as a paradigm for ongoing surveillance to detect emerging malaria parasite resistance to ACTs in countries in sub-Saharan Africa (sSA) where malaria is endemic. With the use of ACTs in Ghana for the past 14 years, there is the growing concern that drug-resistant parasites are evolving and are under selective pressure. This study identified pfk13 polymorphisms in Ghanaian P. falciparum populations from uncomplicated malaria cases reporting at health centers over 9 years (2007 to 2016). Analysis of the pfk13 genetic diversity in parasites from across 3 ecological zones in Ghana uncovered both known and novel mutations, with the majority being nonsynonymous. Although none of the validated mutations for ART resistance were observed, variants with different amino acid substitutions of the codons conferring resistance were observed (R539I instead of T). In addition, some candidate mutations from African isolates were also observed but at low frequencies, such as A578S.
The observed diverse mutations in the Ghanaian isolates corroborate what was observed by Kamau and others in their baseline report of pfk13 polymorphisms in sub-Saharan Africa (12). Although the period of sample collection from the sites was not disclosed, it is noteworthy that the majority of mutations observed in that study were from Ghanaian parasites (Cape Coast and Navrongo). The paper also mentioned the polyclonality of infections observed in samples from Ghana and Kenya, which was apparent by the presence of minor peaks in the DNA sequence chromatograms (12). A similar observation as seen in Ghanaian infections was made in some of the samples analyzed in this study due to the multiplicity of infection (MOI). Depending on the region, the MOIs ranged from 1 to 4 genotypes per infection (27), as such, some samples had more than one pfk13 mutation. This observation is in contrast to the often seen monoclonal infections in the SEA region, where one mutation per sample was observed for the parasites (3). Ocan and others have also indicated that there is a higher diversity of pfk13 mutations in the African region than in the SEA region from a systematic review of the prevalence of the mutations in the two geographic malaria regions (28).
After 2 to 3 years of ACT implementation (2007 to 2008), the wild-type population for pfk13 was approximately 99%, and subsequently, there has been an overall significant decline in the ecological zones. The Guinea savannah (GS) and the coastal savannah (CS) zones had no mutants in the first time period (2007 to 2008), but the mutants were seen from the second time period onwards. This implies that there was an independent emergence of spontaneous mutations with time as observed in Senegal, where no mutations were detected in 2010 but started coming up in the subsequent years (29, 30). There was also an observation of no mutations in the 2012 forest and 2013 to 2014 CS samples, which may be due to the samples available for use (22 and 24, respectively). It is possible that most of the samples with mutants may have been missed due to usage in other studies. However, the overall analysis showed a decline the pfk13 wild-type population with time, suggesting the evolution of the pfk13 C-terminal region with the introduction of ACTs due to the increase in mutations over time. However, there are no prior data to ACT use to confirm this assertion. Three mutations were also observed at three consistent time periods nationally but were seen at 2 time periods for the zones. This is indicative of a possible selection of the mutations with drug pressure, but the frequencies decreased with time. These mutations will be tracked in the subsequent years to investigate the spread in the population.
The observation of mutations at the codons of the ART resistance molecular markers but with different amino acid substitutions is of great interest. The variants of the validated mutations seen include R539I, I543S, R561P, and C580V. The CCF53_62 matrix analysis by Brick and Pizzi (26) which determines the effect of amino acid changes in proteins using bioinformatics indicates that the smaller the stabilization score of the amino acid, the more negative and drastic the change in the protein. From the analysis, R539I, I543S, R561P, and C580V had lower stabilization scores than R539T, I543T, R561H, and C580Y, respectively, which implies an effect on the protein structure and function. The functional roles of the substituted amino acids need to be investigated to reveal any disruption or modification of the domain scaffold and the phenotypic effect on parasite artemisinin resistance. It is worth noting that since the molecular markers of CQ and SP for the SEA region and sSA regions were of different haplotypes, there is a likelihood of a similar scenario with pfk13 mutations.
Differences in the number and distribution of the polymorphisms from the different geographic locations depict the probable role of transmission dynamics as well as and human migratory patterns. High transmission intensity enhances genetic recombination in parasite populations and consequently spontaneous mutations that could be heritable. Although there were mutations that were shared among populations from the different ecological zones, others were zone specific. This implies that the possibility of intracountry gene flow as a result parasite transport with human migratory patterns cannot be overlooked, which has been observed in The Gambia (Alfred Amambua-Ngwa, personal communication). The observation of previously reported mutations in circulating African parasites corroborates an assertion that mutations from the African parasites probably independently evolved locally; however, the possibility of gene flow from the SEA region cannot be overlooked, although there is no evidence of that yet (10).
The data show that after 14 years of the introduction of ACT use in Ghana, there has been a steady increase in the variability of mutations in the pfk13 gene as well as the number of mutants from 2007 to 2016. The pfk13 polymorphisms observed were both SYN and NS mutations which include known mutations from SEA and Africa malaria regions, but in most cases, with a different in amino acid substitution. Novel mutations were also detected in Ghanaian isolates, and three of them were observed consistently at three time periods nationally. In the advent of the unknown markers for ACT resistance in Africa, the observed mutations will be tracked in isolates from ongoing surveillance of antimalarial drug resistance in Ghana for the possible validation of the resistance-conferring mutations in Africa.
MATERIALS AND METHODS
Study sites.A collaboration between the National Malaria Control Program (NMCP) and the Noguchi Memorial Institute for Medical Research (NMIMR) established 10 sentinel sites in the ten regions of Ghana for the surveillance of antimalarial drug efficacy in the country (Fig. 4). The sites include Navrongo (10.8940°N, 1.0921°W), Wa (10.0601°N, 2.5099°W), and Yendi (9.4450°N, 0.0093°W) in the Guinea savannah ecological zone with seasonal malaria transmission, Begoro (6.3916°N, 0.3795°W), Bekwai (6.4532°N, 1.5838°W), Hohoe (7.1519°N, 0.4738°E), Sunyani (7.3349°N, 2.3123°W), and Tarkwa (5.3018°N, 1.9930°W) in the forest ecological zone with perennial malaria transmission, and Accra (5.6037°N, 0.1870°W) and Cape Coast (5.1315°N, 1.2795°W) in the coastal savannah zone with perennial malaria transmission.
Map showing the 10 sentinel sites for surveillance of antimalarial drug resistance in Ghana. (Map is from https://binged.it/31uEjfO, and the URL of the original copyright holder [newafrica.com] is no longer active.)
Study samples.Archived filter paper blood blots collected from uncomplicated malaria patients aged ≤9 years were used. The samples were originally collected during the wet transmission seasons of five time periods (2007 to 2008, 2010 to 2011, 2012, 2013 to 2014, and 2015 to 2016) as determined by the availability of funds. The information on sample collection over the years was previously reported (21, 22, 31–36). The samples were collected before (day 0) and after treatment (with AL or AA) on days 2, 3, 7, 14, 21, and 28 per protocol with microscopically and PCR-detected infection. For this work, only day 0 samples were used. The samples (100 μl blood) were collected on Whatman 3 filter paper (Little Chalfont, UK), stored in plastic bags containing silica gels, and kept at room temperature until use. Mixed infections were not looked at, because one of the criteria for recruitment in our treatment efficacy study is to have solely P. falciparum parasites not in combination with the other species.
Molecular analysis for the detection of pfk13 SNPs in isolates.Malaria parasite DNA was extracted from 1,100 filter paper blood blots by using a QIAamp DNA minikit (Qiagen, Hilden, Germany) according to the manufacturer’s protocol. Nested PCR (nPCR) was performed to amplify the region of pfk13 gene using published protocols (9) with minimal modification. PCR products were sequenced using Sanger sequencing at Macrogen, Netherlands. For quality control purposes, 20 filter paper blood blots of known parasite sequences were sent to Venkatachalam Udhayakumar’s laboratory at the Malaria Centre, Centers for Disease Control and Prevention (CDC), Atlanta, GA, for analysis and confirmation of results.
pfk13 sequence analysis.Sequences were analyzed with the BLAST program (http://blast.ncbi.nlm.nih.gov/) to determine the authenticity of the sequences. Multiple sequences were aligned with MAFFT (EMBL.EBI, Hinxton, Cambridge, UK) using the 3D7 wild-type pfk13 sequence (PF3D7_1343700) as a reference. Consensus sequence editing and single nucleotide polymorphism (SNP) detection was carried out using the CLC Main Workbench 7.9.1 (Qiagen, Aarhus, Denmark) and the Benchling website (San Francisco, CA, USA). Frequencies of mutations were determined for each sample by individual counts. The chi-square test was used to determine the prevalence trends of the wild types over time. The CCF53_62 matrix analysis by Brick and Pizzi (26) was used to determine the effect of amino acid changes in proteins.
Ethics statement.The study protocol was approved by the NMIMR and Naval Medical Research Center’s Institutional Review Boards in compliance with all applicable federal regulations governing the protection of human subjects.
ACKNOWLEDGMENTS
The study was funded by the Armed Forces Health Surveillance Branch (AFHSB) and its Global Emerging Infections Surveillance and Response section (GEIS). The field work for the collection of the samples was partly funded by GEIS and the Global Fund to Fight AIDS, Tuberculosis and Malaria (GFATM)/National Malaria Control Programme (NMCP; Ghana). The funders played no role in the design of the study, sample collection, analysis, and interpretation of results as well as manuscript preparation and submission.
We thank Venkatachalam Udhayakumar and Naomi Lucchi of the Malaria Centre, CDC, Atlanta, GA, for the quality control analysis at their laboratory. We also thank the participants for the various studies under the Surveillance of Antimalarial Drug Resistance in Ghana. We declare no competing interests. The views expressed in this article are those of the authors and do not necessary reflect the official policy or position of the Department of the Navy, Department of Defense, nor the U.S. Government.
Andrew Letizia is a military Service member and Anne Fox is an employee of the U.S. government. This work was prepared as part of their official duties. Title 17, U.S.C., §105 provides that copyright protection under this title is not available for any work of the U.S. Government. Title 17, U.S.C., §101 defines a U.S. Government work as a work prepared by a military service member or employee of the U.S. Government as part of that person’s official duties.
FOOTNOTES
- Received 15 April 2019.
- Returned for modification 14 May 2019.
- Accepted 7 August 2019.
- Accepted manuscript posted online 19 August 2019.
This is a work of the U.S. Government and is not subject to copyright protection in the United States. Foreign copyrights may apply.