Comparative Heterochromatin Profiling Reveals Conserved and Unique Epigenome Signatures Linked to Adaptation and Development of Malaria Parasites

Deciphering the Plasmodium falciparum perinuclear var gene expression site

The protozoan parasite Plasmodium falciparum, responsible for the deadliest form of human malaria, employs antigenic variation via monoallelic expression as a key survival strategy. The selective activation of one out of the 60-member var gene family is key to understanding the parasite's ability to cause severe disease and evade the host immune response. var gene activation is initiated by its relocation to a specialized expression site. While the perinuclear expression site (PES) plays a crucial role in enabling the expression of a single allele, the characteristics of this PES remain largely obscure. Recent breakthroughs in genome editing tools and the discovery of regulatory noncoding RNAs have shed light on this intriguing biological feature, offering significant insights into the mechanisms of pathogen virulence.

Regulation of sexual commitment in malaria parasites - a complex affair

Malaria blood stage parasites commit to either one of two distinct cellular fates while developing within erythrocytes of their mammalian host: they either undergo another round of asexual replication or they differentiate into nonreplicative transmissible gametocytes. Depending on the state of infection, either path may support or impair the ultimate goal of human-to-human transmission via the mosquito vector. Malaria parasites therefore evolved strategies to control investments into asexual proliferation versus gametocyte formation. Recent work provided fascinating molecular insight into shared and unique mechanisms underlying the control and environmental modulation of sexual commitment in the two most widely studied malaria parasite species, Plasmodium falciparum and P. berghei. With this review, we aim at placing these findings into a comparative mechanistic context.

Ready for renascence in mosquito: The regulation of gene expression in Plasmodium sexual development

Malaria remains one of the most perilous vector-borne infectious diseases for humans globally. Sexual gametocyte represents the exclusive stage at which malaria parasites are transmitted from the vertebrate to the Anopheles host. The feasible and effective approach to prevent malaria transmission is by addressing the sexual developmental processes, that is, gametocytogenesis and gametogenesis. Thus, this review will comprehensively cover advances in the regulation of gene expression surrounding the transmissible stages, including epigenetic, transcriptional, and post-transcriptional control.

Expression of the MSPDBL2 antigen in a discrete subset of Plasmodium falciparum schizonts is regulated by GDV1 but may not be linked to sexual commitment

The Plasmodium falciparum merozoite surface protein MSPDBL2 is a polymorphic antigen targeted by acquired immune responses, and normally expressed in only a minority of mature schizonts. The potential relationship of MSPDBL2 to sexual commitment is examined, as variable mspdbl2 transcript levels and proportions of MSPDBL2-positive mature schizonts in clinical isolates have previously correlated with levels of many sexual stage parasite gene transcripts, although not with the master regulator ap2-g. It is demonstrated that conditional overexpression of the gametocyte development protein GDV1, which promotes sexual commitment, also substantially increases the proportion of MSPDBL2-positive schizonts in culture. Conversely, truncation of the gdv1 gene is shown to prevent any expression of MSPDBL2. However, across diverse P. falciparum cultured lines, the variable proportions of MSPDBL2 positivity in schizonts do not correlate significantly with variable gametocyte conversion rates, indicating it is not involved in sexual commitment. Confirming this, examining a line with endogenous hemagglutinin-tagged AP2-G showed that the individual schizonts expressing MSPDBL2 are mostly different from those expressing AP2-G. Using a selection-linked integration system, modified P. falciparum lines were engineered to express an intact or disrupted version of MSPDBL2, showing the protein is not required for sexual commitment or early gametocyte development. Asexual parasite multiplication rates were also not affected by expression of either intact or disrupted MSPDBL2 in a majority of schizonts. Occurring alongside sexual commitment, the role of the discrete MSPDBL2-positive schizont subpopulation requires further investigation in natural infections where it is under immune selection.

IMPORTANCE

Malaria parasites in the blood are remarkably variable, able to switch antigenic targets so they may survive within humans who have already developed specific immune responses. This is one of the challenges in developing vaccines against malaria. MSPDBL2 is a target of naturally acquired immunity expressed in minority proportions of schizonts, the end stages of each 2-day replication cycle in red blood cells which contain merozoites prepared to invade new red blood cells. Results show that the proportion of schizonts expressing MSPDBL2 is positively controlled by the expression of the regulatory gametocyte development protein GDV1. It was previously known that expression of GDV1 leads to increased expression of AP2-G which causes parasites to switch to sexual development, so a surprising finding here is that MSPDBL2-positive parasites are mostly distinct from those that express AP2-G. This discrete antigenic subpopulation of mostly asexual parasites is regulated alongside sexually committed parasites, potentially enabling survival under stress conditions.

Reversible host cell surface remodelling limits immune recognition and maximizes transmission of Plasmodium falciparum gametocytes

Reducing malaria transmission has been a major pillar of control programmes and is considered crucial for achieving malaria elimination. Gametocytes, the transmissible forms of the P. falciparum parasite, arise during the blood stage of the parasite and develop through 5 morphologically distinct stages. Immature gametocytes (stage I-IV) sequester and develop in the extravascular niche of the bone marrow and possibly spleen. Only mature stage V gametocytes re-enter peripheral circulation to be taken up by mosquitoes for successful onward transmission. We have recently shown that immature, but not mature gametocytes are targets of host immune responses and identified putative target surface antigens. We hypothesize that these antigens play a role in gametocyte sequestration and contribute to acquired transmission-reducing immunity. Here we demonstrate that surface antigen expression, serum reactivity by human IgG, and opsonic phagocytosis by macrophages all show similar dynamics during gametocyte maturation, i.e., on in immature and off in mature gametocytes. Moreover, the switch in surface reactivity coincides with reversal in phosphatidylserine (PS) surface exposure, a marker for red blood cell age and clearance. PS is exposed on the surface of immature gametocytes (as well as in late asexual stages) but is removed from the surface in later gametocyte stages (IV-V). Using parasite reverse genetics and drug perturbations, we confirm that parasite protein export into the host cell and phospholipid scramblase activity are required for the observed surface modifications in asexual and sexual P. falciparum stages. These findings suggest that the dynamic surface remodelling allows (i) immature gametocyte sequestration in bone marrow and (ii) mature gametocyte release into peripheral circulation and immune evasion, therefore contributing to mature gametocyte survival in vivo and onward transmission to mosquitoes. Importantly, blocking scramblase activity during gametocyte maturation results in efficient clearance of mature gametocytes, revealing a potential path for transmission blocking interventions. Our studies have important implications for our understanding of parasite biology and form a starting point for novel intervention strategies to simultaneously reduce parasite burden and transmission.

Revisiting the determinants of malaria transmission

Malaria parasites have coevolved with humans over thousands of years, mirroring their migration out of Africa. They persist to this day, despite continuous elimination efforts worldwide. These parasites can adapt to changing environments during infection of human and mosquito, and when expanding the geographical range by switching vector species. Recent studies in the human malaria parasite, Plasmodium falciparum, identified determinants governing the plasticity of sexual conversion rates, sex ratio, and vector competence. Here we summarize the latest literature revealing environmental, epigenetic, and genetic determinants of malaria transmission.

The genetic landscape of origins of replication in P. falciparum

Various origin mapping approaches have enabled genome-wide identification of origins of replication (ORI) in model organisms, but only a few studies have focused on divergent organisms. By employing three complementary approaches we provide a high-resolution map of ORIs in Plasmodium falciparum, the deadliest human malaria parasite. We profiled the distribution of origin of recognition complex (ORC) binding sites by ChIP-seq of two PfORC subunits and mapped active ORIs using NFS and SNS-seq. We show that ORIs lack sequence specificity but are not randomly distributed, and group in clusters. Licensing is biased towards regions of higher GC content and associated with G-quadruplex forming sequences (G4FS). While strong transcription likely enhances firing, active origins are depleted from transcription start sites. Instead, most accumulate in transcriptionally active gene bodies. Single molecule analysis of nanopore reads containing multiple initiation events, which could have only come from individual nuclei, showed a relationship between the replication fork pace and the distance to the nearest origin. While some similarities were drawn with the canonic eukaryote model, the distribution of ORIs in P. falciparum is likely shaped by unique genomic features such as extreme AT-richness-a product of evolutionary pressure imposed by the parasitic lifestyle.

DNA-binding protein PfAP2-P regulates parasite pathogenesis during malaria parasite blood stages

Malaria-associated pathogenesis such as parasite invasion, egress, host cell remodelling and antigenic variation requires concerted action by many proteins, but the molecular regulation is poorly understood. Here we have characterized an essential Plasmodium-specific Apicomplexan AP2 transcription factor in Plasmodium falciparum (PfAP2-P; pathogenesis) during the blood-stage development with two peaks of expression. An inducible knockout of gene function showed that PfAP2-P is essential for trophozoite development, and critical for var gene regulation, merozoite development and parasite egress. Chromatin immunoprecipitation sequencing data collected at timepoints matching the two peaks of pfap2-p expression demonstrate PfAP2-P binding to promoters of genes controlling trophozoite development, host cell remodelling, antigenic variation and pathogenicity. Single-cell RNA sequencing and fluorescence-activated cell sorting revealed de-repression of most var genes in Δpfap2-p parasites. Δpfap2-p parasites also overexpress early gametocyte marker genes, indicating a regulatory role in sexual stage conversion. We conclude that PfAP2-P is an essential upstream transcriptional regulator at two distinct stages of the intra-erythrocytic development cycle.

Sequence-dependent heterochromatin formation in the human malaria parasite Plasmodium falciparum

The human malaria parasite Plasmodium falciparum represses transcription of the gene encoding AP2-G, which is the master regulator of germ cell differentiation, via heterochromatin condensation following histone H3 lysine 9 trimethylation (H3K9me3). Although H3K9me3-marked heterochromatin is typically constitutive and its establishment depends on the RNA interference (RNAi) pathway in fission yeast centromeres, malaria parasites lack molecular members essential for RNAi. We developed a strategy to assess heterochromatin establishment on artificial chromosomes introduced into P. falciparum. We show that a particular DNA sequence in the AP2-G promoter is able to induce de novo H3K9me3 nucleosome deposition. In addition, we also found that the AP2-G promoter contains a distinct element required in maintenance of the repression memory. Thus, we speculate that malaria parasites have evolutionarily acquired a sequence-dependent establishment system of non-constitutive, i.e. facultative, H3K9me3-marked heterochromatin.

Plasmodium falciparum Development from Gametocyte to Oocyst: Insight from Functional Studies

Malaria elimination may never succeed without the implementation of transmission-blocking strategies. The transmission of Plasmodium spp. parasites from the human host to the mosquito vector depends on circulating gametocytes in the peripheral blood of the vertebrate host. Once ingested by the mosquito during blood meals, these sexual forms undergo a series of radical morphological and metabolic changes to survive and progress from the gut to the salivary glands, where they will be waiting to be injected into the vertebrate host. The design of effective transmission-blocking strategies requires a thorough understanding of all the mechanisms that drive the development of gametocytes, gametes, sexual reproduction, and subsequent differentiation within the mosquito. The drastic changes in Plasmodium falciparum shape and function throughout its life cycle rely on the tight regulation of stage-specific gene expression. This review outlines the mechanisms involved in Plasmodium falciparum sexual stage development in both the human and mosquito vector, and zygote to oocyst differentiation. Functional studies unravel mechanisms employed by P. falciparum to orchestrate the expression of stage-specific functional products required to succeed in its complex life cycle, thus providing us with potential targets for developing new therapeutics. These mechanisms are based on studies conducted with various Plasmodium species, including predominantly P. falciparum and the rodent malaria parasites P. berghei. However, the great potential of epigenetics, genomics, transcriptomics, proteomics, and functional genetic studies to improve the understanding of malaria as a disease remains partly untapped because of limitations in studies using human malaria parasites and field isolates.

Sexual differentiation in human malaria parasites is regulated by competition between phospholipid metabolism and histone methylation

For Plasmodium falciparum, the most widespread and virulent malaria parasite that infects humans, persistence depends on continuous asexual replication in red blood cells, while transmission to their mosquito vector requires asexual blood-stage parasites to differentiate into non-replicating gametocytes. This decision is controlled by stochastic derepression of a heterochromatin-silenced locus encoding AP2-G, the master transcription factor of sexual differentiation. The frequency of ap2-g derepression was shown to be responsive to extracellular phospholipid precursors but the mechanism linking these metabolites to epigenetic regulation of ap2-g was unknown. Through a combination of molecular genetics, metabolomics and chromatin profiling, we show that this response is mediated by metabolic competition for the methyl donor S-adenosylmethionine between histone methyltransferases and phosphoethanolamine methyltransferase, a critical enzyme in the parasite's pathway for de novo phosphatidylcholine synthesis. When phosphatidylcholine precursors are scarce, increased consumption of SAM for de novo phosphatidylcholine synthesis impairs maintenance of the histone methylation responsible for silencing ap2-g, increasing the frequency of derepression and sexual differentiation. This provides a key mechanistic link that explains how LysoPC and choline availability can alter the chromatin status of the ap2-g locus controlling sexual differentiation.

Epigenetically regulated RNA-binding proteins signify malaria hypnozoite dormancy

Dormancy enables relapsing malaria parasites, such as Plasmodium vivax and cynomolgi, to survive unfavorable conditions. It is enabled by hypnozoites, parasites remaining quiescent inside hepatocytes before reactivating and establishing blood-stage infection. We integrate omics approaches to explore gene-regulatory mechanisms underlying hypnozoite dormancy. Genome-wide profiling of activating and repressing histone marks identifies a few genes that get silenced by heterochromatin during hepatic infection of relapsing parasites. By combining single-cell transcriptomics, chromatin accessibility profiling, and fluorescent in situ RNA hybridization, we show that these genes are expressed in hypnozoites and that their silencing precedes parasite development. Intriguingly, these hypnozoite-specific genes mainly encode proteins with RNA-binding domains. We hence hypothesize that these likely repressive RNA-binding proteins keep hypnozoites in a developmentally competent but dormant state and that heterochromatin-mediated silencing of the corresponding genes aids reactivation. Exploring the regulation and exact function of these proteins hence could provide clues for targeted reactivation and killing of these latent pathogens.

Plasmodium falciparum gametocytes display global chromatin remodelling during sexual differentiation

BACKGROUND

The protozoan malaria parasite Plasmodium falciparum has a complex life cycle during which it needs to differentiate into multiple morphologically distinct life forms. A key process for transmission of the disease is the development of male and female gametocytes in the human blood, yet the mechanisms determining sexual dimorphism in these haploid, genetically identical sexual precursor cells remain largely unknown. To understand the epigenetic program underlying the differentiation of male and female gametocytes, we separated the two sexual forms by flow cytometry and performed RNAseq as well as comprehensive ChIPseq profiling of several histone variants and modifications.

RESULTS

We show that in female gametocytes the chromatin landscape is globally remodelled with respect to genome-wide patterns and combinatorial usage of histone variants and histone modifications. We identified sex specific differences in heterochromatin distribution, implicating exported proteins and ncRNAs in sex determination. Specifically in female gametocytes, the histone variants H2A.Z/H2B.Z were highly enriched in H3K9me3-associated heterochromatin. H3K27ac occupancy correlated with stage-specific gene expression, but in contrast to asexual parasites this was unlinked to H3K4me3 co-occupancy at promoters in female gametocytes.

CONCLUSIONS

Collectively, we defined novel combinatorial chromatin states differentially organising the genome in gametocytes and asexual parasites and unravelled fundamental, sex-specific differences in the epigenetic code. Our chromatin maps represent an important resource for future understanding of the mechanisms driving sexual differentiation in P. falciparum.

PfSET2 Is Involved in Genome Organization of Var Gene Family in Plasmodium falciparum

The three-dimensional (3D) genome structure of human malaria parasite Plasmodium falciparum is highly organized and plays important roles in regulating coordinated expression patterns of specific genes such as virulence genes which are involved in antigenic variation and immune escape. However, the molecular mechanisms that control 3D genome of the parasite remain elusive. Here, by analyzing genome organization of P. falciparum, we identify high-interacting regions (HIRs) with strong chromatin interactions at telomeres and virulence genes loci. Specifically, HIRs are highly enriched with repressive histone marks (H3K36me3 and H3K9me3) and form the transcriptional repressive center. Deletion of PfSET2, which controls H3K36me3 level, results in marked reduction of both intrachromosomal and interchromosomal interactions for HIRs. Importantly, such chromatin reorganization coordinates with dynamic changes in epigenetic feature in HIRs and transcriptional activation of var genes. Additionally, different cluster of var genes based on the pattern of chromatin interactions show distinct transcriptional activation potential after deletion of PfSET2. Our results uncover a fundamental mechanism that the epigenetic factor PfSET2 controls the 3D organization of heterochromatin to regulate the transcription activities of var genes family in P. falciparum. IMPORTANCE PfSET2 has been reported to play key role in silencing var genes in Plasmodium falciparum, while the underlying molecular mechanisms remain unclear. Here, we provide evidence that PfSET2 is essential to maintain 3D genome organization of heterochromatin region to keep var genes in transcription repressive state. These findings can contribute better understanding of the regulation of high-order chromatin structure in P. falciparum.

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.

High-throughput analysis of the transcriptional patterns of sexual genes in malaria

BACKGROUND

Plasmodium falciparum (Pf) is the leading protozoan causing malaria, the most devastating parasitic disease. To ensure transmission, a small subset of Pf parasites differentiate into the sexual forms (gametocytes). Since the abundance of these essential parasitic forms is extremely low within the human host, little is currently known about the molecular regulation of their sexual differentiation, highlighting the need to develop tools to investigate Pf gene expression during this fundamental mechanism.

METHODS

We developed a high-throughput quantitative Reverse-Transcription PCR (RT-qPCR) platform to robustly monitor Pf transcriptional patterns, in particular, systematically profiling the transcriptional pattern of a large panel of gametocyte-related genes (GRG). Initially, we evaluated the technical performance of the systematic RT-qPCR platform to ensure it complies with the accepted quality standards for: (i) RNA extraction, (ii) cDNA synthesis and (iii) evaluation of gene expression through RT-qPCR. We then used this approach to monitor alterations in gene expression of a panel of GRG upon treatment with gametocytogenesis regulators.

RESULTS

We thoroughly elucidated GRG expression profiles under treatment with the antimalarial drug dihydroartemisinin (DHA) or the metabolite choline over the course of a Pf blood cycle (48 h). We demonstrate that both significantly alter the expression pattern of PfAP2-G, the gametocytogenesis master regulator. However, they also markedly modify the developmental rate of the parasites and thus might bias the mRNA expression. Additionally, we screened the effect of the metabolites lactate and kynurenic acid, abundant in severe malaria, as potential regulators of gametocytogenesis.

CONCLUSIONS

Our data demonstrate that the high-throughput RT-qPCR method enables studying the immediate transcriptional response initiating gametocytogenesis of the parasites from a very low volume of malaria-infected RBC samples. The obtained data expand the current knowledge of the initial alterations in mRNA profiles of GRG upon treatment with reported regulators. In addition, using this method emphasizes that asexual parasite stage composition is a crucial element that must be considered when interpreting changes in GRG expression by RT-qPCR, specifically when screening for novel compounds that could regulate Pf sexual differentiation.

Repurposing of Plasmodium falciparum var genes beyond the blood stage

A commonly observed survival strategy in protozoan parasites is the sequential expression of clonally variant-surface antigens to avoid elimination by the host's immune response. In malaria-causing P. falciparum, the immunovariant erythrocyte-membrane protein-1 (PfEMP1) adhesin family, encoded by var genes, is responsible for both antigenic variation and cytoadherence of infected erythrocytes to the microvasculature. Until recently, the biological function of these variant genes was believed to be restricted to intraerythrocytic developmental stages. With the advent of new technologies, var gene expression has been confirmed in transmission and pre-erythrocytic stages. Here, we discuss how repurposing of var gene expression beyond chronic blood-stage infection may be critical for successful transmission.

Genome-wide profiling of histone modifications in Plasmodium falciparum using CUT&RUN

We recently adapted a CUT&RUN protocol for genome-wide profiling of chromatin modifications in the human malaria parasite Plasmodium Using the step-by-step protocol described below, we were able to generate high-quality profiles of multiple histone modifications using only a small fraction of the cells required for ChIP-seq. Using antibodies against two commonly profiled histone modifications, H3K4me3 and H3K9me3, we show here that CUT&RUN profiling is highly reproducible and closely recapitulates previously published ChIP-seq-based abundance profiles of histone marks. Finally, we show that CUT&RUN requires substantially lower sequencing coverage for accurate profiling compared with ChIP-seq.

Inhibitors of ApiAP2 protein DNA binding exhibit multistage activity against Plasmodium parasites

Plasmodium parasites are reliant on the Apicomplexan AP2 (ApiAP2) transcription factor family to regulate gene expression programs. AP2 DNA binding domains have no homologs in the human or mosquito host genomes, making them potential antimalarial drug targets. Using an in-silico screen to dock thousands of small molecules into the crystal structure of the AP2-EXP (Pf3D7_1466400) AP2 domain (PDB:3IGM), we identified putative AP2-EXP interacting compounds. Four compounds were found to block DNA binding by AP2-EXP and at least one additional ApiAP2 protein. Our top ApiAP2 competitor compound perturbs the transcriptome of P. falciparum trophozoites and results in a decrease in abundance of log₂ fold change > 2 for 50% (46/93) of AP2-EXP target genes. Additionally, two ApiAP2 competitor compounds have multi-stage anti-Plasmodium activity against blood and mosquito stage parasites. In summary, we describe a novel set of antimalarial compounds that interact with AP2 DNA binding domains. These compounds may be used for future chemical genetic interrogation of ApiAP2 proteins or serve as starting points for a new class of antimalarial therapeutics.

Plasmodium parasites are the causative agent of malaria, which resulted in over 600,000 deaths in 2021. Due to resistance arising for every antimalarial therapeutic deployed to date, new drug targets and druggable pathways must be explored. To address this concern, we used a molecular docking screen to predict competitors of DNA binding by the parasite specific family of Apicomplexan AP2 (ApiAP2) transcription factor proteins for testing in vitro and in vivo. We find that ApiAP2 competing compounds have antimalarial activity consistent with the disruption of gene regulation. This work will further our understanding of both the biological role and targetability of parasite transcriptional regulation.

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.

Chromatin structure can introduce systematic biases in genome-wide analyses of Plasmodium falciparum

Background: The maintenance, regulation, and dynamics of heterochromatin in the human malaria parasite, Plasmodium falciparum, has drawn increasing attention due to its regulatory role in mutually exclusive virulence gene expression and the silencing of key developmental regulators. The advent of genome-wide analyses such as chromatin-immunoprecipitation followed by sequencing (ChIP-seq) has been instrumental in understanding chromatin composition; however, even in model organisms, ChIP-seq experiments are susceptible to intrinsic experimental biases arising from underlying chromatin structure. Methods: We performed a control ChIP-seq experiment, re-analyzed previously published ChIP-seq datasets and compared different analysis approaches to characterize biases of genome-wide analyses in P. falciparum. Results: We found that heterochromatic regions in input control samples used for ChIP-seq normalization are systematically underrepresented in regard to sequencing coverage across the P. falciparum genome. This underrepresentation, in combination with a non-specific or inefficient immunoprecipitation, can lead to the identification of false enrichment and peaks across these regions. We observed that such biases can also be seen at background levels in specific and efficient ChIP-seq experiments. We further report on how different read mapping approaches can also skew sequencing coverage within highly similar subtelomeric regions and virulence gene families. To ameliorate these issues, we discuss orthogonal methods that can be used to characterize bona fide chromatin-associated proteins. Conclusions: Our results highlight the impact of chromatin structure on genome-wide analyses in the parasite and the need for caution when characterizing chromatin-associated proteins and features.

Distinct Histone Post-translational Modifications during Plasmodium falciparum Gametocyte Development

Histones are the building units of nucleosomes, which constitute chromatin. Histone post-translational modifications (PTMs) play an essential role in epigenetic gene regulation. The Plasmodium falciparum genome encodes canonical and variant histones and a collection of conserved enzymes for histone PTMs and chromatin remodeling. Herein, we profiled the P. falciparum histone PTMs during the development of gametocytes, the obligatory stage for parasite transmission. Mass spectrometric analysis of histones extracted from the early, middle, and late stages of gametocytes identified 457 unique histone peptides with 90 PTMs, of which 50% were novel. The gametocyte histone PTMs display distinct patterns from asexual stages, with many new methylation sites in histones H3 and H3.3 (e.g., K14, K18, and K37). Quantitative analyses revealed a high abundance of acetylation in H3 and H4, mono-methylation of H3/H3.3 K37, and ubiquitination of H3BK112, suggesting that these PTMs play critical roles in gametocytes. Gametocyte histones also showed extensive and unique combinations of PTMs. These data indicate that the parasite harbors distinct transcription regulation mechanisms during gametocyte development and lay the foundation for further characterization of epigenetic regulation in the life cycle of the malaria parasite.

Adapt or Die: Targeting Unique Transmission-Stage Biology for Malaria Elimination

Plasmodium parasites have a complex life cycle that includes development in the human host as well as the Anopheles vector. Successful transmission of the parasite between its host and vector therefore requires the parasite to balance its investments in asexual replication and sexual reproduction, varying the frequency of sexual commitment to persist within the human host and generate future opportunities for transmission. The transmission window is extended further by the ability of stage V gametocytes to circulate in peripheral blood for weeks, whereas immature stage I to IV gametocytes sequester in the bone marrow and spleen until final maturation. Due to the low gametocyte numbers in blood circulation and with the ease of targeting such life cycle bottlenecks, transmission represents an efficient target for therapeutic intervention. The biological process of Plasmodium transmission is a multistage, multifaceted process and the past decade has seen a much deeper understanding of the molecular mechanisms and regulators involved. Clearly, specific and divergent processes are used during transmission compared to asexual proliferation, which both poses challenges but also opportunities for discovery of transmission-blocking antimalarials. This review therefore presents an update of our molecular understanding of gametocyte and gamete biology as well as the status of transmission-blocking activities of current antimalarials and lead development compounds. By defining the biological components associated with transmission, considerations for the development of new transmission-blocking drugs to target such untapped but unique biology is suggested as an important, main driver for transmission-blocking drug discovery.

PfARID Regulates P. falciparum Malaria Parasite Male Gametogenesis and Female Fertility and Is Critical for Parasite Transmission to the Mosquito Vector

Sexual reproduction of Plasmodium falciparum parasites is critical to the spread of malaria in the human population. The factors that regulate gene expression underlying formation of fertilization-competent gametes, however, remain unknown. Here, we report that P. falciparum expresses a protein with an AT-rich interaction domain (ARID) which, in other organisms, is part of chromatin remodeling complexes. P. falciparum ARID (PfARID) localized to the parasite nucleus and is critical for the formation of male gametes and fertility of female gametes. PfARID gene deletion (Pfarid^–) gametocytes showed downregulation of gene expression important for gametogenesis, antigenic variation, and cell signaling and for parasite development in the mosquito. Our study identifies PfARID as a critical nuclear protein involved in regulating the gene expression landscape of mature gametocytes. This establishes fertility and also prepares the parasite for postfertilization events that are essential for infection of the mosquito vector.

Characterising genome architectures using genome decomposition analysis

Genome architecture describes how genes and other features are arranged in genomes. These arrangements reflect the evolutionary pressures on genomes and underlie biological processes such as chromosomal segregation and the regulation of gene expression. We present a new tool called Genome Decomposition Analysis (GDA) that characterises genome architectures and acts as an accessible approach for discovering hidden features of a genome assembly. With the imminent deluge of high-quality genome assemblies from projects such as the Darwin Tree of Life and the Earth BioGenome Project, GDA has been designed to facilitate their exploration and the discovery of novel genome biology. We highlight the effectiveness of our approach in characterising the genome architectures of single-celled eukaryotic parasites from the phylum Apicomplexa and show that it scales well to large genomes.

The Putative Bromodomain Protein PfBDP7 of the Human Malaria Parasite Plasmodium Falciparum Cooperates With PfBDP1 in the Silencing of Variant Surface Antigen Expression

Epigenetic regulation is a critical mechanism in controlling virulence, differentiation, and survival of the human malaria parasite Plasmodium (P.) falciparum. Bromodomain proteins contribute to this process by binding to acetylated lysine residues of histones and thereby targeting the gene regulatory machinery to gene promoters. A protein complex containing the P. falciparum bromodomain proteins (PfBDP) 1 and PfBDP2 (BDP1/BDP2 core complex) was previously shown to play an essential role for the correct transcription of invasion related genes. Here, we performed a functional characterization of a third component of this complex, which we dubbed PfBDP7, because structural modelling predicted a typical bromodomain fold. We confirmed that PfBDP7 is a nuclear protein that interacts with PfBDP1 at invasion gene promoters in mature schizont stage parasites and contributes to their transcription. Although partial depletion of PfBDP7 showed no significant effect on parasite viability, conditional knock down of either PfBDP7 or PfBDP1 resulted in the de-repression of variant surface antigens (VSA), which are important pathogenicity factors. This de-repression was evident both on mRNA and protein level. To understand the underlying mechanism, we mapped the genome wide binding sites of PfBDP7 by ChIPseq and showed that in early schizonts, PfBDP7 and PfBDP1 are commonly enriched in heterochromatic regions across the gene body of all VSA families, including genes coding for PfEMP1, RIFIN, STEVOR, and PfMC-2TM. This suggests that PfBDP7 and PfBDP1 contribute to the silencing of VSAs by associating with heterochromatin. In conclusion, we identified PfBDP7 as a chromatin binding protein that is a constitutive part of the P. falciparum BDP1/BDP2 core complex and established PfBDP1 and PfBDP7 as novel players in the silencing of heterochromatin regulated virulence gene families of the malaria parasite P. falciparum.

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.

Plasmodium falciparum GBP2 Is a Telomere-Associated Protein That Binds to G-Quadruplex DNA and RNA

In the early-diverging protozoan parasite Plasmodium, few telomere-binding proteins have been identified and several are unique. Plasmodium telomeres, like those of most eukaryotes, contain guanine-rich repeats that can form G-quadruplex structures. In model systems, quadruplex-binding drugs can disrupt telomere maintenance and some quadruplex-binding drugs are potent anti-plasmodial agents. Therefore, telomere-interacting and quadruplex-interacting proteins may offer new targets for anti-malarial therapy. Here, we report that P. falciparum GBP2 is such a protein. It was identified via 'Proteomics of Isolated Chromatin fragments', applied here for the first time in Plasmodium. In vitro, PfGBP2 binds specifically to G-rich telomere repeats in quadruplex form and it can also bind to G-rich RNA. In vivo, PfGBP2 partially colocalises with the known telomeric protein HP1 but is also found in the cytoplasm, probably due to its affinity for RNA. Consistently, its interactome includes numerous RNA-associated proteins. PfGBP2 is evidently a multifunctional DNA/RNA-binding factor in Plasmodium.

The transcriptional regulator HDP1 controls expansion of the inner membrane complex during early sexual differentiation of malaria parasites

Transmission of Plasmodium falciparum and other malaria parasites requires their differentiation from asexual blood stages into gametocytes, the non-replicative sexual stage necessary to infect the mosquito vector. This transition involves changes in gene expression and chromatin reorganization that result in the activation and silencing of stage-specific genes. However, the genomes of malaria parasites have been noted for their limited number of transcriptional and chromatin regulators, and the molecular mediators of these changes remain largely unknown. We recently identified homeodomain protein 1 (HDP1) as a DNA-binding protein, first expressed in gametocytes, that enhances the expression of key genes critical for early sexual differentiation. The discovery of HDP1 marks a new class of transcriptional regulator in malaria parasites outside of the better-characterized ApiAP2 family. Here, using molecular biology, biochemistry and microscopy techniques, we show that HDP1 is essential for gametocyte maturation, facilitating the necessary upregulation of inner membrane complex components during early gametocytogenesis that gives P. falciparum gametocytes their characteristic shape.

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.

OUP accepted manuscript

Heterochromatin-associated gene silencing controls multiple physiological processes in malaria parasites, however, little is known concerning the regulatory network and cis-acting sequences involved in the organization of heterochromatin and how they modulate heterochromatic gene expression. Based on systematic profiling of genome-wide occupancy of eighteen Apicomplexan AP2 transcription factors by ChIP-seq analysis, we identify and characterize eight heterochromatin-associated factors (PfAP2-HFs), which exhibit preferential enrichment within heterochromatic regions but with differential coverage profiles. Although these ApiAP2s target euchromatic gene loci via specific DNA motifs, they are likely integral components of heterochromatin independent of DNA motif recognition. Systematic knockout screenings of ApiAP2 factors coupled with RNA-seq transcriptomic profiling revealed three activators and three repressors of heterochromatic gene expression including four PfAP2-HFs. Notably, expression of virulence genes is either completely silenced or significantly reduced upon the depletion of PfAP2-HC. Integrated multi-omics analyses reveal autoregulation and feed-forward loops to be common features of the ApiAP2 regulatory network, in addition to the occurrence of dynamic interplay between local chromatin structure and ApiAP2s in transcriptional control. Collectively, this study provides a valuable resource describing the genome-wide landscape of the ApiAP2 family and insights into functional divergence and cooperation within this family during the blood-stage development of malaria parasites.

"Reading" a new chapter in protozoan parasite transcriptional regulation

Protozoan parasites continue to cause a significant health and economic burden worldwide. As infectious organisms, they pose unique and difficult challenges due to a level of conservation of critical eukaryotic cellular pathways with their hosts. Gene regulation has been pinpointed as an essential pathway with enough divergence to warrant investigation into therapeutically targeting. Examination of human parasites such as Plasmodium falciparum, Toxoplasma gondii, and kinetoplastids have revealed that epigenetic mechanisms play a key role in their gene regulation. The enzymes involved in adding and removing epigenetic posttranslational modifications (PTMs) have historically been the focus of study. However, the reader proteins that recognize and bind PTMs, initiating recruitment of chromatin-modifying and transcription complexes, are now being realized for their critical role in regulation and their potential as drug targets. In this review, we highlight the current knowledge on epigenetic reader proteins in model parasitic protozoa, focusing on the histone acyl- and methyl-reading domains. With this knowledge base, we compare differences between medically relevant parasites, discuss conceivable functions of these understudied proteins, indicate gaps in knowledge, and provide current progress in drug development.

A unique GCN5 histone acetyltransferase complex controls erythrocyte invasion and virulence in the malaria parasite Plasmodium falciparum

The histone acetyltransferase GCN5-associated SAGA complex is evolutionarily conserved from yeast to human and functions as a general transcription co-activator in global gene regulation. In this study, we identified a divergent GCN5 complex in Plasmodium falciparum, which contains two plant homeodomain (PHD) proteins (PfPHD1 and PfPHD2) and a plant apetela2 (AP2)-domain transcription factor (PfAP2-LT). To dissect the functions of the PfGCN5 complex, we generated parasite lines with either the bromodomain in PfGCN5 or the PHD domain in PfPHD1 deleted. The two deletion mutants closely phenocopied each other, exhibiting significantly reduced merozoite invasion of erythrocytes and elevated sexual conversion. These domain deletions caused dramatic decreases not only in histone H3K9 acetylation but also in H3K4 trimethylation, indicating synergistic crosstalk between the two euchromatin marks. Domain deletion in either PfGCN5 or PfPHD1 profoundly disturbed the global transcription pattern, causing altered expression of more than 60% of the genes. At the schizont stage, these domain deletions were linked to specific down-regulation of merozoite genes involved in erythrocyte invasion, many of which contain the AP2-LT binding motif and are also regulated by AP2-I and BDP1, suggesting targeted recruitment of the PfGCN5 complex to the invasion genes by these specific factors. Conversely, at the ring stage, PfGCN5 or PfPHD1 domain deletions disrupted the mutually exclusive expression pattern of the entire var gene family, which encodes the virulent factor PfEMP1. Correlation analysis between the chromatin state and alteration of gene expression demonstrated that up- and down-regulated genes in these mutants are highly correlated with the silent and active chromatin states in the wild-type parasite, respectively. Collectively, the PfGCN5 complex represents a novel HAT complex with a unique subunit composition including an AP2 transcription factor, which signifies a new paradigm for targeting the co-activator complex to regulate general and parasite-specific cellular processes in this low-branching parasitic protist.

Epigenetic regulation of gene expression plays essential roles in orchestrating the general and parasite-specific cellular pathways in the malaria parasite Plasmodium falciparum. To better understand the epigenetic mechanisms in this parasite, we characterized the histone acetyltransferase GCN5-mediated transcription regulation during intraerythrocytic development of the parasite. Using tandem affinity purification and proteomic characterization, we identified that the PfGCN5-associated complex contains nine core components, including two PHD domain proteins (PfPHD1 and PfPHD2) and an AP2-domain transcription factor, which is divergent from the canonical GCN5 complexes evolutionarily conserved from yeast to human. To understand the functions of the PfGCN5 complex, we performed domain deletions in two subunits of this complex, PfGCN5 and PfPHD1. We found that the two deletion mutants displayed very similar growth phenotypes, including significantly reduced merozoite invasion rates and elevated sexual conversion. These two mutants were associated with dramatic decreases in histone H3K9 acetylation and H3K4 trimethylation, which led to global changes in chromatin states and gene expression. Consistent with the phenotypes, genes significantly affected by the PfGCN5 and PfPHD1 gene disruption include those participating in parasite-specific pathways such as invasion, virulence, and sexual development. In conclusion, this study presents a new model of the PfGCN5 complex for targeting the co-activator complex to regulate general and parasite-specific cellular processes in this low-branching parasitic protist.

CRISPR/Cas9-engineered inducible gametocyte producer lines as a valuable tool for Plasmodium falciparum malaria transmission research

The malaria parasite Plasmodium falciparum replicates inside erythrocytes in the blood of infected humans. During each replication cycle, a small proportion of parasites commits to sexual development and differentiates into gametocytes, which are essential for parasite transmission via the mosquito vector. Detailed molecular investigation of gametocyte biology and transmission has been hampered by difficulties in generating large numbers of these highly specialised cells. Here, we engineer P. falciparum NF54 inducible gametocyte producer (iGP) lines for the routine mass production of synchronous gametocytes via conditional overexpression of the sexual commitment factor GDV1. NF54/iGP lines consistently achieve sexual commitment rates of 75% and produce viable gametocytes that are transmissible by mosquitoes. We also demonstrate that further genetic engineering of NF54/iGP parasites is a valuable tool for the targeted exploration of gametocyte biology. In summary, we believe the iGP approach developed here will greatly expedite basic and applied malaria transmission stage research.

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.

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.

The Transcription Factor PfAP2-O Influences Virulence Gene Transcription and Sexual Development in Plasmodium falciparum

The human malaria parasite Plasmodium falciparum expresses variant PfEMP1 proteins on the infected erythrocyte, which function as ligands for endothelial receptors in capillary vessels, leading to erythrocyte sequestration and severe malaria. The factors that orchestrate the mono-allelic expression of the 45–90 PfEMP1-encoding var genes within each parasite genome are still not fully identified. Here, we show that the transcription factor PfAP2-O influences the transcription of var genes. The temporary knockdown of PfAP2-O leads to a complete loss of var transcriptional memory and a decrease in cytoadherence in CD36 adherent parasites. AP2-O-knocked-down parasites exhibited also significant reductions in transmission through Anopheles mosquitoes. We propose that PfAP2-O is, beside its role in transmission stages, also one of the virulence gene transcriptional regulators and may therefore be exploited as an important target to disrupt severe malaria and block parasite transmission.

The Architectural Factor HMGB1 Is Involved in Genome Organization in the Human Malaria Parasite Plasmodium falciparum

The three-dimensional (3D) genome organization plays a critical role in the regulation of gene expression in eukaryotic organisms. In the unicellular malaria parasite Plasmodium falciparum, the high-order chromosome organization has emerged as an important epigenetic pathway mediating gene expression, particularly for virulence genes, but the related architectural factors and underlying mechanism remain elusive. Herein, we have identified the high-mobility-group protein HMGB1 as a critical architectural factor for maintenance of genome organization in P. falciparum Genome-wide occupancy analysis (chromatin immunoprecipitation sequencing [ChIP-seq]) shows that the HMGB1 protein is recruited mainly to centromeric regions likely via a DNA-binding-independent pathway. Chromosome conformation capture coupled with next-generation sequencing (Hi-C-seq) and 3D modeling analysis show that the loss of HMGB1 disrupts the integrity of centromere/telomere-based chromosome organization accompanied with diminished interaction frequency among centromere clusters. This triggers local chromatin alteration and dysregulated gene expression. Notably, the entire repertoire of the primary virulence genes (var) was completely silenced in the absence of P. falciparum HMGB1 (PfHMGB1). Furthermore, the disrupted nuclear organization was reconstituted by complementation of HMGB1, thereby rescuing the mutually exclusive expression of the var gene family. Collectively, these data demonstrate that the architectural factor HMGB1 is associated with gene expression via mediating the high-order structure of genome organization. This finding not only contributes better understanding of the epigenetic regulation of gene expression but may also provide novel targets for antimalarial strategies.IMPORTANCE Malaria remains a major public health and economic burden currently. The mutually exclusive expression of the virulence genes is associated with the pathogenesis and immune evasion of human malaria parasites in the host. The nuclear architecture provides a well-organized environment for differential gene expression in the nucleus, but the underlying mechanism remains largely unknown. In this study, we have identified the highly conserved high-mobility-group protein HMGB1 as a key architecture regulator involved in virulence gene expression by establishing high-order genome organization in the nucleus of P. falciparum Mechanistic investigation revealed that the specific interaction of HMGB1 and centromeres constructed the precisely organized nuclear architecture, which coordinated with local chromatin structure to control the singular expression of virulence genes. Hence, this protein appears to be a critical architectural regulator for the pathogenesis of malaria infection and may be a new target for the development of an intervention strategy against malaria.

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

The Modular Circuitry of Apicomplexan Cell Division Plasticity

The close-knit group of apicomplexan parasites displays a wide variety of cell division modes, which differ between parasites as well as between different life stages within a single parasite species. The beginning and endpoint of the asexual replication cycles is a 'zoite' harboring the defining apical organelles required for host cell invasion. However, the number of zoites produced per division round varies dramatically and can unfold in several different ways. This plasticity of the cell division cycle originates from a combination of hard-wired developmental programs modulated by environmental triggers. Although the environmental triggers and sensors differ between species and developmental stages, widely conserved secondary messengers mediate the signal transduction pathways. These environmental and genetic input integrate in division-mode specific chromosome organization and chromatin modifications that set the stage for each division mode. Cell cycle progression is conveyed by a smorgasbord of positively and negatively acting transcription factors, often acting in concert with epigenetic reader complexes, that can vary dramatically between species as well as division modes. A unique set of cell cycle regulators with spatially distinct localization patterns insert discrete check points which permit individual control and can uncouple general cell cycle progression from nuclear amplification. Clusters of expressed genes are grouped into four functional modules seen in all division modes: 1. mother cytoskeleton disassembly; 2. DNA replication and segregation (D&S); 3. karyokinesis; 4. zoite assembly. A plug-and-play strategy results in the variety of extant division modes. The timing of mother cytoskeleton disassembly is hard-wired at the species level for asexual division modes: it is either the first step, or it is the last step. In the former scenario zoite assembly occurs at the plasma membrane (external budding), and in the latter scenario zoites are assembled in the cytoplasm (internal budding). The number of times each other module is repeated can vary regardless of this first decision, and defines the modes of cell division: schizogony, binary fission, endodyogeny, endopolygeny.

H3K36 methylation reprograms gene expression to drive early gametocyte development in Plasmodium falciparum

Background

The Plasmodium sexual gametocyte stages are the only transmissible form of the malaria parasite and are thus responsible for the continued transmission of the disease. Gametocytes undergo extensive functional and morphological changes from commitment to maturity, directed by an equally extensive control program. However, the processes that drive the differentiation and development of the gametocyte post-commitment, remain largely unexplored. A previous study reported enrichment of H3K36 di- and tri-methylated (H3K36me2&3) histones in early-stage gametocytes. Using chromatin immunoprecipitation followed by high-throughput sequencing, we identify a stage-specific association between these repressive histone modifications and transcriptional reprogramming that define a stage II gametocyte transition point.

Results

Here, we show that H3K36me2 and H3K36me3 from stage II gametocytes are associated with repression of genes involved in asexual proliferation and sexual commitment, indicating that H3K36me2&3-mediated repression of such genes is essential to the transition from early gametocyte differentiation to intermediate development. Importantly, we show that the gene encoding the transcription factor AP2-G as commitment master regulator is enriched with H3K36me2&3 and actively repressed in stage II gametocytes, providing the first evidence of ap2-g gene repression in post-commitment gametocytes. Lastly, we associate the enhanced potency of the pan-selective Jumonji inhibitor JIB-04 in gametocytes with the inhibition of histone demethylation including H3K36me2&3 and a disruption of normal transcriptional programs.

Conclusions

Taken together, our results provide the first description of an association between global gene expression reprogramming and histone post-translational modifications during P. falciparum early sexual development. The stage II gametocyte-specific abundance of H3K36me2&3 manifests predominantly as an independent regulatory mechanism targeted towards genes that are repressed post-commitment. H3K36me2&3-associated repression of genes is therefore involved in key transcriptional shifts that accompany the transition from early gametocyte differentiation to intermediate development.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13072-021-00393-9.

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.

Investigation of Heterochromatin Protein 1 Function in the Malaria Parasite Plasmodium falciparum Using a Conditional Domain Deletion and Swapping Approach

The human malaria parasite Plasmodium falciparum encodes a single ortholog of heterochromatin protein 1 (PfHP1) that plays a crucial role in the epigenetic regulation of various survival-related processes. PfHP1 is essential for parasite proliferation and the heritable silencing of genes linked to antigenic variation, host cell invasion, and sexual conversion. Here, we employed CRISPR/Cas9-mediated genome editing combined with the DiCre/loxP system to investigate how the PfHP1 chromodomain (CD), hinge domain, and chromoshadow domain (CSD) contribute to overall PfHP1 function. We show that the 76 C-terminal residues are responsible for targeting PfHP1 to the nucleus. Furthermore, we reveal that each of the three functional domains of PfHP1 are required for heterochromatin formation, gene silencing, and mitotic parasite proliferation. Finally, we discovered that the hinge domain and CSD of HP1 are functionally conserved between P. falciparum and P. berghei, a related malaria parasite infecting rodents. In summary, our study provides new insights into PfHP1 function and offers a tool for further studies on epigenetic regulation and life cycle decision in malaria parasites.IMPORTANCE Malaria is caused by unicellular Plasmodium species parasites that repeatedly invade and replicate inside red blood cells. Some blood-stage parasites exit the cell cycle and differentiate into gametocytes that are essential for malaria transmission via the mosquito vector. Epigenetic control mechanisms allow the parasites to alter the expression of surface antigens and to balance the switch between parasite multiplication and gametocyte production. These processes are crucial to establish chronic infection and optimize parasite transmission. Here, we performed a mutational analysis of heterochromatin protein 1 (HP1) in P. falciparum We demonstrate that all three domains of this protein are indispensable for the proper function of HP1 in parasite multiplication, heterochromatin formation, and gene silencing. Moreover, expression of chimeric proteins revealed the functional conservation of HP1 proteins between different Plasmodium species. These results provide new insight into the function and evolution of HP1 as an essential epigenetic regulator of parasite survival.

From Genes to Transcripts, a Tightly Regulated Journey in Plasmodium

Over the past decade, we have witnessed significant progresses in understanding gene regulation in Apicomplexa including the human malaria parasite, Plasmodium falciparum. This parasite possesses the ability to convert in multiple stages in various hosts, cell types, and environments. Recent findings indicate that P. falciparum is talented at using efficient and complementary molecular mechanisms to ensure a tight control of gene expression at each stage of its life cycle. Here, we review the current understanding on the contribution of the epigenome, atypical transcription factors, and chromatin organization to regulate stage conversion in P. falciparum. The adjustment of these regulatory mechanisms occurring during the progression of the life cycle will be extensively discussed.

Dynamic Chromatin Structure and Epigenetics Control the Fate of Malaria Parasites

Multiple hosts and various life cycle stages prompt the human malaria parasite, Plasmodium falciparum, to acquire sophisticated molecular mechanisms to ensure its survival, spread, and transmission to its next host. To face these environmental challenges, increasing evidence suggests that the parasite has developed complex and complementary layers of regulatory mechanisms controlling gene expression. Here, we discuss the recent developments in the discovery of molecular components that contribute to cell replication and differentiation and highlight the major contributions of epigenetics, transcription factors, and nuclear architecture in controlling gene regulation and life cycle progression in Plasmodium spp.

Exploring the virulence gene interactome with CRISPR/dCas9 in the human malaria parasite

Mutually exclusive expression of the var multigene family is key to immune evasion and pathogenesis in Plasmodium falciparum, but few factors have been shown to play a direct role. We adapted a CRISPR-based proteomics approach to identify novel factors associated with var genes in their natural chromatin context. Catalytically inactive Cas9 ("dCas9") was targeted to var gene regulatory elements, immunoprecipitated, and analyzed with mass spectrometry. Known and novel factors were enriched including structural proteins, DNA helicases, and chromatin remodelers. Functional characterization of PfISWI, an evolutionarily divergent putative chromatin remodeler enriched at the var gene promoter, revealed a role in transcriptional activation. Proteomics of PfISWI identified several proteins enriched at the var gene promoter such as acetyl-CoA synthetase, a putative MORC protein, and an ApiAP2 transcription factor. These findings validate the CRISPR/dCas9 proteomics method and define a new var gene-associated chromatin complex. This study establishes a tool for targeted chromatin purification of unaltered genomic loci and identifies novel chromatin-associated factors potentially involved in transcriptional control and/or chromatin organization of virulence genes in the human malaria parasite.

Activity of Epigenetic Inhibitors against Plasmodium falciparum Asexual and Sexual Blood Stages

Earlier genetic and inhibitor studies showed that epigenetic regulation of gene expression is critical for malaria parasite survival in multiple life stages and a promising target for new antimalarials. We therefore evaluated the activity of 350 diverse epigenetic inhibitors against multiple stages of Plasmodium falciparum We observed ≥90% inhibition at 10 μM for 28% of compounds against asexual blood stages and early gametocytes, of which a third retained ≥90% inhibition at 1 μM.

Actin-related protein Arp4 regulates euchromatic gene expression and development through H2A.Z deposition in blood-stage Plasmodium falciparum

BACKGROUND

Malaria caused by Plasmodium spp. is still a major threat to public health globally. The various approaches to developing new antimalarial agents rely on the understanding of the complex regulatory mechanisms of dynamic gene expression in the life-cycle of these malaria parasites. The nuclear members of the evolutionarily conserved actin-related protein nuclear (ARP) superfamily are the major components of nucleosome remodelling complexes. In the human malaria parasite Plasmodium falciparum, bioinformatics analysis has predicted three ARP orthologues: PfArp1, PfArp4 and PfArp6. However, little is known about the biological functions of putative PfArp4. In this study, we aimed to investigate the function and the underlying mechanisms of PfArp4 gene regulation.

METHODS

A conditional gene knockdown approach was adopted by incorporating the glucosamine-inducible glmS ribozyme sequence into the 3' UTR of the PfArp4 and PfArp6 genes. The transgenic parasites PfArp4-Ty1-Ribo, PfArp6-Ty1-Ribo and pL6-PfArp4-Ty1::PfArp6-HA were generated by the CRISPR-Cas9 technique. The knockdown effect in the transgenic parasite was measured by growth curve assay and western blot (WB) analysis. The direct interaction between PfArp4 and PfArp6 was validated by co-IFA and co-IP assays. The euchromatic gene expression mediated through H2A.Z (histone H2A variant) deposition and H3K9ac modification at promoters and regulated by PfArp4, was determined by RNA-seq and ChIP-seq.

RESULTS

The inducible knockdown of PfArp4 inhibited blood-stage development of P. falciparum. PfArp4 and PfArp6 were colocalized in the nucleus of P. falciparum parasites. PfArp4 gene knockdown altered the global transcriptome. PfArp4 protein colocalized with the histone variant H2A.Z and euchromatic marker H3K9ac in intergenic regions. The inducible downregulation of PfArp4 resulted in the depletion of H2A.Z and lower H3K9ac levels at the upstream regions of eukaryotic genes, thereby repressing the transcriptional abundance of H2A.Z-dependent genes.

CONCLUSIONS

Our findings suggest that PfArp4 regulates the cell cycle by controlling H2A.Z deposition and affecting centromere function, contributing to the understanding the complex epigenetic regulation of gene expression and the development of P. falciparum.