Chromatin associated sense and antisense noncoding RNAs are transcribed from the var gene family of virulence genes of the malaria parasite Plasmodium falciparum

A comprehensive epigenome map of Plasmodium falciparum reveals unique mechanisms of transcriptional regulation and identifies H3K36me2 as a global mark of gene suppression

Background

Role of epigenetic mechanisms towards regulation of the complex life cycle/pathogenesis of Plasmodium falciparum, the causative agent of malaria, has been poorly understood. To elucidate stage-specific epigenetic regulation, we performed genome-wide mapping of multiple histone modifications of P. falciparum. Further to understand the differences in transcription regulation in P. falciparum and its host, human, we compared their histone modification profiles.

Results

Our comprehensive comparative analysis suggests distinct mode of transcriptional regulation in malaria parasite by virtue of poised genes and differential histone modifications. Furthermore, analysis of histone modification profiles predicted 562 genes producing anti-sense RNAs and 335 genes having bidirectional promoter activity, which raises the intriguing possibility of RNA-mediated regulation of transcription in P. falciparum. Interestingly, we found that H3K36me2 acts as a global repressive mark and gene regulation is fine tuned by the ratio of activation marks to H3K36me2 in P. falciparum. This novel mechanism of gene regulation is supported by the fact that knockout of SET genes (responsible for H3K36 methylation) leads to up-regulation of genes with highest occupancy of H3K36me2 in wild-type P. falciparum. Moreover, virulence (var) genes are mostly poised and marked by a unique set of activation (H4ac) and repression (H3K9me3) marks, which are mutually exclusive to other Plasmodium housekeeping genes.

Conclusions

Our study reveals unique plasticity in the epigenetic regulation in P. falciparum which can influence parasite virulence and pathogenicity. The observed differences in the histone code and transcriptional regulation in P. falciparum and its host will open new avenues for epigenetic drug development against malaria parasite.

Electronic supplementary material

The online version of this article (doi:10.1186/s13072-015-0029-1) contains supplementary material, which is available to authorized users.

Dual regulatory effects of non-coding GC-rich elements on the expression of virulence genes in malaria parasites

As the primary virulence factor of falciparum malaria, var genes harboring mutually exclusive expression pattern lead to antigenic variation and immune evasion of this pathogen in human host. Although various mechanisms contribute to silence of var genes, little is known of transcriptional activation pathways of a single var gene and maintenance of its active state with other silent var loci. Here, we report a monoallelic expression pattern of the non-coding GC-elements flanking chromosomal internal var genes, and transcript from the active one was required for activation of the var gene in the same array. Meanwhile, GFP reporter assays revealed a repressive effect on the adjacent gene induced by DNA motifs of the insulator-like GC-element, which was linked to heterochromatin subnuclear localization. Taken together, these data for the first time provide experimental evidence of the dual cis- and trans-acting regulatory functions of the GC-elements in both silence and activation of var genes, which would advance our understanding of the complex regulatory network of the virulence gene family in P. falciparum.

Strand-specific RNA sequencing in Plasmodium falciparum malaria identifies developmentally regulated long non-coding RNA and circular RNA

Background

The human malaria parasite Plasmodium falciparum has a complex and multi-stage life cycle that requires extensive and precise gene regulation to allow invasion and hijacking of host cells, transmission, and immune escape. To date, the regulatory elements orchestrating these critical parasite processes remain largely unknown. Yet it is becoming increasingly clear that long non-coding RNAs (lncRNAs) could represent a missing regulatory layer across a broad range of organisms.

Results

To investigate the regulatory capacity of lncRNA in P. falciparum, we harvested fifteen samples from two time-courses. Our sample set profiled 56 h of P. falciparum blood stage development. We then developed and validated strand-specific, non-polyA-selected RNA sequencing methods, and pursued the first assembly of P. falciparum strand-specific transcript structures from RNA sequencing data. This approach enabled the annotation of over one thousand lncRNA transcript models and their comprehensive global analysis: coding prediction, periodicity, stage-specificity, correlation, GC content, length, location relative to annotated transcripts, and splicing. We validated the complete splicing structure of three lncRNAs with compelling properties. Non-polyA-selected deep sequencing also enabled the prediction of hundreds of intriguing P. falciparum circular RNAs, six of which we validated experimentally.

Conclusions

We found that a subset of lncRNAs, including all subtelomeric lncRNAs, strongly peaked in expression during invasion. By contrast, antisense transcript levels significantly dropped during invasion. As compared to neighboring mRNAs, the expression of antisense-sense pairs was significantly anti-correlated during blood stage development, indicating transcriptional interference. We also validated that P. falciparum produces circRNAs, which is notable given the lack of RNA interference in the organism, and discovered that a highly expressed, five-exon antisense RNA is poised to regulate P. falciparum gametocyte development 1 (PfGDV1), a gene required for early sexual commitment events.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1603-4) contains supplementary material, which is available to authorized users.

Antisense long noncoding RNAs regulate var gene activation in the malaria parasite Plasmodium falciparum

The virulence of Plasmodium falciparum, the causative agent of the deadliest form of human malaria, is attributed to its ability to evade human immunity through antigenic variation. These parasites alternate between expression of variable antigens, encoded by members of a multicopy gene family named var. Immune evasion through antigenic variation depends on tight regulation of var gene expression, ensuring that only a single var gene is expressed at a time while the rest of the family is maintained transcriptionally silent. Understanding how a single gene is chosen for activation is critical for understanding mutually exclusive expression but remains a mystery. Here, we show that antisense long noncoding RNAs (lncRNAs) initiating from var introns are associated with the single active var gene at the time in the cell cycle when the single var upstream promoter is active. We demonstrate that these antisense transcripts are incorporated into chromatin, and that expression of these antisense lncRNAs in trans triggers activation of a silent var gene in a sequence- and dose-dependent manner. On the other hand, interference with these lncRNAs using complement peptide nucleic acid molecules down-regulated the active var gene, erased the epigenetic memory, and induced expression switching. Altogether, our data provide evidence that these antisense lncRNAs play a key role in regulating var gene activation and mutually exclusive expression.

The relative ages of eukaryotes and akaryotes

The Last Eukaryote Common Ancestor (LECA) appears to have the genetics required for meiosis, mitosis, nucleus and nuclear substructures, an exon/intron gene structure, spliceosomes, many centres of DNA replication, etc. (and including mitochondria). Most of these features are not generally explained by models for the origin of the Eukaryotic cell based on the fusion of an Archeon and a Bacterium. We find that the term 'prokaryote' is ambiguous and the non-phylogenetic term akaryote should be used in its place because we do not yet know the direction of evolution between eukaryotes and akaryotes. We use the term 'protoeukaryote' for the hypothetical stem group ancestral eukaryote that took up a bacterium as an endosymbiont that formed the mitochondrion. It is easier to make detailed models with a eukaryote to an akaryote transition, rather than vice versa. So we really are at a phylogenetic impasse in not being confident about the direction of change between eukaryotes and akaryotes.

Noncoding RNAs as emerging regulators of Plasmodium falciparum virulence gene expression

The eukaryotic unicellular pathogen Plasmodium falciparum tightly regulates gene expression, both during development and in adaptation to dynamic host environments. This regulation is evident in the mutually exclusive expression of members of clonally variant virulence multigene families. While epigenetic regulators have been selectively identified at active or repressed virulence genes, their specific recruitment remains a mystery. In recent years, noncoding RNAs (ncRNAs) have emerged as lynchpins of eukaryotic gene regulation; by binding to epigenetic regulators, they provide target specificity to otherwise non-specific enzyme complexes. Not surprisingly, there is great interest in understanding the role of ncRNA in P. falciparum, in particular, their contribution to the mutually exclusive expression of virulence genes. The current repertoire of P. falciparum ncRNAs includes, but is not limited to, subtelomeric ncRNAs, virulence gene-associated ncRNAs and natural antisense RNA transcripts. Continued improvement in high-throughput sequencing methods is sure to expand this repertoire. Here, we summarize recent advances in P. falciparum ncRNA biology, with an emphasis on ncRNA-mediated epigenetic modes of gene regulation.

Exonuclease-mediated degradation of nascent RNA silences genes linked to severe malaria

Antigenic variation of the Plasmodium falciparum multicopy var gene family enables parasite evasion of immune destruction by host antibodies. Expression of a particular var subgroup, termed upsA, is linked to the obstruction of blood vessels in the brain and to the pathogenesis of human cerebral malaria. The mechanism determining upsA activation remains unknown. Here we show that an entirely new type of gene silencing mechanism involving an exonuclease-mediated degradation of nascent RNA controls the silencing of genes linked to severe malaria. We identify a novel chromatin-associated exoribonuclease, termed PfRNase II, that controls the silencing of upsA var genes by marking their transcription start site and intron-promoter regions leading to short-lived cryptic RNA. Parasites carrying a deficient PfRNase II gene produce full-length upsA var transcripts and intron-derived antisense long non-coding RNA. The presence of stable upsA var transcripts overcomes monoallelic expression, resulting in the simultaneous expression of both upsA and upsC type PfEMP1 proteins on the surface of individual infected red blood cells. In addition, we observe an inverse relationship between transcript levels of PfRNase II and upsA-type var genes in parasites from severe malaria patients, implying a crucial role of PfRNase II in severe malaria. Our results uncover a previously unknown type of post-transcriptional gene silencing mechanism in malaria parasites with repercussions for other organisms. Additionally, the identification of RNase II as a parasite protein controlling the expression of virulence genes involved in pathogenesis in patients with severe malaria may provide new strategies for reducing malaria mortality.

A var gene upstream element controls protein synthesis at the level of translation initiation in Plasmodium falciparum

Clonally variant protein expression in the malaria parasite Plasmodium falciparum generates phenotypic variability and allows isogenic populations to adapt to environmental changes encountered during blood stage infection. The underlying regulatory mechanisms are best studied for the major virulence factor P. falciparum erythrocyte membrane protein 1 (PfEMP1). PfEMP1 is encoded by the multicopy var gene family and only a single variant is expressed in individual parasites, a concept known as mutual exclusion or singular gene choice. var gene activation occurs in situ and is achieved through the escape of one locus from epigenetic silencing. Singular gene choice is controlled at the level of transcription initiation and var 5' upstream (ups) sequences harbour regulatory information essential for mutually exclusive transcription as well as for the trans-generational inheritance of the var activity profile. An additional level of control has recently been identified for the var2csa gene, where an mRNA element in the 5' untranslated region (5' UTR) is involved in the reversible inhibition of translation of var2csa transcripts. Here, we extend the knowledge on post-transcriptional var gene regulation to the common upsC type. We identified a 5' UTR sequence that inhibits translation of upsC-derived mRNAs. Importantly, this 5' UTR element efficiently inhibits translation even in the context of a heterologous upstream region. Further, we found var 5' UTRs to be significantly enriched in uAUGs which are known to impair the efficiency of protein translation in other eukaryotes. Our findings suggest that regulation at the post-transcriptional level is a common feature in the control of PfEMP1 expression in P. falciparum.

Epigenetic memory takes center stage in the survival strategy of malaria parasites

Malaria parasites run through a complex life cycle in the vertebrate host and mosquito vector. This not only requires tightly controlled mechanisms to govern stage-specific gene expression but also necessitates effective strategies for survival under changing environmental conditions. In recent years, the combination of different -omics approaches and targeted functional studies highlighted that Plasmodium falciparum blood stage parasites use heterochromatin-based gene silencing as a unifying strategy for clonally variant expression of hundreds of genes. In this article, we describe the epigenetic control mechanisms that mediate alternative expression states of genes involved in antigenic variation, nutrient uptake and sexual conversion and discuss the relevance of this strategy for the survival and transmission of malaria parasites.

The var3 genes of Plasmodium falciparum 3D7 strain are differentially expressed in infected erythrocytes

Plasmodium falciparum erythrocyte membrane protein 1 (PfEMP1) is an important virulence factor encoded by a family of 59 var genes, including 56 var genes plus 3 small var3 genes. The var genes are among the most diverse sequences in the P. falciparum genome, but the var3 genes are found conserved in most P. falciparum strains. Previous studies have been mainly focused on the typical var genes, while the biological characteristics of the var3 genes remain unknown. In this study, the three var3 genes, PF3D7_0100300, PF3D7_0600400, and PF3D7_0937600, were found to be transcribed in the erythrocytic stages of P. falciparum, with a peak in the transcription level at 16 h post-invasion, but terminated immediately after 16 h post-invasion. The encoded protein of PF3D7_0600400 could be detected in both the late trophozoite stage and schizont stage, while the encoded proteins of PF3D7_0100300 and PF3D7_0937600 could only be detected in the late trophozoite stage and schizont stage, respectively. Thus, the var3 genes of the P. falciparum 3D7 strain were differentially expressed during the erythrocytic development of the parasite.

Strand-specific RNA-Seq reveals widespread and developmentally regulated transcription of natural antisense transcripts in Plasmodium falciparum

BACKGROUND

Advances in high-throughput sequencing have led to the discovery of widespread transcription of natural antisense transcripts (NATs) in a large number of organisms, where these transcripts have been shown to play important roles in the regulation of gene expression. Likewise, the existence of NATs has been observed in Plasmodium but our understanding towards their genome-wide distribution remains incomplete due to the limited depth and uncertainties in the level of strand specificity of previous datasets.

RESULTS

To gain insights into the genome-wide distribution of NATs in P. falciparum, we performed RNA-ligation based strand-specific RNA sequencing at unprecedented depth. Our data indicate that 78.3% of the genome is transcribed during blood-stage development. Moreover, our analysis reveals significant levels of antisense transcription from at least 24% of protein-coding genes and that while expression levels of NATs change during the intraerythrocytic developmental cycle (IDC), they do not correlate with the corresponding mRNA levels. Interestingly, antisense transcription is not evenly distributed across coding regions (CDSs) but strongly clustered towards the 3'-end of CDSs. Furthermore, for a significant subset of NATs, transcript levels correlate with mRNA levels of neighboring genes.Finally, we were able to identify the polyadenylation sites (PASs) for a subset of NATs, demonstrating that at least some NATs are polyadenylated. We also mapped the PASs of 3443 coding genes, yielding an average 3' untranslated region length of 523 bp.

CONCLUSIONS

Our strand-specific analysis of the P. falciparum transcriptome expands and strengthens the existing body of evidence that antisense transcription is a substantial phenomenon in P. falciparum. For a subset of neighboring genes we find that sense and antisense transcript levels are intricately linked while other NATs appear to be regulated independently of mRNA transcription. Our deep strand-specific dataset will provide a valuable resource for the precise determination of expression levels as it separates sense from antisense transcript levels, which we find to often significantly differ. In addition, the extensive novel data on 3' UTR length will allow others to perform searches for regulatory motifs in the UTRs and help understand post-translational regulation in P. falciparum.

Genome-wide identification and functional annotation of Plasmodium falciparum long noncoding RNAs from RNA-seq data

The life cycle of Plasmodium falciparum is very complex, with an erythrocytic stage that involves the invasion of red blood cells and the survival and growth of the parasite within the host. Over the past several decades, numbers of studies have shown that proteins exported by P. falciparum to the surface of infected red blood cells play a critical role in recognition and interaction with host receptors and are thus essential for the completion of the life cycle of P. falciparum. However, little is known about long noncoding RNAs (lncRNAs). In this study, we designed a computational pipeline to identify new lncRNAs of P. falciparum from published RNA-seq data and analyzed their sequences and expression features. As a result, 164 novel lncRNAs were found. The sequences and expression features of P. falciparum lncRNAs were similar to those of humans and mice: there was a lack of sequence conservation, low expression levels, and high expression coefficient of variance and co-expression with nearby coding sequences in the genome. Next, a coding/noncoding gene co-expression network for P. falciparum was constructed to further annotate the functions of novel and known lncRNAs. In total, the functions of 69 lncRNAs, including 44 novel lncRNAs, were annotated. The main functions of the lncRNAs included metabolic processes, biosynthetic processes, regulation of biological processes, establishment of localization, catabolic processes, cellular component organization, and interspecies interactions between organisms. Our results will provide clues to further the investigation of interactions between human hosts and parasites and the mechanisms of P. falciparum infection.

The role played by alternative splicing in antigenic variability in human endo-parasites

Endo-parasites that affect humans include Plasmodium, the causative agent of malaria, which remains one of the leading causes of death in human beings. Despite decades of research, vaccines to this and other endo-parasites remain elusive. This is in part due to the hyper-variability of the parasites surface proteins. Generally these surface proteins are encoded by a large family of genes, with only one being dominantly expressed at certain life stages. Another layer of complexity can be introduced through the alternative splicing of these surface proteins. The resulting isoforms may differ from each other with regard to cell localisation, substrate affinities and functions. They may even differ in structure to the extent that they are no longer recognised by the host’s immune system. In many cases this leads to changes in the N terminus of these proteins. The geographical localisation of endo-parasitic infections around the tropics and the highest incidences of HIV-1 infection in the same areas, adds a further layer of complexity as parasitic infections affect the host immune system resulting in higher HIV infection rates, faster disease progression, and an increase in the severity of infections and complications in HIV diagnosis. This review discusses some examples of parasite surface proteins that are alternatively spliced in trypanosomes, Plasmodium and the parasitic worm Schistosoma as well as what role alternate splicing may play in the interaction between HIV and these endo-parasites.

Recruitment of PfSET2 by RNA polymerase II to variant antigen encoding loci contributes to antigenic variation in P. falciparum

Histone modifications are important regulators of gene expression in all eukaryotes. In Plasmodium falciparum, these epigenetic marks regulate expression of genes involved in several aspects of host-parasite interactions, including antigenic variation. While the identities and genomic positions of many histone modifications have now been cataloged, how they are targeted to defined genomic regions remains poorly understood. For example, how variant antigen encoding loci (var) are targeted for deposition of unique histone marks is a mystery that continues to perplex the field. Here we describe the recruitment of an ortholog of the histone modifier SET2 to var genes through direct interactions with the C-terminal domain (CTD) of RNA polymerase II. In higher eukaryotes, SET2 is a histone methyltransferase recruited by RNA pol II during mRNA transcription; however, the ortholog in P. falciparum (PfSET2) has an atypical architecture and its role in regulating transcription is unknown. Here we show that PfSET2 binds to the unphosphorylated form of the CTD, a property inconsistent with its recruitment during mRNA synthesis. Further, we show that H3K36me3, the epigenetic mark deposited by PfSET2, is enriched at both active and silent var gene loci, providing additional evidence that its recruitment is not associated with mRNA production. Over-expression of a dominant negative form of PfSET2 designed to disrupt binding to RNA pol II induced rapid var gene expression switching, confirming both the importance of PfSET2 in var gene regulation and a role for RNA pol II in its recruitment. RNA pol II is known to transcribe non-coding RNAs from both active and silent var genes, providing a possible mechanism by which it could recruit PfSET2 to var loci. This work unifies previous reports of histone modifications, the production of ncRNAs, and the promoter activity of var introns into a mechanism that contributes to antigenic variation by malaria parasites.

Chemical modifications to histones, the proteins that serve as the primary units of chromatin, often determine whether specific genes are actively transcribed or condensed into transcriptionally silent regions of the genome. In the malaria parasite Plasmodium falciparum, histone modifications have been shown to play a significant role in controlling gene expression involved in many aspects of their lifecycle, including the complex gene expression patterns associated with antigenic variation. The various histone modifications that are found within the parasite's genome have now been extensively cataloged, and the enzymes that are responsible for adding and removing them have been identified. However, how these enzymes are recruited to specific regions of the genome to coordinate gene expression is not understood. In this paper, we provide the first evidence for recruitment of a unique histone methyltransferase to specific regions of the genome through its tethering to RNA polymerase II. We find that disruption of this interaction results in major changes in expression patterns of genes involved in antigenic variation, demonstrating the importance of regulated recruitment of histone modifiers for parasite biology.

Non-coding RNAs in homeostasis, disease and stress responses: an evolutionary perspective

Cells and organisms are subject to challenges and perturbations in their environment and physiology in all stages of life. The molecular response to such changes, including insulting conditions such as pathogen infections, involves coordinated modulation of gene expression programmes and has not only homeostatic but also ecological and evolutionary importance. Although attention has been primarily focused on signalling pathways and protein networks, non-coding RNAs (ncRNAs), which comprise a significant output of the genomes of prokaryotes and especially eukaryotes, are increasingly implicated in the molecular mechanisms of these responses. Long and short ncRNAs not only regulate development and cell physiology, they are also involved in disease states, including cancers, in host-pathogen interactions, and in a variety of stress responses. Indeed, regulatory RNAs are part of genetically encoded response networks and also underpin epigenetic processes, which are emerging as key mechanisms of adaptation and transgenerational inheritance. Here we present the growing evidence that ncRNAs are intrinsically involved in cellular and organismal adaptation processes, in both robustness and protection to stresses, as well as in mechanisms generating evolutionary change.

Epigenetic regulation of the Plasmodium falciparum genome

Recent research has highlighted some unique aspects of chromatin biology in the malaria parasite Plasmodium falciparum. During its erythrocytic lifecycle P. falciparum maintains its genome primarily as unstructured euchromatin. Indeed there is no clear role for chromatin-mediated silencing of the majority of the developmentally expressed genes in P. falciparum. However discontinuous stretches of heterochromatin are critical for variegated expression of contingency genes that mediate key pathogenic processes in malaria. These range from invasion of erythrocytes and antigenic variation to solute transport and growth adaptation in response to environmental changes. Despite lack of structure within euchromatin the nucleus maintains functional compartments that regulate expression of many genes at the nuclear periphery, particularly genes with clonally variant expression. The typical components of the chromatin regulatory machinery are present in P. falciparum; however, some of these appear to have evolved novel species-specific functions, e.g. the dynamic regulation of histone variants at virulence gene promoters. The parasite also appears to have repeatedly acquired chromatin regulatory proteins through lateral transfer from endosymbionts and from the host. P. falciparum chromatin regulators have been successfully targeted with multiple drugs in laboratory studies; hopefully their functional divergence from human counterparts will allow the development of parasite-specific inhibitors.

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

Background

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

Results

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

Conclusions

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

Parasite epigenetics and immune evasion: lessons from budding yeast

The remarkable ability of many parasites to evade host immunity is the key to their success and pervasiveness. The immune evasion is directly linked to the silencing of the members of extended families of genes that encode for major parasite antigens. At any time only one of these genes is active. Infrequent switches to other members of the gene family help the parasites elude the immune system and cause prolonged maladies. For most pathogens, the detailed mechanisms of gene silencing and switching are poorly understood. On the other hand, studies in the budding yeast Saccharomyces cerevisiae have revealed similar mechanisms of gene repression and switching and have provided significant insights into the molecular basis of these phenomena. This information is becoming increasingly relevant to the genetics of the parasites. Here we summarize recent advances in parasite epigenetics and emphasize the similarities between S. cerevisiae and pathogens such as Plasmodium, Trypanosoma, Candida, and Pneumocystis. We also outline current challenges in the control and the treatment of the diseases caused by these parasites and link them to epigenetics and the wealth of knowledge acquired from budding yeast.

Spleen-dependent regulation of antigenic variation in malaria parasites: Plasmodium knowlesi SICAvar expression profiles in splenic and asplenic hosts

Background

Antigenic variation by malaria parasites was first described in Plasmodium knowlesi, which infects humans and macaque monkeys, and subsequently in P. falciparum, the most virulent human parasite. The schizont-infected cell agglutination (SICA) variant proteins encoded by the SICAvar multigene family in P. knowlesi, and Erythrocyte Membrane Protein-1 (EMP-1) antigens encoded by the var multigene family in P. falciparum, are expressed at the surface of infected erythrocytes, are associated with virulence, and serve as determinants of naturally acquired immunity. A parental P. knowlesi clone, Pk1(A+), and a related progeny clone, Pk1(B+)1+, derived by an invivo induced variant antigen switch, were defined by the expression of distinct SICA variant protein doublets of 210/190 and 205/200 kDa, respectively. Passage of SICA[+] infected erythrocytes through splenectomized rhesus monkeys results in the SICA[-] phenotype, defined by the lack of surface expression and agglutination with variant specific antisera.

Principal Findings

We have investigated SICAvar RNA and protein expression in Pk1(A+), Pk1(B+)1+, and SICA[-] parasites. The Pk1(A+) and Pk1(B+)1+ parasites express different distinct SICAvar transcript and protein repertoires. By comparison, SICA[-] parasites are characterized by a vast reduction in SICAvar RNA expression, the lack of full-length SICAvar transcript signals on northern blots, and correspondingly, the absence of any SICA protein detected by mass spectrometry.

Significance

SICA protein expression may be under transcriptional as well as post-transcriptional control, and we show for the first time that the spleen, an organ central to blood-stage immunity in malaria, exerts an influence on these processes. Furthermore, proteomics has enabled the first in-depth characterization of SICA[+] protein phenotypes and we show that the invivo switch from Pk1(A+) to Pk1(B+)1+ parasites resulted in a complete change in SICA profiles. These results emphasize the importance of studying antigenic variation in the context of the host environment.

PfSETvs methylation of histone H3K36 represses virulence genes in Plasmodium falciparum

The variant antigen Plasmodium falciparum erythrocyte membrane protein 1 (PfEMP1), which is expressed on the surface of P. falciparum-infected red blood cells, is a critical virulence factor for malaria. Each parasite has 60 antigenically distinct var genes that each code for a different PfEMP1 protein. During infection the clonal parasite population expresses only one gene at a time before switching to the expression of a new variant antigen as an immune-evasion mechanism to avoid the host antibody response. The mechanism by which 59 of the 60 var genes are silenced remains largely unknown. Here we show that knocking out the P. falciparum variant-silencing SET gene (here termed PfSETvs), which encodes an orthologue of Drosophila melanogaster ASH1 and controls histone H3 lysine 36 trimethylation (H3K36me3) on var genes, results in the transcription of virtually all var genes in the single parasite nuclei and their expression as proteins on the surface of individual infected red blood cells. PfSETvs-dependent H3K36me3 is present along the entire gene body, including the transcription start site, to silence var genes. With low occupancy of PfSETvs at both the transcription start site of var genes and the intronic promoter, expression of var genes coincides with transcription of their corresponding antisense long noncoding RNA. These results uncover a previously unknown role of PfSETvs-dependent H3K36me3 in silencing var genes in P. falciparum that might provide a general mechanism by which orthologues of PfSETvs repress gene expression in other eukaryotes. PfSETvs knockout parasites expressing all PfEMP1 proteins may also be applied to the development of a malaria vaccine.

Silence, activate, poise and switch! Mechanisms of antigenic variation in Plasmodium falciparum

Phenotypic variation in genetically identical malaria parasites is an emerging topic. Although antigenic variation is only part of a more global parasite strategy to create adaptation through epigenetically controlled transcriptional variability, it is the central mechanism enabling immune evasion and promoting pathogenesis. The var gene family is the best-studied example in a wide range of clonally variant gene families in Plasmodium falciparum. It is unique in its strict selection of a single member for activation, a process termed monoallelic expression. The conceptual advances that have emerged from studying var genes show striking common epigenetic features with many other clonally variant gene families or even single-copy genes that show a variegated expression in parasite populations. However, major mechanistic questions, such as the existence of a potential expression site and the identity of transcription factors or genetic elements driving singular gene choice, are still unanswered. In this review we discuss the recent findings in the molecular processes essential for clonal variation, namely silencing, activation, poising and switching. Integrating findings about all clonally variant gene families and other mutually exclusive expression systems will hopefully drive mechanistic understanding of antigenic variation.

H2A.Z and H2B.Z double-variant nucleosomes define intergenic regions and dynamically occupy var gene promoters in the malaria parasite Plasmodium falciparum

Histone variants are important components of eukaryotic chromatin and can alter chromatin structure to confer specialized functions. H2B variant histones are rare in nature but have evolved independently in the phyla Apicomplexa and Trypanasomatida. Here, we investigate the apicomplexan-specific Plasmodium falciparum histone variant Pf H2B.Z and show that within nucleosomes Pf H2B.Z dimerizes with the H2A variant Pf H2A.Z and that Pf H2B.Z and Pf H2A.Z occupancy correlates in the subset of genes examined. These double-variant nucleosomes also carry common markers of euchromatin like H3K4me3 and histone acetylation. Pf H2B.Z levels are elevated in intergenic regions across the genome, except in the var multigene family, where Pf H2A.Z/Pf H2B.Z double-variant nucleosomes are only enriched in the promoter of the single active var copy and this enrichment is developmentally regulated. Importantly, this pattern seems to be specific for var genes and does not apply to other heterochromatic gene families involved in red blood cell invasion which are also subject to clonal expression. Thus, Pf H2A.Z/Pf H2B.Z double-variant nucleosomes appear to have a highly specific function in the regulation of P. falciparum virulence.

Early gametocytes of the malaria parasite Plasmodium falciparum specifically remodel the adhesive properties of infected erythrocyte surface

In Plasmodium falciparum infections the parasite transmission stages, the gametocytes, mature in 10 days sequestered in internal organs. Recent studies suggest that cell mechanical properties rather than adhesive interactions play a role in sequestration during gametocyte maturation. It remains instead obscure how sequestration is established, and how the earliest sexual stages, morphologically similar to asexual trophozoites, modify the infected erythrocytes and their cytoadhesive properties at the onset of gametocytogenesis. Here, purified P. falciparum early gametocytes were used to ultrastructurally and biochemically analyse parasite-induced modifications on the red blood cell surface and to measure their functional consequences on adhesion to human endothelial cells. This work revealed that stage I gametocytes are able to deform the infected erythrocytes like asexual parasites, but do not modify its surface with adhesive 'knob' structures and associated proteins. Reduced levels of the P. falciparum erythrocyte membrane protein 1 (PfEMP1) adhesins are exposed on the red blood cell surface by these parasites, and the expression of the var gene family, which encodes 50-60 variants of PfEMP1, is dramatically downregulated in the transition from asexual development to gametocytogenesis. Cytoadhesion assays show that such gene expression changes and host cell surface modifications functionally result in the inability of stage I gametocytes to bind the host ligands used by the asexual parasite to bind endothelial cells. In conclusion, these results identify specific differences in molecular and cellular mechanisms of host cell remodelling and in adhesive properties, leading to clearly distinct host parasite interplays in the establishment of sequestration of stage I gametocytes and of asexual trophozoites.

Analysis of single-cell gene transcription by RNA fluorescent in situ hybridization (FISH)

Adhesion of Plasmodium falciparum infected erythrocytes (IE) to human endothelial receptors during malaria infections is mediated by expression of PfEMP1 protein variants encoded by the var genes. The haploid P. falciparum genome harbors approximately 60 different var genes of which only one has been believed to be transcribed per cell at a time during the blood stage of the infection. How such mutually exclusive regulation of var gene transcription is achieved is unclear, as is the identification of individual var genes or sub-groups of var genes associated with different receptors and the consequence of differential binding on the clinical outcome of P. falciparum infections. Recently, the mutually exclusive transcription paradigm has been called into doubt by transcription assays based on individual P. falciparum transcript identification in single infected erythrocytic cells using RNA fluorescent in situ hybridization (FISH) analysis of var gene transcription by the parasite in individual nuclei of P. falciparum IE(1). Here, we present a detailed protocol for carrying out the RNA-FISH methodology for analysis of var gene transcription in single-nuclei of P. falciparum infected human erythrocytes. The method is based on the use of digoxigenin- and biotin- labeled antisense RNA probes using the TSA Plus Fluorescence Palette System(2) (Perkin Elmer), microscopic analyses and freshly selected P. falciparum IE. The in situ hybridization method can be used to monitor transcription and regulation of a variety of genes expressed during the different stages of the P. falciparum life cycle and is adaptable to other malaria parasite species and other organisms and cell types.

Two long non-coding RNAs generated from subtelomeric regions accumulate in a novel perinuclear compartment in Plasmodium falciparum

Chromosome ends have been implicated in the default silencing of clonally variant gene families in the human malaria parasite Plasmodium falciparum. These chromosome regions are organized into heterochromatin, as defined by the presence of a repressive histone H3 lysine 9 trimethylated marker and heterochromatin protein 1. Here, we show that the non-coding subtelomeric region adjacent to virulence genes forms facultative heterochromatin in a cell cycle-dependent manner. We demonstrate that telomere-associated repeat elements (TAREs) and telomeres are transcribed as long non-coding RNAs (lncRNAs) during schizogony. Northern blot assays revealed two classes of lncRNAs: a ~4-kb transcript composed of telomere sequences and a TARE-3 element, and a >6-kb transcript composed of 21-bp repeats from TARE-6. These lncRNAs are transcribed by RNA polymerase II as single-stranded molecules. RNA-FISH analysis showed that these lncRNAs form several nuclear foci during the schizont stage, whereas in the ring stage, they are located in a single perinuclear compartment that does not co-localize with any known nuclear subcompartment. Furthermore, the TARE-6 lncRNA is predicted to form a stable and repetitive hairpin structure that is able to bind histones. Consequently, the characterization of the molecular interactions of these lncRNAs with nuclear proteins may reveal novel modes of gene regulation and nuclear function in P. falciparum.

Malaria var gene expression: keeping up with the neighbors

Plasmodium falciparum PfEMP1 is a malaria virulence protein whose expression is epigenetically regulated. The parasite's ability to express exclusively only one of the sixty var genes that encode PfEMP1 is essential for disease pathogenesis. Two recent papers identify key molecular players in determining whether a var gene is active or silenced (Volz et al., 2012; Zhang et al., 2011).

Evidence for in vitro and in vivo expression of the conserved VAR3 (type 3) plasmodium falciparum erythrocyte membrane protein 1

BACKGROUND

Members of the Plasmodium falciparum erythrocyte membrane protein 1 (PfEMP1) adhesion antigen family are major contributors to the pathogenesis of P. falciparum malaria infections. The PfEMP1-encoding var genes are among the most diverse sequences in nature, but three genes, var1, var2csa and var3 are found conserved in most parasite genomes. The most severe forms of malaria disease are caused by parasites expressing a subset of antigenically conserved PfEMP1 variants. Thus the ubiquitous and conserved VAR3 PfEMP1 is of particular interest to the research field. Evidence of VAR3 expression on the infected erythrocyte surface has never been presented, and var3 genes have been proposed to be transcribed and expressed differently from the rest of the var gene family members.

METHODS

In this study, parasites expressing VAR3 PfEMP1 were generated using anti-VAR3 antibodies and the var transcript and PfEMP1 expression profiles of the generated parasites were investigated. The IgG reactivity by plasma from children living in malaria-endemic Tanzania was tested to parasites and recombinant VAR3 protein. Parasites from hospitalized children were isolated and the transcript level of var3 was investigated.

RESULTS

Var3 is transcribed and its protein product expressed on the surface of infected erythrocytes. The VAR3-expressing parasites were better recognized by children´s IgG than a parasite line expressing a Group B var gene. Two in 130 children showed increased recognition of parasites expressing VAR3 and to the recombinant VAR3 protein after a malaria episode and the isolated parasites showed high levels of var3 transcripts.

CONCLUSIONS

Collectively, the presented data suggest that var3 is transcribed and its protein product expressed on the surface of infected erythrocytes in the same manner as seen for other var genes both in vitro and in vivo. Only very few children exhibit seroconversion to VAR3 following a malaria episode requiring hospitalization, supporting the previous conclusion drawn from var3 transcript analysis of parasites collected from children hospitalized with malaria, that VAR3 is not associated with severe anaemia or cerebral malaria syndromes in children.

Chromatin modifications, epigenetics, and how protozoan parasites regulate their lives

Chromatin structure plays a vital role in epigenetic regulation of protozoan parasite gene expression. Epigenetic gene regulation impacts upon parasite virulence, differentiation and cell-cycle control. Recent work in many laboratories has elucidated the functions of proteins that regulate parasite gene expression by chemical modification of constituent nucleosomes. A major focus of investigation has been the characterization of post-translational modifications (PTMs) of histones and the identification of the enzymes responsible. Despite conserved features and specificity common to all eukaryotes, parasite enzymes involved in chromatin modification have unique functions that regulate unique aspects of parasite biology.

Characterization of the unusual bidirectional ves promoters driving VESA1 expression and associated with antigenic variation in Babesia bovis

Rapid clonal antigenic variation in Babesia bovis involves the variant erythrocyte surface antigen-1 (VESA1) protein expressed on the infected-erythrocyte surface. Because of the significance of this heterodimeric protein for demonstrated mechanisms of parasite survival and virulence, there is a need to understand how expression of the ves multigene family encoding this protein is controlled. As an initial step toward this goal, we present here initial characterization of the ves promoter driving transcription of VESA1a and -1b subunits. A series of transfection constructs containing various sequence elements from the in vivo locus of active ves transcription (LAT) were used to drive expression of the firefly luciferase gene in a dual luciferase-normalized assay. The results of this approach reveal the presence of two bidirectional promoter activities within the 434-bp intergenic region (IGr), influenced by putative regulatory sequences embedded within the flanking ves1α and ves1β genes. Repressor-like effects on the apposing gene were observed for intron 1 of both ves1α and ves1β. This effect is apparently not dependent upon intronic promoter activity and acts only in cis. The expression of genes within the ves family is likely modulated by local elements embedded within ves coding sequences outside the intergenic promoter region in concert with chromatin modifications. These results provide a framework to help us begin to understand gene regulation during antigenic variation in B. bovis.

PfAlbas constitute a new eukaryotic DNA/RNA-binding protein family in malaria parasites

In Plasmodium falciparum, perinuclear subtelomeric chromatin conveys monoallelic expression of virulence genes. However, proteins that directly bind to chromosome ends are poorly described. Here we identify a novel DNA/RNA-binding protein family that bears homology to the archaeal protein Alba (Acetylation lowers binding affinity). We isolated three of the four PfAlba paralogs as part of a molecular complex that is associated with the P. falciparum-specific TARE6 (Telomere-Associated Repetitive Elements 6) subtelomeric region and showed in electromobility shift assays (EMSAs) that the PfAlbas bind to TARE6 repeats. In early blood stages, the PfAlba proteins were enriched at the nuclear periphery and partially co-localized with PfSir2, a TARE6-associated histone deacetylase linked to the process of antigenic variation. The nuclear location changed at the onset of parasite proliferation (trophozoite-schizont), where the PfAlba proteins were also detectable in the cytoplasm in a punctate pattern. Using single-stranded RNA (ssRNA) probes in EMSAs, we found that PfAlbas bind to ssRNA, albeit with different binding preferences. We demonstrate for the first time in eukaryotes that Alba-like proteins bind to both DNA and RNA and that their intracellular location is developmentally regulated. Discovery of the PfAlbas may provide a link between the previously described subtelomeric non-coding RNA and the regulation of antigenic variation.

In Silico Prediction of Evolutionarily Conserved GC-Rich Elements Associated with Antigenic Proteins of Plasmodium falciparum

The Plasmodium falciparum genome being AT-rich, the presence of GC-rich regions suggests functional significance. Evolution imposes selection pressure to retain functionally important coding and regulatory elements. Hence searching for evolutionarily conserved GC-rich, intergenic regions in an AT-rich genome will help in discovering new coding regions and regulatory elements. We have used elevated GC content in intergenic regions coupled with sequence conservation against P. reichenowi, which is evolutionarily closely related to P. falciparum to identify potential sequences of functional importance. Interestingly, ~30% of the GC-rich, conserved sequences were associated with antigenic proteins encoded by var and rifin genes. The majority of sequences identified in the 5′ UTR of var genes are represented by short expressed sequence tags (ESTs) in cDNA libraries signifying that they are transcribed in the parasite. Additionally, 19 sequences were located in the 3′ UTR of rifins and 4 also have overlapping ESTs. Further analysis showed that several sequences associated with var genes have the capacity to encode small peptides. A previous report has shown that upstream peptides can regulate the expression of var genes hence we propose that these conserved GC-rich sequences may play roles in regulation of gene expression.

Genome cartography: charting the apicomplexan genome

Genes reside in particular genomic contexts that can be mapped at many levels. Historically, 'genetic maps' were used primarily to locate genes. Recent technological advances in the determination of genome sequences have made the analysis and comparison of whole genomes possible and increasingly tractable. What do we see if we shift our focus from gene content (the 'inventory' of genes contained within a genome) to the composition and organization of a genome? This review examines what has been learned about the evolution of the apicomplexan genome as well as the significance and impact of genomic location on our understanding of the eukaryotic genome and parasite biology.

A global transcriptional analysis of Plasmodium falciparum malaria reveals a novel family of telomere-associated lncRNAs

BACKGROUND

Mounting evidence suggests a major role for epigenetic feedback in Plasmodium falciparum transcriptional regulation. Long non-coding RNAs (lncRNAs) have recently emerged as a new paradigm in epigenetic remodeling. We therefore set out to investigate putative roles for lncRNAs in P. falciparum transcriptional regulation.

RESULTS

We used a high-resolution DNA tiling microarray to survey transcriptional activity across 22.6% of the P. falciparum strain 3D7 genome. We identified 872 protein-coding genes and 60 putative P. falciparum lncRNAs under developmental regulation during the parasite's pathogenic human blood stage. Further characterization of lncRNA candidates led to the discovery of an intriguing family of lncRNA telomere-associated repetitive element transcripts, termed lncRNA-TARE. We have quantified lncRNA-TARE expression at 15 distinct chromosome ends and mapped putative transcriptional start and termination sites of lncRNA-TARE loci. Remarkably, we observed coordinated and stage-specific expression of lncRNA-TARE on all chromosome ends tested, and two dominant transcripts of approximately 1.5 kb and 3.1 kb transcribed towards the telomere.

CONCLUSIONS

We have characterized a family of 22 telomere-associated lncRNAs in P. falciparum. Homologous lncRNA-TARE loci are coordinately expressed after parasite DNA replication, and are poised to play an important role in P. falciparum telomere maintenance, virulence gene regulation, and potentially other processes of parasite chromosome end biology. Further study of lncRNA-TARE and other promising lncRNA candidates may provide mechanistic insight into P. falciparum transcriptional regulation.

From in vivo to in vitro: dynamic analysis of Plasmodium falciparum var gene expression patterns of patient isolates during adaptation to culture

Plasmodium falciparum erythrocyte membrane protein 1 (PfEMP1), encoded by the var gene family, plays a crucial role in disease virulence through its involvement in binding to various host cellular receptors during infection. Growing evidence suggests that differential expression of the various var subgroups may be involved in parasite virulence. To further explore this issue, we have collected isolates from symptomatic patients in south China-Myanmar border, and characterized their sequence diversity and transcription profiles over time of var gene family, and cytoadherence properties from the time of their initial collection and extending through a two month period of adaptation to culture. Initially, we established a highly diverse, DBLα (4 cysteines) subtype-enriched, but unique local repertoire of var-DBL1α sequences by cDNA cloning and sequencing. Next we observed a rapid transcriptional decline of upsA- and upsB-subtype var genes at ring stage through qRT-PCR assays, and a switching event from initial ICAM-I binding to the CD36-binding activity during the first week of adaptive cultivation in vitro. Moreover, predominant transcription of upsA var genes was observed to be correlated with those isolates that showed a higher parasitemia at the time of collection and the ICAM-1-binding phenotype in culture. Taken together, these data indicate that the initial stage of adaptive process in vitro significantly influences the transcription of virulence-related var subtypes and expression of PfEMP1 variants. Further, the specific upregulation of the upsA var genes is likely linked to the rapid propagation of the parasite during natural infection due to the A-type PfEMP1 variant-mediated growth advantages.

Expression of P. falciparum var genes involves exchange of the histone variant H2A.Z at the promoter

Plasmodium falciparum employs antigenic variation to evade the human immune response by switching the expression of different variant surface antigens encoded by the var gene family. Epigenetic mechanisms including histone modifications and sub-nuclear compartmentalization contribute to transcriptional regulation in the malaria parasite, in particular to control antigenic variation. Another mechanism of epigenetic control is the exchange of canonical histones with alternative variants to generate functionally specialized chromatin domains. Here we demonstrate that the alternative histone PfH2A.Z is associated with the epigenetic regulation of var genes. In many eukaryotic organisms the histone variant H2A.Z mediates an open chromatin structure at promoters and facilitates diverse levels of regulation, including transcriptional activation. Throughout the asexual, intraerythrocytic lifecycle of P. falciparum we found that the P. falciparum ortholog of H2A.Z (PfH2A.Z) colocalizes with histone modifications that are characteristic of transcriptionally-permissive euchromatin, but not with markers of heterochromatin. Consistent with this finding, antibodies to PfH2A.Z co-precipitate the permissive modification H3K4me3. By chromatin-immunoprecipitation we show that PfH2A.Z is enriched in nucleosomes around the transcription start site (TSS) in both transcriptionally active and silent stage-specific genes. In var genes, however, PfH2A.Z is enriched at the TSS only during active transcription in ring stage parasites. Thus, in contrast to other genes, temporal var gene regulation involves histone variant exchange at promoter nucleosomes. Sir2 histone deacetylases are important for var gene silencing and their yeast ortholog antagonises H2A.Z function in subtelomeric yeast genes. In immature P. falciparum parasites lacking Sir2A or Sir2B high var transcription levels correlate with enrichment of PfH2A.Z at the TSS. As Sir2A knock out parasites mature the var genes are silenced, but PfH2A.Z remains enriched at the TSS of var genes; in contrast, PfH2A.Z is lost from the TSS of de-repressed var genes in mature Sir2B knock out parasites. This result indicates that PfH2A.Z occupancy at the active var promoter is antagonized by PfSir2A during the intraerythrocytic life cycle. We conclude that PfH2A.Z contributes to the nucleosome architecture at promoters and is regulated dynamically in active var genes.

Plasmodium falciparum is a protist parasite that causes malaria and kills more than 800,000 people per year. The parasite escapes from the human immune response by antigenic variation through switching between expression of different var genes. These encode different variant antigens that are expressed on the surface of the infected erythrocyte and mediate pathogenic adhesion of the infected erythrocytes to host receptors. Understanding how this process is regulated may lead to the identification of factors that are essential for immune evasion and that could represent novel drug targets. Here, we have identified the parasite's histone variant PfH2A.Z as a novel contributor to the transcriptional regulation of antigenic variation. PfH2A.Z is enriched in the promoter of many genes, but enrichment correlates with gene expression only in var genes. Furthermore we show that PfH2A.Z enrichment in var promoters is antagonised by the var gene silencing factor PfSir2A. These findings further extend our knowledge of the complex mechanisms regulating gene expression in P. falciparum.

Plasmodium falciparum var gene silencing is determined by cis DNA elements that form stable and heritable interactions

Antigenic variation in the human malaria parasite Plasmodium falciparum depends on the transcriptional regulation of the var gene family. In each individual parasite, mRNA is expressed exclusively from 1 var gene out of ∼60, while the rest of the genes are transcriptionally silenced. Both modifications to chromatin structure and DNA regulatory elements associated with each var gene have been implicated in the organization and maintenance of the silent state. Whether silencing is established at the level of entire chromosomal regions via heterochromatin spreading or at the level of individual var promoters through the action of a silencing element within each var intron has been debated. Here, we consider both possibilities, using clonal parasite lines carrying chromosomally integrated transgenes. We confirm a previous finding that the loss of an adjacent var intron results in var promoter activation and further show that transcriptional activation of a var promoter within a cluster does not affect the transcriptional activity of neighboring var promoters. Our results provide more evidence for the hypothesis that var genes are primarily silenced at the level of an individual gene, rather than by heterochromatin spreading. We also tested the intrinsic directionality of an intron's silencing effect on upstream or downstream var promoters. We found that an intron is capable of silencing in either direction and that, once established, a var promoter-intron pair is stably maintained through many generations, suggesting a possible role in epigenetic memory. This study provides insights into the regulation of endogenous var gene clusters.

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

Epigenetic regulatory mechanisms and their enzymes are promising targets for malaria therapeutic intervention; however, the epigenetic component of gene expression in P. falciparum is poorly understood. Dynamic or stable association of epigenetic marks with genomic features provides important clues about their function and helps to understand how histone variants/modifications are used for indexing the Plasmodium epigenome. We describe a novel, linear amplification method for next-generation sequencing (NGS) that allows unbiased analysis of the extremely AT-rich Plasmodium genome. We used this method for high resolution, genome-wide analysis of a histone H2A variant, H2A.Z and two histone H3 marks throughout parasite intraerythrocytic development. Unlike in other organisms, H2A.Z is a constant, ubiquitous feature of euchromatic intergenic regions throughout the intraerythrocytic cycle. The almost perfect colocalisation of H2A.Z with H3K9ac and H3K4me3 suggests that these marks are preferentially deposited on H2A.Z-containing nucleosomes. By performing RNA-seq on 8 time-points, we show that acetylation of H3K9 at promoter regions correlates very well with the transcriptional status whereas H3K4me3 appears to have stage-specific regulation, being low at early stages, peaking at trophozoite stage, but does not closely follow changes in gene expression. Our improved NGS library preparation procedure provides a foundation to exploit the malaria epigenome in detail. Furthermore, our findings place H2A.Z at the cradle of P. falciparum epigenetic regulation by stably defining intergenic regions and providing a platform for dynamic assembly of epigenetic and other transcription related complexes.

Plasmodium falciparum is a unicellular pathogen that is responsible for the most severe form of malaria. Similar to other eukaryotic organisms, its genome is organized into chromosomes by proteins called histones. Modification or replacement of these histones has marked effects on the packaging grade of DNA and instructs the recruitment of protein complexes, thereby regulating essential cellular processes such as gene expression and replication. Here we unveil the genome-wide localization of two histone H3 modifications (K9ac/K4me3) and a histone variant, H2A.Z, during development of the parasite in the human red blood cells. We find that all three epigenetic features are predominantly present in intergenic regions of the P. falciparum genome, suggesting an interconnecting role in regulation of gene expression. H2A.Z levels appear to be largely invariable throughout intraerythrocytic development while placement/removal of the histone marks is dynamic with H3K9ac and H3K4me3 being transcription-coupled and stage-specific, respectively. These observations support a model in which H2A.Z-containing nucleosomes serve to demarcate regulatory regions in the parasite's genome and promote transcription initiation by guiding chromatin modifying and transcription initiating complexes. The findings and methodological developments presented in this paper provide a cornerstone for future epigenome research in eukaryotic pathogens and vital information to understand and to interfere with parasite development and survival.

Var gene expression and human Plasmodium pathogenesis

Plasmodium falciparum is responsible for most of the morbidity and mortality associated with malaria and is unique in its ability to sequester in organ postcapillary venules. Specific host-parasite interactions mediate this phenomenon and the P. falciparum erythrocyte membrane protein 1 is the predominant ligand responsible for adhering to host endothelial receptors. This review focuses on the current knowledge regarding this protein family, evidence for its role in various pathogenic mechanisms and on insights that have been gained in this area from field studies.

Epigenetics in Plasmodium: what do we really know?

In the burgeoning field of Plasmodium gene expression, there are--to borrow some famous words from a former U.S. Secretary of Defense--"known knowns, known unknowns, and unknown unknowns." This is in itself an important achievement, since it is only in the past decade that facts have begun to move from the third category into the first. Nevertheless, much remains in the middle ground of known or suspected "unknowns." It is clear that the malaria parasite controls vital virulence processes such as host cell invasion and cytoadherence at least partly via epigenetic mechanisms, so a proper understanding of epigenetic transcriptional control in this organism should have great clinical relevance. Plasmodium, however, is an obligate intracellular parasite: it operates not in a vacuum but rather in the complicated context of its metazoan hosts. Therefore, as valuable data about the parasite's basic epigenetic machinery begin to emerge, it becomes increasingly important to relate in vitro studies to the situation in vivo. This review will focus upon the challenge of understanding Plasmodium epigenetics in an integrated manner, in the human and insect hosts as well as the petri dish.

Non-coding RNA in apicomplexan parasites

In recent years it has became evident that the transcriptome of most species has little protein-coding capacity and that the abundance of non-coding RNA was previously overlooked. Non-coding RNAs were initially thought to be transcriptional noise, however, a growing number of studies is showing that many of these RNAs have important regulatory functions. Here, we review the progress done in apicomplexan parasites in this rapidly growing field.

Chromatin-mediated epigenetic regulation in the malaria parasite Plasmodium falciparum

Malaria is a major public health problem in many developing countries, with the malignant tertian parasite Plasmodium falciparum causing the most malaria-associated mortality. Extensive research, especially with the advancement of genomics and transfection tools, has highlighted the fundamental importance of chromatin-mediated gene regulation in the developmental program of this early-branching eukaryote. The Plasmodium parasite genomes reveal the existence of both canonical and variant histones that make up the nucleosomes, as well as a full collection of conserved enzymes for chromatin remodeling and histone posttranslational modifications (PTMs). Recent studies have identified a wide array of both conserved and novel histone PTMs in P. falciparum, indicating the presence of a complex and divergent "histone code." Genome-wide analysis has begun to decipher the nucleosome landscape and histone modifications associated with the dynamic organization of chromatin structures during the parasite's life cycle. Focused studies on malaria-specific phenomena such as antigenic variation and red cell invasion pathways shed further light on the involvement of epigenetic mechanisms in these processes. Here we review our current understanding of chromatin-mediated gene regulation in malaria parasites, with specific reference to exemplar studies on antigenic variation and host cell invasion.

Nucleosome landscape and control of transcription in the human malaria parasite

In eukaryotic cells, chromatin reorganizes within promoters of active genes to allow the transcription machinery and various transcription factors to access DNA. In this model, promoter-specific transcription factors bind DNA to initiate the production of mRNA in a tightly regulated manner. In the case of the human malaria parasite, Plasmodium falciparum, specific transcription factors are apparently underrepresented with regards to the size of the genome, and mechanisms underlying transcriptional regulation are controversial. Here, we investigate the modulation of DNA accessibility by chromatin remodeling during the parasite infection cycle. We have generated genome-wide maps of nucleosome occupancy across the parasite erythrocytic cycle using two complementary assays--the formaldehyde-assisted isolation of regulatory elements to extract protein-free DNA (FAIRE) and the MNase-mediated purification of mononucleosomes to extract histone-bound DNA (MAINE), both techniques being coupled to high-throughput sequencing. We show that chromatin architecture undergoes drastic upheavals throughout the parasite's cycle, contrasting with targeted chromatin reorganization usually observed in eukaryotes. Chromatin loosens after the invasion of the red blood cell and then repacks prior to the next cycle. Changes in nucleosome occupancy within promoter regions follow this genome-wide pattern, with a few exceptions such as the var genes involved in virulence and genes expressed at early stages of the cycle. We postulate that chromatin structure and nucleosome turnover control massive transcription during the erythrocytic cycle. Our results demonstrate that the processes driving gene expression in Plasmodium challenge the classical eukaryotic model of transcriptional regulation occurring mostly at the transcription initiation level.

Simultaneous transcription of duplicated var2csa gene copies in individual Plasmodium falciparum parasites

BACKGROUND

Single nucleotide polymorphisms are common in duplicated genes, causing functional preservation, alteration or silencing. The Plasmodium falciparum genes var2csa and Pf332 are duplicated in the haploid genome of the HB3 parasite line. Whereas the molecular function of Pf332 remains to be elucidated, VAR2CSA is known to be the main adhesin in placental parasite sequestration. Sequence variations introduced upon duplication of these genes provide discriminative possibilities to analyze allele-specific transcription with a bearing towards understanding gene dosage impact on parasite biology.

RESULTS

We demonstrate an approach combining real-time PCR allelic discrimination and discriminative RNA-FISH to distinguish between highly similar gene copies in P. falciparum parasites. The duplicated var2csa variants are simultaneously transcribed, both on a population level and intriguingly also in individual cells, with nuclear co-localization of the active genes and corresponding transcripts. This indicates transcriptional functionality of duplicated genes, challenges the dogma of mutually exclusive var gene transcription and suggests mechanisms behind antigenic variation, at least in respect to the duplicated and highly similar var2csa genes.

CONCLUSIONS

Allelic discrimination assays have traditionally been applied to study zygosity in diploid genomes. The assays presented here are instead successfully applied to the identification and evaluation of transcriptional activity of duplicated genes in the haploid genome of the P. falciparum parasite. Allelic discrimination and gene or transcript localization by FISH not only provide insights into transcriptional regulation of genes such as the virulence associated var genes, but also suggest that this sensitive and precise approach could be used for further investigation of genome dynamics and gene regulation.

Clonally variant gene families in Plasmodium falciparum share a common activation factor

The genome of the malaria parasite Plasmodium falciparum contains several multicopy gene families, including var, rifin, stevor and Pfmc-2TM. These gene families undergo expression switching and appear to play a role in antigenic variation. It has recently been shown that forcing parasites to express high copy numbers of transcriptionally active, episomal var promoters led to gradual downregulation and eventual silencing of the entire var gene family, suggesting that a limiting titratable factor plays a role in var gene activation. Through similar experiments using rifin, stevor or Pfmc-2TM episomal promoters we show that promoter titration can be used as a general method to downregulate multicopy gene families in P. falciparum. Additionally, we show that promoter titration with var, rifin, stevor or Pfmc-2TM episomal promoters results in downregulation of expression not only of the family to which the episomal promoter belongs, but also members of the other gene families, suggesting that the var-specific titratable factor previously described is shared by all four families. Further, transcriptionally active promoters from different families colocalize within the same subnuclear expression site, indicating that the role that nuclear architecture plays in var gene regulation also likely applies to the other multicopy gene families of P. falciparum.

Genetics of antigenic variation in Plasmodium falciparum

Malaria caused by the protozoan parasite Plasmodium falciparum is characterized by long-term, persistent infections that can last for many months. The ability of this parasite to avoid clearance by the human immune system is dependent on its capacity to continuously alter the surface exposed antigenic proteins that that are vulnerable to antibody recognition and attack, a process called antigenic variation. Significant work in recent years has contributed to our understanding of the mechanisms underlying this process, including the genes encoding the antigenic proteins and the DNA sequence elements that control their expression. In addition, the epigenetic "marks" that are associated with activation and silencing of individual genes have been extensively characterized. These studies have led to a model that includes multiple layers of regulation that ultimately lead to the tight coordination of expression of the genes responsible for antigenic variation by malaria parasites. Here we review some more recent data that adds additional complexity to our understanding of these regulatory layers.