What exactly does the PfK13 C580Y mutation in Plasmodium falciparum influence?
Parasites & Vectors volume 16, Article number: 421 (2023)
The emergence and spread of artemisinin resistance threaten global malaria control and elimination goals, and encourage research on the mechanisms of drug resistance in malaria parasites. Mutations in Plasmodium falciparum Kelch 13 (PfK13) protein are associated with artemisinin resistance, but the unique or common mechanism which results in this resistance is unclear.
We analyzed the effects of the PfK13 mutation on the transcriptome and proteome of P. falciparum at different developmental stages. Additionally, the number of merozoites, hemozoin amount, and growth of P. falciparum 3D7C580Y and P. falciparum 3D7WT were compared. The impact of iron supplementation on the number of merozoites of P. falciparum 3D7C580Y was also examined.
We found that the PfK13 mutation did not significantly change glycolysis, TCA, pentose phosphate pathway, or oxidative phosphorylation, but did reduce the expression of reproduction- and DNA synthesis-related genes. The reduced number of merozoites, decreased level of hemozoin, and slowed growth of P. falciparum 3D7C580Y were consistent with these changes. Furthermore, adding iron supply could increase the number of the merozoites of P. falciparum 3D7C580Y.
These results revealed that the PfK13 mutation reduced hemoglobin ingestion, leading to artemisinin resistance, likely by decreasing the parasites' requirement for haem and iron. This study helps elucidate the mechanism of artemisinin resistance due to PfK13 mutations.
The emergence and spread of artemisinin resistance endanger worldwide malaria control and elimination goals [1,2,3,4]. Although the mechanism of artemisinin resistance remains a subject of debate, numerous mutations in the Plasmodium falciparum Kelch 13 (PfK13) protein have been associated with artemisinin resistance [5,6,7,8]. Frequently occurring mutations of PfK13 include C580Y, R539T, I543, and Y493Y, with the C580Y mutation becoming increasingly dominant in the Greater Mekong subregion . PfK13 is located in the cytostome (cell mouth), through which the parasite takes up hemoglobin , the uptake and digestion of which likely determine artemisinin activity . PfK13 mutation likely confers artemisinin resistance by dampening hemoglobin endocytosis and artemisinin activation [8, 10]. Although malaria artemisinin resistance is thought to be associated with different PfK13 mutations and relevant genetic and physiological backgrounds [12, 13], not all artemisinin resistance is due to the PfK13 mutation [13,14,15,16,17,18]. Some parasites lacking PfK13 mutations were still parasitaemic on day 3 after treatment with artemisinin derivatives , thus, the PfK13 mutations that reduced hemoglobin endocytosis and led to decreased artemisinin activation were not the sole causes of artemisinin resistance. Malaria parasites that survive with PfK13 mutations show reduced hemoglobin endocytosis and have lower haem and iron (HI) requirements. We recently showed that malaria parasites depend on HI, and that artemisinin likely kills malaria parasites by disturbing their utilization of HI . The high HI requirement of Plasmodium parasites may result in their high sensitivity to artemisinin [19, 20]; thus, there may be a different explanation for artemisinin resistance in Plasmodium with PfK13 mutations. However, regardless of which explanation is better supported, it is essential to clarify how PfK13 mutations affect malaria parasites.
In vitro culture of P. falciparum parasites and synchronization
Plasmodium falciparum 3D7WT and 3D7C580Y were cultured in Roswell Park Memorial Institute (RPMI) 1640 medium with 25 mmol/L HEPES, 0.5% AlbuMAX, 0.2% sodium bicarbonate, 0.2 mmol/L htpoxanthine, and 20 μg/mL gentamicin sulfate at 37 °C in an incubator supplied with 5% CO2 and 5% O2. Fresh O positive human red blood cells from healthy human donors were used to maintain parasitaemia; the parasites were maintained at 2% hematocrit, based on a previous report . The parasites were synchronized by treatment with 5% sorbitol, incubated at 37 °C for 10 min, washed with RPMI 1640, and then returned to normal culture conditions. Mature schizonts were purified using a 40–70% Percoll gradient, washed with RPMI 1640, and cultured for 3 or 6 h with fresh erythrocytes. The cultures were further treated with 5% sorbitol to obtain 0- to 3 or 0- to 6-h ring stages. Parasitaemia was calculated from Giemsa blood smears observed under a light microscope.
Sample preparation and RNA isolation
The cultures were tightly synchronized to obtain 0- to 3-h ring stages. Parasites were harvested at 0, 6, 12, 18, 24, 30, 36, and 42 h post-infection (hpi) for transcriptome analysis. Samples at different stages were collected by centrifugation at 1000 g for 3 min and preserved in TRIzol reagent (Invitrogen, USA). RNA was extracted using phenol and isopropanol precipitation. RNA quantification and quality determination were performed using Nanodrop 2000 (Thermo Fisher Scientific, Waltham, MA).
Library preparation and RNA sequencing
RNA integrity was assessed using an Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA). Libraries were constructed using the TruSeq Stranded mRNA LT Sample Prep Kit (Illumina, San Diego, CA) according to the manufacturer’s instructions. Transcriptome sequencing and analysis were conducted by LC Sciences (Hangzhou, Zhejiang, China). The libraries were sequenced on an Illumina HiSeq X Ten platform, and 150 base pair paired-end reads were generated. About 49.10 million raw reads for each sample were generated. Raw data in FASTQ format were first processed using Trimmomatic, and low-quality reads were removed to obtain the clean reads . About 48.02 million clean reads for each sample were retained for subsequent analyses.
Differentially expressed gene analysis
Clean reads were mapped to the P. falciparum genome by using HISAT2 . The fragments per kilobase of transcript per million mapped reads of each gene was calculated using Cufflinks [24, 25], and the read counts of each gene were obtained using htseq-count . Differential expression analysis was performed using the DESeq (2012) package in R (https://bioconductor.org/packages/release/bioc/html/edgeR.html). The threshold for significant differential expression was set at P- < 0.05, fold change > 2, or fold change < 0.5, and the minimum number of reads for at least one sample had to exceed 100. Gene Ontology (GO) enrichment analysis was performed using R or David software (https://david.ncifcrf.gov). Differentially expressed genes were functionally annotated using PlasmoDB and UniProt.
Sample preparation and mass spectrometry analysis
Plasmodium falciparum 3D7WT and P. falciparum 3D7C580Y were treated to obtain parasites 0–3 h after invasion, as described above, and samples were harvested at 12, 24, and 30 h. Extraction, purification and mass spectrometry analysis of the protein were performed by LC-Bio Technology (Hangzhou, China) through liquid chromatography–tandem mass spectrometry. During the differentially expressed protein analysis, fold change > 2 or fold change < 0.5 served as the screening criterion. The selected differentially expressed proteins were subjected to GO enrichment analysis.
Primer design and real-time polymerase chain reaction analysis
Total RNA from P. falciparum 3D7WT and P. falciparum 3D7C580Y was isolated using TRIzol reagent. Following DNase treatment, total RNA (800 ng) was reverse transcribed into complementary DNA using the ChamQ Universal SYBR qPCR master mix (Vazyme, Nanjing, China) according to the manufacturer’s instructions on a 7500 Real-Time PCR System (Applied Biosystems). We selected genes such as merozoite structure-related genes, merozoite surface protein MSA180 gene, rhoptry protein genes, DNA replication-related genes, DNA replication licensing factor MCM5 gene, and proliferating cell nuclear antigen 1 gene, which showed differential expression in the transcriptomes. Seryl-transfer RNA synthetase served as the internal control, and its expression was measured by detecting the fluorescence from SYBR green dye incorporated into the amplifying target. RNA expression was normalized to the corresponding internal control genes, and relative changes were calculated using the 2−\(\Delta \Delta\)Ct method [27, 28]. The experiment was performed using three independent biological and technical replicates. For primer design, we used Primer BLAST (https://www.ncbi.nlm.nih.gov/tools/primer-blast/). Primers used for quantitative real-time polymerase chain reacion (qPCR) are listed in Additional file 2: Table S1.
Comparison of growth curves
Plasmodium falciparum 3D7WT and P. falciparum 3D7C580Y were synchronized several times to obtain 0- to 3-h ring stages. The initial infection rate was maintained at 0.3% with hematocrit at 2%, after which normal culture conditions were used. The parasites were cultured for 144 h, and infected red blood cells were collected every 48 h. These were then incubated with Hoechst 33342 (catalogue no. C1028; Biyuntian) in the dark at 37 °C for 10 min. After two washes with RPMI 1640, positive cells were sorted using the BD FACS AriaII cell sorter based on a 355-nm (Hoechst 33342) fluorescence signal. Uninfected red blood cells were used as controls. All data were processed using FlowJo software (10.8.1).
Analysis of the number of merozoites
Mature schizonts were synchronized using a 40–70% Percoll gradient and cultured at 37 °C for 3 h with shaking to reduce multiple infections. After synchronization with 5% sorbitol, the parasite was cultured for another 42–45 h. The number of merozoites in the schizonts was counted under light microscopy after Giemsa staining. Data were collected through repeated experiments conducted by three operators.
Comparison of hemozoin content
Plasmodium falciparum 3D7WT and P. falciparum 3D7C580Y were purified using a 40–70% Percoll gradient and treated with 5% sorbitol to obtain 6-h ring stages. Then, P. falciparum 3D7WT and P. falciparum 3D7C580Y were cultured for 36–42 h with parasitaemia adjusted to 0.5–1% by using normal erythrocytes while maintaining 2% hematocrit. Approximately 10 ml of the culture medium was centrifuged at 1000 g for 3 min, and then treated with 0.15% saponin solution on ice for 10 min to lyse the cells. After centrifugation at 12,000 g for 10 min, the supernatant was discarded, and the sediment was washed with 25 mmol/L Tris (pH 7.8) containing 2.5% sodium dodecyl sulphate until the supernatant was clear. The sediment was dissolved in 250 µL of 2.5% sodium dodecyl sulphate buffer and 20 µL of 2.5 mol/L NaOH. The absorbance of hemozoin at 400 nm was measured using a Nanodrop 2000 spectrophotometer to quantify its content .
In vitro iron sucrose assay
Plasmodium falciparum 3D7WT and P. falciparum 3D7C580Y were tightly synchronized to obtain the 0- to 3-h ring stages. Synchronized ring-stage parasites were incubated with 50 µmol/L iron sucrose, and the control group was incubated with an equal volume of RPMI 1640 complete medium. After 90–96 h of incubation, the merozoites were counted under a light microscope following Giemsa staining. Data were obtained from experiments repeated at least three times.
Plasmodium falciparum 3D7WT at 1.5–2.5% parasitaemia was treated with desferrioxamine (DFO) (MedChemExpress, HY-B1625) at a final concentration of 500 nmol/L for 24 h. A portion of the mixture was centrifuged at 1000 g for 5 min to remove the supernatant. After the addition of cells suspended in PBS and Hoechst 33342 (1:100) (Biyuntian), the mixture was incubated at 37 °C for 10 min under continuous mixing, centrifuged at 1000 g for 5 min, and washed three times with phosphate-buffered saline. The positive cell rate was determined using a BD FACS AriaII cell sorter. Infected red blood cells were identified by using 355-nm fluorescence spectroscopy, with uninfected cells serving as controls. All data were processed using FlowJo software (10.8.1).
Graphpad Prism 8.0.2 software was used to perform the statistical analysis. Student's t-test (unpaired) was performed to compare the changes in parasitaemia in the groups. P < 0.05 was considered to indicate statistical significance.
Different effects of PfK13 mutation on the transcriptome and proteome of P. falciparum at different stages
The transcriptomes of P. falciparum 3D7WT and P. falciparum 3D7C580Y were analyzed at 0, 6, 12, 18, 24, 30, 36, and 42 hpi (Fig. 1A–C; Additional file 3: Table S2). Among the 5282 identified genes, 123, 128, 111, 130, 91, 138, 253, and 164 were differentially expressed at various time points post-infection (Additional file 3: Table S2). At 36 and 42 hpi, the downregulated genes (P. falciparum 3D7C580Y vs. P. falciparum 3D7WT) accounted for 96.05% and 95.12% of the differentially expressed genes, respectively (Fig. 1B). Additionally, proteomes of the samples at 12, 24, and 30 hpi were also analyzed. Among all the identified 2558 proteins, 366, 279, and 648 differentially expressed proteins were detected (Additional file 4: Table S3). The downregulated proteins (P. falciparum 3D7C580Y vs. P. falciparum 3D7WT) accounted for 61.20%, 69.18%, and 56.94%, respectively (Figs. 1D, E; Additional file 4: Table S3).
PfK13 mutation especially decreased the expression of DNA synthesis- and merozoite structure-related genes
The characteristic differentially expressed genes between P. falciparum 3D7C580Y and P. falciparum 3D7WT at various time points post-infection indicated differences or changes in main functions at these various time points. Especially at 18 hpi, DNA synthesis-related genes and proteins, such as the minichromosome maintenance protein complex, proliferating cell nuclear antigen complex, and DNA polymerase processivity factor complex were downregulated. At 30 h and 36 hpi, all or most merozoite structure-related genes, such as rhoptry, merozoite dense granules, microneme, and apical complex genes were downregulated (Fig. 2A–C). Proteome analysis also showed that at 24 and 30 hpi, the expression of merozoite structure-related proteins was downregulated (Fig. 2D–F), which supported the results of the transcriptome analysis.
PfK13 mutation did not significantly affect metabolism but did weaken reproduction
To investigate the effect of PfK13 protein mutation on various metabolic processes in Plasmodium, we analyzed the expression of genes associated with glycolysis, the tricarboxylic acid cycle, the pentose phosphate pathway (PPP), oxidative phosphorylation, reproduction, DNA replication, pyrimidine metabolism, fatty acid metabolism, amino acid metabolism, purine metabolism, and redox pathways, including the glutathione and thioredoxin systems (Fig. 3; Additional file 1: Fig. S1), which have often been associated with drug resistance in malaria parasites . The expression of most genes associated with glycolysis, the tricarboxylic acid cycle, the PPP, oxidative phosphorylation, fatty acid metabolism, amino acid metabolism, purine and pyrimidine metabolism, and redox processes was unaffected, that is, most genes were not differentially expressed between P. falciparum 3D7C580Y and P. falciparum 3D7WT. Although a few genes were differentially expressed, some of them were upregulated or downregulated for P. falciparum 3D7C580Y and P. falciparum 3D7WT, and it was difficult to determine if the overall process was being facilitated or inhibited (Fig. 3; Additional file 1: Fig. S1). However, most of the genes associated with merozoites, sporozoites, and gametocytes were downregulated at 30 and 36 hpi (Fig. 3E), and most of the genes associated with DNA replication were downregulated at 18 hpi (Fig. 3G), which is consistent with the findings of a previous study . Some of these results were further validated by real-time PCR, which revealed the close functional relationships of the genes (Fig. 3F, H).
P. falciparum 3D7C580Y had fewer merozoites than P. falciparum 3D7WT
An inhibitory effect of PfK13 mutation on reproduction was revealed. To further confirm this finding, we compared the number of merozoites between P. falciparum 3D7C580Y and P. falciparum 3D7WT. We quantified the number of merozoites in 50 mature schizonts and observed that P. falciparum 3D7C580Y had an average of 16 ± 3 merozoites, whereas P. falciparum 3D7WT had an average of 19 ± 3 (P < 0.0001) (Fig. 4A, B).
Hemozoin content of P. falciparum 3D7C580Y was lower than that of P. falciparum 3D7WT
Hemozoin, which is considered a vector of HI in hemozoin-forming organisms, functions in reproduction [29, 32]. Hemozoin content was found to correlate with the utilization of HI and artemisinin sensitivity . To investigate the effect of the PfK13 mutation on HI utilization, we compared the hemozoin content of an equivalent number of P. falciparum 3D7C580Y and P. falciparum 3D7WT at 36 hpi. The hemozoin content of P. falciparum 3D7WT was 16.02 ± 2.44 nmol and that of P. falciparum 3D7C580Y 10.72 ± 2.13 nmol. The hemozoin content of P. falciparum 3D7C580Y was significantly lower than that of P. falciparum 3D7WT (P < 0.01) (Fig. 4E).
Growth of P. falciparum 3D7C580Y was slower than that of P. falciparum 3D7WT
After three cycles of culture, the mean parasitaemia of P. falciparum 3D7C580Y was 7.35%, which was significantly lower than that observed in P. falciparum 3D7WT (11.2%) (P < 0.01). This suggested that the growth and reproductive capability of P. falciparum 3D7C580Y was lower than that of P. falciparum 3D7WT (Fig. 4F).
Addition of iron increased the number of P. falciparum 3D7C580Y merozoites
The previous findings indicate that P. falciparum 3D7C580Y had a reduced reproductive capacity and hemozoin content compared with P. falciparum 3D7WT. To investigate the effect of iron on the reproduction of P. falciparum 3D7C580Y, we supplemented the culture system with sucrose iron and quantified the number of merozoites after 90–96 h of cultivation. The number of merozoites of P. falciparum 3D7C580Y increased significantly from 16 ± 3 to 19 ± 3 (P < 0.0001) (Fig. 4C, D), which was similar to the increase seen in P. falciparum 3D7WT (Fig. 4D).
Effect of DFO on P. falciparum
Iron supplementation promoted the reproduction of P. falciparum 3D7C580Y. To further elucidate the role of iron in Plasmodium proliferation and development, we treated P. falciparum 3D7 with DFO for 24 h and observed its effect on the parasites. The parasitaemia of P. falciparum 3D7 decreased from 3.21% ± 0.55% to 1.03% ± 0.56% after treatment with DFO (P < 0.01) (Fig. 4G).
The World Health Organization recommends that the presence of persistent parasitaemia on day 3 after treatment with artemisinin should be considered confirmatory for the resistance of P. falciparum to this drug. The phenotypic manifestation of artemisinin-resistant P. falciparum infections is slowed parasite clearance after treatment with artesunate in vivo. This slow-clearance phenotype, defined by a parasite clearance half-life of > 5 h, is attributed to the loss of sensitivity of P. falciparum to artemisinin during the early stage of the intraerythrocytic developmental cycle, the ring stage [33, 34]. Resistance (or tolerance) in cultured parasites is defined as > 1% survival in the ring-stage survival assay, in which parasites at the young ring stage (0–3 hpi) are exposed to 700 nM DHA for 4–6 h, and their survival rate is measured 3 days later . This level of survival is considered to indicate partial resistance to artemisinin or its derivatives. Artemisinin resistance is considered to be associated with the genetic and transcriptomic background or initial transcriptional responses of susceptible parasites to artemisinin [13, 35]. The mechanisms by which PfK13 mutations may lead to resistance are both controversial and interesting. Previous studies have proposed that PfK13 mutations could change various aspects of the parasite’s intraerythrocytic developmental program, by influencing cell-cycle periodicity, unfolded-protein response, protein degradation, vesicular trafficking, DNA replication and repair [8, 30, 31]. Moreover, mitochondrial metabolism or processes including damage sensing and antioxidant properties likely augment the protection of P. falciparum with a PfK13 mutation from artemisinin by helping these parasites to undergo temporary quiescence and then recovery through accelerated growth after drug elimination . PfK13 mutation did not significantly influence the main metabolic processes of glycolysis, the TCA cycle, PPP, fatty acid metabolism, pyrimidine metabolism, purine metabolism, antioxidation, or the electron transport chain in the present study; however, DNA replication decreased at 18 hpi. PfK13 mutation also inhibited the expression of merozoite-related genes, which was accompanied by a reduced number of merozoites, slow growth, and decreased hemozoin content in P. falciparum 3D7C580Y (Fig. 5A, B). Why the PfK13 mutation affected mostly the expression of DNA replication- and merozoite-related genes is of interest.
Superficially, the PfK13 mutation reduces hemoglobin endocytosis, resulting in decreased amino acid supply and nutrient deficiency, ultimately leading to reduced reproductive ability. Simultaneously, reduced HI release diminishes artemisinin activation, thereby contributing to the development of artemisinin resistance . However, it is worth noting that the PfK13 mutation, which resulted in a reduction of HI release, not only diminished the activation of artemisinin but also decreased the requirement for and utilization of HI (Fig. 5C, D). When insufficient HI are provided to parasites, their functions associated with HI utilization are impaired. When the parasites do not die, their related functions and metabolism will change in accordance with the low level of HI utilized. As the level of utilized HI is low, efficient activation of the pentose phosphate pathway (PPP), which synthesizes pentose phosphate and NADPH and can thus promote DNA synthesis and GSH cycling, is hindered. Consequently, DNA synthesis declines and it becomes difficult to maintain the GSH cycle, which ultimately leads to decreased fertility . Notably, parasites that utilize low levels of HI can achieve greater resistance to HI depletion. A potential mechanism to explain this, that artemisinin kills malaria parasites by disrupting their utilization of HI, is another likely reason why the parasites develop resistance to artemisinin . This theory suggests that HI use decreases with a decreased HI requirement, leading to artemisinin-induced parasite resistance to the exhaustion of HI.
HI are crucial for the reproduction and development of malaria parasites [20, 36,37,38,39]. DNA synthesis-related enzymes in particular require an iron–sulfur cluster to form active complexes [40,41,42,43,44,45]. In our experiment, iron supplementation alone restored the number of merozoites of P. falciparum 3D7C580Y to the normal level, which indicated that the PfK13 mutation likely weakens reproductive ability by reducing hemoglobin endocytosis to decrease the HI supply. In fact, PfK13 mutation, on the one hand, can lead to a reduction in the supply of amino acids and HI, and thus delay development and reduce the activation of artemisinin. On the other hand, the parasites can become resistant to artemisinin as a consequence of the PfK13 mutation leading to a reduction in their requirement for and dependence on HI . Determining which is the dominant mechanism, or whether both mechanisms operate concurrently, is of interest.
Given that the parasites were alive after the PfK13 mutation arose, the HI requirement of P. falciparum 3D7C580Y was reduced compared with that of P. falciparum 3D7WT. Our findings revealed the differential expression of some genes associated with glycolysis, the TCA cycle, PPP, fatty acid metabolism, antioxidation, and the electron transport chain between P. falciparum 3D7C580Y and P. falciparum 3D7WT. However, regular expression patterns were not observed, unlike in DNA synthesis-related or reproduction-related genes that were downregulated in P. falciparum 3D7C580Y, which clearly indicated that PfK13 mutation primarily affected HI supply. Krugliak et al.  proposed that the parasite digests up to 65% of the host cell's hemoglobin but utilizes only up to approximately 16% of the amino acids derived from it. This large discrepancy indicates that the parasites prefer HI to amino acids present in the host hemoglobin. The PfK13 mutation may primarily lead to a decrease in the uptake of iron rather than that of amino acids in the parasites.
In the present study, the number of merozoites of P. falciparum 3D7C580Y was restored to a level similar to that of P. falciparum 3D7WT simply by iron supplementation. Furthermore, when iron was removed by using DFO, P. falciparum 3D7WT was inhibited or killed. These results confirm that iron is crucial to the parasite. Therefore, the question arises as to whether PfK13 mutation confers resistance by reducing the activation of artemisinin or by decreasing HI requirements.
In brief, although reduced activation of artemisinin due to a reduction in hemoglobin digestion and a reduced level of HI can logically explain resistance to artemisinin, it is a personalized mechanism that arises as a consequence of PfK13 mutation. It is difficult to base an explanation for artemisinin resistance on other mutations. However, a different mechanism through which the HI requirement of P. falciparum 3D7C580Y is decreased, which results in artemisinin resistance, is useful for elucidating other mutations that are responsible for resistance to artemisinin or identifying organisms that are insensitive to this drug.
Availability of data and materials
All data generated or analyzed during this study are included in the article.
- PfK13 mutation:
Plasmodium falciparum Kelch 13 protein mutation
- TCA cycle:
Tricarboxylic acid cycle
Pentose phosphate pathway
Haem and iron
- MCM proteins:
Mini-chromosome maintenance proteins
Proliferating cell nuclear antigen
- iRBCs :
Infected red blood cells
World health organization
Intraerythrocytic developmental cycle
Nicotinamide adenine dinucleotide phosphate
Lu F, Culleton R, Zhang M, Ramaprasad A, von Seidlein L, Zhou H, et al. Emergence of indigenous artemisinin-resistant Plasmodium falciparum in Africa. N Engl J Med. 2017;376:991–3.
Das S, Saha B, Hati AK, Roy S. Evidence of artemisinin-resistant Plasmodium falciparum malaria in eastern India. N Engl J Med. 2018;379:1962–4.
Ashley EA, Dhorda M, Fairhurst RM, Amaratunga C, Lim P, Suon S, et al. Spread of artemisinin resistance in Plasmodium falciparum malaria. N Engl J Med. 2014;371:411–23. https://doi.org/10.1056/NEJMoa1314981.
Dhorda M, Amaratunga C, Dondorp AM. Artemisinin and multidrug-resistant Plasmodium falciparum—a threat for malaria control and elimination. Curr Opin Infect Dis. 2021;34:432–9. https://doi.org/10.1097/QCO.0000000000000766.
Ariey F, Witkowski B, Amaratunga C, Beghain J, Langlois AC, Khim N, et al. A molecular marker of artemisinin-resistant Plasmodium falciparum malaria. Nature. 2014;505:50–5. https://doi.org/10.1038/nature12876.
Menard D, Khim N, Beghain J, Adegnika AA, Shafiul-Alam M, Amodu O, et al. A worldwide map of Plasmodium falciparum K13-propeller polymorphisms. N Engl J Med. 2016;374:2453–64. https://doi.org/10.1056/NEJMoa1513137.
Straimer J, Gnadig NF, Witkowski B, Amaratunga C, Duru V, Ramadani AP, et al. Drug resistance. K13-propeller mutations confer artemisinin resistance in Plasmodium falciparum clinical isolates. Science. 2015;347:428–31. https://doi.org/10.1126/science.1260867.
Birnbaum J, Scharf S, Schmidt S, Jonscher E, Hoeijmakers WAM, Flemming S, et al. A Kelch13-defined endocytosis pathway mediates artemisinin resistance in malaria parasites. Science. 2020;367:51–9. https://doi.org/10.1126/science.aax4735.
Yoshida N, Yamauchi M, Morikawa R, Hombhanje F, Mita T. Increase in the proportion of Plasmodium falciparum with kelch13 C580Y mutation and decline in pfcrt and pfmdr1 mutant alleles in Papua New Guinea. Malar J. 2021;20:410. https://doi.org/10.1186/s12936-021-03933-6.
Xie SC, Ralph SA, Tilley L. K13, the cytostome, and artemisinin resistance. Trends Parasitol. 2020;36:533–44. https://doi.org/10.1016/j.pt.2020.03.006.
Klonis N, Crespo-Ortiz MP, Bottova I, Abu-Bakar N, Kenny S, Rosenthal PJ, et al. Artemisinin activity against Plasmodium falciparum requires hemoglobin uptake and digestion. Proc Natl Acad Sci U S A. 2011;108:11405–10.
Zaw MT, Emran NA, Lin Z. Updates on k13 mutant alleles for artemisinin resistance in Plasmodium falciparum. J Microbiol Immunol Infect. 2018;51:159–65. https://doi.org/10.1016/j.jmii.2017.06.009.
Zhu L, Tripathi J, Rocamora FM, Miotto O, van der Pluijm R, Voss TS, et al. The origins of malaria artemisinin resistance defined by a genetic and transcriptomic background. Nat Commun. 2018;9:5158. https://doi.org/10.1038/s41467-018-07588-x.
Demas AR, Sharma AI, Wong W, Early AM, Redmond S, Bopp S, et al. Mutations in Plasmodium falciparum actin-binding protein coronin confer reduced artemisinin susceptibility. Proc Natl Acad Sci U S A. 2018;115:12799–804. https://doi.org/10.1073/pnas.1812317115.
Delandre O, Daffe SM, Gendrot M, Diallo MN, Madamet M, Kounta MB, et al. Absence of association between polymorphisms in the pfcoronin and pfk13 genes and the presence of Plasmodium falciparum parasites after treatment with artemisinin derivatives in Senegal. Int J Antimicrob Agents. 2020;56:106190. https://doi.org/10.1016/j.ijantimicag.2020.106190.
Voumbo-Matoumona DF, Akiana J, Madamet M, Kouna LC, Lekana-Douki JB, Pradines B. High prevalence of Plasmodium falciparum antimalarial drug resistance markers in isolates from asymptomatic patients from the Republic of the Congo between 2010 and 2015. J Glob Antimicrob Resist. 2018;14:277–83. https://doi.org/10.1016/j.jgar.2018.08.003.
Rocamora F, Zhu L, Liong KY, Dondorp A, Miotto O, Mok S, et al. Oxidative stress and protein damage responses mediate artemisinin resistance in malaria parasites. PLoS Pathog. 2018;14:e1006930. https://doi.org/10.1371/journal.ppat.1006930.
Mukherjee A, Bopp S, Magistrado P, Wong W, Daniels R, Demas A, et al. Artemisinin resistance without pfkelch13 mutations in Plasmodium falciparum isolates from Cambodia. Malar J. 2017;16:195. https://doi.org/10.1186/s12936-017-1845-5.
Sun J, Guo FJ, Liu C, Yin M, Xu WY, Li YN, et al. Double-kill mechanism of artemisinin against Plasmodium; 23 May 2022, PREPRINT (version 1) available at Research Square. https://doi.org/10.21203/rs.3.rs-1645582/v1.
Si W, Fang C, Liu C, Yin M, Xu W, Li Y, et al. Why is Babesia not killed by artemisinin like Plasmodium? Parasit Vectors. 2023;16:193. https://doi.org/10.1186/s13071-023-05783-4.
Bancells C, Deitsch KW. A molecular switch in the efficiency of translation reinitiation controls expression of var2csa, a gene implicated in pregnancy-associated malaria. Mol Microbiol. 2013;90:472–88. https://doi.org/10.1111/mmi.12379.
Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30:2114–20. https://doi.org/10.1093/bioinformatics/btu170.
Kim D, Langmead B, Salzberg SL. HISAT: a fast spliced aligner with low memory requirements. Nat Methods. 2015;12:357–60. https://doi.org/10.1038/nmeth.3317.
Roberts A, Trapnell C, Donaghey J, Rinn JL, Pachter L. Improving RNA-Seq expression estimates by correcting for fragment bias. Genome Biol. 2011;12:R22. https://doi.org/10.1186/gb-2011-12-3-r22.
Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, et al. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010;28:511–5. https://doi.org/10.1038/nbt.1621.
Anders S, Pyl PT, Huber W. HTSeq–a Python framework to work with high-throughput sequencing data. Bioinformatics. 2015;31:166–9. https://doi.org/10.1093/bioinformatics/btu638.
Shang X, Wang C, Fan Y, Guo G, Wang F, Zhao Y, et al. Genome-wide landscape of ApiAP2 transcription factors reveals a heterochromatin-associated regulatory network during Plasmodium falciparum blood-stage development. Nucleic Acids Res. 2022;50:3413–31. https://doi.org/10.1093/nar/gkac176.
Jiang P, Wang J, Zhu S, Hu C, Lin Y, Pan W. Identification of a Schistosoma japonicum microRNA that suppresses hepatoma cell growth and migration by targeting host FZD4 gene. Front Cell Infect Microbiol. 2022;12:786543. https://doi.org/10.3389/fcimb.2022.786543.
Sun J, Li C, Wang S. Organism-like formation of Schistosoma hemozoin and its function suggest a mechanism for anti-malarial action of artemisinin. Sci Rep. 2016;6:34463. https://doi.org/10.1038/srep34463.
Mok S, Stokes BH, Gnadig NF, Ross LS, Yeo T, Amaratunga C, et al. Artemisinin-resistant K13 mutations rewire Plasmodium falciparum’s intra-erythrocytic metabolic program to enhance survival. Nat Commun. 2021;12:530. https://doi.org/10.1038/s41467-020-20805-w.
Gibbons J, Button-Simons KA, Adapa SR, Li S, Pietsch M, Zhang M, et al. Altered expression of K13 disrupts DNA replication and repair in Plasmodium falciparum. BMC Genomics. 2018;19:849. https://doi.org/10.1186/s12864-018-5207-7.
Sun J, Hu W, Li C. Beyond heme detoxification: a role for hemozoin in iron transport in S. japonicum. Parasitol Res. 2013;112:2983–90. https://doi.org/10.1007/s00436-013-3470-8.
Chotivanich K, Tripura R, Das D, Yi P, Day NP, Pukrittayakamee S, et al. Laboratory detection of artemisinin-resistant Plasmodium falciparum. Antimicrob Agents Chemother. 2014;58:3157–61. https://doi.org/10.1128/AAC.01924-13.
Witkowski B, Khim N, Chim P, Kim S, Ke S, Kloeung N, et al. Reduced artemisinin susceptibility of Plasmodium falciparum ring stages in western Cambodia. Antimicrob Agents Chemother. 2013;57:914–23. https://doi.org/10.1128/AAC.01868-12.
Zhu L, van der Pluijm RW, Kucharski M, Nayak S, Tripathi J, White NJ, et al. Artemisinin resistance in the malaria parasite, Plasmodium falciparum, originates from its initial transcriptional response. Commun Biol. 2022;5:274. https://doi.org/10.1038/s42003-022-03215-0.
Maya-Maldonado K, Cardoso-Jaime V, Gonzalez-Olvera G, Osorio B, Recio-Totoro B, Manrique-Saide P, et al. Mosquito metallomics reveal copper and iron as critical factors for Plasmodium infection. PLoS Negl Trop Dis. 2021;15:e0009509. https://doi.org/10.1371/journal.pntd.0009509.
Goldberg DE, Sigala PA. Plasmodium heme biosynthesis: to be or not to be essential? PLoS Pathog. 2017;13:e1006511. https://doi.org/10.1371/journal.ppat.1006511.
Clark MA, Goheen MM, Fulford A, Prentice AM, Elnagheeb MA, Patel J, et al. Host iron status and iron supplementation mediate susceptibility to erythrocytic stage Plasmodium falciparum. Nat Commun. 2014;5:4446. https://doi.org/10.1038/ncomms5446.
Scholl PF, Tripathi AK, Sullivan DJ. Bioavailable iron and heme metabolism in Plasmodium falciparum. Curr Top Microbiol Immunol. 2005;295:293–324.
Netz DJ, Stith CM, Stumpfig M, Kopf G, Vogel D, Genau HM, et al. Eukaryotic DNA polymerases require an iron-sulfur cluster for the formation of active complexes. Nat Chem Biol. 2011;8:125–32. https://doi.org/10.1038/nchembio.721.
Shi R, Hou W, Wang ZQ, Xu X. Biogenesis of iron-sulfur clusters and their role in DNA metabolism. Front Cell Dev Biol. 2021;9:735678. https://doi.org/10.3389/fcell.2021.735678.
Wu Y, Brosh RM Jr. DNA helicase and helicase-nuclease enzymes with a conserved iron-sulfur cluster. Nucleic Acids Res. 2012;40:4247–60. https://doi.org/10.1093/nar/gks039.
Gari K, Leon Ortiz AM, Borel V, Flynn H, Skehel JM, Boulton SJ. MMS19 links cytoplasmic iron-sulfur cluster assembly to DNA metabolism. Science (New York, NY). 2012;337:243–5. https://doi.org/10.1126/science.1219664.
Makarova AV, Stodola JL, Burgers PM. A four-subunit DNA polymerase zeta complex containing Pol delta accessory subunits is essential for PCNA-mediated mutagenesis. Nucleic Acids Res. 2012;40:11618–26. https://doi.org/10.1093/nar/gks948.
Puig S, Ramos-Alonso L, Romero AM, Martinez-Pastor MT. The elemental role of iron in DNA synthesis and repair. Metallomics. 2017;9:1483–500. https://doi.org/10.1039/c7mt00116a.
Krugliak M, Zhang J, Ginsburg H. Intraerythrocytic Plasmodium falciparum utilizes only a fraction of the amino acids derived from the digestion of host cell cytosol for the biosynthesis of its proteins. Mol Biochem Parasitol. 2002;119:249–56.
We thank CT Fang for their valuable advice and LC Sciences (Hangzhou, Zhejiang, China) and Majorbio (Shanghai, China) for their assistance in conducting the proteomic, transcriptomic, and bioinformatics analyses.
This research was supported by the Innovation Program of Shanghai Municipal Education Commission (201901070007E00017) and National Natural Science Foundation of China (82230077).
The authors declare that they have no conflicts of interest.
Ethical approval and consent to participate
This study was carried out in strict accordance with the recommendations of the Regulations for the Administration of Affairs Concerning Experimental Animals of the State Science and Technology Commission.
Consent for publication
The protocol was approved by the Internal Review Board of Tongji University School of Medicine (TJLAC-017–039).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
. Differential expression of genes associated with pyrimidine metabolism, fatty acid metabolism, amino acid metabolism, purine metabolism, and redox processes between Plasmodium falciparum 3D7C580Y and P. falciparum 3D7WT. A Heat map of differentially expressed genes associated with pyrimidine metabolism (P. falciparum 3D7C580Y/P. falciparum 3D7WT) at different stages, B fatty acid metabolism, C amino acid metabolism, D purine metabolism, and E redox processes.
. Primer sequences of key target genes.
. Analysis of differentially expressed genes at different stages of the parasite.
. Analysis of differentially expressed proteins at 12, 24, and 30 h post-infection.
About this article
Cite this article
Si, W., Zhao, Y., Qin, X. et al. What exactly does the PfK13 C580Y mutation in Plasmodium falciparum influence?. Parasites Vectors 16, 421 (2023). https://doi.org/10.1186/s13071-023-06024-4