Functional analysis of Plasmodium falciparum subpopulations associated with artemisinin resistance in Cambodia

Genomic analysis of Indian isolates of Plasmodium falciparum: Implications for drug resistance and virulence factors

The emergence of drug resistance to frontline treatments such as Artemisinin-based combination therapy (ACT) is a major obstacle to the control and eradication of malaria. This problem is compounded by the inherent genetic variability of the parasites, as many established markers of resistance do not accurately predict the drug-resistant status. There have been reports of declining effectiveness of ACT in the West Bengal and Northeast regions of India, which have traditionally been areas of drug resistance emergence in the country. Monitoring the genetic makeup of a population can help to identify the potential for drug resistance markers associated with it and evaluate the effectiveness of interventions aimed at reducing the spread of malaria. In this study, we performed whole genome sequencing of 53 isolates of Plasmodium falciparum from West Bengal and compared their genetic makeup to isolates from Southeast Asia (SEA) and Africa. We found that the Indian isolates had a distinct genetic makeup compared to those from SEA and Africa, and were more similar to African isolates, with a high prevalence of mutations associated with antigenic variation genes. The Indian isolates also showed a high prevalence of markers of chloroquine resistance (mutations in Pfcrt) and multidrug resistance (mutations in Pfmdr1), but no known mutations associated with artemisinin resistance in the PfKelch13 gene. Interestingly, we observed a novel L152V mutation in PfKelch13 gene and other novel mutations in genes involved in ubiquitination and vesicular transport that have been reported to support artemisinin resistance in the early stages of ACT resistance in the absence of PfKelch13 polymorphisms. Thus, our study highlights the importance of region-specific genomic surveillance for artemisinin resistance and the need for continued monitoring of resistance to artemisinin and its partner drugs.

Cytoprotective autophagy as a pro-survival strategy in ART-resistant malaria parasites

Despite several initiatives to subside the global malaria burden, the spread of artemisinin-resistant parasites poses a big threat to malaria elimination. Mutations in PfKelch13 are predictive of ART resistance, whose underpinning molecular mechanism remains obscure. Recently, endocytosis and stress response pathways such as the ubiquitin-proteasome machinery have been linked to artemisinin resistance. With Plasmodium, however, ambiguity persists regarding a role in ART resistance for another cellular stress defence mechanism called autophagy. Therefore, we investigated whether, in the absence of ART treatment, basal autophagy is augmented in PfK13-R539T mutant ART-resistant parasites and analyzed whether PfK13-R539T endowed mutant parasites with an ability to utilize autophagy as a pro-survival strategy. We report that in the absence of any ART treatment, PfK13-R539T mutant parasites exhibit increased basal autophagy compared to PfK13-WT parasites and respond aggressively through changes in autophagic flux. A clear cytoprotective role of autophagy in parasite resistance mechanism is evident by the observation that a suppression of PI3-Kinase (PI3K) activity (a master autophagy regulator) rendered difficulty in the survival of PfK13-R539T ART-resistant parasites. In conclusion, we now show that higher PI3P levels reported for mutant PfKelch13 backgrounds led to increased basal autophagy that acts as a pro-survival response to ART treatment. Our results highlight PfPI3K as a druggable target with the potential to re-sensitize ART-resistant parasites and identify autophagy as a pro-survival function that modulates ART-resistant parasite growth.

Nutrient Limitation Mimics Artemisinin Tolerance in Malaria

Mounting evidence demonstrates that nutritional environment can alter pathogen drug sensitivity. While the rich media used for in vitro culture contains supraphysiological nutrient concentrations, pathogens encounter a relatively restrictive environment in vivo. We assessed the effect of nutrient limitation on the protozoan parasite that causes malaria and demonstrated that short-term growth under physiologically relevant mild nutrient stress (or "metabolic priming") triggers increased tolerance of a potent antimalarial drug. We observed beneficial effects using both short-term survival assays and longer-term proliferation studies, where metabolic priming increases parasite survival to a level previously defined as resistant (>1% survival). We performed these assessments by either decreasing single nutrients that have distinct roles in metabolism or using a media formulation that simulates the human plasma environment. We determined that priming-induced tolerance was restricted to parasites that had newly invaded the host red blood cell, but the effect was not dependent on genetic background. The molecular mechanisms of this intrinsic effect mimic aspects of genetic tolerance, including translational repression and protein export. This finding suggests that regardless of the impact on survival rates, environmental stress could stimulate changes that ultimately directly contribute to drug tolerance. Because metabolic stress is likely to occur more frequently in vivo compared to the stable in vitro environment, priming-induced drug tolerance has ramifications for how in vitro results translate to in vivo studies. Improving our understanding of how pathogens adjust their metabolism to impact survival of current and future drugs is an important avenue of research to slow the evolution of resistance. IMPORTANCE There is a dire need for effective treatments against microbial pathogens. Yet, the continuing emergence of drug resistance necessitates a deeper knowledge of how pathogens respond to treatments. We have long appreciated the contribution of genetic evolution to drug resistance, but transient metabolic changes that arise in response to environmental factors are less recognized. Here, we demonstrate that short-term growth of malaria parasites in a nutrient-limiting environment triggers cellular changes that lead to better survival of drug treatment. We found that these strategies are similar to those employed by drug-tolerant parasites, which suggests that starvation "primes" parasites to survive and potentially evolve resistance. Since the environment of the human host is relatively nutrient restrictive compared to growth conditions in standard laboratory culture, this discovery highlights the important connections among nutrient levels, protective cellular pathways, and resistance evolution.

Identification of Co-Existing Mutations and Gene Expression Trends Associated With K13-Mediated Artemisinin Resistance in Plasmodium falciparum

Plasmodium falciparum infects millions and kills thousands of people annually the world over. With the emergence of artemisinin and/or multidrug resistant strains of the pathogen, it has become even more challenging to control and eliminate the disease. Multiomics studies of the parasite have started to provide a glimpse into the confounding genetics and mechanisms of artemisinin resistance and identified mutations in Kelch13 (K13) as a molecular marker of resistance. Over the years, thousands of genomes and transcriptomes of artemisinin-resistant/sensitive isolates have been documented, supplementing the search for new genes/pathways to target artemisinin-resistant isolates. This meta-analysis seeks to recap the genetic landscape and the transcriptional deregulation that demarcate artemisinin resistance in the field. To explore the genetic territory of artemisinin resistance, we use genomic single-nucleotide polymorphism (SNP) datasets from 2,517 isolates from 15 countries from the MalariaGEN Network (The Pf3K project, pilot data release 4, 2015) to dissect the prevalence, geographical distribution, and co-existing patterns of genetic markers associated with/enabling artemisinin resistance. We have identified several mutations which co-exist with the established markers of artemisinin resistance. Interestingly, K13-resistant parasites harbor α-ß hydrolase and putative HECT domain-containing protein genes with the maximum number of SNPs. We have also explored the multiple, publicly available transcriptomic datasets to identify genes from key biological pathways whose consistent deregulation may be contributing to the biology of resistant parasites. Surprisingly, glycolytic and pentose phosphate pathways were consistently downregulated in artemisinin-resistant parasites. Thus, this meta-analysis highlights the genetic and transcriptomic features of resistant parasites to propel further exploratory studies in the community to tackle artemisinin resistance.

Mutation in Plasmodium falciparum BTB/POZ domain of K13 protein confers artemisinin resistance

Partial artemisinin resistance, defined in patients as a delayed parasite clearance following artemisinin-based treatment, is conferred by non-synonymous mutations in the Kelch beta-propeller domain of the Plasmodium falciparum k13 (pfk13) gene. Here, we carried out in vitro selection over a one-year period on a West African P. falciparum strain isolated from Kolle (Mali) under a dose-escalating artemisinin regimen. After 18 cycles of sequential drug pressure, the selected parasites exhibited enhanced survival to dihydroartemisinin in the ring-stage survival assay (RSA^0-3h = 9.2%). Sanger and whole-genome sequence analyses identified the PfK13 P413A mutation, localized in the BTB/POZ domain, upstream of the propeller domain. This mutation was sufficient to confer in vitro artemisinin resistance when introduced into the PfK13 coding sequence of the parasite strain Dd2 by CRISPR/Cas9 gene editing. These results together with structural studies of the protein demonstrate that the propeller domain is not the sole in vitro mediator of PfK13-mediated artemisinin resistance, and highlight the importance of monitoring for mutations throughout PfK13.

Integration of population and functional genomics to understand mechanisms of artemisinin resistance in Plasmodium falciparum

Resistance to antimalarial drugs, and in particular to the artemisinin derivatives and their partner drugs, threatens recent progress toward regional malaria elimination and eventual global malaria eradication. Population-level studies utilizing whole-genome sequencing approaches have facilitated the identification of regions of the parasite genome associated with both clinical and in vitro drug-resistance phenotypes. However, the biological relevance of genes identified in these analyses and the establishment of a causal relationship between genotype and phenotype requires functional characterization. Here we examined data from population genomic and transcriptomic studies in the context of data generated from recent functional studies, using a new population genetic approach designed to identify potential favored mutations within the region of a selective sweep (iSAFE). We identified several genes functioning in pathways now known to be associated with artemisinin resistance that were supported in early population genomic studies, as well as potential new drug targets/pathways for further validation and consideration for treatment of artemisinin-resistant Plasmodium falciparum. In addition, we establish the utility of iSAFE in identifying positively-selected mutations in population genomic studies, potentially accelerating the time to functional validation of candidate genes.

In vitro growth competition experiments that suggest consequences of the substandard artemisinin epidemic that may be accelerating drug resistance in P. falciparum malaria

Over the past decade, artemisinin (ART)-combination therapies (ACTs) have shown declining efficacy within Southeast Asia (SEA). These resistance-like phenomena manifest as a delayed clearance phenotype (DCP) in some patients treated with ACTs. ACTs are currently the recommended treatment for P. falciparum infections by the World Health Organization (WHO), and they are our last line of defense to effectively treat all strains of malaria. Acceleration of antimicrobial resistance (AMR) is often theorized to be exacerbated by the use of subtherapeutic dosages of drugs ("substandard" drug), which for ACTs has been well documented over the last decade. Troublingly, in 2017, the WHO estimated that nearly 1 in 10 medical products tested in low- and middle-income countries failed to meet quality standards. We have developed a tissue culture-based approach for testing possible connections between substandard treatment and the spread of ACT resistant blood stage forms of P. falciparum. Via sequencing of pfk13, a molecular marker that is predictive for ART resistance (ARTR), we monitor competition of sensitive vs resistant strains over time and under various conditions and define conditions that favor emergence of ARTR parasites. Our findings help to define the conditions under which substandard drug treatments might favor the proliferation of mutant PfK13-mediated drug resistant strains over drug sensitive.

Optimization of parasite DNA enrichment approaches to generate whole genome sequencing data for Plasmodium falciparum from low parasitaemia samples

Background

Owing to the large amount of host DNA in clinical samples, generation of high-quality Plasmodium falciparum whole genome sequencing (WGS) data requires enrichment for parasite DNA. Enrichment is often achieved by leukocyte depletion of infected blood prior to storage. However, leukocyte depletion is difficult in low-resource settings and limits analysis to prospectively-collected samples. As a result, approaches such as selective whole genome amplification (sWGA) are being used to enrich for parasite DNA. However, sWGA has had limited success in generating reliable sequencing data from low parasitaemia samples. In this study, enzymatic digestion with MspJI prior to sWGA and whole genome sequencing was evaluated to determine whether this approach improved genome coverage compared to sWGA alone. The potential of sWGA to cause amplification bias in polyclonal infections was also examined.

Methods

DNA extracted from laboratory-created dried blood spots was treated with a modification-dependent restriction endonuclease, MspJI, and filtered via vacuum filtration. Samples were then selectively amplified using a previously reported sWGA protocol and subjected to WGS. Genome coverage statistics were compared between the optimized sWGA approach and the previously reported sWGA approach performed in parallel. Differential amplification by sWGA was assessed by comparing WGS data generated from lab-created mixtures of parasite isolates, from the same geographical region, generated with or without sWGA.

Results

MspJI digestion did not enrich for parasite DNA. Samples that underwent vacuum filtration (without MspJI digestion) prior to sWGA had the highest parasite DNA concentration and displayed greater genome coverage compared to MspJI + sWGA and sWGA alone, particularly for low parasitaemia samples. The optimized sWGA (filtration + sWGA) approach was successfully used to generate WGS data from 218 non-leukocyte depleted field samples from Malawi. Sequences from lab-created mixtures of parasites did not show evidence of differential amplification of parasite strains compared to directly sequenced samples.

Conclusion

This optimized sWGA approach is a reliable method to obtain WGS data from non-leukocyte depleted, low parasitaemia samples. The absence of amplification bias in data generated from mixtures of isolates from the same geographic region suggests that this approach can be appropriately used for molecular epidemiological studies.

Origin and Spread of Evolving Artemisinin-Resistant Plasmodium falciparum Malarial Parasites in Southeast Asia

In this review, we provide an epidemiological history of the emergence and ongoing spread of evolving Plasmodium falciparum artemisinin resistance (ARTR). Southeast Asia has been the focal point for emergence and spread of multiple antimalarial drug resistance phenomena, and is once again for evolving ARTR, also known as the "delayed clearance phenotype" (DCP). The five countries most impacted, Cambodia, Thailand, Myanmar, Laos, and Vietnam, each have complex histories of antimalarial drug use over many decades, which have in part molded the use of various artemisinin combination therapies (ACTs) within each country. We catalog the use of ACTs, evolving loss of ACT efficacy, and the frequency of pfk13 mutations (mutations associated with ARTR) in the Greater Mekong Subregion and map the historical spread of ARTR/DCP parasites. These data should assist improved surveillance and deployment of next-generation ACTs.

Strains used in whole organism Plasmodium falciparum vaccine trials differ in genome structure, sequence, and immunogenic potential

BACKGROUND

Plasmodium falciparum (Pf) whole-organism sporozoite vaccines have been shown to provide significant protection against controlled human malaria infection (CHMI) in clinical trials. Initial CHMI studies showed significantly higher durable protection against homologous than heterologous strains, suggesting the presence of strain-specific vaccine-induced protection. However, interpretation of these results and understanding of their relevance to vaccine efficacy have been hampered by the lack of knowledge on genetic differences between vaccine and CHMI strains, and how these strains are related to parasites in malaria endemic regions.

METHODS

Whole genome sequencing using long-read (Pacific Biosciences) and short-read (Illumina) sequencing platforms was conducted to generate de novo genome assemblies for the vaccine strain, NF54, and for strains used in heterologous CHMI (7G8 from Brazil, NF166.C8 from Guinea, and NF135.C10 from Cambodia). The assemblies were used to characterize sequences in each strain relative to the reference 3D7 (a clone of NF54) genome. Strains were compared to each other and to a collection of clinical isolates (sequenced as part of this study or from public repositories) from South America, sub-Saharan Africa, and Southeast Asia.

RESULTS

While few variants were detected between 3D7 and NF54, we identified tens of thousands of variants between NF54 and the three heterologous strains. These variants include SNPs, indels, and small structural variants that fall in regulatory and immunologically important regions, including transcription factors (such as PfAP2-L and PfAP2-G) and pre-erythrocytic antigens that may be key for sporozoite vaccine-induced protection. Additionally, these variants directly contributed to diversity in immunologically important regions of the genomes as detected through in silico CD8+ T cell epitope predictions. Of all heterologous strains, NF135.C10 had the highest number of unique predicted epitope sequences when compared to NF54. Comparison to global clinical isolates revealed that these four strains are representative of their geographic origin despite long-term culture adaptation; of note, NF135.C10 is from an admixed population, and not part of recently formed subpopulations resistant to artemisinin-based therapies present in the Greater Mekong Sub-region.

CONCLUSIONS

These results will assist in the interpretation of vaccine efficacy of whole-organism vaccines against homologous and heterologous CHMI.

Pre-clinical study of iron oxide nanoparticles fortified artesunate for efficient targeting of malarial parasite

BACKGROUND

Artesunate the most potent antimalarial is widely used for the treatment of multidrug-resistant malaria. The antimalarial cytotoxicity of artesunate has been mainly attributed to its selective, irreversible and iron- radical-mediated damage of parasite biomolecules. In the present research, iron oxide nanoparticle fortified artesunate was tested in P. falciparum and in an experimental malaria mouse model for enhancement in the selectivity and toxicity of artesunate towards parasite. Artesunate was fortified with nontoxic biocompatible surface modified iron oxide nanoparticle which is specially designed and synthesized for the sustained pH-dependent release of Fe2+ within the parasitic food vacuole for enhanced ROS spurt.

METHODS

Antimalarial efficacy of Iron oxide nanoparticle fortified artesunate was evaluated in wild type and artemisinin-resistant Plasmodium falciparum (R539T) grown in O + ve human blood and in Plasmodium berghei ANKA infected swiss albino mice. Internalization of nanoparticles, the pH-dependent release of Fe2+, production of reactive oxygen species and parasite biomolecule damage by iron oxide nanoparticle fortified artesunate was studied using various biochemical, biophysical, ultra-structural and fluorescence microscopy. For determining the efficacy of ATA-IONP+ART on resistant parasite ring survival assay was performed.

RESULTS

The nanoparticle fortified artesunate was highly efficient in the 1/8th concentration of artesunate IC50 and led to retarded growth of P. falciparum with significant damage to macromolecules mediated via enhanced ROS production. Similarly, preclinical In vivo studies also signified a radical reduction in parasitemia with ~8-10-fold reduced dosage of artesunate when fortified with iron oxide nanoparticles. Importantly, the ATA-IONP combination was efficacious against artemisinin-resistant parasites.

INTERPRETATION

Surface coated iron-oxide nanoparticle fortified artesunate can be developed into a potent therapeutic agent towards multidrug-resistant and artemisinin-resistant malaria in humans. FUND: This study is supported by the Centre for Study of Complex Malaria in India funded by the National Institute of Health, USA.

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.

Dihydroartemisinin-Ferriprotoporphyrin IX Adduct Abundance in Plasmodium falciparum Malarial Parasites and the Relationship to Emerging Artemisinin Resistance

Previously (Heller, L. E., and Roepe, P. D. Quantification of Free Ferriprotoporphyrin IX Heme and Hemozoin for Artemisinin Sensitive versus Delayed Clearance Phenotype Plasmodium falciparum Malarial Parasites. Biochemistry, DOI: 10.1021/acs.biochem.8b00959, preceding paper in this issue), we quantified free ferriprotoporphyrin IX (FPIX) heme abundance for control versus delayed clearance phenotype (DCP) intraerythrocytic Plasmodium falciparum malarial parasites. Because artemisinin drugs are activated by free FPIX, these data predict that the abundance of long-hypothesized toxic artemisinin drug-FPIX covalent adducts might differ for control versus DCP parasites. If so, this would have important repercussions for understanding the mechanism of the DCP, also known as emerging artemisinin resistance. To test these predictions, we studied in vitro formation of FPIX-dihydroartemisinin (DHA) adducts and then for the first time quantified the abundance of FPIX-DHA adducts formed within live P. falciparum versus the stage of intraerythrocytic development. Using matched isogenic parasite strains, we quantified the adduct for DCP versus control parasite strains and found that mutant PfK13 mediates lower adduct abundance for DCP parasites. The results suggest improved models for the molecular pharmacology of artemisinin-based antimalarial drugs and the molecular mechanism of the DCP.

Mechanisms of artemisinin resistance in Plasmodium falciparum malaria

Artemisinin-based combination therapies (ACTs) have substantially reduced worldwide malaria burden and deaths. But malaria parasites have become resistant to artemisinins. Prior studies suggested two different molecular pathways of artemisinin-resistance. Here we unify recent findings into a single model, where elevation of a lipid, phosphatidylinositol-3- phosphate (PI3P) results in vesicle expansion that increases the engagement with the unfolded protein response (UPR). Vesicle expansion (rather than increasing individual genetic determinants of the UPR) efficiently induces artemisinin resistance likely by promoting ‘proteostasis’ (protein translation coupled to proper protein folding and vesicular remodeling) to mitigate artemisinin-induced proteopathy (death from global abnormal protein-toxicity). Vesicular amplification engages the host red cell, suggesting that artemisinin resistant malaria may also persist by taking advantage of host niches and escaping the immune response.