Protein-based signatures of functional evolution in Plasmodium falciparum

Forces acting on codon bias in malaria parasites

Malaria parasite genomes have a range of codon biases, with Plasmodium falciparum one of the most AT-biased genomes known. We examined the make up of synonymous coding sites and stop codons in the core genomes of representative malaria parasites, showing first that local DNA context influences codon bias similarly across P. falciparum, P. vivax and P. berghei, with suppression of CpG dinucleotides and enhancement of CpC dinucleotides, both within and aross codons. Intense asexual phase gene expression in P. falciparum and P. berghei is associated with increased A3:G3 bias but reduced T3:C3 bias at 2-fold sites, consistent with adaptation of codons to tRNA pools and avoidance of wobble tRNA interactions that potentially slow down translation. In highly expressed genes, the A3:G3 ratio can exceed 30-fold while the T3:C3 ratio can be less than 1, according to the encoded amino acid and subsequent base. Lysine codons (AAA/G) show distinctive behaviour with substantially reduced A3:G3 bias in highly expressed genes, perhaps because of selection against frameshifting when the AAA codon is followed by another adenine. Intense expression is also associated with a strong bias towards TAA stop codons (found in 94% and 89% of highly expressed P. falciparum and P. berghei genes respectively) and a proportional rise in the TAAA stop ‘tetranucleotide’. The presence of these expression-linked effects in the relatively AT-rich malaria parasite species adds weight to the suggestion that AT-richness in the Plasmodium genus might be a fitness adaptation. Potential explanations for the relative lack of codon bias in P. vivax include the distinct features of its lifecycle and its effective population size over evolutionary time.

South-East Asian strains of Plasmodium falciparum display higher ratio of non-synonymous to synonymous polymorphisms compared to African strains

Resistance to frontline anti-malarial drugs, including artemisinin, has repeatedly arisen in South-East Asia, but the reasons for this are not understood. Here we test whether evolutionary constraints on Plasmodium falciparum strains from South-East Asia differ from African strains. We find a significantly higher ratio of non-synonymous to synonymous polymorphisms in P. falciparum from South-East Asia compared to Africa, suggesting differences in the selective constraints on P. falciparum genome in these geographical regions. Furthermore, South-East Asian strains showed a higher proportion of non-synonymous polymorphism at conserved positions, suggesting reduced negative selection. There was a lower rate of mixed infection by multiple genotypes in samples from South-East Asia compared to Africa. We propose that a lower mixed infection rate in South-East Asia reduces intra-host competition between the parasite clones, reducing the efficiency of natural selection. This might increase the probability of fixation of fitness-reducing mutations including drug resistant ones.

The clinical impact of artemisinin resistance in Southeast Asia and the potential for future spread

Artemisinins are the most rapidly acting of currently available antimalarial drugs. Artesunate has become the treatment of choice for severe malaria, and artemisinin-based combination therapies (ACTs) are the foundation of modern falciparum malaria treatment globally. Their safety and tolerability profile is excellent. Unfortunately, Plasmodium falciparum infections with mutations in the ‘K13’ gene, with reduced ring-stage susceptibility to artemisinins, and slow parasite clearance in patients treated with ACTs, are now widespread in Southeast Asia. We review clinical efficacy data from the region (2000–2015) that provides strong evidence that the loss of first-line ACTs in western Cambodia, first artesunate-mefloquine and then DHA-piperaquine, can be attributed primarily to K13 mutated parasites. The ring-stage activity of artemisinins is therefore critical for the sustained efficacy of ACTs; once it is lost, rapid selection of partner drug resistance and ACT failure are inevitable consequences. Consensus methods for monitoring artemisinin resistance are now available. Despite increased investment in regional control activities, ACTs are failing across an expanding area of the Greater Mekong subregion. Although multiple K13 mutations have arisen independently, successful multidrug-resistant parasite genotypes are taking over and threaten to spread to India and Africa. Stronger containment efforts and new approaches to sustaining long-term efficacy of antimalarial regimens are needed to prevent a global malaria emergency.

Artemisinin resistance in Plasmodium falciparum malaria is causing failure of artemisinin-based combination therapies across an expanding area of Southeast Asia, undermining control and elimination efforts. The potential global consequences can only be avoided by new approaches that ensure sustained efficacy for antimalarial regimens in malaria affected populations.

Genomic epidemiology of artemisinin resistant malaria

The current epidemic of artemisinin resistant Plasmodium falciparum in Southeast Asia is the result of a soft selective sweep involving at least 20 independent kelch13 mutations. In a large global survey, we find that kelch13 mutations which cause resistance in Southeast Asia are present at low frequency in Africa. We show that African kelch13 mutations have originated locally, and that kelch13 shows a normal variation pattern relative to other genes in Africa, whereas in Southeast Asia there is a great excess of non-synonymous mutations, many of which cause radical amino-acid changes. Thus, kelch13 is not currently undergoing strong selection in Africa, despite a deep reservoir of variations that could potentially allow resistance to emerge rapidly. The practical implications are that public health surveillance for artemisinin resistance should not rely on kelch13 data alone, and interventions to prevent resistance must account for local evolutionary conditions, shown by genomic epidemiology to differ greatly between geographical regions.

DOI:http://dx.doi.org/10.7554/eLife.08714.001

Malaria is an infectious disease caused by a microscopic parasite called Plasmodium, which is transferred between humans by mosquitos. One species of malaria parasite called Plasmodium falciparum can cause particularly severe and life-threatening forms of the disease. Currently, the most widely used treatment for P. falciparum infections is artemisinin combination therapy, a treatment that combines the drug artemisinin (or a closely related molecule) with another antimalarial drug. However, resistance to artemisinin has started to spread throughout Southeast Asia.

Artemisinin resistance is caused by mutations in a parasite gene called kelch13, and researchers have identified over 20 different mutations in P. falciparum that confer artemisinin resistance. The diversity of mutations involved, and the fact that the same mutation can arise independently in different locations, make it difficult to track the spread of resistance using conventional molecular marker approaches.

Here, Amato, Miotto et al. sequenced the entire genomes of more than 3,000 clinical samples of P. falciparum from Southeast Asia and Africa, collected as part of a global network of research groups called the MalariaGEN Plasmodium falciparum Community Project. Amato, Miotto et al. found that African parasites had independently acquired many of the same kelch13 mutations that are known to cause resistance to artemisinin in Southeast Asia. However the kelch13 mutations seen in Africa remained at low levels in the parasite population, and appeared to be under much less pressure for evolutionary selection than those found in Southeast Asia.

These findings demonstrate that the emergence and spread of resistance to antimalarial drugs does not depend solely on the mutational process, but also on other factors that influence whether the mutations will spread in the population. Understanding how this is affected by different patterns of drug treatments and other environmental conditions will be important in developing more effective strategies for combating malaria.

DOI:http://dx.doi.org/10.7554/eLife.08714.002

K13 mutations and pfmdr1 copy number variation in Plasmodium falciparum malaria in Myanmar

Background

Artemisinin-based combination therapy has been first-line treatment for falciparum malaria in Myanmar since 2005. The wide extent of artemisinin resistance in the Greater Mekong sub-region and the presence of mefloquine resistance at the Myanmar-Thailand border raise concerns over resistance patterns in Myanmar. The availability of molecular markers for resistance to both drugs enables assessment even in remote malaria-endemic areas.

Methods

A total of 250 dried blood spot samples collected from patients with Plasmodium falciparum malarial infection in five malaria-endemic areas across Myanmar were analysed for kelch 13 sequence (k13) and pfmdr1 copy number variation. K13 mutations in the region corresponding to amino acids 210–726 (including the propeller region of the protein) were detected by nested PCR amplification and sequencing, and pfmdr1 copy number variation by real-time PCR. In two sites, a sub-set of patients were prospectively followed up for assessment of day-3 parasite clearance rates after a standard course of artemether-lumefantrine.

Results

K13 mutations and pfmdr1 amplification were successfully analysed in 206 and 218 samples, respectively. Sixty-nine isolates (33.5 %) had mutations within the k13 propeller region with 53 of these (76.8 %) having mutations already known to be associated with artemisinin resistance. F446I (32 isolates) and P574L (15 isolates) were the most common examples. K13 mutation was less common in sites in western border regions (29 of 155 isolates) compared to samples from the east and north (40 of 51 isolates; p < 0.0001). The overall proportion of parasites with multiple pfmdr1 copies (greater than 1.5) was 5.5 %. Seven samples showed both k13 mutation and multiple copies of pfmdr1. Only one of 36 patients followed up after artemether-lumefantrine treatment still had parasites at day 3; molecular analysis indicated wild-type k13 and single copy pfmdr1.

Conclusion

The proportion of P. falciparum isolates with mutations in the propeller region of k13 indicates that artemisinin resistance extends across much of Myanmar. There is a low prevalence of parasites with multiple pfmdr1 copies across the country. The efficacy of artemisinin-based combination therapy containing mefloquine and lumefantrine is, therefore, expected to be high, although regular monitoring of efficacy will be important.

Electronic supplementary material

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

The malaria parasite cation ATPase PfATP4 and its role in the mechanism of action of a new arsenal of antimalarial drugs

The intraerythrocytic malaria parasite, Plasmodium falciparum, maintains a low cytosolic Na(+) concentration and the plasma membrane P-type cation translocating ATPase 'PfATP4' has been implicated as playing a key role in this process. PfATP4 has been the subject of significant attention in recent years as mutations in this protein confer resistance to a growing number of new antimalarial compounds, including the spiroindolones, the pyrazoles, the dihydroisoquinolones, and a number of the antimalarial agents in the Medicines for Malaria Venture's 'Malaria Box'. On exposure of parasites to these compounds there is a rapid disruption of cytosolic Na(+). Whether, and if so how, such chemically distinct compounds interact with PfATP4, and how such interactions lead to parasite death, is not yet clear. The fact that multiple different chemical classes have converged upon PfATP4 highlights its significance as a potential target for new generation antimalarial agents. A spiroindolone (KAE609, now known as cipargamin) has progressed through Phase I and IIa clinical trials with favourable results. In this review we consider the physiological role of PfATP4, summarise the current repertoire of antimalarial compounds for which PfATP4 is implicated in their mechanism of action, and provide an outlook on translation from target identification in the laboratory to patient treatment in the field.

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.

An evolutionary perspective on the kinome of malaria parasites

Malaria parasites belong to an ancient lineage that diverged very early from the main branch of eukaryotes. The approximately 90-member plasmodial kinome includes a majority of eukaryotic protein kinases that clearly cluster within the AGC, CMGC, TKL, CaMK and CK1 groups found in yeast, plants and mammals, testifying to the ancient ancestry of these families. However, several hundred millions years of independent evolution, and the specific pressures brought about by first a photosynthetic and then a parasitic lifestyle, led to the emergence of unique features in the plasmodial kinome. These include taxon-restricted kinase families, and unique peculiarities of individual enzymes even when they have homologues in other eukaryotes. Here, we merge essential aspects of all three malaria-related communications that were presented at the Evolution of Protein Phosphorylation meeting, and propose an integrated discussion of the specific features of the parasite's kinome and phosphoproteome.