A study on combined methylation, phosphorylation, and acetylation of proteins done in a lung cancer model showed a large number of proteins that presented the three PTMs [
120,
121]. Taking a functional perspective, PTMs increase the diversity of functional units of protein origin within a cell. PTMs are still a field with great potential for future exploration in
T. gondii. A crosstalk study with ubiquitination and the other PTMs showed that less than 10% of the proteins overlapped between ubiquitin and arginine methylation, 21% combined acetylation with ubiquitin, 25% of the SUMO proteome was ubiquitinated, while 78% of the phosphorylated proteins are also ubiquitinated [
118]. Among the proteins detected that present combinations of different PTMs, the histone acetyltransferase EP300 and the chaperone Hsp90 [
120] stand out.
4.3.1. Histones
Over several decades of research, it has become evident that various PTMs on histone proteins engage in intricate crosstalk, ultimately influencing the regulation of gene expression. This intricate interplay of histone modifications has been characterized as a "histone code," which specific proteins can interpret, providing in conjunction with the substitution of different histone variants, a means for epigenetic regulation [
122]. Understanding the roles of PTMs is also complicated by the fact that the same writers and readers of marks also write and read PTMs on other proteins [
123].
Protozoan parasites, such as
T. gondii, are no exception to this trend. The use of high-throughput mass spectrometry is enabling the systematic dissection of histone modifications and core markers involved in epigenetic regulation [
124]. Numerous residues have been identified for modification, particularly within the histone tails. However, it is worth noting that various residues within the histone globular domains can also undergo diverse modifications, encompassing processes such as acetylation, phosphorylation, ubiquitination, methylation, and lately crotonylation and 2-hydroxyisobutyrylation [
124]. These last newly described marks have demonstrated their crucial significance in regulating global transcription in mammalian cells. Consequently, in
T. gondii, they may exert a noteworthy influence on the transcription process [
124].
There is a high frequency of lysine residues in histones, and these residues are often modified by PTMs [
86]. Among the histones, H3 and H4 PTMs are widely distributed not only in unicellular organisms, like protozoan parasites, but also in multicellular eukaryotes, and many of the amino acid residues that carry certain PTMs are highly conserved [
123].
Figure 3A shows the different histone marks and their crosstalk in the mammal model as described. Furthermore, histone marks identified in
T. gondii PTM proteomic analysis were highlighted to compare (
Figure 3B). In the figure, it can also see the conservation of possible crosstalk similar to those observed in mammals in punctuated lines, although some new marks are also detected and the crosstalks are marked as questions. As it can be observed, almost all the marks are observed, and although not much research has been conducted in this area in the parasite, their crosstalk can be inferred. As shown in the figure, in
T. gondii, the N-terminal tails of H4 and H3 are highly acetylated and methylated regions [
86]. Interestingly, some of the residues that are acetylated can also be methylated, like lysine 12 (H4K12me) or lysine 16 when not acetylated, was found trimethylated (H4K16me3) [
86]. Another example is H3K9, which was also found methylated and acetylated. While mono methylation and acetylation are active marks, di or tri methylation is often repressive and found in subtelomeres and centromeres [
125]. Acetylation on H4K31 has also been reported in
T. gondii and
P. falciparum and this residue can also be mono-methylated in a mutually exclusive manner [
115,
126]. In
T. gondii, this residue was found to be acetylated at the promoter of nearby active genes associated with H3K4me3 and H3K9ac marks (
Figure 3B). At the same time, it was mono-methylated in the core body of the gene, and it inversely correlated with gene expression, it would be a repressive mark as opposed to H4K31ac [
126]. Some PTMs not found in other species for this histone were found, like H4R23me, in the globular domain and methylation, both di and monomethylation on the C-terminal tail, on R78 [
86].
Another important modification on lysine residues is ubiquitination, but few PTMs occur on the same lysine residues as ubiquitin. One exception is H3K24, which is modified by ubiquitin, acetylation, mono- and trimethylation, implicating it as an important regulatory site [
118]. Apart from this residue, only H3K116 was found to be ubiquitinylated, although other modifications occur at substoichiometric levels only detectable using an enrichment strategy [
118].
There are also few PTMs in H2A.X, but two were found in the C-terminal tail, acetylation of K128 that has not been described, and the phosphorylated S132 [
86], which corresponds to S139p in humans with a conserved function for TgH2A.X in response to DNA damage [
127]. By contrast, H2A.Z displays several PTMs, of which many acetylations in the N-terminal tail stand out. Within the first 40 amino acids, there are 10 lysines (K6, K10, K14, K18, K24, K27, K29, K34, K36, and K37), all of which are acetylated. Additionally, lysine K18 could also be methylated, leading to various PTM combinations with acetylated residues [
86].
A similar situation arises for H2B histones, as canonical H2B possesses only one acetylation on the N-terminus at K4, whereas H2B.Z exhibits multiple modifications in the N-terminal tail, including acetylation on lysines 3, 8, 13, 14, and 18 [
86]. Another study using H2BK120 ubiquitin-specific antibody detected reactivity with
T. gondii H2B and H2B.Z supporting the presence of ubiquitin-histone conjugates in
T. gondii chromatin that could have conserved functions [
118]. Histone H2B.Z is unique to apicomplexan parasites and has been described to conform to a double variant nucleosome (DVN), hyperacetylated in the N-terminal tail and associated with transcriptional activation (
Figure 3C) [
39,
125,
127,
128,
129,
130]. Recent research has examined acetylations in the H2B.Z N-terminal tail lysines and their significance in various biological processes in
T. gondii [
39]. This study revealed that the inability to regulate this N-tail positive charge patch generated when no acetylation is possible, produced a reduced in vitro replication, a heightened differentiation rate, an increased sensitivity to DNA damage, and in an in vivo model produced a complete loss of
T. gondii virulence [
39].
As mentioned above, not much research has been conducted to date on histone crosstalks in
T. gondii. Nevertheless, some conclusions have been made, taking into account results when one mark can be sufficient to elicit a specific biological output while in other cases, multiple marks are required [
131] (
Figure 3B). For example, it has been proposed that gene activation needs the dual signature of H3R17me2 and H3K18 [
132]. In this work, the
T. gondii arginine methyltransferase TgCARM1 was shown to work in concert with the acetylase TgGCN5-A, whose substrate preference is H3K18. Another example is the redundancy or multiplicity of sumoylation sites in the parasite histones suggested by Bougdour et al. [
131].
Also, the transcription-associated PTMs H3K4ac, H3K9ac, and H3K27ac were found to be associated with H4K31ac in euchromatic regions and opposed to the heterochromatic regions revealed by the repressive marks H3K9me3 and H4K20me3 [
126] (
Figure 3B). In another study, the DVN conformed by H2A.Z and H2B.Z, which are both hyperacetylated in the N-terminal tails, was found in association with the activatory mark H3K4me3 [
129]. In turn, H3K4me3 was found along with H3K9ac in most of the
T. gondii genes that are expressed (
Figure 3C) [
131].
Another well-characterized example first in yeast but then in other eukaryotes is the requirement of histone H2B monoubiquitination for proper H3K4 and H3K79 methylation [
122,
133,
134]. Although ubiquitination was detected on TgH2AK119Ub and TgH2BK120Ub [
118], to date it has only been speculated that acetylation and ubiquitination may regulate the differential localization of H2A.Z and H2B.Z on active and silent genes in
T. gondii [
125] and so it is represented as a question in
Figure 3B.
4.3.2. Hsp90 PTM Crosstalk
The Hsp90 chaperone was shown to be simultaneously acetylated and phosphorylated in the lung cancer model [
120]. This protein is central since its main functions are to assist important proteins in biological processes such as DNA replication, regulation of gene expression, proliferation, etc.; in both stressed and non-stressed cells. The Hsp90 protein network includes hundreds of proteins, mainly transcription factors, and kinases, but also proteins associated with chromatin, metabolism, translation, and DNA damage [
135,
136,
137]. Hsp90 can form different complexes, one of them called the Hsp70/Hsp90 cycle. It may include other co-chaperones, the hsp40 type I being the one that initiates the recruitment of the client protein. The Hsp70/Hsp90 cycle is conserved in
T. gondii [
137]. Recently, the proteome of Hsp90 and Hsp40 type I in
T. gondii was analyzed and referred to as Tgj1 [
62]. This study allowed defining a putative PPI network for both chaperones as part of the Hsp70/Hsp90 cycle and in their independent functions between them. Within the Hsp70/Hsp90 cycle, the enriched pathways were translation, cell redox homeostasis, and protein folding. The PPI associated with Hsp90 not related to the Tgj1-Hsp90 axis showed mainly interactors related to protein folding, RNA processing, cell signaling, and transcription. This shows the important role of this chaperone in a wide diversity of biological processes of
T. gondii. The proteomes related to PTMs show that
T. gondii Hsp90 presents phosphorylations, acetylations, ubiquitinylation, and even a monomethyl arginine modification (ToxoDB_TGME49_288380). Backe et al. describe the role of the different PTMs in Hsp90 in the review [
138]. Based on the role of the PTMs present in Hsp90 of other organisms, a proposed role for the
T. gondii Hsp90 PTMs counterpart can be deduced [
61]. Briefly, phosphorylation of Hsp90 in yeast and other species may affect its binding to certain client proteins or even its affinity for ATP. In other cases (e.g., Hsp90 T101), it may promote kinase client activation. Hsp90 phosphorylation in turn is associated with the cell cycle and in some cases promotes cell proliferation [
139,
140]. Therefore, these PTMs could similarly affect
T. gondii Hsp90. Phosphorylation of threonine 5 and 7 on Hsp90 are associated with double-strand break repair [
141], but these modifications were not detected in
T. gondii. The phosphorylation of residues S231 (T220 in
T. gondii) and S263 is related to telomerase activity [
142]. Similarly, acetylation of Hsp90 leads to a loss of affinity for ATP and a decrease in its binding to the client protein and co-chaperones [
138].
T. gondii presents two acetylations in the acetylome (K384 and K559). Acetylation of lysine K384 in human Hsp90 affects the interaction with its client protein, the receptor tyrosine kinase ErbB2 [
143]. Acetylation at K558 (TgK559) could be associated with the export of Hsp90 to the extracellular space [
144]. The ubiquitinylation of Hsp90 seems to be associated only with its degradation. CHIP is one of the co-chaperones associated with the Hsp70/Hsp90 cycle, which negatively regulates the pathway. This CHIP protein is an E3 ubiquitin ligase and binds to Hsp90. However, phosphorylated forms of Hsp90 decrease their association with CHIP, suggesting a crosstalk between both PTMs [
145].