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

Histone H3.3 variant plays a critical role on zygote-to-oocyst development in malaria parasites

The Plasmodium life cycle involves differentiation into multiple morphologically distinct forms, a process regulated by developmental stage-specific gene expression. Histone proteins are involved in epigenetic regulation in eukaryotes, and the histone variant H3.3 plays a key role in the regulation of gene expression and maintenance of genomic integrity during embryonic development in mice. However, the function of H3.3 through multiple developmental stages in Plasmodium remains unknown. To examine the function of H3.3, h3.3-deficient mutants (Δh3.3) were generated in P. berghei. The deletion of h3.3 was not lethal in blood stage parasites, although it had a minor effect of the growth rate in blood stage; however, the in vitro ookinete conversion rate was significantly reduced, and the production of the degenerated form was increased. Regarding the mosquito stage development of Δh3.3, oocysts number was significantly reduced, and no sporozoite production was observed. The h3.3 gene complemented mutant have normal development in mosquito stage producing mature oocysts and salivary glands contained sporozoites, and interestingly, the majority of H3.3 protein was detected in female gametocytes. However, Δh3.3 male and female gametocyte production levels were comparable to the wild-type levels. Transcriptome analysis of Δh3.3 male and female gametocytes revealed the upregulation of several male-specific genes in female gametocytes, suggesting that H3.3 functions as a transcription repressor of male-specific genes to maintain sexual identity in female gametocytes. This study provides new insights into the molecular biology of histone variants H3.3 which plays a critical role on zygote-to-oocyst development in primitive unicellular eukaryotes.

Histone code: a common language and multiple dialects to meet the different developmental requirements of apicomplexan parasites

Apicomplexan parasites have complex life cycles, often requiring transmission between two different hosts, facing periods of dormancy within the host or in the environment to maximize chances of transmission. To support survival in these different conditions, tightly regulated and correctly timed gene expression is critical. The modification of histones and nucleosome composition makes a significant contribution to this regulation, and as eukaryotes, the fundamental mechanisms underlying this process in apicomplexans are similar to those in model eukaryotic organisms. However, single-celled intracellular parasites face unique challenges, and regulation of gene expression at the epigenetic level provides tight control for responses that must often be rapid and robust. Here, we discuss the recent advances in understanding the dynamics of histone modifications across Apicomplexan life cycles and the molecular mechanisms that underlie epigenetic regulation of gene expression to promote parasite life cycle progression, dormancy, and transmission.

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.

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

The Plasmodium falciparum alternative histones Pf H2A.Z and Pf H2B.Z are enriched in the same nucleosomes in intergenic euchromatin but depleted from heterochromatin. They occupy most promoters but are only dynamically associated with expression at var genes. In other organisms, acetylation of H2A.Z is important for its functions in gene expression and chromatin structure. Here, we show that acetylated Pf H2A.Z and Pf H2B.Z are dynamically associated with gene expression at promoters. In addition, acetylated Pf H2A.Z and Pf H2B.Z are antagonized by the sirtuin class III histone deacetylases (HDAC) PfSir2A and B at heterochromatin boundaries and encroach upon heterochromatin in parasites lacking PfSir2A or B. However, the majority of acetylated Pf H2A.Z and Pf H2B.Z are deacetylated by class I or II HDACs. Acetylated Pf H2A.Z and Pf H2B.Z are also dynamically associated with promoter activity of both canonical upstream var gene promoters and var gene introns. These findings suggest that both acetylated Pf H2A.Z and Pf H2B.Z play critical roles in gene expression and contribute to maintenance of chromatin structure at the boundaries of subtelomeric, facultative heterochromatin, critical for the variegated expression of genes that enable rapid adaptation to altered host environments.IMPORTANCEThe 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 GCN5 plays a key role in regulating artemisinin resistance-related stress responses

Plasmodium falciparum causes the most severe malaria and is exposed to various environmental and physiological stresses in the human host. Given that GCN5 plays a critical role in regulating stress responses in model organisms, we aimed to elucidate PfGCN5's function in stress responses in P. falciparum. The protein level of PfGCN5 was substantially induced under three stress conditions [heat shock, low glucose starvation, and dihydroartemisinin, the active metabolite of artemisinin (ART)]. With a TetR-DOZI conditional knockdown (KD) system, we successfully down-regulated PfGCN5 to ~50% and found that KD parasites became more sensitive to all three stress conditions. Transcriptomic analysis via RNA-seq identified ~1,000 up- and down-regulated genes in the wild-type (WT) and KD parasites under these stress conditions. Importantly, DHA induced transcriptional alteration of many genes involved in many aspects of stress responses, which were heavily shared among the altered genes under heat shock and low glucose conditions, including ART-resistance-related genes such as K13 and coronin. Based on the expression pattern between WT and KD parasites under three stress conditions, ~300-400 genes were identified to be involved in PfGCN5-dependent, general, and stress-condition-specific responses with high levels of overlaps among three stress conditions. Notably, using ring-stage survival assay, we found that KD or inhibition of PfGCN5 could sensitize the ART-resistant parasites to the DHA treatment. All these indicate that PfGCN5 is pivotal in regulating general and ART-resistance-related stress responses in malaria parasites, implicating PfGCN5 as a potential target for malaria intervention.

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.

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.

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

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

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.

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.

Structural insights into acetylated histone ligand recognition by the BDP1 bromodomain of Plasmodium falciparum

Plasmodium falciparum requires a two-host system, moving between Anopheles mosquito and humans, to complete its life cycle. To overcome such dynamic growth conditions its histones undergo various post-translational modifications to regulate gene expression. The P. falciparum Bromodomain Protein 1 (PfBDP1) has been shown to interact with acetylated lysine modifications on histone H3 to regulate the expression of invasion-related genes. Here, we investigated the ability of the PfBDP1 bromodomain to interact with acetyllsyine modifications on additional core and variant histones. A crystal structure of the PfBDP1 bromodomain (PfBDP1-BRD) reveals it contains the conserved bromodomain fold, but our comparative analysis between the PfBDP1-BRD and human bromodomain families indicates it has a unique binding mechanism. Solution NMR spectroscopy and ITC binding assays carried out with acetylated histone ligands demonstrate that it preferentially recognizes tetra-acetylated histone H4, and we detected weaker interactions with multi-acetylated H2A.Z in addition to the previously reported interactions with acetylated histone H3. Our findings indicate PfBDP1 may play additional roles in the P. falciparum life cycle, and the distinctive features of its bromodomain binding pocket could be leveraged for the development of new therapeutic agents to help overcome the continuously evolving resistance of P. falciparum against currently available drugs.

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.

Plasmodium 6-Cysteine Proteins: Functional Diversity, Transmission-Blocking Antibodies and Structural Scaffolds

The 6-cysteine protein family is one of the most abundant surface antigens that are expressed throughout the Plasmodium falciparum life cycle. Many members of the 6-cysteine family have critical roles in parasite development across the life cycle in parasite transmission, evasion of the host immune response and host cell invasion. The common feature of the family is the 6-cysteine domain, also referred to as s48/45 domain, which is conserved across Aconoidasida. This review summarizes the current approaches for recombinant expression for 6-cysteine proteins, monoclonal antibodies against 6-cysteine proteins that block transmission and the growing collection of crystal structures that provide insights into the functional domains of this protein family.

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.

Genome-wide localization of histone variants in Toxoplasma gondii implicates variant exchange in stage-specific gene expression

BACKGROUND

Toxoplasma gondii is a protozoan parasite that differentiates from acute tachyzoite stages to latent bradyzoite forms in response to environmental cues that modify the epigenome. We studied the distribution of the histone variants CenH3, H3.3, H2A.X, H2A.Z and H2B.Z, by genome-wide chromatin immunoprecipitation to understand the role of variant histones in developmental transitions of T. gondii parasites.

RESULTS

H3.3 and H2A.X were detected in telomere and telomere associated sequences, whereas H3.3, H2A.X and CenH3 were enriched in centromeres. Histones H2A.Z and H2B.Z colocalize with the transcriptional activation mark H3K4me3 in promoter regions surrounding the nucleosome-free region upstream of the transcription start site. The H2B.Z/H2A.Z histone pair also localizes to the gene bodies of genes that are silent but poised for activation, including bradyzoite stage-specific genes. The majority of H2A.X and H2A.Z/H2B.Z loci do not overlap, consistent with variant histones demarcating specific functional regions of chromatin. The extent of enrichment of H2A.Z/H2B.Z (and H3.3 and H2A.X) within the entire gene (5'UTR and gene body) reflects the timing of gene expression during the cell cycle, suggesting that dynamic turnover of H2B.Z/H2A.Z occurs during the tachyzoite cell cycle. Thus, the distribution of the variant histone H2A.Z/H2B.Z dimer defines active and developmentally silenced regions of the T. gondii epigenome including genes that are poised for expression.

CONCLUSIONS

Histone variants mark functional regions of parasite genomes with the dynamic placement of the H2A.Z/H2B.Z dimer implicated as an evolutionarily conserved regulator of parasite and eukaryotic differentiation.

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.

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

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.

Function of cofactor Akirin2 in the regulation of gene expression in model human Caucasian neutrophil-like HL60 cells

The Akirin family of transcription cofactors are involved throughout the metazoan in the regulation of different biological processes (BPs) such as immunity, interdigital regression, muscle and neural development. Akirin do not have catalytic or DNA-binding capability and exert its regulatory function primarily through interacting proteins such as transcription factors, chromatin remodelers, and RNA-associated proteins. In the present study, we focused on the human Akirin2 regulome and interactome in neutrophil-like model human Caucasian promyelocytic leukemia HL60 cells. Our hypothesis is that metazoan evolved to have Akirin2 functional complements and different Akirin2-mediated mechanisms for the regulation of gene expression. To address this hypothesis, experiments were conducted using transcriptomics, proteomics and systems biology approaches in akirin2 knockdown and wildtype (WT) HL60 cells to characterize Akirin2 gene/protein targets, functional complements and to provide evidence of different mechanisms that may be involved in Akirin2-mediated regulation of gene expression. The results revealed Akirin2 gene/protein targets in multiple BPs with higher representation of immunity and identified immune response genes as candidate Akirin2 functional complements. In addition to linking chromatin remodelers with transcriptional activation, Akirin2 also interacts with histone H3.1 for regulation of gene expression.