iTRAQ-Based Global Phosphoproteomics Reveals Novel Molecular Differences Between Toxoplasma gondii Strains of Different Genotypes

iTRAQ-Based Phosphoproteomic Analysis Exposes Molecular Changes in the Small Intestinal Epithelia of Cats after Toxoplasma gondii Infection

Toxoplasma gondii, an obligate intracellular parasite, has the ability to invade and proliferate within most nucleated cells. The invasion and destruction of host cells by T. gondii lead to significant changes in the cellular signal transduction network. One important post-translational modification (PTM) of proteins is phosphorylation/dephosphorylation, which plays a crucial role in cell signal transmission. In this study, we aimed to investigate how T. gondii regulates signal transduction in definitive host cells. We employed titanium dioxide (TiO2) affinity chromatography to enrich phosphopeptides in the small intestinal epithelia of cats at 10 days post-infection with the T. gondii Prugniuad (Pru) strain and quantified them using iTRAQ technology. A total of 4998 phosphopeptides, 3497 phosphorylation sites, and 1805 phosphoproteins were identified. Among the 705 differentially expressed phosphoproteins (DEPs), 68 were down-regulated and 637 were up-regulated. The bioinformatics analysis revealed that the DE phosphoproteins were involved in various cellular processes, including actin cytoskeleton reorganization, cell necroptosis, and MHC immune processes. Our findings confirm that T. gondii infection leads to extensive changes in the phosphorylation of proteins in the cat intestinal epithelial cells. The results of this study provide a theoretical foundation for understanding the interaction between T. gondii and its definitive host.

Phosphoproteomics Reveal New Candidates in Abnormal Spermatogenesis of Pseudomales in Cynoglossus semilaevis

Phosphorylation is a post-translational modification that contributes to versatile protein functions in spermatogenesis, and the variations they generate usually results in abnormal spermatogenesis or sperm dysfunction. The sex-reversal phenomenon exists in Chinese tongue sole under certain conditions such that individuals with a ZW genotype can acquire a male phenotype and are thus called pseudomales. Pseudomale tongue sole can reach sexual maturity but produce only Z-type sperm, and the Z sperm carries paternal epigenetic information. Whether phosphorylation plays a role in the sperm abnormality of pseudomales is unknown. In this study, a phosphoproteomic analysis was performed to compare protein phosphorylation profiles between pseudomale and male testes. Altogether, we identified 14,253 phosphopeptides matching with 4843 proteins, with 1329 differentially phosphorylated peptides corresponding to 1045 differentially phosphorylated proteins (DPPs). Phosphorylation at 781 sites was upregulated and at 548 sites was downregulated. Four motifs were identified among differentially phosphorylated peptides, which were "SP", "SD", "RxxS", and "TP". Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses suggested that the cell cycle and DNA/RNA processing were significantly enriched with the genes encoding DPPs. To analyze DPP function in depth, a protein-protein interaction network was constructed, and Ran-binding protein 2 was found to play a central role in spermatogenesis by regulating several processes such as the cell cycle, eukaryotic translation, ubiquitination, and minichromosome maintenance. In kinase-associated network analyses, two "mitogen-activated protein kinase (Mapk)-centered" clusters were identified that may account for abnormal spermatogenesis in pseudomales. One cluster was centered on Mapk6, which predominantly regulated the cell cycle by interacting with several cyclin-dependent kinases, and the other was centered on the "testis-expressed kinase 1-like (Tesk1l)/Pim1l-Mapk4l- testis-expressed 14 (Tex14)" kinase cascade, which might contribute to spermatogenesis by regulating β-catenin. Taken together, these data suggested the new candidates involved in pseudomale sperm abnormalities and provided clues to discover the phosphorylated regulatory mechanism underlying tongue sole spermatogenesis.

Transcriptome analysis of a newly established mouse model of Toxoplasma gondii pneumonia

BACKGROUND

Toxoplasmosis is a zoonotic parasitic disease caused by Toxoplasma gondii. Toxoplasma gondii infection of the lungs can lead to severe pneumonia. However, few studies have reported Toxoplasma pneumonia. Most reports were clinical cases due to the lack of a good disease model. Therefore, the molecular mechanisms, development, and pathological damage of Toxoplasma pneumonia remain unclear.

METHODS

A mouse model of Toxoplasma pneumonia was established by nasal infection with T. gondii. The model was evaluated using survival statistics, lung morphological observation, and lung pathology examination by hematoxylin and eosin (H&E) and Evans blue staining at 5 days post-infection (dpi). Total RNA was extracted from the lung tissues of C57BL/6 mice infected with T. gondii RH and TGME49 strains at 5 dpi. Total RNA was subjected to transcriptome analysis by RNA sequencing (RNA-seq) followed by quantitative real-time polymerase chain reaction (qRT-PCR) validation. Transcript enrichment analysis was performed using the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases to assess the biological relevance of differentially expressed transcripts (DETs).

RESULTS

C57BL/6 mice infected with T. gondii via nasal delivery exhibited weight loss, ruffled fur, and respiratory crackles at 5 dpi. The clinical manifestations and lethality of RH strains were more evident than those of TGME49. H&E staining of lung tissue sections from mice infected with T. gondii at 5 dpi showed severe lymphocytic infiltration, pulmonary edema, and typical symptoms of pneumonia. We identified 3167 DETs and 1880 DETs in mice infected with the T. gondii RH and TGME49 strains, respectively, compared with the phosphate-buffered saline (PBS) control group at 5 dpi. GO and KEGG enrichment analyses of DETs showed that they were associated with the immune system and microbial infections. The innate immune, inflammatory signaling, cytokine-mediated signaling, and chemokine signaling pathways displayed high gene enrichment.

CONCLUSION

In this study, we developed a new mouse model for Toxoplasma pneumonia. Transcriptome analysis helped to better understand the molecular mechanisms of the disease. These results provided DETs during acute T. gondii lung infection, which expanded our knowledge of host immune defenses and the pathogenesis of Toxoplasma pneumonia.

Effect of deleting four Toxoplasma gondii calcium-binding EGF domain-containing proteins on parasite replication and virulence

Several calcium-binding proteins including calcium-dependent protein kinases play important roles in several facets of the intracellular infection cycle of the apicomplexan protozoan parasite Toxoplasma gondii. However, the role of the calcium-binding epidermal growth factor (EGF) domain-containing proteins (CBDPs) remains poorly understood. In this study, we examined the functions of four CBDP genes in T. gondii RH strain of type I by generating knock-out strains using CRISPR-Cas9 system. We investigated the ability of mutant strains deficient in CBDP1, CBDP2, CBDP3, or CBDP4 to form plaques, replicate intracellularly, and egress from the host cells. The results showed that no definite differences between any of these four CBDP mutant strains and the wild-type strain in terms of their ability to form plaques, intracellular replication, and egress. Additionally, CBDP mutants did not exhibit any significant attenuated virulence compared to the wild-type strain in mice. The expression profiles of CBDP2-4 genes were conserved among T. gondii strains of different genotypes, life cycle stages, and developmental forms. Whether other CBDP genes play any roles in the pathogenicity of T. gondii strains of different genotypes remains to be elucidated.

Comparative phosphoproteomics reveal new candidates in the regulation of spermatogenesis of Drosophila melanogaster

The most common phenotype induced by the endosymbiont Wolbachia in insects is cytoplasmic incompatibility, where none or fewer progenies can be produced when Wolbachia-infected males mate with uninfected females. This suggests that some modifications are induced in host sperms during spermatogenesis by Wolbachia. To identify the proteins whose phosphorylation states play essential roles in male reproduction in Drosophila melanogaster, we applied isobaric tags for relative and absolute quantitation (iTRAQ)-based proteomic strategy combined with titanium dioxide (TiO₂ ) enrichment to compare the phosphoproteome of Wolbachia-infected with that of uninfected male reproductive systems in D. melanogaster. We identified 182 phosphopeptides, defining 140 phosphoproteins, that have at least a 1.2 fold change in abundance with a P-value of <0.05. Most of the differentially abundant phosphoproteins (DAPPs) were associated with microtubule cytoskeleton organization and spermatid differentiation. The DAPPs included proteins already known to be associated with spermatogenesis, as well as many not previously studied during this process. Six genes coding for DAPPs were knocked down, respectively, in Wolbachia-free fly testes. Among them, Slmap knockdown caused the most severe damage in spermatogenesis, with no mature sperm observed in seminal vesicles. Immunofluorescence staining showed that the formation of individualization complex composed of actin cones was completely disrupted. These results suggest that Wolbachia may induce wide changes in the abundance of phosphorylated proteins which are closely related to male reproduction. By identifying phospho-modulated proteins we also provide a significant candidate set for future studies on their roles in spermatogenesis.

Quantitative phosphoproteomics analyses reveal the regulatory mechanisms related to frozen-thawed sperm capacitation and acrosome reaction in yak (Bos grunniens)

Mammalian spermatozoa are not mature after ejaculation and must undergo additional functional and structural changes within female reproductive tracts to achieve subsequent fertilization, including both capacitation and acrosome reaction (AR), which are dominated by post-translational modifications (PTMs), especially phosphorylation. However, the mechanism of protein phosphorylation during frozen-thawed sperm capacitation and AR has not been well studied. In this study, the phosphoproteomics approach was employed based on tandem mass tag (TMT) labeling combined with liquid chromatography-tandem mass spectrometry (LC-MS/MS) strategy to analyze frozen-thawed sperm in Ashidan yak under three sequential conditions (density gradient centrifugation-based purification, incubation in the capacitation medium and induction of AR processes by the calcium ionophore A23187 treatment). The identification of 1,377 proteins with 5,509 phosphorylation sites revealed changes in phosphorylation levels of sperm-specific proteins involved in regulation of spermatogenesis, sperm motility, energy metabolism, cilium movement, capacitation and AR. Some phosphorylated proteins, such as AKAP3, AKAP4, SPA17, PDMD11, CABYR, PRKAR1A, and PRKAR2A were found to regulate yak sperm capacitation and AR though the cAMP/PKA signaling pathway cascades. Notably, the phosphorylation level of SPA17 at Y156 increased in capacitated sperm, suggesting that it is also a novel functional protein besides AKAPs during sperm capacitation. Furthermore, the results of this study suggested that the phosphorylation of PRKAR1A and PRKAR2A, and the dephosphorylation of CABYR both play key regulatory role in yak sperm AR process. Protein-protein interaction analysis revealed that differentially phosphorylated proteins (AKAP3, AKAP4, FSIP2, PSMD11, CABYR, and TPPP2) related to capacitation and AR process played a key role in protein kinase A binding, sperm motility, reproductive process, cytoskeleton and sperm flagella function. Taken together, these data provide not only a solid foundation for further exploring phosphoproteome of sperm in yak, but an efficient way to identify sperm fertility-related marker phosphorylated proteins.

Global profiling of protein lysine malonylation in Toxoplasma gondii strains of different virulence and genetic backgrounds

Lysine malonylation is a post-translational modification (PTM), which regulates many cellular processes. Limited information is available about the level of lysine malonylation variations between Toxoplasma gondii strains of distinct genetic lineages. Yet, insights into such variations are needed to understand the extent to which lysine malonylation contributes to the differences in the virulence and repertoire of virulence factors between T. gondii genotypes. In this study, we profiled lysine malonylation in T. gondii using quantitative liquid chromatography-tandem mass spectrometry (LC-MS/MS) and immuno-affinity purification. This analysis was performed on three T. gondii strains with distinctive pathogenicity in mice, including RH strain (type I), PRU strain (type II), and VEG strain (type III). In total, 111 differentially malonylated proteins and 152 sites were upregulated, and 17 proteins and 17 sites were downregulated in RH strain versus PRU strain; 50 proteins and 59 sites were upregulated, 50 proteins and 53 sites were downregulated in RH strain versus VEG strain; and 72 proteins and 90 sites were upregulated, and 7 proteins and 8 sites were downregulated in VEG strain versus PRU strain. Differentially malonylated proteins were involved in key processes, such as those mediating the regulation of protein metabolism, stress response, glycolysis, and actin cytoskeleton. These results reveal an association between lysine malonylation and intra-species virulence differences in T. gondii and offer a new resource for elucidating the contribution of lysine malonylation to energy metabolism and virulence in T. gondii.

Lysine malonylation has been shown to play important roles in various biological processes in Toxoplasma gondii. Here, we used quantitative liquid chromatography-tandem mass spectrometry (LC-MS/MS) and immuno-affinity purification to test the hypothesis that lysine malonylation underpins the inter-genotype differences in the virulence of T. gondii. Several up-regulated and down- regulated malonylated proteins were identified in the tachyzoites of RH (type I) strain, PRU (type II) strain, and VEG (type III) strain. Differentially regulated malonylated proteins were enriched in many biological and metabolic pathways, and were found to contribute T. gondii energy metabolism, stress response, and infectivity, suggesting the role of lysine malonylation in the regulation of T. gondii virulence. These findings expand our knowledge of lysine malonylation in T. gondii and provide more insight into the mechanisms mediating the virulence differences between T. gondii strains of different genotypes.

Comparative Phosphoproteomic Analysis of Sporulated Oocysts and Tachyzoites of Toxoplasma gondii Reveals Stage-Specific Patterns

Toxoplasma gondii is an obligate intracellular protozoan of severe threat to humans and livestock, whose life history harbors both gamic and apogamic stages. Chinese 1 (ToxoDB#9) was a preponderant genotype epidemic in food-derived animals and humans in China, with a different pathogenesis from the strains from the other nations of the world. Posttranslational modifications (PTMs) of proteins were critical mediators of the biology, developmental transforms, and pathogenesis of protozoan parasites. The phosphoprotein profiling and the difference between the developmental phases of T. gondii, contributing to development and infectivity, remain unknown. A quantitative phosphoproteomic approach using IBT integrated with TiO₂ affinity chromatography was applied to identify and analyze the difference in the phosphoproteomes between the sporulated oocysts and the tachyzoites of the virulent ToxoDB#9 (PYS) strain of T. gondii. A total of 4058 differential phosphopeptides, consisting of 2597 upregulated and 1461 downregulated phosphopeptides, were characterized between sporulated the oocysts and tachyzoites. Twenty-one motifs extracted from the upregulated phosphopeptides contained 19 serine motifs and 2 threonine motifs (GxxTP and TP), whereas 16 motifs identified from downregulated phosphopeptides included 13 serine motifs and 3 threonine motifs (KxxT, RxxT, and TP). Beyond the traditional kinases, some infrequent classes of kinases, including Ab1, EGFR, INSR, Jak, Src and Syk, were found to be corresponding to motifs from the upregulated and downregulated phosphopeptides. Remarkable functional properties of the differentially expressed phosphoproteins were discovered by GO analysis, KEGG pathway analysis, and STRING analysis. S8GFS8 (DNMT1-RFD domain-containing protein) and S8F5G5 (Histone kinase SNF1) were the two most connected peptides in the kinase-associated network. Out of these, phosphorylated modifications in histone kinase SNF1 have functioned in mitosis and interphase of T. gondii, as well as in the regulation of gene expression relevant to differentiation. Our study discovered a remarkable difference in the abundance of phosphopeptides between the sporulated oocysts and tachyzoites of the virulent ToxoDB#9 (PYS) strain of T. gondii, which may provide a new resource for understanding stage-specific differences in PTMs and may enhance the illustration of the regulatory mechanisms contributing to the development and infectivity of T. gondii.

Global Phosphoproteomics Analysis of IBRS-2 Cells Infected With Senecavirus A

Phosphorylation is a widespread posttranslational modification that regulates numerous biological processes. Viruses can alter the physiological activities of host cells to promote virus particle replication, and manipulating phosphorylation is one of the mechanisms. Senecavirus A (SVA) is the causative agent of porcine idiopathic vesicular disease. Although numerous studies on SVA have been performed, comprehensive phosphoproteomics analysis of SVA infection is lacking. The present study performed a quantitative mass spectrometry-based phosphoproteomics survey of SVA infection in Instituto Biologico-Rim Suino-2 (IBRS-2) cells. Three parallel experiments were performed, and 4,520 phosphosites were quantified on 2,084 proteins. Gene Ontology (GO) functional enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses showed that many phosphorylated proteins were involved in apoptosis and spliceosome pathways, and subcellular structure localization analysis revealed that more than half were located in the nucleus. Motif analysis of proteins with differentially regulated phosphosites showed that proline, aspartic acid, and glutamic acid were the most abundant residues in the serine motif, while proline and arginine were the most abundant in the threonine motif. Forty phosphosites on 27 proteins were validated by parallel reaction monitoring (PRM) phosphoproteomics, and 30 phosphosites in 21 proteins were verified. Nine proteins with significantly altered phosphosites were further discussed, and eight [SRRM2, CDK13, DDX20, DDX21, BAD, ELAVL1, PDZ-binding kinase (PBK), and STAT3] may play a role in SVA infection. Finally, kinase activity prediction showed 10 kinases’ activity was reversed following SVA infection. It is the first phosphoproteomics analysis of SVA infection of IBRS-2 cells, and the results greatly expand our knowledge of SVA infection. The findings provide a basis for studying the interactions of other picornaviruses and their mammalian host cells.

A global phosphoproteomics analysis of adult Fasciola gigantica by LC-MS/MS

Protein phosphorylation plays key roles in a variety of essential cellular processes. Fasciola gigantica is a tropical liver fluke causing hepatobiliary disease fascioliasis, leading to human health threats and heavy economic losses. Although the genome and protein kinases of F. gigantica provided new insights to understand the molecular biology and etiology of this parasite, there is scant knowledge of protein phosphorylation events in F. gigantica. In this study, we characterized the global phosphoproteomics of adult F. gigantica by phosphopeptide enrichment-based LC–MS/MS, a high-throughput analysis to maximize the detection of a large repertoire of phosphoproteins and phosphosites. A total of 1030 phosphopeptides with 1244 phosphosites representing 635 F. gigantica phosphoproteins were identified. The phosphoproteins were involved in a wide variety of biological processes including cellular, metabolic, and single-organism processes. Meanwhile, these proteins were found predominantly in cellular components like membranes and organelles with molecular functions of binding (51.3%) and catalytic activity (40.6%). The KEGG annotation inferred that the most enriched pathways of the phosphoproteins included tight junction, spliceosome, and RNA transport (each one contains 15 identified proteins). Combining the reports in other protozoa and helminths, the phosphoproteins identified in this work play roles in metabolic regulation and signal transduction. To our knowledge, this work performed the first global phosphoproteomics analysis of adult F. gigantica, which provides valuable information for development of intervention strategies for fascioliasis.

A functionally divergent transcription elongation factor 1-like protein in Toxoplasma gondii

Zinc ribbons, one of the largest fold groups among zinc fingers, often include proteins involved in the transcription machinery. Here, we identify and characterize one such zinc ribbon-bearing protein in the apicomplexan parasite Toxoplasma gondii, annotated as putative transcription elongation factor 1 (ELF1), with predicted functions in transcription and chromatin maintenance. We show that this ELF1 homolog, referred to as T. gondii ELF1-like divergent (TgELD), is expressed in both tachyzoite and bradyzoite developmental stages. TgELD associates with the cytoskeleton in the tachyzoites, while it transiently becomes a part of the cyst wall in the early bradyzoites, followed by a cytosolic and peripheral localization in late bradyzoites. TgELD is phosphorylated by a casein kinase 2-like protein, which has potential implications for its localization and function in the parasite.

Global phosphoproteome analysis reveals significant differences between sporulated oocysts of virulent and avirulent strains of Toxoplasma gondii

In this study, the differences in the phosphoproteomic landscape of sporulated oocysts between virulent and avirulent strains of Toxoplasma gondii were examined using a global phosphoproteomics approach. Phosphopeptides from sporulated oocysts of the virulent PYS strain (Chinese ToxoDB#9) and the avirulent PRU strain (type II) were enriched by titanium dioxide (TiO₂) affinity chromatography and quantified using IBT approach. A total of 10,645 unique phosphopeptides, 8181 nonredundant phosphorylation sites and 2792 phosphoproteins were identified. We also detected 4129 differentially expressed phosphopeptides (DEPs) between sporulated oocysts of PYS strain and PRU strain (|log1.5 fold change| > 1 and p < 0.05), including 2485 upregulated and 1644 downregulated phosphopeptides. Motif analysis identified 24 motifs from the upregulated phosphorylated peptides including 22 serine motifs and two threonine motifs (TPE and TP), and 15 motifs from the downregulated phosphorylated peptides including 12 serine motifs and three threonine motifs (TP, RxxT and KxxT) in PYS strain when comparing PYS strain to PRU strain. Several kinases were consistent with motifs of overrepresented phosphopeptides, such as PKA, PKG, CKII, IKK, MAPK, EGFR, INSR, Jak, Syk, Src, Ab1. GO enrichment, KEGG pathway analysis and STRING analysis revealed DEPs significantly enriched in many biological processes and pathways. Kinase related network analysis showed that AGC kinase was the most connected kinase peptide. Our findings reveal significant difference in phosphopeptide profiles of sporulated oocysts between virulent and avirulent T. gondii strains, providing new resources for further elucidation of the mechanisms underpinning the virulence of T. gondii.

iTRAQ-Based Phosphoproteomic Analysis of Toxoplasma gondii Tachyzoites Provides Insight Into the Role of Phosphorylation for its Invasion and Egress

The invasion and egress are two key steps in lytic cycle vital to the propagation of Toxoplasma gondii infection, and phosphorylation is believed to play important roles in these processes. However, the phosphoproteome of T. gondii at these two stages has not been characterized. In this study, we profiled the phosphoproteome of tachyzoites at the stages of “just invading” (JI) and “prior to egress” (PE) based on iTRAQ quantitative analysis, in which a total of 46 phosphopeptides, 42 phosphorylation sites, and 38 phosphoproteins were detected. In the comparison of PE vs. JI, 10 phosphoproteins were detected with their phosphorylation level significantly changed, and four of them were demonstrated to be significantly down-regulated at the transcriptional level. Bioinformatic analysis of these identified phosphoproteins suggested that phosphorylation-mediated modulation of protein function was employed to regulate the pathway of toxoplasmosis and metabolism and cellular processes correlated with tachyzoite’s binding, location, and metabolism, and thus play vital roles in the parasite lytic cycle. Moreover, cytoskeletal network (CN)-associated Inner Membrane Complex (IMC1, IMC4, IMC6 and IMC12), Intravascular Network (IVN)-related GRAs (GRA2, GRA3, GRA7 and GRA12), and Parasitophorous Vacuole Membrane (PVM)-localized ROP5 were shown to be enriched at the central nodes in the protein interaction network generated by bioinformatic analysis, in which the phosphorylation level of IMC4, GRA2, GRA3, and GRA12 were found to be significantly regulated. This study revealed the main cellular processes and key phosphoproteins crucial for the invasion and egress of T. gondii, which will provide new insights into the developmental biology of T. gondii in vitro and contribute to the understanding of pathogen-host interaction from the parasite perspective.

Toxoplasma gondii GRA60 is an effector protein that modulates host cell autonomous immunity and contributes to virulence

Toxoplasma gondii infects virtually any nucleated cell and resides inside a non-phagocytic vacuole surrounded by a parasitophorous vacuolar membrane (PVM). Pivotal to the restriction of T. gondii dissemination upon infection in murine cells is the recruitment of immunity regulated GTPases (IRGs) and guanylate binding proteins (GBPs) to the PVM that leads to pathogen elimination. The virulent T. gondii type I RH strain secretes a handful of effectors including the dense granule protein GRA7, the serine-threonine kinases ROP17 and ROP18, and a pseudo-kinase ROP5, that synergistically inhibit the recruitment of IRGs to the PVM. Here, we characterise GRA60, a novel dense granule effector, which localises to the vacuolar space and PVM and contributes to virulence of RH in mice, suggesting a role in the subversion of host cell defence mechanisms. Members of the host cell IRG defence system Irgb10 and Irga6 are recruited to the PVM of RH parasites lacking GRA60 as observed previously for the avirulent RHΔrop5 mutant, with RH preventing such recruitment. Deletion of GRA60 in RHΔrop5 leads to a recruitment of IRGs comparable to the single knockouts. GRA60 therefore represents a novel parasite effector conferring resistance to IRGs in type I parasites, and found associated to ROP18, a member of the virulence complex.

Label-Free Quantitative Acetylome Analysis Reveals Toxoplasma gondii Genotype-Specific Acetylomic Signatures

Distinct genotypic and pathogenic differences exist between Toxoplasma gondii genotypes. For example, genotype I is highly virulent, whereas genotype II and genotype III are less virulent. Moreover, Chinese 1 genotype (ToxoDB#9) is also virulent. Here, we compare the acetylomes of genotype 1 (RH strain) and Chinese 1 genotype (ToxoDB#9, PYS strain) of T. gondii. Using mass spectrometry enriched for acetylated peptides, we found a relationship between the levels of protein acetylation and parasite genotype-specific virulence. Notably, lysine acetylation was the largest (458 acetylated proteins) in RH strain, followed by PYS strain (188 acetylated proteins), whereas only 115 acetylated proteins were detected in PRU strain. Our analysis revealed four, three, and four motifs in RH strain, PRU strain and PYS strain, respectively. Three conserved sequences around acetylation sites, namely, xxxxxK^AcHxxxx, xxxxxK^AcFxxxx, and xxxxGK^AcSxxxx, were detected in the acetylome of the three strains. However, xxxxxK^AcNxxxx (asparagine) was found in RH and PYS strains but was absent in PRU strain. Our analysis also identified 15, 3, and 26 differentially expressed acetylated proteins in RH strain vs. PRU strain, PRU strain vs. PYS strain and PYS strain vs. RH strain, respectively. KEGG pathway analysis showed that a large proportion of the acetylated proteins are involved in metabolic processes. Pathways for the biosynthesis of secondary metabolites, biosynthesis of antibiotics and microbial metabolism in diverse environments were featured in the top five enriched pathways in all three strains. However, acetylated proteins from the virulent strains (RH and PYS) were more enriched in the pyruvate metabolism pathway compared to acetylated proteins from PRU strain. Increased levels of histone-acetyl-transferase and glycyl-tRNA synthase were detected in RH strain compared to PRU strain and PYS strain. Both enzymes play roles in stress tolerance and proliferation, key features in the parasite virulence. These findings reveal novel insight into the acetylomic profiles of major T. gondii genotypes and provide a new important resource for further investigations of the roles of the acetylated parasite proteins in the modulation of the host cell response to the infection of T. gondii.