Transcriptional variation in the malaria parasite Plasmodium falciparum

The Plasmodium falciparum histone methyltransferase PfSET10 is dispensable for the regulation of antigenic variation and gene expression in blood-stage parasites

The malaria parasite Plasmodium falciparum employs antigenic variation of the virulence factor P. falciparum erythrocyte membrane protein 1 (PfEMP1) to escape adaptive immune responses during blood infection. Antigenic variation of PfEMP1 occurs through epigenetic switches in the mutually exclusive expression of individual members of the multi-copy var gene family. var genes are located in perinuclear clusters of transcriptionally inactive heterochromatin. Singular var gene activation is linked to locus repositioning into a dedicated zone at the nuclear periphery and deposition of histone 3 lysine 4 di-/trimethylation (H3K4me2/3) and H3K9 acetylation marks in the promoter region. While previous work identified the putative H3K4-specific methyltransferase PfSET10 as an essential enzyme and positive regulator of var gene expression, a recent study reported conflicting data. Here, we used iterative genome editing to engineer a conditional PfSET10 knockout line tailored to study the function of PfSET10 in var gene regulation. We demonstrate that PfSET10 is not required for mutually exclusive var gene expression and switching. We also show that PfSET10 is dispensable not only for asexual parasite proliferation but also for sexual conversion and gametocyte differentiation. Furthermore, comparative RNA-seq experiments revealed that PfSET10 plays no obvious role in regulating gene expression during asexual parasite development and gametocytogenesis. Interestingly, however, PfSET10 shows different subnuclear localization patterns in asexual and sexual stage parasites and female-specific expression in mature gametocytes. In summary, our work confirms in detail that PfSET10 is not involved in regulating var gene expression and is not required for blood-stage parasite viability, indicating PfSET10 may be important for life cycle progression in the mosquito vector or during liver stage development.IMPORTANCEThe malaria parasite Plasmodium falciparum infects hundreds of millions of people every year. To survive and proliferate in the human bloodstream, the parasites need to escape recognition by the host's immune system. To achieve this, P. falciparum can change the expression of surface antigens via a process called antigenic variation. This fascinating survival strategy is based on infrequent switches in the expression of single members of the var multigene family. Previous research reported conflicting results on the role of the epigenetic regulator PfSET10 in controlling mutually exclusive var gene expression and switching. Here, we unequivocally demonstrate that PfSET10 is neither required for antigenic variation nor the expression of any other proteins during blood-stage infection. This information is critical in directing our attention toward exploring alternative molecular mechanisms underlying the control of antigenic variation and investigating the function of PfSET10 in other life cycle stages.

Heterochromatin dynamics during the initial stages of sexual development in Plasmodium falciparum

Asexual replication of Plasmodium falciparum in the human blood results in exponential parasite growth and causes all clinical symptoms of malaria. However, at each round of the replicative cycle, some parasites convert into sexual precursors called gametocytes, which develop through different stages until they become infective to mosquito vectors. The genome-wide distribution of heterochromatin, a type of chromatin generally refractory to gene expression, is identical at all asexual blood stages, but is altered in stage II/III and more mature gametocytes. However, it is not known if these changes occur concomitantly with sexual conversion or at a later time during gametocyte development. Using a transgenic line in which massive sexual conversion can be conditionally induced, we show that the genome-wide distribution of heterochromatin at the initial stages of sexual development (i.e., sexual rings and stage I gametocytes) is almost identical to asexual blood stages, and major changes do not occur until stage II/III. However, we found that at loci with heterochromatin alterations, transcriptional changes associated with sexual development typically precede, rather than follow, changes in heterochromatin occupancy.

Comparative analysis of peripheral whole blood transcriptome from asymptomatic carriers reveals upregulation of subsets of surface proteins implicated in Plasmodium falciparum phenotypic plasticity

The molecular mechanism underlying Plasmodium falciparum's persistence in the asymptomatic phase of infection remains largely unknown. However, large-scale shifts in the parasites' gene expression during asymptomatic infections may enhance phenotypic plasticity, maximizing their fitness and leading to the persistence of the asymptomatic infections. To uncover these mechanisms, we aimed to identify parasite genetic factors implicated in asymptomatic infections through whole transcriptome analysis. We analyzed publicly available transcriptome datasets containing asymptomatic malaria (ASM), uncomplicated malaria (SM), and malaria-naïve (NSM) samples from 35 subjects for differentially expressed genes (DEGs) and long noncoding RNAs. Our analysis identified 755 and 1773 DEGs in ASM vs SM and NSM, respectively. These DEGs revealed sets of genes coding for proteins of unknown functions (PUFs) upregulated in ASM vs SM and ASM, suggesting their role in underlying fundamental molecular mechanisms during asymptomatic infections. Upregulated genes in ASM vs SM revealed a subset of 24 clonal variant genes (CVGs) involved in host-parasite and symbiotic interactions and modulation of the symbiont of host erythrocyte aggregation pathways. Moreover, we identified 237 differentially expressed noncoding RNAs in ASM vs SM, of which 11 were found to interact with CVGs, suggesting their possible role in regulating the expression of CVGs. Our results suggest that P. falciparum utilizes phenotypic plasticity as an adaptive mechanism during asymptomatic infections by upregulating clonal variant genes, with long noncoding RNAs possibly playing a crucial role in their regulation. Thus, our study provides insights into the parasites' genetic factors that confer a fitness advantage during asymptomatic infections.

Epigenetic regulation as a therapeutic target in the malaria parasite Plasmodium falciparum

Over the past thirty years, epigenetic regulation of gene expression has gained increasing interest as it was shown to be implicated in illnesses ranging from cancers to parasitic diseases. In the malaria parasite, epigenetics was shown to be involved in several key steps of the complex life cycle of Plasmodium, among which asexual development and sexual commitment, but also in major biological processes like immune evasion, response to environmental changes or DNA repair. Because epigenetics plays such paramount roles in the Plasmodium parasite, enzymes involved in these regulating pathways represent a reservoir of potential therapeutic targets. This review focuses on epigenetic regulatory processes and their effectors in the malaria parasite, as well as the inhibitors of epigenetic pathways and their potential as new anti-malarial drugs. Such types of drugs could be formidable tools that may contribute to malaria eradication in a context of widespread resistance to conventional anti-malarials.

Detection of Developmental Asexual Stage-Specific RNA Editing Events in Plasmodium falciparum 3D7 Malaria Parasite

Transcriptional variation has been studied but post-transcriptional modification due to RNA editing has not been investigated in Plasmodium. We investigated developmental stage-specific RNA editing in selected genes in Plasmodium falciparum 3D7. We detected extensive amination- and deamination-type RNA editing at 8, 16, 24, 32, 40, and 46 h in tightly synchronized Plasmodium. Most of the editing events were observed in 8 and 16 h ring-stage parasites. Extensive A-to-G deamination-type editing was detected more during the 16 h ring stage (25%) than the 8 h ring stage (20%). Extensive U-to-C amination-type editing was detected more during the 16 h ring stage (31%) than the 8 h ring stage (22%). In 28S, rRNA editing converted the loop structure to the stem structure. The hemoglobin binding activity of PF3D7_0216900 was also altered due to RNA editing. Among the expressed 28S rRNA genes, PF3D7_0532000 and PF3D7_0726000 expression was higher. Increased amounts of the transcripts of these two genes were found, particularly PF3D7_0726000 in the ring stage and PF3D7_0532000 in the trophozoite and schizont stages. Adenosine deaminase (ADA) expression did not correlate with the editing level. This first experimental report of RNA editing will help to identify the editing machinery that might be useful for antimalarial drug discovery and malaria control.

The many paths to artemisinin resistance in Plasmodium falciparum

Emerging resistance against artemisinin (ART) poses a major challenge in controlling malaria. Parasites with mutations in PfKelch13, the major marker for ART resistance, are known to reduce hemoglobin endocytosis, induce unfolded protein response (UPR), elevate phosphatidylinositol-3-phosphate (PI3P) levels, and stimulate autophagy. Nonetheless, PfKelch13-independent resistance is also reported, indicating extensive complementation by reconfiguration in the parasite metabolome and transcriptome. These findings implicate that there may not be a single 'universal identifier' of ART resistance. This review sheds light on the molecular, transcriptional, and metabolic pathways associated with ART resistance, while also highlighting the interplay between cellular heterogeneity, environmental stress, and ART sensitivity.

The global transcriptome of Plasmodium falciparum mid-stage gametocytes (stages II-IV) appears largely conserved and gametocyte-specific gene expression patterns vary in clinical isolates

Our overall understanding of the developmental biology of malaria parasites has been greatly enhanced by recent advances in transcriptomic analysis. However, most of these investigations rely on laboratory strains (LS) that were adapted into in vitro culture many years ago, and the transcriptomes of clinical isolates (CI) circulating in human populations have not been assessed. In this study, RNA-seq was used to compare the global transcriptome of mid-stage gametocytes derived from three short-term cultured CI, with gametocytes derived from the NF54 reference laboratory strain. The core transcriptome appeared to be consistent between CI- and LS-derived gametocyte preparations, but some important differences were also observed. A majority of gametocyte-specific genes (43/53) appear to have relatively higher expression in CI-derived gametocytes than in LS-derived gametocytes, but a K-means clustering analysis showed that genes involved in flagellum- and microtubule-based processes (movement/motility) were more abundant in both groups, albeit with some differences between them. In addition, gametocytes from one CI described as CI group II gametocytes (CI:GGII) showed gene expression variation in the form of reduced gametocyte-specific gene expression compared to the other two CI-derived gametocytes (CI gametocyte group I, CI:GGI), although the mixed developmental stages used in our study is a potential confounder, only partially mitigated by the inclusion of multiple replicates for each CI. Overall, our study suggests that there may be subtle differences in the gene expression profiles of mid-stage gametocytes from CI relative to the NF54 reference strain of Plasmodium falciparum. Thus, it is necessary to deploy gametocyte-producing clinical parasite isolates to fully understand the diversity of gene expression strategies that may occur during the sequestered development of parasite sexual stages. IMPORTANCE Maturing gametocytes of Plasmodium falciparum are known to sequester away from peripheral circulation into the bone marrow until they are mature. Blocking gametocyte sequestration can prevent malaria transmission from humans to mosquitoes, but most studies aim to understand gametocyte development utilizing long-term adapted laboratory lines instead of clinical isolates. This is a particular issue for our understanding of the sexual stages, which are known to decrease rapidly during adaptation to long-term culture, meaning that many LS are unable to produce transmissible gametocytes. Using RNA-seq, we investigated the global transcriptome of mid-stage gametocytes derived from three clinical isolates and a reference strain (NF54). This identified important differences in gene expression profiles between immature gametocytes of CI and the NF54 reference strain of P. falciparum, suggesting increased investment in gametocytogenesis in clinical isolates. Our transcriptomic data highlight the use of clinical isolates in studying the morphological, cellular features and molecular biology of gametocytes.

Cultivation of Asexual Intraerythrocytic Stages of Plasmodium falciparum

Successfully developed in 1976, the continuous in vitro culture of Plasmodium falciparum has many applications in the field of malaria research. It has become an important experimental model that directly uses a human pathogen responsible for a high prevalence of morbidity and mortality in many parts of the world and is a major source of biological material for immunological, biochemical, molecular, and pharmacological studies. Until present, the basic techniques described by Trager and Jensen and Haynes et al. remain unchanged in many malaria research laboratories. Nonetheless, different factors, including culture media, buffers, serum substitutes and supplements, sources of erythrocytes, and conditions of incubation (especially oxygen concentration), have been modified by different investigators to adapt the original technique in their laboratories or enhance the in vitro growth of the parasites. The possible effects and benefits of these modifications for the continuous cultivation of asexual intraerythrocytic stages of P. falciparum, as well as future challenges in developing a serum-free cultivation system and axenic cultures, are discussed.

The role of long noncoding RNAs in malaria parasites

The human malaria parasites, including Plasmodium falciparum, persist as a major cause of global morbidity and mortality. The recent stalling of progress toward malaria elimination substantiates a need for novel interventions. Controlled gene expression is central to the parasite's numerous life cycle transformations and adaptation. With few specific transcription factors (TFs) identified, crucial roles for chromatin states and epigenetics in parasite transcription have become evident. Although many chromatin-modifying enzymes are known, less is known about which factors mediate their impacts on transcriptional variation. Like those of higher eukaryotes, long noncoding RNAs (lncRNAs) have recently been shown to have integral roles in parasite gene regulation. This review aims to summarize recent developments and key findings on the role of lncRNAs in P. falciparum.

The exception that proves the rule: Virulence gene expression at the onset of Plasmodium falciparum blood stage infections

Controlled human malaria infections (CHMI) are a valuable tool to study parasite gene expression in vivo under defined conditions. In previous studies, virulence gene expression was analyzed in samples from volunteers infected with the Plasmodium falciparum (Pf) NF54 isolate, which is of African origin. Here, we provide an in-depth investigation of parasite virulence gene expression in malaria-naïve European volunteers undergoing CHMI with the genetically distinct Pf 7G8 clone, originating in Brazil. Differential expression of var genes, encoding major virulence factors of Pf, PfEMP1s, was assessed in ex vivo parasite samples as well as in parasites from the in vitro cell bank culture that was used to generate the sporozoites (SPZ) for CHMI (Sanaria PfSPZ Challenge (7G8)). We report broad activation of mainly B-type subtelomeric located var genes at the onset of a 7G8 blood stage infection in naïve volunteers, mirroring the NF54 expression study and suggesting that the expression of virulence-associated genes is generally reset during transmission from the mosquito to the human host. However, in 7G8 parasites, we additionally detected a continuously expressed single C-type variant, Pf7G8_040025600, that was most highly expressed in both pre-mosquito cell bank and volunteer samples, suggesting that 7G8, unlike NF54, maintains expression of some previously expressed var variants during transmission. This suggests that in a new host, the parasite may preferentially express the variants that previously allowed successful infection and transmission. Trial registration: ClinicalTrials.gov - NCT02704533; 2018-004523-36.

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 rising tide of parasite transcriptomics propels pathogen biology

Twenty years ago, the first transcriptome of the intraerythrocytic developmental cycle of the malaria parasite Plasmodium falciparum was published in PLOS Biology. Since then, transcriptomics studies have transformed the study of parasite biology.

Detection of positive selection acting on protein surfaces at the whole-genome scale in the human malaria parasite Plasmodium falciparum

The host-parasite evolutionary arms race is a fundamental process with medical implications. During this process, the host develops parasite resistance, and the parasite develops host immune evasion strategies. Thus, this process accelerates relevant protein evolution. This study test hypothesizes that proteins subject to sequence evolution structural constraints play a crucial role and that these constraints hinder the modification of such proteins in this process. These hypotheses were tested using Plasmodium falciparum model and evaluated protein structures predicted for the entire proteome by the AlphaFold method. Based on dN/dS test results and P. falciparum and P. reichenowi comparisons, the presented approach identified proteins subject to purifying selection acting on the whole sequence and buried residues (dN < dS) and positive selection on nonburied residues. Of the 26 proteins, some known antigens (ring-exported protein 3, RAP protein, erythrocyte binding antigen-140, and protein P47) targeted by the host immune system are promising vaccine candidates. The set also contained 11 enzymes, including FIKK kinase, which modifies host proteins. This set was compared with genes for which the dN/dS test suggested that positive selection acts on the whole gene (i.e., dN > dS). The present study found that such genes encode enzymes and antigenic vaccine candidates less frequently than genes for which evolution is not subject to selection constraints and positive selection acts on only exposed residues. The analysis was repeated comparing P. falciparum with P. alderi, which is more distantly related. The study discusses the potential implications of the presented methodology for rational vaccine design and the parasitology and evolutionary biology fields.

Malian children infected with Plasmodium ovale and Plasmodium falciparum display very similar gene expression profiles

Plasmodium parasites caused 241 million cases of malaria and over 600,000 deaths in 2020. Both P. falciparum and P. ovale are endemic to Mali and cause clinical malaria, with P. falciparum infections typically being more severe. Here, we sequenced RNA from nine pediatric blood samples collected during infections with either P. falciparum or P. ovale, and characterized the host and parasite gene expression profiles. We found that human gene expression varies more between individuals than according to the parasite species causing the infection, while parasite gene expression profiles cluster by species. Additionally, we characterized DNA polymorphisms of the parasites directly from the RNA-seq reads and found comparable levels of genetic diversity in both species, despite dramatic differences in prevalence. Our results provide unique insights into host-pathogen interactions during malaria infections and their variations according to the infecting Plasmodium species, which will be critical to develop better elimination strategies against all human Plasmodium parasites.

Plasmodium falciparum Sexual Commitment Rate Variation among Clinical Isolates and Diverse Laboratory-Adapted Lines

Asexual blood-stage malaria parasites must produce sexual progeny to infect mosquitoes. It is important to understand the scope and causes of intraspecific variation in sexual commitment rates, particularly for the major human parasite P. falciparum. First, two alternative assay methods of measuring sexual commitment were compared to test a genetically modified P. falciparum line with elevated commitment rates inducible by overexpression of GDV1. The methods yielded correlated measurements with higher sensitivity and precision being achieved by one employing detection of the early gametocyte differentiation marker Pfs16. Thus, this was used to survey a diverse range of parasite lines and test each in multiple biological replicate assays in a serum-free medium supplemented with Albumax. There were differences among six recent clinical isolates from Ghana in their mean rates of sexual commitment per cycle, ranging from 3.3% to 12.2%. Among 13 diverse long-term laboratory-adapted lines, mean sexual commitment rates for most ranged from 4.7% to 13.4%, a few had lower rates with means from 0.3 to 1.6%, and one with a nonfunctional ap2-g gene always showed zero commitment. Among a subset of lines tested for the effects of exogenous choline to suppress commitment, there were significant differences. As expected, there was no effect in a line that had lost the gdv1 gene and that had generally low commitment, whereas the others showed quantitatively variable but significant responses to choline, suggesting potential trait variation. The results indicated the value of performing multiple replicate assays for understanding the variation of this key reproductive trait that likely affects transmission. IMPORTANCE Only sexual-stage malaria parasites are transmitted from human blood to mosquitoes. Thus, it is vital to understand variations in sexual commitment rates because these may be modifiable or susceptible to blocking. Two different methods of commitment rate measurement were first compared, demonstrating higher sensitivity and precision by the detection of an early differentiation marker, which was subsequently used to survey diverse lines. Clinical isolates from Ghana showed significant variation in mean per-cycle commitment rates and variation among biological replicates. Laboratory-adapted lines of diverse origins had a wider range with most being within the range observed for the clinical isolates, while a minority consistently had lower or zero rates. There was quantitative variation in the effects when adding choline to suppress commitment, indicating differing responsiveness of parasites to this environmental modification. Performing multiple assay replicates and comparisons of diverse isolates was important to understand this trait and its potential effects on transmission.

Naturally Acquired Kelch13 Mutations in Plasmodium falciparum Strains Modulate In Vitro Ring-Stage Artemisinin-Based Drug Tolerance and Parasite Survival in Response to Hyperoxia

The ring-stage survival assay was utilized to assess the impact of physiological hyperoxic stress on dihydroartemisinin (DHA) tolerance for a panel of Plasmodium falciparum strains with and without Kelch13 mutations. Strains without naturally acquired Kelch13 mutations or the postulated genetic background associated with delayed parasite clearance time demonstrated reduced proliferation under hyperoxic conditions in the subsequent proliferation cycle. Dihydroartemisinin tolerance in three isolates with naturally acquired Kelch13 mutations but not two genetically manipulated laboratory strains was modulated by in vitro hyperoxic stress exposure of early-ring-stage parasites in the cycle before drug exposure. Reduced parasite tolerance to additional derivatives, including artemisinin, artesunate, and OZ277, was observed within the second proliferation cycle. OZ439 and epoxomicin completely prevented parasite survival under both hyperoxia and normoxic in vitro culture conditions, highlighting the unique relationship between DHA tolerance and Kelch13 mutation-associated genetic background. IMPORTANCE Artemisinin-based combination therapy (ACT) for treating malaria is under intense scrutiny following treatment failures in the Greater Mekong subregion of Asia. This is further compounded by the potential for extensive loss of life if treatment failures extend to the African continent. Although Plasmodium falciparum has become resistant to all antimalarial drugs, artemisinin "resistance" does not present in the same way as resistance to other antimalarial drugs. Instead, a partial resistance or tolerance is demonstrated, associated with the parasite's genetic profile and linked to a molecular marker referred to as K13. It is suggested that parasites may have adapted to drug treatment, as well as the presence of underlying population health issues such as hemoglobinopathies, and/or environmental pressures, resulting in parasite tolerance to ACT. Understanding parasite evolution and control of artemisinin tolerance will provide innovative approaches to mitigate the development of artemisinin tolerance and thereby artemisinin-based drug treatment failure and loss of life globally to malaria infections.

Analysis of Plasmodium falciparum myosin B ATPase activity and structure in complex with the calmodulin-like domain of its light chain MLC-B

Myosin B (MyoB) is a class 14 myosin expressed in all invasive stages of the malaria parasite, Plasmodium falciparum. It is not associated with the glideosome complex that drives motility and invasion of host cells. During red blood cell invasion, MyoB remains at the apical tip of the merozoite but is no longer observed once invasion is completed. MyoB is not essential for parasite survival, but when it is knocked out, merozoites are delayed in the initial stages of red blood cell invasion, giving rise to a growth defect that correlates with reduced invasion success. Therefore, further characterization is needed to understand how MyoB contributes to parasite invasion. Here, we have expressed and purified functional MyoB with the help of parasite-specific chaperones Hsp90 and Unc45, characterized its binding to actin and its known light chain MLC-B using biochemical and biophysical methods and determined its low-resolution structure in solution using small angle X-ray scattering. In addition to MLC-B, we found that four other putative regulatory light chains bind to the MyoB IQ2 motif in vitro. The purified recombinant MyoB adopted the overall shape of a myosin, exhibited actin-activated ATPase activity, and moved actin filaments in vitro. Additionally, we determined that the ADP release rate was faster than the ATP turnover number, and thus, does not appear to be rate limiting. This, together with the observed high affinity to actin and the specific localization of MyoB, may point toward a role in tethering and/or force sensing during early stages of invasion.

Epigenetic and Epitranscriptomic Gene Regulation in Plasmodium falciparum and How We Can Use It against Malaria

Malaria, caused by Plasmodium parasites, is still one of the biggest global health challenges. P. falciparum is the deadliest species to humans. In this review, we discuss how this parasite develops and adapts to the complex and heterogenous environments of its two hosts thanks to varied chromatin-associated and epigenetic mechanisms. First, one small family of transcription factors, the ApiAP2 proteins, functions as master regulators of spatio-temporal patterns of gene expression through the parasite life cycle. In addition, chromatin plasticity determines variable parasite cell phenotypes that link to parasite growth, virulence and transmission, enabling parasite adaptation within host conditions. In recent years, epitranscriptomics is emerging as a new regulatory layer of gene expression. We present evidence of the variety of tRNA and mRNA modifications that are being characterized in Plasmodium spp., and the dynamic changes in their abundance during parasite development and cell fate. We end up outlining that new biological systems, like the mosquito model, to decipher the unknowns about epigenetic mechanisms in vivo; and novel methodologies, to study the function of RNA modifications; are needed to discover the Achilles heel of the parasite. With this new knowledge, future strategies manipulating the epigenetics and epitranscriptomic machinery of the parasite have the potential of providing new weapons against malaria.

Pervasive sequence-level variation in the transcriptome of Plasmodium falciparum

Single-nucleotide variations (SNVs) in RNA, arising from co- and post-transcriptional phenomena including transcription errors and RNA-editing, are well studied in a range of organisms. In the malaria parasite Plasmodium falciparum, stage-specific and non-specific gene-expression variations accompany the parasite's array of developmental and morphological phenotypes over the course of its complex life cycle. However, the extent, rate and effect of sequence-level variation in the parasite's transcriptome are unknown. Here, we report the presence of pervasive, non-specific SNVs in the P. falciparum transcriptome. SNV rates for a gene were correlated to gene length (r[Formula: see text]0.65-0.7) but not to the AT-content of that gene. Global SNV rates for the P. falciparum lines we used, and for publicly available P. vivax and P. falciparum clinical isolate datasets, were of the order of 10-3 per base, ∼10× higher than rates we calculated for bacterial datasets. These variations may reflect an intrinsic transcriptional error rate in the parasite, and RNA editing may be responsible for a subset of them. This seemingly characteristic property of the parasite may have implications for clinical outcomes and the basic biology and evolution of P. falciparum and parasite biology more broadly. We anticipate that our study will prompt further investigations into the exact sources, consequences and possible adaptive roles of these SNVs.

RhopH2 and RhopH3 export enables assembly of the RhopH complex on P. falciparum-infected erythrocyte membranes

RhopH complexes consists of Clag3, RhopH2 and RhopH3 and are essential for growth of Plasmodium falciparum inside infected erythrocytes. Proteins are released from rhoptry organelles during merozoite invasion and trafficked to the surface of infected erythrocytes and enable uptake of nutrients. RhopH3, unlike other RhopH proteins, is required for parasite invasion, suggesting some cellular processes RhopH proteins function as single players rather than a complex. We show the RhopH complex has not formed during merozoite invasion. Clag3 is directly released into the host cell cytoplasm, whilst RhopH2 and RhopH3 are released into the nascent parasitophorous vacuole. Export of RhopH2 and RhopH3 from the parasitophorous vacuole into the infected erythrocyte cytoplasm enables assembly of Clag3/RhopH2/RhopH3 complexes and incorporation into the host cell membrane concomitant with activation of nutrient uptake. This suggests compartmentalisation prevents premature channel assembly before intact complex is assembled at the host cell membrane.

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.

Histone Modification Landscapes as a Roadmap for Malaria Parasite Development

Plasmodium falciparum remains the deadliest parasite species in the world, responsible for 229 million cases of human malaria in 2019. The ability of the P. falciparum parasite to progress through multiple life cycle stages and thrive in diverse host and vector species hinges on sophisticated mechanisms of epigenetic regulation of gene expression. Emerging evidence indicates such epigenetic control exists in concentric layers, revolving around core histone post-translational modification (PTM) landscapes. Here, we provide a necessary update of recent epigenome research in malaria parasites, focusing specifically on the ability of dynamic histone PTM landscapes to orchestrate the divergent development and differentiation pathways in P. falciparum parasites. In addition to individual histone PTMs, we discuss recent findings that imply functional importance for combinatorial PTMs in P. falciparum parasites, representing an operational histone code. Finally, this review highlights the remaining gaps and provides strategies to address these to obtain a more thorough understanding of the histone modification landscapes that are at the center of epigenetic regulation in human malaria parasites.

Malaria parasites do respond to heat

The capacity of malaria parasites to respond to changes in their environment at the transcriptional level has been the subject of debate, but recent evidence has unambiguously demonstrated that Plasmodium spp. can produce adaptive transcriptional responses when exposed to some specific types of stress. These include metabolic conditions and febrile temperature. The Plasmodium falciparum protective response to thermal stress is similar to the response in other organisms, but it is regulated by a transcription factor evolutionarily unrelated to the conserved transcription factor that drives the heat shock (HS) response in most eukaryotes. Of the many genes that change expression during HS, only a subset constitutes an authentic response that contributes to parasite survival.

Contribution of Transcriptome to Elucidate the Biology of Plasmodium spp

In the present review, we discuss some of the new technologies that have been applied to elucidate how Plasmodium spp escape from the immune system and subvert the host physiology to orchestrate the regulation of its biological pathways. Our manuscript describes how techniques such as microarray approaches, RNA-Seq and single-cell RNA sequencing have contributed to the discovery of transcripts and changed the concept of gene expression regulation in closely related malaria parasite species. Moreover, the text highlights the contributions of high-throughput RNA sequencing for the current knowledge of malaria parasite biology, physiology, vaccine target and the revelation of new players in parasite signaling.

The 3-phosphoinositide-dependent protein kinase 1 is an essential upstream activator of protein kinase A in malaria parasites

Cyclic adenosine monophosphate (cAMP)-dependent protein kinase A (PKA) signalling is essential for the proliferation of Plasmodium falciparum malaria blood stage parasites. The mechanisms regulating the activity of the catalytic subunit PfPKAc, however, are only partially understood, and PfPKAc function has not been investigated in gametocytes, the sexual blood stage forms that are essential for malaria transmission. By studying a conditional PfPKAc knockdown (cKD) mutant, we confirm the essential role for PfPKAc in erythrocyte invasion by merozoites and show that PfPKAc is involved in regulating gametocyte deformability. We furthermore demonstrate that overexpression of PfPKAc is lethal and kills parasites at the early phase of schizogony. Strikingly, whole genome sequencing (WGS) of parasite mutants selected to tolerate increased PfPKAc expression levels identified missense mutations exclusively in the gene encoding the parasite orthologue of 3-phosphoinositide-dependent protein kinase-1 (PfPDK1). Using targeted mutagenesis, we demonstrate that PfPDK1 is required to activate PfPKAc and that T189 in the PfPKAc activation loop is the crucial target residue in this process. In summary, our results corroborate the importance of tight regulation of PfPKA signalling for parasite survival and imply that PfPDK1 acts as a crucial upstream regulator in this pathway and potential new drug target.

Single-Cell RNA Sequencing Reveals Cellular Heterogeneity and Stage Transition under Temperature Stress in Synchronized Plasmodium falciparum Cells

The malaria parasite has a complex life cycle exhibiting phenotypic and morphogenic variations in two different hosts by existing in heterogeneous developmental states. To investigate this cellular heterogeneity of the parasite within the human host, we performed single-cell RNA sequencing of synchronized Plasmodium cells under control and temperature treatment conditions. Using the Malaria Cell Atlas (https://www.sanger.ac.uk/science/tools/mca) as a guide, we identified 9 subtypes of the parasite distributed across known intraerythrocytic stages. Interestingly, temperature treatment results in the upregulation of the AP2-G gene, the master regulator of sexual development in a small subpopulation of the parasites. Moreover, we identified a heterogeneous stress-responsive subpopulation (clusters 5, 6, and 7 [∼10% of the total population]) that exhibits upregulation of stress response pathways under normal growth conditions. We also developed an online exploratory tool that will provide new insights into gene function under normal and temperature stress conditions. Thus, our study reveals important insights into cell-to-cell heterogeneity in the parasite population under temperature treatment that will be instrumental toward a mechanistic understanding of cellular adaptation and population dynamics in Plasmodium falciparum. IMPORTANCE The malaria parasite has a complex life cycle exhibiting phenotypic variations in two different hosts accompanied by cell-to-cell variability that is important for stress tolerance, immune evasion, and drug resistance. To investigate cellular heterogeneity determined by gene expression, we performed single-cell RNA sequencing (scRNA-seq) of about 12,000 synchronized Plasmodium cells under physiologically relevant normal (37°C) and temperature stress (40°C) conditions phenocopying the cyclic bouts of fever experienced during malarial infection. In this study, we found that parasites exhibit transcriptional heterogeneity in an otherwise morphologically synchronized culture. Also, a subset of parasites is continually committed to gametocytogenesis and stress-responsive pathways. These observations have important implications for understanding the mechanisms of drug resistance generation and vaccine development against the malaria parasite.

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.

Plasmodium falciparum DDX3X is a nucleocytoplasmic protein and requires N-terminal for DNA helicase activity

Malaria is a haemato-protozoan disease which causes thousands of deaths every year. Due to the alarming increase of drug resistant strains of P. falciparum, malaria is now becoming more deadly. Helicases are the most important components of the cellular machinery without which cells are unable to survive. The importance of helicases has been proven in variety of organisms. In this study we have reported detailed biochemical characterization of human homologue of DDX3X from Plasmodium falciparum (PfDDX3X). Our study revealed that PfDDX3X is ATP- dependent DNA helicase whereas in human host it is ATP-dependent RNA helicase. We show that N-terminal is essential for its activity and it is present in nucleus and cytoplasm in intraerythrocytic developmental stages of P. falciparum 3D7 strain. Also, it is highly expressed in the schizont stage of P. falciparum 3D7strain. The present study suggests that a protein can perform different functions in different systems. The present study will help to understand the basic biology of malaria parasite P. falciparum.

The private life of malaria parasites: Strategies for sexual reproduction

Malaria parasites exhibit a complex lifecycle, requiring extensive asexual replication in the liver and blood of the vertebrate host, and in the haemocoel of the insect vector. Yet, they must also undergo a single round of sexual reproduction, which occurs in the vector's midgut upon uptake of a blood meal. Sexual reproduction is obligate for infection of the vector and thus, is essential for onwards transmission to new hosts. Sex in malaria parasites involves several bottlenecks in parasite number, making the stages involved attractive targets for blocking disease transmission. Malaria parasites have evolved a suite of adaptations ("strategies") to maximise the success of sexual reproduction and transmission, which could undermine transmission-blocking interventions. Yet, understanding parasite strategies may also reveal novel opportunities for such interventions. Here, we outline how evolutionary and ecological theories, developed to explain reproductive strategies in multicellular taxa, can be applied to explain two reproductive strategies (conversion rate and sex ratio) expressed by malaria parasites within the vertebrate host.

Identification of a protein unique to the genus Plasmodium that contains a WD40 repeat domain and extensive low-complexity sequence

Proteins containing WD40 domains play important roles in the formation of multiprotein complexes. Little is known about WD40 proteins in the malaria parasite. This report contains the initial description of a WD40 protein that is unique to the genus Plasmodium and possibly closely related genera. The N-terminal portion of this protein consists of seven WD40 repeats that are highly conserved in all Plasmodium species. Following the N-terminal region is a central region that is conserved within the major Plasmodium clades, such as parasites of great apes, monkeys, rodents, and birds, but partially conserved across all Plasmodium species. This central region contains extensive low-complexity sequence and is predicted to have a disordered structure. Proteins with disordered structure generally function in molecular interactions. The C-terminal region is semi-conserved across all Plasmodium species and has no notable features. This WD40 repeat protein likely functions in some aspect of parasite biology that is unique to Plasmodium and this uniqueness makes the protein a possible target for therapeutic intervention.

Transcriptome assembly and differential gene expression of the invasive avian malaria parasite Plasmodium relictum in Hawai'i

The malaria parasite Plasmodium relictum (lineage GRW4) was introduced less than a century ago to the native avifauna of Hawai'i, where it has since caused major declines of endemic bird populations. One of the native bird species that is frequently infected with GRW4 is the Hawai'i 'amakihi (Chlorodrepanis virens). To achieve a better understanding of the transcriptional activities of this virulent parasite, we performed a controlled challenge experiment of 15 'amakihi that were infected with GRW4. Blood samples containing malaria parasites were collected at two time points (intermediate and peak infection stages) from host individuals that were either experimentally infected by mosquitoes or inoculated with infected blood. We then used RNA sequencing to assemble a high-quality blood transcriptome of P. relictum GRW4, allowing us to quantify parasite expression levels inside individual birds. We found few significant differences (one to two transcripts) in GRW4 expression levels between host infection stages and between inoculation methods. However, 36 transcripts showed differential expression levels among all host individuals, indicating a potential presence of host-specific gene regulation across hosts. To reduce the extinction risk of the remaining native bird species in Hawai'i, genetic resources of the local Plasmodium lineage are needed to enable further molecular characterization of this parasite. Our newly built Hawaiian GRW4 transcriptome assembly, together with analyses of the parasite's transcriptional activities inside the blood of Hawai'i 'amakihi, can provide us with important knowledge on how to combat this deadly avian disease in the future.

The PfAP2-G2 transcription factor is a critical regulator of gametocyte maturation

Differentiation from asexual blood stages to mature sexual gametocytes is required for the transmission of malaria parasites. Here, we report that the ApiAP2 transcription factor, PfAP2-G2 (PF3D7_1408200) plays a critical role in the maturation of Plasmodium falciparum gametocytes. PfAP2-G2 binds to the promoters of a wide array of genes that are expressed at many stages of the parasite life cycle. Interestingly, we also find binding of PfAP2-G2 within the gene body of almost 3,000 genes, which strongly correlates with the location of H3K36me3 and several other histone modifications as well as Heterochromatin Protein 1 (HP1), suggesting that occupancy of PfAP2-G2 in gene bodies may serve as an alternative regulatory mechanism. Disruption of pfap2-g2 does not impact asexual development, but the majority of sexual parasites are unable to mature beyond stage III gametocytes. The absence of pfap2-g2 leads to overexpression of 28% of the genes bound by PfAP2-G2 and none of the PfAP2-G2 bound genes are downregulated, suggesting that it is a repressor. We also find that PfAP2-G2 interacts with chromatin remodeling proteins, a microrchidia (MORC) protein, and another ApiAP2 protein (PF3D7_1139300). Overall our data demonstrate that PfAP2-G2 establishes an essential gametocyte maturation program in association with other chromatin-related proteins.

Plasmodium falciparum DDX17 is an RNA helicase crucial for parasite development

Malaria is one of the major global health concerns still prevailing in this 21st century. Even the effect of artemisinin combination therapies (ACT) have declined and causing more mortality across the globe. Therefore, it is important to understand the basic biology of malaria parasite in order to find novel drug targets. Helicases play important role in nucleic acid metabolism and are components of cellular machinery in various organisms. In this manuscript we have performed the biochemical characterization of homologue of DDX17 from Plasmodium falciparum (PfDDX17). Our results show that PfDDX17 is an active RNA helicase and uses mostly ATP for its function. The qRT-PCR experiment results suggest that PfDDX17 is highly expressed in the trophozoite stage and it is localised mainly in the cytoplasm and in infected RBC (iRBC) membrane mostly in the trophozoite stage. The dsRNA knockdown study suggests that PfDDX17 is important for cell cycle progression. These studies report the biochemical functions of PfDDX17 helicase and further augment the fundamental knowledge about helicase families of P. falciparum.

•

Biochemical characterization of homologue of DDX17 from Plasmodium falciparum (PfDDX17) is presented.

•

Results show that PfDDX17 is an active RNA helicase and uses mostly ATP for its function.

•

Results also suggest that PfDDX17 is highly expressed in the trophozoite stage.

•

dsRNA knockdown study revealed that PfDDX17 is important for cell cycle progression.

The ApiAP2 factor PfAP2-HC is an integral component of heterochromatin in the malaria parasite Plasmodium falciparum

Malaria parasites undergo a complex life cycle in the human host and the mosquito vector. The ApiAP2 family of DNA-binding proteins plays a dominant role in parasite development and life cycle progression. Most ApiAP2 factors studied to date act as transcription factors regulating stage-specific gene expression. Here, we characterized an ApiAP2 factor in Plasmodium falciparum that we termed PfAP2-HC. We demonstrate that PfAP2-HC specifically binds to heterochromatin throughout the genome. Intriguingly, PfAP2-HC does not bind DNA in vivo and recruitment of PfAP2-HC to heterochromatin is independent of its DNA-binding domain but strictly dependent on heterochromatin protein 1. Furthermore, our results suggest that PfAP2-HC functions neither in the regulation of gene expression nor in heterochromatin formation or maintenance. In summary, our findings reveal PfAP2-HC as a core component of heterochromatin in malaria parasites and identify unexpected properties and substantial functional divergence among the members of the ApiAP2 family of regulatory proteins.

•

The ApiAP2 factor AP2-HC is a core component of heterochromatin in malaria parasites

•

Binding of AP2-HC to heterochromatin strictly depends on heterochromatin protein 1

•

The AP2 DNA-binding domain of AP2-HC is dispensable for heterochromatin association

Borgheti-Cardoso L, Albertazzi L, Dimova R, Fernàndez-Busquets X. The ESCRT-III machinery participates in the production of extracellular vesicles and protein export during Plasmodium falciparum infection

Infection with Plasmodium falciparum enhances extracellular vesicle (EV) production in parasitized red blood cells (pRBCs), an important mechanism for parasite-to-parasite communication during the asexual intraerythrocytic life cycle. The endosomal sorting complex required for transport (ESCRT), and in particular the ESCRT-III sub-complex, participates in the formation of EVs in higher eukaryotes. However, RBCs have lost the majority of their organelles through the maturation process, including an important reduction in their vesicular network. Therefore, the mechanism of EV production in P. falciparum-infected RBCs remains to be elucidated. Here we demonstrate that P. falciparum possesses a functional ESCRT-III machinery activated by an alternative recruitment pathway involving the action of PfBro1 and PfVps32/PfVps60 proteins. Additionally, multivesicular body formation and membrane shedding, both reported mechanisms of EV production, were reconstituted in the membrane model of giant unilamellar vesicles using the purified recombinant proteins. Moreover, the presence of PfVps32, PfVps60 and PfBro1 in EVs purified from a pRBC culture was confirmed by super-resolution microscopy and dot blot assays. Finally, disruption of the PfVps60 gene led to a reduction in the number of the produced EVs in the KO strain and affected the distribution of other ESCRT-III components. Overall, our results increase the knowledge on the underlying molecular mechanisms during malaria pathogenesis and demonstrate that ESCRT-III P. falciparum proteins participate in EV production.

Malaria is a disease caused by Plasmodium parasites that is still a leading cause of death in many low-income countries, and for which currently available therapeutic strategies are not succeeding in its control, let alone eradication. An interesting feature observed after Plasmodium invasion is the increase of extracellular vesicles (EVs) generated by parasitized red blood cells (pRBCs), which lack a vesicular trafficking that would explain EV production. Here, by combining different approaches, we demonstrated the participation of the endosomal sorting complex required for transport (ESCRT) machinery from Plasmodium falciparum in the production of EVs in pRBCs. Moreover, we were able to detect ESCRT-III proteins adjacent to the membrane of the host and in EVs purified from a pRBC culture, which shows the export of these proteins and their participation in EV production. Finally, the disruption of an ESCRT-III associated gene, Pfvps60, led to a significant reduction in the amount of EVs. Altogether, these results confirm ESCRT-III participation in EV production and provide novel information on the P. falciparum protein export mechanisms, which can be used for the development of new therapeutic strategies against malaria, based on the disruption of EV formation and trafficking.

Identification of Three Novel Plasmodium Factors Involved in Ookinete to Oocyst Developmental Transition

Plasmodium falciparum malaria remains a major cause of global morbidity and mortality, mainly in sub-Saharan Africa. The numbers of new malaria cases and deaths have been stable in the last years despite intense efforts for disease elimination, highlighting the need for new approaches to stop disease transmission. Further understanding of the parasite transmission biology could provide a framework for the development of such approaches. We phenotypically and functionally characterized three novel genes, PIMMS01, PIMMS57, and PIMMS22, using targeted disruption of their orthologs in the rodent parasite Plasmodium berghei. PIMMS01 and PIMMS57 are specifically and highly expressed in ookinetes, while PIMMS22 transcription starts already in gametocytes and peaks in sporozoites. All three genes show strong phenotypes associated with the ookinete to oocyst transition, as their disruption leads to very low numbers of oocysts and complete abolishment of transmission. PIMMS22 has a secondary essential function in the oocyst. Our results enrich the molecular understanding of the parasite-vector interactions and identify PIMMS01, PIMMS57, and PIMMS22 as new targets of transmission blocking interventions.

A suitable RNA preparation methodology for whole transcriptome shotgun sequencing harvested from Plasmodium vivax-infected patients

Plasmodium vivax is a world-threatening human malaria parasite, whose biology remains elusive. The unavailability of in vitro culture, and the difficulties in getting a high number of pure parasites makes RNA isolation in quantity and quality a challenge. Here, a methodological outline for RNA-seq from P. vivax isolates with low parasitemia is presented, combining parasite maturation and enrichment with efficient RNA extraction, yielding ~ 100 pg.µL⁻¹ of RNA, suitable for SMART-Seq Ultra-Low Input RNA library and Illumina sequencing. Unbiased coding transcriptome of ~ 4 M reads was achieved for four patient isolates with ~ 51% of transcripts mapped to the P. vivax P01 reference genome, presenting heterogeneous profiles of expression among individual isolates. Amongst the most transcribed genes in all isolates, a parasite-staged mixed repertoire of conserved parasite metabolic, membrane and exported proteins was observed. Still, a quarter of transcribed genes remain functionally uncharacterized. In parallel, a P. falciparum Brazilian isolate was also analyzed and 57% of its transcripts mapped against IT genome. Comparison of transcriptomes of the two species revealed a common trophozoite-staged expression profile, with several homologous genes being expressed. Collectively, these results will positively impact vivax research improving knowledge of P. vivax biology.

Using scRNA-seq to Identify Transcriptional Variation in the Malaria Parasite Ookinete Stage

The crossing of the mosquito midgut epithelium by the malaria parasite motile ookinete form represents the most extreme population bottleneck in the parasite life cycle and is a prime target for transmission blocking strategies. However, we have little understanding of the clonal variation that exists in a population of ookinetes in the vector, partially because the parasites are difficult to access and are found in low numbers. Within a vector, variation may result as a response to specific environmental cues or may exist independent of those cues as a potential bet-hedging strategy. Here we use single-cell RNA-seq to profile transcriptional variation in Plasmodium berghei ookinetes across different vector species, and between and within individual midguts. We then compare our results to low-input transcriptomes from individual Anopheles coluzzii midguts infected with the human malaria parasite Plasmodium falciparum. Although the vast majority of transcriptional changes in ookinetes are driven by development, we have identified candidate genes that may be responding to environmental cues or are clonally variant within a population. Our results illustrate the value of single-cell and low-input technologies in understanding clonal variation of parasite populations.

A systems-level gene regulatory network model for Plasmodium falciparum

Many of the gene regulatory processes of Plasmodium falciparum, the deadliest malaria parasite, remain poorly understood. To develop a comprehensive guide for exploring this organism's gene regulatory network, we generated a systems-level model of P. falciparum gene regulation using a well-validated, machine-learning approach for predicting interactions between transcription regulators and their targets. The resulting network accurately predicts expression levels of transcriptionally coherent gene regulatory programs in independent transcriptomic data sets from parasites collected by different research groups in diverse laboratory and field settings. Thus, our results indicate that our gene regulatory model has predictive power and utility as a hypothesis-generating tool for illuminating clinically relevant gene regulatory mechanisms within P. falciparum. Using the set of regulatory programs we identified, we also investigated correlates of artemisinin resistance based on gene expression coherence. We report that resistance is associated with incoherent expression across many regulatory programs, including those controlling genes associated with erythrocyte-host engagement. These results suggest that parasite populations with reduced artemisinin sensitivity are more transcriptionally heterogenous. This pattern is consistent with a model where the parasite utilizes bet-hedging strategies to diversify the population, rendering a subpopulation more able to navigate drug treatment.

Silent control: microbial plant pathogens evade host immunity without coding sequence changes

Both animals and plants have evolved a robust immune system to surveil and defeat invading pathogenic microbes. Evasion of host immune surveillance is the key for pathogens to initiate successful infection. To evade the host immunity, plant pathogens evolved a variety of strategies such as masking themselves from host immune recognitions, blocking immune signaling transductions, reprogramming immune responses and adapting to immune microenvironmental changes. Gain of new virulence genes, sequence and structural variations enables plant pathogens to evade host immunity through changes in the genetic code. However, recent discoveries demonstrated that variations at the transcriptional, post-transcriptional, post-translational and glycome level enable pathogens to cope with the host immune system without coding sequence changes. The biochemical modification of pathogen associated molecular patterns and silencing of effector genes emerged as potent ways for pathogens to hide from host recognition. Altered processing in mRNA activities provide pathogens with resilience to microenvironment changes. Importantly, these hiding variants are directly or indirectly modulated by catalytic enzymes or enzymatic complexes and cannot be revealed by classical genomics alone. Unveiling these novel host evasion mechanisms in plant pathogens enables us to better understand the nature of plant disease and pinpoints strategies for rational diseases management in global food protection.

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.

Analysis of Plasmodium vivax schizont transcriptomes from field isolates reveals heterogeneity of expression of genes involved in host-parasite interactions

Plasmodium vivax gene regulation remains difficult to study due to the lack of a robust in vitro culture method, low parasite densities in peripheral circulation and asynchronous parasite development. We adapted an RNA-seq protocol “DAFT-seq” to sequence the transcriptome of four P. vivax field isolates that were cultured for a short period ex vivo before using a density gradient for schizont enrichment. Transcription was detected from 78% of the PvP01 reference genome, despite being schizont-enriched samples. This extensive data was used to define thousands of 5′ and 3′ untranslated regions, some of which overlapped with neighbouring transcripts, and to improve the gene models of 352 genes, including identifying 20 novel gene transcripts. This dataset has also significantly increased the known amount of heterogeneity between P. vivax schizont transcriptomes from individual patients. The majority of genes found to be differentially expressed between the isolates lack Plasmodium falciparum homologs and are predicted to be involved in host-parasite interactions, with an enrichment in reticulocyte binding proteins, merozoite surface proteins and exported proteins with unknown function. An improved understanding of the diversity within P. vivax transcriptomes will be essential for the prioritisation of novel vaccine targets.

N6-Adenosine methylation on mRNA is recognized by YTH2 domain protein of human malaria parasite Plasmodium falciparum

Background

Plasmodium falciparum exhibits high translational plasticity during its development in RBCs, yet the regulation at the post-transcriptional level is not well understood. The N⁶-methyl adenosine (m6A) is an important epigenetic modification primarily present on mRNA that controls the levels of transcripts and efficiency of translation in eukaryotes. Recently, the dynamics of m6A on mRNAs at all three developmental stages of P. falciparum in RBCs have been profiled; however, the proteins that regulate the m6A containing mRNAs in the parasites are unknown.

Results

Using sequence analysis, we computationally identified that the P. falciparum genome encodes two putative YTH (YT521-B Homology) domain-containing proteins, which could potentially bind to m6A containing mRNA. We developed a modified methylated RNA immunoprecipitation (MeRIP) assay using PfYTH2 and find that it binds selectively to m6A containing transcripts. The PfYTH2 has a conserved aromatic amino acid cage that forms the methyl-binding pocket. Through site-directed mutagenesis experiments and molecular dynamics simulations, we show that F98 residue is important for m6A binding on mRNA. Fluorescence depolarization assay confirmed that PfYTH2 binds to methylated RNA oligos with high affinity. Further, MeRIP sequencing data revealed that PfYTH2 has more permissive sequence specificity on target m6A containing mRNA than other known eukaryotic YTH proteins. Taken together, here we identify and characterize PfYTH2 as the major protein that could regulate m6A containing transcripts in P. falciparum.

Conclusion

Plasmodium spp. lost the canonical m6A-specific demethylases in their genomes, however, the YTH domain-containing proteins seem to be retained. This study presents a possibility that the YTH proteins are involved in post-transcriptional control in P. falciparum, and might orchestrate the translation of mRNA in various developmental stages of P. falciparum. This is perhaps the first characterization of the methyl-reading function of YTH protein in any parasites.

Conditional expression of PfAP2-G for controlled massive sexual conversion in Plasmodium falciparum

Malaria transmission requires that some asexual parasites convert into sexual forms termed gametocytes. The initial stages of sexual development, including sexually committed schizonts and sexual rings, remain poorly characterized, mainly because they are morphologically identical to their asexual counterparts and only a small subset of parasites undergo sexual development. Here, we describe a system for controlled sexual conversion in the human malaria parasite Plasmodium falciparum, based on conditional expression of the PfAP2-G transcription factor. Using this system, ~90 percent of the parasites converted into sexual forms upon induction, enabling the characterization of committed and early sexual stages without further purification. We characterized sexually committed schizonts and sexual rings at the transcriptomic and phenotypic levels, which revealed down-regulation of genes involved in solute transport upon sexual commitment, among other findings. The new inducible lines will facilitate the study of early sexual stages at additional levels, including multiomic characterization and drug susceptibility assays.

Disruption of the Plasmodium falciparum Life Cycle through Transcriptional Reprogramming by Inhibitors of Jumonji Demethylases

Little is known about the role of the three Jumonji C (JmjC) enzymes in Plasmodium falciparum (Pf). Here, we show that JIB-04 and other established inhibitors of mammalian JmjC histone demethylases kill asexual blood stage parasites and are even more potent at blocking gametocyte development and gamete formation. In late stage parasites, JIB-04 increased levels of trimethylated lysine residues on histones, suggesting the inhibition of P. falciparum Jumonji demethylase activity. These epigenetic defects coincide with deregulation of invasion, cell motor, and sexual development gene programs, including gene targets coregulated by the PfAP2-I transcription factor and chromatin-binding factor, PfBDP1. Mechanistically, we demonstrate that PfJmj3 converts 2-oxoglutarate to succinate in an iron-dependent manner consistent with mammalian Jumonji enzymes, and this catalytic activity is inhibited by JIB-04 and other Jumonji inhibitors. Our pharmacological studies of Jumonji activity in the malaria parasite provide evidence that inhibition of these enzymatic activities is detrimental to the parasite.

The second life of Plasmodium in the mosquito host: gene regulation on the move

Malaria parasites face dynamically changing environments and strong selective constraints within human and mosquito hosts. To survive such hostile and shifting conditions, Plasmodium switches transcriptional programs during development and has evolved mechanisms to adjust its phenotype through heterogeneous patterns of gene expression. In vitro studies on culture-adapted isolates have served to set the link between chromatin structure and functional gene expression. Yet, experimental evidence is limited to certain stages of the parasite in the vertebrate, i.e. blood, while the precise mechanisms underlying the dynamic regulatory landscapes during development and in the adaptation to within-host conditions remain poorly understood. In this review, we discuss available data on transcriptional and epigenetic regulation in Plasmodium mosquito stages in the context of sporogonic development and phenotypic variation, including both bet-hedging and environmentally triggered direct transcriptional responses. With this, we advocate the mosquito offers an in vivo biological model to investigate the regulatory networks, transcription factors and chromatin-modifying enzymes and their modes of interaction with regulatory sequences, which might be responsible for the plasticity of the Plasmodium genome that dictates stage- and cell type-specific blueprints of gene expression.