Fitness Costs and the Rapid Spread of kelch13-C580Y Substitutions Conferring Artemisinin Resistance

A Plasmodium falciparum genetic cross reveals the contributions of pfcrt and plasmepsin II/III to piperaquine drug resistance

Piperaquine (PPQ) is widely used in combination with dihydroartemisinin as a first-line treatment against malaria. Multiple genetic drivers of PPQ resistance have been reported, including mutations in the Plasmodium falciparum chloroquine resistance transporter (pfcrt) and increased copies of plasmepsin II/III (pm2/3). We generated a cross between a Cambodia-derived multidrug-resistant KEL1/PLA1 lineage isolate (KH004) and a drug-susceptible Malawian parasite (Mal31). Mal31 harbors a wild-type (3D7-like) pfcrt allele and a single copy of pm2/3, while KH004 has a chloroquine-resistant (Dd2-like) pfcrt allele with an additional G367C substitution and multiple copies of pm2/3. We recovered 104 unique recombinant parasites and examined a targeted set of progeny representing all possible combinations of variants at pfcrt and pm2/3. We performed a detailed analysis of competitive fitness and a range of PPQ susceptibility phenotypes with these progenies, including PPQ survival assay, area under the dose response curve, and a limited point IC50. We find that inheritance of the KH004 pfcrt allele is required for reduced PPQ sensitivity, whereas copy number variation in pm2/3 further decreases susceptibility but does not confer resistance in the absence of additional mutations in pfcrt. A deep investigation of genotype-phenotype relationships demonstrates that progeny clones from experimental crosses can be used to understand the relative contributions of pfcrt, pm2/3, and parasite genetic background to a range of PPQ-related traits. Additionally, we find that the resistance phenotype associated with parasites inheriting the G367C substitution in pfcrt is consistent with previously validated PPQ resistance mutations in this transporter.IMPORTANCEResistance to piperaquine, used in combination with dihydroartemisinin, has emerged in Cambodia and threatens to spread to other malaria-endemic regions. Understanding the causal mutations of drug resistance and their impact on parasite fitness is critical for surveillance and intervention and can also reveal new avenues to limiting the evolution and spread of drug resistance. An experimental genetic cross is a powerful tool for pinpointing the genetic determinants of key drug resistance and fitness phenotypes and has the distinct advantage of quantifying the effects of naturally evolved genetic variation. Our study was strengthened since the full range of copies of KH004 pm2/3 was inherited among the progeny clones, allowing us to directly test the role of the pm2/3 copy number on resistance-related phenotypes in the context of a unique pfcrt allele. Our multigene model suggests an important role for both loci in the evolution of this multidrug-resistant parasite lineage.

Evaluation of Dihydroartemisinin-Piperaquine Efficacy and Molecular Markers in Uncomplicated Falciparum Patients: A Study across Binh Phuoc and Dak Nong, Vietnam

Background and Objectives: Malaria continues to be a significant global health challenge. The efficacy of artemisinin-based combination therapies (ACTs) has declined in many parts of the Greater Mekong Subregion, including Vietnam, due to the spread of resistant malaria strains. This study was conducted to assess the efficacy of the Dihydroartemisinin (DHA)-Piperaquine (PPQ) regimen in treating uncomplicated falciparum malaria and to conduct molecular surveillance of antimalarial drug resistance in Binh Phuoc and Dak Nong provinces. Materials and Methods: The study included 63 uncomplicated malaria falciparum patients from therapeutic efficacy studies (TES) treated following the WHO treatment guidelines (2009). Molecular marker analysis was performed on all 63 patients. Methods encompassed Sanger sequencing for pfK13 mutations and quantitative real-time PCR for the pfpm2 gene. Results: This study found a marked decrease in the efficacy of the DHA-PPQ regimen, with an increased rate of treatment failures at two study sites. Genetic analysis revealed a significant presence of pfK13 mutations and pfpm2 amplifications, indicating emerging resistance to artemisinin and its partner drug. Conclusions: The effectiveness of the standard DHA-PPQ regimen has sharply declined, with rising treatment failure rates. This decline necessitates a review and possible revision of national malaria treatment guidelines. Importantly, molecular monitoring and clinical efficacy assessments together provide a robust framework for understanding and addressing detection drug resistance in malaria.

tRNA modification reprogramming contributes to artemisinin resistance in Plasmodium falciparum

Plasmodium falciparum artemisinin (ART) resistance is driven by mutations in kelch-like protein 13 (PfK13). Quiescence, a key aspect of resistance, may also be regulated by a yet unidentified epigenetic pathway. Transfer RNA modification reprogramming and codon bias translation is a conserved epitranscriptomic translational control mechanism that allows cells to rapidly respond to stress. We report a role for this mechanism in ART-resistant parasites by combining tRNA modification, proteomic and codon usage analyses in ring-stage ART-sensitive and ART-resistant parasites in response to drug. Post-drug, ART-resistant parasites differentially hypomodify mcm5s2U on tRNA and possess a subset of proteins, including PfK13, that are regulated by Lys codon-biased translation. Conditional knockdown of the terminal s2U thiouridylase, PfMnmA, in an ART-sensitive parasite background led to increased ART survival, suggesting that hypomodification can alter the parasite ART response. This study describes an epitranscriptomic pathway via tRNA s2U reprogramming that ART-resistant parasites may employ to survive ART-induced stress.

Emergence, transmission dynamics and mechanisms of artemisinin partial resistance in malaria parasites in Africa

Malaria, mostly due to Plasmodium falciparum infection in Africa, remains one of the most important infectious diseases in the world. Standard treatment for uncomplicated P. falciparum malaria is artemisinin-based combination therapy (ACT), which includes a rapid-acting artemisinin derivative plus a longer-acting partner drug, and standard therapy for severe P. falciparum malaria is intravenous artesunate. The efficacy of artemisinins and ACT has been threatened by the emergence of artemisinin partial resistance in Southeast Asia, mediated principally by mutations in the P. falciparum Kelch 13 (K13) protein. High ACT treatment failure rates have occurred when resistance to partner drugs is also seen. Recently, artemisinin partial resistance has emerged in Rwanda, Uganda and the Horn of Africa, with independent emergences of different K13 mutants in each region. In this Review, we summarize our current knowledge of artemisinin partial resistance and focus on the emergence of resistance in Africa, including its epidemiology, transmission dynamics and mechanisms. At present, the clinical impact of emerging resistance in Africa is unclear and most available evidence suggests that the efficacies of leading ACTs remain excellent, but there is an urgent need to better appreciate the extent of the problem and its consequences for the treatment and control of malaria.

Critical interdependencies between Plasmodium nutrient flux and drugs

Nutrient import and waste efflux are critical dependencies for intracellular Plasmodium falciparum parasites. Nutrient transport proteins are often lineage specific and can provide unique targets for antimalarial drug development. P. falciparum nutrient transport pathways can be a double-edged sword for the parasite, not only mediating the import of nutrients and excretion of waste products but also providing an access route for drugs. Here we briefly summarise the nutrient acquisition pathways of intracellular P. falciparum blood-stage parasites and then highlight how these pathways influence many aspects relevant to antimalarial drugs, resulting in complex and often underappreciated interdependencies.

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.

Fitness Costs of pfhrp2 and pfhrp3 Deletions Underlying Diagnostic Evasion in Malaria Parasites

BACKGROUND

Rapid diagnostic tests based on detection of histidine-rich proteins (HRPs) are widely used for malaria diagnosis, but parasites carrying pfhrp deletions can evade detection and are increasing in frequency in some countries. Models aim to predict conditions under which pfhrp2 and/or pfhrp3 deletions will increase, but a key parameter-the fitness cost of deletions-is unknown.

METHODS

We removed pfhrp2 and/or pfhrp3 from a Malawian parasite clone using gene editing approaches) and measured fitness costs by conducting pairwise competition experiments.

RESULTS

We observed significant fitness costs of 0.087 ± 0.008 (1 standard error) per asexual cycle for pfhrp2 deletion and 0.113 ± 0.008 for the pfhrp2/3 double deletion, relative to the unedited progenitor parasite. Selection against deletions is strong and comparable to that resulting from drug resistance mutations.

CONCLUSIONS

Prior modeling suggested that diagnostic selection may drive increased frequency of pfhrp deletions only when fitness costs are mild. Our experiments show that costs of pfhrp deletions are higher than these thresholds, but modeling and empirical results can be reconciled if the duration of infection is short. These results may inform future modeling to understand why pfhrp2/3 deletions are increasing in some locations (Ethiopia and Eritrea) but not in others (Mekong region).

Barcoding Genetically Distinct Plasmodium falciparum Strains for Comparative Assessment of Fitness and Antimalarial Drug Resistance

The repeated emergence of antimalarial drug resistance in Plasmodium falciparum, including to the current frontline antimalarial artemisinin, is a perennial problem for malaria control. Next-generation sequencing has greatly accelerated the identification of polymorphisms in resistance-associated genes but has also highlighted the need for more sensitive and accurate laboratory tools to profile current and future antimalarials and to quantify the impact of drug resistance acquisition on parasite fitness. The interplay of fitness and drug response is of fundamental importance in understanding why particular genetic backgrounds are better at driving the evolution of drug resistance in natural populations, but the impact of parasite fitness landscapes on the epidemiology of drug resistance has typically been laborious to accurately quantify in the lab, with assays being limited in accuracy and throughput. Here we present a scalable method to profile fitness and drug response of genetically distinct P. falciparum strains with well-described sensitivities to several antimalarials. We leverage CRISPR/Cas9 genome-editing and barcode sequencing to track unique barcodes integrated into a nonessential gene (pfrh3). We validate this approach in multiplex competitive growth assays of three strains with distinct geographical origins. Furthermore, we demonstrate that this method can be a powerful approach for tracking artemisinin response as it can identify an artemisinin resistant strain within a mix of multiple parasite lines, suggesting an approach for scaling the laborious ring-stage survival assay across libraries of barcoded parasite lines. Overall, we present a novel high-throughput method for multiplexed competitive growth assays to evaluate parasite fitness and drug response. IMPORTANCE The complex interplay between antimalarial resistance and parasite fitness has important implications for understanding the development and spread of drug resistance alleles and the impact of genetic background on transmission. One limitation with current methodologies to measure parasite fitness is the ability to scale this beyond simple head-to-head competition experiments between a wildtype control line and test line, with a need for a scalable approach that allows tracking of parasite growth in complex mixtures. In our study, we have used CRISPR editing to insert unique DNA barcodes into a safe-harbor genomic locus to tag multiple parasite strains and use next-generation sequencing to read out strain dynamics. We observe inherent fitness differences between the strains, as well as sensitive modulation of responses to challenge with clinically relevant antimalarials, including artemisinin.

Leveraging genome editing to functionally evaluate Plasmodium diversity

The ambitious goal of malaria elimination requires an in-depth understanding of the parasite's biology to counter the growing threat of antimalarial resistance and immune evasion. Timely assessment of the functional impact of antigenic diversity in the early stages of vaccine development will be critical for achieving the goal of malaria control, elimination, and ultimately eradication. Recent advances in targeted genome editing enabled the functional validation of resistance-associated markers in Plasmodium falciparum, the deadliest malaria-causing pathogen and strain-specific immune neutralization. This review explores recent advances made in leveraging genome editing to aid the functional evaluation of Plasmodium diversity and highlights how these techniques can assist in prioritizing both therapeutic and vaccine candidates.

A Malaria Parasite Cross Reveals Genetic Determinants of Plasmodium falciparum Growth in Different Culture Media

What genes determine in vitro growth and nutrient utilization in asexual blood-stage malaria parasites? Competition experiments between NF54, clone 3D7, a lab-adapted African parasite, and a recently isolated Asian parasite (NHP4026) reveal contrasting outcomes in different media: 3D7 outcompetes NHP4026 in media containing human serum, while NHP4026 outcompetes 3D7 in media containing AlbuMAX, a commercial lipid-rich bovine serum formulation. To determine the basis for this polymorphism, we conducted parasite genetic crosses using humanized mice and compared genome-wide allele frequency changes in three independent progeny populations cultured in media containing human serum or AlbuMAX. This bulk segregant analysis detected three quantitative trait loci (QTL) regions [on chromosome (chr) 2 containing aspartate transaminase AST; chr 13 containing EBA-140; and chr 14 containing cysteine protease ATG4] linked with differential growth in serum or AlbuMAX in each of the three independent progeny pools. Selection driving differential growth was strong (s = 0.10 – 0.23 per 48-hour lifecycle). We conducted validation experiments for the strongest QTL on chr 13: competition experiments between ΔEBA-140 and 3D7 wildtype parasites showed fitness reversals in the two medium types as seen in the parental parasites, validating this locus as the causative gene. These results (i) demonstrate the effectiveness of bulk segregant analysis for dissecting fitness traits in P. falciparum genetic crosses, and (ii) reveal intimate links between red blood cell invasion and nutrient composition of growth media. Use of parasite crosses combined with bulk segregant analysis will allow systematic dissection of key nutrient acquisition/metabolism and red blood cell invasion pathways in P. falciparum.

Nutrient Limitation Magnifies Fitness Costs of Antimalarial Drug Resistance Mutations

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

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.

Comparative analysis of targeted next-generation sequencing for Plasmodium falciparum drug resistance markers

Well-defined molecular resistance markers are available for a range of antimalarial drugs, and molecular surveillance is increasingly important for monitoring antimalarial drug resistance. Different genotyping platforms are available, but these have not been compared in detail. We compared Targeted Amplicon Deep sequencing (TADs) using Ion Torrent PGM with Illumina MiSeq for the typing of antimalarial drug resistance genes. We developed and validated protocols to type the molecular resistance markers pfcrt, pfdhfr, pfdhps, pfmdr1, pfkelch, and pfcytochrome b, in Plasmodium falciparum for the Ion Torrent PGM and Illumina MiSeq sequencing platforms. With P. falciparum 3D7 and K1 as reference strains, whole blood samples (N = 20) and blood spots from Rapid Diagnostic Test (RDT) samples (N = 5) from patients with uncomplicated falciparum malaria from Ubon Ratchathani were assessed on both platforms and compared for coverage (average reads per amplicon), sequencing accuracy, variant accuracy, false positive rate, false negative rate, and alternative allele detection, with conventional Sanger sequencing as the reference method for SNP calling. Both whole blood and RDT samples could be successfully sequenced using the Ion Torrent PGM and Illumina MiSeq platforms. Coverage of reads per amplicon was higher with Illumina MiSeq (28,886 reads) than with Ion Torrent PGM (1754 reads). In laboratory generated artificial mixed infections, the two platforms could detect the minor allele down to 1% density at 500X coverage. SNPs calls from both platforms were in complete agreement with conventional Sanger sequencing. The methods can be multiplexed with up to 96 samples per run, which reduces cost by 86% compared to conventional Sanger sequencing. Both platforms, using the developed TAD protocols, provide an accurate method for molecular surveillance of drug resistance markers in P. falciparum, but Illumina MiSeq provides higher coverage than Ion Torrent PGM.

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.

Genetic analysis of praziquantel response in schistosome parasites implicates a transient receptor potential channel

[Figure: see text].

The parasitophorous vacuole nutrient channel is critical for drug access in malaria parasites and modulates the artemisinin resistance fitness cost

Intraerythrocytic malaria parasites proliferate bounded by a parasitophorous vacuolar membrane (PVM). The PVM contains nutrient permeable channels (NPCs) conductive to small molecules, but their relevance for parasite growth for individual metabolites is largely untested. Here we show that growth-relevant levels of major carbon and energy sources pass through the NPCs. Moreover, we find that NPCs are a gate for several antimalarial drugs, highlighting their permeability properties as a critical factor for drug design. Looking into NPC-dependent amino acid transport, we find that amino acid shortage is a reason for the fitness cost in artemisinin-resistant (ART^R) parasites and provide evidence that NPC upregulation to increase amino acids acquisition is a mechanism of ART^R parasites in vitro and in human infections to compensate this fitness cost. Hence, the NPCs are important for nutrient and drug access and reveal amino acid deprivation as a critical constraint in ART^R parasites.

Increase in the proportion of Plasmodium falciparum with kelch13 C580Y mutation and decline in pfcrt and pfmdr1 mutant alleles in Papua New Guinea

Background

The C580Y mutation in the Plasmodium falciparum kelch13 gene is the most commonly observed variant in artemisinin-resistant isolates in the Greater Mekong Subregion (GMS). Until 2017, it had not been identified outside the GMS, except for Guyana/Amazonia. In 2017, three parasites carrying the C580Y mutation were identified in Papua New Guinea (PNG). As the C580Y allele rapidly spread in the GMS, there is concern that this mutant is now spreading in PNG.

Methods

In 2020, a cross-sectional survey was conducted at two clinics in Wewak, PNG. Symptomatic patients infected with P. falciparum were treated with artemether plus lumefantrine following a national treatment policy. Blood samples were obtained before treatment, and polymorphisms in kelch13, pfcrt, and pfmdr1 were determined. Parasite positivity was examined on day 3. The results were compared with those of previous studies conducted in 2002, 2003, and 2016–2018.

Results

A total of 94 patients were included in this analysis. The proportion of C580Y was significantly increased (2.2% in 2017, 5.7% in 2018, and 6.4% in 2020; p = 4.2 × 10^–3). A significant upward trend was observed in the wild-type proportion for pfcrt (1.9% in 2016 to 46.7% in 2020; p = 8.9 × 10^–16) and pfmdr1 (59.5% in 2016 to 91.4% in 2020; p = 2.3 × 10^–6). Among 27 patients successfully followed on day 3, including three with C580Y infections, none showed positive parasitaemia.

Conclusions

Under the conditions of significant increases in pfcrt K76 and pfmdr1 N86 alleles in PNG, the increase in kelch13 C580Y mutants may be a warning indicator of the emergence of parasites resistant to the currently used first-line treatment regimen of artemether plus lumefantrine. Therefore, nationwide surveillance of molecular markers for drug resistance and assessment of its therapeutic effects are important.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12936-021-03933-6.

Fitness of sulfadoxine-resistant Plasmodium berghei harboring a single mutation in dihydropteroate synthase (DHPS)

Genetic changes conferring drug resistance are generally believed to impose fitness costs to pathogens in the absence of the drug. However, the fitness of resistant parasites against sulfadoxine/pyrimethamine has been inconclusive in Plasmodium falciparum. This is because resistance is conferred by the complex combination of mutations in dihydropteroate synthase (dhps) and dihydrofolate reductase (dhfr), which makes it difficult to separately assess the extent and magnitude of the costs imposed by mutations in dhps and dhfr. To assess the fitness costs imposed by sulfadoxine resistance alone, we generated a transgenic rodent malaria parasite, P. berghei clone harboring an A394G mutation in dhps (PbDHPS-A394G), corresponding to the causative mutation for sulfadoxine resistance in P. falciparum (PfDHPS-A437G). A four-day suppressive test confirmed that the PbDHPS-A394G clone was resistant to sulfadoxine. PbDHPS-A394G and wild-type clones showed similar growth rates and gametocyte production. This observation was confirmed in competitive experiments in which PbDHPS-A394G and wild-type clones were co-infected into mice to directly assess the survival competition between them. In the mosquitoes, there were no significant differences in oocyst production between PbDHPS-A394G and wild-type. These results indicate that the PbDHPS-A394G mutation alters the parasites to sulfadoxine resistance but may not impose fitness disadvantages during the blood stages in mice and oocyst formation in mosquitoes. These results partly explain the persistence of the PfDHPS-A437G mutant in the natural parasite populations.

Artemisinin and multidrug-resistant Plasmodium falciparum - a threat for malaria control and elimination

PURPOSE OF REVIEW

Artemisinin-based combination therapies (ACTs) are globally the first-line treatment for uncomplicated falciparum malaria and new compounds will not be available within the next few years. Artemisinin-resistant Plasmodium falciparum emerged over a decade ago in the Greater Mekong Subregion (GMS) and, compounded by ACT partner drug resistance, has caused significant ACT treatment failure. This review provides an update on the epidemiology, and mechanisms of artemisinin resistance and approaches to counter multidrug-resistant falciparum malaria.

RECENT FINDINGS

An aggressive malaria elimination programme in the GMS has helped prevent the spread of drug resistance to neighbouring countries. However, parasites carrying artemisinin resistance-associated mutations in the P. falciparum Kelch13 gene (pfk13) have now emerged independently in multiple locations elsewhere in Asia, Africa and South America. Notably, artemisinin-resistant infections with parasites carrying the pfk13 R561H mutation have emerged and spread in Rwanda.

SUMMARY

Enhancing the geographic coverage of surveillance for resistance will be key to ensure prompt detection of emerging resistance in order to implement effective countermeasures without delay. Treatment strategies designed to prevent the emergence and spread of multidrug resistance must be considered, including deployment of triple drug combination therapies and multiple first-line therapies.

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.

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.

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 K13 mutations in Africa and Asia impact artemisinin resistance and parasite fitness

The emergence of mutant K13-mediated artemisinin (ART) resistance in Plasmodium falciparum malaria parasites has led to widespread treatment failures across Southeast Asia. In Africa, K13-propeller genotyping confirms the emergence of the R561H mutation in Rwanda and highlights the continuing dominance of wild-type K13 elsewhere. Using gene editing, we show that R561H, along with C580Y and M579I, confer elevated in vitro ART resistance in some African strains, contrasting with minimal changes in ART susceptibility in others. C580Y and M579I cause substantial fitness costs, which may slow their dissemination in high-transmission settings, in contrast with R561H that in African 3D7 parasites is fitness neutral. In Cambodia, K13 genotyping highlights the increasing spatio-temporal dominance of C580Y. Editing multiple K13 mutations into a panel of Southeast Asian strains reveals that only the R561H variant yields ART resistance comparable to C580Y. In Asian Dd2 parasites C580Y shows no fitness cost, in contrast with most other K13 mutations tested, including R561H. Editing of point mutations in ferredoxin or mdr2, earlier associated with resistance, has no impact on ART susceptibility or parasite fitness. These data underline the complex interplay between K13 mutations, parasite survival, growth and genetic background in contributing to the spread of ART resistance.

Immune selection suppresses the emergence of drug resistance in malaria parasites but facilitates its spread

Although drug resistance in Plasmodium falciparum typically evolves in regions of low transmission, resistance spreads readily following introduction to regions with a heavier disease burden. This suggests that the origin and the spread of resistance are governed by different processes, and that high transmission intensity specifically impedes the origin. Factors associated with high transmission, such as highly immune hosts and competition within genetically diverse infections, are associated with suppression of resistant lineages within hosts. However, interactions between these factors have rarely been investigated and the specific relationship between adaptive immunity and selection for resistance has not been explored. Here, we developed a multiscale, agent-based model of Plasmodium parasites, hosts, and vectors to examine how host and parasite dynamics shape the evolution of resistance in populations with different transmission intensities. We found that selection for antigenic novelty (“immune selection”) suppressed the evolution of resistance in high transmission settings. We show that high levels of population immunity increased the strength of immune selection relative to selection for resistance. As a result, immune selection delayed the evolution of resistance in high transmission populations by allowing novel, sensitive lineages to remain in circulation at the expense of the spread of a resistant lineage.

In contrast, in low transmission settings, we observed that resistant strains were able to sweep to high population prevalence without interference. Additionally, we found that the relationship between immune selection and resistance changed when resistance was widespread. Once resistance was common enough to be found on many antigenic backgrounds, immune selection stably maintained resistant parasites in the population by allowing them to proliferate, even in untreated hosts, when resistance was linked to a novel epitope. Our results suggest that immune selection plays a role in the global pattern of resistance evolution.

Drug resistance in the malaria parasite, Plasmodium falciparum, presents an ongoing public health challenge, but aspects of its evolution are poorly understood. Although antimalarial resistance is common worldwide, it can typically be traced to just a handful of evolutionary origins. Counterintuitively, although Sub Saharan Africa bears 90% of the global malaria burden, resistance typically originates in regions where transmission intensity is low. In high transmission regions, infections are genetically diverse, and hosts have significant standing adaptive immunity, both of which are known to suppress the frequency of resistance within infections. However, interactions between immune-driven selection, transmission intensity, and resistance have not been investigated. Using a multiscale, agent-based model, we found that high transmission intensity slowed the evolution of resistance via its effect on host population immunity. High host immunity strengthened selection for antigenic novelty, interfering with selection for resistance and allowing sensitive lineages to suppress resistant lineages in untreated hosts. However, once resistance was common in the circulating parasite population, immune selection maintained it in the population at a high prevalence. Our findings provide a novel explanation for observations about the origin of resistance and suggest that adaptive immunity is a critical component of selection.

Balanced impacts of fitness and drug pressure on the evolution of PfMDR1 polymorphisms in Plasmodium falciparum

Background

Anti-malarial drug resistance may be limited by decreased fitness in resistant parasites. Important contributors to resistance are mutations in the Plasmodium falciparum putative drug transporter PfMDR1.

Methods

Impacts on in vitro fitness of two common PfMDR1 polymorphisms, N86Y, which is associated with sensitivity to multiple drugs, and Y184F, which has no clear impact on drug sensitivity, were evaluated to study associations between resistance mediators and parasite fitness, measured as relative growth in competitive culture experiments. NF10 P. falciparum lines engineered to represent all PfMDR1 N86Y and Y184F haplotypes were co-cultured for 40 days, and the genetic make-up of the cultures was characterized every 4 days by pyrosequencing. The impacts of culture with anti-malarials on the growth of different haplotypes were also assessed. Lastly, the engineering of P. falciparum containing another common polymorphism, PfMDR1 D1246Y, was attempted.

Results

Co-culture results were as follows. With wild type (WT) Y184 fixed (N86/Y184 vs. 86Y/Y184), parasites WT and mutant at 86 were at equilibrium. With mutant 184 F fixed (N86/184F vs. 86Y/184F), mutants at 86 overgrew WT. With WT N86 fixed (N86/Y184 vs. N86/184F), WT at 184 overgrew mutants. With mutant 86Y fixed (86Y/Y184 vs. 86Y/184F), WT and mutant at 86 were at equilibrium. Parasites with the double WT were in equilibrium with the double mutant, but 86Y/Y184 overgrew N86/184F. Overall, WT N86/mutant 184F parasites were less fit than parasites with all other haplotypes. Parasites engineered for another mutation, PfMDR1 1246Y, were unstable in culture, with reversion to WT over time. Thus, the N86 WT is stable when accompanied by the Y184 WT, but incurs a fitness cost when accompanied by mutant 184F. Culturing in the presence of chloroquine favored 86Y mutant parasites and in the presence of lumefantrine favored N86 WT parasites; piperaquine had minimal impact.

Conclusions

These results are consistent with those for Ugandan field isolates, suggest reasons for varied haplotypes, and highlight the interplay between drug pressure and fitness that is guiding the evolution of resistance-mediating haplotypes in P. falciparum.

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.

Transmission of Artemisinin-Resistant Malaria Parasites to Mosquitoes under Antimalarial Drug Pressure

Resistance to artemisinin-based combination therapy (ACT) in the Plasmodium falciparum parasite is threatening to reverse recent gains in reducing global deaths from malaria. While resistance manifests as delayed parasite clearance in patients, the phenotype can only spread geographically via the sexual stages and mosquito transmission. In addition to their asexual killing properties, artemisinin and its derivatives sterilize sexual male gametocytes. Whether resistant parasites overcome this sterilizing effect has not, however, been fully tested.

Resistance to artemisinin-based combination therapy (ACT) in the Plasmodium falciparum parasite is threatening to reverse recent gains in reducing global deaths from malaria. While resistance manifests as delayed parasite clearance in patients, the phenotype can only spread geographically via the sexual stages and mosquito transmission. In addition to their asexual killing properties, artemisinin and its derivatives sterilize sexual male gametocytes. Whether resistant parasites overcome this sterilizing effect has not, however, been fully tested. Here, we analyzed P. falciparum clinical isolates from the Greater Mekong Subregion, each demonstrating delayed clinical clearance and known resistance-associated polymorphisms in the Kelch13 (PfK13^var) gene. As well as demonstrating reduced asexual sensitivity to drug, certain PfK13^var isolates demonstrated a marked reduction in sensitivity to artemisinin in an in vitro male gamete formation assay. Importantly, this same reduction in sensitivity was observed when the most resistant isolate was tested directly in mosquito feeds. These results indicate that, under artemisinin drug pressure, while sensitive parasites are blocked, resistant parasites continue transmission. This selective advantage for resistance transmission could favor acquisition of additional host-specificity or polymorphisms affecting partner drug sensitivity in mixed infections. Favored resistance transmission under ACT coverage could have profound implications for the spread of multidrug-resistant malaria beyond Southeast Asia.

Plasmodium berghei K13 Mutations Mediate In Vivo Artemisinin Resistance That Is Reversed by Proteasome Inhibition

Recent successes in malaria control have been seriously threatened by the emergence of Plasmodium falciparum parasite resistance to the frontline artemisinin drugs in Southeast Asia. P. falciparum artemisinin resistance is associated with mutations in the parasite K13 protein, which associates with a delay in the time required to clear the parasites upon drug treatment. Gene editing technologies have been used to validate the role of several candidate K13 mutations in mediating P. falciparum artemisinin resistance in vitro under laboratory conditions. Nonetheless, the causal role of these mutations under in vivo conditions has been a matter of debate. Here, we have used CRISPR/Cas9 gene editing to introduce K13 mutations associated with artemisinin resistance into the related rodent-infecting parasite, Plasmodium berghei. Phenotyping of these P. berghei K13 mutant parasites provides evidence of their role in mediating artemisinin resistance in vivo, which supports in vitro artemisinin resistance observations. However, we were unable to introduce some of the P. falciparum K13 mutations (C580Y and I543T) into the corresponding amino acid residues, while other introduced mutations (M476I and R539T equivalents) carried pronounced fitness costs. Our study provides evidence of a clear causal role of K13 mutations in modulating susceptibility to artemisinins in vitro and in vivo using the well-characterized P. berghei model. We also show that inhibition of the P. berghei proteasome offsets parasite resistance to artemisinins in these mutant lines.

The recent emergence of Plasmodium falciparum parasite resistance to the first line antimalarial drug artemisinin is of particular concern. Artemisinin resistance is primarily driven by mutations in the P. falciparum K13 protein, which enhance survival of early ring-stage parasites treated with the artemisinin active metabolite dihydroartemisinin in vitro and associate with delayed parasite clearance in vivo. However, association of K13 mutations with in vivo artemisinin resistance has been problematic due to the absence of a tractable model. Herein, we have employed CRISPR/Cas9 genome editing to engineer selected orthologous P. falciparum K13 mutations into the K13 gene of an artemisinin-sensitive Plasmodium berghei rodent model of malaria. Introduction of the orthologous P. falciparum K13 F446I, M476I, Y493H, and R539T mutations into P. berghei K13 yielded gene-edited parasites with reduced susceptibility to dihydroartemisinin in the standard 24-h in vitro assay and increased survival in an adapted in vitro ring-stage survival assay. Mutant P. berghei K13 parasites also displayed delayed clearance in vivo upon treatment with artesunate and achieved faster recrudescence upon treatment with artemisinin. Orthologous C580Y and I543T mutations could not be introduced into P. berghei, while the equivalents of the M476I and R539T mutations resulted in significant growth defects. Furthermore, a Plasmodium-selective proteasome inhibitor strongly synergized dihydroartemisinin action in these P. berghei K13 mutant lines, providing further evidence that the proteasome can be targeted to overcome artemisinin resistance. Taken together, our findings provide clear experimental evidence for the involvement of K13 polymorphisms in mediating susceptibility to artemisinins in vitro and, most importantly, under in vivo conditions.

Characterization and mitigation of fragmentation enzyme-induced dual stranded artifacts

High-throughput short-read sequencing relies on fragmented DNA for optimal sampling of input nucleic acid. Several vendors now offer proprietary enzyme cocktails as a cheaper and more streamlined method of fragmentation when compared to acoustic shearing. We have discovered that these enzymes induce the formation of library molecules containing regions of nearby DNA from opposite strands. Sequencing reads derived from these molecules can lead to artifact-derived variant calls appearing at variant allele frequencies <5%. We present Fragmentation Artifact Detection and Elimination (FADE), software to remove these artifacts from mapped reads and mitigate artifact-related effects on downstream analysis. We find that the artifacts principally affect downstream analyses that are sensitive to a 1-3% artifact bias in the sequencing reads, such as targeted resequencing and rare variant discovery.

Humanized Mice and the Rebirth of Malaria Genetic Crosses

The first experimental crosses carried out with the human malaria parasite Plasmodium falciparum played a key role in determining the genetic loci responsible for drug resistance, virulence, invasion, growth rate, and transmission. These crosses relied on splenectomized chimpanzees to complete the liver stage of the parasite's life cycle and the subsequent transition to asexual blood stage culture followed by cloning of recombinant progeny in vitro. Crosses can now be routinely carried out using human-liver-chimeric mice infused with human erythrocytes to generate hundreds of unique recombinant progeny for genetic linkage mapping, bulk segregant analysis, and high-throughput 'omics readouts. The high number of recombinant progeny should allow for unprecedented power and efficiency in the execution of a systems genetics approach to study P. falciparum biology.

Sand Fly Studies Predict Transmission Potential of Drug-resistant Leishmania

Leishmania parasites have the capacity to rapidly adapt to changing environments in their digenetic life cycle which alternates between a vertebrate and an invertebrate host. Emergence of resistance following drug exposure can evoke phenotypic alterations that affect several aspects of parasite fitness in both hosts. Current studies of the impact of resistance are mostly limited to interactions with the mammalian host and characterization of in vitro parasite growth and differentiation. Development in the vector and transmission capacity have been largely ignored. This review reflects on the impact of drug resistance on its spreading potential with specific focus on the use of the sand fly infection model to evaluate parasite development in the vector and the ensuing transmission potential of drug-resistant phenotypes.

Local emergence in Amazonia of Plasmodium falciparum k13 C580Y mutants associated with in vitro artemisinin resistance

Antimalarial drug resistance has historically arisen through convergent de novo mutations in Plasmodium falciparum parasite populations in Southeast Asia and South America. For the past decade in Southeast Asia, artemisinins, the core component of first-line antimalarial therapies, have experienced delayed parasite clearance associated with several pfk13 mutations, primarily C580Y. We report that mutant pfk13 has emerged independently in Guyana, with genome analysis indicating an evolutionary origin distinct from Southeast Asia. Pfk13 C580Y parasites were observed in 1.6% (14/854) of samples collected in Guyana in 2016-2017. Introducing pfk13 C580Y or R539T mutations by gene editing into local parasites conferred high levels of in vitro artemisinin resistance. In vitro growth competition assays revealed a fitness cost associated with these pfk13 variants, potentially explaining why these resistance alleles have not increased in frequency more quickly in South America. These data place local malaria control efforts at risk in the Guiana Shield.

A Kelch13-defined endocytosis pathway mediates artemisinin resistance in malaria parasites

Artemisinin and its derivatives (ARTs) are the frontline drugs against malaria, but resistance is jeopardizing their effectiveness. ART resistance is mediated by mutations in the parasite's Kelch13 protein, but Kelch13 function and its role in resistance remain unclear. In this study, we identified proteins located at a Kelch13-defined compartment. Inactivation of eight of these proteins, including Kelch13, rendered parasites resistant to ART, revealing a pathway critical for resistance. Functional analysis showed that these proteins are required for endocytosis of hemoglobin from the host cell. Parasites with inactivated Kelch13 or a resistance-conferring Kelch13 mutation displayed reduced hemoglobin endocytosis. ARTs are activated by degradation products of hemoglobin. Hence, reduced activity of Kelch13 and its interactors diminishes hemoglobin endocytosis and thereby ART activation, resulting in parasite resistance.

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.

Genetic mapping of fitness determinants across the malaria parasite Plasmodium falciparum life cycle

Determining the genetic basis of fitness is central to understanding evolution and transmission of microbial pathogens. In human malaria parasites (Plasmodium falciparum), most experimental work on fitness has focused on asexual blood stage parasites, because this stage can be easily cultured, although the transmission of malaria requires both female Anopheles mosquitoes and vertebrate hosts. We explore a powerful approach to identify the genetic determinants of parasite fitness across both invertebrate and vertebrate life-cycle stages of P. falciparum. This combines experimental genetic crosses using humanized mice, with selective whole genome amplification and pooled sequencing to determine genome-wide allele frequencies and identify genomic regions under selection across multiple lifecycle stages. We applied this approach to genetic crosses between artemisinin resistant (ART-R, kelch13-C580Y) and ART-sensitive (ART-S, kelch13-WT) parasites, recently isolated from Southeast Asian patients. Two striking results emerge: we observed (i) a strong genome-wide skew (>80%) towards alleles from the ART-R parent in the mosquito stage, that dropped to ~50% in the blood stage as selfed ART-R parasites were selected against; and (ii) repeatable allele specific skews in blood stage parasites with particularly strong selection (selection coefficient (s) ≤ 0.18/asexual cycle) against alleles from the ART-R parent at loci on chromosome 12 containing MRP2 and chromosome 14 containing ARPS10. This approach robustly identifies selected loci and has strong potential for identifying parasite genes that interact with the mosquito vector or compensatory loci involved in drug resistance.

Cutting back malaria: CRISPR/Cas9 genome editing of Plasmodium

CRISPR/Cas9 approaches are revolutionizing our ability to perform functional genomics across a wide range of organisms, including the Plasmodium parasites that cause malaria. The ability to deliver single point mutations, epitope tags and gene deletions at increased speed and scale is enabling our understanding of the biology of these complex parasites, and pointing to potential new therapeutic targets. In this review, we describe some of the biological and technical considerations for designing CRISPR-based experiments, and discuss potential future developments that broaden the applications for CRISPR/Cas9 interrogation of the malaria parasite genome.

Gene silencing through RNAi and antisense Vivo-Morpholino increases the efficacy of pyrethroids on larvae of Anopheles stephensi

Background

Insecticides are still at the core of insect pest and vector control programmes. Several lines of evidence indicate that ABC transporters are involved in detoxification processes against insecticides, including permethrin and other pyrethroids. In particular, the ABCG4 gene, a member of the G subfamily, has consistently been shown to be up-regulated in response to insecticide treatments in the mosquito malaria vector Anopheles stephensi (both adults and larvae).

Methods

To verify the actual involvement of this transmembrane protein in the detoxification process of permethrin, bioassays on larvae of An. stephensi, combining the insecticide with a siRNA, specifically designed for the inhibition of ABCG4 gene expression were performed. Administration to larvae of the same siRNA, labeled with a fluorescent molecule, was effected to investigate the systemic distribution of the inhibitory RNA into the larval bodies. Based on siRNA results, similar experiments using antisense Vivo-Morpholinos (Vivo-MOs) were effected. These molecules, compared to siRNA, are expected to guarantee a higher stability in environmental conditions and in the insect gut, and present thus a higher potential for future in-field applications.

Results

Bioassays using two different concentrations of siRNA, associated with permethrin, led to an increase of larval mortality, compared with results with permethrin alone. These outcomes confirm that ABCG4 transporter plays a role in the detoxification process against the selected insecticide. Moreover, after fluorescent labelling, it was shown the systemic dissemination of siRNA in different body districts of An. stephensi larvae, which suggest a potential systemic effect of the molecule. At the same time, results of Vivo-MO experiments were congruent with those obtained using siRNA, thus confirming the potential of ABCG4 inhibition as a strategy to increase permethrin susceptibility in mosquitoes. For the first time, Vivo-MOs were administered in water to larvae, with evidence for a biological effect.

Conclusions

Targeting ABCG4 gene for silencing through both techniques resulted in an increased pyrethroid efficacy. These results open the way toward the possibility to exploit ABCG4 inhibition in the context of integrated programmes for the control An. stephensi mosquitoes and malaria transmission.

Pairwise growth competitions identify relative fitness relationships among artemisinin resistant Plasmodium falciparum field isolates

Background

Competitive outcomes between co-infecting malaria parasite lines can reveal fitness disparities in blood stage growth. Blood stage fitness costs often accompany the evolution of drug resistance, with the expectation that relatively fitter parasites will be more likely to spread in populations. With the recent emergence of artemisinin resistance, it is important to understand the relative competitive fitness of the metabolically active asexual blood stage parasites. Genetically distinct drug resistant parasite clones with independently evolved sets of mutations are likely to vary in asexual proliferation rate, contributing to their chance of transmission to the mosquito vector.

Methods

An optimized in vitro 96-well plate-based protocol was used to quantitatively measure-head-to-head competitive fitness during blood stage development between seven genetically distinct field isolates from a hotspot of emerging artemisinin resistance and the laboratory strain, NF54. These field isolates were isolated from patients in Southeast Asia carrying different alleles of kelch13 and included both artemisinin-sensitive and artemisinin-resistant isolates. Fluorescent labeled microsatellite markers were used to track the relative densities of each parasite throughout the co-growth period of 14–60 days. All-on-all competitions were conducted for the panel of eight parasite lines (28 pairwise competitions) to determine their quantitative competitive fitness relationships.

Results

Twenty-eight pairwise competitive growth outcomes allowed for an unambiguous ranking among a set of seven genetically distinct parasite lines isolated from patients in Southeast Asia displaying a range of both kelch13 alleles and clinical clearance times and a laboratory strain, NF54. This comprehensive series of assays established the growth relationships among the eight parasite lines. Interestingly, a clinically artemisinin resistant parasite line that carries the wild-type form of kelch13 outcompeted all other parasites in this study. Furthermore, a kelch13 mutant line (E252Q) was competitively more fit without drug than lines with other resistance-associated kelch13 alleles, including the C580Y allele that has expanded to high frequencies under drug pressure in Southeast Asian resistant populations.

Conclusions

This optimized competitive growth assay can be employed for assessment of relative growth as an index of fitness during the asexual blood stage growth between natural lines carrying different genetic variants associated with artemisinin resistance. Improved understanding of the fitness costs of different parasites proliferating in human blood and the role different resistance mutations play in the context of specific genetic backgrounds will contribute to an understanding of the potential for specific mutations to spread in populations, with the potential to inform targeted strategies for malaria therapy.

Evolution and Genetic Diversity of the k13 Gene Associated with Artemisinin Delayed Parasite Clearance in Plasmodium falciparum

Mutations in the Plasmodium falciparum k13 (Pfk13) gene are linked to delayed parasite clearance in response to artemisinin-based combination therapies (ACTs) in Southeast Asia. To explore the evolutionary rate and constraints acting on this gene, k13 orthologs from species sharing a recent common ancestor with P. falciparum and Plasmodium vivax were analyzed. These comparative studies were followed by genetic polymorphism analyses within P. falciparum using 982 complete Pfk13 sequences from public databases and new data obtained by next-generation sequencing from African and Haitian isolates. Although k13 orthologs evolve at heterogeneous rates, the gene was conserved across the genus, with only synonymous substitutions being found at residues where mutations linked to the delayed parasite clearance phenotype have been reported. This suggests that those residues were under constraint from undergoing nonsynonymous changes during evolution of the genus. No fixed nonsynonymous differences were found between Pfk13 and its orthologs in closely related species found in African apes. This indicates that all nonsynonymous substitutions currently found in Pfk13 are younger than the time of divergence between P. falciparum and its closely related species. At the population level, no mutations linked to delayed parasite clearance were found in our samples from Africa and Haiti. However, there is a high number of single Pfk13 mutations segregating in P. falciparum populations, and two predominant alleles are distributed worldwide. This pattern is discussed in terms of how changes in the efficacy of natural selection, affected by population expansion, may have allowed for the emergence of mutations tolerant to ACTs.

Detection of low-density Plasmodium falciparum infections using amplicon deep sequencing

BACKGROUND

Deep sequencing of targeted genomic regions is becoming a common tool for understanding the dynamics and complexity of Plasmodium infections, but its lower limit of detection is currently unknown. Here, a new amplicon analysis tool, the Parallel Amplicon Sequencing Error Correction (PASEC) pipeline, is used to evaluate the performance of amplicon sequencing on low-density Plasmodium DNA samples. Illumina-based sequencing of two Plasmodium falciparum genomic regions (CSP and SERA2) was performed on two types of samples: in vitro DNA mixtures mimicking low-density infections (1-200 genomes/μl) and extracted blood spots from a combination of symptomatic and asymptomatic individuals (44-653,080 parasites/μl). Three additional analysis tools-DADA2, HaplotypR, and SeekDeep-were applied to both datasets and the precision and sensitivity of each tool were evaluated.

RESULTS

Amplicon sequencing can contend with low-density samples, showing reasonable detection accuracy down to a concentration of 5 Plasmodium genomes/μl. Due to increased stochasticity and background noise, however, all four tools showed reduced sensitivity and precision on samples with very low parasitaemia (< 5 copies/μl) or low read count (< 100 reads per amplicon). PASEC could distinguish major from minor haplotypes with an accuracy of 90% in samples with at least 30 Plasmodium genomes/μl, but only 61% at low Plasmodium concentrations (< 5 genomes/μl) and 46% at very low read counts (< 25 reads per amplicon). The four tools were additionally used on a panel of extracted parasite-positive blood spots from natural malaria infections. While all four identified concordant patterns of complexity of infection (COI) across four sub-Saharan African countries, COI values obtained for individual samples differed in some cases.

CONCLUSIONS

Amplicon deep sequencing can be used to determine the complexity and diversity of low-density Plasmodium infections. Despite differences in their approach, four state-of-the-art tools resolved known haplotype mixtures with similar sensitivity and precision. Researchers can therefore choose from multiple robust approaches for analysing amplicon data, however, error filtration approaches should not be uniformly applied across samples of varying parasitaemia. Samples with very low parasitaemia and very low read count have higher false positive rates and call for read count thresholds that are higher than current default recommendations.

Global Spread of Mutant PfCRT and Its Pleiotropic Impact on Plasmodium falciparum Multidrug Resistance and Fitness

The global spread of Plasmodium falciparum chloroquine resistance transporter (PfCRT) variant haplotypes earlier caused the widespread loss of chloroquine (CQ) efficacy. In Asia, novel PfCRT mutations that emerged on the Dd2 allelic background have recently been implicated in high-level resistance to piperaquine, and N326S and I356T have been associated with genetic backgrounds in which resistance emerged to artemisinin derivatives. By analyzing large-scale genome sequencing data, we report that the predominant Asian CQ-resistant Dd2 haplotype is undetectable in Africa. Instead, the GB4 and previously unexplored Cam783 haplotypes predominate, along with wild-type, drug-sensitive PfCRT that has reemerged as the major haplotype. To interrogate how these alleles impact drug susceptibility, we generated pfcrt-modified isogenic parasite lines spanning the mutational interval between GB4 and Dd2, which includes Cam783 and involves amino acid substitutions at residues 326 and 356. Relative to Dd2, the GB4 and Cam783 alleles were observed to mediate lower degrees of resistance to CQ and the first-line drug amodiaquine, while resulting in higher growth rates. These findings suggest that differences in growth rates, a surrogate of parasite fitness, influence selection in the context of African infections that are frequently characterized by high transmission rates, mixed infections, increased immunity, and less recourse to treatment. We also observe that the Asian Dd2 allele affords partial protection against piperaquine yet does not directly impact artemisinin efficacy. Our results can help inform the regional recommendations of antimalarials, whose activity is influenced by and, in certain cases, enhanced against select PfCRT variant haplotypes.IMPORTANCE Our study defines the allelic distribution of pfcrt, an important mediator of multidrug resistance in Plasmodium falciparum, in Africa and Asia. We leveraged whole-genome sequence analysis and gene editing to demonstrate how current drug combinations can select different allelic variants of this gene and shape region-specific parasite population structures. We document the ability of PfCRT mutations to modulate parasite susceptibility to current antimalarials in dissimilar, pfcrt allele-specific ways. This study underscores the importance of actively monitoring pfcrt genotypes to identify emerging patterns of multidrug resistance and help guide region-specific treatment options.

Artemisinin-Resistant Malaria as a Global Catastrophic Biological Threat

The global spread of artemisinin resistance brings with it the threat of incurable malaria. Already, this disease threatens over 219 million lives per year and causes 5-6% losses in GDP in endemic areas, even with current advances in prevention and treatment. This chapter discusses the currently tenuous position we are in globally, and the impact that could be seen if artemisinin treatment is lost, whether due to the unchecked spread of K13 mutations or poor global investment in treatment and prevention advances. Artemisinin is the backbone of current ACT treatment programs and severe malarial treatment; without it, the success of future malaria eradication programs will be in jeopardy.