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

To gain insights into differences in the virulence among strains at the post-translational level, we conducted a quantitative analysis of the phosphoproteome profile of 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. 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 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 genotypes, which explain some of the mechanisms that contribute to the virulence heterogeneity of this parasite.