Genetic diversity and chloroquine selective sweeps in Plasmodium falciparum

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.

Two decades of molecular surveillance in Senegal reveal rapid changes in known drug resistance mutations over time

BACKGROUND

Drug resistance in Plasmodium falciparum is a major threat to malaria control efforts. Pathogen genomic surveillance could be invaluable for monitoring current and emerging parasite drug resistance.

METHODS

Data from two decades (2000-2020) of continuous molecular surveillance of P. falciparum parasites from Senegal were retrospectively examined to assess historical changes in malaria drug resistance mutations. Several known drug resistance markers and their surrounding haplotypes were profiled using a combination of single nucleotide polymorphism (SNP) molecular surveillance and whole genome sequence based population genomics.

RESULTS

This dataset was used to track temporal changes in drug resistance markers whose timing correspond to historically significant events such as the withdrawal of chloroquine (CQ) and the introduction of sulfadoxine-pyrimethamine (SP) in 2003. Changes in the mutation frequency at Pfcrt K76T and Pfdhps A437G coinciding with the 2014 introduction of seasonal malaria chemoprevention (SMC) in Senegal were observed. In 2014, the frequency of Pfcrt K76T increased while the frequency of Pfdhps A437G declined. Haplotype-based analyses of Pfcrt K76T showed that this rapid increase was due to a recent selective sweep that started after 2014.

DISCUSSION (CONCLUSION)

The rapid increase in Pfcrt K76T is troubling and could be a sign of emerging amodiaquine (AQ) resistance in Senegal. Emerging AQ resistance may threaten the future clinical efficacy of artesunate-amodiaquine (ASAQ) and AQ-dependent SMC chemoprevention. These results highlight the potential of molecular surveillance for detecting rapid changes in parasite populations and stress the need to monitor the effectiveness of AQ as a partner drug for artemisinin-based combination therapy (ACT) and for chemoprevention.

Temporal and spatial dynamics of Plasmodium falciparum clonal lineages in Guyana

Plasmodium parasites, the causal agents of malaria, are eukaryotic organisms that obligately undergo sexual recombination within mosquitoes. In low transmission settings, parasites recombine with themselves, and the clonal lineage is propagated rather than broken up by outcrossing. We investigated whether stochastic/neutral factors drive the persistence and abundance of Plasmodium falciparum clonal lineages in Guyana, a country with relatively low malaria transmission, but the only setting in the Americas in which an important artemisinin resistance mutation (pfk13 C580Y) has been observed. We performed whole genome sequencing on 1,727 Plasmodium falciparum samples collected from infected patients across a five-year period (2016-2021). We characterized the relatedness between each pair of monoclonal infections (n = 1,409) through estimation of identity-by-descent (IBD) and also typed each sample for known or candidate drug resistance mutations. A total of 160 multi-isolate clones (mean IBD ≥ 0.90) were circulating in Guyana during the study period, comprising 13 highly related clusters (mean IBD ≥ 0.40). In the five-year study period, we observed a decrease in frequency of a mutation associated with artemisinin partner drug (piperaquine) resistance (pfcrt C350R) and limited co-occurence of pfcrt C350R with duplications of plasmepsin 2/3, an epistatic interaction associated with piperaquine resistance. We additionally observed 61 nonsynonymous substitutions that increased markedly in frequency over the study period as well as a novel pfk13 mutation (G718S). However, P. falciparum clonal dynamics in Guyana appear to be largely driven by stochastic factors, in contrast to other geographic regions, given that clones carrying drug resistance polymorphisms do not demonstrate enhanced persistence or higher abundance than clones carrying polymorphisms of comparable frequency that are unrelated to resistance. The use of multiple artemisinin combination therapies in Guyana may have contributed to the disappearance of the pfk13 C580Y mutation.

Sanger sequencing and deconvolution of polyclonal infections: a quantitative approach to monitor drug-resistant Plasmodium falciparum

BACKGROUND

Anti-malarial drug resistance in Plasmodium falciparum is a major public health problem in malaria-endemic regions. Although various technical improvements in sequencing methods have been introduced to identify SNPs, the conventional approach with current tools does not discriminate mixed infections, and thus can be improved for more sensitive surveillance of anti-malarial resistance to better inform control strategies.

METHODS

We developed a computational approach for deconvolution of chromatograms generated by standard Sanger sequencing of PCR amplicons in order to quantify molecular marker variants of anti-malarial drug resistance genes [Plasmodium falciparum dihydropteorate synthase (Pfdhps) and P. falciparum dihydrofolate reductase (Pfdhfr)]. We validated this computational approach using mixtures of V1/S and FCR3 at varying proportions between 0 and 100%, then applied it to field samples collected in Doneguebougou, Mali in 2018. We determined the mean fraction of resistance alleles in individual samples, as well as the prevalence of infections carrying resistant parasites.

FINDINGS

We observed a highly significant correlation between the predicted and measured proportions of V1/S and FCR3 alleles in mixed laboratory samples (all p < 0.001). Among field samples, the mean fraction of resistant Pfdhps alleles was 4.7% 431V, 95.9% 436F/A, 49.9% 437G, 0.0% 540E, 1.2% 581G and 1.5% 613S/T; corresponding prevalences were 50.0%, 100%, 72.5%, 0.0%, 25.0%, and 12.5%, respectively. The mean fraction of resistant Pfdhfr alleles was 0.6% 16V, 11.1% 50R, 89.0% 51I, 98.3% 59R, 74.7% 108T/N, 8.6% 140L and 8.7% 164L; corresponding prevalences were 12.5%, 75.0%, 100%, 100%, 100%, 50.0%, and 28.6%, respectively. We identified two new point mutations on the Pfdhps gene at codons D484T and D545N.

INTERPRETATION

Computational deconvolution of sequencing chromatograms can discriminate varying proportions of antimalarial drug-sensitive versus -resistant alleles. This cost-effective and quantitative variant-sequencing approach will be useful for population-based surveys that characterize mixed infections at the individual level to survey known and unknown mutations in P. falciparum drug-resistance genes.

FUNDING

This work was supported by the Division of Intramural Research of the National Institute of Allergy and Infectious Diseases, National Institutes of Health (NIH). HM was supported by the African Postdoctoral Training Initiative (APTI) Fellowship program jointly managed by the US NIH, The African Academy of Sciences (AAS) and Bill & Melinda Gates Foundation (BMGF); Grant Reference Number: APTI-18-01.

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.

The problem of antimalarial resistance and its implications for drug discovery

INTRODUCTION

Malaria remains a devastating infectious disease with hundreds of thousands of casualties each year. Antimalarial drug resistance has been a threat to malaria control and elimination for many decades and is still of concern today. Despite the continued effectiveness of current first-line treatments, namely artemisinin-based combination therapies, the emergence of drug-resistant parasites in Southeast Asia and even more alarmingly the occurrence of resistance mutations in Africa is of great concern and requires immediate attention.

AREAS COVERED

A comprehensive overview of the mechanisms underlying the acquisition of drug resistance in Plasmodium falciparum is given. Understanding these processes provides valuable insights that can be harnessed for the development and selection of novel antimalarials with reduced resistance potential. Additionally, strategies to mitigate resistance to antimalarial compounds on the short term by using approved drugs are discussed.

EXPERT OPINION

While employing strategies that utilize already approved drugs may offer a prompt and cost-effective approach to counter antimalarial drug resistance, it is crucial to recognize that only continuous efforts into the development of novel antimalarial drugs can ensure the successful treatment of malaria in the future. Incorporating resistance propensity assessment during this developmental process will increase the likelihood of effective and enduring malaria treatments.

Temporal and spatial dynamics of Plasmodium falciparum clonal lineages in Guyana

Plasmodium parasites, the causal agents of malaria, are eukaryotic organisms that obligately undergo sexual recombination within mosquitoes. However, in low transmission settings where most mosquitoes become infected with only a single parasite clone, parasites recombine with themselves, and the clonal lineage is propagated rather than broken up by outcrossing. We investigated whether stochastic/neutral factors drive the persistence and abundance of Plasmodium falciparum clonal lineages in Guyana, a country with relatively low malaria transmission, but the only setting in the Americas in which an important artemisinin resistance mutation (pfk13 C580Y) has been observed. To investigate whether this clonality was potentially associated with the persistence and spatial spread of the mutation, we performed whole genome sequencing on 1,727 Plasmodium falciparum samples collected from infected patients across a five-year period (2016-2021). We characterized the relatedness between each pair of monoclonal infections (n=1,409) through estimation of identity by descent (IBD) and also typed each sample for known or candidate drug resistance mutations. A total of 160 clones (mean IBD ≥ 0.90) were circulating in Guyana during the study period, comprising 13 highly related clusters (mean IBD ≥ 0.40). In the five-year study period, we observed a decrease in frequency of a mutation associated with artemisinin partner drug (piperaquine) resistance (pfcrt C350R) and limited co-occurence of pfcrt C350R with duplications of plasmepsin 2/3, an epistatic interaction associated with piperaquine resistance. We additionally report polymorphisms exhibiting evidence of selection for drug resistance or other phenotypes and reported a novel pfk13 mutation (G718S) as well as 61 nonsynonymous substitutions that increased markedly in frequency. However, P. falciparum clonal dynamics in Guyana appear to be largely driven by stochastic factors, in contrast to other geographic regions. The use of multiple artemisinin combination therapies in Guyana may have contributed to the disappearance of the pfk13 C580Y mutation.

Role of Seasonal Importation and Random Genetic Drift on Selection for Drug-Resistant Genotypes of Plasmodium falciparum in High Transmission Settings

Historically Plasmodium falciparum has followed a pattern of drug resistance first appearing in low transmission settings before spreading to high transmission settings. Several features of low-transmission regions are hypothesized as explanations: higher chance of symptoms and treatment seeking, better treatment access, less within-host competition among clones, and lower rates of recombination. Here, we test whether importation of drug-resistant parasites is more likely to lead to successful emergence and establishment in low-transmission or high-transmission periods of the same epidemiological setting, using a spatial, individual-based stochastic model of malaria and drug-resistance evolution calibrated for Burkina Faso. Upon controlling for the timing of importation of drug-resistant genotypes and examination of key model variables, we found that drug-resistant genotypes imported during the low transmission season were, (1) more susceptible to stochastic extinction due to the action of random genetic drift, and (2) more likely to lead to establishment of drug resistance when parasites are able to survive early stochastic loss due to drift. This implies that rare importation events are more likely to lead to establishment if they occur during a high-transmission season, but that constant importation (e.g., neighboring countries with high levels of resistance) may produce a greater risk during low-transmission periods.

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 (DHA) as a first-line treatment against malaria parasites. 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 multi-drug resistant KEL1/PLA1 lineage isolate (KH004) and a drug susceptible parasite isolated in Malawi (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 four copies of pm2/3. We recovered 104 unique recombinant progeny and examined a targeted set of progeny representing all possible combinations of variants at pfcrt and pm2/3 for detailed analysis of competitive fitness and a range of PPQ susceptibility phenotypes, including PPQ survival assay (PSA), area under the dose-response curve (AUC), and a limited point IC50 (LP-IC50). We find that inheritance of the KH004 pfcrt allele is required for PPQ resistance, whereas copy number variation in pm2/3 further enhances resistance but does not confer resistance in the absence of PPQ-R-associated mutations in pfcrt. Deeper 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 and confirm the critical role of the PfCRT G367C substitution in PPQ resistance.

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.

Chloroquine resistance evolution in Plasmodium falciparum is mediated by the putative amino acid transporter AAT1

Malaria parasites break down host haemoglobin into peptides and amino acids in the digestive vacuole for export to the parasite cytoplasm for growth: interrupting this process is central to the mode of action of several antimalarial drugs. Mutations in the chloroquine (CQ) resistance transporter, pfcrt, located in the digestive vacuole membrane, confer CQ resistance in Plasmodium falciparum, and typically also affect parasite fitness. However, the role of other parasite loci in the evolution of CQ resistance is unclear. Here we use a combination of population genomics, genetic crosses and gene editing to demonstrate that a second vacuolar transporter plays a key role in both resistance and compensatory evolution. Longitudinal genomic analyses of the Gambian parasites revealed temporal signatures of selection on a putative amino acid transporter (pfaat1) variant S258L, which increased from 0% to 97% in frequency between 1984 and 2014 in parallel with the pfcrt1 K76T variant. Parasite genetic crosses then identified a chromosome 6 quantitative trait locus containing pfaat1 that is selected by CQ treatment. Gene editing demonstrated that pfaat1 S258L potentiates CQ resistance but at a cost of reduced fitness, while pfaat1 F313S, a common southeast Asian polymorphism, reduces CQ resistance while restoring fitness. Our analyses reveal hidden complexity in CQ resistance evolution, suggesting that pfaat1 may underlie regional differences in the dynamics of resistance evolution, and modulate parasite resistance or fitness by manipulating the balance between both amino acid and drug transport.

Forward genetic analysis of monensin and diclazuril resistance in Eimeria tenella

Worldwide distributed coccidiosis is caused by infection of both Eimeria species and Cystoisospora in the host intestine and causes huge economic losses to the livestock industry, especially the poultry industry. The control of such diseases relies mainly on chemoprophylaxis with anticoccidials, which has led to a very common drug resistance in this field. However, the genetic mechanisms underlying resistance to many anticoccidial drugs remain unknown. In this study, strains of E. tenella resistant to 250 mg/kg monensin were generated and characterized. Forward genetic approaches based on pooled genome sequencing, including experimental evolution and linkage group selection, were used to locate candidate targets responsible for resistance to monensin and diclazuril in E. tenella. A total of 16 nonsynonymous mutants in protein-coding genes were identified in monensin-resistant strains, and two genomic regions with strong selection signals were also detected in diclazuril-resistant strains. Our study reveals the genetic characterization of the experimental evolution and linkage group selection in Eimeria species, and also provides important information that contributes to the understanding of the molecular mechanism of drug resistance in coccidia.

Plasmodium falciparum Chloroquine-pfcrt Resistant Haplotypes in Brazilian Endemic Areas Four Decades after CQ Withdrawn

(1) Background: Malaria is a public health problem worldwide. Despite global efforts to control it, antimalarial drug resistance remains a great challenge. In 2009, our team identified, for the first time in Brazil, chloroquine (CQ)-susceptible Plasmodium falciparum parasites in isolates from the Brazilian Amazon. The present study extends those observations to include survey samples from 2010 to 2018 from the Amazonas and Acre states for the purpose of tracking pfcrt molecular changes in P. falciparum parasites. (2) Objective: to investigate SNPs in the P. falciparum gene associated with chemoresistance to CQ (pfcrt). (3) Methods: Sixty-six P. falciparum samples from the Amazonas and Acre states were collected from 2010 to 2018 in patients diagnosed at the Reference Research Center for Treatment and Diagnosis of Malaria (CPD-Mal/Fiocruz), FMT-HVD and Acre Health Units. These samples were subjected to PCR and DNA Sanger sequencing to identify mutations in pfcrt (C72S, M74I, N75E, and K76T). (4) Results: Of the 66 P. falciparum samples genotyped for pfcrt, 94% carried CQ-resistant genotypes and only 4 showed a CQ pfcrt sensitive-wild type genotype, i.e., 1 from Barcelos and 3 from Manaus. (5) Conclusion: CQ-resistant P. falciparum populations are fixed, and thus, CQ cannot be reintroduced in malaria falciparum therapy.

Two decades of molecular surveillance in Senegal reveal changes in known drug resistance mutations associated with historical drug use and seasonal malaria chemoprevention

Drug resistance in Plasmodium falciparum is a major threat to malaria control efforts. We analyzed data from two decades (2000-2020) of continuous molecular surveillance of P. falciparum parasite strains in Senegal to determine how historical changes in drug administration policy may have affected parasite evolution. We profiled several known drug resistance markers and their surrounding haplotypes using a combination of single nucleotide polymorphism (SNP) molecular surveillance and whole-genome sequence (WGS) based population genomics. We observed rapid changes in drug resistance markers associated with the withdrawal of chloroquine and introduction of sulfadoxine-pyrimethamine in 2003. We also observed a rapid increase in Pfcrt K76T and decline in Pfdhps A437G starting in 2014, which we hypothesize may reflect changes in resistance or fitness caused by seasonal malaria chemoprevention (SMC). Parasite populations evolve rapidly in response to drug use, and SMC preventive efficacy should be closely monitored.

Drug Resistance Molecular Markers of Plasmodium falciparum and Severity of Malaria in Febrile Children in the Sentinel Site for Malaria Surveillance of Melen in Gabon: Additional Data from the Plasmodium Diversity Network African Network

The objective of this study was to analyze the relationship between the frequency of artemisinin-based combination (ACT) drug resistance molecular markers and clinical forms of P. falciparum malaria and parasitemia. A cross-sectional study was carried out between January and April 2014 at the Operational Clinical Research Unit of Melen in febrile children aged 12 to 240 months with a Plasmodium sp. infection. A total of 3 mL of peripheral blood collected from an EDTA tube was used for leukocyte depletion. DNA mutation detection was performed by next generation sequencing (NGS). A total of 1075 patients were screened for malaria. Among them, 384 had a Plasmodium infection. P. falciparum mono-infection was found in 98.9% of the patients. Pfcrt-326T mutation was found in all isolates, while 37.9% had Pfmdr2-484I mutant allele. The highest median parasite densities were found in patients infected by parasites carrying the CVIET haplotype of the Pfcrt gene. The different genetic profiles found here, and their variations according to clinical and biological signs of severe malaria, are additional arguments for the surveillance of P. falciparum strains.

Resolving drug selection and migration in an inbred South American Plasmodium falciparum population with identity-by-descent analysis

The human malaria parasite Plasmodium falciparum is globally widespread, but its prevalence varies significantly between and even within countries. Most population genetic studies in P. falciparum focus on regions of high transmission where parasite populations are large and genetically diverse, such as sub-Saharan Africa. Understanding population dynamics in low transmission settings, however, is of particular importance as these are often where drug resistance first evolves. Here, we use the Pacific Coast of Colombia and Ecuador as a model for understanding the population structure and evolution of Plasmodium parasites in small populations harboring less genetic diversity. The combination of low transmission and a high proportion of monoclonal infections means there are few outcrossing events and clonal lineages persist for long periods of time. Yet despite this, the population is evolutionarily labile and has successfully adapted to changes in drug regime. Using newly sequenced whole genomes, we measure relatedness between 166 parasites, calculated as identity by descent (IBD), and find 17 distinct but highly related clonal lineages, six of which have persisted in the region for at least a decade. This inbred population structure is captured in more detail with IBD than with other common population structure analyses like PCA, ADMIXTURE, and distance-based trees. We additionally use patterns of intra-chromosomal IBD and an analysis of haplotypic variation to explore past selection events in the region. Two genes associated with chloroquine resistance, crt and aat1, show evidence of hard selective sweeps, while selection appears soft and/or incomplete at three other key resistance loci (dhps, mdr1, and dhfr). Overall, this work highlights the strength of IBD analyses for studying parasite population structure and resistance evolution in regions of low transmission, and emphasizes that drug resistance can evolve and spread in small populations, as will occur in any region nearing malaria elimination.

Genomic approaches for monitoring transmission dynamics of malaria: A case for malaria molecular surveillance in Sub-Saharan Africa

Transmission dynamics is an important indicator for malaria control and elimination. As we move closer to eliminating malaria in Sub-Saharan Africa (sSA), transmission indices with higher resolution (genomic approaches) will complement our current measurements of transmission. Most of the present programmatic knowledge of malaria transmission patterns are derived from assessments of epidemiologic and clinical data, such as case counts, parasitological estimates of parasite prevalence, and Entomological Inoculation Rates (EIR). However, to eliminate malaria from endemic areas, we need to track changes in the parasite population and how they will impact transmission. This is made possible through the evolving field of genomics and genetics, as well as the development of tools for more in-depth studies on the diversity of parasites and the complexity of infections, among other topics. If malaria elimination is to be achieved globally, country-specific elimination activities should be supported by parasite genomic data from regularly collected blood samples for diagnosis, surveillance and possibly from other programmatic interventions. This presents a unique opportunity to track the spread of malaria parasites and shed additional light on intervention efficacy. In this review, various genetic techniques are highlighted along with their significance for an enhanced understanding of transmission patterns in distinct topological settings throughout Sub-Saharan Africa. The importance of these methods and their limitations in malaria surveillance to guide control and elimination strategies, are explored.

In vitro susceptibility profile of Plasmodium falciparum clinical isolates from Ghana to antimalarial drugs and polymorphisms in resistance markers

Drug resistance in Plasmodium falciparum compromises the effectiveness of antimalarial therapy. This study aimed to evaluate the extent of drug resistance in parasites obtained from international travelers returning from Ghana to guide the management of malaria cases. Eighty-two clinical parasite isolates were obtained from patients returning from Ghana in 2016–2018, of which 29 were adapted to continuous in vitro culture. Their geometric mean IC₅₀ values to a panel of 11 antimalarial drugs, assessed using the standard SYBR Green-I drug sensitivity assay, were 2.1, 3.8, 1.0, 2.7, 17.2, 4.6, 8.3, 8.3, 19.6, 55.1, and 11,555 nM for artemether, artesunate, dihydroartemisinin, lumefantrine, mefloquine, piperaquine, naphthoquine, pyronaridine, chloroquine, quinine, and pyrimethamine, respectively. Except for chloroquine and pyrimethamine, the IC₅₀ values for other tested drugs were below the resistance threshold. The mean ring-stage survival assay value was 0.8%, with four isolates exceeding 1%. The mean piperaquine survival assay value was 2.1%, all below 10%. Mutations associated with chloroquine resistance (pfcrt K76T and pfmdr1 N86Y) were scarce, consistent with the discontinuation of chloroquine a decade ago. Instead, the pfmdr1 86N-184F-1246D haplotype was predominant, suggesting selection by the extensive use of artemether-lumefantrine. No mutations in the pfk13 propeller domain were detected. The pfdhfr/pfdhps quadruple mutant IRNGK associated with resistance to sulfadoxine-pyrimethamine reached an 82% prevalence. In addition, five isolates had pfgch1 gene amplification but, intriguingly, increased susceptibilities to pyrimethamine. This study showed that parasites originating from Ghana were susceptible to artemisinins and the partner drugs of artemisinin-based combination therapies. Genotyping drug resistance genes identified the signature of selection by artemether-lumefantrine. Parasites showed substantial levels of resistance to the antifolate drugs. Continuous resistance surveillance is necessary to guide timely changes in drug policy.

Artemisinin resistance and malaria elimination: Where are we now?

The emergence of artemisinin resistance is a major obstacle to the global malaria eradication/elimination programs. Artemisinin is a very fast-acting antimalarial drug and is the most important drug in the treatment of severe and uncomplicated malaria. For the treatment of acute uncomplicated falciparum malaria, artemisinin derivatives are combined with long half-life partner drugs and widely used as artemisinin-based combination therapies (ACTs). Some ACTs have shown decreased efficacy in the Southeast Asian region. Fortunately, artemisinin has an excellent safety profile and resistant infections can still be treated successfully by modifying the ACT. This review describes the pharmacological properties of ACTs, mechanisms of artemisinin resistance and the potential changes needed in the treatment regimens to overcome resistance. The suggested ACT modifications are extension of the duration of the ACT course, alternating use of different ACT regimens, and addition of another antimalarial drug to the standard ACTs (Triple-ACT). Furthermore, a malaria vaccine (e.g., RTS,S vaccine) could be added to mass drug administration (MDA) campaigns to enhance the treatment efficacy and to prevent further artemisinin resistance development. This review concludes that artemisinin remains the most important antimalarial drug, despite the development of drug-resistant falciparum malaria.

Efficacy of three anti-malarial regimens for uncomplicated Plasmodium falciparum malaria in Cambodia, 2009-2011: a randomized controlled trial and brief review

BACKGROUND

Anti-malarial resistance remains an important public health challenge in Cambodia. The effectiveness of three therapies for uncomplicated falciparum malaria was evaluated in Oddar Meanchey province in Northern Cambodia from 2009 to 2011.

METHODS

In this randomized, open-label, parallel group-controlled trial, 211 subjects at least 5 years old with uncomplicated falciparum malaria were treated with 3 days of directly observed therapy: 63 received artesunate-mefloquine (AS/MQ), 77 received dihydroartemisinin-piperaquine (DHA/PPQ), and 71 received atovaquone-proguanil (ATQ/PG). The subjects were followed for 42 days or until recurrent parasitaemia. Genotyping of msp1, msp2, and glurp among individual parasite isolates distinguished recrudescence from reinfection. Pfmdr1 copy number was measured by real-time PCR and half-maximal parasite inhibitory concentrations (IC50) were measured in vitro by 48-h isotopic hypoxanthine incorporation assay.

RESULTS

The per-protocol PCR-adjusted efficacy (95% confidence interval) at 42 days was 80.6% (70.8-90.5%) for AS/MQ, 97.2% (93.3-100%) for DHA/PPQ, and 92.9% (86.1-99.6%) for ATQ/PG. On day 3, 57.9% remained parasitaemic in the AS/MQ and DHA/PPQ arms. At baseline, 46.9% had microscopic Plasmodium falciparum gametocytaemia. Both recurrences in the DHA/PPQ arm lost Pfmdr1 copy number amplification at recrudescence. All four recurrences in the ATQ/PG arm were wild-type for cytochrome bc1. One subject withdrew from the ATQ/PG arm due to drug allergy.

CONCLUSIONS

This study was conducted at the epicentre of substantial multi-drug resistance that emerged soon thereafter. Occurring early in the national transition from AS/MQ to DHA/PPQ, both DHA/PPQ and ATQ/PG had acceptable efficacy against uncomplicated falciparum malaria. However, efficacy of AS/MQ was only 80% with apparent mefloquine resistance based on elevated Pfmdr1 copy number and IC50. By 2009, there was already significant evidence of artemisinin resistance not previously reported at the Northern Cambodia-Thai border. This study suggests the basis for early development of significant DHA/PPQ failures within 3 years of introduction. Artemisinin resistance likely occurred on the Northern border concurrently with that reported along the Western border in Pailin. Trial registration This legacy trial was conducted prior to International Committee of Medical Journal Editors' requirements for preregistration on ClinicalTrials.gov. The full protocol has been provided.

Sources of transcription variation in Plasmodium falciparum

Variation in transcript abundance can contribute to both short-term environmental response and long-term evolutionary adaptation. Most studies are designed to assess differences in mean transcription levels and do not consider other potentially important and confounding sources of transcriptional variation. Detailed quantification of variation sources will improve our ability to detect and identify the mechanisms that contribute to genome-wide transcription changes that underpin adaptive responses. To quantify innate levels of expression variation, we measured mRNA levels for more than 5000 genes in the malaria parasite, Plasmodium falciparum, among clones derived from two parasite strains across biologically and experimentally replicated batches. Using a mixed effects model, we partitioned the total variation among four sources - between strain, within strain, environmental batch effects, and stochastic noise. We found 646 genes with significant variation attributable to at least one of these sources. These genes were categorized by their predominant variation source and further examined using gene ontology enrichment analysis to associate function with each source of variation. Genes with environmental batch effect and within strain transcript variation may contribute to phenotypic plasticity, while genes with between strain variation may contribute to adaptive responses and processes that lead to parasite strain-specific survival under varied conditions.

Malaria chemoprevention and drug resistance: a review of the literature and policy implications

Chemoprevention strategies reduce malaria disease and death, but the efficacy of anti-malarial drugs used for chemoprevention is perennially threatened by drug resistance. This review examines the current impact of chemoprevention on the emergence and spread of drug resistant malaria, and the impact of drug resistance on the efficacy of each of the chemoprevention strategies currently recommended by the World Health Organization, namely, intermittent preventive treatment in pregnancy (IPTp); intermittent preventive treatment in infants (IPTi); seasonal malaria chemoprevention (SMC); and mass drug administration (MDA) for the reduction of disease burden in emergency situations. While the use of drugs to prevent malaria often results in increased prevalence of genetic mutations associated with resistance, malaria chemoprevention interventions do not inevitably lead to meaningful increases in resistance, and even high rates of resistance do not necessarily impair chemoprevention efficacy. At the same time, it can reasonably be anticipated that, over time, as drugs are widely used, resistance will generally increase and efficacy will eventually be lost. Decisions about whether, where and when chemoprevention strategies should be deployed or changed will continue to need to be made on the basis of imperfect evidence, but practical considerations such as prevalence patterns of resistance markers can help guide policy recommendations.

Identification of polymorphisms in genes associated with drug resistance in Plasmodium falciparum isolates from school-age children in Kinshasa, Democratic Republic of Congo

BACKGROUND

The emergence and spread of Plasmodium falciparum parasites resistant to antimalarial drugs constitutes an obstacle to malaria control and elimination. This study aimed to identify the prevalence of polymorphisms in pfk13, pfmdr1, pfdhfr, pfdhps and pfcrt genes in isolates from asymptomatic and symptomatic school-age children in Kinshasa.

METHODS

Nested-PCR followed by sequencing was performed for the detection of pfk13, pfmdr1, pfdhfr, pfdhps and pfcrt polymorphisms.

RESULTS

Two mutations in pfk13, C532S and Q613E were identified in the Democratic Republic of Congo for the first time. The prevalence of the drug-resistance associated mutations pfcrt K76T, pfdhps K540E and pfmdr1 N86Y was low, being 27%, 20% and 9%, respectively.

CONCLUSION

We found a low prevalence of genetic markers associated with chloroquine and sulfadoxine-pyrimethamine resistance in Kinshasa. Furthermore, no mutations previously associated with resistance against artemisinin and is derivatives were observed in the pfK13 gene. These findings support the continued use of ACTs and IPTp-SP. Continuous molecular monitoring of antimalarial resistance markers is recommended.

Population-level genome-wide STR discovery and validation for population structure and genetic diversity assessment of Plasmodium species

Short tandem repeats (STRs) are highly informative genetic markers that have been used extensively in population genetics analysis. They are an important source of genetic diversity and can also have functional impact. Despite the availability of bioinformatic methods that permit large-scale genome-wide genotyping of STRs from whole genome sequencing data, they have not previously been applied to sequencing data from large collections of malaria parasite field samples. Here, we have genotyped STRs using HipSTR in more than 3,000 Plasmodium falciparum and 174 Plasmodium vivax published whole-genome sequence data from samples collected across the globe. High levels of noise and variability in the resultant callset necessitated the development of a novel method for quality control of STR genotype calls. A set of high-quality STR loci (6,768 from P. falciparum and 3,496 from P. vivax) were used to study Plasmodium genetic diversity, population structures and genomic signatures of selection and these were compared to genome-wide single nucleotide polymorphism (SNP) genotyping data. In addition, the genome-wide information about genetic variation and other characteristics of STRs in P. falciparum and P. vivax have been available in an interactive web-based R Shiny application PlasmoSTR (https://github.com/bahlolab/PlasmoSTR).

Antimalarial Natural Products

Natural products have made a crucial and unique contribution to human health, and this is especially true in the case of malaria, where the natural products quinine and artemisinin and their derivatives and analogues, have saved millions of lives. The need for new drugs to treat malaria is still urgent, since the most dangerous malaria parasite, Plasmodium falciparum, has become resistant to quinine and most of its derivatives and is becoming resistant to artemisinin and its derivatives. This volume begins with a short history of malaria and follows this with a summary of its biology. It then traces the fascinating history of the discovery of quinine for malaria treatment and then describes quinine's biosynthesis, its mechanism of action, and its clinical use, concluding with a discussion of synthetic antimalarial agents based on quinine's structure. The volume then covers the discovery of artemisinin and its development as the source of the most effective current antimalarial drug, including summaries of its synthesis and biosynthesis, its mechanism of action, and its clinical use and resistance. A short discussion of other clinically used antimalarial natural products leads to a detailed treatment of other natural products with significant antiplasmodial activity, classified by compound type. Although the search for new antimalarial natural products from Nature's combinatorial library is challenging, it is very likely to yield new antimalarial drugs. The chapter thus ends by identifying over ten natural products with development potential as clinical antimalarial agents.

Genetic surveillance for monitoring the impact of drug use on Plasmodium falciparum populations

The use of antimalarial drugs is an effective strategy in the fight against malaria. However, selection of drug resistant parasites is a constant threat to the continued use of this approach. Antimalarial drugs are used not only to treat infections but also as part of population-level strategies to reduce malaria transmission toward elimination. While there is strong evidence that the ongoing use of antimalarial drugs increases the risk of the emergence and spread of drug-resistant parasites, it is less clear how population-level use of drug-based interventions like seasonal malaria chemoprevention (SMC) or mass drug administration (MDA) may contribute to drug resistance or loss of drug efficacy. Critical to sustained use of drug-based strategies for reducing the burden of malaria is the surveillance of population-level signals related to transmission reduction and resistance selection. Here we focus on Plasmodium falciparum and discuss the genetic signatures of a parasite population that are correlated with changes in transmission and related to drug pressure and resistance as a result of drug use. We review the evidence for MDA and SMC contributing to malaria burden reduction and drug resistance selection and examine the use and impact of these interventions in Senegal. Throughout we consider best strategies for ongoing surveillance of both population and resistance signals in the context of different parasite population parameters. Finally, we propose a roadmap for ongoing surveillance during population-level drug-based interventions to reduce the global malaria burden.

Comparative Analysis of Plasmodium falciparum Genotyping via SNP Detection, Microsatellite Profiling, and Whole-Genome Sequencing

Research efforts to combat antimalarial drug resistance rely on quick, robust, and sensitive methods to genetically characterize Plasmodium falciparum parasites. We developed a single-nucleotide polymorphism (SNP)-based genotyping method that can assess 33 drug resistance-conferring SNPs in dhfr, dhps, pfmdr1, pfcrt, and k13 in nine PCRs, performed directly from P. falciparum cultures or infected blood. We also optimized multiplexed fragment analysis and gel electrophoresis-based microsatellite typing methods using a set of five markers that can distinguish 12 laboratory strains of diverse geographical and temporal origin. We demonstrate how these methods can be applied to screen for the multidrug-resistant KEL1/PLA1/PfPailin (KelPP) lineage that has been sweeping across the Greater Mekong Subregion, verify parasite in vitro SNP-editing, identify novel recombinant genetic cross progeny, or cluster strains to infer their geographical origins. Results were compared with Illumina-based whole-genome sequence analysis that provides the most detailed sequence information but is cost-prohibitive. These adaptable, simple, and inexpensive methods can be easily implemented into routine genotyping of P. falciparum parasites in both laboratory and field settings.

Polymeric Nanoparticle Based Diagnosis and Nanomedicine for Treatment and Development of Vaccines for Cerebral Malaria: A Review on Recent Advancement

Cerebral malaria occurs due to Plasmodium falciparum infection, which causes 228 million infections and 450,000 deaths worldwide every year. African people are mostly affected with nearly 91% cases, of which 86% are pregnant women and infants. India and Brazil are the other two countries severely suffering from malaria endemicity. Commonly used drugs have severe side effects, and unfortunately no suitable vaccine is available in the market today. In this line, this review is focused on polymeric nanomaterials and nanocapsules that can be used for the development of effective diagnostic strategies, nanomedicines, and vaccines in the management of cerebral malaria. Further, this review will help scientists and medical professionals by updating the status on the development stages of polymeric nanoparticle based diagnostics, nanomedicines, and vaccines and strategies to eradicate cerebral malaria. In addition to this, the predominant focus of this review is antimalarial agents based on polymer nanomedicines that are currently in the preclinical and clinical trial stages, and potential developments are suggested as well. This review further will have an important social and commercial impact worldwide for the development of polymeric nanomedicines and strategies for the treatment of cerebral malaria.

Arterolane-piperaquine-mefloquine versus arterolane-piperaquine and artemether-lumefantrine in the treatment of uncomplicated Plasmodium falciparum malaria in Kenyan children: a single-centre, open-label, randomised, non-inferiority trial

BACKGROUND

Triple antimalarial combination therapies combine potent and rapidly cleared artemisinins or related synthetic ozonides, such as arterolane, with two, more slowly eliminated partner drugs to reduce the risk of resistance. We aimed to assess the safety, tolerability, and efficacy of arterolane-piperaquine-mefloquine versus arterolane-piperaquine and artemether-lumefantrine for the treatment of uncomplicated falciparum malaria in Kenyan children.

METHODS

In this single-centre, open-label, randomised, non-inferiority trial done in Kilifi County Hospital, Kilifi, coastal Kenya, children with uncomplicated Plasmodium falciparum malaria were recruited. Eligible patients were aged 2-12 years and had an asexual parasitaemia of 5000-250 000 parasites per μL. The exclusion criteria included the presence of an acute illness other than malaria, the inability to tolerate oral medications, treatment with an artemisinin derivative in the previous 7 days, a known hypersensitivity or contraindication to any of the study drugs, and a QT interval corrected for heart rate (QTc interval) longer than 450 ms. Patients were randomly assigned (1:1:1), by use of blocks of six, nine, and 12, and opaque, sealed, and sequentially numbered envelopes, to receive either arterolane-piperaquine, arterolane-piperaquine-mefloquine, or artemether-lumefantrine. Laboratory staff, but not the patients, the patients' parents or caregivers, clinical or medical officers, nurses, or trial statistician, were masked to the intervention groups. For 3 days, oral artemether-lumefantrine was administered twice daily (target dose 5-24 mg/kg of bodyweight of artemether and 29-144 mg/kg of bodyweight of lumefantrine), and oral arterolane-piperaquine (arterolane dose 4 mg/kg of bodyweight; piperaquine dose 20 mg/kg of bodyweight) and oral arterolane-piperaquine-mefloquine (mefloquine dose 8 mg/kg of bodyweight) were administered once daily. All patients received 0·25 mg/kg of bodyweight of oral primaquine at hour 24. All patients were admitted to Kilifi County Hospital for at least 3 consecutive days and followed up at day 7 and, thereafter, weekly for up to 42 days. The primary endpoint was 42-day PCR-corrected efficacy, defined as the absence of treatment failure in the first 42 days post-treatment, of arterolane-piperaquine-mefloquine versus artemether-lumefantrine, and, along with safety, was analysed in the intention-to-treat population, which comprised all patients who received at least one dose of a study drug. The 42-day PCR-corrected efficacy of arterolane-piperaquine-mefloquine versus arterolane-piperaquine was an important secondary endpoint and was also analysed in the intention-to-treat population. The non-inferiority margin for the risk difference between treatments was -7%. The study is registered in ClinicalTrials.gov, NCT03452475, and is completed.

FINDINGS

Between March 7, 2018, and May 2, 2019, 533 children with P falciparum were screened, of whom 217 were randomly assigned to receive either arterolane-piperaquine (n=73), arterolane-piperaquine-mefloquine (n=72), or artemether-lumefantrine (n=72) and comprised the intention-to-treat population. The 42-day PCR-corrected efficacy after treatment with arterolane-piperaquine-mefloquine (100%, 95% CI 95-100; 72/72) was non-inferior to that after treatment with artemether-lumefantrine (96%, 95% CI 88-99; 69/72; risk difference 4%, 95% CI 0-9; p=0·25). The 42-day PCR-corrected efficacy of arterolane-piperaquine-mefloquine was non-inferior to that of arterolane-piperaquine (100%, 95% CI 95-100; 73/73; risk difference 0%). Vomiting rates in the first hour post-drug administration were significantly higher in patients treated with arterolane-piperaquine (5%, 95% CI 2-9; ten of 203 drug administrations; p=0·0013) or arterolane-piperaquine-mefloquine (5%, 3-9; 11 of 209 drug administrations; p=0·0006) than in patients treated with artemether-lumefantrine (1%, 0-2; three of 415 drug administrations). Upper respiratory tract complaints (n=26 for artemether-lumefantrine; n=19 for arterolane-piperaquine-mefloquine; n=23 for arterolane-piperaquine), headache (n=13; n=4; n=5), and abdominal pain (n=7; n=5; n=5) were the most frequently reported adverse events. There were no deaths.

INTERPRETATION

This study shows that arterolane-piperaquine-mefloquine is an efficacious and safe treatment for uncomplicated falciparum malaria in children and could potentially be used to prevent or delay the emergence of antimalarial resistance.

FUNDING

UK Department for International Development, The Wellcome Trust, The Bill & Melinda Gates Foundation, Sun Pharmaceutical Industries.

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.

Fixation of pfcrt chloroquine resistance alleles in Plasmodium falciparum clinical isolates collected from unrest tribal agencies of Pakistan

Plasmodium falciparum resistance to Chloroquine (CQ) is a significant cause of mortality and morbidity worldwide. There is a paucity of documented data on the prevalence of CQ-resistant mutant haplotypes of Pfcrt and Pfmdr1 genes from malaria-endemic war effected Federally Administered Tribal Areas of Pakistan. The objective of this study was to investigate the prevalence of P. falciparum CQ-resistance in this area. Clinical isolates were collected between May 2017 and May 2018 from North Waziristan and South Waziristan agencies of Federally Administrated Trial Area. Subsequently, Giemsa-stained blood smears were examined to detect Plasmodium falciparum. Extraction of malarial DNA was done from microscopy positive P. falciparum samples, and P. falciparum infections were confirmed by nested PCR (targeting Plasmodium small subunit ribosomal ribonucleic acid (ssrRNA) genes). All PCR confirmed P. falciparum samples were sequenced by pyrosequencing to find out mutation in Pfcrt gene at codon K76T and in pfmdr1 at codons N86Y, Y184F, N1042D, and D1246Y. Out of 121 microscopies positive P. falciparum cases, 109 samples were positive for P. falciparum by nested PCR. Pfcrt K76T mutation was found in 96% of isolates, Pfmdr1 N86Y mutation was observed in 20%, and 11% harboured Y184F mutation. All samples were wild type for Pfmdr1 codon N1042D and D1246Y. In the FATA, Pakistan, the frequency of resistant allele 76T remained high despite the removal of CQ. However, current findings of the study suggest complete fixation of P. falciparum CQ-resistant genotype in the study area.

Current and emerging strategies to combat antimalarial resistance

INTRODUCTION

Since the spread of chloroquine resistance in Plasmodium falciparum in the 1960s, recommendations have been made on how to respond to antimalarial resistance. Only with the advent of artemisinin partial resistance were large scale efforts made in the Greater Mekong Subregion to carry out recommendations in a coordinated and well-funded manner. Independent emergence of parasites partially resistant to artemisinins has now been reported in Rwanda.

AREAS COVERED

We reviewed past recommendations and activities to respond to resistance as well as the research ongoing into new ways to stop or delay the spread of resistant parasites.

EXPERT OPINION

Inadequate information limits the options and support for a strong, coordinated response to artemisinin partial resistance in Africa, making better phenotypic and genotypic surveillance a priority. A response to resistance needs to address factors that may have hastened the emergence and could speed the spread, including overuse of drugs and lack of access to quality treatment. New ways to use the existing treatments in the response to resistance such as multiple first-lines are currently impeded by the limited number of drugs available.

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.

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.

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.

An open dataset of Plasmodium falciparum genome variation in 7,000 worldwide samples

MalariaGEN is a data-sharing network that enables groups around the world to work together on the genomic epidemiology of malaria. Here we describe a new release of curated genome variation data on 7,000 Plasmodium falciparum samples from MalariaGEN partner studies in 28 malaria-endemic countries. High-quality genotype calls on 3 million single nucleotide polymorphisms (SNPs) and short indels were produced using a standardised analysis pipeline. Copy number variants associated with drug resistance and structural variants that cause failure of rapid diagnostic tests were also analysed.  Almost all samples showed genetic evidence of resistance to at least one antimalarial drug, and some samples from Southeast Asia carried markers of resistance to six commonly-used drugs. Genes expressed during the mosquito stage of the parasite life-cycle are prominent among loci that show strong geographic differentiation. By continuing to enlarge this open data resource we aim to facilitate research into the evolutionary processes affecting malaria control and to accelerate development of the surveillance toolkit required for malaria elimination.

An open dataset of Plasmodium falciparum genome variation in 7,000 worldwide samples

MalariaGEN is a data-sharing network that enables groups around the world to work together on the genomic epidemiology of malaria. Here we describe a new release of curated genome variation data on 7,000 Plasmodium falciparum samples from MalariaGEN partner studies in 28 malaria-endemic countries. High-quality genotype calls on 3 million single nucleotide polymorphisms (SNPs) and short indels were produced using a standardised analysis pipeline. Copy number variants associated with drug resistance and structural variants that cause failure of rapid diagnostic tests were also analysed.  Almost all samples showed genetic evidence of resistance to at least one antimalarial drug, and some samples from Southeast Asia carried markers of resistance to six commonly-used drugs. Genes expressed during the mosquito stage of the parasite life-cycle are prominent among loci that show strong geographic differentiation. By continuing to enlarge this open data resource we aim to facilitate research into the evolutionary processes affecting malaria control and to accelerate development of the surveillance toolkit required for malaria elimination.

Theory of change: Drama and arts-based community engagement for malaria research and elimination in Cambodia

Background: Across the Greater Mekong Sub-region, malaria persists in isolated communities along international borders. Arts and drama have been used to reach to communities in Cambodia to engage them in malaria research, prevention and control. The “Village Drama Against Malaria” (VDAM) project was conducted in north eastern and western Cambodia: Stung Treng; Battambang and Pailin provinces during 2016 to 2019.  In total, VDAM reached 55 rural villages, 2,378 student participants and 43,502 audience members.

Methods: This article presents the results of two stakeholder-led evaluation workshops in which participants collaboratively developed theories of change to better understand the potential and actual impact of arts and drama-based activities on malaria in these communities. The workshops had a particular focus on identifying areas for monitoring and evaluation so that impact can be measured. Workshop participants included village malaria workers, community leaders, professional and student drama performers, and representatives from the local health authorities and the national malaria control programme.

Results: Five broad areas were identified as relevant for monitoring and evaluation: logistical and practical challenges; embeddedness and reach of engagement; health knowledge and confidence of young people; effectiveness of communications; impact on malaria. These areas align well with the monitoring and evaluation conducted to date and point to additional opportunities for data collection.

Conclusions: The findings from these workshops will inform future engagement strategies, for example, we may engage a smaller number of young people but over a longer period and more in-depth.

Theory of change: Drama and arts-based community engagement for malaria research and elimination in Cambodia

Background: Across the Greater Mekong Sub-region, malaria persists in isolated communities along international borders. Arts and drama have been used to reach to communities in Cambodia to engage them in malaria research, prevention and control. The "Village Drama Against Malaria" (VDAM) project was conducted in north eastern and western Cambodia: Stung Treng; Battambang and Pailin provinces during 2016 to 2019.  In total, VDAM reached 55 rural villages, 2,378 student participants and 43,502 audience members. Methods: This article presents the results of two stakeholder-led evaluation workshops in which participants collaboratively developed theories of change to better understand the potential and actual impact of arts and drama-based activities on malaria in these communities. The workshops had a particular focus on identifying areas for monitoring and evaluation so that impact can be measured. Workshop participants included village malaria workers, community leaders, professional and student drama performers, and representatives from the local health authorities and the national malaria control programme. Results: Five broad areas were identified as relevant for monitoring and evaluation: logistical and practical challenges; embeddedness and reach of engagement; health knowledge and confidence of young people; effectiveness of communications; impact on malaria. These areas align well with the monitoring and evaluation conducted to date and point to additional opportunities for data collection. Conclusions: The findings from these workshops will inform future engagement strategies, for example, we may engage a smaller number of young people but over a longer period and more in-depth.

In Vitro Susceptibility of Plasmodium falciparum Isolates from the China-Myanmar Border Area to Piperaquine and Association with Candidate Markers

Plasmodium falciparum from the Greater Mekong subregion has evolved resistance to the artemisinin-based combination therapy dihydroartemisinin and the partner drug piperaquine. To monitor the potential westward spread or independent evolution of piperaquine resistance, we evaluated the in vitro susceptibility of 120 P. falciparum isolates collected at the China-Myanmar border during 2007-2016. The parasite isolates displayed a relatively wide range of piperaquine susceptibility estimates. While 56.7% of the parasites showed bimodal drug response curves, all but five generated area-under-the-curve (AUC) estimates consistent with a susceptible phenotype. Using the piperaquine survival assay (PSA), 5.6% parasites showed reduced susceptibility. Of note, parasites from 2014-2016 showed the highest AUC value and the highest proportion with a bimodal curve, suggesting falling effectiveness in these later years. Unsupervised K-mean analysis of the combined data assigned parasites into three clusters and identified significant correlations between IC₅₀, IC_90, and AUC values. No parasites carried the E415G mutation in a putative exo-nuclease, new mutations in PfCRT, or amplification of the plasmepsin 2/3 genes, suggesting mechanisms of reduced piperaquine susceptibility that differ from those described in other countries of the region. The association of increased AUC, IC₅₀, and IC₉₀ values with major PfK13 mutations (F446I and G533S) suggests that piperaquine resistance may evolve in these PfK13 genetic backgrounds. Additionally, the Pfmdr1 F1226Y mutation was associated with significantly higher PSA values. Further elucidation of piperaquine resistance mechanisms and continuous surveillance are warranted.

Association of Plasmodium falciparum kelch13 R561H genotypes with delayed parasite clearance in Rwanda: an open-label, single-arm, multicentre, therapeutic efficacy study

BACKGROUND

Partial artemisinin resistance is suspected if delayed parasite clearance (ie, persistence of parasitaemia on day 3 after treatment initiation) is observed. Validated markers of artemisinin partial resistance in southeast Asia, Plasmodium falciparum kelch13 (Pfkelch13) R561H and P574L, have been reported in Rwanda but no association with parasite clearance has been observed. We aimed to establish the efficacy of artemether-lumefantrine and genetic characterisation of Pfkelch13 alleles and their association with treatment outcomes.

METHODS

This open-label, single-arm, multicentre, therapeutic efficacy study was done in 2018 in three Rwandan sites: Masaka, Rukara, and Bugarama. Children aged 6-59 months with P falciparum monoinfection and fever were eligible and treated with a 3-day course of artemether-lumefantrine. Treatment response was monitored for 28 days using weekly microscopy screenings of blood samples for P falciparum. Mutations in Pfkelch13 and P falciparum multidrug resistance-1 (Pfmdr1) genes were characterised in parasites collected from enrolled participants. Analysis of flanking microsatellites surrounding Pfkelch13 was done to define the origins of the R561H mutations. The primary endpoint was PCR-corrected parasitological cure on day 28, as per WHO protocol.

FINDINGS

228 participants were enrolled and 224 (98·2%) reached the study endpoint. PCR-corrected efficacies were 97·0% (95% CI 88-100) in Masaka, 93·8% (85-98) in Rukara, and 97·2% (91-100) in Bugarama. Pfkelch13 R561H mutations were present in 28 (13%) of 218 pre-treatment samples and P574L mutations were present in two (1%) pre-treatment samples. 217 (90%) of the 240 Pfmdr1 haplotypes observed in the pretreatment samples, had either the NFD (N86Y, Y184F, D1246Y) or NYD haplotype. Eight (16%) of 51 participants in Masaka and 12 (15%) of 82 participants in Rukara were microscopically positive 3 days after treatment initiation, which was associated with pre-treatment presence of Pfkelch13 R561H in Masaka (p=0·0005). Genetic analysis of Pfkelch13 R561H mutations suggest their common ancestry and local origin in Rwanda.

INTERPRETATION

We confirm evidence of emerging artemisinin partial resistance in Rwanda. Although artemether-lumefantrine remains efficacious, vigilance for decreasing efficacy, further characterisation of artemisinin partial resistance, and evaluation of additional antimalarials in Rwanda should be considered.

FUNDING

The US President's Malaria Initiative.

TRANSLATION

For the French translation of the abstract see Supplementary Materials section.

Historical trends and new surveillance of Plasmodium falciparum drug resistance markers in Angola

BACKGROUND

Plasmodium falciparum resistance to chloroquine (CQ) and sulfadoxine-pyrimethamine (SP) has historically posed a major threat to malaria control throughout the world. The country of Angola officially replaced CQ with artemisinin-based combination therapy (ACT) as a first-line treatment in 2006, but malaria cases and deaths have recently been rising. Many classic resistance mutations are relevant for the efficacy of currently available drugs, making it important to continue monitoring their frequency in Angola.

METHODS

Plasmodium falciparum DNA was sampled from the blood of 50 hospital patients in Cabinda, Angola from October-December of 2018. Each infection was genotyped for 13 alleles in the genes crt, mdr1, dhps, dhfr, and kelch13, which are collectively involved in resistance to six common anti-malarials. To compare frequency patterns over time, P. falciparum genotype data were also collated from studies published from across Angola in the last two decades.

RESULTS

The two most important alleles for CQ resistance, crt 76T and mdr1 86Y, were found at respective frequencies of 71.4% and 6.5%. Historical data suggest that mdr1 N86 has been steadily replacing 86Y throughout Angola in the last decade, while the frequency of crt 76T has been more variable across studies. Over a third of new samples from Cabinda were 'quintuple mutants' for SP resistance in dhfr/dhps, with a sixth mutation at dhps A581G present at 9.6% frequency. The markers dhfr 51I, dhfr 108N, and dhps 437G have been nearly fixed in Angola since the early 2000s, whereas dhfr 59R may have risen to high frequency more recently. Finally, no non-synonymous polymorphisms were detected in kelch13, which is involved in artemisinin resistance in Southeast Asia.

CONCLUSIONS

Genetic markers of P. falciparum resistance to CQ are likely declining in frequency in Angola, consistent with the official discontinuation of CQ in 2006. The high frequency of multiple genetic markers of SP resistance is consistent with the continued public and private use of SP. In the future, more complete haplotype data from mdr1, dhfr, and dhps will be critical for understanding the changing efficacy of multiple anti-malarial drugs. These data can be used to support effective drug policy decisions in Angola.

Treatment of COVID-19 with Chloroquine: Implication for Malaria Chemotherapy Using ACTs in Disease Endemic Countries

Based on reports of parasite resistance and on World Health Organization recommendation, chloroquine was replaced with the artemisinin-based combination therapies (ACTs) as the first choice of drugs for the treatment of uncomplicated malaria. Disuse of chloroquine led to restoration of drug-sensitive parasite to some extent in certain countries. Ever since chloroquine and hydroxychloroquine were touted as potential treatment for coronavirus disease 2019 (COVID-19), there has been a dramatic surge in demand for the drugs. Even in areas where chloroquine is proscribed, there has been an unexpected increase in demand and supply of the drug. This situation is quite worrying as the indiscriminate use of chloroquine may produce drug-resistant parasites which may impact negatively on the efficacy of amodiaquine due to cross-resistance. Amodiaquine is a partner drug in one of the ACTs and in some of the drugs used for intermittent preventive treatment. We herein discuss the consequences of the escalated use of chloroquine in the management of COVID-19 on chemotherapy or chemoprevention of malaria and offer an advice. We speculate that parasite strains resistant to chloroquine will escalate due to the increased and indiscriminate use of the drug and consequently lead to cross-resistance with amodiaquine which is present in some drug schemes aforementioned. Under the circumstance, the anticipated hope of reverting to the use of the 'resurrected chloroquine' to manage malaria in future is likely to diminish. The use of chloroquine and its derivatives for the management of COVID-19 should be controlled.

Restricted genetic heterogeneity of the Plasmodium vivax transmission-blocking vaccine (TBV) candidate Pvs48/45 in a low transmission setting: Implications for the Plasmodium vivax malaria vaccine development

Plasmodium vivax is the most widespread malaria species parasitizing humans outside Africa, with approximately 100 million cases reported per year. Most human cases of P. vivax are asymptomatic with low parasitemia, making active case detection-based elimination programme challenging and less effective. Despite the widespread distribution of P. vivax, no effective vaccines are currently available. Transmission blocking vaccines have recently emerged as potential vaccine candidates to reduce transmission rates to below the essential levels required for the maintenance of the parasite life cycle. Here, we demonstrated that P. vivax was the predominant species found in a malaria-endemic area, although P. vivax/P. falciparum co-infections were also common. Through genomic sequence analysis and neighbor-joining algorithms, we demonstrated limited genetic heterogeneity in the P. vivax transmission-blocking vaccine candidate Pvs48/45 among clinical isolates of P. vivax. Restricted genetic polymorphism occurred at both nucleotide and amino acid levels. The most frequent mutation was A → G at nucleotide position 77 (46.7%), whereas the least frequent was C → T at nucleotide position 1230 (3.3%). The occurrence of single nucleotide polymorphisms (SNPs) distribution at 6/8 positions (75%) led to changes in amino acid sequences in the Pvs48/45 loci, whereas 2/8 (25%) of SNPs resulted in no amino acid sequence variations. Consistently, the nucleotide diversity in the Pvs48/45 locus among the P. vivax population studied was extremely low (π = 0.000525). Changes in amino acid sequences in the Pvs48/45 protein did not result in substantial conformational modifications in the tertiary structures of these proteins. Unveiling the population genetic structure and genetic heterogeneity of vaccine target antigens are necessary for rational design of transmission-blocking antibody vaccines and to monitor the vaccine efficacy in clinical trials in endemic areas for malaria.

Association of variants in IL1B, TLR9, TREM1, IL10RA, and CD3G and Native American ancestry on malaria susceptibility in Colombian populations

Host genetics is an influencing factor in the manifestation of infectious diseases. In this study, the association of mild malaria with 28 variants in 16 genes previously reported in other populations and/or close to ancestry-informative markers (AIMs) selected was evaluated in an admixed 736 Colombian population sample. Additionally, the effect of genetic ancestry on phenotype expression was explored. For this purpose, the ancestral genetic composition of Turbo and El Bagre was determined. A higher Native American ancestry trend was found in the population with lower malaria susceptibility [odds ratio (OR) = 0.416, 95% confidence interval (95% CI) = 0.234-0.740, P = 0.003]. Three AIMs presented significant associations with the disease phenotype (MID1752, MID921, and MID1586). The first two were associated with greater malaria susceptibility (D/D, OR = 2.23, 95% CI = 1.06-4.69, P = 0.032 and I/D-I/I, OR = 2.14, 95% CI = 1.18-3.87, P = 0.011, respectively), and the latter has a protective effect on the appearance of malaria (I/I, OR = 0.18, 95% CI = 0.08-0.40, P < 0.0001). After adjustment by age, sex, municipality, and genetic ancestry, genotype association analysis showed evidence of association with malaria susceptibility for variants in or near IL1B, TLR9, TREM1, IL10RA, and CD3G genes: rs1143629-IL1B (G/A-A/A, OR = 0.41, 95% CI = 0.21-0.78, P = 0.0051), rs352139-TLR9 (T/T, OR = 0.28, 95% CI = 0.11-0.72, P = 0.0053), rs352140-TLR9 (C/C, OR = 0.41, 95% CI = 0.20-0.87, P = 0.019), rs2234237-TREM1 (T/A-A/A, OR = 0.43, 95% CI = 0.23-0.79, P = 0.0056), rs4252246-IL10RA (C/A-A/A, OR = 2.11, 95% CI = 1.18-3.75, P = 0.01), and rs1561966-CD3G (A/A, OR = 0.20, 95% CI = 0.06-0.69, P = 0.0058). The results showed the participation of genes involved in immunological processes and suggested an effect of ancestral genetic composition over the traits analyzed. Compared to the paisa population (Antioquia), Turbo and El Bagre showed a strong decrease in European ancestry and an increase in African and Native American ancestries. Also, a novel association of two single nucleotide polymorphisms with malaria susceptibility was identified in this study.

Emergence of artemisinin-resistant Plasmodium falciparum with kelch13 C580Y mutations on the island of New Guinea

The rapid and aggressive spread of artemisinin-resistant Plasmodium falciparum carrying the C580Y mutation in the kelch13 gene is a growing threat to malaria elimination in Southeast Asia, but there is no evidence of their spread to other regions. We conducted cross-sectional surveys in 2016 and 2017 at two clinics in Wewak, Papua New Guinea (PNG) where we identified three infections caused by C580Y mutants among 239 genotyped clinical samples. One of these mutants exhibited the highest survival rate (6.8%) among all parasites surveyed in ring-stage survival assays (RSA) for artemisinin. Analyses of kelch13 flanking regions, and comparisons of deep sequencing data from 389 clinical samples from PNG, Indonesian Papua and Western Cambodia, suggested an independent origin of the Wewak C580Y mutation, showing that the mutants possess several distinctive genetic features. Identity by descent (IBD) showed that multiple portions of the mutants' genomes share a common origin with parasites found in Indonesian Papua, comprising several mutations within genes previously associated with drug resistance, such as mdr1, ferredoxin, atg18 and pnp. These findings suggest that a P. falciparum lineage circulating on the island of New Guinea has gradually acquired a complex ensemble of variants, including kelch13 C580Y, which have affected the parasites' drug sensitivity. This worrying development reinforces the need for increased surveillance of the evolving parasite populations on the island, to contain the spread of resistance.

Expansion of a Specific Plasmodium falciparum PfMDR1 Haplotype in Southeast Asia with Increased Substrate Transport

Global efforts to eliminate malaria depend on the continued success of artemisinin-based combination therapies (ACTs) that target Plasmodium asexual blood-stage parasites. Resistance to ACTs, however, has emerged, creating the need to define the underlying mechanisms. Mutations in the P. falciparum multidrug resistance protein 1 (PfMDR1) transporter constitute an important determinant of resistance. Applying gene editing tools combined with an analysis of a public database containing thousands of parasite genomes, we show geographic selection and expansion of a pfmdr1 gene amplification encoding the N86/184F haplotype in Southeast Asia. Parasites expressing this PfMDR1 variant possess a higher transport capacity that modulates their responses to antimalarials. These data could help tailor and optimize antimalarial drug usage in different regions where malaria is endemic by taking into account the regional prevalence of pfmdr1 polymorphisms.

Artemisinin-based combination therapies (ACTs) have been vital in reducing malaria mortality rates since the 2000s. Their efficacy, however, is threatened by the emergence and spread of artemisinin resistance in Southeast Asia. The Plasmodium falciparum multidrug resistance protein 1 (PfMDR1) transporter plays a central role in parasite resistance to ACT partner drugs through gene copy number variations (CNV) and/or single nucleotide polymorphisms (SNPs). Using genomic epidemiology, we show that multiple pfmdr1 copies encoding the N86 and 184F haplotype are prevalent across Southeast Asia. Applying genome editing tools on the Southeast Asian Dd2 strain and using a surrogate assay to measure transporter activity in infected red blood cells, we demonstrate that parasites harboring multicopy N86/184F PfMDR1 have a higher Fluo-4 transport capacity compared with those expressing the wild-type N86/Y184 haplotype. Multicopy N86/184F PfMDR1 is also associated with decreased parasite susceptibility to lumefantrine. These findings provide evidence of the geographic selection and expansion of specific multicopy PfMDR1 haplotypes associated with multidrug resistance in Southeast Asia.

Kelch 13-propeller polymorphisms in Plasmodium falciparum from Jazan region, southwest Saudi Arabia

BACKGROUND

Artemisinin-based combination therapy (ACT) is recommended at the initial phase for treatment of Plasmodium falciparum, to reduce morbidity and mortality in all countries where malaria is endemic. Polymorphism in portions of P. falciparum gene encoding kelch (K13)-propeller domains is associated with delayed parasite clearance after ACT. Of about 124 different non-synonymous mutations, 46 have been identified in Southeast Asia (SEA), 62 in sub-Saharan Africa (SSA) and 16 in both the regions. This is the first study designed to analyse the prevalence of polymorphism in the P. falciparum k13-propeller domain in the Jazan region of southwest Saudi Arabia, where malaria is endemic.

METHODS

One-hundred and forty P. falciparum samples were collected from Jazan region of southwest Saudi Arabia at three different times: 20 samples in 2011, 40 samples in 2016 and 80 samples in 2020 after the implementation of ACT. Plasmodium falciparum kelch13 (k13) gene DNA was extracted, amplified, sequenced, and analysed using a basic local alignment search tool (BLAST).

RESULTS

This study obtained 51 non-synonymous (NS) mutations in three time groups, divided as follows: 6 single nucleotide polymorphisms (SNPs) '11.8%' in samples collected in 2011 only, 3 (5.9%) in 2011and 2016, 5 (9.8%) in 2011 and 2020, 5 (9.8%) in 2016 only, 8 (15.7%) in 2016 and 2020, 14 (27.5%) in 2020 and 10 (19.6%) in all the groups. The BLAST revealed that the 2011 isolates were genetically closer to African isolates (53.3%) than Asian ones (46.7%). Interestingly, this proportion changed completely in 2020, to become closer to Asian isolates (81.6%) than to African ones (18.4%).

CONCLUSIONS

Despite the diversity of the identified mutations in the k13-propeller gene, these data did not report widespread artemisinin-resistant polymorphisms in the Jazan region where these samples were collected. Such a process would be expected to increase frequencies of mutations associated with the resistance of ACT.