Expansion of artemisinin partial resistance mutations and lack of histidine rich protein-2 and -3 deletions in Plasmodium falciparum infections from Rukara, Rwanda

BACKGROUND

Emerging artemisinin partial resistance and diagnostic resistance are a threat to malaria control in Africa. Plasmodium falciparum kelch13 (k13) propeller-domain mutations that confer artemisinin partial resistance have emerged in Africa. k13-561H was initially described at a frequency of 7.4% from Masaka in 2014-2015, but not present in nearby Rukara. By 2018, 19.6% of isolates in Masaka and 22% of isolates in Rukara contained the mutation. Longitudinal monitoring is essential to inform control efforts. In Rukara, an assessment was conducted to evaluate recent k13-561H prevalence changes, as well as other key mutations. Prevalence of hrp2/3 deletions was also assessed.

METHODS

Samples collected in Rukara in 2021 were genotyped for key artemisinin and partner drug resistance mutations using molecular inversion probe assays and for hrp2/3 deletions using qPCR.

RESULTS

Clinically validated k13 artemisinin partial resistance mutations continue to increase in prevalence with the overall level of mutant infections reaching 32% in Rwanda. The increase appears to be due to the rapid emergence of k13-675V (6.4%, 6/94 infections), previously not observed, rather than continued expansion of 561H (23.5% 20/85). Mutations to partner drugs and other anti-malarials were variable, with high levels of multidrug resistance 1 (mdr1) N86 (95.5%) associated with lumefantrine decreased susceptibility and dihydrofolate reductase (dhfr) 164L (24.7%) associated with a high level of antifolate resistance, but low levels of amodiaquine resistance polymorphisms with chloroquine resistance transporter (crt) 76T: at 6.1% prevalence. No hrp2 or hrp3 gene deletions associated with diagnostic resistance were found.

CONCLUSIONS

Increasing prevalence of artemisinin partial resistance due to k13-561H and the rapid expansion of k13-675V is concerning for the longevity of artemisinin effectiveness in the region. False negative RDT results do not appear to be an issue with no hrp2 or hpr3 deletions detected. Continued molecular surveillance in this region and surrounding areas is needed to follow artemisinin partial resistance and provide early detection of partner drug resistance, which would likely compromise control and increase malaria morbidity and mortality in East Africa.

Examining the Early Distribution of the Artemisinin-Resistant Plasmodium falciparum kelch13 R561H Mutation in Areas of Higher Transmission in Rwanda

Background

Artemisinin resistance mutations in Plasmodium falciparum kelch13 (Pfk13) have begun to emerge in Africa, with Pfk13-R561H being the first reported in Rwanda in 2014, but limited sampling left questions about its early distribution and origin.

Methods

We genotyped P. falciparum positive dried blood spot (DBS) samples from a nationally representative 2014-2015 Rwanda Demographic Health Surveys (DHS) HIV study. DBS were subsampled from DHS sampling clusters with >15% P. falciparum prevalence, as determined by rapid testing or microscopy done during the DHS study (n clusters = 67, n samples = 1873).

Results

We detected 476 parasitemias among 1873 residual blood spots from a 2014-2015 Rwanda Demographic Health Survey. We sequenced 351 samples: 341/351 were wild-type (97.03% weighted), and 4 samples (1.34% weighted) harbored R561H that were significantly spatially clustered. Other nonsynonymous mutations found were V555A (3), C532W (1), and G533A (1).

Conclusions

Our study better defines the early distribution of R561H in Rwanda. Previous studies only observed the mutation in Masaka as of 2014, but our study indicates its presence in higher-transmission regions in the southeast of the country at that time.

Molecular markers of dihydroartemisinin-piperaquine resistance in northwestern Thailand

BACKGROUND

Dihydroartemisinin-piperaquine (DHA-PPQ) combination therapy is the current first-line treatment for Plasmodium falciparum malaria in Thailand. Since its introduction in 2015, resistance to this drug combination has emerged in the eastern part of the Greater Mekong Subregion including the eastern part of Thailand near Cambodia. This study aimed to assess whether the resistance genotypes have arisen the western part of country.

METHODS

Fifty-seven P. falciparum-infected blood samples were collected in Tak province of northwestern Thailand between 2013 and 2019. Resistance to DHA was examined through the single nucleotide polymorphisms (SNPs) of kelch13. PPQ resistance was examined through the copy number plasmepsin-2 and the SNPs of Pfcrt.

RESULTS

Among the samples whose kelch13 were successfully sequenced, approximately half (31/55; 56%) had mutation associated with artemisinin resistance, including G533S (23/55; 42%), C580Y (6/55; 11%), and G538V (2/55; 4%). During the study period, G533S mutation appeared and increased from 20% (4/20) in 2014 to 100% (9/9) in 2019. No plasmepsin-2 gene amplification was observed, but one sample (1/54) had the Pfcrt F145I mutation previously implicated in PPQ resistance.

CONCLUSIONS

Kelch13 mutation was common in Tak Province in 2013-2019. A new mutation G533S emerged in 2014 and rose to dominance in 2019. PPQ resistance marker Pfcrt F145I was also detected in 2019. Continued surveillance of treatment efficacy and drug resistance markers is warranted.

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.

Decreased K13 Abundance Reduces Hemoglobin Catabolism and Proteotoxic Stress, Underpinning Artemisinin Resistance

Increased tolerance of Plasmodium falciparum to front-line artemisinin antimalarials (ARTs) is associated with mutations in Kelch13 (K13), although the precise role of K13 remains unclear. Here, we show that K13 mutations result in decreased expression of this protein, while mislocalization of K13 mimics resistance-conferring mutations, pinpointing partial loss of function of K13 as the relevant molecular event. K13-GFP is associated with ∼170 nm diameter doughnut-shaped structures at the parasite periphery, consistent with the location and dimensions of cytostomes. Moreover, the hemoglobin-peptide profile of ring-stage parasites is reduced when K13 is mislocalized. We developed a pulse-SILAC approach to quantify protein turnover and observe less disruption to protein turnover following ART exposure when K13 is mislocalized. Our findings suggest that K13 regulates digestive vacuole biogenesis and the uptake/degradation of hemoglobin and that ART resistance is mediated by a decrease in heme-dependent drug activation, less proteotoxicity, and increased survival of parasite ring stages.

Identification of Mutations in Antimalarial Resistance Gene Kelch13 from plasmodium falciparum Isolates in Kano, Nigeria

Malaria control relies on first-line treatments that use artemisinin-combination therapies (ACT). Unfortunately, mutations in the plasmodium falciparum kelch13 gene result in delayed parasite clearance. Research on what is causing ACT failure is non-existent in northwestern Nigeria. Thus, the presence of mutations in kelch13 in P. falciparum isolates from Kano, Nigeria was investigated in this study. Microscopic examination of 154 blood samples obtained from patients revealed a high prevalence of P. falciparum infection (114 positive individuals, slide positivity rate = 74.03%). The 114 patients were administered Cartef® (ACT) and out of the 50 patients that returned for the 14-day follow up, 11 were positive for P. falciparum (slide positivity rate = 22%). On day 0, 80 samples out of 114 and 11 samples on day 14 (91 out of 125 microscopy-positive samples) were positive with Plasmodium according to the PCR of cytochrome oxidase I, which corresponds to 72.8%. A fragment of the kelch13 gene encompassing the propeller domains was sequenced in 49 samples, alongside samples of the susceptible strain pf_3D7. Low polymorphism was observed, suggesting a lack of selection on this gene, and only six mutations (Glu433Gly, Phe434Ile, Phe434Ser, Ile684Asn, Ile684Thr and Glu688Lys) were found. The epidemiologic impact of these mutations and their potential role in ACT resistance needs to be investigated further.

The extended recovery ring-stage survival assay provides a superior association with patient clearance half-life and increases throughput

BACKGROUND

Tracking and understanding artemisinin resistance is key for preventing global setbacks in malaria eradication efforts. The ring-stage survival assay (RSA) is the current gold standard for in vitro artemisinin resistance phenotyping. However, the RSA has several drawbacks: it is relatively low throughput, has high variance due to microscopy readout, and correlates poorly with the current benchmark for in vivo resistance, patient clearance half-life post-artemisinin treatment. Here a modified RSA is presented, the extended Recovery Ring-stage Survival Assay (eRRSA), using 15 cloned patient isolates from Southeast Asia with a range of patient clearance half-lives, including parasite isolates with and without kelch13 mutations.

METHODS

Plasmodium falciparum cultures were synchronized with single layer Percoll during the schizont stage of the intraerythrocytic development cycle. Cultures were left to reinvade to early ring-stage and parasitaemia was quantified using flow cytometry. Cultures were diluted to 2% haematocrit and 0.5% parasitaemia in a 96-well plate to start the assay, allowing for increased throughput and decreased variability between biological replicates. Parasites were treated with 700 nM of dihydroartemisinin or 0.02% dimethyl sulfoxide (DMSO) for 6 h, washed three times in drug-free media, and incubated for 66 or 114 h, when samples were collected and frozen for PCR amplification. A SYBR Green-based quantitative PCR method was used to quantify the fold-change between treated and untreated samples.

RESULTS

15 cloned patient isolates from Southeast Asia with a range of patient clearance half-lives were assayed using the eRRSA. Due to the large number of pyknotic and dying parasites at 66 h post-exposure (72 h sample), parasites were grown for an additional cell cycle (114 h post-exposure, 120 h sample), which drastically improved correlation with patient clearance half-life compared to the 66 h post-exposure sample. A Spearman correlation of - 0.8393 between fold change and patient clearance half-life was identified in these 15 isolates from Southeast Asia, which is the strongest correlation reported to date.

CONCLUSIONS

eRRSA drastically increases the efficiency and accuracy of in vitro artemisinin resistance phenotyping compared to the traditional RSA, which paves the way for extensive in vitro phenotyping of hundreds of artemisinin resistant parasites.

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.

Absence of kelch13 artemisinin resistance markers but strong selection for lumefantrine-tolerance molecular markers following 18 years of artemisinin-based combination therapy use in Mpumalanga Province, South Africa (2001-2018)

Background

The ability of Plasmodium falciparum parasites to develop resistance to widely used anti-malarials threatens malaria control and elimination efforts. Regular drug efficacy monitoring is essential for ensuring effective treatment policies. In low transmission settings where therapeutic efficacy studies are often not feasible, routine surveillance for molecular markers associated with anti-malarial resistance provides an alternative for the early detection of emerging resistance. Such a longitudinal survey of changes in the prevalence of selected molecular markers of resistance was conducted in the malaria-endemic regions of Mpumalanga Province, South Africa, where malaria elimination at a district-level is being pursued.

Methods

Molecular analyses to determine the prevalence of alleles associated with resistance to lumefantrine (mdr86N, crt76K and mdr1 copy number variation) and sulfadoxine–pyrimethamine (dhfr triple, dhps double, SP quintuple) were conducted between 2001 and 2018, while artemisinin resistance markers (kelch13 mutations) were assessed only in 2018.

Results

Parasite DNA was successfully amplified from 1667/2393 (70%) of malaria-positive rapid diagnostic tests routinely collected at primary health care facilities. No artemisinin resistance-associated kelch13 mutations nor amplification of the mdr1 gene copy number associated with lumefantrine resistance were observed. However, prevalence of both the mdr86N and crt76K alleles increased markedly over the study period, with all isolates collected in 2018 carrying these markers. SP quintuple mutation prevalence increased steadily from 14% in 2001 to 96% in 2018. Mixed alleles at any of the codons assessed were rare by 2018.

Conclusion

No kelch13 mutations confirmed or suspected to be associated with artemisinin resistance were identified in 2018. Although parasites carrying the mdr86N and crt76K alleles associated with reduced lumefantrine susceptibility were strongly selected for over the study period, nearing fixation by 2018, the marker for lumefantrine resistance, namely increased mdr1 copy number, was not observed in this study. The increase in mdr86N and crt76K allele prevalence together with intense regional artemether–lumefantrine drug pressure, raises concern regarding the sustained artemether–lumefantrine efficacy. Regular, rigorous anti-malarial resistance marker surveillance across all three South African malaria-endemic provinces to inform case management is recommended.

Emergence and Spread of kelch13 Mutations Associated with Artemisinin Resistance in Plasmodium falciparum Parasites in 12 Thai Provinces from 2007 to 2016

Artemisinin-based combination therapy (ACT) is the most effective and widely used treatment for uncomplicated Plasmodium falciparum malaria and is a cornerstone for malaria control and prevention globally. Resistance to artemisinin derivatives has been confirmed in the Greater Mekong Subregion (GMS) and manifests as slow parasite clearance in patients and reduced ring stage susceptibility to artemisinins in survival assays. The P. falciparumkelch13 gene mutations associated with artemisinin-resistant parasites are now widespread in the GMS. We genotyped 277 samples collected during an observational study from 2012 to 2016 from eight provinces in Thailand to identify P. falciparum kelch13 mutations. The results were combined with previously reported genotyping results from Thailand to construct a map illustrating the evolution of P. falciparum kelch13 mutations from 2007 to 2016 in that country. Different mutant alleles were found in strains with different geographical origins. The artemisinin resistance-conferring Y493H and R539T mutations were detected mainly in eastern Thailand (bordering Cambodia), while P574L was found only in western Thailand and R561H only in northwestern Thailand. The C580Y mutation was found across the entire country and was nearing fixation along the Thai-Cambodia border. Overall, the prevalence of artemisinin resistance mutations increased over the last 10 years across Thailand, especially along the Thai-Cambodia border. Molecular surveillance and therapeutic efficacy monitoring should be intensified in the region to further assess the extent and spread of artemisinin resistance.

Molecular Epidemiology of Plasmodium falciparum kelch13 Mutations in Senegal Determined by Using Targeted Amplicon Deep Sequencing

The emergence of Plasmodium falciparum resistance to artemisinin in Southeast Asia threatens malaria control and elimination activities worldwide. Multiple polymorphisms in the P. falciparum kelch gene found in chromosome 13 (Pfk13) have been associated with artemisinin resistance. Surveillance of potential drug resistance loci within a population that may emerge under increasing drug pressure is an important public health activity. In this context, P. falciparum infections from an observational surveillance study in Senegal were genotyped using targeted amplicon deep sequencing (TADS) for Pfk13 polymorphisms. The results were compared to previously reported Pfk13 polymorphisms from around the world. A total of 22 Pfk13 propeller domain polymorphisms were identified in this study, of which 12 have previously not been reported. Interestingly, of the 10 polymorphisms identified in the present study that were also previously reported, all had a different amino acid substitution at these codon positions. Most of the polymorphisms were present at low frequencies and were confined to single isolates, suggesting they are likely transient polymorphisms that are part of naturally evolving parasite populations. The results of this study underscore the need to identify potential drug resistance loci existing within a population, which may emerge under increasing drug pressure.

Genomic epidemiology of artemisinin resistant malaria

The current epidemic of artemisinin resistant Plasmodium falciparum in Southeast Asia is the result of a soft selective sweep involving at least 20 independent kelch13 mutations. In a large global survey, we find that kelch13 mutations which cause resistance in Southeast Asia are present at low frequency in Africa. We show that African kelch13 mutations have originated locally, and that kelch13 shows a normal variation pattern relative to other genes in Africa, whereas in Southeast Asia there is a great excess of non-synonymous mutations, many of which cause radical amino-acid changes. Thus, kelch13 is not currently undergoing strong selection in Africa, despite a deep reservoir of variations that could potentially allow resistance to emerge rapidly. The practical implications are that public health surveillance for artemisinin resistance should not rely on kelch13 data alone, and interventions to prevent resistance must account for local evolutionary conditions, shown by genomic epidemiology to differ greatly between geographical regions.

DOI:http://dx.doi.org/10.7554/eLife.08714.001

Malaria is an infectious disease caused by a microscopic parasite called Plasmodium, which is transferred between humans by mosquitos. One species of malaria parasite called Plasmodium falciparum can cause particularly severe and life-threatening forms of the disease. Currently, the most widely used treatment for P. falciparum infections is artemisinin combination therapy, a treatment that combines the drug artemisinin (or a closely related molecule) with another antimalarial drug. However, resistance to artemisinin has started to spread throughout Southeast Asia.

Artemisinin resistance is caused by mutations in a parasite gene called kelch13, and researchers have identified over 20 different mutations in P. falciparum that confer artemisinin resistance. The diversity of mutations involved, and the fact that the same mutation can arise independently in different locations, make it difficult to track the spread of resistance using conventional molecular marker approaches.

Here, Amato, Miotto et al. sequenced the entire genomes of more than 3,000 clinical samples of P. falciparum from Southeast Asia and Africa, collected as part of a global network of research groups called the MalariaGEN Plasmodium falciparum Community Project. Amato, Miotto et al. found that African parasites had independently acquired many of the same kelch13 mutations that are known to cause resistance to artemisinin in Southeast Asia. However the kelch13 mutations seen in Africa remained at low levels in the parasite population, and appeared to be under much less pressure for evolutionary selection than those found in Southeast Asia.

These findings demonstrate that the emergence and spread of resistance to antimalarial drugs does not depend solely on the mutational process, but also on other factors that influence whether the mutations will spread in the population. Understanding how this is affected by different patterns of drug treatments and other environmental conditions will be important in developing more effective strategies for combating malaria.

DOI:http://dx.doi.org/10.7554/eLife.08714.002

Plasmodium falciparum kelch 13: a potential molecular marker for tackling artemisinin-resistant malaria parasites

Although artemisinin combination therapies have been deployed as a first-line treatment for uncomplicated malaria in almost all endemic countries, artemisinin-resistant parasites have emerged and have gradually spread across the Greater Mekong subregions. There is growing concern that the resistant parasites may migrate to or emerge indigenously in sub-Saharan Africa, which might provoke a global increase in malaria-associated morbidity and mortality. Therefore, development of molecular markers that enable identification of artemisinin resistance with high sensitivity is urgently required to combat this issue. In 2014, a potential artemisinin-resistance responsible gene, Plasmodium falciparum kelch13, was discovered. Here, we review the genetic features of P. falciparum kelch13 and discuss its related resistant mechanisms and potential as a molecular marker.