Experimental evolution, genetic analysis and genome re-sequencing reveal the mutation conferring artemisinin resistance in an isogenic lineage of malaria parasites

Genomics and integrated systems biology in Plasmodium falciparum: a path to malaria control and eradication

The first draft of the human malaria parasite's genome was released in 2002. Since then, the malaria scientific community has witnessed a steady embrace of new and powerful functional genomic studies. Over the years, these approaches have slowly revolutionized malaria research and enabled the comprehensive, unbiased investigation of various aspects of the parasite's biology. These genome-wide analyses delivered a refined annotation of the parasite's genome, delivered a better knowledge of its RNA, proteins and metabolite derivatives, and fostered the discovery of new vaccine and drug targets. Despite the positive impacts of these genomic studies, most research and investment still focus on protein targets, drugs and vaccine candidates that were known before the publication of the parasite genome sequence. However, recent access to next-generation sequencing technologies, along with an increased number of genome-wide applications, is expanding the impact of the parasite genome on biomedical research, contributing to a paradigm shift in research activities that may possibly lead to new optimized diagnosis and treatments. This review provides an update of Plasmodium falciparum genome sequences and an overview of the rapid development of genomics and system biology applications that have an immense potential of creating powerful tools for a successful malaria eradication campaign.

A major genome region underlying artemisinin resistance in malaria

Evolving resistance to artemisinin-based compounds threatens to derail attempts to control malaria. Resistance has been confirmed in western Cambodia and has recently emerged in western Thailand, but is absent from neighboring Laos. Artemisinin resistance results in reduced parasite clearance rates (CRs) after treatment. We used a two-phase strategy to identify genome region(s) underlying this ongoing selective event. Geographical differentiation and haplotype structure at 6969 polymorphic single-nucleotide polymorphisms (SNPs) in 91 parasites from Cambodia, Thailand, and Laos identified 33 genome regions under strong selection. We screened SNPs and microsatellites within these regions in 715 parasites from Thailand, identifying a selective sweep on chromosome 13 that shows strong association (P = 10(-6) to 10(-12)) with slow CRs, illustrating the efficacy of targeted association for identifying the genetic basis of adaptive traits.

Quantitative genome re-sequencing defines multiple mutations conferring chloroquine resistance in rodent malaria

Background

Drug resistance in the malaria parasite Plasmodium falciparum severely compromises the treatment and control of malaria. A knowledge of the critical mutations conferring resistance to particular drugs is important in understanding modes of drug action and mechanisms of resistances. They are required to design better therapies and limit drug resistance.

A mutation in the gene (pfcrt) encoding a membrane transporter has been identified as a principal determinant of chloroquine resistance in P. falciparum, but we lack a full account of higher level chloroquine resistance. Furthermore, the determinants of resistance in the other major human malaria parasite, P. vivax, are not known. To address these questions, we investigated the genetic basis of chloroquine resistance in an isogenic lineage of rodent malaria parasite P. chabaudi in which high level resistance to chloroquine has been progressively selected under laboratory conditions.

Results

Loci containing the critical genes were mapped by Linkage Group Selection, using a genetic cross between the high-level chloroquine-resistant mutant and a genetically distinct sensitive strain. A novel high-resolution quantitative whole-genome re-sequencing approach was used to reveal three regions of selection on chr11, chr03 and chr02 that appear progressively at increasing drug doses on three chromosomes. Whole-genome sequencing of the chloroquine-resistant parent identified just four point mutations in different genes on these chromosomes. Three mutations are located at the foci of the selection valleys and are therefore predicted to confer different levels of chloroquine resistance. The critical mutation conferring the first level of chloroquine resistance is found in aat1, a putative aminoacid transporter.

Conclusions

Quantitative trait loci conferring selectable phenotypes, such as drug resistance, can be mapped directly using progressive genome-wide linkage group selection. Quantitative genome-wide short-read genome resequencing can be used to reveal these signatures of drug selection at high resolution. The identities of three genes (and mutations within them) conferring different levels of chloroquine resistance generate insights regarding the genetic architecture and mechanisms of resistance to chloroquine and other drugs. Importantly, their orthologues may now be evaluated for critical or accessory roles in chloroquine resistance in human malarias P. vivax and P. falciparum.

A high-coverage artificial chromosome library for the genome-wide screening of drug-resistance genes in malaria parasites

The global spread of drug-resistant parasites is a serious problem for the treatment of malaria. Although identifying drug-resistance genes is crucial for the efforts against resistant parasites, an effective approach has not yet been developed. Here, we report a robust method for identifying resistance genes from parasites by using a Plasmodium artificial chromosome (PAC). Large genomic DNA fragments (10–50 kb) from the drug-resistant rodent malaria parasite Plasmodium berghei were ligated into the PAC and directly introduced into the drug-sensitive (i.e., wild-type) parasite by electroporation, resulting in a PAC library that encompassed the whole genomic sequence of the parasite. Subsequently, the transformed parasites that acquired resistance were selected by screening with the drug, and the resistance gene in the PAC was successfully identified. Furthermore, the drug-resistance gene was identified from a PAC library that was made from the pyrimethamine-resistant parasite Plasmodium chabaudi, further demonstrating the utility of our method. This method will promote the identification of resistance genes and contribute to the global fight against drug-resistant parasites.

Multiple mutations in heterogeneous miltefosine-resistant Leishmania major population as determined by whole genome sequencing

Background

Miltefosine (MF) is the first oral compound used in the chemotherapy against leishmaniasis. Since the mechanism of action of this drug and the targets of MF in Leishmania are unclear, we generated in a step-by-step manner Leishmania major promastigote mutants highly resistant to MF. Two of the mutants were submitted to a short-read whole genome sequencing for identifying potential genes associated with MF resistance.

Methods/Principal Findings

Analysis of the genome assemblies revealed several independent point mutations in a P-type ATPase involved in phospholipid translocation. Mutations in two other proteins—pyridoxal kinase and α-adaptin like protein—were also observed in independent mutants. The role of these proteins in the MF resistance was evaluated by gene transfection and gene disruption and both the P-type ATPase and pyridoxal kinase were implicated in MF susceptibility. The study also highlighted that resistance can be highly heterogeneous at the population level with individual clones derived from this population differing both in terms of genotypes but also susceptibility phenotypes.

Conclusions/Significance

Whole genome sequencing was used to pinpoint known and new resistance markers associated with MF resistance in the protozoan parasite Leishmania. The study also demonstrated the polyclonal nature of a resistant population with individual cells with varying susceptibilities and genotypes.

Leishmania spp. are parasitic protozoa responsible for a spectrum of diseases known as leishmaniasis. There are few drugs available for the treatment of these diseases, and miltefosine is the first oral drug used in treatment of visceral leishmaniasis, a form of the disease that can be lethal if not treated. In this study, we seek to understand the mechanism of action and identify targets of the drug by generating promastigote mutants highly resistant to miltefosine. Two independent mutants were submitted to short read whole genome sequencing. Genome analysis of these mutants has permitted us to identify point mutations in three genes (P-type ATPase, pyridoxal kinase and α-adaptin like protein) that were also present in other independent miltefosine resistant mutants. Some of the new genes identified here could be useful as potential markers for miltefosine resistance in Leishmania. Moreover, our approach has permitted us to highlight that resistance can be highly heterogeneous at the population level with individual clones derived from this population differing both in terms of genotypes but also susceptibility phenotypes. This may have practical applications while studying resistance.

Recent clinical and molecular insights into emerging artemisinin resistance in Plasmodium falciparum

PURPOSE OF REVIEW

Artemisinin-based combination therapies (ACTs) have been deployed globally with remarkable success for more than 10 years without having lost their malaria treatment efficacy. However, recent reports from the Thai-Cambodian border reveal evidence of emerging resistance to artemisinins. The latest published clinical and molecular findings are summarized herein.

RECENT FINDINGS

Clinical studies have identified delayed parasite clearance time as the most robust marker of artemisinin resistance. Resistance has only been documented from South-east Asia and has been observed in isolates that show no significant decrease in drug susceptibility in vitro. Genetic investigations have yet to uncover robust molecular markers. In-vitro studies have identified parasite quiescence or dormancy mechanisms that protect early 'ring-stage' intra-erythrocytic parasites against short-term artemisinin exposure. This might be achieved by reducing the rate of hemoglobin degradation, important for artemisinin bioactivation.

SUMMARY

Should ACTs fail, no suitable alternatives exist as first-line treatments of P. falciparum malaria. Intensified efforts are essential to monitor the spread of resistance, define therapeutic and operational strategies to counter its impact, and understand its molecular basis. Success in these areas is critical to ensuring that recent gains in reducing the burden of malaria are not lost.

Phenotypic and genotypic analysis of in vitro-selected artemisinin-resistant progeny of Plasmodium falciparum

Emergence of artemisinin resistance in Cambodia highlights the importance of characterizing resistance to this class of drugs. Previously, intermediate levels of resistance in Plasmodium falciparum were generated in vitro for artelinic acid (AL) and artemisinin (QHS). Here we expanded on earlier selection efforts to produce levels of clinically relevant concentrations, and the resulting lines were characterized genotypically and phenotypically. Recrudescence assays determined the ability of resistant and parent lines to recover following exposure to clinically relevant levels of drugs. Interestingly, the parent clone (D6) tolerated up to 1,500 ng/ml QHS, but the resistant parasite, D6.QHS340×3, recovered following exposure to 2,400 ng/ml QHS. Resistant D6, W2, and TM91c235 parasites all exhibited elevated 50% inhibitory concentrations (IC(50)s) to multiple artemisinin drugs, with >3-fold resistance to QHS and AL; however, the degree of resistance obtained with standard methods was remarkably less than expected for parasite lines that recovered from 2,400-ng/ml drug pressure. A novel assay format with radiolabeled hypoxanthine demonstrated a greater degree of resistance in vitro than the standard SYBR green method. Analysis of merozoite number in resistant parasites found D6 and TM91c235 resistant progeny had significantly fewer merozoites than parent strains, whereas W2 resistant progeny had significantly more. Amplification of pfmdr1 increased proportionately to the increased drug levels tolerated by W2 and TM91c235, but not in resistant D6. In summary, we define the artemisinin resistance phenotype as a decrease in susceptibility to artemisinins along with the ability to recover from drug-induced dormancy following supraclinical concentrations of the drug.

Artemisinin resistance in Plasmodium falciparum is associated with an altered temporal pattern of transcription

BACKGROUND

Artemisinin resistance in Plasmodium falciparum malaria has emerged in Western Cambodia. This is a major threat to global plans to control and eliminate malaria as the artemisinins are a key component of antimalarial treatment throughout the world. To identify key features associated with the delayed parasite clearance phenotype, we employed DNA microarrays to profile the physiological gene expression pattern of the resistant isolates.

RESULTS

In the ring and trophozoite stages, we observed reduced expression of many basic metabolic and cellular pathways which suggests a slower growth and maturation of these parasites during the first half of the asexual intraerythrocytic developmental cycle (IDC). In the schizont stage, there is an increased expression of essentially all functionalities associated with protein metabolism which indicates the prolonged and thus increased capacity of protein synthesis during the second half of the resistant parasite IDC. This modulation of the P. falciparum intraerythrocytic transcriptome may result from differential expression of regulatory proteins such as transcription factors or chromatin remodeling associated proteins. In addition, there is a unique and uniform copy number variation pattern in the Cambodian parasites which may represent an underlying genetic background that contributes to the resistance phenotype.

CONCLUSIONS

The decreased metabolic activities in the ring stages are consistent with previous suggestions of higher resilience of the early developmental stages to artemisinin. Moreover, the increased capacity of protein synthesis and protein turnover in the schizont stage may contribute to artemisinin resistance by counteracting the protein damage caused by the oxidative stress and/or protein alkylation effect of this drug. This study reports the first global transcriptional survey of artemisinin resistant parasites and provides insight to the complexities of the molecular basis of pathogens with drug resistance phenotypes in vivo.

Genomewide scan reveals amplification of mdr1 as a common denominator of resistance to mefloquine, lumefantrine, and artemisinin in Plasmodium chabaudi malaria parasites

Multidrug-resistant Plasmodium falciparum malaria parasites pose a threat to effective drug control, even to artemisinin-based combination therapies (ACTs). Here we used linkage group selection and Solexa whole-genome resequencing to investigate the genetic basis of resistance to component drugs of ACTs. Using the rodent malaria parasite P. chabaudi, we analyzed the uncloned progeny of a genetic backcross between the mefloquine-, lumefantrine-, and artemisinin-resistant mutant AS-15MF and a genetically distinct sensitive clone, AJ, following drug treatment. Genomewide scans of selection showed that parasites surviving each drug treatment bore a duplication of a segment of chromosome 12 (translocated to chromosome 04) present in AS-15MF. Whole-genome resequencing identified the size of the duplicated segment and its position on chromosome 4. The duplicated fragment extends for ∼393 kbp and contains over 100 genes, including mdr1, encoding the multidrug resistance P-glycoprotein homologue 1. We therefore show that resistance to chemically distinct components of ACTs is mediated by the same genetic mutation, highlighting a possible limitation of these therapies.

A research agenda for malaria eradication: basic science and enabling technologies

Today's malaria control efforts are limited by our incomplete understanding of the biology of Plasmodium and of the complex relationships between human populations and the multiple species of mosquito and parasite. Research priorities include the development of in vitro culture systems for the complete life cycle of P. falciparum and P. vivax and the development of an appropriate liver culture system to study hepatic stages. In addition, genetic technologies for the manipulation of Plasmodium need to be improved, the entire parasite metabolome needs to be characterized to identify new druggable targets, and improved information systems for monitoring the changes in epidemiology, pathology, and host-parasite-vector interactions as a result of intensified control need to be established to bridge the gap between bench, preclinical, clinical, and population-based sciences.

How can we identify parasite genes that underlie antimalarial drug resistance?

This article outlines genome-scale approaches that can be used to identify mutations in malaria (Plasmodium) parasites that underlie drug resistance and contribute to treatment failure. These approaches include genetic mapping by linkage or genome-wide association studies, drug selection and characterization of resistant mutants, and the identification of genome regions under strong recent selection. While these genomic approaches can identify candidate resistance loci, genetic manipulation is needed to demonstrate causality. We therefore also describe the growing arsenal of available transfection approaches for direct incrimination of mutations suspected to play a role in resistance. Our intention is both to review past progress and highlight promising approaches for future investigations.

Plasmodium sensitivity to artemisinins: magic bullets hit elusive targets

Artemisinins are efficacious antimalarial drugs widely employed as first-line treatment in endemic countries under the form of combined therapies. Different molecular modes of action have been postulated to explain the parasiticidal effect of these compounds; however, none has been unequivocally accepted, and their physiological relevance is still questioned. Similarly, no definite genetic determinant of Plasmodium sensitivity to artemisinins has been identified so far. A better understanding of the mode of action of artemisinins and the genetic basis of laboratory-induced or field-observed altered susceptibility is crucial for malaria control. In this review different models of artemisinins' molecular action are briefly presented, focusing on recent advances, and the evidence of potential association between various gene polymorphisms and artemisinin resistance is comprehensively reviewed.

Whole genome re-sequencing identifies a mutation in an ABC transporter (mdr2) in a Plasmodium chabaudi clone with altered susceptibility to antifolate drugs

In malaria parasites, mutations in two genes of folate biosynthesis encoding dihydrofolate reductase (dhfr) and dihydropteroate synthase (dhps) modify responses to antifolate therapies which target these enzymes. However, the involvement of other genes which modify the availability of exogenous folate, for example, has been proposed. Here, we used short-read whole-genome re-sequencing to determine the mutations in a clone of the rodent malaria parasite, Plasmodium chabaudi, which has altered susceptibility to both sulphadoxine and pyrimethamine. This clone bears a previously identified S106N mutation in dhfr and no mutation in dhps. Instead, three additional point mutations in genes on chromosomes 2, 13 and 14 were identified. The mutated gene on chromosome 13 (mdr2 K392Q) encodes an ABC transporter. Because Quantitative Trait Locus analysis previously indicated an association of genetic markers on chromosome 13 with responses to individual and combined antifolates, MDR2 is proposed to modulate antifolate responses, possibly mediated by the transport of folate intermediates.