The nucleosome landscape of Plasmodium falciparum reveals chromatin architecture and dynamics of regulatory sequences

Nascent RNA sequencing reveals mechanisms of gene regulation in the human malaria parasite Plasmodium falciparum

Gene expression in Plasmodium falciparum is tightly regulated to ensure successful propagation of the parasite throughout its complex life cycle. The earliest transcriptomics studies in P. falciparum suggested a cascade of transcriptional activity over the course of the 48-hour intraerythrocytic developmental cycle (IDC); however, the just-in-time transcriptional model has recently been challenged by findings that show the importance of post-transcriptional regulation. To further explore the role of transcriptional regulation, we performed the first genome-wide nascent RNA profiling in P. falciparum. Our findings indicate that the majority of genes are transcribed simultaneously during the trophozoite stage of the IDC and that only a small subset of genes is subject to differential transcriptional timing. RNA polymerase II is engaged with promoter regions prior to this transcriptional burst, suggesting that Pol II pausing plays a dominant role in gene regulation. In addition, we found that the overall transcriptional program during gametocyte differentiation is surprisingly similar to the IDC, with the exception of relatively small subsets of genes. Results from this study suggest that further characterization of the molecular players that regulate stage-specific gene expression and Pol II pausing will contribute to our continuous search for novel antimalarial drug targets.

GT-rich promoters can drive RNA pol II transcription and deposition of H2A.Z in African trypanosomes

Genome‐wide transcription studies are revealing an increasing number of “dispersed promoters” that, unlike “focused promoters”, lack well‐conserved sequence motifs and tight regulation. Dispersed promoters are nevertheless marked by well‐defined chromatin structures, suggesting that specific sequence elements must exist in these unregulated promoters. Here, we have analyzed regions of transcription initiation in the eukaryotic parasite Trypanosoma brucei, in which RNA polymerase II transcription initiation occurs over broad regions without distinct promoter motifs and lacks regulation. Using a combination of site‐specific and genome‐wide assays, we identified GT‐rich promoters that can drive transcription and promote the targeted deposition of the histone variant H2A.Z in a genomic context‐dependent manner. In addition, upon mapping nucleosome occupancy at high resolution, we find nucleosome positioning to correlate with RNA pol II enrichment and gene expression, pointing to a role in RNA maturation. Nucleosome positioning may thus represent a previously unrecognized layer of gene regulation in trypanosomes. Our findings show that even highly dispersed, unregulated promoters contain specific DNA elements that are able to induce transcription and changes in chromatin structure.

Malaria Epigenetics

Organisms with identical genome sequences can show substantial differences in their phenotypes owing to epigenetic changes that result in different use of their genes. Epigenetic regulation of gene expression plays a key role in the control of several fundamental processes in the biology of malaria parasites, including antigenic variation and sexual differentiation. Some of the histone modifications and chromatin-modifying enzymes that control the epigenetic states of malaria genes have been characterized, and their functions are beginning to be unraveled. The fundamental principles of epigenetic regulation of gene expression appear to be conserved between malaria parasites and model eukaryotes, but important peculiarities exist. Here, we review the current knowledge of malaria epigenetics and discuss how it can be exploited for the development of new molecular markers and new types of drugs that may contribute to malaria eradication efforts.

The Role of Chromatin Structure in Gene Regulation of the Human Malaria Parasite

The human malaria parasite, Plasmodium falciparum, depends on a coordinated regulation of gene expression for development and propagation within the human host. Recent developments suggest that gene regulation in the parasite is largely controlled by epigenetic mechanisms. Here, we discuss recent advancements contributing to our understanding of the mechanisms controlling gene regulation in the parasite, including nucleosome landscape, histone modifications, and nuclear architecture. In addition, various processes involved in regulation of parasite-specific genes and gene families are examined. Finally, we address the use of epigenetic processes as targets for novel antimalarial therapies. Collectively, these topics highlight the unique biology of P. falciparum, and contribute to our understanding of mechanisms regulating gene expression in this deadly parasite.

Epigenetic landscapes underlining global patterns of gene expression in the human malaria parasite, Plasmodium falciparum

The dynamic chromatin landscape displaying combinatorial complexity of the epigenome impacts gene expression that underlies many events of differentiation and cell cycle progression. In the past few years, epigenetic mechanisms have emerged as important processes involved in the tight gene regulation in malaria parasites, Plasmodium spp. Focusing predominantly on Plasmodium falciparum, the species associated with the most severe form of the disease, many advances have been made in our understanding of the interaction between transcriptional regulation and epigenetic mechanisms as the pivotal processes in regulating life cycle progression, host parasite interactions and parasite adaptation to the host environment. This review focuses on the epigenome and its effect on transcriptional regulation in P. falciparum, highlighting its unique, evolutionary diverse features.

Malaria parasites possess a telomere repeat-binding protein that shares ancestry with transcription factor IIIA

Telomere repeat-binding factors (TRFs) are essential components of the molecular machinery that regulates telomere function. TRFs are widely conserved across eukaryotes and bind duplex telomere repeats via a characteristic MYB-type domain. Here, we identified the telomere repeat-binding protein PfTRZ in the malaria parasite Plasmodium falciparum, a member of the Alveolate phylum for which TRFs have not been described so far. PfTRZ lacks an MYB domain and binds telomere repeats via a C₂H₂-type zinc finger domain instead. In vivo, PfTRZ binds with high specificity to the telomeric tract and to interstitial telomere repeats upstream of subtelomeric virulence genes. Conditional depletion experiments revealed that PfTRZ regulates telomere length homeostasis and is required for efficient cell cycle progression. Intriguingly, we found that PfTRZ also binds to and regulates the expression of 5S rDNA genes. Combined with detailed phylogenetic analyses, our findings identified PfTRZ as a remote functional homologue of the basic transcription factor TFIIIA, which acquired a new function in telomere maintenance early in the apicomplexan lineage. Our work sheds unexpected new light on the evolution of telomere repeat-binding proteins and paves the way for dissecting the presumably divergent mechanisms regulating telomere functionality in one of the most deadly human pathogens.

Epigenetic regulation of Plasmodium falciparum clonally variant gene expression during development in Anopheles gambiae

P. falciparum phenotypic plasticity is linked to the variant expression of clonal multigene families such as the var genes. We have examined changes in transcription and histone modifications that occur during sporogonic development of P. falciparum in the mosquito host. All var genes are silenced or transcribed at low levels in blood stages (gametocyte/ring) of the parasite in the human host. After infection of mosquitoes, a single var gene is selected for expression in the oocyst, and transcription of this gene increases dramatically in the sporozoite. The same PF3D7_1255200 var gene was activated in 4 different experimental infections. Transcription of this var gene during parasite development in the mosquito correlates with the presence of low levels of H3K9me3 at the binding site for the PF3D7_1466400 AP2 transcription factor. This chromatin state in the sporozoite also correlates with the expression of an antisense long non-coding RNA (lncRNA) that has previously been shown to promote var gene transcription during the intraerythrocytic cycle in vitro. Expression of both the sense protein-coding transcript and the antisense lncRNA increase dramatically in sporozoites. The findings suggest a complex process for the activation of a single particular var gene that involves AP2 transcription factors and lncRNAs.

Plasmodium falciparum Nucleosomes Exhibit Reduced Stability and Lost Sequence Dependent Nucleosome Positioning

The packaging and organization of genomic DNA into chromatin represents an additional regulatory layer of gene expression, with specific nucleosome positions that restrict the accessibility of regulatory DNA elements. The mechanisms that position nucleosomes in vivo are thought to depend on the biophysical properties of the histones, sequence patterns, like phased di-nucleotide repeats and the architecture of the histone octamer that folds DNA in 1.65 tight turns. Comparative studies of human and P. falciparum histones reveal that the latter have a strongly reduced ability to recognize internal sequence dependent nucleosome positioning signals. In contrast, the nucleosomes are positioned by AT-repeat sequences flanking nucleosomes in vivo and in vitro. Further, the strong sequence variations in the plasmodium histones, compared to other mammalian histones, do not present adaptations to its AT-rich genome. Human and parasite histones bind with higher affinity to GC-rich DNA and with lower affinity to AT-rich DNA. However, the plasmodium nucleosomes are overall less stable, with increased temperature induced mobility, decreased salt stability of the histones H2A and H2B and considerable reduced binding affinity to GC-rich DNA, as compared with the human nucleosomes. In addition, we show that plasmodium histone octamers form the shortest known nucleosome repeat length (155bp) in vitro and in vivo. Our data suggest that the biochemical properties of the parasite histones are distinct from the typical characteristics of other eukaryotic histones and these properties reflect the increased accessibility of the P. falciparum genome.

Nucleosomes are not positioned randomly on DNA but on preferential sites with respect to the underlying DNA sequence. Histones belong to the most conserved eukaryotic proteins, as sequence dependent nucleosome positioning is an essential regulatory feature of nucleosomes, determining the accessibility of regulatory factors to DNA. We determined the biochemical properties of plasmodium histones and show that they are distinct from human forms, explaining the accessible chromatin structure of P. falciparum. Amino acid exchanges in the histones do not present an adaption to the AT-rich genome, but rather reduce the binding affinity to GC-rich DNA sequences, resulting in rather unstable nucleosomes with labile H2A and H2B, requiring extra-nucleosomal positioning signals to keep them on place. Plasmodium chromatin exhibits the shortest nucleosome spacing known to date potentially inhibiting the formation of higher order structures and maintaining chromatin accessible.

Comparative Bioinformatics Analysis of Transcription Factor Genes Indicates Conservation of Key Regulatory Domains among Babesia bovis, Babesia microti, and Theileria equi

Apicomplexa tick-borne hemoparasites, including Babesia bovis, Babesia microti, and Theileria equi are responsible for bovine and human babesiosis and equine theileriosis, respectively. These parasites of vast medical, epidemiological, and economic impact have complex life cycles in their vertebrate and tick hosts. Large gaps in knowledge concerning the mechanisms used by these parasites for gene regulation remain. Regulatory genes coding for DNA binding proteins such as members of the Api-AP2, HMG, and Myb families are known to play crucial roles as transcription factors. Although the repertoire of Api-AP2 has been defined and a HMG gene was previously identified in the B. bovis genome, these regulatory genes have not been described in detail in B. microti and T. equi. In this study, comparative bioinformatics was used to: (i) identify and map genes encoding for these transcription factors among three parasites' genomes; (ii) identify a previously unreported HMG gene in B. microti; (iii) define a repertoire of eight conserved Myb genes; and (iv) identify AP2 correlates among B. bovis and the better-studied Plasmodium parasites. Searching the available transcriptome of B. bovis defined patterns of transcription of these three gene families in B. bovis erythrocyte stage parasites. Sequence comparisons show conservation of functional domains and general architecture in the AP2, Myb, and HMG proteins, which may be significant for the regulation of common critical parasite life cycle transitions in B. bovis, B. microti, and T. equi. A detailed understanding of the role of gene families encoding DNA binding proteins will provide new tools for unraveling regulatory mechanisms involved in B. bovis, B. microti, and T. equi life cycles and environmental adaptive responses and potentially contributes to the development of novel convergent strategies for improved control of babesiosis and equine piroplasmosis.

Nucleotide diversity analysis highlights functionally important genomic regions

We analyzed functionality and relative distribution of genetic variants across the complete Oryza sativa genome, using the 40 million single nucleotide polymorphisms (SNPs) dataset from the 3,000 Rice Genomes Project (http://snp-seek.irri.org), the largest and highest density SNP collection for any higher plant. We have shown that the DNA-binding transcription factors (TFs) are the most conserved group of genes, whereas kinases and membrane-localized transporters are the most variable ones. TFs may be conserved because they belong to some of the most connected regulatory hubs that modulate transcription of vast downstream gene networks, whereas signaling kinases and transporters need to adapt rapidly to changing environmental conditions. In general, the observed profound patterns of nucleotide variability reveal functionally important genomic regions. As expected, nucleotide diversity is much higher in intergenic regions than within gene bodies (regions spanning gene models), and protein-coding sequences are more conserved than untranslated gene regions. We have observed a sharp decline in nucleotide diversity that begins at about 250 nucleotides upstream of the transcription start and reaches minimal diversity exactly at the transcription start. We found the transcription termination sites to have remarkably symmetrical patterns of SNP density, implying presence of functional sites near transcription termination. Also, nucleotide diversity was significantly lower near 3′ UTRs, the area rich with regulatory regions.

Recent advances in malaria genomics and epigenomics

Malaria continues to impose a significant disease burden on low- and middle-income countries in the tropics. However, revolutionary progress over the last 3 years in nucleic acid sequencing, reverse genetics, and post-genome analyses has generated step changes in our understanding of malaria parasite (Plasmodium spp.) biology and its interactions with its host and vector. Driven by the availability of vast amounts of genome sequence data from Plasmodium species strains, relevant human populations of different ethnicities, and mosquito vectors, researchers can consider any biological component of the malarial process in isolation or in the interactive setting that is infection. In particular, considerable progress has been made in the area of population genomics, with Plasmodium falciparum serving as a highly relevant model. Such studies have demonstrated that genome evolution under strong selective pressure can be detected. These data, combined with reverse genetics, have enabled the identification of the region of the P. falciparum genome that is under selective pressure and the confirmation of the functionality of the mutations in the kelch13 gene that accompany resistance to the major frontline antimalarial, artemisinin. Furthermore, the central role of epigenetic regulation of gene expression and antigenic variation and developmental fate in P. falciparum is becoming ever clearer. This review summarizes recent exciting discoveries that genome technologies have enabled in malaria research and highlights some of their applications to healthcare. The knowledge gained will help to develop surveillance approaches for the emergence or spread of drug resistance and to identify new targets for the development of antimalarial drugs and perhaps vaccines.

H3.3 demarcates GC-rich coding and subtelomeric regions and serves as potential memory mark for virulence gene expression in Plasmodium falciparum

Histones, by packaging and organizing the DNA into chromatin, serve as essential building blocks for eukaryotic life. The basic structure of the chromatin is established by four canonical histones (H2A, H2B, H3 and H4), while histone variants are more commonly utilized to alter the properties of specific chromatin domains. H3.3, a variant of histone H3, was found to have diverse localization patterns and functions across species but has been rather poorly studied in protists. Here we present the first genome-wide analysis of H3.3 in the malaria-causing, apicomplexan parasite, P. falciparum, which revealed a complex occupancy profile consisting of conserved and parasite-specific features. In contrast to other histone variants, PfH3.3 primarily demarcates euchromatic coding and subtelomeric repetitive sequences. Stable occupancy of PfH3.3 in these regions is largely uncoupled from the transcriptional activity and appears to be primarily dependent on the GC-content of the underlying DNA. Importantly, PfH3.3 specifically marks the promoter region of an active and poised, but not inactive antigenic variation (var) gene, thereby potentially contributing to immune evasion. Collectively, our data suggest that PfH3.3, together with other histone variants, indexes the P. falciparum genome to functionally distinct domains and contribute to a key survival strategy of this deadly pathogen.

The nuclear envelope and gene organization in parasitic protozoa: Specializations associated with disease

The parasitic protozoa Trypanosoma brucei and Plasmodium falciparum are lethal human parasites that have developed elegant strategies of immune evasion by antigenic variation. Despite the vast evolutionary distance between the two taxa, both parasites employ strict monoallelic expression of their membrane proteins, variant surface glycoproteins in Trypanosomes and the var, rif and stevor genes in Plasmodium, in order to evade their host's immune system. Additionally, both telomeric location and epigenetic controls are prominent features of these membrane proteins. As such, telomeres, chromatin structure and nuclear organization all contribute to control of gene expression and immune evasion. Here, we discuss the importance of epigenetics and sub-nuclear context for the survival of these disease-causing parasites.

Integrated transcriptomic and proteomic analyses of P. falciparum gametocytes: molecular insight into sex-specific processes and translational repression

Sexual differentiation of malaria parasites into gametocytes in the vertebrate host and subsequent gamete fertilization in mosquitoes is essential for the spreading of the disease. The molecular processes orchestrating these transitions are far from fully understood. Here, we report the first transcriptome analysis of male and female Plasmodium falciparum gametocytes coupled with a comprehensive proteome analysis. In male gametocytes there is an enrichment of proteins involved in the formation of flagellated gametes; proteins involved in DNA replication, chromatin organization and axoneme formation. On the other hand, female gametocytes are enriched in proteins required for zygote formation and functions after fertilization; protein-, lipid- and energy-metabolism. Integration of transcriptome and proteome data revealed 512 highly expressed maternal transcripts without corresponding protein expression indicating large scale translational repression in P. falciparum female gametocytes for the first time. Despite a high degree of conservation between Plasmodium species, 260 of these ‘repressed transcripts’ have not been previously described. Moreover, for some of these genes, protein expression is only reported in oocysts and sporozoites indicating that repressed transcripts can be partitioned into short- and long-term storage. Finally, these data sets provide an essential resource for identification of vaccine/drug targets and for further mechanistic studies.