ABSTRACT: To gain insights into differences in the virulence among T. gondii strains at the post-translational level, we conducted a quantitative analysis of the phosphoproteome profile of T. gondii strains belonging to three different genotypes. Phosphopeptides from three strains, type I (RH strain), type II (PRU strain) and ToxoDB#9 (PYS strain), were enriched by titanium dioxide (TiO2) affinity chromatography and quantified using iTRAQ technology. A total of 1,441 phosphopeptides, 1,250 phosphorylation sites and 759 phosphoproteins were detected. In addition, 392, 298, and 436 differentially expressed phosphoproteins (DEPs) were identified in RH strain when comparing RH/PRU strains, in PRU strain when comparing PRU/PYS strains, and in PYS strain when comparing PYS/RH strains, respectively. Functional characterization of the DEPs using GO, KEGG, and STRING analyses revealed marked differences between the three strains. In silico kinase substrate motif analysis of the DEPs revealed three (RxxS, SxxE, and SxxxE), three (RxxS, SxxE, and SP), and five (SxxE, SP, SxE, LxRxxS, and RxxS) motifs in RH strain when comparing RH/PRU strains, in PRU strain when comparing PRU/PYS, and in PYS strain when comparing PYS/RH strains, respectively. This suggests that multiple overrepresented protein kinases including PKA, PKG, CKII, IKK, and MAPK could be involved in such a difference between T. gondii strains. Kinase associated network analysis showed that ROP5, ROP16, and cell-cycle-associated protein kinase CDK were the most connected kinase peptides. Our data reveal significant changes in the abundance of phosphoproteins between T. gondii genotypes, which explain some of the mechanisms that contribute to the virulence heterogeneity of this parasite.