Epigenetic Regulation and Chromatin Remodeling in Malaria Parasites.

Plasmodium falciparum, the human malaria parasite, infects two hosts and various cell types, inducing distinct morphological and physiological changes in the parasite in response to different environmental conditions. These variations required the parasite to adapt and develop elaborate molecular mechanisms to ensure its spread and transmission. Recent findings have significantly improved our understanding of the regulation of gene expression in P. falciparum. Here, we provide an up-to-date overview of technologies used to highlight the transcriptomic adjustments occurring in the parasite throughout its life cycle. We also emphasize the complementary and complex epigenetic mechanisms regulating gene expression in malaria parasites. This review concludes with an outlook on the chromatin architecture, the remodeling systems, and how this 3D genome organization is critical in various biological processes.

The Toxoplasma gondii F-Box Protein L2 Functions as a Repressor of Stage Specific Gene Expression.

Toxoplasma gondii is a foodborne pathogen that can cause severe and life-threatening infections in fetuses and immunocompromised patients. Felids are its only definitive hosts, and a wide range of animals, including humans, serve as intermediate hosts. When the transmissible bradyzoite stage is orally ingested by felids, they transform into merozoites that expand asexually, ultimately generating millions of gametes for the parasite sexual cycle. However, bradyzoites in intermediate hosts differentiate exclusively to disease-causing tachyzoites, which rapidly disseminate throughout the host. Though tachyzoites are well-studied, the molecular mechanisms governing transitioning between developmental stages are poorly understood. Each parasite stage can be distinguished by a characteristic transcriptional signature, with one signature being repressed during the other stages. Switching between stages require substantial changes in the proteome, which is achieved in part by ubiquitination. F-box proteins mediate protein poly-ubiquitination by recruiting substrates to SKP1, Cullin-1, F-Box protein E3 ubiquitin ligase (SCF-E3) complexes. We have identified an F-box protein named Toxoplasma gondii F-Box Protein L2 (TgFBXL2), which localizes to distinct perinucleolar sites. TgFBXL2 is stably engaged in an SCF-E3 complex that is surprisingly also associated with a COP9 signalosome complex that negatively regulates SCF-E3 function. At the cellular level, TgFBXL2-depleted parasites are severely defective in centrosome replication and daughter cell development. Most remarkable, RNAseq data show that TgFBXL2 conditional depletion induces the expression of stage-specific genes including a large cohort of genes necessary for sexual commitment. Together, these data suggest that TgFBXL2 is a latent guardian of stage specific gene expression in Toxoplasma and poised to remove conflicting proteins in response to an unknown trigger of development.

Genome-wide functional analysis reveals key roles for kinesins in the mammalian and mosquito stages of the malaria parasite life cycle.

Kinesins are microtubule (MT)-based motors important in cell division, motility, polarity, and intracellular transport in many eukaryotes. However, they are poorly studied in the divergent eukaryotic pathogens Plasmodium spp., the causative agents of malaria, which manifest atypical aspects of cell division and plasticity of morphology throughout the life cycle in both mammalian and mosquito hosts. Here, we describe a genome-wide screen of Plasmodium kinesins, revealing diverse subcellular locations and functions in spindle assembly, axoneme formation, and cell morphology. Surprisingly, only kinesin-13 is essential for growth in the mammalian host while the other 8 kinesins are required during the proliferative and invasive stages of parasite transmission through the mosquito vector. In-depth analyses of kinesin-13 and kinesin-20 revealed functions in MT dynamics during apical cell polarity formation, spindle assembly, and axoneme biogenesis. These findings help us to understand the importance of MT motors and may be exploited to discover new therapeutic interventions against malaria.

Dynamic Chromatin Structure and Epigenetics Control the Fate of Malaria Parasites.

Multiple hosts and various life cycle stages prompt the human malaria parasite, Plasmodium falciparum, to acquire sophisticated molecular mechanisms to ensure its survival, spread, and transmission to its next host. To face these environmental challenges, increasing evidence suggests that the parasite has developed complex and complementary layers of regulatory mechanisms controlling gene expression. Here, we discuss the recent developments in the discovery of molecular components that contribute to cell replication and differentiation and highlight the major contributions of epigenetics, transcription factors, and nuclear architecture in controlling gene regulation and life cycle progression in Plasmodium spp.

Third-generation sequencing revises the molecular karyotype for Toxoplasma gondii and identifies emerging copy number variants in sexual recombinants.

Toxoplasma gondii is a useful model for intracellular parasitism given its ease of culture in the laboratory and genomic resources. However, as for many other eukaryotes, the T. gondii genome contains hundreds of sequence gaps owing to repetitive and/or unclonable sequences that disrupt the assembly process. Here, we use the Oxford Nanopore Minion platform to generate near-complete de novo genome assemblies for multiple strains of T. gondii and its near relative, N. caninum We significantly improved T. gondii genome contiguity (average N50 of ∼6.6 Mb) and added ∼2 Mb of newly assembled sequence. For all of the T. gondii strains that we sequenced (RH, ME49, CTG, II×III progeny clones CL13, S27, S21, S26, and D3X1), the largest contig ranged in size between 11.9 and 12.1 Mb in size, which is larger than any previously reported T. gondii chromosome, and found to be due to a consistent fusion of Chromosomes VIIb and VIII. These data were validated by mapping existing T. gondii ME49 Hi-C data to our assembly, providing parallel lines of evidence that the T. gondii karyotype consists of 13, rather than 14, chromosomes. By using this technology, we also resolved hundreds of tandem repeats of varying lengths, including in well-known host-targeting effector loci like rhoptry protein 5 (ROP5) and ROP38 Finally, when we compared T. gondii with N. caninum, we found that although the 13-chromosome karyotype was conserved, extensive, previously unappreciated chromosome-scale rearrangements had occurred in T. gondii and N. caninum since their most recent common ancestry.

An Enantiospecific Synthesis of Isoneoamphilectane Confirms Its Strained Tricyclic Structure.

We describe a total synthesis of the rare isocyanoterpene natural product isoneoamphilectane and two of its unnatural diastereomers. The significantly strained ring system of the reported natural product─along with a hypothesis about a biosynthetic relationship to related family members─inspired us to consider a potential misassignment in the structure's relative configuration. As a result, we initially targeted two less strained, more accessible, stereoisomers of the reported natural product. When these compounds failed to exhibit spectroscopic data that matched those of isoneoamphilectane, we embarked on a synthesis of the originally proposed strained structure via an approach that hinged on a challenging cis-to-trans decalone epimerization. Ultimately, we implemented a novel cyclic sulfite pinacol-type rearrangement to generate the strained ring system. Additional features of this work include the application of a stereocontrolled Mukaiyama-Michael addition of an acyclic silylketene acetal, an unusual intramolecular alkoxide-mediated regioselective elimination, and an HAT-mediated alkene hydroazidation to forge the C-N bond of the tertiary isonitrile. Throughout this work, our synthetic planning was heavily guided by computational analyses to inform on key issues of stereochemical control.

Three-dimensional chromatin in infectious disease-A role for gene regulation and pathogenicity?

The recent Coronavirus Disease 2019 pandemic has once again reminded us the importance of understanding infectious diseases. One important but understudied area in infectious disease research is the role of nuclear architecture or the physical arrangement of the genome in the nucleus in controlling gene regulation and pathogenicity. Recent advances in research methods, such as Genome-wide chromosome conformation capture using high-throughput sequencing (Hi-C), have allowed for easier analysis of nuclear architecture and chromosomal reorganization in both the infectious disease agents themselves as well as in their host cells. This review will discuss broadly on what is known about nuclear architecture in infectious disease, with an emphasis on chromosomal reorganization, and briefly discuss what steps are required next in the field.

The Arabidopsis PHD-finger protein EDM2 has multiple roles in balancing NLR immune receptor gene expression.

Plant NLR-type receptors serve as sensitive triggers of host immunity. Their expression has to be well-balanced, due to their interference with various cellular processes and dose-dependency of their defense-inducing activity. A genetic "arms race" with fast-evolving pathogenic microbes requires plants to constantly innovate their NLR repertoires. We previously showed that insertion of the COPIA-R7 retrotransposon into RPP7 co-opted the epigenetic transposon silencing signal H3K9me2 to a new function promoting expression of this Arabidopsis thaliana NLR gene. Recruitment of the histone binding protein EDM2 to COPIA-R7-associated H3K9me2 is required for optimal expression of RPP7. By profiling of genome-wide effects of EDM2, we now uncovered additional examples illustrating effects of transposons on NLR gene expression, strongly suggesting that these mobile elements can play critical roles in the rapid evolution of plant NLR genes by providing the "raw material" for gene expression mechanisms. We further found EDM2 to have a global role in NLR expression control. Besides serving as a positive regulator of RPP7 and a small number of other NLR genes, EDM2 acts as a suppressor of a multitude of additional NLR genes. We speculate that the dual functionality of EDM2 in NLR expression control arose from the need to compensate for fitness penalties caused by high expression of some NLR genes by suppression of others. Moreover, we are providing new insights into functional relationships of EDM2 with its interaction partner, the RNA binding protein EDM3/AIPP1, and its target gene IBM1, encoding an H3K9-demethylase.

Real-time dynamics of Plasmodium NDC80 reveals unusual modes of chromosome segregation during parasite proliferation.

Eukaryotic cell proliferation requires chromosome replication and precise segregation to ensure daughter cells have identical genomic copies. Species of the genus Plasmodium, the causative agents of malaria, display remarkable aspects of nuclear division throughout their life cycle to meet some peculiar and unique challenges to DNA replication and chromosome segregation. The parasite undergoes atypical endomitosis and endoreduplication with an intact nuclear membrane and intranuclear mitotic spindle. To understand these diverse modes of Plasmodium cell division, we have studied the behaviour and composition of the outer kinetochore NDC80 complex, a key part of the mitotic apparatus that attaches the centromere of chromosomes to microtubules of the mitotic spindle. Using NDC80-GFP live-cell imaging in Plasmodium berghei, we observe dynamic spatiotemporal changes during proliferation, including highly unusual kinetochore arrangements during sexual stages. We identify a very divergent candidate for the SPC24 subunit of the NDC80 complex, previously thought to be missing in Plasmodium, which completes a canonical, albeit unusual, NDC80 complex structure. Altogether, our studies reveal the kinetochore to be an ideal tool to investigate the non-canonical modes of chromosome segregation and cell division in Plasmodium.

Measuring significant changes in chromatin conformation with ACCOST.

Chromatin conformation assays such as Hi-C cannot directly measure differences in 3D architecture between cell types or cell states. For this purpose, two or more Hi-C experiments must be carried out, but direct comparison of the resulting Hi-C matrices is confounded by several features of Hi-C data. Most notably, the genomic distance effect, whereby contacts between pairs of genomic loci that are proximal along the chromosome exhibit many more Hi-C contacts that distal pairs of loci, dominates every Hi-C matrix. Furthermore, the form that this distance effect takes often varies between different Hi-C experiments, even between replicate experiments. Thus, a statistical confidence measure designed to identify differential Hi-C contacts must accurately account for the genomic distance effect or risk being misled by large-scale but artifactual differences. ACCOST (Altered Chromatin COnformation STatistics) accomplishes this goal by extending the statistical model employed by DEseq, re-purposing the 'size factors,' which were originally developed to account for differences in read depth between samples, to instead model the genomic distance effect. We show via analysis of simulated and real data that ACCOST provides unbiased statistical confidence estimates that compare favorably with competing methods such as diffHiC, FIND and HiCcompare. ACCOST is freely available with an Apache license at https://bitbucket.org/noblelab/accost.

Comparative 3D genome organization in apicomplexan parasites.

The positioning of chromosomes in the nucleus of a eukaryotic cell is highly organized and has a complex and dynamic relationship with gene expression. In the human malaria parasite Plasmodium falciparum, the clustering of a family of virulence genes correlates with their coordinated silencing and has a strong influence on the overall organization of the genome. To identify conserved and species-specific principles of genome organization, we performed Hi-C experiments and generated 3D genome models for five Plasmodium species and two related apicomplexan parasites. Plasmodium species mainly showed clustering of centromeres, telomeres, and virulence genes. In P. falciparum, the heterochromatic virulence gene cluster had a strong repressive effect on the surrounding nuclear space, while this was less pronounced in Plasmodium vivax and Plasmodium berghei, and absent in Plasmodium yoelii In Plasmodium knowlesi, telomeres and virulence genes were more dispersed throughout the nucleus, but its 3D genome showed a strong correlation with gene expression. The Babesia microti genome showed a classical Rabl organization with colocalization of subtelomeric virulence genes, while the Toxoplasma gondii genome was dominated by clustering of the centromeres and lacked virulence gene clustering. Collectively, our results demonstrate that spatial genome organization in most Plasmodium species is constrained by the colocalization of virulence genes. P. falciparum and P. knowlesi, the only two Plasmodium species with gene families involved in antigenic variation, are unique in the effect of these genes on chromosome folding, indicating a potential link between genome organization and gene expression in more virulent pathogens.

Unraveling the 3D genome of human malaria parasites.

The chromosomes within the eukaryotic cell nucleus are highly dynamic and adopt complex hierarchical structures. Understanding how this three-dimensional (3D) nuclear architectureaffects gene regulation, cell cycle progression and disease pathogenesis are important biological questions in development and disease. Recently, many genome-wide technologies including chromosome conformation capture (3C) and 3C-based methodologies (4C, 5C, and Hi-C) have been developed to investigate 3D chromatin structure. In this review, we introduce 3D genome methodologies, with a focus on their application for understanding the nuclear architecture of the human malaria parasite, Plasmodium falciparum. An increasing amount of evidence now suggests that gene regulation in the parasite is largely regulated by epigenetic mechanisms and nuclear reorganization. Here, we explore the 3D genome architecture of P. falciparum, including local and global chromatin structure. In addition, molecular components important for maintaining 3D chromatin organization including architectural proteins and long non-coding RNAs are discussed. Collectively, these studies contribute to our understanding of how the plasticity of 3D genome architecture regulates gene expression and cell cycle progression in this deadly parasite.

Plasmodium kinesin-8X associates with mitotic spindles and is essential for oocyst development during parasite proliferation and transmission.

Kinesin-8 proteins are microtubule motors that are often involved in regulation of mitotic spindle length and chromosome alignment. They move towards the plus ends of spindle microtubules and regulate the dynamics of these ends due, at least in some species, to their microtubule depolymerization activity. Plasmodium spp. exhibit an atypical endomitotic cell division in which chromosome condensation and spindle dynamics in the different proliferative stages are not well understood. Genome-wide shared orthology analysis of Plasmodium spp. revealed the presence of two kinesin-8 motor proteins, kinesin-8X and kinesin-8B. Here we studied the biochemical properties of kinesin-8X and its role in parasite proliferation. In vitro, kinesin-8X has motility and depolymerization activities like other kinesin-8 motors. To understand the role of Plasmodium kinesin-8X in cell division, we used fluorescence-tagging and live cell imaging to define its location, and gene targeting to analyse its function, during all proliferative stages of the rodent malaria parasite P. berghei life cycle. The results revealed a spatio-temporal involvement of kinesin-8X in spindle dynamics and an association with both mitotic and meiotic spindles and the putative microtubule organising centre (MTOC). Deletion of the kinesin-8X gene revealed a defect in oocyst development, confirmed by ultrastructural studies, suggesting that this protein is required for oocyst development and sporogony. Transcriptome analysis of Δkinesin-8X gametocytes revealed modulated expression of genes involved mainly in microtubule-based processes, chromosome organisation and the regulation of gene expression, supporting a role for kinesin-8X in cell division. Kinesin-8X is thus required for parasite proliferation within the mosquito and for transmission to the vertebrate host.

The Arabidopsis RRM domain protein EDM3 mediates race-specific disease resistance by controlling H3K9me2-dependent alternative polyadenylation of RPP7 immune receptor transcripts.

The NLR-receptor RPP7 mediates race-specific immunity in Arabidopsis. Previous screens for enhanced downy mildew (edm) mutants identified the co-chaperone SGT1b (EDM1) and the PHD-finger protein EDM2 as critical regulators of RPP7. Here, we describe a third edm mutant compromised in RPP7 immunity, edm3. EDM3 encodes a nuclear-localized protein featuring an RNA-recognition motif. Like EDM2, EDM3 promotes histone H3 lysine 9 dimethylation (H3K9me2) at RPP7. Global profiling of H3K9me2 showed EDM3 to affect this silencing mark at a large set of loci. Importantly, both EDM3 and EDM2 co-associate in vivo with H3K9me2-marked chromatin and transcripts at a critical proximal polyadenylation site of RPP7, where they suppress proximal transcript polyadeylation/termination. Our results highlight the complexity of plant NLR gene regulation, and establish a functional and physical link between a histone mark and NLR-transcript processing.

Predicting gene expression in the human malaria parasite Plasmodium falciparum using histone modification, nucleosome positioning, and 3D localization features.

Empirical evidence suggests that the malaria parasite Plasmodium falciparum employs a broad range of mechanisms to regulate gene transcription throughout the organism's complex life cycle. To better understand this regulatory machinery, we assembled a rich collection of genomic and epigenomic data sets, including information about transcription factor (TF) binding motifs, patterns of covalent histone modifications, nucleosome occupancy, GC content, and global 3D genome architecture. We used these data to train machine learning models to discriminate between high-expression and low-expression genes, focusing on three distinct stages of the red blood cell phase of the Plasmodium life cycle. Our results highlight the importance of histone modifications and 3D chromatin architecture in Plasmodium transcriptional regulation and suggest that AP2 transcription factors may play a limited regulatory role, perhaps operating in conjunction with epigenetic factors.

Design and tests of prospective property predictions for novel antimalarial 2-aminopropylaminoquinolones.

There is a pressing need to improve the efficiency of drug development, and nowhere is that need more clear than in the case of neglected diseases like malaria. The peculiarities of pyrimidine metabolism in Plasmodium species make inhibition of dihydroorotate dehydrogenase (DHODH) an attractive target for antimalarial drug design. By applying a pair of complementary quantitative structure-activity relationships derived for inhibition of a truncated, soluble form of the enzyme from Plasmodium falciparum (s-PfDHODH) to data from a large-scale phenotypic screen against cultured parasites, we were able to identify a class of antimalarial leads that inhibit the enzyme and abolish parasite growth in blood culture. Novel analogs extending that class were designed and synthesized with a goal of improving potency as well as the general pharmacokinetic and toxicological profiles. Their synthesis also represented an opportunity to prospectively validate our in silico property predictions. The seven analogs synthesized exhibited physicochemical properties in good agreement with prediction, and five of them were more active against P. falciparum growing in blood culture than any of the compounds in the published lead series. The particular analogs prepared did not inhibit s-PfDHODH in vitro, but advanced biological assays indicated that other examples from the class did inhibit intact PfDHODH bound to the mitochondrial membrane. The new analogs, however, killed the parasites by acting through some other, unidentified mechanism 24-48 h before PfDHODH inhibition would be expected to do so.

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

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

Concise Synthesis of the Antiplasmodial Isocyanoterpene 7,20-Diisocyanoadociane.

The flagship member of the antiplasmodial isocyanoterpenes, 7,20-diisocyanoadociane (DICA), was synthesized from dehydrocryptone in 10 steps, and in 13 steps from commercially available material. Our previous formal synthesis was reengineered, leveraging only productive transformations to deliver DICA in fewer than half the number of steps of our original effort. Important contributions, in addition to the particularly concise strategy, include a solution to the problem of axial nucleophilic methylation of a late-stage cyclohexanone, and the first selective synthesis and antiplasmodial evaluation of the DICA stereoisomer with both isonitriles equatorial.

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

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

Multiple dimensions of epigenetic gene regulation in the malaria parasite Plasmodium falciparum: gene regulation via histone modifications, nucleosome positioning and nuclear architecture in P. falciparum.

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