Origin of robustness in generating drug-resistant malaria parasites

Evolutionary biology of antimalarial drug resistance: Understanding of the evolutionary dynamics

Anti-malarial drug resistance poses a significant challenge to global malaria control efforts, necessitating a deeper understanding of the evolutionary dynamics underlying the emergence and spread of resistance. This study explores how evolutionary theory provides a framework for elucidating the molecular mechanisms and genetic variation within parasite populations that drive resistance evolution. Drawing on recent research findings, we discuss the role of natural selection, genetic diversity, and fitness costs in shaping resistance phenotypes. Additionally, we highlight the implications of evolutionary insights for antimalarial drug policy, treatment guidelines, and future research directions. By integrating evolutionary biology principles with molecular epidemiology, this review aims to inform strategies for combating antimalarial drug resistance and advancing malaria treatment efforts. Using evolutionary theory to understand the dynamics of antimalarial drug resistance at the molecular level, we explored the influence of genetic variation within parasite populations on the emergence and spread of resistance. Antimalarial drug resistance poses a formidable challenge to global malaria control. By applying evolutionary theory to understand the dynamics of resistance emergence and spread at the molecular level, researchers can develop more effective strategies for surveillance, prevention, and treatment of drug-resistant malaria.

Assessing fitness costs in malaria parasites: a comprehensive review and implications for drug resistance management

Artemisinin-based combination therapy (ACT) remains a broadly effective anti-malarial drug combination, but the emergence of resistance is threatening its effectiveness. Limiting the spread of these drug-resistant parasites and delaying the emergence of resistance in new areas are of high priority. Understanding the evolution of resistance relies on discerning the fitness costs and benefits associated with resistance mutations. If the cost associated with resistance in an untreated host is sufficiently large relative to the benefit of resistance in a treated host, then the spread of resistance can be mitigated by ensuring sufficient hosts free from that active pharmaceutical ingredient. There is no straightforward way to measure these fitness costs, and each approach that has been used has its limitations. Here, the evidence of fitness costs as measured using field data, animal models, and in vitro models is reviewed for three of the main current or past first-line treatments for malaria: chloroquine (CQ), sulfadoxine-pyrimethamine (SP), and artemisinin derivatives (ART). Despite the difficulties of assessing fitness costs, there is a good amount of evidence of fitness costs in drug-resistant Plasmodium falciparum parasites. The most persuasive evidence comes from resistance reversal observed following the cessation of the use of chloroquine. Comparable evidence cannot be obtained for SP- and ART-resistant parasites, due to the absence of complete cessation of these drugs in the field. Data from in vitro and animal models are variable. While fitness costs are often observed, their presence is not universal across all resistant strains. The extent and nature of these fitness costs can vary greatly depending on the specific genetic factors involved and the ecological context in which the parasites evolve. As a result, it is essential to avoid making broad generalizations about the prevalence or impact of fitness costs in drug-resistant malaria parasites. Focusing on fitness costs as a vulnerability in resistant parasites can guide their evolutionary trajectory towards minimizing their fitness. By accurately predicting these costs, efforts to extend the effectiveness of anti-malarials can be enhanced, limiting resistance evolution and advancing malaria control and elimination goals.

folA thyA knockout E. coli as a suitable surrogate model for evaluation of antifolate sensitivity against PfDHFR-TS

A new superior bacteria complementation model was achieved for testing antifolate compounds and investigating antifolate resistance in the dihydrofolate reductase (DHFR) enzyme of the malaria parasite. Earlier models depended on the addition of trimethoprim (TMP) to chemically suppress the host Escherichia coli (Ec) DHFR function. However, incomplete suppression of EcDHFR and potential interference of antibiotics needed to maintain plasmids for complementary gene expression can complicate the interpretations. To overcome such limitations, the folA (F) and thyA (T) genes were genetically knocked out (Δ) in E. coli BL21(DE3). The resulting EcΔFΔT cells were thymidine auxotroph where thymidine supplementation or functional complementation with heterologous DHFR-thymidylate synthase (TS) is needed to restore the loss of gene functions. When tested against pyrimethamine (PYR) and its analogs designed to target Plasmodium falciparum (Pf) DHFR-TS, the 50 % inhibitory concentration values obtained from EcΔFΔT surrogates expressing wildtype (PfTM4) or double mutant (PfK1) DHFR-TS showed strong correlations to the results obtained from the standard in vitro P. falciparum growth inhibition assay. Interestingly, while TMP had little effect on the susceptibility to PYR and analogs in EcΔFΔT expressing PfDHFR-TS, it hypersensitized the chemically knockdown E. coli BL21(DE3) expressing PfTM4 DHFR-TS but desensitized the one carrying PfK1 DHFR-TS. The low intrinsic expression level of PfTM4 in E. coli BL21(DE3) by western blot analysis may explain the hypersensitive to antifolates of chemical knockdown bacteria surrogate. These results demonstrated the usefulness of EcΔFΔT surrogate as a new tool for antifolate antimalarial screening with potential application for investigation of antifolate resistance mechanism.

Strong positive selection biases identity-by-descent-based inferences of recent demography and population structure in Plasmodium falciparum

Malaria genomic surveillance often estimates parasite genetic relatedness using metrics such as Identity-By-Decent (IBD), yet strong positive selection stemming from antimalarial drug resistance or other interventions may bias IBD-based estimates. In this study, we use simulations, a true IBD inference algorithm, and empirical data sets from different malaria transmission settings to investigate the extent of this bias and explore potential correction strategies. We analyze whole genome sequence data generated from 640 new and 3089 publicly available Plasmodium falciparum clinical isolates. We demonstrate that positive selection distorts IBD distributions, leading to underestimated effective population size and blurred population structure. Additionally, we discover that the removal of IBD peak regions partially restores the accuracy of IBD-based inferences, with this effect contingent on the population's background genetic relatedness and extent of inbreeding. Consequently, we advocate for selection correction for parasite populations undergoing strong, recent positive selection, particularly in high malaria transmission settings.

Direct long read visualization reveals metabolic interplay between two antimalarial drug targets

Increases in the copy number of large genomic regions, termed genome amplification, are an important adaptive strategy for malaria parasites. Numerous amplifications across the Plasmodium falciparum genome contribute directly to drug resistance or impact the fitness of this protozoan parasite. During the characterization of parasite lines with amplifications of the dihydroorotate dehydrogenase (DHODH) gene, we detected increased copies of an additional genomic region that encompassed 3 genes (~5 kb) including GTP cyclohydrolase I (GCH1 amplicon). While this gene is reported to increase the fitness of antifolate resistant parasites, GCH1 amplicons had not previously been implicated in any other antimalarial resistance context. Here, we further explored the association between GCH1 and DHODH copy number. Using long read sequencing and single read visualization, we directly observed a higher number of tandem GCH1 amplicons in parasites with increased DHODH copies (up to 9 amplicons) compared to parental parasites (3 amplicons). While all GCH1 amplicons shared a consistent structure, expansions arose in 2-unit steps (from 3 to 5 to 7, etc copies). Adaptive evolution of DHODH and GCH1 loci was further bolstered when we evaluated prior selection experiments; DHODH amplification was only successful in parasite lines with pre-existing GCH1 amplicons. These observations, combined with the direct connection between metabolic pathways that contain these enzymes, lead us to propose that the GCH1 locus is beneficial for the fitness of parasites exposed to DHODH inhibitors. This finding highlights the importance of studying variation within individual parasite genomes as well as biochemical connections of drug targets as novel antimalarials move towards clinical approval.

Strong Positive Selection Biases Identity-By-Descent-Based Inferences of Recent Demography and Population Structure in Plasmodium falciparum

Malaria genomic surveillance often estimates parasite genetic relatedness using metrics such as Identity-By-Decent (IBD). Yet, strong positive selection stemming from antimalarial drug resistance or other interventions may bias IBD-based estimates. In this study, we utilized simulations, a true IBD inference algorithm, and empirical datasets from different malaria transmission settings to investigate the extent of such bias and explore potential correction strategies. We analyzed whole genome sequence data generated from 640 new and 4,026 publicly available Plasmodium falciparum clinical isolates. Our findings demonstrated that positive selection distorts IBD distributions, leading to underestimated effective population size and blurred population structure. Additionally, we discovered that the removal of IBD peak regions partially restored the accuracy of IBD-based inferences, with this effect contingent on the population's background genetic relatedness. Consequently, we advocate for selection correction for parasite populations undergoing strong, recent positive selection, particularly in high malaria transmission settings.

Prevalence of malaria resistance-associated mutations in Plasmodium falciparum circulating in 2017–2018, Bo, Sierra Leone

Introduction

In spite of promising medical, sociological, and engineering strategies and interventions to reduce the burden of disease, malaria remains a source of significant morbidity and mortality, especially among children in sub-Saharan Africa. In particular, progress in the development and administration of chemotherapeutic agents is threatened by evolved resistance to most of the antimalarials currently in use, including artemisinins.

Methods

This study analyzed the prevalence of mutations associated with antimalarial resistance in Plasmodium falciparum from 95 clinical samples collected from individuals with clinically confirmed malaria at a hospital in Bo, Sierra Leone between May 2017 and December 2018. The combination of polymerase chain reaction amplification and subsequent high throughput DNA sequencing was used to determine the presence of resistance-associated mutations in five P. falciparum genes – pfcrt, pfmdr1, pfdhfr, pfdhps and pfkelch13. The geographic origin of parasites was assigned using mitochondrial sequences.

Results

Relevant mutations were detected in the pfcrt (22%), pfmdr1 (>58%), pfdhfr (100%) and pfdhps (>80%) genes while no resistance-associated mutations were found in the pfkelch13 gene. The mitochondrial barcodes were consistent with a West African parasite origin with one exception indicating an isolate imported from East Africa.

Discussion

Detection of the pfmdr1 NFSND haplotype in 50% of the samples indicated the increasing prevalence of strains with elevated tolerance to artemeter + lumefantrine (AL) threatening the combination currently used to treat uncomplicated malaria in Sierra Leone. The frequency of mutations linked to resistance to antifolates suggests widespread resistance to the drug combination used for intermittent preventive treatment during pregnancy.

Prevalence of Mutations in the Pfdhfr, Pfdhps, and Pfmdr1 Genes of Malarial Parasites Isolated from Symptomatic Patients in Dogondoutchi, Niger

The effectiveness of artemisinin-based combination therapies (ACTs) depends not only on that of artemisinin but also on that of partner molecules. This study aims to evaluate the prevalence of mutations in the Pfdhfr, Pfdhps, and Pfmdr1 genes from isolates collected during a clinical study. Plasmodium genomic DNA samples extracted from symptomatic malaria patients from Dogondoutchi, Niger, were sequenced by the Sanger method to determine mutations in the Pfdhfr (codons 51, 59, 108, and 164), Pfdhps (codons 436, 437, 540, 581, and 613), and Pfmdr1 (codons 86, 184, 1034, and 1246) genes. One hundred fifty-five (155) pre-treatment samples were sequenced for the Pfdhfr, Pfdhps, and Pfmdr1 genes. A high prevalence of mutations in the Pfdhfr gene was observed at the level of the N51I (84.97%), C59R (92.62%), and S108N (97.39%) codons. The key K540E mutation in the Pfdhps gene was not observed. Only one isolate was found to harbor a mutation at codon I431V. The most common mutation on the Pfmdr1 gene was Y184F in 71.43% of the mutations found, followed by N86Y in 10.20%. The triple-mutant haplotype N51I/C59R/S108N (IRN) was detected in 97% of the samples. Single-mutant (ICS and NCN) and double-mutant (IRS, NRN, and ICN) haplotypes were prevalent at 97% and 95%, respectively. Double-mutant haplotypes of the Pfdhps (581 and 613) and Pfmdr (86 and 184) were found in 3% and 25.45% of the isolates studied, respectively. The study focused on the molecular analysis of the sequencing of the Pfdhfr, Pfdhps, and Pfmdr1 genes. Although a high prevalence of mutations in the Pfdhfr gene have been observed, there is a lack of sulfadoxine pyrimethamine resistance. There is a high prevalence of mutation in the Pfmdr184 codon associated with resistance to amodiaquine. These data will be used by Niger’s National Malaria Control Program to better monitor the resistance of Plasmodium to partner molecules in artemisinin-based combination therapies.

Nutrient Limitation Magnifies Fitness Costs of Antimalarial Drug Resistance Mutations

Drug resistance mutations tend to disrupt key physiological processes and frequently carry fitness costs, which are a central determinant of the rate of spread of these mutations in natural populations. Head-to-head competition assays provide a standard approach to measuring fitness for malaria parasites. These assays typically use a standardized culture medium containing RPMI 1640, which has a 1.4- to 5.5-fold higher concentration of amino acids than human blood. In this rich medium, we predict that fitness costs will be underestimated because resource competition is weak. We tested this prediction using an artemisinin-sensitive parasite edited to contain kelch-C580Y or R561H mutations conferring resistance to artemisinin or synonymous control mutations. We examined the impact of these single amino acid mutations on fitness, using replicated head-to-head competition experiments conducted in media containing (i) normal RPMI, (ii) modified RPMI with reduced amino acid concentration, (iii) RPMI containing only isoleucine, or (iv) 3-fold diluted RPMI. We found a significant 1.3- to 1.4-fold increase in fitness costs measured in modified and isoleucine-only media relative to normal media, while fitness costs were 2.5-fold higher in diluted media. We conclude that fitness costs are strongly affected by media composition and will be significantly underestimated in normal RPMI. Several components differed between media, including pABA and sodium bicarbonate concentrations, so we cannot directly determine which is responsible. Elevated fitness costs in nature will limit spread of artemisinin (ART) resistance but will also promote evolution of compensatory mutations that restore fitness and can be exploited to maximize selection in laboratory experiments.

Evidence for the early emergence of piperaquine-resistant Plasmodium falciparum malaria and modeling strategies to mitigate resistance

Multidrug-resistant Plasmodium falciparum parasites have emerged in Cambodia and neighboring countries in Southeast Asia, compromising the efficacy of first-line antimalarial combinations. Dihydroartemisinin + piperaquine (PPQ) treatment failure rates have risen to as high as 50% in some areas in this region. For PPQ, resistance is driven primarily by a series of mutant alleles of the P. falciparum chloroquine resistance transporter (PfCRT). PPQ resistance was reported in China three decades earlier, but the molecular driver remained unknown. Herein, we identify a PPQ-resistant pfcrt allele (China C) from Yunnan Province, China, whose genotypic lineage is distinct from the PPQ-resistant pfcrt alleles currently observed in Cambodia. Combining gene editing and competitive growth assays, we report that PfCRT China C confers moderate PPQ resistance while re-sensitizing parasites to chloroquine (CQ) and incurring a fitness cost that manifests as a reduced rate of parasite growth. PPQ transport assays using purified PfCRT isoforms, combined with molecular dynamics simulations, highlight differences in drug transport kinetics and in this transporter’s central cavity conformation between China C and the current Southeast Asian PPQ-resistant isoforms. We also report a novel computational model that incorporates empirically determined fitness landscapes at varying drug concentrations, combined with antimalarial susceptibility profiles, mutation rates, and drug pharmacokinetics. Our simulations with PPQ-resistant or -sensitive parasite lines predict that a three-day regimen of PPQ combined with CQ can effectively clear infections and prevent the evolution of PfCRT variants. This work suggests that including CQ in combination therapies could be effective in suppressing the evolution of PfCRT-mediated multidrug resistance in regions where PPQ has lost efficacy.

The recent emergence of Plasmodium falciparum parasite resistance to the antimalarial drug piperaquine (PPQ) has contributed to frequent treatment failures across Southeast Asia, originating in Cambodia. Here, we show that earlier reports of PPQ resistance in Yunnan Province, China could be explained by the unique China C variant of the P. falciparum chloroquine resistance transporter PfCRT. Gene-edited parasites show a loss of fitness and parasite resensitization to the chemically related former first-line antimalarial chloroquine, while acquiring PPQ resistance via drug efflux. Molecular features of drug resistance were examined using biochemical assays to measure mutant PfCRT-mediated drug transport and molecular dynamics simulations with the recently solved PfCRT structure to assess changes in the central drug-binding cavity. We also describe a new computational model that incorporates parasite mutation rates, fitness costs, antimalarial susceptibilities, and drug pharmacological profiles to predict how infections with parasite strains expressing distinct PfCRT variants can evolve and be selected in response to different drug pressures and regimens. Simulations predict that a three-day regimen of PPQ plus chloroquine would be fully effective at preventing recrudescence of drug-resistant infections.

GTP Cyclohydrolase I as a Potential Drug Target: New Insights into Its Allosteric Modulation via Normal Mode Analysis

Guanosine triphosphate (GTP) cyclohydrolase I (GCH1) catalyzes the conversion of GTP into dihydroneopterin triphosphate (DHNP). DHNP is the first intermediate of the folate de novo biosynthesis pathway in prokaryotic and lower eukaryotic microorganisms and the tetrahydrobiopterin (BH4) biosynthesis pathway in higher eukaryotes. The de novo folate biosynthesis provides essential cofactors for DNA replication, cell division, and synthesis of key amino acids in rapidly replicating pathogen cells, such as Plasmodium falciparum (P. falciparum), a causative agent of malaria. In eukaryotes, the product of the BH4 biosynthesis pathway is essential for the production of nitric oxide and several neurotransmitter precursors. An increased copy number of the malaria parasite P. falciparum GCH1 gene has been reported to influence antimalarial antifolate drug resistance evolution, whereas mutations in the human GCH1 are associated with neuropathic and inflammatory pain disorders. Thus, GCH1 stands as an important and attractive drug target for developing therapeutics. The GCH1 intrinsic dynamics that modulate its activity remains unclear, and key sites that exert allosteric effects across the structure are yet to be elucidated. This study employed the anisotropic network model to analyze the intrinsic motions of the GCH1 structure alone and in complex with its regulatory partner protein. We showed that the GCH1 tunnel-gating mechanism is regulated by a global shear motion and an outward expansion of the central five-helix bundle. We further identified hotspot residues within sites of structural significance for the GCH1 intrinsic allosteric modulation. The obtained results can provide a solid starting point to design novel antineuropathic treatments for humans and novel antimalarial drugs against the malaria parasite P. falciparum GCH1 enzyme.

Decreased Susceptibility to Dihydrofolate Reductase Inhibitors Associated With Genetic Polymorphisms in Ugandan Plasmodium falciparum Isolates

BACKGROUND

The Plasmodium falciparum dihydrofolate reductase (PfDHFR) inhibitors pyrimethamine and cycloguanil (the active metabolite of proguanil) have important roles in malaria chemoprevention, but drug resistance challenges their efficacies. A new compound, P218, was designed to overcome resistance, but drug-susceptibility data for P falciparum field isolates are limited.

METHODS

We studied ex vivo PfDHFR inhibitor susceptibilities of 559 isolates from Tororo and Busia districts, Uganda, from 2016 to 2020, sequenced 383 isolates, and assessed associations between genotypes and drug-susceptibility phenotypes.

RESULTS

Median half-maximal inhibitory concentrations (IC50s) were 42 100 nM for pyrimethamine, 1200 nM for cycloguanil, 13000 nM for proguanil, and 0.6 nM for P218. Among sequenced isolates, 3 PfDHFR mutations, 51I (100%), 59R (93.7%), and 108N (100%), were very common, as previously seen in Uganda, and another mutation, 164L (12.8%), had moderate prevalence. Increasing numbers of mutations were associated with decreasing susceptibility to pyrimethamine, cycloguanil, and P218, but not proguanil, which does not act directly against PfDHFR. Differences in P218 susceptibilities were modest, with median IC50s of 1.4 nM for parasites with mixed genotype at position 164 and 5.7 nM for pure quadruple mutant (51I/59R/108N/164L) parasites.

CONCLUSIONS

Resistance-mediating PfDHFR mutations were common in Ugandan isolates, but P218 retained excellent activity against mutant parasites.

Chloroquine and Sulfadoxine-Pyrimethamine Resistance in Sub-Saharan Africa-A Review

Malaria is a great concern for global health and accounts for a large amount of morbidity and mortality, particularly in Africa, with sub-Saharan Africa carrying the greatest burden of the disease. Malaria control tools such as insecticide-treated bed nets, indoor residual spraying, and antimalarial drugs have been relatively successful in reducing the burden of malaria; however, sub-Saharan African countries encounter great challenges, the greatest being antimalarial drug resistance. Chloroquine (CQ) was the first-line drug in the 20th century until it was replaced by sulfadoxine-pyrimethamine (SP) as a consequence of resistance. The extensive use of these antimalarials intensified the spread of resistance throughout sub-Saharan Africa, thus resulting in a loss of efficacy for the treatment of malaria. SP was replaced by artemisinin-based combination therapy (ACT) after the emergence of resistance toward SP; however, the use of ACTs is now threatened by the emergence of resistant parasites. The decreased selective pressure on CQ and SP allowed for the reintroduction of sensitivity toward those antimalarials in regions of sub-Saharan Africa where they were not the primary drug for treatment. Therefore, the emergence and spread of antimalarial drug resistance should be tracked to prevent further spread of the resistant parasites, and the re-emergence of sensitivity should be monitored to detect the possible reappearance of sensitivity in sub-Saharan Africa.

Changing Prevalence of Potential Mediators of Aminoquinoline, Antifolate, and Artemisinin Resistance Across Uganda

BACKGROUND

In Uganda, artemether-lumefantrine is recommended for malaria treatment and sulfadoxine-pyrimethamine for chemoprevention during pregnancy, but drug resistance may limit efficacies.

METHODS

Genetic polymorphisms associated with sensitivities to key drugs were characterized in samples collected from 16 sites across Uganda in 2018 and 2019 by ligase detection reaction fluorescent microsphere, molecular inversion probe, dideoxy sequencing, and quantitative polymerase chain reaction assays.

RESULTS

Considering transporter polymorphisms associated with resistance to aminoquinolines, the prevalence of Plasmodium falciparum chloroquine resistance transporter (PfCRT) 76T decreased, but varied markedly between sites (0-46% in 2018; 0-23% in 2019); additional PfCRT polymorphisms and plasmepsin-2/3 amplifications associated elsewhere with resistance to piperaquine were not seen. For P. falciparum multidrug resistance protein 1, in 2019 the 86Y mutation was absent at all sites, the 1246Y mutation had prevalence ≤20% at 14 of 16 sites, and gene amplification was not seen. Considering mutations associated with high-level sulfadoxine-pyrimethamine resistance, prevalences of P. falciparum dihydrofolate reductase 164L (up to 80%) and dihydropteroate synthase 581G (up to 67%) were high at multiple sites. Considering P. falciparum kelch protein propeller domain mutations associated with artemisinin delayed clearance, prevalence of the 469Y and 675V mutations has increased at multiple sites in northern Uganda (up to 23% and 41%, respectively).

CONCLUSIONS

We demonstrate concerning spread of mutations that may limit efficacies of key antimalarial drugs.

New glycine derived peptides bearing benzenesulphonamide as an antiplasmodial agent

In the tropics, malaria is among the most serious infectious diseases in developing countries. The discovery of the artemesinin antimalarial drug not too long ago was a major breakthrough in the effort to combat the malaria disease.

Genetic diversity of the Plasmodium falciparum GTP-cyclohydrolase 1, dihydrofolate reductase and dihydropteroate synthetase genes reveals new insights into sulfadoxine-pyrimethamine antimalarial drug resistance

Plasmodium falciparum parasites resistant to antimalarial treatments have hindered malaria disease control. Sulfadoxine-pyrimethamine (SP) was used globally as a first-line treatment for malaria after wide-spread resistance to chloroquine emerged and, although replaced by artemisinin combinations, is currently used as intermittent preventive treatment of malaria in pregnancy and in young children as part of seasonal malaria chemoprophylaxis in sub-Saharan Africa. The emergence of SP-resistant parasites has been predominantly driven by cumulative build-up of mutations in the dihydrofolate reductase (pfdhfr) and dihydropteroate synthetase (pfdhps) genes, but additional amplifications in the folate pathway rate-limiting pfgch1 gene and promoter, have recently been described. However, the genetic make-up and prevalence of those amplifications is not fully understood. We analyse the whole genome sequence data of 4,134 P. falciparum isolates across 29 malaria endemic countries, and reveal that the pfgch1 gene and promoter amplifications have at least ten different forms, occurring collectively in 23% and 34% in Southeast Asian and African isolates, respectively. Amplifications are more likely to be present in isolates with a greater accumulation of pfdhfr and pfdhps substitutions (median of 1 additional mutations; P<0.00001), and there was evidence that the frequency of pfgch1 variants may be increasing in some African populations, presumably under the pressure of SP for chemoprophylaxis and anti-folate containing antibiotics used for the treatment of bacterial infections. The selection of P. falciparum with pfgch1 amplifications may enhance the fitness of parasites with pfdhfr and pfdhps substitutions, potentially threatening the efficacy of this regimen for prevention of malaria in vulnerable groups. Our work describes new pfgch1 amplifications that can be used to inform the surveillance of SP drug resistance, its prophylactic use, and future experimental work to understand functional mechanisms.

Pyrimethamine Elicits Antitumor Effects on Prostate Cancer by Inhibiting the p38-NF-κB Pathway

Since incurable castration-resistant prostate cancer (CRPC) inevitably develops following treatment with androgen deprivation therapy, there is an urgent need to devise new therapeutic strategies to treat this cancer. Pyrimethamine, an FDA-approved antimalarial drug, is known to exert an antitumor activity in various types of human cancer cells. However, whether pyrimethamine can inhibit prostate cancer is not well established. Hence, the present study aimed to characterize the mechanism of action of pyrimethamine on prostate cancer. We investigated the potential effect of pyrimethamine on cell proliferation, cell cycle, and apoptosis in metastatic DU145 and PC3 prostate cancer cells. We found that pyrimethamine inhibited cell proliferation, induced cell cycle arrest in the S phase, and promoted cell apoptosis of prostate cells in vitro; it also suppressed tumor growth in xenograft models. In addition, we observed that pyrimethamine suppressed prostate cancer growth by inhibiting the p38-NF-κB axis in vitro and in vivo. Thus, this study demonstrates that pyrimethamine is a novel p38 inhibitor that can exert antiproliferative and proapoptotic effects in prostate cancer by affecting cell cycle and intrinsic apoptotic signaling, thereby providing a novel strategy for using pyrimethamine in CRPC treatment.

Role of a Concentration Gradient in Malaria Drug Resistance Evolution: A Combined within- and between-Hosts Modelling Approach

Resistance to antimalarial drugs is currently a growing public health problem, resulting in more cases with treatment failure. Although previous studies suggested that a concentration gradient facilitates the antibiotic resistance evolution in bacteria, no attempt has been made to investigate the roles of a concentration gradient in malaria drug resistance. Unlike the person-to-person mode of transmission of bacteria, the malaria parasites need to switch back and forth between the human and mosquito hosts to complete the life cycle and to spread the resistant alleles. Here we developed a stochastic combined within- and between-hosts evolutionary dynamics model specific to malaria parasites in order to investigate the influence of an antimalarial concentration gradient on the evolutionary dynamics of malaria drug resistance. Every stage of malaria development in both human and mosquito hosts are individually modelled using the tau-leaping algorithm. We found that the concentration gradient can accelerate antimalarial resistance evolution. The gain in resistance evolution was improved by the increase in the parasite mutation rate and the mosquito biting rate. In addition, even though the rate of resistance evolution is not sensitive to the changes in parasite reduction ratios (PRRs) of antimalarial drugs, the probability of finding the antimalarial drug resistant parasites decreases when the PRR increases.

The genomic architecture of antimalarial drug resistance

Plasmodium falciparum and Plasmodium vivax, the two protozoan parasite species that cause the majority of cases of human malaria, have developed resistance to nearly all known antimalarials. The ability of malaria parasites to develop resistance is primarily due to the high numbers of parasites in the infected person's bloodstream during the asexual blood stage of infection in conjunction with the mutability of their genomes. Identifying the genetic mutations that mediate antimalarial resistance has deepened our understanding of how the parasites evade our treatments and reveals molecular markers that can be used to track the emergence of resistance in clinical samples. In this review, we examine known genetic mutations that lead to resistance to the major classes of antimalarial medications: the 4-aminoquinolines (chloroquine, amodiaquine and piperaquine), antifolate drugs, aryl amino-alcohols (quinine, lumefantrine and mefloquine), artemisinin compounds, antibiotics (clindamycin and doxycycline) and a napthoquinone (atovaquone). We discuss how the evolution of antimalarial resistance informs strategies to design the next generation of antimalarial therapies.

Covalent Plasmodium falciparum-selective proteasome inhibitors exhibit a low propensity for generating resistance in vitro and synergize with multiple antimalarial agents

Therapeutics with novel modes of action and a low risk of generating resistance are urgently needed to combat drug-resistant Plasmodium falciparum malaria. Here, we report that the peptide vinyl sulfones WLL-vs (WLL) and WLW-vs (WLW), highly selective covalent inhibitors of the P. falciparum proteasome, potently eliminate genetically diverse parasites, including K13-mutant, artemisinin-resistant lines, and are particularly active against ring-stage parasites. Selection studies reveal that parasites do not readily acquire resistance to WLL or WLW and that mutations in the β2, β5 or β6 subunits of the 20S proteasome core particle or in components of the 19S proteasome regulatory particle yield only <five-fold decreases in parasite susceptibility. This result compares favorably against previously published non-covalent inhibitors of the Plasmodium proteasome that can select for resistant parasites with >hundred-fold decreases in susceptibility. We observed no cross-resistance between WLL and WLW. Moreover, most mutations that conferred a modest loss of parasite susceptibility to one inhibitor significantly increased sensitivity to the other. These inhibitors potently synergized multiple chemically diverse classes of antimalarial agents, implicating a shared disruption of proteostasis in their modes of action. These results underscore the potential of targeting the Plasmodium proteasome with covalent small molecule inhibitors as a means of combating multidrug-resistant malaria.

The spread of artemisinin-resistant Plasmodium falciparum malaria across Southeast Asia creates an imperative to develop new treatment options with compounds that are not susceptible to existing mechanisms of antimalarial drug resistance. Recent work has identified the P. falciparum proteasome as a promising drug target. Here, we report potent antimalarial activity of highly selective vinyl sulfone-conjugated peptide proteasome inhibitors, including against artemisinin-resistant P. falciparum early ring-stage parasites that are traditionally difficult to treat. Unlike many advanced antimalarial candidates, these covalent proteasome inhibitors do not readily select for resistance. Selection studies with cultured parasites reveal infrequent and minor decreases in susceptibility resulting from point mutations in components of the 26S proteasome, which we model using cryo-electron microscopy-based structural data. No parasites were observed to be cross-resistant to both compounds; in fact, partial resistance to one compound often created hypersensitivity to the other. We also document potent synergy between these covalent proteasome inhibitors and multiple classes of antimalarial agents, including dihydroartemisinin, the clinical candidate OZ439, and the parasite transmission-blocking agent methylene blue. Proteasome inhibitors have significant promise as components of novel combination therapies to treat multidrug-resistant malaria.

Association of mutations in the Plasmodium falciparum Kelch13 gene (Pf3D7_1343700) with parasite clearance rates after artemisinin-based treatments-a WWARN individual patient data meta-analysis

BACKGROUND

Plasmodium falciparum infections with slow parasite clearance following artemisinin-based therapies are widespread in the Greater Mekong Subregion. A molecular marker of the slow clearance phenotype has been identified: single genetic changes within the propeller region of the Kelch13 protein (pfk13; Pf3D7_1343700). Global searches have identified almost 200 different non-synonymous mutant pfk13 genotypes. Most mutations occur at low prevalence and have uncertain functional significance. To characterize the impact of different pfk13 mutations on parasite clearance, we conducted an individual patient data meta-analysis of the associations between parasite clearance half-life (PC_1/2) and pfk13 genotype based on a large set of individual patient records from Asia and Africa.

METHODS

A systematic literature review following the PRISMA protocol was conducted to identify studies published between 2000 and 2017 which included frequent parasite counts and pfk13 genotyping. Four databases (Ovid Medline, PubMed, Ovid Embase, and Web of Science Core Collection) were searched. Eighteen studies (15 from Asia, 2 from Africa, and one multicenter study with sites on both continents) met inclusion criteria and were shared. Associations between the log transformed PC_1/2 values and pfk13 genotype were assessed using multivariable regression models with random effects for study site.

RESULTS

Both the pfk13 genotypes and the PC_1/2 were available from 3250 (95%) patients (n = 3012 from Asia (93%), n = 238 from Africa (7%)). Among Asian isolates, all pfk13 propeller region mutant alleles observed in five or more specific isolates were associated with a 1.5- to 2.7-fold longer geometric mean PC_1/2 compared to the PC_1/2 of wild type isolates (all p ≤ 0.002). In addition, mutant allele E252Q located in the P. falciparum region of pfk13 was associated with 1.5-fold (95%CI 1.4-1.6) longer PC_1/2. None of the isolates from four countries in Africa showed a significant difference between the PC_1/2 of parasites with or without pfk13 propeller region mutations. Previously, the association of six pfk13 propeller mutant alleles with delayed parasite clearance had been confirmed. This analysis demonstrates that 15 additional pfk13 alleles are associated strongly with the slow-clearing phenotype in Southeast Asia.

CONCLUSION

Pooled analysis associated 20 pfk13 propeller region mutant alleles with the slow clearance phenotype, including 15 mutations not confirmed previously.

Overexpression of plasmepsin II and plasmepsin III does not directly cause reduction in Plasmodium falciparum sensitivity to artesunate, chloroquine and piperaquine

Artemisinin derivatives and their partner drugs in artemisinin combination therapies (ACTs) have played a pivotal role in global malaria mortality reduction during the last two decades. The loss of artemisinin efficacy due to evolving drug-resistant parasites could become a serious global health threat. Dihydroartemisinin-piperaquine is a well tolerated and generally highly effective ACT. The implementation of a partner drug in ACTs is critical in the control of emerging artemisinin resistance. Even though artemisinin is highly effective in parasite clearance, it is labile in the human body. A partner drug is necessary for killing the remaining parasites when the pulses of artemisinin have ceased. A population of Plasmodium falciparum parasites in Cambodia and adjacent countries has become resistant to piperaquine. Increased copy number of the genes encoding the haemoglobinases Plasmepsin II and Plasmepsin III has been linked with piperaquine resistance by genome-wide association studies and in clinical trials, leading to the use of increased plasmepsin II/plasmepsin III copy number as a molecular marker for piperaquine resistance. Here we demonstrate that overexpression of plasmepsin II and plasmepsin III in the 3D7 genetic background failed to change the susceptibility of P. falciparum to artemisinin, chloroquine and piperaquine by both a standard dose-response analysis and a piperaquine survival assay. Whilst plasmepsin copy number polymorphism is currently implemented as a molecular surveillance resistance marker, further studies to discover the molecular basis of piperaquine resistance and potential epistatic interactions are needed.

Exploration of the Plasmodium falciparum Resistome and Druggable Genome Reveals New Mechanisms of Drug Resistance and Antimalarial Targets

Plasmodium parasites, the causative agent of malaria infections, rapidly evolve drug resistance and escape detection by the human immune response via the incredible mutability of its genome. Understanding the genetic mechanisms by which Plasmodium parasites develop antimalarial resistance is essential to understanding why most drugs fail in the clinic and designing the next generation of therapies. A systematic genomic analysis of 262 Plasmodium falciparum clones with stable in vitro resistance to 37 diverse compounds with potent antimalarial activity was undertaken with the main goal of identifying new drug targets. Despite several challenges inherent to this method of in vitro drug resistance generation followed by whole genome sequencing, the study was able to identify a likely drug target or resistance gene for every compound for which resistant parasites could be generated. Known and novel P falciparum resistance mediators were discovered along with several new promising antimalarial drug targets. Surprisingly, gene amplification events contributed to one-third of the drug resistance acquisition events. The study can serve as a model for drug discovery and resistance analyses in other similar microbial pathogens amenable to in vitro culture.

Amplification of GTP-cyclohydrolase 1 gene in Plasmodium falciparum isolates with the quadruple mutant of dihydrofolate reductase and dihydropteroate synthase genes in Ghana

Sulfadoxine-pyrimethamine (SP) is used as malaria chemoprophylaxis for pregnant women and children in Ghana. Plasmodium falciparum resistance to SP is linked to mutations in the dihydropteroate synthase gene (pfdhps), dihydrofolate reductase gene (pfdhfr) and amplification of GTP cyclohydrolase 1 (pfgch1) gene. The pfgch1 duplication is associated with pfdhfr L164, a crucial mutant for high level pyrimethamine resistance which is rare in Ghana. The presence of amplified pfgch1 in Ghanaian isolates could be an indicator of the evolution of the L164 mutant. This study therefore determined the pfgch1 copy number variations and SP resistance mutations in clinical isolates from Ghana. One hundred and ninety-two (192) blood samples collected from children aged ≤14 years with uncomplicated malaria in 2013-14 and 2015-16 were used. Quantitative real-time polymerase chain reaction (qRT-PCR) was used to detect the pfgch1 copy number and nested PCR-Sanger sequencing used to detect mutations in pfdhps and pfdhfr genes. Twelve parasites (6.3%) harbored double copies of the pfgch1 gene out of the 192 samples. Of the 12, 75% had the pfdhfr I51-R59-N108, 92% had the pfdhps G437 mutant, 8% had the pfdhps E540 and 67% had the SP resistance haplotype IRNG. No L164 was detected in samples with amplified pfgch1. The rare T108 mutant associated with cycloguanil resistance showed predominance (60%) over N108 in the 2015-16 isolates. The observation of parasites with increased copy number of pfgch1 gene is indicative of the future evolution of the rare quadruple pfdhfr mutant, I51-R59-N108-L164, in Ghanaian parasites. Mutant pfdhps isolates also had increased gch1 copy number suggestive that it may also facilitate sulphadoxine resistance. The selection of parasites with pfgch1 gene amplification will enhance the sustenance and persistence of parasites with SP resistance in the country. Policy makers need to begin the search for a replacement chemoprophylaxis drug for malaria vulnerable groups in Ghana.

Fitness Loss under Amino Acid Starvation in Artemisinin-Resistant Plasmodium falciparum Isolates from Cambodia

Artemisinin is the most rapidly effective drug for Plasmodium falciparum malaria treatment currently in clinical use. Emerging artemisinin-resistant parasites pose a great global health risk. At present, the level of artemisinin resistance is still relatively low with evidence pointing towards a trade-off between artemisinin resistance and fitness loss. Here we show that artemisinin-resistant P. falciparum isolates from Cambodia manifested fitness loss, showing fewer progenies during the intra-erythrocytic developmental cycle. The loss in fitness was exacerbated under the condition of low exogenous amino acid supply. The resistant parasites failed to undergo maturation, whereas their drug-sensitive counterparts were able to complete the erythrocytic cycle under conditions of amino acid deprivation. The artemisinin-resistant phenotype was not stable, and loss of the phenotype was associated with changes in the expression of a putative target, Exp1, a membrane glutathione transferase. Analysis of SNPs in haemoglobin processing genes revealed associations with parasite clearance times, suggesting changes in haemoglobin catabolism may contribute to artemisinin resistance. These findings on fitness and protein homeostasis could provide clues on how to contain emerging artemisinin-resistant parasites.

How to combat emerging artemisinin resistance: Lessons from "The Three Little Pigs"

It is rare to come across an Aesop’s fable in respectable journals. It might catch scientists outside the malaria field by surprise to learn that the famous story of “The Boy Who Cried Wolf” has been repeatedly compared to the threat from artemisinin-resistant malaria parasites, including the two latest reports on the rise of a specific haplotype in Cambodia and Thailand, sensationally dubbed “Super Malaria” by the media [1, 2]. The comparison to a children’s tale should not negate the fact that malaria drug resistance is one of the most pressing threats to the global public health community. Here, the findings leading to this contentious discourse will be delineated in order to provide a perspective. Possible solutions will be presented to stimulate further research and discussion to solve one of the greatest public health challenges of our lifetime.

Antimalarial drug resistance: linking Plasmodium falciparum parasite biology to the clinic

The global adoption of artemisinin-based combination therapies (ACTs) in the early 2000s heralded a new era in effectively treating drug-resistant Plasmodium falciparum malaria. However, several Southeast Asian countries have now reported the emergence of parasites that have decreased susceptibility to artemisinin (ART) derivatives and ACT partner drugs, resulting in increasing rates of treatment failures. Here we review recent advances in understanding how antimalarials act and how resistance develops, and discuss new strategies for effectively combatting resistance, optimizing treatment and advancing the global campaign to eliminate malaria.

Assessment of copy number variation in genes related to drug resistance in Plasmodium vivax and Plasmodium falciparum isolates from the Brazilian Amazon and a systematic review of the literature

Background

Parasite resistance to anti-malarials represents a great obstacle for malaria elimination. The majority of studies have investigated the association between single-nucleotide polymorphisms (SNPs) and drug resistance; however, it is becoming clear that the copy number variation (CNV) is also associated with this parasite phenotype. To provide a baseline for molecular surveillance of anti-malarial drug resistance in the Brazilian Amazon, the present study characterized the genetic profile of both markers in the most common genes associated with drug resistance in Plasmodium falciparum and Plasmodium vivax isolates. Additionally, these data were compared to data published elsewhere applying a systematic review of the literature published over a 20-year time period.

Methods

The genomic DNA of 67 patients infected by P. falciparum and P. vivax from three Brazilian States was obtained between 2002 and 2012. CNV in P. falciparum multidrug resistance gene-1 (pfmdr1), GTP cyclohydrolase 1 (pfgch1) and P. vivax multidrug resistance gene-1 (pvmdr1) were assessed by real-time PCR assays. SNPs in the pfmdr1 and pfcrt genes were assessed by PCR–RFLP. A literature search for studies that analysed CNP in the same genes of P. falciparum and P. vivax was conducted between May 2014 and March 2017 across four databases.

Results

All analysed samples of P. falciparum carried only one copy of pfmdr1 or pfgch1. Although the pfcrt K76T polymorphism, a determinant of CQ resistance, was present in all samples genotyped, the pfmdr1 N86Y was absent. For P. vivax isolates, an amplification rate of 20% was found for the pvmdr1 gene. The results of the study are in agreement with the low amplification rates for pfmdr1 gene evidenced in the Americas and Africa, while higher rates have been described in Southeast Asia. For P. vivax, very low rates of amplification for pvmdr1 have been described worldwide, with exceptions in French Guiana, Cambodia, Thailand and Brazil.

Conclusions

The present study was the first to evaluate gch1 CNV in P. falciparum isolates from Brazil, showing an absence of amplification of this gene more than 20 years after the withdrawal of the Brazilian antifolates therapeutic scheme. Furthermore, the rate of pvmdr1 amplification was significantly higher than that previously reported for isolates circulating in Northern Brazil.

Electronic supplementary material

The online version of this article (doi:10.1186/s12936-017-1806-z) contains supplementary material, which is available to authorized users.

Identifying antimalarial compounds targeting dihydrofolate reductase-thymidylate synthase (DHFR-TS) by chemogenomic profiling

The mode of action of many antimalarial drugs is unknown. Chemogenomic profiling is a powerful method to address this issue. This experimental approach entails disruption of gene function and phenotypic screening for changes in sensitivity to bioactive compounds. Here, we describe the application of reverse genetics for chemogenomic profiling in Plasmodium. Plasmodium falciparum parasites harbouring a transgenic insertion of the glmS ribozyme downstream of the dihydrofolate reductase-thymidylate synthase (DHFR-TS) gene were used for chemogenomic profiling of antimalarial compounds to identify those which target DHFR-TS. DHFR-TS expression can be attenuated by exposing parasites to glucosamine. Parasites with attenuated DHFR-TS expression were significantly more sensitive to antifolate drugs known to target DHFR-TS. In contrast, no change in sensitivity to other antimalarial drugs with different modes of action was observed. Chemogenomic profiling was performed using the Medicines for Malaria Venture (Switzerland) Malaria Box compound library, and two compounds were identified as novel DHFR-TS inhibitors. We also tested the glmS ribozyme in Plasmodium berghei, a rodent malaria parasite. The expression of reporter genes with downstream glmS ribozyme could be attenuated in transgenic parasites comparable with that obtained in P. falciparum. The chemogenomic profiling method was applied in a P. berghei line expressing a pyrimethamine-resistant Toxoplasma gondii DHFR-TS reporter gene under glmS ribozyme control. Parasites with attenuated expression of this gene were significantly sensitised to antifolates targeting DHFR-TS, but not other drugs with different modes of action. In conclusion, these data show that the glmS ribozyme reverse genetic tool can be applied for identifying primary targets of antimalarial compounds in human and rodent malaria parasites.

Interrogating the Plasmodium Sporozoite Surface: Identification of Surface-Exposed Proteins and Demonstration of Glycosylation on CSP and TRAP by Mass Spectrometry-Based Proteomics

Malaria parasite infection is initiated by the mosquito-transmitted sporozoite stage, a highly motile invasive cell that targets hepatocytes in the liver for infection. A promising approach to developing a malaria vaccine is the use of proteins located on the sporozoite surface as antigens to elicit humoral immune responses that prevent the establishment of infection. Very little of the P. falciparum genome has been considered as potential vaccine targets, and candidate vaccines have been almost exclusively based on single antigens, generating the need for novel target identification. The most advanced malaria vaccine to date, RTS,S, a subunit vaccine consisting of a portion of the major surface protein circumsporozoite protein (CSP), conferred limited protection in Phase III trials, falling short of community-established vaccine efficacy goals. In striking contrast to the limited protection seen in current vaccine trials, sterilizing immunity can be achieved by immunization with radiation-attenuated sporozoites, suggesting that more potent protection may be achievable with a multivalent protein vaccine. Here, we provide the most comprehensive analysis to date of proteins located on the surface of or secreted by Plasmodium falciparum salivary gland sporozoites. We used chemical labeling to isolate surface-exposed proteins on sporozoites and identified these proteins by mass spectrometry. We validated several of these targets and also provide evidence that components of the inner membrane complex are in fact surface-exposed and accessible to antibodies in live sporozoites. Finally, our mass spectrometry data provide the first direct evidence that the Plasmodium surface proteins CSP and TRAP are glycosylated in sporozoites, a finding that could impact the selection of vaccine antigens.

Malaria remains one of the most important infectious diseases in the world, responsible for an estimated 500 million new cases and 600,000 deaths annually. The etiologic agents of the disease are protozoan parasites of the genus Plasmodium that have a complex cycle between mosquito and mammalian hosts. Though all clinical symptoms are attributable to the blood stages, it is only by attacking the transmission stages that we can make an impact on the economic and health burdens of malaria. Infection is initiated when mosquitoes inoculate sporozoites into the skin as they probe for blood. Sporozoites must locate blood vessels and enter the circulation to reach the liver where they invade and grow in hepatocytes. The inoculum is low and these early stages of infection are asymptomatic. Though the small amounts of material available for study has made large scale -omics studies difficult, killing the parasite at this stage would prevent infection and block downstream transmission to mosquitoes, thus preventing spread of disease. Here we use state-of-the-art biochemistry tools to identify the proteins on the sporozoite surface and find that two of the most studied proteins, CSP and TRAP, have post-translational modifications. These studies will aid investigations into the novel biology of sporozoites and importantly, significantly expand the pool of potential vaccine candidates.

Within-host competition and drug resistance in the human malaria parasite Plasmodium falciparum

Infections with the malaria parasite Plasmodium falciparum typically comprise multiple strains, especially in high-transmission areas where infectious mosquito bites occur frequently. However, little is known about the dynamics of mixed-strain infections, particularly whether strains sharing a host compete or grow independently. Competition between drug-sensitive and drug-resistant strains, if it occurs, could be a crucial determinant of the spread of resistance. We analysed 1341 P. falciparum infections in children from Angola, Ghana and Tanzania and found compelling evidence for competition in mixed-strain infections: overall parasite density did not increase with additional strains, and densities of individual chloroquine-sensitive (CQS) and chloroquine-resistant (CQR) strains were reduced in the presence of competitors. We also found that CQR strains exhibited low densities compared with CQS strains (in the absence of chloroquine), which may underlie observed declines of chloroquine resistance in many countries following retirement of chloroquine as a first-line therapy. Our observations support a key role for within-host competition in the evolution of drug-resistant malaria. Malaria control and resistance-management efforts in high-transmission regions may be significantly aided or hindered by the effects of competition in mixed-strain infections. Consideration of within-host dynamics may spur development of novel strategies to minimize resistance while maximizing the benefits of control measures.

Combinatorial Genetic Modeling of pfcrt-Mediated Drug Resistance Evolution in Plasmodium falciparum

The emergence of drug resistance continuously threatens global control of infectious diseases, including malaria caused by the protozoan parasite Plasmodium falciparum. A critical parasite determinant is the P. falciparum chloroquine resistance transporter (PfCRT), the primary mediator of chloroquine (CQ) resistance (CQR), and a pleiotropic modulator of susceptibility to several first-line artemisinin-based combination therapy partner drugs. Aside from the validated CQR molecular marker K76T, P. falciparum parasites have acquired at least three additional pfcrt mutations, whose contributions to resistance and fitness have been heretofore unclear. Focusing on the quadruple-mutant Ecuadorian PfCRT haplotype Ecu1110 (K76T/A220S/N326D/I356L), we genetically modified the pfcrt locus of isogenic, asexual blood stage P. falciparum parasites using zinc-finger nucleases, producing all possible combinations of intermediate pfcrt alleles. Our analysis included the related quintuple-mutant PfCRT haplotype 7G8 (Ecu1110 + C72S) that is widespread throughout South America and the Western Pacific. Drug susceptibilities and in vitro growth profiles of our combinatorial pfcrt-modified parasites were used to simulate the mutational trajectories accessible to parasites as they evolved CQR. Our results uncover unique contributions to parasite drug resistance and growth for mutations beyond K76T and predict critical roles for the CQ metabolite monodesethyl-CQ and the related quinoline-type drug amodiaquine in driving mutant pfcrt evolution. Modeling outputs further highlight the influence of parasite proliferation rates alongside gains in drug resistance in dictating successful trajectories. Our findings suggest that P. falciparum parasites have navigated constrained pfcrt adaptive landscapes by means of probabilistically rare mutational bursts that led to the infrequent emergence of pfcrt alleles in the field.

Comparative genome-wide analysis and evolutionary history of haemoglobin-processing and haem detoxification enzymes in malarial parasites

Background

Malaria parasites have evolved a series of intricate mechanisms to survive and propagate within host red blood cells. Intra-erythrocytic parasitism requires these organisms to digest haemoglobin and detoxify iron-bound haem. These tasks are executed by haemoglobin-specific proteases and haem biocrystallization factors that are components of a large multi-subunit complex. Since haemoglobin processing machineries are functionally and genetically linked to the modes of action and resistance mechanisms of several anti-malarial drugs, an understanding of their evolutionary history is important for drug development and drug resistance prevention.

Methods

Maximum likelihood trees of genetic repertoires encoding haemoglobin processing machineries within Plasmodium species, and with the representatives of Apicomplexan species with various host tropisms, were created. Genetic variants were mapped onto existing three-dimensional structures. Genome-wide single nucleotide polymorphism data were used to analyse the selective pressure and the effect of these mutations at the structural level.

Results

Recent expansions in the falcipain and plasmepsin repertoires are unique to human malaria parasites especially in the Plasmodium falciparum and P. reichenowi lineage. Expansion of haemoglobin-specific plasmepsins occurred after the separation event of Plasmodium species, but the other members of the plasmepsin family were evolutionarily conserved with one copy for each sub-group in every Apicomplexan species. Haemoglobin-specific falcipains are separated from invasion-related falcipain, and their expansions within one specific locus arose independently in both P. falciparum and P. vivax lineages. Gene conversion between P. falciparum falcipain 2A and 2B was observed in artemisinin-resistant strains. Comparison between the numbers of non-synonymous and synonymous mutations suggests a strong selective pressure at falcipain and plasmepsin genes. The locations of amino acid changes from non-synonymous mutations mapped onto protein structures revealed clusters of amino acid residues in close proximity or near the active sites of proteases.

Conclusion

A high degree of polymorphism at the haemoglobin processing genes implicates an imposition of selective pressure. The identification in recent years of functional redundancy of haemoglobin-specific proteases makes them less appealing as potential drug targets, but their expansions, especially in the human malaria parasite lineages, unequivocally point toward their functional significance during the independent and repetitive adaptation events in malaria parasite evolutionary history.

Electronic supplementary material

The online version of this article (doi:10.1186/s12936-016-1097-9) contains supplementary material, which is available to authorized users.

Genome-wide transcriptome profiling reveals functional networks involving the Plasmodium falciparum drug resistance transporters PfCRT and PfMDR1

BACKGROUND

The acquisition of multidrug resistance by Plasmodium falciparum underscores the need to understand the underlying molecular mechanisms so as to counter their impact on malaria control. For the many antimalarials whose mode of action relates to inhibition of heme detoxification inside infected erythrocytes, the digestive vacuole transporters PfCRT and PfMDR1 constitute primary resistance determinants.

RESULTS

Using gene expression microarrays over the course of the parasite intra-erythrocytic developmental cycle, we compared the transcriptomic profiles between P. falciparum strains displaying mutant or wild-type pfcrt or varying in pfcrt or pfmdr1 expression levels. To account for differences in the time of sampling, we developed a computational method termed Hypergeometric Analysis of Time Series, which combines Fast Fourier Transform with a modified Gene Set Enrichment Analysis. Our analysis revealed coordinated changes in genes involved in protein catabolism, translation initiation and DNA/RNA metabolism. We also observed differential expression of genes with a role in transport or coding for components of the digestive vacuole. Interestingly, a global comparison of all profiled transcriptomes uncovered a tight correlation between the transcript levels of pfcrt and pfmdr1, extending to dozens of other genes, suggesting an intricate regulatory balance in order to maintain optimal physiological processes.

CONCLUSIONS

This study provides insight into the mechanisms by which P. falciparum adjusts to the acquisition of mutations or gene amplification in key transporter loci that mediate drug resistance. Our results implicate several biological pathways that may be differentially regulated to compensate for impaired transporter function and alterations in parasite vacuole physiology.

Population Structure Shapes Copy Number Variation in Malaria Parasites

If copy number variants (CNVs) are predominantly deleterious, we would expect them to be more efficiently purged from populations with a large effective population size (Ne) than from populations with a small Ne. Malaria parasites (Plasmodium falciparum) provide an excellent organism to examine this prediction, because this protozoan shows a broad spectrum of population structures within a single species, with large, stable, outbred populations in Africa, small unstable inbred populations in South America and with intermediate population characteristics in South East Asia. We characterized 122 single-clone parasites, without prior laboratory culture, from malaria-infected patients in seven countries in Africa, South East Asia and South America using a high-density single-nucleotide polymorphism/CNV microarray. We scored 134 high-confidence CNVs across the parasite exome, including 33 deletions and 102 amplifications, which ranged in size from <500 bp to 59 kb, as well as 10,107 flanking, biallelic single-nucleotide polymorphisms. Overall, CNVs were rare, small, and skewed toward low frequency variants, consistent with the deleterious model. Relative to African and South East Asian populations, CNVs were significantly more common in South America, showed significantly less skew in allele frequencies, and were significantly larger. On this background of low frequency CNV, we also identified several high-frequency CNVs under putative positive selection using an FST outlier analysis. These included known adaptive CNVs containing rh2b and pfmdr1, and several other CNVs (e.g., DNA helicase and three conserved proteins) that require further investigation. Our data are consistent with a significant impact of genetic structure on CNV burden in an important human pathogen.

The molecular basis of antifolate resistance in Plasmodium falciparum: looking beyond point mutations

Drugs that target the folate-synthesis pathway have a long history of effectiveness against a variety of pathogens. As antimalarials, the antifolates were safe and well tolerated, but resistance emerged quickly and has persisted even with decreased drug pressure. The primary determinants of resistance in Plasmodium falciparum are well-described point mutations in the enzymes dihydropteroate synthase and dihydrofolate reductase targeted by the combination sulfadoxine-pyrimethamine. Recent work has highlighted the contributions of additional parasite adaptation to antifolate resistance. In fact, the evolution of antifolate-resistant parasites is multifaceted and complex. Gene amplification of the first enzyme in the parasite folate synthesis pathway, GTP-cyclohydrolase, is strongly associated with resistant parasites and potentially contributes to persistence of resistant parasites. Further understanding of how parasites adjust flux through the folate pathway is important to the further development of alternative agents targeting this crucial synthesis pathway.

Effective temperature of mutations

Biological macromolecules experience two seemingly very different types of noise acting on different time scales: (i) point mutations corresponding to changes in molecular sequence and (ii) thermal fluctuations. Examining the secondary structures of a large number of microRNA precursor sequences and model lattice proteins, we show that the effects of single point mutations are statistically indistinguishable from those of an increase in temperature by a few tens of kelvins. The existence of such an effective mutational temperature establishes a quantitative connection between robustness to genetic (mutational) and environmental (thermal) perturbations.

Biochemical and functional characterization of Plasmodium falciparum GTP cyclohydrolase I

Background

Antifolates are currently in clinical use for malaria preventive therapy and treatment. The drugs kill the parasites by targeting the enzymes in the de novo folate pathway. The use of antifolates has now been limited by the spread of drug-resistant mutations. GTP cyclohydrolase I (GCH1) is the first and the rate-limiting enzyme in the folate pathway. The amplification of the gch1 gene found in certain Plasmodium falciparum isolates can cause antifolate resistance and influence the course of antifolate resistance evolution. These findings showed the importance of P. falciparum GCH1 in drug resistance intervention. However, little is known about P. falciparum GCH1 in terms of kinetic parameters and functional assays, precluding the opportunity to obtain the key information on its catalytic reaction and to eventually develop this enzyme as a drug target.

Methods

Plasmodium falciparum GCH1 was cloned and expressed in bacteria. Enzymatic activity was determined by the measurement of fluorescent converted neopterin with assay validation by using mutant and GTP analogue. The genetic complementation study was performed in ∆folE bacteria to functionally identify the residues and domains of P. falciparum GCH1 required for its enzymatic activity. Plasmodial GCH1 sequences were aligned and structurally modeled to reveal conserved catalytic residues.

Results

Kinetic parameters and optimal conditions for enzymatic reactions were determined by the fluorescence-based assay. The inhibitor test against P. falciparum GCH1 is now possible as indicated by the inhibitory effect by 8-oxo-GTP. Genetic complementation was proven to be a convenient method to study the function of P. falciparum GCH1. A series of domain truncations revealed that the conserved core domain of GCH1 is responsible for its enzymatic activity. Homology modelling fits P. falciparum GCH1 into the classic Tunnelling-fold structure with well-conserved catalytic residues at the active site.

Conclusions

Functional assays for P. falciparum GCH1 based on enzymatic activity and genetic complementation were successfully developed. The assays in combination with a homology model characterized the enzymatic activity of P. falciparum GCH1 and the importance of its key amino acid residues. The potential to use the assay for inhibitor screening was validated by 8-oxo-GTP, a known GTP analogue inhibitor.