Spontaneous mutations in the Plasmodium falciparum sarcoplasmic/ endoplasmic reticulum Ca2+-ATPase (PfATP6) gene among geographically widespread parasite populations unexposed to artemisinin-based combination therapies

Geographic Plasmodium falciparum sarcoplasmic reticulum Ca2+ ATPase (PfSERCA) genotype diversity in India

Plasmodium falciparum sarcoplasmic reticulum Ca²⁺ ATPase (PfSERCA) is sarcoplasmic reticulum membrane bound transporter to regulate cytosol Ca²⁺ ions. Ca²⁺ act as secondary messenger and play important role in differentiation of parasite during its life cycle. Present study is epidemiological surveillance of PfSERCA (Pf3D7_0106300) gene fragment harboring 263, 402, 431 codon to look for its single nucleotide polymorphism which is well documented to be associated with Artemisinin tolerance. Filter paper with finger pricked blood samples for Plasmodium falciparum infected uncomplicated malaria patients were obtained for region as diverse as down the longitude from east to west of India i.e. Mizoram, Tripura, Meghalaya, Jharkhand, Odhisa. There observed no mutation for codon 263 at all study sites. Mizoram showed highest PfSERCA diversity with well known SNPs of L402 V, E431 K, A438 V and novel mutations as well i.e. A338 V, S357Y, S379Y. Tripura reported highest proportion of Plasmodium isolates (18.5%) with E431 K single nucleotide polymorphism. Moving towards the west i.e. Meghalaya, Jharkhand, Odhisa showed no occurrence of most prevalent PfSERCA 431, 402 polymorphism worldwide but some novel mutations and its haplotypes. In present study, significantly increased proportion of novel PfSERCA polymorphism among children suggests the susceptibility of these Plasmodium falciparum strains to acquired immunity. Mizoram, sharing open international border with south east asia, demonstrated highest PfSERCA diversity. Spatial PfSERCA diversity from far north east India to moving towards west implies its association with antimalarial susceptibility.

Molecular markers of anti-malarial drug resistance in Central, West and East African children with severe malaria

Background

The Plasmodium falciparum multidrug resistance 1 (PfMDR1), P. falciparum Ca²⁺-ATPase (PfATP6) and Kelch-13 propeller domain (PfK13) loci are molecular markers of parasite susceptibility to anti-malarial drugs. Their frequency distributions were determined in the isolates collected from children with severe malaria originating from three African countries.

Methods

Samples from 287 children with severe malaria [(Gabon: n = 114); (Ghana: n = 89); (Kenya: n = 84)] were genotyped for pfmdr1, pfatp6 and pfk13 loci by DNA sequencing and assessing pfmdr1 copy number variation (CNV) by real-time PCR.

Results

Pfmdr1-N86Y mutation was detected in 48, 10 and 10% in Lambaréné, Kumasi and Kisumu, respectively. At codon 184, the prevalence of the mutation was 73% in Lambaréné, 63% in Kumasi and 49% Kisumu. The S1034C and N1042D variants were absent at all three sites, while the frequency of the D1246Y mutation was 1, 3 and 13% in Lambaréné, Kumasi and Kisumu, respectively. Isolates with two pfmdr1 gene copy number predominantly harboured the N86Y wild-type allele and were mostly found in Kumasi (10%) (P < 0.0001). Among the main pfmdr1 haplotypes (NFD, NYD and YFD), NYD was associated with highest parasitaemia (P = 0.04). At the pfatp6 locus, H243Y and A623E mutations were observed at very low frequency at all three sites. The prevalence of the pfatp6 E431K variant was 6, 18 and 17% in Lambaréné, Kumasi and Kisumu, respectively. The L263E and S769N mutations were absent in all isolates. The pfk13 variants associated with artemisinin resistance in Southeast Asia were not observed. Eleven novel substitutions in the pfk13 locus occurring at low frequency were observed.

Conclusions

Artemisinins are still highly efficacious in large malaria-endemic regions though declining efficacy has occurred in Southeast Asia. The return of chloroquine-sensitive strains following the removal of drug pressure is observed. However, selection of wild-type alleles in the multidrug-resistance gene and the increased gene copy number is associated with reduced lumefantrine sensitivity. This study indicates a need to constantly monitor drug resistance to artemisinin in field isolates from malaria-endemic countries.

Molecular surveillance of Plasmodium falciparum drug resistance in the Republic of Congo: four and nine years after the introduction of artemisinin-based combination therapy

BACKGROUND

Resistance to anti-malarial drugs hinders efforts on malaria elimination and eradication. Following the global spread of chloroquine-resistant parasites, the Republic of Congo adopted artemisinin-based combination therapy (ACT) in 2006 as a first-line treatment for uncomplicated malaria. To assess the impacts after implementation of ACT, a molecular surveillance for anti-malarial drug resistance was conducted in Congo 4 and 9 years after the introduction of ACT.

METHODS

Blood samples of 431 febrile children aged 1-10 years were utilized from two previous studies conducted in 2010 (N = 311) and 2015 (N = 120). All samples were screened for malaria parasites using nested PCR. Direct sequencing was used to determine the frequency distribution of genetic variants in the anti-malarial drug-resistant Plasmodium falciparum genes (Pfcrt, Pfmdr1, Pfatp6, Pfk13) in malaria-positive isolates.

RESULTS

One-hundred and nineteen (N = 70 from 2010 and N = 49 from 2015) samples were positive for P. falciparum. A relative decrease in the proportion of chloroquine-resistant haplotype (CVIET) from 100% in 2005, 1 year before the introduction and implementation of ACT in 2006, to 98% in 2010 to 71% in 2015 was observed. Regarding the multidrug transporter gene, a considerable reduction in the frequency of the mutations N86Y (from 73 to 27%) and D1246Y (from 22 to 0%) was observed. However, the prevalence of the Y184F mutation remained stable (49% in 2010 compared to 54% in 2015). Isolates carrying the Pfatp6 H243Y was 25% in 2010 and this frequency was reduced to null in 2015. None of the parasites harboured the Pfk13 mutations associated with prolonged artemisinin clearance in Southeast Asia. Nevertheless, 13 new Pfk13 variants are reported among the investigated isolates.

CONCLUSION

The implementation of ACT has led to the decline in prevalence of chloroquine-resistant parasites in the Republic of Congo. However, the constant prevalence of the PfMDR1 Y184F mutation, associated with lumefantrine susceptibility, indicate a selective drug pressure still exists. Taken together, this study could serve as the basis for epidemiological studies monitoring the distribution of molecular markers of artemisinin resistance in the Republic of Congo.

Evaluation of artemether-lumefantrine efficacy in the treatment of uncomplicated malaria and its association with pfmdr1, pfatpase6 and K13-propeller polymorphisms in Luanda, Angola

Background

Drug resistance in Plasmodiumfalciparum has posed an obstacle to effective treatment and challenges many malaria control programmes in endemic areas. In Angola, until 2003, chloroquine (CQ) was used as first-line therapy for uncomplicated malaria. It was replaced initially by amodiaquine and, in 2006, by artemisinin-based combination therapy (ACT) with artemether-lumefantrine (AL, Coartem^®). Efficacy study of ACT, conducted in Angola between 2004 and 2005, showed a baseline efficacy of ≈99 %.

Methods

103 malaria patients were enrolled according to WHO proceedings. Patients were followed up with clinical and parasitological evaluations for 28 days, parasite density and identification was evaluated by microscopy, the pfmsp2 were genotyped by nested-PCR, to distinguish parasite recrudescence from new infections; the polymorphisms at codons 86 and 1246 of pfmdr1 gene, and 769 of pfatp6 gene were assessed by PCR–RFLP and sequencing for pfk13-propeller genotype.

Results

The cure rate was 91.3 %. The obtained results showed that from 103 patients, 12.6 % (n = 13) still had parasitaemia 1 day after the treatment was finished. On day 0, of the 94 evaluated samples, wild-type alleles were identified in 73.4 % (n = 69) for pfmdr1 N86Y position and only one sample carried the mutant allele (Y) for pfmdr1 1246; 14 % of these samples showed increased pfmdr1 copy number; 100 % (n = 21) had wild-type allele of k13 gene in all the studied positions.

Discussion

These results showed changes in parasite profile susceptibility to AL in comparison to the baseline data from 2002 to 2004 and on the genotyping characteristics; the clinical outcome after treatment with AL did not link a particular genotype with treatment failure; observed changes do not provide sufficient evidence for a treatment policy change, but they suggest that a carefully monitoring is needed in this area.

Ecotope-based entomological surveillance and molecular xenomonitoring of multidrug resistant malaria parasites in anopheles vectors

The emergence and spread of multidrug resistant (MDR) malaria caused by Plasmodium falciparum or Plasmodium vivax have become increasingly important in the Greater Mekong Subregion (GMS). MDR malaria is the heritable and hypermutable property of human malarial parasite populations that can decrease in vitro and in vivo susceptibility to proven antimalarial drugs as they exhibit dose-dependent drug resistance and delayed parasite clearance time in treated patients. MDR malaria risk situations reflect consequences of the national policy and strategy as this influences the ongoing national-level or subnational-level implementation of malaria control strategies in endemic GMS countries. Based on our experience along with current literature review, the design of ecotope-based entomological surveillance (EES) and molecular xenomonitoring of MDR falciparum and vivax malaria parasites in Anopheles vectors is proposed to monitor infection pockets in transmission control areas of forest and forest fringe-related malaria, so as to bridge malaria landscape ecology (ecotope and ecotone) and epidemiology. Malaria ecotope and ecotone are confined to a malaria transmission area geographically associated with the infestation of Anopheles vectors and particular environments to which human activities are related. This enables the EES to encompass mosquito collection and identification, salivary gland DNA extraction, Plasmodium- and species-specific identification, molecular marker-based PCR detection methods for putative drug resistance genes, and data management. The EES establishes strong evidence of Anopheles vectors carrying MDR P. vivax in infection pockets epidemiologically linked with other data obtained during which a course of follow-up treatment of the notified P. vivax patients receiving the first-line treatment was conducted. For regional and global perspectives, the EES would augment the epidemiological surveillance and monitoring of MDR falciparum and vivax malaria parasites in hotspots or suspected areas established in most endemic GMS countries implementing the National Malaria Control Programs, in addition to what is guided by the World Health Organization.

Plasmodium falciparum: genetic diversity and complexity of infections in an isolated village in western Thailand

Genetic diversity of Plasmodium falciparum is intimately associated with morbidity, mortality and malaria control strategies. It is therefore imperative to study genetic makeup and population structure of this parasite in endemic areas. In Kong Mong Tha, an isolated village in western Thailand, the majority of P. falciparum infections are asymptomatic. In this study we investigated complexity of infections and single nucleotide polymorphisms (SNPs) in the P. falciparum population of Kong Mong Tha, and compared results with those previously obtained from Mae Sod, in northwestern Thailand, where the majority of infections were symptomatic. Using PCR-based determination of the 5' merozoite surface protein 1 gene (msp1) recombinant types, we found that 39% of 59 P. falciparum isolates from Kong Mong Tha had multiple 5' recombinant types with a mean number of 1.54. These values were much lower than those obtained from Mae Sod: 96% for multiple infections and with a mean number of 3.61. Analysis of full-length sequences of two housekeeping genes, the P-type Ca(2+)-transporting ATPase gene (n=33) plus adenylosuccinate lyase gene (n=33), and three vaccine candidate antigen genes, msp1 (n=26), the circumsporozoite protein gene, csp (n=30) and the apical membrane antigen 1 gene, ama 1 (n=32), revealed that in all of these genes within-population SNP diversity was at similar levels between Kong Mong Tha and Mae Sod, suggesting that the extent of MOI and clinical manifestations of malaria are not strongly associated with genetic diversity. Additionally, we did not detect significant genetic differentiation between the two parasite populations, as estimated by the Wright's fixation index of inter-population variance in allele frequencies, suggesting that gene flow prevented the formation of population structuring. Thus, this study highlights unique features of P. falciparum populations in Thailand. The implications of these finding are discussed.

Plasmodium falciparum mitochondrial genetic diversity exhibits isolation-by-distance patterns supporting a sub-Saharan African origin

The geographical distribution of single nucleotide polymorphism (SNP) in the mitochondrial genome of the human malaria parasite Plasmodium falciparum was investigated. We identified 88 SNPs in 516 isolates from seven parasite populations in Africa, Southeast Asia and Oceania. Analysis of the SNPs postulated a sub-Saharan African origin and recovered a strong negative correlation between within-population SNP diversity and geographic distance from the putative African origin over Southeast Asia and Oceania. These results are consistent with those previously obtained for nuclear genome-encoded housekeeping genes, indicating that the pattern of inheritance does not substantially affect the geographical distribution of SNPs.

The interplay between drug resistance and fitness in malaria parasites

Controlling the spread of antimalarial drug resistance, especially resistance of Plasmodium falciparum to artemisinin-based combination therapies, is a high priority. Available data indicate that, as with other microorganisms, the spread of drug-resistant malaria parasites is limited by fitness costs that frequently accompany resistance. Resistance-mediating polymorphisms in malaria parasites have been identified in putative drug transporters and in target enzymes. The impacts of these polymorphisms on parasite fitness have been characterized in vitro and in animal models. Additional insights have come from analyses of samples from clinical studies, both evaluating parasites under different selective pressures and determining the clinical consequences of infection with different parasites. With some exceptions, resistance-mediating polymorphisms lead to malaria parasites that, compared with wild type, grow less well in culture and in animals, and are replaced by wild type when drug pressure diminishes in the clinical setting. In some cases, the fitness costs of resistance may be offset by compensatory mutations that increase virulence or changes that enhance malaria transmission. However, not enough is known about effects of resistance mediators on parasite fitness. A better appreciation of the costs of fitness-mediating mutations will facilitate the development of optimal guidelines for the treatment and prevention of malaria.

Antimalarial screening via large-scale purification of Plasmodium falciparum Ca2+-ATPase 6 and in vitro studies

The most severe form of human malaria is caused by the parasite Plasmodium falciparum. Despite the current need, there is no effective vaccine and parasites are becoming resistant to most of the antimalarials available. Therefore, there is an urgent need to discover new drugs from targets that have not yet suffered from drug pressure with the aim of overcoming the problem of new emerging resistance. Membrane transporters, such as P. falciparum Ca(2+)-ATPase 6 (PfATP6), the P. falciparum sarcoplasmic/endoplasmic reticulum Ca(2+)-ATPase (SERCA), have been proposed as potentially good antimalarial targets. The present investigation focuses on: (a) the large-scale purification of PfATP6 for maintenance of its enzymatic activity; (b) screening for PfATP6 inhibitors from a compound library; and (c) the selection of the best inhibitors for further tests on P. falciparum growth in vitro. We managed to heterologously express in yeast and purify an active form of PfATP6 as previously described, although in larger amounts. In addition to some classical SERCA inhibitors, a chemical library of 1680 molecules was screened. From these, we selected a pool of the 20 most potent inhibitors of PfATP6, presenting half maximal inhibitory concentration values in the range 1-9 μm. From these, eight were chosen for evaluation of their effect on P. falciparum growth in vitro, and the best compound presented a half maximal inhibitory concentration of ~ 2 μm. We verified the absence of an inhibitory effect of most of the compounds on mammalian SERCA1a, representing a potential advantage in terms of human toxicity. The present study describes a multidisciplinary approach allowing the selection of promising PfATP6-specific inhibitors with good antimalarial activity.

Genetic diversity and lack of artemisinin selection signature on the Plasmodium falciparum ATP6 in the Greater Mekong Subregion

The recent detection of clinical Artemisinin (ART) resistance manifested as delayed parasite clearance in the Cambodia-Thailand border area raises a serious concern. The mechanism of ART resistance is not clear; but the P. falciparum sarco/endoplasmic reticulum Ca²⁺-ATPase (PfSERCA or PfATP6) has been speculated to be the target of ARTs and thus a potential marker for ART resistance. Here we amplified and sequenced pfatp6 gene (∼3.6 Kb) in 213 samples collected after 2005 from the Greater Mekong Subregion, where ART drugs have been used extensively in the past. A total of 24 single nucleotide polymorphisms (SNPs), including 8 newly found in this study and 13 nonsynonymous, were identified. However, these mutations were either uncommon or also present in other geographical regions with limited ART use. None of the mutations were suggestive of directional selection by ARTs. We further analyzed pfatp6 from a worldwide collection of 862 P. falciparum isolates in 19 populations from Asia, Africa, South America and Oceania, which include samples from regions prior to and after deployments ART drugs. A total of 71 SNPs were identified, resulting in 106 nucleotide haplotypes. Similarly, many of the mutations were continent-specific and present at frequencies below 5%. The most predominant and perhaps the ancestral haplotype occurred in 441 samples and was present in 16 populations from Asia, Africa, and Oceania. The 3D7 haplotype found in 54 samples was the second most common haplotype and present in nine populations from all four continents. Assessment of the selection strength on pfatp6 in the 19 parasite populations found that pfatp6 in most of these populations was under purifying selection with an average d_N/d_S ratio of 0.333. Molecular evolution analyses did not detect significant departures from neutrality in pfatp6 for most populations, challenging the suitability of this gene as a marker for monitoring ART resistance.

Characterization of the Plasmodium falciparum sarcoplasmic/endoplasmic reticulum Ca2+-ATPase gene in samples from Equatorial Guinea before implementation of artemisinin-based combination therapy

Plasmodium falciparum resistance to the primary drugs used for treatment of malaria has become the main obstacle to malaria control. Artemisinin combination therapies are the current treatment strategy, and it has been suggested that resistance to artemisinin derivatives may be related to mutations in the Plasmodium falciparum sarcoplasmic-endoplasmic reticulum Ca(2+)-ATPase ortholog of the mammalian sarco-endoplasmic reticulum Ca(2+) ATPase gene, known as the pfatp6 gene. Thus, the purpose of this study was to determine the prevalence of single-nucleotide polymorphisms (SNPs) in pfatp6. The presence of different SNPs was detected by polymerase chain reaction amplification of the pfatp6 gene, and then sequencing to identify all possible alleles of the gene. A total of 20 SNPs were detected, including eight SNPs that have not been previously described: K481R in Malabo; R801H on Annobon Island; and the synonymous SNPs a141t, c1788t, a2211g, t2739g, a2760c, and g2836a. The genotypic profile of pfatp6 in samples from Equatorial Guinea, may be a useful epidemiologic tool for monitoring local efficacy of artemisinin combination therapies.

Molecular assessment of atpase6 mutations associated with artemisinin resistance among unexposed and exposed Plasmodium falciparum clinical isolates to artemisinin-based combination therapy

Background

Artemisinin-based combination therapy (ACT) is the mainstay of global efforts for treatment of Plasmodium falciparum malaria, but decline in its efficacy is the most important obstacle towards malaria control and elimination. Therefore, the present molecular analysis provides information on putative mutations associated with artemisinin resistance in P. falciparum clinical population unexposed and exposed to artesunate 4 years after adoption of ACT as the first-line anti-malarial therapy in Iran.

Methods

In this study, blood samples (n = 226) were collected from uncomplicated P. falciparum-infected patients from different health centers of Chabahar district in Sistan and Baluchistan province in the south-eastern part of Iran, during 2003 to 2010. All collected isolates were analysed for putative candidate mutations (TTA) L263E (GAA), (GAA) E431K (AAA), (GCA) A623E (GAA) and (AGT) S769N (AAT) of pfatpase6 gene using nested PCR/RFLP, followed by sequencing. Furthermore, the gene copy number was assessed by real-time quantitative PCR (RT-qPCR) in the presence of SYBR green.

Results

Neither the pfatpase6 L263E nor the A623E mutation was detected among all examined isolates. The E431K mutation was found in 23% of the analysed samples unexposed to ACT; however, it was detected in 17.8% (34/191) of P. falciparum isolates exposed to artesunate after 2007. High frequency of this single nucleotide polymorphisms (SNP) (overall 18.6%) among both examined groups (X² test, P>0.05) indicated that this SNP should be considered as an unrelated mutation to artemisinin resistance. In contrast, S769N mutation was not detected in unexposed isolates; however, it was found in 2.6% (5/191), four years after introduction of ACT in this malaria setting. Also, detected SNPs were not significantly frequent in both unexposed and exposed examined isolates (X² test, P> 0.05). Investigation in the copy number of pfatpase6 gene revealed a similar number of copy (n = 1) as in an isolate sensitive to artemisinin.

Conclusion

Taken together, the results suggest, in particular, that pfatpase6 S769N gene needs more consideration for its possible association with artesunate resistance among P. falciparum isolates.

Mechanisms of in vitro resistance to dihydroartemisinin in Plasmodium falciparum

The recent reports of artemisinin (ART) resistance in the Thai-Cambodian border area raise a serious concern on the long-term efficacy of ARTs. To elucidate the resistance mechanisms, we performed in vitro selection with dihydroartemisinin (DHA) and obtained two parasite clones from Dd2 with more than 25-fold decrease in susceptibility to DHA. The DHA-resistant clones were more tolerant of stressful growth conditions and more resistant to several commonly used antimalarial drugs than Dd2. The result is worrisome as many of the drugs are currently used as ART partners in malaria control. This study showed that the DHA resistance is not limited to ring stage, but also occurred in trophozoites and schizonts. Microarray and biochemical analyses revealed pfmdr1 amplification, elevation of the antioxidant defence network, and increased expression of many chaperones in the DHA-resistant parasites. Without drug pressure, the DHA-resistant parasites reverted to sensitivity in approximately 8 weeks, accompanied by de-amplification of pfmdr1 and reduced antioxidant activities. The parallel decrease and increase in pfmdr1 copy number and antioxidant activity and the up and down of DHA sensitivity strongly suggest that pfmdr1 and antioxidant defence play a role in in vitro resistance to DHA, providing potential molecular markers for ART resistance.

Questions over high frequency of mutant PfATP6 haplotypes in traveller isolates

A recent paper in Malaria Journal suggests that a high proportion of Plasmodium falciparum isolates found in travellers returning from a range of African countries carry the PfATP6 A623E S769N haplotype, and that this genotype is associated with artemether resistance. Such a finding would represent a substantial departure from the extensive literature reporting these individual mutations to be very rare, with the double mutation never documented. The number of isolates screened to obtain these double mutants is unstated, but highly relevant, not least because selection of isolates could have introduced significant confounders, such as timing of in vitro testing. An additional concern relates to the location of sequencing primers used to assess these positions. In the absence of clear information on these fundamental questions it would be appropriate to treat the findings with caution.

In vitro sensitivity of Plasmodium falciparum from China-Myanmar border area to major ACT drugs and polymorphisms in potential target genes

Drug resistance has always been one of the most important impediments to global malaria control. Artemisinin resistance has recently been confirmed in the Greater Mekong Subregion (GMS) and efforts for surveillance and containment are intensified. To determine potential mechanisms of artemisinin resistance and monitor the emergence and spread of resistance in other regions of the GMS, we investigated the in vitro sensitivity of 51 culture-adapted parasite isolates from the China-Myanmar border area to four drugs. The 50% inhibitory concentrations (IC₅₀s) of dihydroartemisinin, mefloquine and lumefantrine were clustered in a relatively narrow, 3- to 6-fold range, whereas the IC₅₀ range of artesunate was 12-fold. We assessed the polymorphisms of candidate resistance genes pfcrt, pfmdr1, pfATP6, pfmdr6 and pfMT (a putative metabolite/drug transporter). The K76T mutation in pfcrt reached fixation in the study parasite population, whereas point mutations in pfmdr1 and pfATP6 had low levels of prevalence. In addition, pfmdr1 gene amplification was not detected. None of the mutations in pfmdr1 and pfATP6 was associated significantly with in vitro sensitivity to artemisinin derivatives. The ABC transporter gene pfmdr6 harbored two point mutations, two indels, and number variations in three simple repeats. Only the length variation in a microsatellite repeat appeared associated with altered sensitivity to dihydroartemisinin. The PfMT gene had two point mutations and one codon deletion; the I30N and N496- both reached high levels of prevalence. However, none of the SNPs or haplotypes in PfMT were correlated significantly with resistance to the four tested drugs. Compared with other parasite populations from the GMS, our studies revealed drastically different genotype and drug sensitivity profiles in parasites from the China-Myanmar border area, where artemisinins have been deployed extensively for over 30 years.

Artemether resistance in vitro is linked to mutations in PfATP6 that also interact with mutations in PfMDR1 in travellers returning with Plasmodium falciparum infections

Background

Monitoring resistance phenotypes for Plasmodium falciparum, using in vitro growth assays, and relating findings to parasite genotype has proved particularly challenging for the study of resistance to artemisinins.

Methods

Plasmodium falciparum isolates cultured from 28 returning travellers diagnosed with malaria were assessed for sensitivity to artemisinin, artemether, dihydroartemisinin and artesunate and findings related to mutations in pfatp6 and pfmdr1.

Results

Resistance to artemether in vitro was significantly associated with a pfatp6 haplotype encoding two amino acid substitutions (pfatp6 A623E and S769N; (mean IC₅₀ (95% CI) values of 8.2 (5.7 – 10.7) for A623/S769 versus 623E/769 N 13.5 (9.8 – 17.3) nM with a mean increase of 65%; p = 0.012). Increased copy number of pfmdr1 was not itself associated with increased IC₅₀ values for artemether, but when interactions between the pfatp6 haplotype and increased copy number of pfmdr1 were examined together, a highly significant association was noted with IC₅₀ values for artemether (mean IC₅₀ (95% CI) values of 8.7 (5.9 – 11.6) versus 16.3 (10.7 – 21.8) nM with a mean increase of 87%; p = 0.0068). Previously described SNPs in pfmdr1 are also associated with differences in sensitivity to some artemisinins.

Conclusions

These findings were further explored in molecular modelling experiments that suggest mutations in pfatp6 are unlikely to affect differential binding of artemisinins at their proposed site, whereas there may be differences in such binding associated with mutations in pfmdr1. Implications for a hypothesis that artemisinin resistance may be exacerbated by interactions between PfATP6 and PfMDR1 and for epidemiological studies to monitor emerging resistance are discussed.

Lack of association of the S769N mutation in Plasmodium falciparum SERCA (PfATP6) with resistance to artemisinins

The recent emergence of artemisinin (ART) resistance in Plasmodium falciparum in western Cambodia, manifested as delayed parasite clearance, is a big threat to the long-term efficacy of this family of antimalarial drugs. Among the multiple candidate genes associated with ART resistance in P. falciparum, the sarcoplasmic/endoplasmic reticulum Ca(2+)-ATPase PfATP6 has been postulated as a specific target of ARTs. The PfATP6 gene harbors multiple single-nucleotide polymorphisms in field parasite populations, and S769N has been associated with decreased sensitivity to artemether in parasite populations from French Guiana. In this study, we used an allelic exchange strategy to engineer parasite lines carrying the S769N mutations in P. falciparum strain 3D7 and evaluated whether introduction of this mutation modulated parasite sensitivity to ART derivatives. Using three transgenic lines carrying the 769N mutation and two transgenic lines carrying the wild-type 769S as controls, we found that S769N did not affect PfATP6 gene expression. We compared the sensitivities of these parasite lines to three ART derivatives, artemether, artesunate, and dihydroartemisinin, in 18 biological experiments and detected no significant effect of the S769N mutation on parasite response to these ART derivatives. This study provides further evidence for the lack of association of PfATP6 with ART resistance.

Protein-based signatures of functional evolution in Plasmodium falciparum

Background

It has been known for over a decade that Plasmodium falciparum proteins are enriched in non-globular domains of unknown function. The potential for these regions of protein sequence to undergo high levels of genetic drift provides a fundamental challenge to attempts to identify the molecular basis of adaptive change in malaria parasites.

Results

Evolutionary comparisons were undertaken using a set of forty P. falciparum metabolic enzyme genes, both within the hominid malaria clade (P. reichenowi) and across the genus (P. chabaudi). All genes contained coding elements highly conserved across the genus, but there were also a large number of regions of weakly or non-aligning coding sequence. These displayed remarkable levels of non-synonymous fixed differences within the hominid malaria clade indicating near complete release from purifying selection (dN/dS ratio at residues non-aligning across genus: 0.64, dN/dS ratio at residues identical across genus: 0.03). Regions of low conservation also possessed high levels of hydrophilicity, a marker of non-globularity. The propensity for such regions to act as potent sources of non-synonymous genetic drift within extant P. falciparum isolates was confirmed at chromosomal regions containing genes known to mediate drug resistance in field isolates, where 150 of 153 amino acid variants were located in poorly conserved regions. In contrast, all 22 amino acid variants associated with drug resistance were restricted to highly conserved regions. Additional mutations associated with laboratory-selected drug resistance, such as those in PfATPase4 selected by spiroindolone, were similarly restricted while mutations in another calcium ATPase (PfSERCA, a gene proposed to mediate artemisinin resistance) that reach significant frequencies in field isolates were located exclusively in poorly conserved regions consistent with genetic drift.

Conclusion

Coding sequences of malaria parasites contain prospectively definable domains subject to neutral or nearly neutral evolution on a scale that appears unrivalled in biology. This distinct evolutionary landscape has potential to confound analytical methods developed for other genera. Against this tide of genetic drift, polymorphisms mediating functional change stand out to such an extent that evolutionary context provides a useful signal for identifying the molecular basis of drug resistance in malaria parasites, a finding that is of relevance to both genome-wide and candidate gene studies in this genus.