Genetic loci associated with delayed clearance of Plasmodium falciparum following artemisinin treatment in Southeast Asia

A Biotinylated cpFIT-PNA Platform for the Facile Detection of Drug Resistance to Artemisinin in Plasmodium falciparum

The evolution of drug resistance to many antimalarial drugs in the lethal strain of malaria (Plasmodium falciparum) has been a great concern over the past 50 years. Among these drugs, artemisinin has become less effective for treating malaria. Indeed, several P. falciparum variants have become resistant to this drug, as elucidated by specific mutations in the pfK13 gene. This study presents the development of a diagnostic kit for the detection of a common point mutation in the pfK13 gene of P. falciparum, namely, the C580Y point mutation. FIT-PNAs (forced-intercalation peptide nucleic acid) are DNA mimics that serve as RNA sensors that fluoresce upon hybridization to their complementary RNA. Herein, FIT-PNAs were designed to sense the C580Y single nucleotide polymorphism (SNP) and were conjugated to biotin in order to bind these molecules to streptavidin-coated plates. Initial studies with synthetic RNA were conducted to optimize the sensing system. In addition, cyclopentane-modified PNA monomers (cpPNAs) were introduced to improve FIT-PNA sensing. Lastly, total RNA was isolated from red blood cells infected with P. falciparum (WT strain - NF54-WT or mutant strain - NF54-C580Y). Streptavidin plates loaded with either FIT-PNA or cpFIT-PNA were incubated with the total RNA. A significant difference in fluorescence for mutant vs WT total RNA was found only for the cpFIT-PNA probe. In summary, this study paves the way for a simple diagnostic kit for monitoring artemisinin drug resistance that may be easily adapted to malaria endemic regions.

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.

Increase of Plasmodium falciparum parasites carrying lumefantrine-tolerance molecular markers and lack of South East Asian pfk13 artemisinin-resistance mutations in samples collected from 2013 to 2016 in Côte d'Ivoire

One of the major obstacles to malaria elimination in the world is the resistance in Plasmodium falciparum to most antimalarial drugs. This study aimed to estimate the prevalence of molecular markers of antimalarial drugs resistance in Côte d'Ivoire. Samples were collected from 2013 to 2016 from asymptomatic and symptomatic subjects in Abengourou, Abidjan, Grand Bassam, and San Pedro. A total of 704 participants aged between 1 year and 65 years (Mean age: 9 years ± 7.7) were enrolled. All the dried filter paper blood spots were genotyped by sequencing. Plasmodium falciparum kelch propeller domain 13 (pfk13) gene were analyzed for all the samples, while 344 samples were examined for Plasmodium falciparum multi-drug resistance 1 (pfmdr1). Overall, the success rate of molecular tests was 98.8% (340/344), 99.1% (341/344), and 94.3% (664/704) for pfmdr1 N86Y, pfmdr1 Y184F, and pfk13 genes respectively. Molecular analysis revealed twenty (5.9%; 20/340) and 219 (64.2%; 219/341) mutant alleles for pfmdr1 86Y and pfmdr1 184 F, respectively. Twenty-nine mutations in pfk13 gene (4.4%; 29/664) with 2.7% (18/664) of non-synonymous mutations was found. None of the mutations previously described in South East Asia (SEA) involved in P. falciparum resistance to artemisinin derivatives were observed in this study. According to year of collection, a decrease of the prevalence of pfk13 mutation (from 3.6 to 1.8%) and pfmdr1 N86Y mutation (from 8.5 to 4.5%) and an increase of mutant allele of pfmdr1 Y184F proportion (from 39.8 to 66.4%) were found. Comparing to previous studies in the country, this study showed an increase in lumefantrine tolerance of P. falciparum strains. This demonstrates the importance of establishing a strong system for molecular surveillance of malaria in Côte d'Ivoire.

Fused Enzyme Glucose-6-Phosphate Dehydrogenase::6-Phosphogluconolactonase (G6PD::6PGL) as a Potential Drug Target in Giardia lamblia, Trichomonas vaginalis, and Plasmodium falciparum

Several microaerophilic parasites such as Giardia lamblia, Trichomonas vaginalis, and Plasmodium falciparum are major disease-causing organisms and are responsible for spreading infections worldwide. Despite significant progress made in understanding the metabolism and molecular biology of microaerophilic parasites, chemotherapeutic treatment to control it has seen limited progress. A current proposed strategy for drug discovery against parasitic diseases is the identification of essential key enzymes of metabolic pathways associated with the parasite's survival. In these organisms, glucose-6-phosphate dehydrogenase::6-phosphogluconolactonase (G6PD:: 6PGL), the first enzyme of the pentose phosphate pathway (PPP), is essential for its metabolism. Since G6PD:: 6PGL provides substrates for nucleotides synthesis and NADPH as a source of reducing equivalents, it could be considered an anti-parasite drug target. This review analyzes the anaerobic energy metabolism of G. lamblia, T. vaginalis, and P. falciparum, with a focus on glucose metabolism through the pentose phosphate pathway and the significance of the fused G6PD:: 6PGL enzyme as a therapeutic target in the search for new drugs.

Genomic complexity of parasites and vectors challenges malaria control in Southeast Asia

Southeast Asia is a uniquely complex region of malaria transmission that maintains an astounding level of species diversity among potential malaria vectors and also generates drug-resistant and quickly diverging populations of malaria parasites. All five human malaria species circulate in Southeast Asia with over 50 Anopheles species that vary in their ability to transmit these pathogens. The intricate relationships of these parasites and vectors are not well-understood. Human activity in Southeast Asian countries has created an increasingly fragmented landscape, bringing humans and mosquitoes into more frequent contact, sustaining malaria transmission in a region where few control tools are effective. Genomic shifts at the species, population, and individual level in parasites and vectors introduce variation that has produced drug- and insecticide resistance. The goal of this review is to highlight genomic studies of Southeast Asian malaria parasites and vectors that demonstrate how diversity in these organisms presents unique challenges and opportunities for global malaria control and eradication efforts.

Artemisinin resistance in P. falciparum: probing the interacting partners of Kelch13 protein in parasite

OBJECTIVES

Artemisinin (ART) resistance in Plasmodium is threatening the artemisinin combination therapies-the first line of defence against malaria. ART resistance has been established to be mediated by the Plasmodium Kelch13 (PfK13) protein. For the crucial role of PfK13 in multiple pathways of the Plasmodium life cycle and ART resistance, it is imperative that we investigate its interacting partners.

METHODS

We recombinantly expressed PfK13-p (Bric a brac/Poxvirus and zinc finger and propeller domains), generating anti-PfK13-p antibodies to perform co-immunoprecipitation assays and probed PfK13 interacting partners. Surface plasmon resonance and pull-down assays were performed to establish physical interactions of representative proteins with PfK13-p.

RESULTS

The co-immunoprecipitation assays identified 17 proteins with distinct functions in the parasite life cycle- protein folding, cellular metabolism, and protein binding and invasion. In addition to the overlap with previously identified proteins, our study identified 10 unique proteins. Fructose-biphosphate aldolase and heat shock protein 70 demonstrated strong biophysical interaction with PfK13-p, with KD values of 6.6 µM and 7.6 µM, respectively. Additionally, Plasmodium merozoite surface protein 1 formed a complex with PfK13-p, which is evident from the pull-down assay.

CONCLUSION

This study adds to our knowledge of the PfK13 protein in mediating ART resistance by identifying new PfK13 interacting partners. Three representative proteins-fructose-biphosphate aldolase, heat shock protein 70, and merozoite surface protein 1-demonstrated clear evidence of biophysical interactions with PfK13-p. However, elucidation of the functional relevance of these physical interactions are crucial in context of PfK13 role in ART resistance.

Molecular surveillance of Plasmodium falciparum drug-resistance markers in Vietnam using multiplex amplicon sequencing (2000-2016)

Emergence and spread of Plasmodium falciparum resistance to artemisinin-based combination therapies (ACT) is a major challenge for Greater Mekong Subregion countries in their goal to eliminate malaria by 2030. Tools to efficiently monitor drug resistance beyond resource-demanding therapeutic efficacy studies are necessary. A custom multiplex amplicon sequencing assay based on Illumina technology was designed to target the marker of partial resistance to artemisinin (K13), five candidate modulators of artemisinin resistance, the marker of resistance to chloroquine (crt), and four neutral microsatellite loci. The assay was used to genotype 635 P. falciparum-positive blood samples collected across seven provinces of Vietnam and one of Cambodia between 2000 and 2016. Markers of resistance to artemisinin partner-drugs piperaquine (copy number of plasmepsin-2) and mefloquine (copy number of multidrug-resistance 1) were determined by qPCR. Parasite population structure was further assessed using a 101-SNP barcode. Validated mutations of artemisinin partial resistance in K13 were found in 48.1% of samples, first detection was in 2000, and by 2015 prevalence overcame > 50% in Central Highlands and Binh Phuoc province. K13-C580Y variant became predominant country-wide, quickly replacing an outbreak of K13-I543T in Central Highlands. Mutations in candidate artemisinin resistance modulator genes paralleled the trends of K13 mutants, whereas resistance to piperaquine and mefloquine remained low (≈ 10%) by 2015-2016. Genomic tools applied to malaria surveillance generate comprehensive information on dynamics of drug resistance and population structure and reflect drug efficacy profiles from in vivo studies.

Malaria Genomics, Vaccine Development, and Microbiome

Recent advances in malaria genetics and genomics have transformed many aspects of malaria research in areas of molecular evolution, epidemiology, transmission, host-parasite interaction, drug resistance, pathogenicity, and vaccine development. Here, in addition to introducing some background information on malaria parasite biology, parasite genetics/genomics, and genotyping methods, we discuss some applications of genetic and genomic approaches in vaccine development and in studying interactions with microbiota. Genetic and genomic data can be used to search for novel vaccine targets, design an effective vaccine strategy, identify protective antigens in a whole-organism vaccine, and evaluate the efficacy of a vaccine. Microbiota has been shown to influence disease outcomes and vaccine efficacy; studying the effects of microbiota in pathogenicity and immunity may provide information for disease control. Malaria genetics and genomics will continue to contribute greatly to many fields of malaria research.

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.

Toxins from Animal Venoms as a Potential Source of Antimalarials: A Comprehensive Review

Malaria is an infectious disease caused by Plasmodium spp. and it is mainly transmitted to humans by female mosquitoes of the genus Anopheles. Malaria is an important global public health problem due to its high rates of morbidity and mortality. At present, drug therapies and vector control with insecticides are respectively the most commonly used methods for the treatment and control of malaria. However, several studies have shown the resistance of Plasmodium to drugs that are recommended for the treatment of malaria. In view of this, it is necessary to carry out studies to discover new antimalarial molecules as lead compounds for the development of new medicines. In this sense, in the last few decades, animal venoms have attracted attention as a potential source for new antimalarial molecules. Therefore, the aim of this review was to summarize animal venom toxins with antimalarial activity found in the literature. From this research, 50 isolated substances, 4 venom fractions and 7 venom extracts from animals such as anurans, spiders, scorpions, snakes, and bees were identified. These toxins act as inhibitors at different key points in the biological cycle of Plasmodium and may be important in the context of the resistance of Plasmodium to currently available antimalarial drugs.

Population-based genomic study of Plasmodium vivax malaria in seven Brazilian states and across South America

Background

Brazil is a unique and understudied setting for malaria, with complex foci of transmission associated with human and environmental conditions. An understanding of the population genomic diversity of P. vivax parasites across Brazil can support malaria control strategies.

Methods

Through whole genome sequencing of P. vivax isolates across 7 Brazilian states, we use population genomic approaches to compare genetic diversity within country (n = 123), continent (6 countries, n = 315) and globally (26 countries, n = 885).

Findings

We confirm that South American isolates are distinct, have more ancestral populations than the other global regions, with differentiating mutations in genes under selective pressure linked to antimalarial drugs (pvmdr1, pvdhfr-ts) and mosquito vectors (pvcrmp3, pvP45/48, pvP47). We demonstrate Brazil as a distinct parasite population, with signals of selection including ABC transporter (PvABCI3) and PHIST exported proteins.

Interpretation

Brazil has a complex population structure, with evidence of P. simium infections and Amazonian parasites separating into multiple clusters. Overall, our work provides the first Brazil-wide analysis of P. vivax population structure and identifies important mutations, which can inform future research and control measures.

Funding

AI is funded by an MRC LiD PhD studentship. TGC is funded by the Medical Research Council (Grant no. MR/M01360X/1, MR/N010469/1, MR/R025576/1, MR/R020973/1 and MR/X005895/1). SC is funded by Medical Research Council UK grants (MR/M01360X/1, MR/R025576/1, MR/R020973/1 and MR/X005895/1) and Bloomsbury SET (ref. CCF17-7779). FN is funded by The Shloklo Malaria Research Unit - part of the Mahidol Oxford Research Unit, supported by the Wellcome Trust (Grant no. 220211). ARSB is funded by São Paulo Research Foundation - FAPESP (Grant no. 2002/09546-1). RLDM is funded by Brazilian National Council for Scientific and Technological Development - CNPq (Grant no. 302353/2003-8 and 471605/2011-5); CRFM is funded by FAPESP (Grant no. 2020/06747-4) and CNPq (Grant no. 302917/2019-5 and 408636/2018-1); JGD is funded by FAPESP fellowships (2016/13465-0 and 2019/12068-5) and CNPq (Grant no. 409216/2018-6).

Molecular detection of drug-resistant Plasmodium falciparum mutants in Assam

Background & objectives

The spread of drug-resistant Plasmodium falciparum (Pf) poses a serious threat to the control and elimination of malaria. The objective of this study was to detect the molecular biomarkers of antimalarial drug resistance in Pf in patients visiting a tertiary care hospital in Assam.

Methods

Malaria was first detected in fever cases using microscopy and a rapid diagnostic test (RDT), and then confirmed using PCR. Pf chloroquine resistance transporter (Pfcrt), Pf multidrug resistance-1 (Pfmdr-1), and single-nucleotide polymorphisms linked to delayed parasite clearance after treatment with artemisinin MAL 10-688956 and MAL 13-1718319 and Kelch-13 propeller (PfK-13) genes were evaluated by PCR-restriction fragment length polymorphism (RFLP).

Results

Sixty nine cases of malaria were found among 300 cases of fever. Of these, 54 were positive for Pf, 47 of which were confirmed by PCR. Pfcrt-K76T mutation was seen in 96.6 per cent and Pfmdr1-N86Y mutation in 84.2 per cent of cases. Mutation was not detected in MAL10 and MAL13 genes. Sequence analysis of Kelch-13 gene showed the presence of a novel mutation at amino acid position 675. Statistically, no significant association was found between the molecular biomarkers and demographic profile, clinical presentation and outcome of the cases.

Interpretation & conclusions

Molecular surveillance is essential to assess the therapeutic efficacy of the drugs against circulating Pf isolates in Assam which are found to be highly resistant to CQ. The role of the new mutation found in the Kelch-13 gene in the development of artemisinin resistance in Assam needs to be thoroughly monitored in future research.

Artemisinin resistance in the malaria parasite, Plasmodium falciparum, originates from its initial transcriptional response

The emergence and spread of artemisinin-resistant Plasmodium falciparum, first in the Greater Mekong Subregion (GMS), and now in East Africa, is a major threat to global malaria elimination ambitions. To investigate the artemisinin resistance mechanism, transcriptome analysis was conducted of 577 P. falciparum isolates collected in the GMS between 2016–2018. A specific artemisinin resistance-associated transcriptional profile was identified that involves a broad but discrete set of biological functions related to proteotoxic stress, host cytoplasm remodelling, and REDOX metabolism. The artemisinin resistance-associated transcriptional profile evolved from initial transcriptional responses of susceptible parasites to artemisinin. The genetic basis for this adapted response is likely to be complex.

Transcriptomic analysis of isolates from the malaria parasite (Plasmodium falciparum) in the Greater Mekong Subregion of Southeast Asia identifies gene expression patterns that are correlated with resistance to a common anti-malaria drug, artemisinin.

Reduction in Plasmodium falciparum Pfk13 and pfg377 allele diversity through time in southern Vietnam

Background

Plasmodium falciparum has acquired resistance to artemisinin in Southeast Asia, with mutations in the P. falciparum Kelch-13 (Pfk13) gene associated with the resistance phenotype. The widespread use of Artemisinin-based combination therapy (ACT)s in Southeast Asia has led to the selection and spread of parasites carrying mutations in Pfk13. We characterised the allele diversity of Pfk13 and pfg377, an artemisinin-resistance neutral polymorphic gene, in parasite DNA extracted human blood from in southern Vietnam in 2003, 2012, 2015 and 2018.

Method

This study was conducted in Bu Gia Map commune, Binh Phuoc province, Vietnam, from May 2018 to January 2019. Twenty-four samples from 2018 to 2019, 30 from 2003, 24 from 2012 and 32 from 2015 were analysed. Malaria-infected human blood was collected by finger-prick and used for molecular analysis. A nested-PCR targeting the small subunit ribosomal RNA gene was used for Plasmodium species identification, followed by amplification and nucleotide sequencing of Pfk13 and region 3 of pfg377. Archived blood samples collected in the same region in 2012 and 2015 were also analysed as above for comparison.

Results

The genetic diversity of Pfk13 and pfg377 was lower in 2018–2019 compared to 2012 and 2015. The number of distinct Pfk13 mutants decreased from three in 2012 and 2015, P553L, V568G and C580Y, to one, C580Y in 2018–2019. In 2018–2019, the frequency of C580Y mutant strains was 71% (17/24 isolates). All samples were wild type in 2003. In 2012 and 2015, there were single-strain infections as well as co-infections with two mutant strains or with mutant and wild strains, whereas there were no co-infections in 2018. pfg377 allele diversity decreased from five alleles in 2012 to two alleles in 2018–2019.

Conclusion

The genetic diversity of P. falciparum was reduced at the two genetic loci surveyed in this study, Pfk13 and pfg377. In the case of the former gene, we observed an increase in the prevalence of parasites carrying the C580Y gene, known to confer reduced susceptibility to ACTs. The reduction in the diversity of pfg377 may be linked to the clonal expansion of parasite strains carrying the C580Y mutation, leading to an overall reduction in parasite genetic diversity across the population.

Supplementary Information

The online version contains supplementary material available at 10.1186/s41182-022-00409-4.

Machine learning for artemisinin resistance in malaria treatment across in vivo-in vitro platforms

Drug resistance has been rapidly evolving with regard to the first-line malaria treatment, artemisinin-based combination therapies. It has been an open question whether predictive models for this drug resistance status can be generalized across in vivo-in vitro transcriptomic measurements. In this study, we present a model that predicts artemisinin treatment resistance developed with transcriptomic information of Plasmodium falciparum. We demonstrated the robustness of this model across in vivo clearance rate and in vitro IC50 measurement and based on different microarray and data processing modalities. The validity of the algorithm is further supported by its first placement in the DREAM Malaria challenge. We identified transcription biomarkers to artemisinin treatment resistance that can predict artemisinin resistance and are conserved in their expression modules. This is a critical step in the research of malaria treatment, as it demonstrated the potential of a platform-robust, personalized model for artemisinin resistance using molecular biomarkers.

•

Artemisinin resistance can be predicted from transcriptomes by machine learning

•

Our model can be transferred between in vivo and in vitro and different platforms

•

We identified top transcription biomarkers of artemisinin resistance

Impact of Sickle Cell Trait Hemoglobin on the Intraerythrocytic Transcriptional Program of Plasmodium falciparum

Sickle-trait hemoglobin (HbAS) confers nearly complete protection from severe, life-threatening falciparum malaria in African children. Despite this clear protection, the molecular mechanisms by which HbAS confers these protective phenotypes remain incompletely understood. As a forward genetic screen for aberrant parasite transcriptional responses associated with parasite neutralization in HbAS red blood cells (RBCs), we performed comparative transcriptomic analyses of Plasmodium falciparum in normal (HbAA) and HbAS erythrocytes during both in vitro cultivation of reference parasite strains and naturally occurring P. falciparum infections in Malian children with HbAA or HbAS. During in vitro cultivation, parasites matured normally in HbAS RBCs, and the temporal expression was largely unperturbed of the highly ordered transcriptional program that underlies the parasite’s maturation throughout the intraerythrocytic development cycle (IDC). However, differential expression analysis identified hundreds of transcripts aberrantly expressed in HbAS, largely occurring late in the IDC. Surprisingly, transcripts encoding members of the Maurer’s clefts were overexpressed in HbAS despite impaired parasite protein export in these RBCs, while parasites in HbAS RBCs underexpressed transcripts associated with the endoplasmic reticulum and those encoding serine repeat antigen proteases that promote parasite egress. Analyses of P. falciparum transcriptomes from 32 children with uncomplicated malaria identified stage-specific differential expression: among infections composed of ring-stage parasites, only cyclophilin 19B was underexpressed in children with HbAS, while trophozoite-stage infections identified a range of differentially expressed transcripts, including downregulation in HbAS of several transcripts associated with severe malaria in collateral studies. Collectively, our comparative transcriptomic screen in vitro and in vivo indicates that P. falciparum adapts to HbAS by altering its protein chaperone and folding machinery, oxidative stress response, and protein export machinery. Because HbAS consistently protects from severe P. falciparum, modulation of these responses may offer avenues by which to neutralize P. falciparum parasites.

IMPORTANCE Sickle-trait hemoglobin (HbAS) confers nearly complete protection from severe, life-threatening malaria, yet the molecular mechanisms that underlie HbAS protection from severe malaria remain incompletely understood. Here, we used transcriptome sequencing (RNA-seq) to measure the impact of HbAS on the blood-stage transcriptome of Plasmodium falciparum in in vitro time series experiments and in vivo samples from natural infections. Our in vitro time series data reveal that, during its blood stage, P. falciparum’s gene expression in HbAS is impacted primarily through alterations in the abundance of gene products as opposed to variations in the timing of gene expression. Collectively, our in vitro and in vivo data indicate that P. falciparum adapts to HbAS by altering its protein chaperone and folding machinery, oxidative stress response, and protein export machinery. Due to the persistent association of HbAS and protection from severe disease, these processes that are modified in HbAS may offer strategies to neutralize P. falciparum.

Single-cell genome sequencing of protozoan parasites

Despite considerable genetic variation within hosts, most parasite genome sequencing studies focus on bulk samples composed of millions of cells. Analysis of bulk samples is biased toward the dominant genotype, concealing cell-to-cell variation and rare variants. To tackle this, single-cell sequencing approaches have been developed and tailored to specific host-parasite systems. These are allowing the genetic diversity and kinship in complex parasite populations to be deciphered and for de novo genetic variation to be captured. Here, we outline the methodologies being used for single-cell sequencing of parasitic protozoans, such as Plasmodium and Leishmania spp., and how these tools are being applied to understand parasite biology.

Plasmodium falciparum resistance to ACTs: Emergence, mechanisms, and outlook

Emergence and spread of resistance in Plasmodium falciparum to the frontline treatment artemisinin-based combination therapies (ACTs) in the epicenter of multidrug resistance of Southeast Asia threaten global malaria control and elimination. Artemisinin (ART) resistance (or tolerance) is defined clinically as delayed parasite clearance after treatment with an ART drug. The resistance phenotype is restricted to the early ring stage and can be measured in vitro using a ring-stage survival assay. ART resistance is associated with mutations in the propeller domain of the Kelch family protein K13. As a pro-drug, ART is activated primarily by heme, which is mainly derived from hemoglobin digestion in the food vacuole. Activated ARTs can react promiscuously with a wide range of cellular targets, disrupting cellular protein homeostasis. Consistent with this mode of action for ARTs, the molecular mechanisms of K13-mediated ART resistance involve reduced hemoglobin uptake/digestion and increased cellular stress response. Mutations in other genes such as AP-2μ (adaptor protein-2 μ subunit), UBP-1 (ubiquitin-binding protein-1), and Falcipain 2a that interfere with hemoglobin uptake and digestion also increase resistance to ARTs. ART resistance has facilitated the development of resistance to the partner drugs, resulting in rapidly declining ACT efficacies. The molecular markers for resistance to the partner drugs are mostly associated with point mutations in the two food vacuole membrane transporters PfCRT and PfMDR1, and amplification of pfmdr1 and the two aspartic protease genes plasmepsin 2 and 3. It has been observed that mutations in these genes can have opposing effects on sensitivities to different partner drugs, which serve as the principle for designing triple ACTs and drug rotation. Although clinical ACT resistance is restricted to Southeast Asia, surveillance for drug resistance using in vivo clinical efficacy, in vitro assays, and molecular approaches is required to prevent or slow down the spread of resistant parasites.

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.

The newly discovered role of endocytosis in artemisinin resistance

Artemisinin and its derivatives (ART) are the cornerstone of malaria treatment as part of artemisinin combination therapy (ACT). However, reduced susceptibility to artemisinin as well as its partner drugs threatens the usefulness of ACTs. Single point mutations in the parasite protein Kelch13 (K13) are necessary and sufficient for the reduced sensitivity of malaria parasites to ART but several alternative mechanisms for this resistance have been proposed. Recent work found that K13 is involved in the endocytosis of host cell cytosol and indicated that this is the process responsible for resistance in parasites with mutated K13. These studies also identified a series of further proteins that act together with K13 in the same pathway, including previously suspected resistance proteins such as UBP1 and AP-2μ. Here, we give a brief overview of artemisinin resistance, present the recent evidence of the role of endocytosis in ART resistance and discuss previous hypotheses in light of this new evidence. We also give an outlook on how the new insights might affect future research.

Genomic Approaches to Drug Resistance in Malaria

Although the last two decades have seen a substantial decline in malaria incidence and mortality due to the use of insecticide-treated bed nets and artemisinin combination therapy, the threat of drug resistance is a constant obstacle to sustainable malaria control. Given that patients can die quickly from this disease, public health officials and doctors need to understand whether drug resistance exists in the parasite population, as well as how prevalent it is so they can make informed decisions about treatment. As testing for drug efficacy before providing treatment to malaria patients is impractical, researchers need molecular markers of resistance that can be more readily tracked in parasite populations. To this end, much work has been done to unravel the genetic underpinnings of drug resistance in Plasmodium falciparum. The aim of this review is to provide a broad overview of common genomic approaches that have been used to discover the alleles that drive drug response phenotypes in the most lethal human malaria parasite.

Molecular Mechanisms of Drug Resistance in Plasmodium falciparum Malaria

Understanding and controlling the spread of antimalarial resistance, particularly to artemisinin and its partner drugs, is a top priority. Plasmodium falciparum parasites resistant to chloroquine, amodiaquine, or piperaquine harbor mutations in the P. falciparum chloroquine resistance transporter (PfCRT), a transporter resident on the digestive vacuole membrane that in its variant forms can transport these weak-base 4-aminoquinoline drugs out of this acidic organelle, thus preventing these drugs from binding heme and inhibiting its detoxification. The structure of PfCRT, solved by cryogenic electron microscopy, shows mutations surrounding an electronegative central drug-binding cavity where they presumably interact with drugs and natural substrates to control transport. P. falciparum susceptibility to heme-binding antimalarials is also modulated by overexpression or mutations in the digestive vacuole membrane-bound ABC transporter PfMDR1 (P. falciparum multidrug resistance 1 transporter). Artemisinin resistance is primarily mediated by mutations in P. falciparum Kelch13 protein (K13), a protein involved in multiple intracellular processes including endocytosis of hemoglobin, which is required for parasite growth and artemisinin activation. Combating drug-resistant malaria urgently requires the development of new antimalarial drugs with novel modes of action.

Plasmodium falciparum kelch 13 Mutations, 9 Countries in Africa, 2014-2018

The spread of drug resistance to antimalarial treatments poses a serious public health risk globally. To combat this risk, molecular surveillance of drug resistance is imperative. We report the prevalence of mutations in the Plasmodium falciparum kelch 13 propeller domain associated with partial artemisinin resistance, which we determined by using Sanger sequencing samples from patients enrolled in therapeutic efficacy studies from 9 sub-Saharan countries during 2014-2018. Of the 2,865 samples successfully sequenced before treatment (day of enrollment) and on the day of treatment failure, 29 (1.0%) samples contained 11 unique nonsynonymous mutations and 83 (2.9%) samples contained 27 unique synonymous mutations. Two samples from Kenya contained the S522C mutation, which has been associated with delayed parasite clearance; however, no samples contained validated or candidate artemisinin-resistance mutations.

Artemisinin susceptibility in the malaria parasite Plasmodium falciparum: propellers, adaptor proteins and the need for cellular healing

Studies of the susceptibility of Plasmodium falciparum to the artemisinin family of antimalarial drugs provide a complex picture of partial resistance (tolerance) associated with increased parasite survival in vitro and in vivo. We present an overview of the genetic loci that, in mutant form, can independently elicit parasite tolerance. These encode Kelch propeller domain protein PfK13, ubiquitin hydrolase UBP-1, actin filament-organising protein Coronin, also carrying a propeller domain, and the trafficking adaptor subunit AP-2μ. Detailed studies of these proteins and the functional basis of artemisinin tolerance in blood-stage parasites are enabling a new synthesis of our understanding to date. To guide further experimental work, we present two major conclusions. First, we propose a dual-component model of artemisinin tolerance in P. falciparum comprising suppression of artemisinin activation in early ring stage by reducing endocytic haemoglobin capture from host cytosol, coupled with enhancement of cellular healing mechanisms in surviving cells. Second, these two independent requirements limit the likelihood of development of complete artemisinin resistance by P. falciparum, favouring deployment of existing drugs in new schedules designed to exploit these biological limits, thus extending the useful life of current combination therapies.

Prevalence of pfk13 and pfmdr1 polymorphisms in Bounkiling, Southern Senegal

Background

Delayed Plasmodium falciparum parasite clearance has been associated with Single Nucleotide Polymorphisms (SNPs) in the kelch protein propeller domain (coded by pfk13 gene). SNPs in the Plasmodium falciparum multidrug resistance gene 1 (pfmdr1) are associated with multi-drug resistance including the combination artemether-lumefantrine. To our knowledge, this is the first work providing information on the prevalence of k13-propeller and pfmdr1 mutations from Sédhiou, a region in the south of Senegal.

Methods

147 dried blood spots on filter papers were collected from symptomatic patients attending a hospital located in Bounkiling City, Sédhiou Region, Southern Senegal. All samples were collected between 2015–2017 during the malaria transmission season. Specific regions of the gene pfk13 and pfmdr1 were analyzed using PCR amplification and Sanger sequencing.

Results

The majority of parasites (92.9%) harboured the pfk13 wild type sequence and 6 samples harboured synonymous changes. Regarding pfmdr1, wild-type alleles represented the majority except at codon 184. Overall, prevalence of 86Y was 11.9%, 184F was 56.3% and 1246Y was 1.5%. The mutant allele 184F decreased from 73.7% in 2015 to 40.7% in 2017. The prevalence of haplotype NFD decreased from 71.4% in 2015 to 20.8% in 2017.

Conclusions

This study provides the first description of pfk13 and pfmdr1 genes variations in Bounkiling, a city in the Sédhiou Region of Senegal, contributing to closing the gap of information on anti-malaria drug resistance molecular markers in southern Senegal.

Multidrug-Resistant Plasmodium falciparum Parasites in the Central Highlands of Vietnam Jeopardize Malaria Control and Elimination Strategies

Plasmodium falciparum resistance to dihydroartemisinin-piperaquine has spread through the Greater Mekong Subregion to southwestern Vietnam. In 2018 to 2019, we collected 127 P. falciparum isolates from Dak Nong (36), Dak Lak (55), Gia Lai (13), and Kon Tum (23) provinces in Vietnam's Central Highlands and found parasites bearing the Pfkelch13 C580Y mutation and multiple plasmepsin 2/3 genes (mean prevalence, 17.9%; range, 4.3% to 27.8%), conferring resistance to dihydroartemisinin-piperaquine. This information is important for drug policy decisions in Vietnam.

Antimalarial Drug Resistance and Implications for the WHO Global Technical Strategy

PURPOSE OF REVIEW

Five years have passed since the World Health Organization released its Global Technical Strategy for Malaria (GTS). In that time, progress against malaria has plateaued. This review focuses on the implications of antimalarial drug resistance for the GTS and how interim progress in parasite genomics and antimalarial pharmacology offer a bulwark against it.

RECENT FINDINGS

For the first time, drug resistance-conferring genes have been identified and validated before their global expansion in malaria parasite populations. More efficient methods for their detection and elaboration have been developed, although low-density infections and polyclonality remain a nuisance to be solved. Clinical trials of alternative regimens for multidrug-resistant malaria have delivered promising results. New agents continue down the development pipeline, while a nascent infrastructure in sub-Saharan Africa for conducting phase I trials and trials of transmission-blocking agents has come to fruition after years of preparation.

SUMMARY

These and other developments can help inform the GTS as the world looks ahead to the next two decades of its implementation. To remain ahead of the threat that drug resistance poses, wider application of genomic-based surveillance and optimization of existing and forthcoming antimalarial drugs are essential.

Genomic and Transcriptomic Analysis for Identification of Genes and Interlinked Pathways Mediating Artemisinin Resistance in Leishmania donovani

Current therapy for visceral leishmaniasis (VL), compromised by drug resistance, toxicity, and high cost, demands for more effective, safer, and low-cost drugs. Artemisinin has been found to be an effectual drug alternative in experimental models of leishmaniasis. Comparative genome and transcriptome analysis of in vitro-adapted artesunate-resistant (K133AS-R) and -sensitive wild-type (K133WT) Leishmania donovani parasites was carried out using next-generation sequencing and single-color DNA microarray technology, respectively, to identify genes and interlinked pathways contributing to drug resistance. Whole-genome sequence analysis of K133WT vs. K133AS-R parasites revealed substantial variation among the two and identified 240 single nucleotide polymorphisms (SNPs), 237 insertion deletions (InDels), 616 copy number variations (CNVs) (377 deletions and 239 duplications), and trisomy of chromosome 12 in K133AS-R parasites. Transcriptome analysis revealed differential expression of 208 genes (fold change ≥ 2) in K133AS-R parasites. Functional categorization and analysis of modulated genes of interlinked pathways pointed out plausible adaptations in K133AS-R parasites, such as (i) a dependency on lipid and amino acid metabolism for generating energy, (ii) reduced DNA and protein synthesis leading to parasites in the quiescence state, and (iii) active drug efflux. The upregulated expression of cathepsin-L like protease, amastin-like surface protein, and amino acid transporter and downregulated expression of the gene encoding ABCG2, pteridine receptor, adenylatecyclase-type receptor, phosphoaceylglucosamine mutase, and certain hypothetical proteins are concordant with genomic alterations suggesting their potential role in drug resistance. The study provided an understanding of the molecular basis linked to artemisinin resistance in Leishmania parasites, which may be advantageous for safeguarding this drug for future use.

Parasite Viability as a Superior Measure of Antimalarial Drug Activity in Humans

BACKGROUND

Artemisinin derivatives are the leading class of antimalarial drugs due to their rapid onset of action and rapid clearance of circulating parasites. The parasite clearance half-life measures the rate of loss of parasites from blood after treatment, and this is currently used to assess antimalarial activity of novel agents and to monitor resistance. However, a number of recent studies have challenged the use of parasite clearance to measure drug activity, arguing that many circulating parasites may be nonviable.

METHODS

Plasmodium falciparum-infected subjects (n = 10) in a malaria volunteer infection study were administered a single dose of artesunate (2 mg/kg). Circulating parasite concentration was assessed by means of quantitative polymerase chain reaction (qPCR). Parasite viability after artesunate administration was estimated by mathematical modeling of the ex vivo growth of parasites collected from subjects.

RESULTS

We showed that in artemisinin-sensitive infection, viable parasites declined to <0.1% of baseline within 8 hours after artesunate administration, while the total number of circulating parasites measured with quantitative polymerase chain reaction remained unchanged. In artemisinin-resistant infections over the same interval, viable parasites declined to 51.4% (standard error of the mean, 4.6%) of baseline.

CONCLUSIONS

These results demonstrate that in vivo drug activity of artesunate is faster than is indicated by the parasite clearance half-life.

An insight into the recent development of the clinical candidates for the treatment of malaria and their target proteins

Malaria is an endemic disease, prevalent in tropical and subtropical regions which cost half of million deaths annually. The eradication of malaria is one of the global health priority nevertheless, current therapeutic efforts seem to be insufficient due to the emergence of drug resistance towards most of the available drugs, even first-line treatment ACT, unavailability of the vaccine, and lack of drugs with a new mechanism of action. Intensification of antimalarial research in recent years has resulted into the development of single dose multistage therapeutic agents which has advantage of overcoming the antimalarial drug resistance. The present review explored the current progress in the development of new promising antimalarials against prominent target proteins that have the potential to be a clinical candidate. Here, we also reviewed different aspects of drug resistance and highlighted new drug candidates that are currently in a clinical trial or clinical development, along with a few other molecules with excellent antimalarial activity overs ACTs. The summarized scientific value of previous approaches and structural features of antimalarials related to the activity are highlighted that will be helpful for the development of next-generation antimalarials.

Polymorphism analysis of propeller domain of k13 gene in Plasmodium ovale curtisi and Plasmodium ovale wallikeri isolates original infection from Myanmar and Africa in Yunnan Province, China

Background

Eighteen imported ovale malaria cases imported from Myanmar and various African countries have been reported in Yunnan Province, China from 2013 to 2018. All of them have been confirmed by morphological examination and 18S small subunit ribosomal RNA gene (18S rRNA) based PCR in YNRL. Nevertheless, the subtypes of Plasmodium ovale could not be identified based on 18S rRNA gene test, thus posing challenges on its accurate diagnosis. To help establish a more sensitive and specific method for the detection of P. ovale genes, this study performs sequence analysis on k13-propeller polymorphisms in P. ovale.

Methods

Dried blood spots (DBS) from ovale malaria cases were collected from January 2013 to December 2018, and the infection sources were confirmed according to epidemiological investigation. DNA was extracted, and the coding region (from 206th aa to 725th aa) in k13 gene propeller domain was amplified using nested PCR. Subsequently, the amplified products were sequenced and compared with reference sequence to obtain CDS. The haplotypes and mutation loci of the CDS were analysed, and the spatial structure of the amino acid peptide chain of k13 gene propeller domain was predicted by SWISS-MODEL.

Results

The coding region from 224th aa to 725th aa of k13 gene from P. ovale in 83.3% of collected samples (15/18) were amplified. Three haplotypes were observed in 15 samples, and the values of Ka/Ks, nucleic acid diversity index (π) and expected heterozygosity (He) were 3.784, 0.0095, and 0.4250. Curtisi haplotype, Wallikeri haplotype, and mutant type accounted for 73.3% (11/15), 20.0% (3/15), and 6.7% (1/15). The predominant haplotypes of P. ovale curtisi were determined in all five Myanmar isolates. Of the ten African isolates, six were identified as P. o. curtisi, three were P. o. wallikeri and one was mutant type. Base substitutions between the sequences of P. o. curtisi and P. o. wallikeri were determined at 38 loci, such as c.711. Moreover, the A > T base substitution at c.1428 was a nonsynonymous mutation, resulting in amino acid variation of T476S in the 476th position. Compared with sequence of P. o. wallikeri, the double nonsynonymous mutations of G > A and A > T at the sites of c.1186 and c.1428 leads to the variations of D396N and T476S for the 396th and 476th amino acids positions. For P. o. curtisi and P. o. wallikeri, the peptide chains in the coding region from 224th aa to 725th aa of k13 gene merely formed a monomeric spatial model, whereas the double-variant peptide chains of D396N and T476S formed homodimeric spatial model.

Conclusion

The propeller domain of k13 gene in the P. ovale isolates imported into Yunnan Province from Myanmar and Africa showed high differentiation. The sequences of Myanmar-imported isolates belong to P. o. curtisi, while the sequences of African isolates showed the sympatric distribution from P. o. curtisi, P. o. wallikeri and mutant isolates. The CDS with a double base substitution formed a dimeric spatial model to encode the peptide chain, which is completely different from the monomeric spatial structure to encode the peptide chain from P. o. curtisi and P. o. wallikeri.

Plasmodium falciparum Artemisinin Resistance: The Effect of Heme, Protein Damage, and Parasite Cell Stress Response

Despite a significant decline in morbidity and mortality over the last two decades, in 2018 there were 228 million reported cases of malaria and 405000 malaria-related deaths. Artemisinin, the cornerstone of artemisinin-based combination therapies, is the most potent drug in the antimalarial armamentarium against falciparum malaria. Heme-mediated activation of artemisinin and its derivatives results in widespread parasite protein alkylation, which is thought to lead to parasite death. Alarmingly, cases of decreased artemisinin efficacy have been widely detected across Cambodia and in neighboring countries, and a few cases have been reported in the Guiana Shield, India, and Africa. The grim prospect of widespread artemisinin resistance propelled a concerted effort to understand the mechanisms of artemisinin action and resistance. The identification of genetic markers and the knowledge of molecular mechanisms underpinning artemisinin resistance allow prospective surveillance and inform future drug development strategies, respectively. Here, we highlight recent advances in our understanding of how parasite vesicle trafficking, hemoglobin digestion, and cell stress responses contribute to artemisinin resistance.

Surveillance of genetic markers associated with Plasmodium falciparum resistance to artemisinin-based combination therapy in Pakistan, 2018-2019

BACKGROUND

The spread of artemisinin resistance in the Greater Mekong Subregion of Southeast Asia poses a significant threat for current anti-malarial treatment guidelines globally. The aim of this study was to assess the current prevalence of molecular markers of drug resistance in Plasmodium falciparum in the four provinces with the highest malaria burden in Pakistan, after introducing artemether-lumefantrine as first-line treatment in 2017.

METHODS

Samples were collected during routine malaria surveillance in Punjab, Sindh, Baluchistan, and Khyber Pakhtunkhwa provinces of Pakistan between January 2018 and February 2019. Plasmodium falciparum infections were confirmed by rapid diagnostic test or microscopy. Plasmodium falciparum positive isolates (n = 179) were screened by Sanger sequencing for single nucleotide polymorphisms (SNPs) in the P. falciparum kelch 13 (pfk13) propeller domain and in P. falciparum coronin (pfcoronin). SNPs in P. falciparum multidrug resistance 1 (pfmdr1) N86Y, Y184F, D1246Y and P. falciparum chloroquine resistance transporter (pfcrt) K76T were genotyped by PCR-restriction fragment length polymorphism.

RESULTS

No artemisinin resistance associated SNPs were identified in the pfk13 propeller domain or in pfcoronin. The pfmdr1 N86, 184F, D1246 and pfcrt K76 alleles associated with reduced lumefantrine sensitivity were present in 83.8% (150/179), 16.9% (29/172), 100.0% (173/173), and 8.4% (15/179) of all infections, respectively. The chloroquine resistance associated pfcrt 76T allele was present in 98.3% (176/179) of infections.

CONCLUSION

This study provides an update on the current prevalence of molecular markers associated with reduced P. falciparum sensitivity to artemether and/or lumefantrine in Pakistan, including a first baseline assessment of polymorphisms in pfcoronin. No mutations associated with artemisinin resistance were observed in pfk13 or pfcoronin. However, the prevalence of the pfmdr1 N86 and D1246 alleles, that have been associated with decreased susceptibility to lumefantrine, remain high. Although clinical and molecular data suggest that the current malaria treatment guidelines for P. falciparum are presently effective in Pakistan, close monitoring for artemisinin and lumefantrine resistance will be critical to ensure early detection and enhanced containment of emerging ACT resistance spreading across from Southeast Asia.

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.

Role of Plasmodium falciparum Kelch 13 Protein Mutations in P. falciparum Populations from Northeastern Myanmar in Mediating Artemisinin Resistance

Artemisinin resistance has emerged in Southeast Asia, endangering the substantial progress in malaria elimination worldwide. It is associated with mutations in the PfK13 protein, but how PfK13 mediates artemisinin resistance is not completely understood. Here we used a new antibody against PfK13 to show that the PfK13 protein is expressed in all stages of the asexual intraerythrocytic cycle as well as in gametocytes and is partially localized in the endoplasmic reticulum. By introducing four PfK13 mutations into the 3D7 strain and reverting these mutations in field parasite isolates, we determined the impacts of these mutations identified in the parasite populations from northern Myanmar on the ring stage using the in vitro ring survival assay. The introduction of the N458Y mutation into the 3D7 background significantly increased the survival rates of the ring-stage parasites but at the cost of the reduced fitness of the parasites. Introduction of the F446I mutation, the most prevalent PfK13 mutation in northern Myanmar, did not result in a significant increase in ring-stage survival after exposure to dihydroartemisinin (DHA), but these parasites showed extended ring-stage development. Further, parasites with the F446I mutation showed only a marginal loss of fitness, partially explaining its high frequency in northern Myanmar. Conversely, reverting all these mutations, except for the C469Y mutation, back to their respective wild types reduced the ring-stage survival of these isolates in response to in vitro DHA treatment.

Mutations in the Plasmodium falciparum Kelch 13 (PfK13) protein are associated with artemisinin resistance. PfK13 is essential for asexual erythrocytic development, but its function is not known. We tagged the PfK13 protein with green fluorescent protein in P. falciparum to study its expression and localization in asexual and sexual stages. We used a new antibody against PfK13 to show that the PfK13 protein is expressed ubiquitously in both asexual erythrocytic stages and gametocytes and is localized in punctate structures, partially overlapping an endoplasmic reticulum marker. We introduced into the 3D7 strain four PfK13 mutations (F446I, N458Y, C469Y, and F495L) identified in parasites from the China-Myanmar border area and characterized the in vitro artemisinin response phenotypes of the mutants. We found that all the parasites with the introduced PfK13 mutations showed higher survival rates in the ring-stage survival assay (RSA) than the wild-type (WT) control, but only parasites with N458Y displayed a significantly higher RSA value (26.3%) than the WT control. After these PfK13 mutations were reverted back to the WT in field parasite isolates, all revertant parasites except those with the C469Y mutation showed significantly lower RSA values than their respective parental isolates. Although the 3D7 parasites with introduced F446I, the predominant PfK13 mutation in northern Myanmar, did not show significantly higher RSA values than the WT, they had prolonged ring-stage development and showed very little fitness cost in in vitro culture competition assays. In comparison, parasites with the N458Y mutations also had a prolonged ring stage and showed upregulated resistance pathways in response to artemisinin, but this mutation produced a significant fitness cost, potentially leading to their lower prevalence in the Greater Mekong subregion.

Malaria Molecular Epidemiology: An Evolutionary Genetics Perspective

Malaria is a vector-borne disease that involves multiple parasite species in a variety of ecological settings. However, the parasite species causing the disease, the prevalence of subclinical infections, the emergence of drug resistance, the scale-up of interventions, and the ecological factors affecting malaria transmission, among others, are aspects that vary across areas where malaria is endemic. Such complexities have propelled the study of parasite genetic diversity patterns in the context of epidemiologic investigations. Importantly, molecular studies indicate that the time and spatial distribution of malaria cases reflect epidemiologic processes that cannot be fully understood without characterizing the evolutionary forces shaping parasite population genetic patterns. Although broad in scope, this review in the Microbiology Spectrum Curated Collection: Advances in Molecular Epidemiology highlights the need for understanding population genetic concepts when interpreting parasite molecular data. First, we discuss malaria complexity in terms of the parasite species involved. Second, we describe how molecular data are changing our understanding of malaria incidence and infectiousness. Third, we compare different approaches to generate parasite genetic information in the context of epidemiologically relevant questions related to malaria control. Finally, we describe a few Plasmodium genomic studies as evidence of how these approaches will provide new insights into the malaria disease dynamics. *This article is part of a curated collection.

Overexpression of the HECT ubiquitin ligase PfUT prolongs the intraerythrocytic cycle and reduces invasion efficiency of Plasmodium falciparum

The glms ribozyme system has been used as an amenable tool to conditionally control expression of genes of interest. It is generally assumed that insertion of the ribozyme sequence does not affect expression of the targeted gene in the absence of the inducer glucosamine-6-phosphate, although experimental support for this assumption is scarce. Here, we report the unexpected finding that integration of the glms ribozyme sequence in the 3′ untranslated region of a gene encoding a HECT E3 ubiquitin ligase, termed Plasmodium falciparum ubiquitin transferase (PfUT), increased steady state RNA and protein levels 2.5-fold in the human malaria parasite P. falciparum. Overexpression of pfut resulted in an S/M phase-associated lengthening of the parasite’s intraerythrocytic developmental cycle and a reduced merozoite invasion efficiency. The addition of glucosamine partially restored the wild type phenotype. Our study suggests a role of PfUT in controlling cell cycle progression and merozoite invasion. Our study further raises awareness regarding unexpected effects on gene expression when inserting the glms ribozyme sequence into a gene locus.

Profiles of Kelch mutations in Plasmodium falciparum across South Asia and their implications for tracking drug resistance

Artemisinin-based combination therapy (ACT) offers highly successful treatment of malaria. Emergence and spread of Plasmodium falciparum (Pf) parasites with decreased susceptibility to ACT in South-East Asia has caused concern worldwide. The current accepted criteria to assess artemisinin (ART) resistance relies upon data on treatment failure, delayed parasite clearance at day 3 (DPC3), parasite clearance half-life (PCHL) and in-vitro/ex-vivo ring stage survival assays (RSAs). Interestingly, some studies suggest that DPC3 does not provide a distinct separation between ART sensitive/resistant strains, and RSA differences may also be inconclusive. More recently, recrudescence of ART treated Pf, independent of the presence of Kelch 13 (K13) mutation (C580Y), has been reported in the monkey malaria model suggesting that genes other than K13 like coronin, dhps, dhfr, crt, mdr1 and plasmepsin1 may contribute towards ACT failure. Here we have collated the distribution of K13 mutants from Pf strains in South Asia. A total of fifty Pf-K13 mutations have been studied for ART resistance in South Asia of which nine have been validated while eleven are potentials for ART resistance. The remaining thirty K13 mutations have been reported from various locations in South Asia but lack corroborative clinical data on ART resistance/ACT failure. Of the fifty, fourteen K13 mutations have been identified in India including four novel mutations (S549Y, G625R, N657H, D702N). Structural mapping of these K13 mutations does not offer any coherent explanation for their contribution towards ART resistance as they are scattered in the K13 structure. Thus, K13 mutations likely provide only a partial synopsis, and we propose that all suspect cases of ACT failure be assessed by: 1) DPC3, 2) PCHL, 3) in-vitro/ex-vivo RSAs and 4) GWAS data in an effort to annotate the resistance status of the parasites. These efforts may help in surveillance and containment of ART resistance/ACT failure in South Asia.

K13-propeller gene polymorphisms of Plasmodium falciparum and the therapeutic effect of artesunate among migrant workers returning to Guangxi, China (2014-2017)

BACKGROUND

The resistance of Plasmodium falciparum to artemisinin has been identified in Asia and some parts of Africa. The drug resistance of P. falciparum will be an obstacle to the successful elimination of malaria by 2025. Whole-genome sequencing of the artemisinin-resistant parasite line revealed mutations on the k13 gene associated with drug resistance in P. falciparum. To understand the artemisinin resistance of the imported P. falciparum cases from Africa, the mutations in the k13 gene in parasites from imported malaria cases in Guangxi Province were detected and the treatment efficiency of artesunate monotherapy was observed.

METHODS

DNA was extracted from 319 blood samples from migrant workers with P. falciparum infection who returned to their hometown in Guangxi Province from Africa between 2014 and 2017. The k13-propeller gene was amplified by nested PCR, and sequencing, gene mutation frequency and geographic difference of imported P. falciparum cases were analysed by comparison with the wild-type strain. Of 319 patients, 158 were P. falciparum-infected and were treated with intravenous injection of artesunate and were observed, including the time of asexual stage clearance and the dose of artesunate used.

RESULTS

Of the 319 P. falciparum samples, 12 samples had the k13-propeller mutation, and 11 point mutations were detected; 5 were non-synonymous mutations (T474I, A481T, A578S, V603E, G665S) and were not associated with artemisinin resistance. The clinical treatment observation showed that the median (IQR) dose of artesunate for peripheral blood parasite asexual stage clearance was 407.55 (360-510) mg, and the D3 parasite clearance rate was 70.25%, including the five k13-propeller mutations of P. falciparum. After 7 days of treatment, 98.73% of cases were cleared. Two cases were treated with artemisinin for 8 days with a 960-mg dose to completely clear the asexual parasite, but they did not have a mutation in the k13 gene.

CONCLUSIONS

Five mutations of the k13-propeller gene in 319 P. falciparum samples from patients returning from Africa were identified. The frequency of the k13-propeller mutants was low, and the mutations were not strongly associated with artemisinin resistance. The median (IQR) dose of artesunate monotherapy in actual clinical treatment to remove asexual parasite stages was 407.55 (360-510) mg, equivalent to D3-D4. Some P. falciparum cases without a k13-propeller mutation showed obvious delayed clearance of the parasite from peripheral blood. Trial registration The diagnosis of malaria and the treatment of malaria-infected patients are the routine work of Centres for Disease Control and Prevention. Information on the patients was conveyed with the patient's approval, and the research aim, methods, risks and benefits of the study were explained in detail to the patients.

A systematic review on anti-malarial drug discovery and antiplasmodial potential of green synthesis mediated metal nanoparticles: overview, challenges and future perspectives

BACKGROUND

The recent emergence in Southeast Asia of artemisinin resistance poses major threats to malaria control and elimination globally. Green nanotechnologies can constitute interesting tools for discovering anti-malarial medicines. This systematic review focused on the green synthesis of metal nanoparticles as potential source of new antiplasmodial drugs.

METHODS

Seven electronic database were used following the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) guidelines.

RESULTS

A total of 17 papers were included in the systematic review. 82.4% of the studies used plant leaves to produce nanoparticles (NPs) while three studies used microorganisms, including bacteria and fungi. Silver was the main metal precursor for the synthesis of NPs. The majority of studies obtained nanoparticles spherical in shape, with sizes ranging between 4 and 65 nm, and reported no or little cytotoxic effect of the NPs. Results based on 50% inhibitory concentration (IC50) varied between studies but, in general, could be divided into three NP categories; (i) those more effective than positive controls, (ii) those more effective than corresponding plant extracts and, (iii) those less effective than the positive controls or plant extracts.

CONCLUSIONS

This study highlights the high antiplasmodial potential of green-synthesized metal nanoparticles thereby underscoring the possibility to find and develop new anti-malarial drugs based on green synthesis approaches. However, the review also highlights the need for extensive in vitro and in vivo studies to confirm their safety in humans and the elucidation of the mechanism of action.

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.

The Diversity of the Plasmodium falciparum K13 Propeller Domain Did Not Increase after Implementation of Artemisinin-Based Combination Therapy in Uganda

Artemisinin-based combination therapies (ACTs) are the standard of care to treat uncomplicated falciparum malaria. However, resistance to artemisinins, defined as delayed parasite clearance after therapy, has emerged in Southeast Asia, and the spread of resistance to sub-Saharan Africa could have devastating consequences. Artemisinin resistance has been associated in Southeast Asia with multiple nonsynonymous single nucleotide polymorphisms (NS-SNPs) in the propeller domain of the gene encoding the Plasmodium falciparum K13 protein (K13PD). Some K13PD NS-SNPs have been seen in Africa, but the relevance of these mutations is unclear. To assess whether ACT use has selected for specific K13PD mutations, we compared the K13PD genetic diversity in clinical isolates collected before and after the implementation of ACT use from seven sites across Uganda. We detected K13PD NS-SNPs in 16 of 683 (2.3%) clinical isolates collected between 1999 and 2004 and in 26 of 716 (3.6%) isolates collected between 2012 and 2016 (P = 0.16), representing a total of 29 different polymorphisms at 27 codons. Individual NS-SNPs were usually detected only once, and none were found in more than 0.7% of the isolates. Three SNPs (C469F, P574L, and A675V) associated with delayed clearance in Southeast Asia were seen in samples collected between 2012 and 2016, each in a single isolate. No differences in diversity following implementation of ACT use were found at any of the seven sites, nor was there evidence of selective pressures acting on the locus. Our results suggest that selection by ACTs is not impacting on K13PD diversity in Uganda.

K-13 propeller gene polymorphisms isolated between 2014 and 2017 from Cameroonian Plasmodium falciparum malaria patients

The emergence of artemisinin-resistant parasites since the late 2000s at the border of Cambodia and Thailand poses serious threats to malaria control globally, particularly in Africa which bears the highest malaria transmission burden. This study aimed to obtain reliable data on the current state of the kelch13 molecular marker for artemisinin resistance in Plasmodium falciparum in Cameroon. DNA was extracted from the dried blood spots collected from epidemiologically distinct endemic areas in the Center, Littoral and North regions of Cameroon. Nested PCR products from the Kelch13-propeller gene were sequenced and analyzed on an ABI 3730XL automatic sequencer. Of 219 dried blood spots, 175 were sequenced successfully. We identified six K13 mutations in 2.9% (5/175) of samples, including 2 non-synonymous, the V589I allele had been reported in Africa already and one new allele E612K had not been reported yet. These two non-synonymous mutations were uniquely found in parasites from the Littoral region. One sample showed two synonymous mutations within the kelch13 gene. We also observed two infected samples with mixed K13 mutant and K13 wild-type infection. Taken together, our data suggested the circulation of the non-synonymous K13 mutations in Cameroon. Albeit no mutations known to be associated with parasite clearance delays in the study population, there is need for continuous surveillance for earlier detection of resistance as long as ACTs are used and scaled up in the community.

Genetic affinities of an eradicated European Plasmodium falciparum strain

Malaria was present in most of Europe until the second half of the 20th century, when it was eradicated through a combination of increased surveillance and mosquito control strategies, together with cross-border and political collaboration. Despite the severe burden of malaria on human populations, it remains contentious how the disease arrived and spread in Europe. Here, we report a partial Plasmodium falciparum nuclear genome derived from a set of antique medical slides stained with the blood of malaria-infected patients from Spain’s Ebro Delta, dating to the 1940s. Our analyses of the genome of this now eradicated European P. falciparum strain confirms stronger phylogeographical affinity to present-day strains in circulation in central south Asia, rather than to those in Africa. This points to a longitudinal, rather than a latitudinal, spread of malaria into Europe. In addition, this genome displays two derived alleles in the pfmrp1 gene that have been associated with drug resistance. Whilst this could represent standing variation in the ancestral P. falciparum population, these mutations may also have arisen due to the selective pressure of quinine treatment, which was an anti-malarial drug already in use by the time the sample we sequenced was mounted on a slide.

In vitro susceptibility of Plasmodium falciparum isolates from the China-Myanmar border area to artemisinins and correlation with K13 mutations

Mutations in the Kelch domain of the K13 gene (PF3D7_1343700) were previously associated with artemisinin resistance in Plasmodium falciparum. This study followed the dynamics of the K13 polymorphisms in P. falciparum parasites from the China-Myanmar border area obtained in 2007-2016, and their in vitro sensitivities to artesunate (AS) and dihydroartemisinin (DHA). The 50% effective concentration (EC5072h) values of 133 culture-adapted field isolates to AS and DHA, measured by the conventional 72 h SYBR Green I-based assay, varied significantly among the parasites from different years; all were significantly higher than that of the reference strain 3D7. Compared with parasites from 2007 to 2008, ring survival rates almost doubled in parasites obtained in later years. Sequencing the full-length K13 genes identified 11 point mutations present in 85 (63.9%) parasite isolates. F446I was the predominant (55/133) variant, and its frequency was increased from 17.6% (3/17) in 2007 to 55.9% (19/34) in 2014-2016. No wild-type (WT) Kelch domain sequences were found in the 34 samples obtained from 2014 to 2016. In the 2014-2016 samples, a new mutation (G533S) appeared and reached 44.1% (15/34). Collectively, parasites with the Kelch domain mutations (after amino acid 440) had significantly higher ring survival rates than the WT parasites. Individually, F446I, G533S and A676D showed significantly higher ring survival rates than the WT. Although the drug sensitivity phenotypes measured by the RSA6h and EC5072h assays may be intrinsically linked to the in vivo clinical efficacy data, the values determined by these two assays were not significantly correlated. This study identified the trend of K13 mutations in parasite populations from the China-Myanmar border area, confirmed an overall correlation of Kelch domain mutations with elevated ring-stage survival rates, and emphasized the importance of monitoring the evolution and spread of parasites with reduced artemisinin sensitivity along the malaria elimination course.

Plasmodium Genomics and Genetics: New Insights into Malaria Pathogenesis, Drug Resistance, Epidemiology, and Evolution

Protozoan Plasmodium parasites are the causative agents of malaria, a deadly disease that continues to afflict hundreds of millions of people every year. Infections with malaria parasites can be asymptomatic, with mild or severe symptoms, or fatal, depending on many factors such as parasite virulence and host immune status. Malaria can be treated with various drugs, with artemisinin-based combination therapies (ACTs) being the first-line choice. Recent advances in genetics and genomics of malaria parasites have contributed greatly to our understanding of parasite population dynamics, transmission, drug responses, and pathogenesis. However, knowledge gaps in parasite biology and host-parasite interactions still remain. Parasites resistant to multiple antimalarial drugs have emerged, while advanced clinical trials have shown partial efficacy for one available vaccine. Here we discuss genetic and genomic studies of Plasmodium biology, host-parasite interactions, population structures, mosquito infectivity, antigenic variation, and targets for treatment and immunization. Knowledge from these studies will advance our understanding of malaria pathogenesis, epidemiology, and evolution and will support work to discover and develop new medicines and vaccines.

Kelch 13 propeller gene polymorphism among Plasmodium falciparum isolates in Lagos, Nigeria: Molecular Epidemiologic Study

OBJECTIVE

To assess polymorphism in Kelch 13 gene of Plasmodium falciparum isolates in Lagos, Nigeria.

METHODS

195 Plasmodium falciparum-positive dried blood spots collected from individuals that accessed diagnostic care at some health facilities and during community surveys across several Local Government Areas of Lagos State, Nigeria, were investigated for the presence of mutations in the K13 gene by nested polymerase chain reaction (PCR) using haplotype-specific probes and sequencing.

RESULTS

Three mutant genotypes of K13 gene were observed: A578S in 0.5%, D464N in 0.5% and Q613H in 1.5%. The frequency of K13 polymorphism was 3.1%, while the remaining parasite population had the wild K13 propeller genes.

CONCLUSION

No validated Kelch 13 polymorphism associated with artemisinin resistance was seen among P. falciparum isolates from Lagos, Nigeria. As no clinical study was done, this could not be correlated with artemisinin sensitivity.

Ozonide Antimalarial Activity in the Context of Artemisinin-Resistant Malaria

The ozonides are one of the most advanced drug classes in the antimalarial development pipeline and were designed to improve on limitations associated with current front-line artemisinin-based therapies. Like the artemisinins, the pharmacophoric peroxide bond of ozonides is essential for activity, and it appears that these antimalarials share a similar mode of action, raising the possibility of cross-resistance. Resistance to artemisinins is associated with Plasmodium falciparum mutations that allow resistant parasites to escape short-term artemisinin-mediated damage (elimination half-life ~1 h). Importantly, some ozonides (e.g., OZ439) have a sustained in vivo drug exposure profile, providing a major pharmacokinetic advantage over the artemisinin derivatives. Here, we describe recent progress made towards understanding ozonide antimalarial activity and discuss ozonide utility within the context of artemisinin resistance.

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.