New Assays to Characterise Growth-Related Phenotypes of Plasmodium falciparum Reveal Variation in Density-Dependent Growth Inhibition between Parasite Lines

Beyond the MEP Pathway: A novel kinase required for prenol utilization by malaria parasites

A proposed treatment for malaria is a combination of fosmidomycin and clindamycin. Both compounds inhibit the methylerythritol 4-phosphate (MEP) pathway, the parasitic source of farnesyl and geranylgeranyl pyrophosphate (FPP and GGPP, respectively). Both FPP and GGPP are crucial for the biosynthesis of several essential metabolites such as ubiquinone and dolichol, as well as for protein prenylation. Dietary prenols, such as farnesol (FOH) and geranylgeraniol (GGOH), can rescue parasites from MEP inhibitors, suggesting the existence of a missing pathway for prenol salvage via phosphorylation. In this study, we identified a gene in the genome of P. falciparum, encoding a transmembrane prenol kinase (PolK) involved in the salvage of FOH and GGOH. The enzyme was expressed in Saccharomyces cerevisiae, and its FOH/GGOH kinase activities were experimentally validated. Furthermore, conditional knockout parasites (Δ-PolK) were created to investigate the biological importance of the FOH/GGOH salvage pathway. Δ-PolK parasites were viable but displayed increased susceptibility to fosmidomycin. Their sensitivity to MEP inhibitors could not be rescued by adding prenols. Additionally, Δ-PolK parasites lost their capability to utilize prenols for protein prenylation. Experiments using culture medium supplemented with whole/delipidated human plasma in transgenic parasites revealed that human plasma has components that can diminish the effectiveness of fosmidomycin. Mass spectrometry tests indicated that both bovine supplements used in culture and human plasma contain GGOH. These findings suggest that the FOH/GGOH salvage pathway might offer an alternate source of isoprenoids for malaria parasites when de novo biosynthesis is inhibited. This study also identifies a novel kind of enzyme related to isoprenoid metabolism.

High-content imaging as a tool to quantify and characterize malaria parasites

In 2021, Plasmodium falciparum was responsible for 619,000 reported malaria-related deaths. Resistance has been detected to every clinically used antimalarial, urging the development of novel antimalarials with uncompromised mechanisms of actions. High-content imaging allows researchers to collect and quantify numerous phenotypic properties at the single-cell level, and machine learning-based approaches enable automated classification and clustering of cell populations. By combining these technologies, we developed a method capable of robustly differentiating and quantifying P. falciparum asexual blood stages. These phenotypic properties also allow for the quantification of changes in parasite morphology. Here, we demonstrate that our analysis can be used to quantify schizont nuclei, a phenotype that previously had to be enumerated manually. By monitoring stage progression and quantifying parasite phenotypes, our method can discern stage specificity of new compounds, thus providing insight into the compound's mode of action.

Patterns of Heterochromatin Transitions Linked to Changes in the Expression of Plasmodium falciparum Clonally Variant Genes

The survival of malaria parasites in the changing human blood environment largely depends on their ability to alter gene expression by epigenetic mechanisms. The active state of Plasmodium falciparum clonally variant genes (CVGs) is associated with euchromatin characterized by the histone mark H3K9ac, whereas the silenced state is characterized by H3K9me3-based heterochromatin. Expression switches are linked to euchromatin-heterochromatin transitions, but these transitions have not been characterized for the majority of CVGs. To define the heterochromatin distribution patterns associated with the alternative transcriptional states of CVGs, we compared H3K9me3 occupancy at a genome-wide level among several parasite subclones of the same genetic background that differed in the transcriptional state of many CVGs. We found that de novo heterochromatin formation or the complete disruption of a heterochromatin domain is a relatively rare event, and for the majority of CVGs, expression switches can be explained by the expansion or retraction of heterochromatin domains. We identified different modalities of heterochromatin changes linked to transcriptional differences, but despite this complexity, heterochromatin distribution patterns generally enable the prediction of the transcriptional state of specific CVGs. We also found that in some subclones, several var genes were simultaneously in an active state. Furthermore, the heterochromatin levels in the putative regulatory region of the gdv1 antisense noncoding RNA, a regulator of sexual commitment, varied between parasite lines with different sexual conversion rates. IMPORTANCE The malaria parasite P. falciparum is responsible for more than half a million deaths every year. P. falciparum clonally variant genes (CVGs) mediate fundamental host-parasite interactions and play a key role in parasite adaptation to fluctuations in the conditions of the human host. The expression of CVGs is regulated at the epigenetic level by changes in the distribution of a type of chromatin called heterochromatin. Here, we describe at a genome-wide level the changes in the heterochromatin distribution associated with the different transcriptional states of CVGs. Our results also reveal a likely role for heterochromatin at a particular locus in determining the parasite investment in transmission to mosquitoes. Additionally, this data set will enable the prediction of the transcriptional state of CVGs from epigenomic data, which is important for the study of parasite adaptation to the conditions of the host in natural malaria infections.

A heat-shock response regulated by the PfAP2-HS transcription factor protects human malaria parasites from febrile temperatures

Periodic fever is a characteristic clinical feature of human malaria, but how parasites survive febrile episodes is not known. Although the genomes of Plasmodium species encode a full set of chaperones, they lack the conserved eukaryotic transcription factor HSF1, which activates the expression of chaperones following heat shock. Here, we show that PfAP2-HS, a transcription factor in the ApiAP2 family, regulates the protective heat-shock response in Plasmodium falciparum. PfAP2-HS activates the transcription of hsp70-1 and hsp90 at elevated temperatures. The main binding site of PfAP2-HS in the entire genome coincides with a tandem G-box DNA motif in the hsp70-1 promoter. Engineered parasites lacking PfAP2-HS have reduced heat-shock survival and severe growth defects at 37 °C but not at 35 °C. Parasites lacking PfAP2-HS also have increased sensitivity to imbalances in protein homeostasis (proteostasis) produced by artemisinin, the frontline antimalarial drug, or the proteasome inhibitor epoxomicin. We propose that PfAP2-HS contributes to the maintenance of proteostasis under basal conditions and upregulates specific chaperone-encoding genes at febrile temperatures to protect the parasite against protein damage.

Expression Patterns of Plasmodium falciparum Clonally Variant Genes at the Onset of a Blood Infection in Malaria-Naive Humans

Clonally variant genes (CVGs) play fundamental roles in the adaptation of Plasmodium falciparum to fluctuating conditions of the human host. However, their expression patterns under the natural conditions of the blood circulation have been characterized in detail for only a few specific gene families. Here, we provide a detailed characterization of the complete P. falciparum transcriptome across the full intraerythrocytic development cycle (IDC) at the onset of a blood infection in malaria-naive human volunteers. We found that the vast majority of transcriptional differences between parasites obtained from the volunteers and the parental parasite line maintained in culture occurred in CVGs. In particular, we observed a major increase in the transcript levels of most genes of the pfmc-2tm and gbp families and of specific genes of other families, such as phist, hyp10, rif, or stevor, in addition to previously reported changes in var and clag3 gene expression. Increased transcript levels of individual pfmc-2tm, rif, and stevor genes involved activation in small subsets of parasites. Large transcriptional differences correlated with changes in the distribution of heterochromatin, confirming their epigenetic nature. Furthermore, the similar expression of several CVGs between parasites collected at different time points along the blood infection suggests that the epigenetic memory for multiple CVG families is lost during transmission stages, resulting in a reset of their transcriptional state. Finally, the CVG expression patterns observed in a volunteer likely infected by a single sporozoite suggest that new epigenetic patterns are established during liver stages.

Transcriptional Analysis of Tightly Synchronized Plasmodium falciparum Intraerythrocytic Stages by RT-qPCR

In Plasmodium falciparum, the parasite responsible for the most severe forms of human malaria, many fundamental processes are controlled at the transcriptional level. Studies on diverse aspects of basic parasite biology as well as molecular epidemiology studies often rely on the ability to accurately measure transcript levels, but this is complicated by the cyclic expression patterns of the majority of malaria parasite genes. Here, we provide a complete workflow to measure transcript levels in P. falciparum intraerythrocytic blood stages, overcoming the confounding factors that are commonly encountered. The method described covers all the steps from synchronization of parasite cultures to reverse transcriptase quantitative PCR (RT-qPCR) analysis.

Artemisinin exposure at the ring or trophozoite stage impacts Plasmodium falciparum sexual conversion differently

Malaria transmission is dependent on the formation of gametocytes in the human blood. The sexual conversion rate, the proportion of asexual parasites that convert into gametocytes at each multiplication cycle, is variable and reflects the relative parasite investment between transmission and maintaining the infection. The impact of environmental factors such as drugs on sexual conversion rates is not well understood. We developed a robust assay using gametocyte-reporter parasite lines to accurately measure the impact of drugs on sexual conversion rates, independently from their gametocytocidal activity. We found that exposure to subcurative doses of the frontline antimalarial drug dihydroartemisinin (DHA) at the trophozoite stage resulted in a ~ fourfold increase in sexual conversion. In contrast, no increase was observed when ring stages were exposed or in cultures in which sexual conversion was stimulated by choline depletion. Our results reveal a complex relationship between antimalarial drugs and sexual conversion, with potential public health implications.

A modified two-color flow cytometry assay to quantify in-vitro reinvasion and determine invasion phenotypes at low Plasmodium falciparum parasitemia

Invasion of human red blood cells (RBCs) by Plasmodium parasites is a crucial yet poorly characterised phenotype. Two-color flow cytometry (2cFCM) promises to be a very sensitive and high throughput method for phenotyping parasite invasion. However, current protocols require high (~1.0%) parasitemia for assay set-up and need to be adapted for low parasitemia samples, which are becoming increasingly common in low transmission settings. Background fluorescence from nuclei-containing uninfected RBCs and high autologous reinvasion rates (merozoite invasion of donor uninfected RBCs present at 50% assay volume) are some of the limitations to the method's sensitivity to enumerate low parasitemia (<0.5%) with nucleic acid-based stains. Here, we describe modifications for plating unlabeled donor to labeled target RBCs per assay well and for gating parasitemia, that produces accurate quantifications of low reinvasion parasitemia. Plasmodium falciparum 3D7, Dd2 and field isolates at various low and high parasitemia (0.05%-2.0%) were used to set-up SyBr Green 1-based 2cFCM invasion assays. Target RBCs were labeled with CTFR proliferation dye. We show that this dye combination allowed for efficient parasite invasion into target RBCs and that a 1:3 ratio of unlabeled to labeled RBCs per assay greatly skewed autologous reinvasion (p < 0.001). Accuracy of quantifying reinvasion was limited to an assay parasitemia of 0.02% with minimal background interference. Invasion inhibition by enzymatic treatments increased averagely by 10% (p<0.05) across the entire parasitemia range. The effect was greater for samples with <0.5% parasitemia. Overall, a more sensitive method for phenotyping invasion of low P. falciparum parasitemia is described.

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.

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

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

Reporter lines based on the gexp02 promoter enable early quantification of sexual conversion rates in the malaria parasite Plasmodium falciparum

Transmission of malaria parasites from humans to mosquito vectors requires that some asexual parasites differentiate into sexual forms termed gametocytes. The balance between proliferation in the same host and conversion into transmission forms can be altered by the conditions of the environment. The ability to accurately measure the rate of sexual conversion under different conditions is essential for research addressing the mechanisms underlying sexual conversion, and to assess the impact of environmental factors. Here we describe new Plasmodium falciparum transgenic lines with genome-integrated constructs in which a fluorescent reporter is expressed under the control of the promoter of the gexp02 gene. Using these parasite lines, we developed a sexual conversion assay that shortens considerably the time needed for an accurate determination of sexual conversion rates, and dispenses the need to add chemicals to inhibit parasite replication. Furthermore, we demonstrate that gexp02 is expressed specifically in sexual parasites, with expression starting as early as the sexual ring stage, which makes it a candidate marker for circulating sexual rings in epidemiological studies.

Revisiting the initial steps of sexual development in the malaria parasite Plasmodium falciparum

Human to vector transmission of malaria requires that some blood-stage parasites abandon asexual growth and convert into non-replicating sexual forms called gametocytes. The initial steps of gametocytogenesis remain largely uncharacterized. Here, we study this part of the malaria life cycle in Plasmodium falciparum using PfAP2-G, the master regulator of sexual conversion, as a marker of commitment. We demonstrate the existence of PfAP2-G-positive sexually committed parasite stages that precede the previously known committed schizont stage. We also found that sexual conversion can occur by two different routes: the previously described route in which PfAP2-G-expressing parasites complete a replicative cycle as committed forms before converting into gametocytes upon re-invasion, or a direct route with conversion within the same cycle as initial PfAP2-G expression. The latter route is linked to early PfAP2-G expression in ring stages. Reanalysis of published single-cell RNA-sequencing (RNA-seq) data confirmed the presence of both routes. Consistent with these results, using plaque assays we observed that, in contrast to the prevailing model, many schizonts produced mixed plaques containing both asexual parasites and gametocytes. Altogether, our results reveal unexpected features of the initial steps of sexual development and extend the current view of this part of the malaria life cycle.

The Plasmodium falciparum transcriptome in severe malaria reveals altered expression of genes involved in important processes including surface antigen-encoding var genes

Within the human host, the malaria parasite Plasmodium falciparum is exposed to multiple selection pressures. The host environment changes dramatically in severe malaria, but the extent to which the parasite responds to-or is selected by-this environment remains unclear. From previous studies, the parasites that cause severe malaria appear to increase expression of a restricted but poorly defined subset of the PfEMP1 variant, surface antigens. PfEMP1s are major targets of protective immunity. Here, we used RNA sequencing (RNAseq) to analyse gene expression in 44 parasite isolates that caused severe and uncomplicated malaria in Papuan patients. The transcriptomes of 19 parasite isolates associated with severe malaria indicated that these parasites had decreased glycolysis without activation of compensatory pathways; altered chromatin structure and probably transcriptional regulation through decreased histone methylation; reduced surface expression of PfEMP1; and down-regulated expression of multiple chaperone proteins. Our RNAseq also identified novel associations between disease severity and PfEMP1 transcripts, domains, and smaller sequence segments and also confirmed all previously reported associations between expressed PfEMP1 sequences and severe disease. These findings will inform efforts to identify vaccine targets for severe malaria and also indicate how parasites adapt to-or are selected by-the host environment in severe malaria.

Multiplication rate variation in the human malaria parasite Plasmodium falciparum

It is important to understand intrinsic variation in asexual blood stage multiplication rates of the most virulent human malaria parasite, Plasmodium falciparum. Here, multiplication rates of long-term laboratory adapted parasite clones and new clinical isolates were measured, using a newly standardised assay of growth from low starting density in replicate parallel cultures with erythrocytes from multiple different donors, across multiple cycles. Multiplication rates of long-term established clones were between 7.6 and 10.5 fold per 48 hours, with clone Dd2 having a higher rate than others (clones 3D7, HB3 and D10). Parasite clone-specific growth was then analysed in co-culture assays with all possible heterologous pairwise combinations. This showed that co-culture of different parasites did not affect their replication rates, indicating that there were no suppressive interactions operating between parasites. Multiplication rates of eleven new clinical isolates were measured after a few weeks of culture, and showed a spectrum of replication rates between 2.3 and 6.0 fold per 48 hours, the entire range being lower than for the long-term laboratory adapted clones. Multiplication rate estimates remained stable over time for several isolates tested repeatedly up to three months after culture initiation, indicating considerable persistence of this important trait variation.