Genome-wide localization of histone variants in Toxoplasma gondii implicates variant exchange in stage-specific gene expression

Histone code: a common language and multiple dialects to meet the different developmental requirements of apicomplexan parasites

Apicomplexan parasites have complex life cycles, often requiring transmission between two different hosts, facing periods of dormancy within the host or in the environment to maximize chances of transmission. To support survival in these different conditions, tightly regulated and correctly timed gene expression is critical. The modification of histones and nucleosome composition makes a significant contribution to this regulation, and as eukaryotes, the fundamental mechanisms underlying this process in apicomplexans are similar to those in model eukaryotic organisms. However, single-celled intracellular parasites face unique challenges, and regulation of gene expression at the epigenetic level provides tight control for responses that must often be rapid and robust. Here, we discuss the recent advances in understanding the dynamics of histone modifications across Apicomplexan life cycles and the molecular mechanisms that underlie epigenetic regulation of gene expression to promote parasite life cycle progression, dormancy, and transmission.

Mechanisms of adaptation and evolution in Toxoplasma gondii

Toxoplasma has high host flexibility, infecting all nucleated cells of mammals and birds. This implies that during its infective process the parasite must constantly adapt to different environmental situations, which in turn leads to modifications in its metabolism, regulation of gene transcription, translation of mRNAs and stage specific factors. There are conserved pathways that support these adaptations, which we aim to elucidate in this review. We begin by exploring the widespread epigenetic mechanisms and transcription regulators, continue with the supportive role of Heat Shock Proteins (Hsp), the translation regulation, stress granules, and finish with the emergence of contingency genes in highly variable genomic domains, such as subtelomeres. Within epigenetics, the discovery of a new histone variant of the H2B family (H2B.Z), contributing to T. gondii virulence and differentiation, but also gene expression regulation and its association with the metabolic state of the parasite, is highlighted. Associated with the regulation of gene expression are transcription factors (TFs). An overview of the main findings on TF and development is presented. We also emphasize the role of Hsp90 and Tgj1 in T. gondii metabolic fitness and the regulation of protein translation. Translation regulation is also highlighted as a mechanism for adaptation to conditions encountered by the parasite as well as stress granules containing mRNA and proteins generated in the extracellular tachyzoite. Another important aspect in evolution and adaptability are the subtelomeres because of their high variability and gene duplication rate. Toxoplasma possess multigene families of membrane proteins and contingency genes that are associated with different metabolic stresses. Among them parasite differentiation and environmental stresses stand out, including those that lead tachyzoite to bradyzoite conversion. Finally, we are interested in positioning protozoa as valuable evolution models, focusing on research related to the Extended Evolutionary Synthesis, based on models recently generated, such as extracellular adaptation and ex vivo cyst recrudescence.

MYST regulates DNA repair and forms a NuA4-like complex in the malaria parasite Plasmodium falciparum

Histone lysine acetyltransferase MYST-associated NuA4 complex is conserved from yeast to humans and plays key roles in cell cycle regulation, gene transcription, and DNA replication/repair. Here, we identified a Plasmodium falciparum MYST-associated complex, PfNuA4, which contains 11 of the 13 conserved NuA4 subunits. Reciprocal pulldowns using PfEAF2, a shared component between the NuA4 and SWR1 complexes, not only confirmed the PfNuA4 complex but also identified the PfSWR1 complex, a histone remodeling complex, although their identities are low compared to the homologs in yeast or humans. Notably, both H2A.Z/H2B.Z were associated with the PfSWR1 complex, indicating that this complex is involved in the deposition of H2A.Z/H2B.Z, the variant histone pair that is enriched in the activated promoters. Overexpression of PfMYST resulted in earlier expression of genes involved in cell cycle regulation, DNA replication, and merozoite invasion, and upregulation of the genes related to antigenic variation and DNA repair. Consistently, PfMYST overexpression led to high basal phosphorylated PfH2A (γ-PfH2A), the mark of DNA double-strand breaks, and conferred protection against genotoxic agent methyl methanesulfonate (MMS), X-rays, and artemisinin, the first-line antimalarial drug. In contrast, the knockdown of PfMYST caused a delayed parasite recovery upon MMS treatment. MMS induced the gradual disappearance of PfMYST in the cytoplasm and concomitant accumulation of PfMYST in the nucleus, suggesting cytoplasm-nucleus shuttling of PfMYST. Meanwhile, PfMYST colocalized with the γ-PfH2A, indicating PfMYST was recruited to the DNA damage sites. Collectively, PfMYST plays critical roles in cell cycle regulation, gene transcription, and DNA replication/DNA repair in this low-branching parasitic protist.IMPORTANCEUnderstanding gene regulation and DNA repair in malaria parasites is critical for identifying targets for antimalarials. This study found PfNuA4, a PfMYST-associated, histone modifier complex, and PfSWR1, a chromatin remodeling complex in malaria parasite Plasmodium falciparum. These complexes are divergent due to the low identities compared to their homologs from yeast and humans. Furthermore, overexpression of PfMYST resulted in substantial transcriptomic changes, indicating that PfMYST is involved in regulating the cell cycle, antigenic variation, and DNA replication/repair. Consistently, PfMYST was found to protect against DNA damage caused by the genotoxic agent methyl methanesulfonate, X-rays, and artemisinin, the first-line antimalarial drug. Additionally, DNA damage led to the relocation of cytoplasmic PfMYST to the nucleus and colocalization of PfMYST with γ-PfH2A, the mark of DNA damage. In summary, this study demonstrated that the PfMYST complex has critical functions in regulating cell cycle, antigenic variation, and DNA replication/DNA repair in P. falciparum.

Proteomics Applications in Toxoplasma gondii: Unveiling the Host-Parasite Interactions and Therapeutic Target Discovery

Toxoplasma gondii, a protozoan parasite with the ability to infect various warm-blooded vertebrates, including humans, is the causative agent of toxoplasmosis. This infection poses significant risks, leading to severe complications in immunocompromised individuals and potentially affecting the fetus through congenital transmission. A comprehensive understanding of the intricate molecular interactions between T. gondii and its host is pivotal for the development of effective therapeutic strategies. This review emphasizes the crucial role of proteomics in T. gondii research, with a specific focus on host-parasite interactions, post-translational modifications (PTMs), PTM crosstalk, and ongoing efforts in drug discovery. Additionally, we provide an overview of recent advancements in proteomics techniques, encompassing interactome sample preparation methods such as BioID (BirA*-mediated proximity-dependent biotin identification), APEX (ascorbate peroxidase-mediated proximity labeling), and Y2H (yeast two hybrid), as well as various proteomics approaches, including single-cell analysis, DIA (data-independent acquisition), targeted, top-down, and plasma proteomics. Furthermore, we discuss bioinformatics and the integration of proteomics with other omics technologies, highlighting its potential in unraveling the intricate mechanisms of T. gondii pathogenesis and identifying novel therapeutic targets.

The P. falciparum alternative histones Pf H2A.Z and Pf H2B.Z are dynamically acetylated and antagonized by PfSir2 histone deacetylases at heterochromatin boundaries

The Plasmodium falciparum alternative histones Pf H2A.Z and Pf H2B.Z are enriched in the same nucleosomes in intergenic euchromatin but depleted from heterochromatin. They occupy most promoters but are only dynamically associated with expression at var genes. In other organisms, acetylation of H2A.Z is important for its functions in gene expression and chromatin structure. Here, we show that acetylated Pf H2A.Z and Pf H2B.Z are dynamically associated with gene expression at promoters. In addition, acetylated Pf H2A.Z and Pf H2B.Z are antagonized by the sirtuin class III histone deacetylases (HDAC) PfSir2A and B at heterochromatin boundaries and encroach upon heterochromatin in parasites lacking PfSir2A or B. However, the majority of acetylated Pf H2A.Z and Pf H2B.Z are deacetylated by class I or II HDACs. Acetylated Pf H2A.Z and Pf H2B.Z are also dynamically associated with promoter activity of both canonical upstream var gene promoters and var gene introns. These findings suggest that both acetylated Pf H2A.Z and Pf H2B.Z play critical roles in gene expression and contribute to maintenance of chromatin structure at the boundaries of subtelomeric, facultative heterochromatin, critical for the variegated expression of genes that enable rapid adaptation to altered host environments.IMPORTANCEThe malaria parasite Plasmodium falciparum relies on variant expression of members of multi-gene families as a strategy for environmental adaptation to promote parasite survival and pathogenesis. These genes are located in transcriptionally silenced DNA regions. A limited number of these genes escape gene silencing, and switching between them confers variant fitness on parasite progeny. Here, we show that PfSir2 histone deacetylases antagonize DNA-interacting acetylated alternative histones at the boundaries between active and silent DNA. This finding implicates acetylated alternative histones in the mechanism regulating P. falciparum variant gene silencing and thus malaria pathogenesis. This work also revealed that acetylation of alternative histones at promoters is dynamically associated with promoter activity across the genome, implicating acetylation of alternative histones in gene regulation genome wide. Understanding mechanisms of gene regulation in P. falciparum may aid in the development of new therapeutic strategies for malaria, which killed 619,000 people in 2021.

Histone variant H2B.Z acetylation is necessary for maintenance of Toxoplasma gondii biological fitness

Through regulation of DNA packaging, histone proteins are fundamental to a wide array of biological processes. A variety of post-translational modifications (PTMs), including acetylation, constitute a proposed histone code that is interpreted by "reader" proteins to modulate chromatin structure. Canonical histones can be replaced with variant versions that add an additional layer of regulatory complexity. The protozoan parasite Toxoplasma gondii is unique among eukaryotes in possessing a novel variant of H2B designated H2B.Z. The combination of PTMs and the use of histone variants are important for gene regulation in T. gondii, offering new targets for drug development. In this work, T. gondii parasites were generated in which the 5 N-terminal acetylatable lysines in H2B.Z were mutated to either alanine (c-Myc-A) or arginine (c-Myc-R). The c-Myc-A mutant displayed no phenotype over than a mild defect in its ability to kill mice. The c-Myc-R mutant presented an impaired ability to grow and an increase in differentiation to latent bradyzoites. The c-Myc-R mutant was also more sensitive to DNA damage, displayed no virulence in mice, and provided protective immunity against future infection. While nucleosome composition was unaltered, key genes were abnormally expressed during in vitro bradyzoite differentiation. Our results show that regulation of the N-terminal positive charge patch of H2B.Z is important for these processes. We also show that acetylated N-terminal H2B.Z interacts with some unique proteins compared to its unacetylated counterpart; the acetylated peptide pulled down proteins associated with chromosome maintenance/segregation and cell cycle, suggesting a link between H2B.Z acetylation status and mitosis.

Babesia duncani multi-omics identifies virulence factors and drug targets

Babesiosis is a malaria-like disease in humans and animals that is caused by Babesia species, which are tick-transmitted apicomplexan pathogens. Babesia duncani causes severe to lethal infection in humans, but despite the risk that this parasite poses as an emerging pathogen, little is known about its biology, metabolic requirements or pathogenesis. Unlike other apicomplexan parasites that infect red blood cells, B. duncani can be continuously cultured in vitro in human erythrocytes and can infect mice resulting in fulminant babesiosis and death. We report comprehensive, detailed molecular, genomic, transcriptomic and epigenetic analyses to gain insights into the biology of B. duncani. We completed the assembly, 3D structure and annotation of its nuclear genome, and analysed its transcriptomic and epigenetics profiles during its asexual life cycle stages in human erythrocytes. We used RNA-seq data to produce an atlas of parasite metabolism during its intraerythrocytic life cycle. Characterization of the B. duncani genome, epigenome and transcriptome identified classes of candidate virulence factors, antigens for diagnosis of active infection and several attractive drug targets. Furthermore, metabolic reconstitutions from genome annotation and in vitro efficacy studies identified antifolates, pyrimethamine and WR-99210 as potent inhibitors of B. duncani to establish a pipeline of small molecules that could be developed as effective therapies for the treatment of human babesiosis.

Histone variant H2B.Z acetylation is necessary for maintenance of Toxoplasma gondii biological fitness

Through regulation of DNA packaging, histone proteins are fundamental to a wide array of biological processes. A variety of post-translational modifications (PTMs), including acetylation, constitute a proposed histone code that is interpreted by "reader" proteins to modulate chromatin structure. Canonical histones can be replaced with variant versions that add an additional layer of regulatory complexity. The protozoan parasite Toxoplasma gondii is unique among eukaryotes in possessing a novel variant of H2B designated H2B.Z. The combination of PTMs and the use of histone variants is important for gene regulation in T. gondii, offering new targets for drug development. In this work, T. gondii parasites were generated in which the 5 N-terminal acetylatable lysines in H2B.Z were mutated to either alanine (c-Myc-A) or arginine (c-Myc-R). c-Myc-A mutant only displayed a mild effect in its ability to kill mice. c-Myc-R mutant presented an impaired ability to grow and an increase in differentiation to latent bradyzoites. This mutant line was also more sensitive to DNA damage, displayed no virulence in mice, and provided protective immunity against future infection. While nucleosome composition was unaltered, key genes were abnormally expressed during in vitro bradyzoite differentiation. Our results show that the N-terminal positive charge patch of H2B.Z is important for these procceses. Pull down assays with acetylated N-terminal H2B.Z peptide and unacetylated one retrieved common and differential interactors. Acetylated peptide pulled down proteins associated with chromosome maintenance/segregation and cell cycle, opening the question of a possible link between H2B.Z acetylation status and mitosis.

Previously Unidentified Histone H1-Like Protein Is Involved in Cell Division and Ribosome Biosynthesis in Toxoplasma gondii

Chromatin dynamics can regulate all DNA-dependent processes. Access to DNA within chromatin is orchestrated mainly by histones and their posttranslational modifications (PTMs). Like other eukaryotes, the apicomplexan parasite Toxoplasma gondii encodes four canonical histones and five histone variants. In contrast, the linker histone (H1) has never been identified in apicomplexan parasites. In other eukaryotes, histone H1 compacts the chromatin by linking the nucleosome and increasing the DNA compaction. H1 is a multifunctional protein and can be involved in different steps of DNA metabolism or associated with protein complexes related to distinct biological processes. We have identified a novel protein in T. gondii ("TgH1-like") that, although lacking the globular domain of mammalian H1, is remarkably like the H1-like proteins of bacteria and trypanosomatids. Our results demonstrate that TgH1-like is a nuclear protein associated with chromatin and other histones. Curiously, TgH1-like is also in the nucleolus and associated with ribosomal proteins, indicating a versatile function in this parasite. Although knockout of the tgh1-like gene does not affect the cell cycle, it causes endopolygeny and asynchronous division. Interestingly, mutation of posttranslationally modified amino acids results in defects in cell division like those in the Δtgh1-like mutant, showing that these sites are important for protein function. Furthermore, in the bradyzoite stage, this protein is expressed only in dividing parasites, reinforcing its importance in cell division. Indeed, the absence of TgH1-like decreases compaction of peripheral chromatin, confirming its role in the chromatin modulation in T. gondii. IMPORTANCE Histone H1, or linker histone, is an important protein that binds to the nucleosome, aiding chromatin compaction. Here, we characterize for the first time a linker histone in T. gondii, named TgH1-like. It is a small and basic protein that corresponds only to the C-terminal portion of the human H1 but is similar to histone H1 from trypanosomatids and bacteria. TgH1-like is located in the nucleus, interacts with nucleosome histones, and acts in chromatin structure and cell division. Our findings show for the first time the presence of a histone H1 protein in an apicomplexan parasite and will provide new insights into cell division and chromatin dynamics in T. gondii and related parasites.

Genomic insights into host and parasite interactions during intracellular infection by Toxoplasma gondii

To gain insights into the molecular interactions of an intracellular pathogen and its host cell, we studied the gene expression and chromatin states of human fibroblasts infected with the Apicomplexan parasite Toxoplasma gondii. We show a striking activation of host cell genes that regulate a number of cellular processes, some of which are protective of the host cell, others likely to be advantageous to the pathogen. The simultaneous capture of host and parasite genomic information allowed us to gain insights into the regulation of the T. gondii genome. We show how chromatin accessibility and transcriptional profiling together permit novel annotation of the parasite's genome, including more accurate mapping of known genes and the identification of new genes and cis-regulatory elements. Motif analysis reveals not only the known T. gondii AP2 transcription factor-binding site but also a previously-undiscovered candidate TATA box-containing motif at one-quarter of promoters. By inferring the transcription factor and upstream cell signaling responses involved in the host cell, we can use genomic information to gain insights into T. gondii's perturbation of host cell physiology. Our resulting model builds on previously-described human host cell signalling responses to T. gondii infection, linked to induction of specific transcription factors, some of which appear to be solely protective of the host cell, others of which appear to be co-opted by the pathogen to enhance its own survival.

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.