Telomeric heterochromatin in Plasmodium falciparum

Single-cell transcriptomics reveal transcriptional programs underlying male and female cell fate during Plasmodium falciparum gametocytogenesis

The Plasmodium falciparum life cycle includes obligate transition between a human and mosquito host. Gametocytes are responsible for transmission from the human to the mosquito vector where gamete fusion followed by meiosis occurs. To elucidate how male and female gametocytes differentiate in the absence of sex chromosomes, we perform FACS-based cell enrichment of a P. falciparum gametocyte reporter line followed by single-cell RNA-seq. In our analyses we define the transcriptional programs and predict candidate driver genes underlying male and female development, including genes from the ApiAP2 family of transcription factors. A motif-driven, gene regulatory network analysis indicates that AP2-G5 specifically modulates male development. Additionally, genes linked to the inner membrane complex, involved in morphological changes, are uniquely expressed in the female lineage. The transcriptional programs of male and female development detailed herein allow for further exploration of the evolution of sex in eukaryotes and provide targets for future development of transmission blocking therapies.

Telomere maintenance in African trypanosomes

Telomere maintenance is essential for genome integrity and chromosome stability in eukaryotic cells harboring linear chromosomes, as telomere forms a specialized structure to mask the natural chromosome ends from DNA damage repair machineries and to prevent nucleolytic degradation of the telomeric DNA. In Trypanosoma brucei and several other microbial pathogens, virulence genes involved in antigenic variation, a key pathogenesis mechanism essential for host immune evasion and long-term infections, are located at subtelomeres, and expression and switching of these major surface antigens are regulated by telomere proteins and the telomere structure. Therefore, understanding telomere maintenance mechanisms and how these pathogens achieve a balance between stability and plasticity at telomere/subtelomere will help develop better means to eradicate human diseases caused by these pathogens. Telomere replication faces several challenges, and the "end replication problem" is a key obstacle that can cause progressive telomere shortening in proliferating cells. To overcome this challenge, most eukaryotes use telomerase to extend the G-rich telomere strand. In addition, a number of telomere proteins use sophisticated mechanisms to coordinate the telomerase-mediated de novo telomere G-strand synthesis and the telomere C-strand fill-in, which has been extensively studied in mammalian cells. However, we recently discovered that trypanosomes lack many telomere proteins identified in its mammalian host that are critical for telomere end processing. Rather, T. brucei uses a unique DNA polymerase, PolIE that belongs to the DNA polymerase A family (E. coli DNA PolI family), to coordinate the telomere G- and C-strand syntheses. In this review, I will first briefly summarize current understanding of telomere end processing in mammals. Subsequently, I will describe PolIE-mediated coordination of telomere G- and C-strand synthesis in T. brucei and implication of this recent discovery.

Architecture, Chromatin and Gene Organization of Toxoplasma gondii Subtelomeres

Subtelomeres (ST) are chromosome regions that separate telomeres from euchromatin and play relevant roles in various biological processes of the cell. While their functions are conserved, ST structure and genetic compositions are unique to each species. This study aims to identify and characterize the subtelomeric regions of the 13 Toxoplasma gondii chromosomes of the Me49 strain. Here, STs were defined at chromosome ends based on poor gene density. The length of STs ranges from 8.1 to 232.4 kbp, with a gene density of 0.049 genes/kbp, lower than the Me49 genome (0.15 kbp). Chromatin organization showed that H3K9me3, H2A.X, and H3.3 are highly enriched near telomeres and the 5' end of silenced genes, decaying in intensity towards euchromatin. H3K4me3 and H2A.Z/H2B.Z are shown to be enriched in the 5' end of the ST genes. Satellite DNA was detected in almost all STs, mainly the sat350 family and a novel satellite named sat240. Beyond the STs, only short dispersed fragments of sat240 and sat350 were found. Within STs, there were 12 functional annotated genes, 59 with unknown functions (Hypothetical proteins), 15 from multigene FamB, and 13 from multigene family FamC. Some genes presented low interstrain synteny associated with the presence of satellite DNA. Orthologues of FamB and FamC were also detected in Neospora caninum and Hammondia hammondi. A re-analysis of previous transcriptomic data indicated that ST gene expression is strongly linked to the adaptation to different situations such as extracellular passage (evolve and resequencing study) and changes in metabolism (lack of acetyl-CoA cofactor). In conclusion, the ST region of the T. gondii chromosomes was defined, the STs genes were determined, and it was possible to associate them with high interstrain plasticity and a role in the adaptability of T. gondii to environmental changes.

Proteomic analysis of Plasmodium falciparum histone deacetylase 1 complex proteins

Plasmodium falciparum histone deacetylases (PfHDACs) are an important class of epigenetic regulators that alter protein lysine acetylation, contributing to regulation of gene expression and normal parasite growth and development. PfHDACs are therefore under investigation as drug targets for malaria. Despite this, our understanding of the biological roles of these enzymes is only just beginning to emerge. In higher eukaryotes, HDACs function as part of multi-protein complexes and act on both histone and non-histone substrates. Here, we present a proteomics analysis of PfHDAC1 immunoprecipitates, identifying 26 putative P. falciparum complex proteins in trophozoite-stage asexual intraerythrocytic parasites. The co-migration of two of these (P. falciparum heat shock proteins 70-1 and 90) with PfHDAC1 was validated using Blue Native PAGE combined with Western blot. These data provide a snapshot of possible PfHDAC1 interactions and a starting point for future studies focused on elucidating the broader function of PfHDACs in Plasmodium parasites.

Modeling of Plasmodium falciparum Telomerase Reverse Transcriptase Ternary Complex: Repurposing of Nucleoside Analog Inhibitors

The Plasmodium falciparum telomerase reverse transcriptase (PfTERT) is a ribonucleoprotein that assists the maintenance of the telomeric ends of chromosomes by reverse transcription of its own RNA subunit. It represents an attractive therapeutic target for eradication of the plasmodial parasite at the asexual liver stage. Automated modeling using MUSTER and knowledge-based techniques were used to obtain a three-dimensional model of the active site of reverse transcriptase domain of PfTERT, which is responsible for catalyzing the addition of incoming dNTPs to the growing DNA strand in presence of divalent magnesium ions. Further, the ternary complex of the active site of PfTERT bound to a DNA-RNA duplex was also modeled using Haddock server and represents the functional form of the enzyme. Initially, established nucleoside analog inhibitors of PfTERT, AZTTP, and ddGTP were docked in the modeled binding site of the PfTERT ternary complex using AutoDock v4.2. Subsequently, docking studies were carried out with 14 approved nucleoside analog inhibitors. Docking studies predicted that floxuridine, gemcitabine, stavudine, and vidarabine have high affinity for the PfTERT ternary complex. Further analysis on the basis of known side effects led us to propose repositioning of vidarabine as a suitable drug candidate for inhibition of PfTERT.

Genomic characterization of variable surface antigens reveals a telomere position effect as a prerequisite for RNA interference-mediated silencing in Paramecium tetraurelia

Antigenic or phenotypic variation is a widespread phenomenon of expression of variable surface protein coats on eukaryotic microbes. To clarify the mechanism behind mutually exclusive gene expression, we characterized the genetic properties of the surface antigen multigene family in the ciliate Paramecium tetraurelia and the epigenetic factors controlling expression and silencing. Genome analysis indicated that the multigene family consists of intrachromosomal and subtelomeric genes; both classes apparently derive from different gene duplication events: whole-genome and intrachromosomal duplication. Expression analysis provides evidence for telomere position effects, because only subtelomeric genes follow mutually exclusive transcription. Microarray analysis of cultures deficient in Rdr3, an RNA-dependent RNA polymerase, in comparison to serotype-pure wild-type cultures, shows cotranscription of a subset of subtelomeric genes, indicating that the telomere position effect is due to a selective occurrence of Rdr3-mediated silencing in subtelomeric regions. We present a model of surface antigen evolution by intrachromosomal gene duplication involving the maintenance of positive selection of structurally relevant regions. Further analysis of chromosome heterogeneity shows that alternative telomere addition regions clearly affect transcription of closely related genes. Consequently, chromosome fragmentation appears to be of crucial importance for surface antigen expression and evolution. Our data suggest that RNAi-mediated control of this genetic network by trans-acting RNAs allows rapid epigenetic adaptation by phenotypic variation in combination with long-term genetic adaptation by Darwinian evolution of antigen genes.

IMPORTANCE

Alternating surface protein structures have been described for almost all eukaryotic microbes, and a broad variety of functions have been described, such as virulence factors, adhesion molecules, and molecular camouflage. Mechanisms controlling gene expression of variable surface proteins therefore represent a powerful tool for rapid phenotypic variation across kingdoms in pathogenic as well as free-living eukaryotic microbes. However, the epigenetic mechanisms controlling synchronous expression and silencing of individual genes are hardly understood. Using the ciliate Paramecium tetraurelia as a (epi)genetic model, we showed that a subtelomeric gene position effect is associated with the selective occurrence of RNAi-mediated silencing of silent surface protein genes, suggesting small interfering RNA (siRNA)-mediated epigenetic cross talks between silent and active surface antigen genes. Our integrated genomic and molecular approach discloses the correlation between gene position effects and siRNA-mediated trans-silencing, thus providing two new parameters for regulation of mutually exclusive gene expression and the genomic organization of variant gene families.

Characterization of Toxoplasma gondii subtelomeric-like regions: identification of a long-range compositional bias that is also associated with gene-poor regions

BACKGROUND

Chromosome ends are composed of telomeric repeats and subtelomeric regions, which are patchworks of genes interspersed with repeated elements. Although chromosome ends display similar arrangements in different species, their sequences are highly divergent. In addition, these regions display a particular nucleosomal composition and bind specific factors, therefore producing a special kind of heterochromatin. Using data from currently available draft genomes we have characterized these putative Telomeric Associated Sequences in Toxoplasma gondii.

RESULTS

An all-vs-all pairwise comparison of T. gondii assembled chromosomes revealed the presence of conserved regions of ∼ 30 Kb located near the ends of 9 of the 14 chromosomes of the genome of the ME49 strain. Sequence similarity among these regions is ∼ 70%, and they are also highly conserved in the GT1 and VEG strains. However, they are unique to Toxoplasma with no detectable similarity in other Apicomplexan parasites. The internal structure of these sequences consists of 3 repetitive regions separated by high-complexity sequences without annotated genes, except for a gene from the Toxoplasma Specific Family. ChIP-qPCR experiments showed that nucleosomes associated to these sequences are enriched in histone H4 monomethylated at K20 (H4K20me1), and the histone variant H2A.X, suggesting that they are silenced sequences (heterochromatin). A detailed characterization of the base composition of these sequences, led us to identify a strong long-range compositional bias, which was similar to that observed in other genomic silenced fragments such as those containing centromeric sequences, and was negatively correlated to gene density.

CONCLUSIONS

We identified and characterized a region present in most Toxoplasma assembled chromosomes. Based on their location, sequence features, and nucleosomal markers we propose that these might be part of subtelomeric regions of T. gondii. The identified regions display a unique trinucleotide compositional bias, which is shared (despite the lack of any detectable sequence similarity) with other silenced sequences, such as those making up the chromosome centromeres. We also identified other genomic regions with this compositional bias (but no detectable sequence similarity) that might be functionally similar.

Analysis of the spatial and temporal arrangement of transcripts over intergenic regions in the human malarial parasite Plasmodium falciparum

Background

The ability of the human malarial parasite Plasmodium falciparum to invade, colonise and multiply within diverse host environments, as well as to manifest its virulence within the human host, are activities tightly linked to the temporal and spatial control of gene expression. Yet, despite the wealth of high throughput transcriptomic data available for this organism there is very little information regarding the location of key transcriptional landmarks or their associated cis-acting regulatory elements. Here we provide a systematic exploration of the size and organisation of transcripts within intergenic regions to yield surrogate information regarding transcriptional landmarks, and to also explore the spatial and temporal organisation of transcripts over these poorly characterised genomic regions.

Results

Utilising the transcript data for a cohort of 105 genes we demonstrate that the untranscribed regions of mRNA are large and apportioned predominantly to the 5′ end of the open reading frame. Given the relatively compact size of the P. falciparum genome, we suggest that whilst transcriptional units are likely to spatially overlap, temporal co-transcription of adjacent transcriptional units is actually limited. Critically, the size of intergenic regions is directly dependent on the orientation of the two transcriptional units arrayed over them, an observation we extend to an analysis of the complete sequences of twelve additional organisms that share moderately compact genomes.

Conclusions

Our study provides a theoretical framework that extends our current understanding of the transcriptional landscape across the P. falciparum genome. Demonstration of a consensus gene-spacing rule that is shared between P. falciparum and ten other moderately compact genomes of apicomplexan parasites reveals the potential for our findings to have a wider impact across a phylum that contains many organisms important to human and veterinary health.

Telomeres, tethers and trypanosomes

Temporal and spatial organization of the nucleus is critical for the control of transcription, mRNA processing and the assembly of ribosomes. This includes the occupancy of specific territories by mammalian chromosomes, the presence of subnuclear compartments such as the nucleolus and Cajal bodies and the division of chromatin between active and inactive states. These latter are commonly associated with the location of DNA within euchromatin and heterochromatin respectively; critically these distinctions arise through modifications to chromatin-associated proteins, including histones, as well as the preferential localization of heterochromatin at the nuclear periphery. Most research on nuclear organization has focused on metazoa and fungi; however, recent technical advances have made more divergent eukaryotes accessible to study, with some surprising results. For example, the organization of heterochromatin is mediated in metazoan nuclei in large part by lamins, the prototypical intermediate filament proteins. Despite the presence of heterochromatin, detected both biochemically and by EM in most eukaryotic organisms, until this year lamins were thought to be restricted to metazoan taxa, and the proteins comprising the lamina in other lineages were unknown. Recent work indicates the presence of lamin orthologs in amoeba, while trypanosomatids possess a large coiled-coil protein, NUP-1, that performs functions analogous to lamins. These data indicate that the presence of a nuclear lamina is substantially more widespread than previously thought, with major implications for the evolution of eukaryotic gene expression mechanisms. We discuss these and other recent findings on the organization of nuclei in diverse organisms, and the implications of these findings for the evolutionary origin of eukaryotes.

Sirtuins of parasitic protozoa: in search of function(s)

The SIR2 family of NAD⁺-dependent protein deacetylases, collectively called sirtuins, has been of central interest due to their proposed roles in life-span regulation and ageing. Sirtuins are one group of environment sensors of a cell interpreting external information and orchestrating internal responses at the sub-cellular level, through participation in gene regulation mechanisms. Remarkably conserved across all kingdoms of life SIR2 proteins in several protozoan parasites appear to have both conserved and intriguing unique functions. This review summarises our current knowledge of the members of the sirtuin families in Apicomplexa, including Plasmodium, and other protozoan parasites such as Trypanosoma and Leishmania. The wide diversity of processes regulated by SIR2 proteins makes them targets worthy of exploitation in anti-parasitic therapies.

Transfection of rodent malaria parasites

Gene manipulation is an invaluable tool to investigate and understand the biology of an organism. Although this technology has been applied to both the human and rodent malarial parasites (RMP), Plasmodium berghei in particular offers a more robust system due to a higher and more efficient transformation rate. Here, we describe a comprehensive transfection and selection protocol using P. berghei including a variant negative selection protocol administering 5-fluorocytosine to the animals in drinking water. Additionally, we discuss and assess the latest advances in gene manipulation technologies developed in RMP to gain a better understanding of Plasmodium biology.

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.