Evelien M. Bunnik

Three-dimensional modeling of the P. falciparum genome during the erythrocytic cycle reveals a strong connection between genome architecture and gene expression

The development of the human malaria parasite Plasmodium falciparum is controlled by coordinated changes in gene expression throughout its complex life cycle, but the corresponding regulatory mechanisms are incompletely understood. To study the relationship between genome architecture and gene regulation in Plasmodium, we assayed the genome architecture of P. falciparum at three time points during its erythrocytic (asexual) cycle. Using chromosome conformation capture coupled with next-generation sequencing technology (Hi-C), we obtained high-resolution chromosomal contact maps, which we then used to construct a consensus three-dimensional genome structure for each time point. We observed strong clustering of centromeres, telomeres, ribosomal DNA, and virulence genes, resulting in a complex architecture that cannot be explained by a simple volume exclusion model. Internal virulence gene clusters exhibit domain-like structures in contact maps, suggesting that they play an important role in the genome architecture. Midway during the erythrocytic cycle, at the highly transcriptionally active trophozoite stage, the genome adopts a more open chromatin structure with increased chromosomal intermingling. In addition, we observed reduced expression of genes located in spatial proximity to the repressive subtelomeric center, and colocalization of distinct groups of parasite-specific genes with coordinated expression profiles. Overall, our results are indicative of a strong association between the P. falciparum spatial genome organization and gene expression. Understanding the molecular processes involved in genome conformation dynamics could contribute to the discovery of novel antimalarial strategies.

Polysome profiling reveals translational control of gene expression in the human malaria parasite Plasmodium falciparum

BackgroundIn eukaryotic organisms, gene expression is regulated at multiple levels during the processes of transcription and translation. The absence of a tight regulatory network for transcription in the human malaria parasite suggests that gene expression may largely be controlled at post-transcriptional and translational levels.ResultsIn this study, we compare steady-state mRNA and polysome-associated mRNA levels of Plasmodium falciparum at different time points during its asexual cell cycle. For more than 30% of its genes, we observe a delay in peak transcript abundance in the polysomal fraction as compared to the steady-state mRNA fraction, suggestive of strong translational control. Our data show that key regulatory mechanisms could include inhibitory activity of upstream open reading frames and translational repression of the major virulence gene family by intronic transcripts. In addition, we observe polysomal mRNA-specific alternative splicing events and widespread transcription of non-coding transcripts.ConclusionsThese different layers of translational regulation are likely to contribute to a complex network that controls gene expression in this eukaryotic pathogen. Disrupting the mechanisms involved in such translational control could provide novel anti-malarial strategies.

The multifunctional autophagy pathway in the human malaria parasite, Plasmodium falciparum

Autophagy is a catabolic pathway typically induced by nutrient starvation to recycle amino acids, but can also function in removing damaged organelles. In addition, this pathway plays a key role in eukaryotic development. To date, not much is known about the role of autophagy in apicomplexan parasites and more specifically in the human malaria parasite Plasmodium falciparum. Comparative genomic analysis has uncovered some, but not all, orthologs of autophagy-related (ATG) genes in the malaria parasite genome. Here, using a genome-wide in silico analysis, we confirmed that ATG genes whose products are required for vesicle expansion and completion are present, while genes involved in induction of autophagy and cargo packaging are mostly absent. We subsequently focused on the molecular and cellular function of P. falciparum ATG8 (PfATG8), an autophagosome membrane marker and key component of the autophagy pathway, throughout the parasite asexual and sexual erythrocytic stages. In this context, we showed that PfATG8 has a distinct and atypical role in parasite development. PfATG8 localized in the apicoplast and in vesicles throughout the cytosol during parasite development. Immunofluorescence assays of PfATG8 in apicoplast-minus parasites suggest that PfATG8 is involved in apicoplast biogenesis. Furthermore, treatment of parasite cultures with bafilomycin A1 and chloroquine, both lysosomotropic agents that inhibit autophagosome and lysosome fusion, resulted in dramatic morphological changes of the apicoplast, and parasite death. Furthermore, deep proteomic analysis of components associated with PfATG8 indicated that it may possibly be involved in ribophagy and piecemeal microautophagy of the nucleus. Collectively, our data revealed the importance and specificity of the autophagy pathway in the malaria parasite and offer potential novel therapeutic strategies.

Reconstructing the Dynamics of HIV Evolution within Hosts from Serial Deep Sequence Data

At the early stage of infection, human immunodeficiency virus (HIV)-1 predominantly uses the CCR5 coreceptor for host cell entry. The subsequent emergence of HIV variants that use the CXCR4 coreceptor in roughly half of all infections is associated with an accelerated decline of CD4+ T-cells and rate of progression to AIDS. The presence of a ‘fitness valley’ separating CCR5- and CXCR4-using genotypes is postulated to be a biological determinant of whether the HIV coreceptor switch occurs. Using phylogenetic methods to reconstruct the evolutionary dynamics of HIV within hosts enables us to discriminate between competing models of this process. We have developed a phylogenetic pipeline for the molecular clock analysis, ancestral reconstruction, and visualization of deep sequence data. These data were generated by next-generation sequencing of HIV RNA extracted from longitudinal serum samples (median 7 time points) from 8 untreated subjects with chronic HIV infections (Amsterdam Cohort Studies on HIV-1 infection and AIDS). We used the known dates of sampling to directly estimate rates of evolution and to map ancestral mutations to a reconstructed timeline in units of days. HIV coreceptor usage was predicted from reconstructed ancestral sequences using the geno2pheno algorithm. We determined that the first mutations contributing to CXCR4 use emerged about 16 (per subject range 4 to 30) months before the earliest predicted CXCR4-using ancestor, which preceded the first positive cell-based assay of CXCR4 usage by 10 (range 5 to 25) months. CXCR4 usage arose in multiple lineages within 5 of 8 subjects, and ancestral lineages following alternate mutational pathways before going extinct were common. We observed highly patient-specific distributions and time-scales of mutation accumulation, implying that the role of a fitness valley is contingent on the genotype of the transmitted variant.

Genome-wide Mapping of DNA Methylation in the Human Malaria Parasite Plasmodium falciparum

Cytosine DNA methylation is an epigenetic mark in most eukaryotic cells that regulates numerous processes, including gene expression and stress responses. We performed a genome-wide analysis of DNA methylation in the human malaria parasite Plasmodium falciparum. We mapped the positions of methylated cytosines and identified a single functional DNA methyltransferase (Plasmodium falciparum DNA methyltransferase; PfDNMT) that may mediate these genomic modifications. These analyses revealed that the malaria genome is asymmetrically methylated and shares common features with undifferentiated plant and mammalian cells. Notably, core promoters are hypomethylated, and transcript levels correlate with intraexonic methylation. Additionally, there are sharp methylation transitions at nucleosome and exon-intron boundaries. These data suggest that DNA methylation could regulate virulence gene expression and transcription elongation. Furthermore, the broad range of action of DNA methylation and the uniqueness of PfDNMT suggest that the methylation pathway is a potential target for antimalarial strategies.

Multiple dimensions of epigenetic gene regulation in the malaria parasite Plasmodium falciparum

Plasmodium falciparum is the most deadly human malarial parasite, responsible for an estimated 207 million cases of disease and 627,000 deaths in 2012. Recent studies reveal that the parasite actively regulates a large fraction of its genes throughout its replicative cycle inside human red blood cells and that epigenetics plays an important role in this precise gene regulation. Here, we discuss recent advances in our understanding of three aspects of epigenetic regulation in P. falciparum: changes in histone modifications, nucleosome occupancy and the three‐dimensional genome structure. We compare these three aspects of the P. falciparum epigenome to those of other eukaryotes, and show that large‐scale compartmentalization is particularly important in determining histone decomposition and gene regulation in P. falciparum. We conclude by presenting a gene regulation model for P. falciparum that combines the described epigenetic factors, and by discussing the implications of this model for the future of malaria research.

The mRNA-bound proteome of the human malaria parasite Plasmodium falciparum

BackgroundGene expression is controlled at multiple levels, including transcription, stability, translation, and degradation. Over the years, it has become apparent that Plasmodium falciparum exerts limited transcriptional control of gene expression, while at least part of Plasmodium’s genome is controlled by post-transcriptional mechanisms. To generate insights into the mechanisms that regulate gene expression at the post-transcriptional level, we undertook complementary computational, comparative genomics, and experimental approaches to identify and characterize mRNA-binding proteins (mRBPs) in P. falciparum.ResultsClose to 1000 RNA-binding proteins are identified by hidden Markov model searches, of which mRBPs encompass a relatively large proportion of the parasite proteome as compared to other eukaryotes. Several abundant mRNA-binding domains are enriched in apicomplexan parasites, while strong depletion of mRNA-binding domains involved in RNA degradation is observed. Next, we experimentally capture 199 proteins that interact with mRNA during the blood stages, 64 of which with high confidence. These captured mRBPs show a significant overlap with the in silico identified candidate RBPs (pu2009<u20090.0001). Among the experimentally validated mRBPs are many known translational regulators active in other stages of the parasite’s life cycle, such as DOZI, CITH, PfCELF2, Musashi, and PfAlba1–4. Finally, we also detect several proteins with an RNA-binding domain abundant in Apicomplexans (RAP domain) that is almost exclusively found in apicomplexan parasites.ConclusionsCollectively, our results provide the most complete comparative genomics and experimental analysis of mRBPs in P. falciparum. A better understanding of these regulatory proteins will not only give insight into the intricate parasite life cycle but may also provide targets for novel therapeutic strategies.

DNA-encoded nucleosome occupancy is associated with transcription levels in the human malaria parasite Plasmodium falciparum

BackgroundIn eukaryotic organisms, packaging of DNA into nucleosomes controls gene expression by regulating access of the promoter to transcription factors. The human malaria parasite Plasmodium falciparum encodes relatively few transcription factors, while extensive nucleosome remodeling occurs during its replicative cycle in red blood cells. These observations point towards an important role of the nucleosome landscape in regulating gene expression. However, the relation between nucleosome positioning and transcriptional activity has thus far not been explored in detail in the parasite.ResultsHere, we analyzed nucleosome positioning in the asexual and sexual stages of the parasite’s erythrocytic cycle using chromatin immunoprecipitation of MNase-digested chromatin, followed by next-generation sequencing. We observed a relatively open chromatin structure at the trophozoite and gametocyte stages, consistent with high levels of transcriptional activity in these stages. Nucleosome occupancy of genes and promoter regions were subsequently compared to steady-state mRNA expression levels. Transcript abundance showed a strong inverse correlation with nucleosome occupancy levels in promoter regions. In addition, AT-repeat sequences were strongly unfavorable for nucleosome binding in P. falciparum, and were overrepresented in promoters of highly expressed genes.ConclusionsThe connection between chromatin structure and gene expression in P. falciparum shares similarities with other eukaryotes. However, the remarkable nucleosome dynamics during the erythrocytic stages and the absence of a large variety of transcription factors may indicate that nucleosome binding and remodeling are critical regulators of transcript levels. Moreover, the strong dependency between chromatin structure and DNA sequence suggests that the P. falciparum genome may have been shaped by nucleosome binding preferences. Nucleosome remodeling mechanisms in this deadly parasite could thus provide potent novel anti-malarial targets.

Dynamic and Combinatorial Landscape of Histone Modifications during the Intraerythrocytic Developmental Cycle of the Malaria Parasite

A major obstacle in understanding the complex biology of the malaria parasite remains to discover how gene transcription is controlled during its life cycle. Accumulating evidence indicates that the parasites epigenetic state plays a fundamental role in gene expression and virulence. Using a comprehensive and quantitative mass spectrometry approach, we determined the global and dynamic abundance of histones and their covalent post-transcriptional modifications throughout the intraerythrocytic developmental cycle of Plasmodium falciparum. We detected a total of 232 distinct modifications, of which 160 had never been detected in Plasmodium and 88 had never been identified in any other species. We further validated over 10% of the detected modifications and their expression patterns by multiple reaction monitoring assays. In addition, we uncovered an unusual chromatin organization with parasite-specific histone modifications and combinatorial dynamics that may be directly related to transcriptional activity, DNA replication, and cell cycle progression. Overall, our data suggest that the malaria parasite has a unique histone modification signature that correlates with parasite virulence.

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

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