H2A.Z demarcates intergenic regions of the plasmodium falciparum epigenome that are dynamically marked by H3K9ac and H3K4me3

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.

DNA sequence and chromatin differentiate sequence-specific transcription factor binding in the human malaria parasite Plasmodium falciparum

Development of the malaria parasite, Plasmodium falciparum, is regulated by a limited number of sequence-specific transcription factors (TFs). However, the mechanisms by which these TFs recognize genome-wide binding sites is largely unknown. To address TF specificity, we investigated the binding of two TF subsets that either bind CACACA or GTGCAC DNA sequence motifs and further characterized two additional ApiAP2 TFs, PfAP2-G and PfAP2-EXP, which bind unique DNA motifs (GTAC and TGCATGCA). We also interrogated the impact of DNA sequence and chromatin context on P. falciparum TF binding by integrating high-throughput in vitro and in vivo binding assays, DNA shape predictions, epigenetic post-translational modifications, and chromatin accessibility. We found that DNA sequence context minimally impacts binding site selection for paralogous CACACA-binding TFs, while chromatin accessibility, epigenetic patterns, co-factor recruitment, and dimerization correlate with differential binding. In contrast, GTGCAC-binding TFs prefer different DNA sequence context in addition to chromatin dynamics. Finally, we determined that TFs that preferentially bind divergent DNA motifs may bind overlapping genomic regions due to low-affinity binding to other sequence motifs. Our results demonstrate that TF binding site selection relies on a combination of DNA sequence and chromatin features, thereby contributing to the complexity of P. falciparum gene regulatory mechanisms.

Understanding the complex formation of falstatin; an endogenous macromolecular inhibitor of falcipains

Proteolytic activity constitutes a fundamental process essential for the survival of the malaria parasite and is thus highly regulated. Falstatin, a protease inhibitor of Plasmodium falciparum, tightly regulates the activity of cysteine hemoglobinases, falcipain-2 and 3 (FP2, FP3), by inhibiting FP2 through a single surface exposed loop. However, the multimeric nature of falstatin and its interaction with FP2 remained unexplored. Here we report that the N-terminal falstatin region is highly disordered, and needs chaperone activity (heat-shock protein 70, HSP70) for its folding. Protein-protein interaction assays showed a significant interaction between falstatin and HSP70. Further, characterization of the falstatin multimer through a series of biophysical techniques identified the formation of a falstatin decamer, which was extremely thermostable. Computational analysis of the falstatin decamer showed the presence of five falstatin dimers, with each dimer aligned in a head-to-tail orientation. Further, the falstatin C-terminal region was revealed to be primarily involved in the oligomerization process. Stoichiometric analysis of the FP2-falstatin multimer showed the formation of a heterooligomeric complex in a 1:1 ratio, with the participation of ten subunits of each protein. Taken together, our results report a novel protease-inhibitor complex and strengthens our understanding of the regulatory mechanisms of major plasmodium hemoglobinases.

Incorporation of the histone variant H2A.Z counteracts gene silencing mediated by H3K27 trimethylation in Fusarium fujikuroi

BACKGROUND

Fusarium fujikuroi is a pathogen of rice causing diverse disease symptoms such as 'bakanae' or stunting, most likely due to the production of various natural products (NPs) during infection. Fusaria have the genetic potential to synthesize a plethora of these compounds with often diverse bioactivity. The capability to synthesize NPs exceeds the number of those being produced by far, implying a gene regulatory network decisive to induce production. One such regulatory layer is the chromatin structure and chromatin-based modifications associated with it. One prominent example is the exchange of histones against histone variants such as the H2A variant H2A.Z. Though H2A.Z already is well studied in several model organisms, its regulatory functions are not well understood. Here, we used F. fujikuroi as a model to explore the role of the prominent histone variant FfH2A.Z in gene expression within euchromatin and facultative heterochromatin.

RESULTS

Through the combination of diverse '-omics' methods, we show the global distribution of FfH2A.Z and analyze putative crosstalks between the histone variant and two prominent histone marks, i.e., H3K4me3 and H3K27me3, important for active gene transcription and silencing, respectively. We demonstrate that, if FfH2A.Z is positioned at the + 1-nucleosome, it poises chromatin for gene transcription, also within facultative heterochromatin. Lastly, functional characterization of FfH2A.Z overexpression and depletion mutants revealed that FfH2A.Z is important for wild type-like fungal development and secondary metabolism.

CONCLUSION

In this study, we show that the histone variant FfH2A.Z is a mark of positive gene transcription and acts independently of the chromatin state most likely through the stabilization of the + 1-nucleosome. Furthermore, we demonstrate that FfH2A.Z depletion does not influence the establishment of both H3K27me3 and H3K4me3, thus indicating no crosstalk between FfH2A.Z and both histone marks. These results highlight the manifold functions of the histone variant FfH2A.Z in the phytopathogen F. fujikuroi, which are distinct regarding gene transcription and crosstalk with the two prominent histone marks H3K27me3 and H3K4me3, as proposed for other model organisms.

Histone globular domain epigenetic modifications: The regulators of chromatin dynamics in malaria parasite

Plasmodium species adapt a complex lifecycle with multiple phenotypes to survive inside various cell types of humans and mosquitoes. Stage-specific gene expression in the developmental stages of parasites is tightly controlled in Plasmodium species; however, the underlying mechanisms have yet to be explored. Genome organization and gene expression for each stage of the malaria parasite need to be better characterized. Recent studies indicated that epigenetic modifications of histone proteins play a vital role in chromatin plasticity. Like other eukaryotes, Plasmodium species N-terminal tail modifications form a distinct "histone code," which creates the docking sites for histone reader proteins, including gene activator/repressor complexes, to regulate gene expression. The emerging research findings shed light on various unconventional epigenetic changes in histone proteins' core/globular domain regions, which might contribute to the chromatin organization in different developmental stages of the malaria parasite. The malaria parasite lost many transcription factors during evolution, and it is proposed that the nature of local chromatin structure essentially regulates the stage-specific gene expression. This review highlights recent discoveries of unconventional histone globular domain epigenetic modifications and their functions in regulating chromatin structure dynamics in various developmental stages of malaria parasites.

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.

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.

The P. falciparum alternative histones Pf H2A.Z and Pf H2B.Z are dynamically acetylated and antagonized by PfSir2 histone deacetylases at heterochromatin boundaries

IMPORTANCE

The malaria parasite Plasmodium falciparum relies on variant expression of members of multi-gene families as a strategy for environmental adaptation to promote parasite survival and pathogenesis. These genes are located in transcriptionally silenced DNA regions. A limited number of these genes escape gene silencing, and switching between them confers variant fitness on parasite progeny. Here, we show that PfSir2 histone deacetylases antagonize DNA-interacting acetylated alternative histones at the boundaries between active and silent DNA. This finding implicates acetylated alternative histones in the mechanism regulating P. falciparum variant gene silencing and thus malaria pathogenesis. This work also revealed that acetylation of alternative histones at promoters is dynamically associated with promoter activity across the genome, implicating acetylation of alternative histones in gene regulation genome wide. Understanding mechanisms of gene regulation in P. falciparum may aid in the development of new therapeutic strategies for malaria, which killed 619,000 people in 2021.

Plasmodium falciparum utilizes pyrophosphate to fuel an essential proton pump in the ring stage and the transition to trophozoite stage

During asexual growth and replication cycles inside red blood cells, the malaria parasite Plasmodium falciparum primarily relies on glycolysis for energy supply, as its single mitochondrion performs little or no oxidative phosphorylation. Post merozoite invasion of a host red blood cell, the ring stage lasts approximately 20 hours and was traditionally thought to be metabolically quiescent. However, recent studies have shown that the ring stage is active in several energy-costly processes, including gene transcription, protein translation, protein export, and movement inside the host cell. It has remained unclear whether a low glycolytic flux alone can meet the energy demand of the ring stage over a long period post invasion. Here, we demonstrate that the metabolic by-product pyrophosphate (PPi) is a critical energy source for the development of the ring stage and its transition to the trophozoite stage. During early phases of the asexual development, the parasite utilizes Plasmodium falciparum vacuolar pyrophosphatase 1 (PfVP1), an ancient pyrophosphate-driven proton pump, to export protons across the parasite plasma membrane. Conditional deletion of PfVP1 leads to a delayed ring stage that lasts nearly 48 hours and a complete blockage of the ring-to-trophozoite transition before the onset of parasite death. This developmental arrest can be partially rescued by an orthologous vacuolar pyrophosphatase from Arabidopsis thaliana, but not by the soluble pyrophosphatase from Saccharomyces cerevisiae, which lacks proton pumping activities. Since proton-pumping pyrophosphatases have been evolutionarily lost in human hosts, the essentiality of PfVP1 suggests its potential as an antimalarial drug target. A drug target of the ring stage is highly desired, as current antimalarials have limited efficacy against this stage.

Comparative spatial proteomics of Plasmodium-infected erythrocytes

Plasmodium parasites contribute to one of the highest global infectious disease burdens. To achieve this success, the parasite has evolved a range of specialized subcellular compartments to extensively remodel the host cell for its survival. The information to fully understand these compartments is likely hidden in the so far poorly characterized Plasmodium species spatial proteome. To address this question, we determined the steady-state subcellular location of more than 12,000 parasite proteins across five different species by extensive subcellular fractionation of erythrocytes infected by Plasmodium falciparum, Plasmodium knowlesi, Plasmodium yoelii, Plasmodium berghei, and Plasmodium chabaudi. This comparison of the pan-species spatial proteomes and their expression patterns indicates increasing species-specific proteins associated with the more external compartments, supporting host adaptations and post-transcriptional regulation. The spatial proteome offers comprehensive insight into the different human, simian, and rodent Plasmodium species, establishing a powerful resource for understanding species-specific host adaptation processes in the parasite.

Genome-wide binding sites of Plasmodium falciparum mini chromosome maintenance protein MCM6 show new insights into parasite DNA replication

Multiple rounds of DNA replication take place in various stages of the life cycle in the human malaria parasite Plasmodium falciparum. Previous bioinformatics analysis has shown the presence of putative Autonomously Replicating Sequence (ARS) like sequences in the Plasmodium genome. However, the actual sites and frequency of replication origins in the P. falciparum genome based on experimental data still remain elusive. Minichromosome maintenance (MCM) proteins are recruited by the Origin recognition complex (ORC) to the origins of replication in eukaryotes including P. falciparum. We used PfMCM6 for chromatin immunoprecipitation followed by sequencing (ChIP-seq) in the quest for identification of putative replication origins in the parasite. PfMCM6 DNA binding sites annotation revealed high enrichment at exon regions. This is contrary to higher eukaryotes that show an inclination of origin sites towards transcriptional start sites. ChIP-seq results were further validated by ChIP-qPCR results as well as nascent strand abundance assay at the selected PfMCM6 enriched sites that also showed preferential binding of PfORC1 suggesting potential of these sites as origin sites. Further, PfMCM6 ChIP-seq data showed a positive correlation with previously published histone H4K8Ac genome-wide binding sites but not with H3K9Ac sites suggesting epigenetic control of replication initiation sites in the parasites. Overall, our data show the genome-wide distribution of PfMCM6 binding sites with their potential as replication origins in this deadly human pathogen that not only broadens our knowledge of parasite DNA replication and its unique biology, it may help to find new avenues for intervention processes.

Epigenetic Regulation and Chromatin Remodeling in Malaria Parasites

Plasmodium falciparum, the human malaria parasite, infects two hosts and various cell types, inducing distinct morphological and physiological changes in the parasite in response to different environmental conditions. These variations required the parasite to adapt and develop elaborate molecular mechanisms to ensure its spread and transmission. Recent findings have significantly improved our understanding of the regulation of gene expression in P. falciparum. Here, we provide an up-to-date overview of technologies used to highlight the transcriptomic adjustments occurring in the parasite throughout its life cycle. We also emphasize the complementary and complex epigenetic mechanisms regulating gene expression in malaria parasites. This review concludes with an outlook on the chromatin architecture, the remodeling systems, and how this 3D genome organization is critical in various biological processes.

Decoding the impact of nuclear organization on antigenic variation in parasites

Antigenic variation as a strategy to evade the host adaptive immune response has evolved in divergent pathogens. Antigenic variation involves restricted, and often mutually exclusive, expression of dominant antigens and a periodic switch in antigen expression during infection. In eukaryotes, nuclear compartmentalization, including three-dimensional folding of the genome and physical separation of proteins in compartments or condensates, regulates mutually exclusive gene expression and chromosomal translocations. In this Review, we discuss the impact of nuclear organization on antigenic variation in the protozoan pathogens Trypanosoma brucei and Plasmodium falciparum. In particular, we highlight the relevance of nuclear organization in both mutually exclusive antigen expression and genome stability, which underlie antigenic variation.

A type II protein arginine methyltransferase regulates merozoite invasion in Plasmodium falciparum

Protein arginine methyltransferases (PRMTs) regulate many important cellular processes, such as transcription and RNA processing in model organisms but their functions in human malaria parasites are not elucidated. Here, we characterize PfPRMT5 in Plasmodium falciparum, which catalyzes symmetric dimethylation of histone H3 at R2 (H3R2me2s) and R8, and histone H4 at R3 in vitro. PfPRMT5 disruption results in asexual stage growth defects primarily due to lower invasion efficiency of the merozoites. Transcriptomic analysis reveals down-regulation of many transcripts related to invasion upon PfPRMT5 disruption, in agreement with H3R2me2s being an active chromatin mark. Genome-wide chromatin profiling detects extensive H3R2me2s marking of genes of different cellular processes, including invasion-related genes in wildtype parasites and PfPRMT5 disruption leads to the depletion of H3R2me2s. Interactome studies identify the association of PfPRMT5 with invasion-related transcriptional regulators such as AP2-I, BDP1, and GCN5. Furthermore, PfPRMT5 is associated with the RNA splicing machinery, and PfPRMT5 disruption caused substantial anomalies in RNA splicing events, including those for invasion-related genes. In summary, PfPRMT5 is critical for regulating parasite invasion and RNA splicing in this early-branching eukaryote.

Histone modification analysis reveals common regulators of gene expression in liver and blood stage merozoites of Plasmodium parasites

Gene expression in malaria parasites is subject to various layers of regulation, including histone post-translational modifications (PTMs). Gene regulatory mechanisms have been extensively studied during the main developmental stages of Plasmodium parasites inside erythrocytes, from the ring stage following invasion to the schizont stage leading up to egress. However, gene regulation in merozoites that mediate the transition from one host cell to the next is an understudied area of parasite biology. Here, we sought to characterize gene expression and the corresponding histone PTM landscape during this stage of the parasite lifecycle through RNA-seq and ChIP-seq on P. falciparum blood stage schizonts, merozoites, and rings, as well as P. berghei liver stage merozoites. In both hepatic and erythrocytic merozoites, we identified a subset of genes with a unique histone PTM profile characterized by a region of H3K4me3 depletion in their promoter. These genes were upregulated in hepatic and erythrocytic merozoites and rings, had roles in protein export, translation, and host cell remodeling, and shared a DNA motif. These results indicate that similar regulatory mechanisms may underlie merozoite formation in the liver and blood stages. We also observed that H3K4me2 was deposited in gene bodies of gene families encoding variant surface antigens in erythrocytic merozoites, which may facilitate switching of gene expression between different members of these families. Finally, H3K18me and H2K27me were uncoupled from gene expression and were enriched around the centromeres in erythrocytic schizonts and merozoites, suggesting potential roles in the maintenance of chromosomal organization during schizogony. Together, our results demonstrate that extensive changes in gene expression and histone landscape occur during the schizont-to-ring transition to facilitate productive erythrocyte infection. The dynamic remodeling of the transcriptional program in hepatic and erythrocytic merozoites makes this stage attractive as a target for novel anti-malarial drugs that may have activity against both the liver and blood stages.

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.

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.

A coordinated transcriptional switching network mediates antigenic variation of human malaria parasites

Malaria parasites avoid immune clearance through their ability to systematically alter antigens exposed on the surface of infected red blood cells. This is accomplished by tightly regulated transcriptional control of individual members of a large, multicopy gene family called var and is the key to both the virulence and chronic nature of malaria infections. Expression of var genes is mutually exclusive and controlled epigenetically, however how large populations of parasites coordinate var gene switching to avoid premature exposure of the antigenic repertoire is unknown. Here, we provide evidence for a transcriptional network anchored by a universally conserved gene called var2csa that coordinates the switching process. We describe a structured switching bias that shifts overtime and could shape the pattern of var expression over the course of a lengthy infection. Our results provide an explanation for a previously mysterious aspect of malaria infections and shed light on how parasites possessing a relatively small repertoire of variant antigen-encoding genes can coordinate switching events to limit antigen exposure, thereby maintaining chronic infections.

Transcriptomic complexity of the human malaria parasite Plasmodium falciparum revealed by long-read sequencing

The Plasmodium falciparum human malaria parasite genome is incompletely annotated and does not accurately represent the transcriptomic diversity of this species. To address this need, we performed long-read transcriptomic sequencing. 5' capped mRNA was enriched from samples of total and nuclear-fractionated RNA from intra-erythrocytic stages and converted to cDNA library. The cDNA libraries were sequenced on PacBio and Nanopore long-read platforms. 12,495 novel isoforms were annotated from the data. Alternative 5' and 3' ends represent the majority of isoform events among the novel isoforms, with retained introns being the next most common event. The majority of alternative 5' ends correspond to genomic regions with features similar to those of the reference transcript 5' ends. However, a minority of alternative 5' ends showed markedly different features, including locations within protein-coding regions. Alternative 3' ends showed similar features to the reference transcript 3' ends, notably adenine-rich termination signals. Distinguishing features of retained introns could not be observed, except for a tendency towards shorter length and greater GC content compared with spliced introns. Expression of antisense and retained intron isoforms was detected at different intra-erythrocytic stages, suggesting developmental regulation of these isoform events. To gain insights into the possible functions of the novel isoforms, their protein-coding potential was assessed. Variants of P. falciparum proteins and novel proteins encoded by alternative open reading frames suggest that P. falciparum has a greater proteomic repertoire than the current annotation. We provide a catalog of annotated transcripts and encoded alternative proteins to support further studies on gene and protein regulation of this pathogen.

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 log2 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.

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.

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.

A Leak-Free Inducible CRISPRi/a System for Gene Functional Studies in Plasmodium falciparum

By fusing catalytically dead Cas9 (dCas9) to active domains of histone deacetylase (Sir2a) or acetyltransferase (GCN5), this CRISPR interference/activation (CRISPRi/a) system allows gene regulation at the transcriptional level without causing permanent changes in the parasite genome. However, the constitutive expression of dCas9 poses a challenge for studying essential genes, which may lead to adaptive changes in the parasite, masking the true phenotypes. Here, we developed a leak-free inducible CRISPRi/a system by integrating the DiCre/loxP regulon to allow the expression of dCas9-GCN5/-Sir2a upon transient induction with rapamycin, which allows convenient transcriptional regulation of a gene of interest by introducing a guide RNA targeting its transcription start region. Using eight genes that are either silent or expressed from low to high levels during asexual erythrocytic development, we evaluated the robustness and versatility of this system in the asexual parasites. For most genes analyzed, this inducible CRISPRi/a system led to 1.5- to 3-fold up-or downregulation of the target genes at the mRNA level. Alteration in the expression of PfK13 and PfMYST resulted in altered sensitivities to artemisinin. For autophagy-related protein 18, an essential gene related to artemisinin resistance, a >2-fold up- or downregulation was obtained by inducible CRISPRi/a, leading to growth retardation. For the master regulator of gametocytogenesis, PfAP2-G, a >10-fold increase of the PfAP2-G transcripts was obtained by CRISPRa, resulting in >4-fold higher gametocytemia in the induced parasites. Additionally, inducible CRISPRi/a could also regulate gene expression in gametocytes. This inducible epigenetic regulation system offers a fast way of studying gene functions in Plasmodium falciparum. IMPORTANCE Understanding the fundamental biology of malaria parasites through functional genetic/genomic studies is critical for identifying novel targets for antimalarial development. Conditional knockout/knockdown systems are required to study essential genes in the haploid blood stages of the parasite. In this study, we developed an inducible CRISPRi/a system via the integration of DiCre/loxP. We evaluated the robustness and versatility of this system by activating or repressing eight selected genes and achieved up- and downregulation of the targeted genes located in both the euchromatin and heterochromatin regions. This system offers the malaria research community another tool for functional genetic studies.

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.

Shared Mechanisms for Mutually Exclusive Expression and Antigenic Variation by Protozoan Parasites

Cellular decision-making at the level of gene expression is a key process in the development and evolution of every organism. Variations in gene expression can lead to phenotypic diversity and the development of subpopulations with adaptive advantages. A prime example is the mutually exclusive activation of a single gene from within a multicopy gene family. In mammals, this ranges from the activation of one of the two immunoglobulin (Ig) alleles to the choice in olfactory sensory neurons of a single odorant receptor (OR) gene from a family of more than 1,000. Similarly, in parasites like Trypanosoma brucei, Giardia lamblia or Plasmodium falciparum, the process of antigenic variation required to escape recognition by the host immune system involves the monoallelic expression of vsg, vsp or var genes, respectively. Despite the importance of this process, understanding how this choice is made remains an enigma. The development of powerful techniques such as single cell RNA-seq and Hi-C has provided new insights into the mechanisms these different systems employ to achieve monoallelic gene expression. Studies utilizing these techniques have shown how the complex interplay between nuclear architecture, physical interactions between chromosomes and different chromatin states lead to single allele expression. Additionally, in several instances it has been observed that high-level expression of a single gene is preceded by a transient state where multiple genes are expressed at a low level. In this review, we will describe and compare the different strategies that organisms have evolved to choose one gene from within a large family and how parasites employ this strategy to ensure survival within their hosts.

A dynamic and combinatorial histone code drives malaria parasite asexual and sexual development

Histone posttranslational modifications (PTMs) frequently co-occur on the same chromatin domains or even in the same molecule. It is now established that these “histone codes” are the result of cross talk between enzymes that catalyze multiple PTMs with univocal readout as compared with these PTMs in isolation. Here, we performed a comprehensive identification and quantification of histone codes of the malaria parasite, Plasmodium falciparum. We used advanced quantitative middle-down proteomics to identify combinations of PTMs in both the proliferative, asexual stages and transmissible, sexual gametocyte stages of P. falciparum. We provide an updated, high-resolution compendium of 77 PTMs on H3 and H3.3, of which 34 are newly identified in P. falciparum. Coexisting PTMs with unique stage distinctions were identified, indicating that many of these combinatorial PTMs are associated with specific stages of the parasite life cycle. We focused on the code H3R17me2K18acK23ac for its unique presence in mature gametocytes; chromatin proteomics identified a gametocyte-specific SAGA-like effector complex including the transcription factor AP2-G2, which we tied to this specific histone code, as involved in regulating gene expression in mature gametocytes. Ultimately, this study unveils previously undiscovered histone PTMs and their functional relationship with coexisting partners. These results highlight that investigating chromatin regulation in the parasite using single histone PTM assays might overlook higher-order gene regulation for distinct proliferation and differentiation processes.

•

First middle-down chromatin proteomics compendium of the malaria parasite, Plasmodium falciparum.

•

Novel histone PTMs (including arginine methylation) in both asexual parasites and transmissible gametocytes.

•

Histone PTM cross talk is dynamic life cycle stage stratified.

•

Gametocytes rely on histone PTM connectivity to allow onward transmission.

•

AP2-G2 is an important effector of H3K18acK23ac in mature gametocytes.

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.

An Uninvited Seat at the Dinner Table: How Apicomplexan Parasites Scavenge Nutrients from the Host

Obligate intracellular parasites have evolved a remarkable assortment of strategies to scavenge nutrients from the host cells they parasitize. Most apicomplexans form a parasitophorous vacuole (PV) within the invaded cell, a replicative niche within which they survive and multiply. As well as providing a physical barrier against host cell defense mechanisms, the PV membrane (PVM) is also an important site of nutrient uptake that is essential for the parasites to sustain their metabolism. This means nutrients in the extracellular milieu are separated from parasite metabolic machinery by three different membranes, the host plasma membrane, the PVM, and the parasite plasma membrane (PPM). In order to facilitate nutrient transport from the extracellular environment into the parasite itself, transporters on the host cell membrane of invaded cells can be modified by secreted and exported parasite proteins to maximize uptake of key substrates to meet their metabolic demand. To overcome the second barrier, the PVM, apicomplexan parasites secrete proteins contained in the dense granules that remodel the vacuole and make the membrane permissive to important nutrients. This bulk flow of host nutrients is followed by a more selective uptake of substrates at the PPM that is operated by specific transporters of this third barrier. In this review, we recapitulate and compare the strategies developed by Apicomplexa to scavenge nutrients from their hosts, with particular emphasis on transporters at the parasite plasma membrane and vacuolar solute transporters on the parasite intracellular digestive organelle.

Identifying transcript 5' capped ends in Plasmodium falciparum

Background

The genome of the human malaria parasite Plasmodium falciparum is poorly annotated, in particular, the 5' capped ends of its mRNA transcripts. New approaches are needed to fully catalog P. falciparum transcripts for understanding gene function and regulation in this organism.

Methods

We developed a transcriptomic method based on next-generation sequencing of complementary DNA (cDNA) enriched for full-length fragments using eIF4E, a 5' cap-binding protein, and an unenriched control. DNA sequencing adapter was added after enrichment of full-length cDNA using two different ligation protocols. From the mapped sequence reads, enrichment scores were calculated for all transcribed nucleotides and used to calculate P-values of 5' capped nucleotide enrichment. Sensitivity and accuracy were increased by combining P-values from replicate experiments. Data were obtained for P. falciparum ring, trophozoite and schizont stages of intra-erythrocytic development.

Results

5' capped nucleotide signals were mapped to 17,961 non-overlapping P. falciparum genomic intervals. Analysis of the dominant 5' capped nucleotide in these genomic intervals revealed the presence of two groups with distinctive epigenetic features and sequence patterns. A total of 4,512 transcripts were annotated as 5' capped based on the correspondence of 5' end with 5' capped nucleotide annotated from full-length cDNA data.

Discussion

The presence of two groups of 5' capped nucleotides suggests that alternative mechanisms may exist for producing 5' capped transcript ends in P. falciparum. The 5' capped transcripts that are antisense, outside of, or partially overlapping coding regions may be important regulators of gene function in P. falciparum.

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.

Peculiarities of Plasmodium falciparum Gene Regulation and Chromatin Structure

The highly complex life cycle of the human malaria parasite, Plasmodium falciparum, is based on an orchestrated and tightly regulated gene expression program. In general, eukaryotic transcription regulation is determined by a combination of sequence-specific transcription factors binding to regulatory DNA elements and the packaging of DNA into chromatin as an additional layer. The accessibility of regulatory DNA elements is controlled by the nucleosome occupancy and changes of their positions by an active process called nucleosome remodeling. These epigenetic mechanisms are poorly explored in P. falciparum. The parasite genome is characterized by an extraordinarily high AT-content and the distinct architecture of functional elements, and chromatin-related proteins also exhibit high sequence divergence compared to other eukaryotes. Together with the distinct biochemical properties of nucleosomes, these features suggest substantial differences in chromatin-dependent regulation. Here, we highlight the peculiarities of epigenetic mechanisms in P. falciparum, addressing chromatin structure and dynamics with respect to their impact on transcriptional control. We focus on the specialized chromatin remodeling enzymes and discuss their essential function in P. falciparum gene regulation.

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

Dynamic association of the H3K64 trimethylation mark with genes encoding exported proteins in Plasmodium falciparum

Epigenetic modifications have emerged as critical regulators of virulence genes and stage-specific gene expression in Plasmodium falciparum. However, the specific roles of histone core epigenetic modifications in regulating the stage-specific gene expression are not well understood. In this study, we report an unconventional trimethylation at lysine 64 on histone 3 (H3K64me3) and characterize its functional relevance in P. falciparum. We show that PfSET4 and PfSET5 proteins of P. falciparum methylate H3K64 and that they prefer the nucleosome as a substrate over free histone 3 proteins. Structural analysis of PfSET5 revealed that it interacts with the nucleosome as a dimer. The H3K64me3 mark is dynamic, being enriched in the ring and trophozoite stages and drastically reduced in the schizont stages. Stage-specific global chromatin immunoprecipitation –sequencing analysis of the H3K64me3 mark revealed the selective enrichment of this methyl mark on the genes of exported family proteins in the ring and trophozoite stages and a significant reduction of the same in the schizont stages. Collectively, our data identify a novel epigenetic mark that is associated with the subset of genes encoding for exported proteins, which may regulate their expression in different stages of P. falciparum.

Plasmodium Reproduction, Cell Size, and Transcription: How to Cope With Increasing DNA Content?

Plasmodium, the unicellular parasite that causes malaria, evolved a highly unusual mode of reproduction. During its complex life cycle, invasive or transmissive stages alternate with proliferating stages, where a single parasite can produce tens of thousands of progeny. In the clinically relevant blood stage of infection, the parasite replicates its genome up to thirty times and forms a multinucleated cell before daughter cells are assembled. Thus, within a single cell cycle, Plasmodium develops from a haploid to a polypoid cell, harboring multiple copies of its genome. Polyploidy creates several biological challenges, such as imbalances in genome output, and cells can respond to this by changing their size and/or alter the production of RNA species and protein to achieve expression homeostasis. However, the effects and possible adaptations of Plasmodium to the massively increasing DNA content are unknown. Here, we revisit and embed current Plasmodium literature in the context of polyploidy and propose potential mechanisms of the parasite to cope with the increasing gene dosage.

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.

Identification of long regulatory elements in the genome of Plasmodium falciparum and other eukaryotes

Long regulatory elements (LREs), such as CpG islands, polydA:dT tracts or AU-rich elements, are thought to play key roles in gene regulation but, as opposed to conventional binding sites of transcription factors, few methods have been proposed to formally and automatically characterize them. We present here a computational approach named DExTER (Domain Exploration To Explain gene Regulation) dedicated to the identification of candidate LREs (cLREs) and apply it to the analysis of the genomes of P. falciparum and other eukaryotes. Our analyses show that all tested genomes contain several cLREs that are somewhat conserved along evolution, and that gene expression can be predicted with surprising accuracy on the basis of these long regions only. Regulation by cLREs exhibits very different behaviours depending on species and conditions. In P. falciparum and other Apicomplexan organisms as well as in Dictyostelium discoideum, the process appears highly dynamic, with different cLREs involved at different phases of the life cycle. For multicellular organisms, the same cLREs are involved in all tissues, but a dynamic behavior is observed along embryonic development stages. In P. falciparum, whose genome is known to be strongly depleted of transcription factors, cLREs are predictive of expression with an accuracy above 70%, and our analyses show that they are associated with both transcriptional and post-transcriptional regulation signals. Moreover, we assessed the biological relevance of one LRE discovered by DExTER in P. falciparum using an in vivo reporter assay. The source code (python) of DExTER is available at https://gite.lirmm.fr/menichelli/DExTER.

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.

Histone acetyltransferase PfGCN5 regulates stress responsive and artemisinin resistance related genes in Plasmodium falciparum

Plasmodium falciparum has evolved resistance to almost all front-line drugs including artemisinin, which threatens malaria control and elimination strategies. Oxidative stress and protein damage responses have emerged as key players in the generation of artemisinin resistance. In this study, we show that PfGCN5, a histone acetyltransferase, binds to the stress-responsive genes in a poised state and regulates their expression under stress conditions. Furthermore, we show that upon artemisinin exposure, genome-wide binding sites for PfGCN5 are increased and it is directly associated with the genes implicated in artemisinin resistance generation like BiP and TRiC chaperone. Interestingly, expression of genes bound by PfGCN5 was found to be upregulated during stress conditions. Moreover, inhibition of PfGCN5 in artemisinin-resistant parasites increases the sensitivity of the parasites to artemisinin treatment indicating its role in drug resistance generation. Together, these findings elucidate the role of PfGCN5 as a global chromatin regulator of stress-responses with a potential role in modulating artemisinin drug resistance and identify PfGCN5 as an important target against artemisinin-resistant parasites.

Phosphorylation of the Canonical Histone H2A Marks Foci of Damaged DNA in Malaria Parasites

Plasmodium falciparum is the deadliest human parasite that causes malaria when it reaches the bloodstream and begins proliferating inside red blood cells, where the parasites are particularly prone to DNA damage. The molecular mechanisms that allow these pathogens to maintain their genome integrity under such conditions are also the driving force for acquiring genome plasticity that enables them to create antigenic variation and become resistant to essentially all available drugs.

Plasmodium falciparum parasites proliferate within circulating red blood cells and are responsible for the deadliest form of human malaria. These parasites are exposed to numerous intrinsic and external sources that could cause DNA damage; therefore, they have evolved efficient mechanisms to protect their genome integrity and allow them to proliferate under such conditions. In higher eukaryotes, double-strand breaks rapidly lead to phosphorylation of the core histone variant H2A.X, which marks the site of damaged DNA. We show that in P. falciparum that lacks the H2A.X variant, the canonical P. falciparum H2A (PfH2A) is phosphorylated on serine 121 upon exposure to sources of DNA damage. We further demonstrate that phosphorylated PfH2A is recruited to foci of damaged chromatin shortly after exposure to sources of damage, while the nonphosphorylated PfH2A remains spread throughout the nucleoplasm. In addition, we found that PfH2A phosphorylation is dynamic and that over time, as the parasite activates the repair machinery, this phosphorylation is removed. Finally, we demonstrate that these phosphorylation dynamics could be used to establish a novel and direct DNA repair assay in P. falciparum.

IMPORTANCEPlasmodium falciparum is the deadliest human parasite that causes malaria when it reaches the bloodstream and begins proliferating inside red blood cells, where the parasites are particularly prone to DNA damage. The molecular mechanisms that allow these pathogens to maintain their genome integrity under such conditions are also the driving force for acquiring genome plasticity that enables them to create antigenic variation and become resistant to essentially all available drugs. However, mechanisms of DNA damage response and repair have not been extensively studied for these parasites. The paper addresses our recent discovery that P. falciparum that lacks the histone variant H2A.X phosphorylates its canonical core histone PfH2A in response to exposure to DNA damage. The process of DNA repair in Plasmodium was mostly studied indirectly. Our findings enabled us to establish a direct DNA repair assay for P. falciparum similar to assays that are widely used in model organisms.

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.

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.

Elucidation of DNA Repair Function of PfBlm and Potentiation of Artemisinin Action by a Small-Molecule Inhibitor of RecQ Helicase

Malaria continues to be a serious threat to humankind not only because of the morbidity and mortality associated with the disease but also due to the huge economic burden that it imparts. Resistance to all available drugs and the unavailability of an effective vaccine cry for an urgent discovery of newer drug targets.

Artemisinin (ART)-based combination therapies are recommended as first- and second-line treatments for Plasmodium falciparum malaria. Here, we investigated the impact of the RecQ inhibitor ML216 on the repair of ART-mediated damage in the genome of P. falciparum. PfBLM and PfWRN were identified as members of the RecQ helicase family in P. falciparum. However, the role of these RecQ helicases in DNA double-strand break (DSB) repair in this parasite has not been explored. Here, we provide several lines of evidence to establish the involvement of PfBlm in DSB repair in P. falciparum. First, we demonstrate that PfBlm interacts with two well-characterized DSB repair proteins of this parasite, namely, PfRad51 and PfalMre11. Second, we found that PfBLM expression was upregulated in response to DNA-damaging agents. Third, through yeast complementation studies, we demonstrated that PfBLM could complement the DNA damage sensitivity of a Δsgs1 mutant of Saccharomyces cerevisiae, in contrast to the helicase-dead mutant PfblmK83R. Finally, we observe that the overexpression of PfBLM induces resistance to DNA-damaging agents and offers a survival advantage to the parasites. Most importantly, we found that the RecQ inhibitor ML216 inhibits the repair of DSBs and thereby renders parasites more sensitive to ART. Such synergism between ART and ML216 actions was observed for both drug-sensitive and multidrug-resistant strains of P. falciparum. Taken together, these findings establish the implications of PfBlm in the Plasmodium DSB repair pathway and provide insights into the antiparasitic activity of the ART-ML216 combination.

IMPORTANCE Malaria continues to be a serious threat to humankind not only because of the morbidity and mortality associated with the disease but also due to the huge economic burden that it imparts. Resistance to all available drugs and the unavailability of an effective vaccine cry for an urgent discovery of newer drug targets. Here, we uncovered a role of the PfBlm helicase in Plasmodium DNA double-strand break repair and established that the parasitic DNA repair mechanism can be targeted to curb malaria. The small-molecule inhibitor of PfBlm tested in this study acts synergistically with two first-line malaria drugs, artemisinin (ART) and chloroquine, in both drug-sensitive and multidrug-resistant strains of P. falciparum, thus qualifying this chemical as a potential partner in ART-based combination therapy. Additionally, the identification of this new specific inhibitor of the Plasmodium homologous recombination (HR) mechanism will now allow us to investigate the role of HR in Plasmodium biology.

Protein Modification Characteristics of the Malaria Parasite Plasmodium falciparum and the Infected Erythrocytes

Malaria elimination is still pending on the development of novel tools that rely on a deep understanding of parasite biology. Proteins of all living cells undergo myriad posttranslational modifications (PTMs) that are critical to multifarious life processes. An extensive proteome-wide dissection revealed a fine PTM map of most proteins in both Plasmodium falciparum, the causative agent of severe malaria, and the infected red blood cells. More than two-thirds of proteins of the parasite and its host cell underwent extensive and dynamic modification throughout the erythrocytic developmental stage. PTMs critically modulate the virulence factors involved in the host-parasite interaction and pathogenesis. Furthermore, P. falciparum stabilized the supporting proteins of erythrocyte origin by selective demodification. Collectively, our multiple omic analyses, apart from having furthered a deep understanding of the systems biology of P. falciparum and malaria pathogenesis, provide a valuable resource for mining new antimalarial targets.

Histone modifications associated with gene expression and genome accessibility are dynamically enriched at Plasmodium falciparum regulatory sequences

Background

The malaria parasite Plasmodium falciparum has an unusually euchromatic genome with poorly conserved positioning of nucleosomes in intergenic sequences and poorly understood mechanisms of gene regulation. Variant histones and histone modifications determine nucleosome stability and recruit trans factors, but their combinatorial contribution to gene regulation is unclear.

Results

Here, we show that the histone H3 acetylations H3K18ac and H3K27ac and the variant histone Pf H2A.Z are enriched together at regulatory sites upstream of genes. H3K18ac and H3K27ac together dynamically mark regulatory regions of genes expressed during the asexual life cycle. In contrast, H3K4me1 is depleted in intergenic sequence and dynamically depleted upstream of expressed genes. The temporal pattern of H3K27ac and H3K18ac enrichment indicates that they accumulate during S phase and mitosis and are retained at regulatory sequences until at least G1 phase and after cessation of expression of the cognate genes. We integrated our ChIPseq data with existing datasets to show that in schizont stages H3K18ac, H3K27ac and Pf H2A.Z colocalise with the transcription factor PfAP2-I and the bromodomain protein PfBDP1 and are enriched at stably positioned nucleosomes within regions of exposed DNA at active transcriptional start sites. Using transient transfections we showed that sequences enriched with colocalised H3K18ac, H3K27ac and Pf H2A.Z possess promoter activity in schizont stages, but no enhancer-like activity.

Conclusions

The dynamic H3 acetylations define P. falciparum regulatory sequences and contribute to gene activation. These findings expand the knowledge of the chromatin landscape that regulates gene expression in P. falciparum.

Plasmodium falciparum replication factor C subunit 1 is involved in genotoxic stress response

About half the world's population is at risk of malaria, with Plasmodium falciparum malaria being responsible for the most malaria related deaths globally. Antimalarial drugs such as chloroquine and artemisinin are directed towards the proliferating intra-erythrocytic stages of the parasite, which is responsible for all the clinical symptoms of the disease. These antimalarial drugs have been reported to function via multiple pathways, one of which induces DNA damage via the generation of free radicals and reactive oxygen species. An urgent need to understand the mechanistic details of drug response and resistance is highlighted by the decreasing clinical efficacy of the front line drug, Artemisinin. The replication factor C subunit 1 is an important component of the DNA replication machinery and DNA damage response mechanism. Here we show the translocation of PfRFC1 from an intranuclear localisation to the nuclear periphery, indicating an orchestrated progression of distinct patterns of replication in the developing parasites. PfRFC1 responds to genotoxic stress via elevated protein levels in soluble and chromatin bound fractions. Reduction of PfRFC1 protein levels upon treatment with antimalarials suggests an interplay of replication, apoptosis and DNA repair pathways leading to cell death. Additionally, mislocalisation of the endogenously tagged protein confirmed its essential role in parasites' replication and DNA repair. This study provides key insights into DNA replication, DNA damage response and cell death in P. falciparum.

The human malaria parasite genome is configured into thousands of coexpressed linear regulatory units

The human malaria parasite Plasmodium falciparum thrives in radically different host environments in mosquitoes and humans, with only a limited set of transcription factors. The nature of regulatory elements or their target genes in the P. falciparum genome remains elusive. Here, we found that this eukaryotic parasite uses an efficient way to maximally use genetic and epigenetic regulation to form regulatory units (RUs) during blood infections. Genes located in the same RU tend to have the same pattern of expression over time and are associated with open chromatin along regulatory elements. To precisely define and quantify these RUs, a novel hidden Markov model was developed to capture the regulatory structure in a genome-wide fashion by integrating expression and epigenetic evidence. We successfully identified thousands of RUs and cross-validated with previous findings. We found more genes involved in red blood cell (RBC) invasion located in the same RU as the PfAP2-I (AP2-I) transcription factor, demonstrating that AP2-I is responsible for regulating RBC invasion. Our study has provided a regulatory mechanism for a compact eukaryotic genome and offers new insights into the in vivo transcriptional regulation of the P. falciparum intraerythrocytic stage.

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.