RESEARCH ARTICLE Temporal and thermal profiling of the Toxoplasma proteome implicates parasite Protein Phosphatase 1 in the regulation of Ca2+- responsive pathways Alice L Herneisen1,2, Zhu- Hong Li3, Alex W Chan1,3, Silvia NJ Moreno3, Sebastian Lourido1,2* 1Whitehead Institute for Biomedical Research, Cambridge, United States; 2Biology Department, Massachusetts Institute of Technology, Cambridge, United States; 3Center for Tropical and Emerging Global Diseases, University of Georgia, Athens, United States Abstract Apicomplexan parasites cause persistent mortality and morbidity worldwide through diseases including malaria, toxoplasmosis, and cryptosporidiosis. Ca2+ signaling pathways have been repurposed in these eukaryotic pathogens to regulate parasite- specific cellular processes governing the replicative and lytic phases of the infectious cycle, as well as the transition between them. Despite the presence of conserved Ca2+- responsive proteins, little is known about how specific signaling elements interact to impact pathogenesis. We mapped the Ca2+- responsive proteome of the model apicomplexan Taxoplasma gondii via time- resolved phosphoproteomics and thermal proteome profiling. The waves of phosphoregulation following PKG activation and stimulated Ca2+ release corroborate known physiological changes but identify specific proteins operating in these pathways. Thermal profiling of parasite extracts identified many expected Ca2+- responsive proteins, such as parasite Ca2+- dependent protein kinases. Our approach also identified numerous Ca2+- responsive proteins that are not predicted to bind Ca2+, yet are critical components of the parasite *For correspondence: signaling network. We characterized protein phosphatase 1 (PP1) as a Ca2+- responsive enzyme that lourido@wi.mit.edu relocalized to the parasite apex upon Ca2+ store release. Conditional depletion of PP1 revealed that Competing interest: See page the phosphatase regulates Ca2+ uptake to promote parasite motility. PP1 may thus be partly respon- 32 sible for Ca2+-r egulated serine/threonine phosphatase activity in apicomplexan parasites. Funding: See page 33 Preprinted: 25 May 2022 Received: 17 May 2022 Editor's evaluation Accepted: 17 August 2022 Herneisen et al., provide a comprehensive and thorough exploration of Ca 2+ responsive changes in Published: 17 August 2022 the Toxoplasma proteome and the resulting phosphorylation events during the transition from intra- cellular residing parasites to egress from the host cell. Furthermore, a novel temperature stability Reviewing Editor: Dominique profiling method of all proteins responding to Ca2+ concentration with a change in stability is a novel Soldati- Favre, University of Geneva, Switzerland applicable tool that here is used to map Ca 2+-r esponsive proteins in the parasites. They provide a compelling analysis of the complex data and carefully validate their findings using genetics and cell Copyright Herneisen et al. biology. This work is of the highest quality in the field. This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that Introduction the original author and source Apicomplexan parasites cause persistent mortality and morbidity worldwide through diseases are credited. including malaria, toxoplasmosis, and cryptosporidiosis (Havelaar et al., 2015). The phylum member Herneisen et al. eLife 2022;11:e80336. DOI: https://doi.org/10.7554/eLife.80336 1 of 41 Research article Cell Biology | Microbiology and Infectious Disease Toxoplasma gondii alone infects >2 billion people. As obligate intracellular pathogens, apicomplexans are exquisitely tuned to transduce environmental signals into programs of motility, replication, and quiescence responsible for parasite pathogenesis and spread (Bisio and Soldati- Favre, 2019). Ca2+ signals and their downstream effectors are a part of the signaling cascade that impacts almost every cellular function (Lourido and Moreno, 2015). Signaling begins with release of Ca2+ from intracellular stores or influx through plasma membrane channels, resulting in diverse downstream events central to parasite virulence, including secretion of adhesive proteins, motility, and invasion into and egress from host cells. Together, these cellular processes orchestrate a dramatic transition from the replicative to the kinetic phase of the life cycle that allows parasites to spread to new host cells. Signal-t ransducing components downstream of Ca2+ release are largely unknown yet are likely essential for apicomplexan viability and virulence (Lourido and Moreno, 2015; Nagamune and Sibley, 2006). Ca2+ can change a protein’s state by direct binding or indirect effects such as triggering post- translational modification (PTM), interaction with other proteins, or relocalization. Indirect effects are fundamental to the propagation and amplification of signals across the Ca2+- regulated network, and in most organisms, they are largely mediated by three classes of Ca2+-b inding proteins: Ca2+-r egulated kinases, Ca2+- regulated phosphatases, and calmodulin (CaM) and related proteins (Villalobo et al., 2019). Genomic searches for canonical Ca2+- binding domains in apicomplexans have identified several individual proteins involved in transducing and effectuating Ca2+ signals (Farrell et al., 2012; Huet et al., 2018; McCoy et al., 2017), like kinases, phosphatases, and transporters (Hortua Triana et al., 2018; Lourido et al., 2012; Lourido et al., 2010; Luo et al., 2005; Márquez- Nogueras et al., 2021). However, many of the key signaling elements involved in fundamental Ca2+ responses—including the channels responsible for its stimulated release—are either missing from apicomplexan genomes or have diverged beyond recognition, suggesting that eukaryotic pathogens evolved novel pathways for Ca2+ mobilization and transduction (Billker et al., 2009; Lourido and Moreno, 2015). In apicomplexans, Ca2+- dependent protein kinases (CDPKs) have garnered the most attention as the only known Ca2+- regulated kinases in the phylum (Billker et al., 2009). These kinases possess intrinsic Ca2+- binding sites and do not rely on CaM- like mammalian CaMKs. In T. gondii, several of these CDPKs trigger parasite motility (Lourido et  al., 2012; Lourido et  al., 2010; McCoy et  al., 2012; Smith et al., 2022; Treeck et al., 2014). The roles of TgCDPK7 in replication and TgCDPK2 in amylopectin granule formation also implicate Ca2+ signaling in cellular functions outside the kinetic phase of the parasite lytic cycle (Morlon-G uyot et  al., 2014; Bansal et  al., 2021; Uboldi et  al., 2015). Dephosphorylation has garnered comparatively little attention in these parasites (Yang and Arrizabalaga, 2017). The roles of the prototypical CaM and the Ca2+/CaM- dependent phosphatase calcineurin have only been phenotypically examined, and their client proteins remain largely unknown (Paul et al., 2015; Philip and Waters, 2015). Although key players are conserved and essential across the Apicomplexa, no systematic efforts have been undertaken to globally map the Ca2+ signaling pathways of these pathogens. Ca2+ signaling pathways have been repurposed in apicomplexans to regulate parasite-s pecific cellular processes governing the transition between the replicative and kinetic phases of the infectious cycle (Brown et al., 2020; Pace et al., 2020). Despite the presence of conserved Ca2+-r esponsive proteins (Lourido and Moreno, 2015; Nagamune and Sibley, 2006), uncovering the Ca2+ signaling architecture of apicomplexans demands a reevaluation of the entire network to understand how specific signaling elements interact to impact pathogenesis. We present an atlas of Ca2+- regulated proteins in the model apicomplexan T. gondii, assembled from high-d imensional proteomic datasets. The physiological changes associated with stimulated motility in the asexual stages of parasites have been characterized for decades (Blader et al., 2015). Our approach identified at once hundreds of molecular components underpinning these processes. We find numerous Ca2+-r esponsive proteins that are not predicted to bind Ca2+, yet operate at critical junctures in the parasite signaling network. From this analysis, the protein phosphatase PP1 emerges as an unanticipated Ca2+-r esponsive protein. Herneisen et al. eLife 2022;11:e80336. DOI: https://doi.org/10.7554/eLife.80336 2 of 41 Research article Cell Biology | Microbiology and Infectious Disease Results Sub-minute phosphoproteomics reveals the topology of Ca2+- dependent signaling processes We can emulate the endogenous signaling pathways that mediate the kinetic phase of the lytic cycle by treating isolated parasites with zaprinast, which stimulates PKG activation by inhibiting phospho- diesterases that degrade cGMP (Brown et  al., 2016; Lourido et  al., 2012). Zaprinast-s timulated motility occurs rapidly and in defined sequence in apicomplexans: initially, an increase in cGMP acti- vates parasite protein kinase G (PKG), which phosphorylates substrates and stimulates Ca2+ release from internal stores (Brown et al., 2020; Lourido and Moreno, 2015). PKG likely performs functions that extend beyond regulating Ca2+ stores, for example by mobilizing diacylglycerol and phosphatidic acid (Lourido et al., 2012; Brown et al., 2017; Bullen et al., 2016; Bisio et al., 2019); however, the use or phosphodiesterase inhibitors like zaprinast allows us to stimulate endogenous Ca2+ release without flooding the cell with Ca2+, as is the case with ionophores. Ca2+- dependent protein kinases, such as TgCDPK1 and TgCDPK3, synergize with PKG to effectuate microneme secretion and parasite motility (Brown et al., 2016; Lourido et al., 2012; Lourido et al., 2010; McCoy et al., 2012). This signaling cascade is active within seconds of cGMP elevation; however, existing T. gondii phosphopro- teomes compare changes at a single time point following Ca2+- ionophore stimulation (Treeck et al., 2014), providing only a snapshot of diverging signaling states. Here, we add kinetic resolution to these signaling pathways. We quantified dynamic changes in the phosphoproteome within a minute of stimulation with zaprinast and thus activation of the cGMP/Ca2+ pathway. We collected five timepoints in the 60 s following stimulation (0, 5, 10, 30, and 60 s), as well as three DMSO- treated matched timepoints (0, 10, and 30 s), in biological duplicate (Figure 1A). Using TMTpro labeling methods (Li et  al., 2020), we multiplexed 16  samples, allowing us to analyze a complete time course, with replicates and controls, in a single MS experiment. Phosphopeptides were enriched from the rest of the sample using sequential metal- oxide affinity chromatography, which maximizes phosphopeptide capture (Tsai et al., 2014). Our experiments quantified 4,055 parasite proteins, none of which exhibited more than a twofold change in abundance in the 60 s following stimulation (Figure 1—figure supplement 1). Given the paucity of known phosphoreglatory interactions in apicomplexans compared to other organisms (Weiss et al., 2020), we employed several analysis approaches to maximize the identifi- cation of changing phosphosites. We first calculated phosphoregulation scores by summing peptide abundances of vehicle (DMSO) and zaprinast-t reated samples, taking their ratios, and standardizing the values with a modified Z score (Figure 1B). From a phosphoproteome of 11,755 unique peptides with quantification values (belonging to 2,690 phosphoproteins), 839 phosphopeptides increased in abundance three modified Z scores above the median, whereas 154 decreased 1.8 Z scores below the median. Principal component analysis on the significant peptides distinguished the agonist treatment and time- course kinetics in the two principal components accounting for the greatest variability in the data (Figure 1C and Figure 1—figure supplement 1). Kinetically resolved clusters reveal regulatory subnetworks during zaprinast stimulation We leveraged the kinetic resolution of our comprehensive phosphoproteomics datasets to identify subregulatory networks. A Gaussian mixture-m odel clustering algorithm (Invergo et al., 2017) heuris- tically resolved four clusters for phosphopeptides arising from zaprinast treatment: three clusters increasing with different kinetics, and one decreasing (Figure 1D). On average, the 173 phosphopep- tides belonging to cluster 1 increased sharply in abundance within 5 s of treatment and continued to increase for the remainder of the time course (Figure 1D), suggesting that they belonged to the first wave of phosphoregulation. This cluster was enriched for phosphoproteins associated with phos- phodiesterase activity, phospholipid binding, and Ca2+ binding (Figure 1E), including PDE1, PDE2, PI4,5K, PI3,4K, PI- PLC, a phosphatidylinositol- 3,4,5-t riphosphate 5-p hosphatase, a putative Sec14, TgCDPK2A and TgCDPK7, and PPM2B (Table 1 and Figure 1F). Peptides in clusters 2 and 3 (173 and 527, respectively) increased more gradually and exhibited lower fold- changes than cluster 1 sites (Figure 1D). Cluster 2 was notably enriched in proteins func- tioning in transport of monovalent ions and lipids, cyclase activity, and Ca2+ binding (Figure 1E). This Herneisen et al. eLife 2022;11:e80336. DOI: https://doi.org/10.7554/eLife.80336 3 of 41 Research article Cell Biology | Microbiology and Infectious Disease A B C DMSO zaprinast DMSO MS TMT16plex replicate 1 zaprinast 3 replicate 2 SMOAC proteome 2 time (s): 0 5 10 30 60 phosphoproteome 2 3 summed abundance vehicle PC1 (72.5%) D F G 2 2 1 PDE1 n 32 1 12 1 1 PDE21 2 3 2 ATPase2 PDE2 GC AAP2 SCE1 3 KinesinB 6 3 GAPM3 AAP54 4 ERK7 2 CGP RNG2 4 TKL2 1 2 3 seconds post-zaprinast treatment 5time (s) clusters represented E 411 434 5 509 12 81 87 61 6 3 75 69 9 61 0 4 3 26 96 88 24 75 62 29 90 95 78 25 73 89 15 0700 05 05 04 08 002 52 52 820 0 0 099 152 8 8 3 0 3 3 0 6 8 8 6 5 6 8 5 2 70 0 0 0 0 0 0 0 0 16 22 40 05 05 08 15 15 19 22 30 42 42 46 60 05 04 –log (p) gene ratio O: :0 :0 :0 :0 0 0 0 0O O O O O: O: O: O: O:0 0 :00 :00 :00 :01 :00 :00 :00 :00 :00 :00 :00 :00 :00 :00 :00 :00 :00 :00 G G G G G G G G G G GO GO GO GO GO GO GO GO GO GO GO GO GO GO GO GO GO GO 2.5 1 2 3 1.5 4 as e ing ng se se ng e l l i i n ne se se e l se el ter or m er ter m er er or se er er or ity ity ste r nd nd era a t i i t rol ind an anh h tid a cla n n t an r lya an o iva i s rt r is rt rt t t t e b b s d b c c p cy ch n ch sp ct ha n pos p o an o o ula ol a po r po r ula cti v cti v di id n o die hy n o n + pe m yg e d ran e a ec ran ran s ch ns p p ns eg yd r ns ns eg r a a ho p r r e p oli 2+ i ho ter 2+ i tio Ka pe siu -ox at e t s m t t me tra tra se h a a e s h a p s g c y s g id Pa ve on ion ve on n a te r r r t s n t n t se or ina ho p s ie -t a s li p T ti ti ti ti tio TP es io l io at a ns p k e p ho s C ho d M ine o t oru ed G yla ca yla ca ca G c d a h a AP d p p ic e p h l r r d i t pic or ple e os tr M oti tid e or st d p p ho ho le n h u m h le eo sp h cy ate s uiv ho -c o sp spp u p ga sp -co rip uc cl ho ac t as e ph o ho -co o r ho p e in p as e e-t c-n u p + - r, r, s P sid ycl i ic- n a2 AT P rte rte Pa AT leo c yc l C po po ATs nu c ',5 '- c ran n s 3 2+ t n t ra Ca io Figure 1. Phosphoregulation triggered by Ca2+ release. (A) Schematic of the sub-m inute phosphoproteomics experiments with the Ca2+ signaling agonist zaprinast. (B) The summed abundances of unique phosphopeptides during zaprinast or vehicle (DMSO) treatment. The abundance ratios were transformed into a modified Z score and were used to threshold increasing (Z>3; blue) or decreasing (Z<–1.8; orange) phosphopeptides. (C) Principal component analysis of phosphopeptides identified as significantly changing. Symbols follow the schematic in A. (D) Gaussian mixture- model- based clustering of phosphopeptides changing during zaprinast treatment. Solid lines show the median relative abundance of each cluster. Opaque lines show the individual phosphopeptides belonging to each cluster. (E) GO terms enriched among phosphopeptides changing with zaprinast treatment, grouped by cluster. Gene ratio is the proportion of proteins with the indicated GO term divided by the total number of proteins belonging to each cluster. Significance was determined with a hypergeometric test; only GO terms with p<0.05 are shown. Redundant GO terms were removed. (F) Examples of phosphopeptides belonging to each cluster. (G) The number of clusters each phosphoprotein belongs to plotted against the number of changing phosphopeptides belonging to each protein. Gene names or IDs indicate proteins discussed in the text. The online version of this article includes the following figure supplement(s) for figure 1: Figure supplement 1. Metrics describing the zaprinast-d ependent phosphoproteome. set includes two putative Ca2+-a ctivated K+ channels and the sodium- hydrogen exchangers NHE1 and NHE3 (Table 1). ATPase2 and the guanylyl cyclase in this cluster both have putative phospholipid- transporting ATPase domains. PI4,5K, PI3,4K, and TgCDPK2A, noted in cluster 1, have additional phosphoregulatory sites belonging to cluster 2. Clusters 2 and 3 are enriched for phosphoproteins with metal transporter activity, including ATPases such as TgA1, the putative copper transporter CuTP, a putative P5B- ATPase, a putative E1- E2 ATPase, and TgATP4. As observed in other T. gondii Herneisen et al. eLife 2022;11:e80336. DOI: https://doi.org/10.7554/eLife.80336 4 of 41 log2 change rel. DMSO t cluster log2 change rel. DMSO t summed abundance zap peptides per protein PC2 (7.3%) Research article Cell Biology | Microbiology and Infectious Disease Table 1. Gene IDs of proteins discussed in the text. NS, not significant (using thresholds defined in the text). ND, not detected. TR, temperature range. CR, concentration range. Phospho Thermal Gene ID Description in text Reference cluster profiling TGGT1_202540 PDE1 Jia et al., 2017 1, 3 NS TGGT1_293000 PDE2 Jia et al., 2017; Moss et al., 2022; Vo et al., 2020 1, 2, 4 NS TGGT1_245730 PI4,5K Garcia et al., 2017 1, 2 NS TGGT1_276170 PI3,4K Garcia et al., 2017 1, 2 NS TGGT1_248830 PI-P LC Bullen et al., 2016; Fang et al., 2006 1, 3 NS Phosphatidylinositol- 3,4,5- triphosphate TGGT1_288800 5- phosphatase 1 NS TGGT1_254390 Putative Sec14 1 NS TGGT1_206590 CDPK2A Billker et al., 2009 1, 2, 3 CR TGGT1_228750 CDPK7 Bansal et al., 2021 1, 3 TR TGGT1_267100 PPM2B Yang et al., 2019; Yang and Arrizabalaga, 2017 1 NS TGGT1_238995 Ca2+- activated K+ channel 2 ND TGGT1_273380 Ca2+- activated K+ channel 2, 3 ND TGGT1_259200B Na+/H+ exchanger Arrizabalaga et al., 2004 2, 3 NS TGGT1_305180 Na+/H+ exchanger Francia et al., 2011 2, 3, 4 NS TGGT1_245510 ATPase2 Chen et al., 2021 2 ND TGGT1_254370 Guanylyl cyclase Bisio et al., 2019; Brown and Sibley, 2018 1, 2, 3, 4 NS TGGT1_312100 Calcium ATPase TgA1 Luo et al., 2005; Luo et al., 2001 3 ND TGGT1_201150 Copper transporter CuTP Kenthirapalan et al., 2014 3 NS TGGT1_318460 Putative P5B- ATPase Møller et al., 2008 2 NS TGGT1_289070 Putative E1- E2 ATPase 3 NS TGGT1_278660 TgATP4 Lehane et al., 2019 3, 4 NS TGGT1_226020 MFS transporter 3 NS TGGT1_230570 MFS transporter 3 NS TGGT1_253700 MFS transporter 3, 4 NS TGGT1_257530 Tyrosine transporter ApiAT5- 3 Parker et al., 2019; Wallbank et al., 2019 3, 4 NS TGGT1_292110 Formate transporter TgFNT2 Erler et al., 2018; Zeng et al., 2021 3 NS TGGT1_270865 Adenylyl cyclase Brown and Sibley, 2018; Jia et al., 2017 2, 3 NS TGGT1_238390 Unique guanylyl cyclase organizer UGO Bisio et al., 2019 2, 3 NS ARO-i nteracting protein (adenylyl cyclase TGGT1_309190 organizer) AIP Mueller et al., 2016; Mueller et al., 2013 1, 3 NS TGGT1_273560 Divergent kinesin KinesinB Leung et al., 2017 2, 3 NS TGGT1_201230 Divergent kinesin Wickstead et al., 2010 3 ND TGGT1_247600 Dynein light chain 3 NS TGGT1_255190 MyoC Frénal et al., 2017 3 NS TGGT1_278870 MyoF Heaslip et al., 2016; Jacot et al., 2013 1, 2, 3 NS TGGT1_257470 MyoJ Frénal et al., 2014 3 CR TGGT1_213325 Uncharacterized TBC domain protein 3 ND Table 1 continued on next page Herneisen et al. eLife 2022;11:e80336. DOI: https://doi.org/10.7554/eLife.80336 5 of 41 Research article Cell Biology | Microbiology and Infectious Disease Table 1 continued Phospho Thermal Gene ID Description in text Reference cluster profiling TGGT1_221710 Uncharacterized TBC domain protein 3 NS TGGT1_237280 Uncharacterized TBC domain protein 1, 3, 4 NS TGGT1_274130 Uncharacterized TBC domain protein 2, 3, 4 NS TGGT1_289820 Uncharacterized TBC domain protein 3 NS TGGT1_206690 GAPM2B Harding et al., 2019 3 NS TGGT1_271970 GAPM3 Harding et al., 2019 3 NS TGGT1_233010 ERK7 O’Shaughnessy et al., 2020 3, 4 CR TGGT1_234970 Tyrosine kinase- like protein TgTLK2 Varberg et al., 2018 4 NS Putative K+ voltage- gated channel TGGT1_202900 complex subunit 4 NS TGGT1_228200 Vacuolar (H+)- ATPase G subunit 1, 4 NS TGGT1_233130 Putative nucleoside transporter 4 NS TGGT1_258700 MFS family transporter 4 ND TGGT1_269260 SCE1 McCoy et al., 2017 1, 2, 3 NS TGGT1_295850 AAP2 Engelberg et al., 2020 1, 2, 4 CR TGGT1_319900 AAP5 Engelberg et al., 2020 1, 2, 3, 4 NS TGGT1_227000 Apical polar ring protein Koreny et al., 2021 1, 3, 4 ND TGGT1_244470 RNG2 Katris et al., 2014 1, 2, 3 ND TGGT1_292950 Apical cap protein AC10 Back et al., 2020; Tosetti et al., 2020 1, 2, 3 NS TGGT1_240380 Conoid gliding protein CGP Li et al., 2022 1, 3, 4 NS TGGT1_216620 Ca2+ influx channel with EF hands Chang et al., 2019 1, 3, 4 NS TGGT1_246930 Calmodulin- like protein CAM1 Long et al., 2017b ND TR TGGT1_262010 ACalmodulin- like protein CAM2 Long et al., 2017b ND TR, CR TGGT1_216080 apical lysine methyltransferase (AKMT) Heaslip et al., 2011 NS TR TGGT1_270690 DrpC Heredero-B ermejo et al., 2019; Melatti et al., 2019 NS TR TGGT1_201880 TgQCR9 ND TR Hayward et al., 2021; Huet et al., 2018; Mühleip TGGT1_204400 ATP synthase subunit alpha et al., 2021; Salunke et al., 2018; Seidi et al., 2018 NS TR Huet et al., 2018; Mühleip et al., 2021; Salunke TGGT1_231910 ATP synthase subunit gamma et al., 2018; Seidi et al., 2018 ND TR Huet et al., 2018; Mühleip et al., 2021; Salunke TGGT1_208440 ATP synthase subunit 8/ASAP- 15 et al., 2018; Seidi et al., 2018 ND TR, CR Huet et al., 2018; Mühleip et al., 2021; Salunke TGGT1_215610 ATP synthase subunit f/ICAP11/ASAP- 10 et al., 2018; Seidi et al., 2018 ND TR ATP synthase- associated protein Huet et al., 2018; Mühleip et al., 2021; Salunke TGGT1_263080 ASAP- 18/ATPTG14 et al., 2018; Seidi et al., 2018 ND TR, CR TGGT1_246540 ATP synthase- associated protein ATPTG1 Mühleip et al., 2021 NS TR TGGT1_249240 CaM Paul et al., 2015 ND CR TGGT1_213800 CnB Paul et al., 2015 NS CR TGGT1_227800 Eps15 Birnbaum et al., 2020; Chern et al., 2021 NS CR TGGT1_269442 ELC1 Nebl et al., 2011 ND CR Table 1 continued on next page Herneisen et al. eLife 2022;11:e80336. DOI: https://doi.org/10.7554/eLife.80336 6 of 41 Research article Cell Biology | Microbiology and Infectious Disease Table 1 continued Phospho Thermal Gene ID Description in text Reference cluster profiling TGGT1_297470 MLC1 Gaskins et al., 2004 ND CR TGGT1_297470 MLC5 Graindorge et al., 2016 ND CR TGGT1_315780 MLC7 Graindorge et al., 2016 ND CR TGGT1_226030 PKA-C 1 Jia et al., 2017; Uboldi et al., 2018 NS CR TGGT1_242070 PKA- R Jia et al., 2017; Uboldi et al., 2018 NS CR TGGT1_210830 Putative RIO1 kinase NS CR TGGT1_310700 PP1 Paul et al., 2020; Zeeshan et al., 2021 ND CR TGGT1_207910 Calcium-h ydrogen exchanger TgCAX Guttery et al., 2013 NS CR TGGT1_311080 Apicoplast two-p ore channel TgTPC Li et al., 2021 NS CR TGGT1_204050 Subtilisin 1 SUB1 NS CR TGGT1_206490 Metacaspase 1 with a C2 domain Li et al., 2015 ND CR TGGT1_310810 Ca2+- activated apyrase ND CR TGGT1_321650 RON13 Lentini et al., 2021 ND CR Uncharacterized metal- binding protein TGGT1_286710 with zinc fingers ND CR Uncharacterized metal- binding protein TGGT1_309290 with HD domain ND CR TGGT1_225690 Apical cap protein AC7 Chen et al., 2015 2, 3 NS TGGT1_234250 CHP interacting protein CIP1 Long et al., 2017a NS ND Ca2+- stimulated phosphoproteomes, at later time points small- molecule transporters are phosphory- lated, including MFS transporters, ApiAT5- 3, and TgFNT2. The guanylyl and adenylyl cyclases are also extensively modified, as are the cyclase organizers UGO and AIP. Cluster 3, which represents a latter wave of phosphoregulation, is uniquely enriched in proteins functioning in subcellular remodeling, vesicle trafficking, and glideosome activity (Figure 1E). This set includes two divergent kinesins (KinesinB and TGGT1_201230), a dynein light chain, and the myosin motors MyoC, MyoF, and MyoJ (Table  1). Several uncharacterized GTPase regulators are phos- phorylated later in the zaprinast response, including putative ARF1 activators (TGGT1_225310 and TGGT1_266830) as well as five uncharacterized TBC domain proteins. The glideosome- associated membrane proteins GAPM2B and GAPM3, which link the IMC, alveolin network, and microtubules, are also dynamically phosphorylated. Cluster 4, the only cluster characterized by decreasing phosphorylation, was functionally enriched in phosphoproteins involved in transporter and MAP kinase activity. ERK7 was dephosphorylated within 30 s of zaprinast treatment, whereas TgTLK2 was dephosphorylated only at the final 60 s time point. ERK7 regulates conoid and cytoskeletal stability during cell division, with secondary func- tions in parasite egress, motility, and invasion (Back et al., 2020; Dos Santos Pacheco et al., 2021; O’Shaughnessy et al., 2020). TLK2 regulates parasite replication (Smith et al., 2022; Varberg et al., 2018). Several phosphoproteins functioning in ion transport belonged to cluster 4, including a puta- tive K+ voltage- gated channel complex subunit, a putative Vacuolar (H+)-A TPase G subunit, TgATP4, and the guanylyl cyclase (Table 1), which has a P-t ype ATPase domain with unknown ion specificity. Small- molecule transporters included a putative nucleoside transporter, ApiAT5- 3, and a MFS family transporter. Cluster 4 was the only class of peptides that was not functionally enriched in Ca2+ binding proteins (Figure 1E). We identified phosphoproteins with peptides belonging to several clusters (Figure  1G). Such proteins may have multiple phosphosites regulated with different kinetics by the same enzyme; or by Herneisen et al. eLife 2022;11:e80336. DOI: https://doi.org/10.7554/eLife.80336 7 of 41 Research article Cell Biology | Microbiology and Infectious Disease different enzymes that alight upon the target at varying spatiotemporal scales. For example, SCE1 and TGGT1_309910 have phosphopeptides belonging to all three increasing clusters, likely resulting from phosphorylation by TgCDPK3 and PKG, respectively. Indeed, SCE1 was implied to be a TgCDPK3 target through a genetic suppressor screen (McCoy et al., 2017). TGGT1_309910 is the ortholog of Plasmodium falciparum Pf ICM1, a PKG substrate identified through proteomic interaction studies (Balestra et al., 2021), although it remains uncharacterized in T. gondii. By contrast, several proteins belong to both increasing and decreasing clusters, likely targets of both kinases and phosphatases. We consider these proteins candidate signaling platforms. Several such proteins localize to discrete domains of the apical complex, including the guanylyl cyclase; apical annuli proteins AAP2 and AAP5; apical polar ring proteins TGGT1_227000 and RNG2; and the apical cap protein AC10 and a recently identified conoid gliding protein CGP (Table 1). Along with the guanylyl cyclase, the cAMP- specific phosphodiesterase PDE2 and a Ca2+ influx channel with EF hands exhibit peptides belonging to both increasing and decreasing clusters. Our phosphoproteome thus identifies candidate mediators of the feedback between cyclic- nucleotide and Ca2+ signaling. The waves of phosphoregulation largely corroborate the sequence of physiological events observed following zaprinast treatment. The first cluster likely includes the most proximal targets of PKG. Previous PKG-d ependent phosphoproteomes from the related apicomplexan P. falciparum simi- larly implicate PI-P LC as a substrate of the kinase (Alam et al., 2015; Balestra et al., 2021; Brochet et al., 2014). The product of PI-P LC activity, IP3, stimulates Ca2+ store release in parasites (Garcia et al., 2017). In turn, phosphoproteins in clusters 2, 3 and 4 likely contain the targets of Ca2+-r egulated kinases and phosphatases. Thermal profiling identifies Ca2+-dependent shifts in T. gondii protein stability The phosphoproteome identified numerous dynamic changes in response to Ca2+ release. However, this information alone is not sufficient to infer the enzymes responsible for phosphorylation, or which events may be functionally relevant. Thermal proteome profiling (TPP) has been used to identify small molecule–target interactions in living cells and cell extracts (Dai et al., 2019; Mateus et al., 2020; Savitski et al., 2014). TPP operates on the premise that ligand binding induces a thermal stability shift, stabilizing, or destabilizing proteins that change conformationally in response to the ligand, and such changes in stability can be quantified by MS (Dziekan et al., 2020; Reinhard et al., 2015). Cells are incubated with different concentrations of ligand and heated, causing thermal denaturation of proteins. The soluble protein is extracted and quantified with multiplexed, quantitative methods, giving rise to thousands of thermal denaturation profiles. Proteins engaging the ligand are identified by their concentration-d ependent thermal shift. We previously used this method to identify the target of the antiparasitic compound ENH1 (Herneisen et al., 2020; Herneisen and Lourido, 2021) and to measure changes to the T. gondii proteome when depleting the mitochondrial protein DegP2 (Harding et al., 2020). In a conceptual leap, we reasoned that TPP could also detect Ca2+- responsive proteins if parasite extracts were exposed to different concentrations of Ca2+, allowing us to system- atically identify the protein components of signaling pathways on the basis of biochemical interactions with Ca2+ and its effectors. Intracellular free Ca2+ levels span three orders of magnitude, from low nanomolar in the cytoplasm to high micromolar in organelles like the ER (Lourido and Moreno, 2015). To measure protein thermal stability at precisely defined [Ca2+]free, we combined crude parasite extracts with calibrated Ca2+ buffers in a solution mimicking the ionic composition of the cytoplasm. As a proof of principle, we measured the thermal stability of the Ca2+-d ependent protein kinase TgCDPK1, as the conformational changes of this enzyme have been structurally characterized (Ingram et al., 2015; Wernimont et al., 2010). Parasite lysates were adjusted to 10 different concentrations of free Ca2+ and heated to 58 °C, based on prior estimates of the melting temperature of TgCDPK1 (Herneisen et al., 2020). As measured by immunoblot band intensity, TgCDPK1 was strongly stabilized by Ca2+ (Figure 2A), suggesting that our experimental system is sensitive to Ca2+-d ependent stability changes. The calculated EC50 using this approach was in the low µM range, consistent with studies using recombinant enzymes (Ingram et al., 2015; Wernimont et al., 2010). The effect of Ca2+ on the global thermostability of the proteome has not been assessed. There- fore, we first generated thermal profiles of the T. gondii proteome without or with 10 µM Ca2+, which Herneisen et al. eLife 2022;11:e80336. DOI: https://doi.org/10.7554/eLife.80336 8 of 41 Research article Cell Biology | Microbiology and Infectious Disease A C EGTA 10 µM Ca2+ ECDPK1 ICAP110.4 TGGT1_214340 EF hand protein TUB1 apicomplexan proteinCDPK7 1.00 EC50 = 2.4 µM AKMT 0.75 1.4 DrpC 0.3 CAM2 1.2 CAM1 median + 2MAD0.50 TgQCR9 0.25 1.0 0.8 0.2 0.00 median CDPK1 0.6 10-1 102 103 0.4 0.1 [Ca2+] (µM) 0.2 0.0 B 0.0 parasite lysates 0 500 1000 1500 2000 2500 °C protein rank [Ca2+] F experiment 1 2 1.5 EGTA 10 µM Ca2+ heat 1.0isolate soluble 0.5 fraction °C DNA Pol β CDPK7 CAM2 (spin/filter) D 1.25 TMT10plex 1.00 ED MS 0.75 EGTA 0.50 °C AUC 10 µM Ca2+ 0.25 255900 206320 AKMT 40 50 60 40 50 60 40 50 60 EC50 40 45 50 55 60 65 70(HFX or FAIMS) [Ca2+] temperature (°C)Tm (°C) Figure 2. Thermal profiling identifies proteins that change stability in response to Ca2+. (A) Thermal shift assays can detect Ca2+-d ependent stability of CDPK1 in extracts. Parasite lysates were combined with 10 concentrations of Ca2+ spanning the nanomolar to micromolar range. After denaturation at 58 °C, the soluble fraction was separated by SDS- PAGE and probed for CDPK1. Band intensity was normalized to the no- Ca2+ control and scaled. Points in shades of gray represent two different replicates. A dose- response curve was calculated for the mean abundances. (B) Schematic of the thermal profiling workflow. In the temperature- range experiment, parasite lysates were combined with EGTA or 10 µM [Ca2+]free and heated at 10 temperatures spanning 37–67°C. In the concentration-r ange experiment, parasite lysates were combined with 10 different [Ca2+]free (nM–mM range) and heated at 50, 54, or 58 °C. Temperature- range shifts were quantified by the Euclidean distance (ED) score, a weighted ratio of thermal stability differences between treatments and replicates. Concentration- range shifts were summarized by pEC50, area under the curve (AUC), and goodness of fit (R2). (C) Heat map of protein thermal stability relative to the lowest temperature (37 °C) in 0 or 10 µM Ca2+. The mean relative abundance at each temperature was calculated for 2381 proteins. Proteins are plotted in the same order in both treatments. (D) Raincloud plots summarizing the distribution of Tm in lysates with EGTA (gray) or 10 µM [Ca2+]free (blue). The average melting temperatures of proteins identified in two replicates were plotted. (E) Proteins rank-o rdered by euclidean distance score quantifying the Ca2+- dependent shift in thermal stability. Solid and dotted lines represent the median ED score and two modified Z scores above the median, respectively. Highlighted proteins have EF hand domains (blue) or are conserved in apicomplexans (pink). (F) Thermal profiles of individual proteins: DNA polymerase β (TGGT1_233820); the EF hand domain-c ontaining proteins CDPK7 (TGGT1_228750) and the calmodulin- like protein CAM2 (TGGT1_262010); potential Ca2+- leak channels TGGT1_255900 and TGGT1_206320; and AKMT (TGGT1_216080). The online version of this article includes the following source data and figure supplement(s) for figure 2: Source data 1. This file contains the source data that was quantified to make the graph presented in Figure 2. Source data 2. This file contains the source data that was quantified to make the graph presented in Figure 2. Source data 3. This file contains the annotated source data that was quantified to make the graph presented in Figure 2. Figure supplement 1. Extended data for thermal profiling experiments. is representative of the resting and stimulated Ca2+ concentrations of cell cytoplasm (Lourido and Moreno, 2015). A thermal challenge between 37 and 67°C induced denatured aggregates, which were separated from stable proteins by ultracentrifugation. The soluble fraction was digested and labeled with isobaric mass tags, pooled, fractionated, and analyzed with an orbitrap mass spectrom- eter (Figure 2B). We detected 3,754 proteins, of which 2,381 yielded thermal profiles for both condi- tions in both replicates (Figure  2C and Figure  2—figure supplement 1); the remaining proteins were not detected in all experiments. The median melting temperatures—the temperatures at which proteins are 50% denatured—were 54.4 and 54.7°C in the lysates without and with Ca2+, respectively. The distribution of melting temperatures was largely overlapping in the two conditions (Figure 2D), Herneisen et al. eLife 2022;11:e80336. DOI: https://doi.org/10.7554/eLife.80336 9 of 41 relative abundance abundance proteins relative abundance euclidian distance score relative abundance Research article Cell Biology | Microbiology and Infectious Disease suggesting that Ca2+- dependent changes in protein stability were restricted to a subset of proteins. We additionally calculated an area under the curve (AUC) metric by numerical integration using the trapezoidal rule (Herneisen and Lourido, 2021) to compare the stabilities of proteins with atypical melting behavior (Figure 2—figure supplement 1), such as components of the tubulin cytoskeleton or parasite conoid. To discover proteins with Ca2+- dependent stability shifts in the initial temperature range experiment, we rank-o rdered proteins by Euclidean distance (ED) scores (Dziekan et al., 2020) quantifying the shift in thermal profiles with and without Ca2+ (Figure 2B and E). The majority of proteins, such as DNA polymerase β (TGGT1_233820), exhibited similar melting behavior in both conditions (Figure 2F). Our analysis identified as Ca2+- responsive parasite-s pecific proteins with EF hands, including TgCDPK7 and the calmodulin- like proteins CAM1 and CAM2 (Table 1, Figure 2F and Figure 2—figure supple- ment 1). The ED metric identified the stability changes of both CAM proteins despite a lack of typical melting behavior, supporting the use of this statistic. Membrane proteins, including potential Ca2+- leak channels localizing to the ER (Barylyuk et al., 2020; TGGT1_255900 and 206320), also exhibited thermal shifts (Figure 2F). Several proteins specific to the apicomplexan parasite phylum exhibited Ca2+- regulation (Barylyuk et al., 2020; Sidik et al., 2016a), including an apical lysine methyltrans- ferase (AKMT) that relocalizes during Ca2+-s timulated egress Heaslip et al., 2011; DrpC, which regu- lates the stability of parasite organelles; a hypothetical protein (TGGT1_214340) with no annotated domains (Figure 2—figure supplement 1); and enzyme subunits involved in cellular metabolism. The data also inform hypotheses about Ca2+ homeostasis and energetics in the parasite mito- chondrion and apicoplast. A divergent subunit of the mitochondrial complex III, TgQCR9, as well as components of the ATP synthase complex, were destabilized by Ca2+ (Figure 2—figure supplement 1). These include the ATP synthase subunits α, γ, 8/ASAP- 15, and f/ICAP11/ASAP-1 0; and the ATP synthase- associated proteins ASAP-1 8/ATPTG14 and ATPTG1 (Table 1). TGGT1_209950, a conserved alveolate thioredoxin- like protein suggested to localize to the apicoplast by spatial proteomics (Bary- lyuk et al., 2020), was similarly destabilized by Ca2+ (Figure 2—figure supplement 1). The apicom- plexan mitochondrion has been reported to sense cytosolic Ca2+ fluctuations (Gazarini and Garcia, 2004), and the apicoplast may take up Ca2+ through direct contact sites with the ER (Li et al., 2021). Furthermore, redox signals were recently reported to induce parasite Ca2+ signaling and motility (Alves et al., 2021; Stommel et al., 1997). Our resource identifies several proteins that may couple ion homeostasis and cellular metabolism. The recent structure of the apicomplexan ATP synthase suggests that all of the Ca2+- responsive elements identified here, including subunits restricted to apicomplexan parasites, protrude into the mitochondrial matrix (Mühleip et  al., 2021), which has elevated Ca2+ (Giorgi et al., 2018). In the case of TGGT1_209950, plants and algae have Ca2+- sensing thioredoxins (calredoxins) in chloroplasts (Hochmal et al., 2016) with Ca2+-r egulated activity mediated by EF hands (Charoenwattanasatien et al., 2020). TGGT1_209950 lacks EF- hands; however, struc- tural modeling based on primary sequence and predicted secondary features (Kelley et al., 2015; Meier and Söding, 2015) suggested similarity to calsequestrin, which lacks a structured Ca2+- binding motif. Our approach has thus identified candidate Ca2+-r esponsive proteins in metabolic organelles that would have been missed by bioinformatic searches; although, additional experiments are needed to determine the interactions each of these candidates has with Ca2+. Determining the specificity and sensitivity of Ca2+-responsive proteins The temperature range experiment has the advantage of generating complete thermal stability profiles, but does not inform the relative affinities of the Ca2+- dependent stability change. To address this gap, we examined protein stability across 10 Ca2+ concentrations (EGTA and 10 nM to 1 mM). Based on the temperature- range experiments, we selected thermal challenge temperatures of 50, 54, and 58 °C to target protein with a range of thermal stabilities in these dose- response experiments (Figure  2—figure supplement 1). We hypothesized that Ca2+- responsive proteins would exhibit sigmoidal dose- response trends in thermal stability, similarly to TgCDPK1 (Figure 2A). We performed four independent concentration- range experiments on two mass spectrometers using different sepa- ration methods (ultracentrifugation or filtration) to capture different types of aggregates. The concentration range thermal- profiling experiments provide insight into the magnitude of the Ca2+- dependent change (AUC parameter) and its dose-d ependency (EC50) (Figure  2B). Clear responses were observed for several known Ca2+-b inding enzymes with EF hands, including calmodulin, Herneisen et al. eLife 2022;11:e80336. DOI: https://doi.org/10.7554/eLife.80336 10 of 41 Research article Cell Biology | Microbiology and Infectious Disease calcineurin, and several parasite Ca2+-d ependent protein kinases (CDPKs; Figure 3A). These parame- ters can be used to classify proteins as stabilized (AUC >1; e.g., TgCDPK1, TgCDPK2A, and TgCDPK3) or destabilized (AUC <1; e.g., calmodulin and calcineurin B) by Ca2+ (Figure 3A). Furthermore, EC50 measurements may inform specific hypotheses about a protein’s involvement in signaling or homeo- stasis, based on the Ca2+ concentration at which the protein is predicted to change. Our results generated a de novo catalog of proteins with Ca2+- dependent stability. Of the 4,403 proteins identified with sufficient quantification values for curve fitting, 228 were classified as Ca2+- responsive based on exhibiting a dose-r esponse trend (R2) greater than 0.8 with a stability change (AUC) of two modified Z scores from the median (Figure 3B, Figure 3—figure supplement 1, and Supplementary file 2). Functional enrichment showed that such Ca2+- responsive proteins were signifi- cantly enriched in Ca2+- binding functions, protein kinase activity, and metal affinity (Figure 3C). We examined the 40 EF hand domain–containing proteins detected in our mass spectrometry experi- ments (Figure 3D), 25 of which exhibited dose-r esponsive behavior. Signaling and trafficking proteins were tuned to respond at lower Ca2+ concentrations: TgCDPK2A and CaM responded robustly with nM EC50; TgCDPK1, TgCDPK3, CnB, and Eps15 exhibited low µM EC50. In the cases of TgCDPK1, CaM, and CnB, these values are consistent with the Ca2+ activation constants determined from studies with recombinant proteins (Feng and Stemmer, 2001; Feng and Stemmer, 1999; Wernimont et al., 2010). By contrast, myosin-a ssociated subunits responded at higher levels: Myosin light chains (ELC1, MLC1, MLC5, and MLC7) and the calmodulin- like proteins CAM1 and CAM2 responded above 10 µM Ca2+ (Table 1). Two additional proteins (TGGT1_255660 and TGGT1_259710) consistently responded with µM EC50 values but have not been functionally characterized (Figure  3—figure supplement 1). The remaining EF hand domain–containing proteins may bind Ca2+ outside of the experimental concentration range; or may not undergo detectable stability changes upon Ca2+ binding (Mateus et al., 2020). Several signaling enzymes lacking canonical Ca2+- binding domains exhibited a dose-r esponse trend (R2 >0.8) in our thermal profiling experiments (Figure 3B and E). The PKA regulatory and cata- lytic subunits (TGGT1_242070 and TGGT1_226030) responded robustly to Ca2+ with 10–100 µM EC50 values in several experiments. A putative RIO1 kinase (TGGT1_210830) also responded consistently to Ca2+, although the magnitude of the change was small (Figure  3—figure supplement 1). The atypical kinase ERK7 was destabilized at high Ca2+ concentrations (Figure 3—figure supplement 1), similarly to other Ca2+- responsive proteins detected at the parasite apex. Rhoptry and dense granule kinases had EC50 values in between 100 µM and 1 mM, consistent with the high concentration of Ca2+ in the secretory pathway (Figure 3E). We also searched for Ca2+- responsive phosphatases. The known Ca2+- regulated phosphatase subunit, CnB, was destabilized (Figure 3F). Protein Phosphatase 1 (PP1) responded to Ca2+ consistently (Figure 3F and 5A), although the catalytic subunit is not an intrinsic Ca2+ sensor. Our resource thus places proteins without previously characterized Ca2+ responses within the broader Ca2+ signaling network. Our search and annotation of Ca2+- responsive proteins also identified candidates that link Ca2+ and ion homeostasis. This list includes several transporters and channels (Table 1): TgCAX, a protein with structural homology to LDL receptors (TGGT1_245610), two proteins with structural homology to the mitochondrial Ca2+ uniporter (TGGT1_257040 and TGGT1_211710), and an apicoplast two-p ore channel. Several Ca2+- responsive hydrolases were also found, including SUB1, a metacaspase with a C2 domain, and a Ca2+- activated apyrase. These proteins may function during the intracellular, repli- cative phase of the lytic cycle, for which the function of Ca2+ is relatively unexplored. Several divalent cations regulate protein structure and activity in addition to Ca2+, including Mg2+ and Zn2+. Such metal- binding proteins might appear Ca2+- responsive in our thermal profiling approach. We buffered our parasite extracts with excess Mg2+ (1 mM) to mitigate non- specific changes caused by the concentration of divalent cations. Nevertheless, to determine whether our approach identified additional metal- binding proteins, we compared the EC50 values of candidates that bind different diva- lent cations, based on the presence of Interpro domains and through manual annotation (Figure 3— figure supplement 1). The EC values of putative Mg2+50 - binding proteins were mostly in the mM range, which may result from nonspecific displacement of Mg2+ by Ca2+. However, a subset of Mg2+- binding proteins were stabilized by micromolar Ca2+–for instance, glucosephosphate-m utase (GPM1) family proteins TGGT1_285980  A and B. The GPM1 ortholog in Paramecium, parafusin (PFUS), is involved in Ca2+- regulated exocytosis (Subramanian and Satir, 1992) A PFUS homolog in T. gondii Herneisen et al. eLife 2022;11:e80336. DOI: https://doi.org/10.7554/eLife.80336 11 of 41 Research article Cell Biology | Microbiology and Infectious Disease A stabilized destabilized 301440 CDPK1 206590 CDPK2A 305860 CDPK3 213800 CnB 249240 CaM 1.6 1.75EC50 = 28 µM EC50 = 70 nM EC50 = 12 µM 1.25 EC = 1.5 µM 1.0050 ECexperiment 1 2 50 = 28 nM 2.0 50°C 1.4 1.50 1.00 0.75 54°C 58°C 1.6 1.25 0.75 1.2 0.50 1.2 0.501.00 0.25 1.0 0.25 10-1 102 103 10-1 102 103 10-1 102 103 10-1 102 103 10-1 102 103 [Ca2+] (µM) B Ca2+ affinity C 3 centrifugation 226030 PKA C1 3 filtration GO:00055092+ 227800 Eps15 Ca bindingnM 2 242070 PKA R GO:0004674286710 zinc finger 2 protein S/T kinase 309290 hyp. 1 310700 PP1 GO:0004672 321650 RON13 1 protein kinase GO:0016773 OH-phosphotransferase µM 0 0 GO:0016301 kinase –log10(p) set size−1 −1 GO:0016772 phosphotransfer 5 100200 −2 −2 GO:0046872 4 metal ion binding 3 300 mM GO:0043169 2 −3 −3 cation binding 0 1 2 3 0 1 2 3 0 5 10 15 20 25 AUC AUC fold enrichment D cutoff AUC R2 experiment 0 0 count 1 1 54°C 2 2 58°C 3 3 50°C 54°C 4 4 58°C 5 5 –log10(EC50) µM3 2 E F 1 0 −1 −2 experiment count −3 54°C 58°C 50°C 54°C 58°C Figure 3. Thermal profiling identifies anticipated and unexplored Ca2+-r esponsive proteins. (A) Mass spectrometry- derived thermal profiles of EF hand- containing proteins stabilized or destabilized by Ca2+. Relative abundance is calculated relative to the protein abundance at 0 µM Ca2+. EC50 is the median of the EC50 values of the curves displayed on the plots. (B) The magnitude of Ca2+- dependent stabilization (AUC) plotted against the sensitivity (pEC50) for protein abundances exhibiting a dose- response trend with an R2 >0.8. Point size is scaled to R2. Summary parameters for the different separation methods (ultracentrifugation or filtration) are plotted separately. Colors identify candidates with Ca2+- responsive behavior validated in Figure 4. (C) Gene ontology (GO) terms enriched among candidate Ca2+- responsive proteins (AUC greater than two modified Z scores and R2 dose- response >0.8). Fold enrichment is the frequency of Ca2+- responsive proteins in the set relative to the frequency of the GO term in the population of detected proteins. Significance was determined with a hypergeometric test; only GO terms with p<0.05 are shown. (D–F) EF hand domain proteins (D), protein kinases (E), and protein phosphatases (F) detected in the thermal profiling mass spectrometry datasets. The top rows indicate if a protein passed the AUC cutoff (orange) or R2 cutoff (blue) for dose- response behavior. The opacity of the band represents the number of experiments in which the protein exhibited the behavior (out of five). The five rows below summarize the pEC50 of each experiment in which the protein exhibited a dose- response trend with R2 >0.8. Kinases are loosely grouped as CDPK’s (included as a reference), non- rhoptry kinases, and secretory pathway kinases. The online version of this article includes the following figure supplement(s) for figure 3: Figure supplement 1. Extended analysis of thermal profiling experiments. Herneisen et al. eLife 2022;11:e80336. DOI: https://doi.org/10.7554/eLife.80336 12 of 41 Exp2 Exp1 Exp2 Exp1 –log10(EC50) µM CDPK1 TGGT1_301440 TGGT1_200400 relative abundance CDPK2A TGGT1_206590 CDPK2A TGGT1_206590 CDPK3 TGGT1_305860 CnB TGGT1_213800 TGGT1_216620 CDPK9 TGGT1_217600 PKA C1 TGGT1_226030 CDPK6 TGGT1_218720 PKA R TGGT1_242070 MSC family TGGT1_219650 RIO1 TGGT1_210830 TGGT1_220140 ERK7 TGGT1_233010 CDPK5 TGGT1_224950 TLK3 TGGT1_253860 TGGT1_225490 TGGT1_260030 CAM3 TGGT1_226040 CDK TGGT1_304970 Eps15 TGGT1_227800 TGGT1_240640 CDPK7 TGGT1_228750 TGGT1_245580 TGGT1_229480 CMCG TGGT1_263070 RON11 TGGT1_230350 GSK TGGT1_265330 PI4,5K TGGT1_230490 NIMA TGGT1_292140 TGGT1_231210 PKG TGGT1_311360 CAM1 TGGT1_246930 PRP4 TGGT1_313180 centrin 1 TGGT1_247230 PI-PLC TGGT1_248830 ROP5 TGGT1_308090 calmodulin TGGT1_249240 ROP8 TGGT1_215775 centrin 2 TGGT1_250340 ROP11 TGGT1_227810 TGGT1_251850 ROP16 TGGT1_262730 TGGT1_255660 ROP17 TGGT1_258580 MLC1 TGGT1_257680 ROP18 TGGT1_205250 TGGT1_259710 ROP20 TGGT1_258230 CEN3 TGGT1_260670 ROP25 TGGT1_202780 CAM2 TGGT1_262010 ROP26 TGGT1_211260 ELC1 TGGT1_269442 ROP28 TGGT1_258370 EFPP TGGT1_269460 ROP34 TGGT1_240090 TGGT1_280480 ROP35 TGGT1_304740 CDPK8 TGGT1_292055 ROP39 TGGT1_262050 CDPK4A TGGT1_295760 ROP40 TGGT1_291960 MLC2 TGGT1_297470 RON11 TGGT1_321650 CDPK1 TGGT1_301440 ELC2 TGGT1_305050 CDPK3 TGGT1_305860 MLC5 TGGT1_311260 CnB TGGT1_213800 MLC7 TGGT1_315780 PP1 TGGT1_310700 TGGT1_315860 Research article Cell Biology | Microbiology and Infectious Disease has been observed to co-l ocalize with micronemes in the apical third of the cell (Matthiesen et al., 2001). Therefore, the thermal stabilization of GPM1 could arise from Ca2+- dependent modifications to the enzyme (Subramanian et al., 1994; Subramanian and Satir, 1992). Validation of Ca2+-dependent thermal stability We selected five proteins exhibiting consistent Ca2+- responsive behavior by MS for validation: an EF hand domain–containing protein (Eps15), two kinases (RON13 and PKA-C 1), and two uncharacter- ized putative metal-b inding proteins (TGGT1_286710 and TGGT1_309290). The first three candidates were selected for potential involvement in dynamic Ca2+- regulated processes. Eps15 (TGGT1_227800) was recently shown to mediate endocytosis in P. falciparum (Birnbaum et  al., 2020) and local- ized to puncta bridging the inner membrane complex (IMC) and cytoskeleton in T. gondii (Chern et al., 2021). PKA-C 1 is thought to antagonize Ca2+ signaling in T. gondii (Jia et al., 2017; Uboldi et al., 2018). RON13 is a secretory-p athway kinase that phosphorylates substrates in apicomplexan invasion-a ssociated organelles called rhoptries (Lentini et al., 2021). The function of the remaining two proteins is unknown, although TGGT1_286710 contains a zinc finger domain and TGGT1_309290, annotated as a hypothetical protein, contains an HD/PDEase domain—both contribute to parasite fitness in cell culture (Sidik et al., 2016a). We appended epitope tags to the endogenous locus of each candidate. These Ca2+-r esponsive proteins localized to diverse structures (Figure 4B). PKA- C1 resided at the parasite pellicle, as previ- ously reported (Jia et al., 2017; Uboldi et al., 2018). Eps15 concentrated at puncta in the apical third of the parasite. TGGT1_286710 and TGGT1_309290 localized to the parasite nucleus, with TGGT1_286710 exhibiting additional staining in the residual body. Finally, RON13 appeared to stain the rhoptries (Lentini et al., 2021). To validate the Ca2+- dependent stability of individual candidates, we prepared parasite extracts as described earlier, but relied on an immunoblot readout instead of MS. In these five cases, stabilization of the candidates was confirmed in multiple biological repli- cates. Several controls (TUB1, MIC2, GAP45, and SAG1) were stable across all Ca2+ concentrations tested (Figure 4 and Figure 4—figure supplement 1). In the case of PKA- C1, Eps15, and both puta- tive metal- binding proteins, the immunoblot experiments revealed an even higher EC50 than was measured in the MS experiments. We conclude that the stability changes detected by our global proteomics methods are robust, although the precise features of Ca2+-d ependent stability may differ based on the method used to assess them. A PP1 holoenzyme performs Ca2+-responsive functions required for parasite spread Our orthogonal proteomics approaches map Ca2+-r esponsive phosphorylation and thermal stability. Unexpectedly, the catalytic subunit of Protein Phosphatase 1 (PP1, TGGT1_310700) exhibited consis- tent stabilization by Ca2+ in our mass spectrometry experiments (Figure  3F and Figure  5A). The contribution of phosphatases to Ca2+ signaling in apicomplexans is poorly understood (Yang and Arrizabalaga, 2017) and the Ca2+- responsive behavior of PP1 has not been reported in other eukary- otes. The function of the phosphatase has not been directly studied in T. gondii, although PP1 inhibi- tors have been shown to block invasion of host cells (Delorme et al., 2002), and recent experiments in P. falciparum suggest that PP1 is required for the merozoite egress- to-i nvasion transition (Paul et al., 2020). Intriguingly, PP1 was recently observed to relocalize to the apical complex in highly motile Plasmodium berghei ookinetes (Zeeshan et al., 2021), suggesting that the enzyme may serve apicomplexan- specific, Ca2+-r esponsive functions in remodeling the parasite phosphoproteome. To track localization during the T. gondii lytic cycle, we tagged the endogenous C terminus of PP1 with an mNeonGreen (mNG) fluorophore and Ty epitope. We confirmed the Ca2+- dependent stability of PP1 by immunoblotting (Figure 5B). Live microscopy revealed a diverse array of PP1 localizations in parasites (Figure 5C). In accordance with imaging in Plasmodium (Zeeshan et al., 2021), PP1- mNG was distributed diffusely in the cytoplasm, as well as in foci resembling the nucleus, centrosome, and in some parasites, the periphery. These diverse localizations may arise from the association of PP1 with distinct regulatory subunits forming different functional holoenzymes, as characterized in metazoans (Brautigan and Shenolikar, 2018). To determine whether PP1 exhibits Ca2+-d ependent relocaliza- tion, we treated parasites with zaprinast and A23187. The PP1 signal intensity increased at the apical cap and pellicle following zaprinast stimulation (Figure 5C and Video 1). However, Ca2+ ionophore Herneisen et al. eLife 2022;11:e80336. DOI: https://doi.org/10.7554/eLife.80336 13 of 41 Research article Cell Biology | Microbiology and Infectious Disease Figure 4. Validation of Ca2+-d ependent thermal stability. (A) Mass spectrometry- derived thermal profiles of the candidates, as in Figure 3A. (B) Immunofluorescence images of fixed intracellular parasites expressing the indicated proteins with C- terminal epitopes at endogenous loci. Hoechst and anti-C DPK1 were used as counterstains in the merged image. Green arrowheads highlight an example of TGGT1_286710 residual body staining. In the case of PKA C1/R, the stain of the R subunit is shown, as both subunits colocalize. (C) Immunoblot- derived thermal profiles of the candidates. Colors correspond to two independent replicates. Uncropped blots are shown in the Figure 4—figure supplement 1. The online version of this article includes the following source data and figure supplement(s) for figure 4: Source data 1. This file contains the source data that was quantified to make the graph presented in Figure 4. Source data 2. This file contains the source data that was quantified to make the graph presented in Figure 4. Source data 3. This file contains the source data that was quantified to make the graph presented in Figure 4. Source data 4. This file contains the source data that was quantified to make the graph presented in Figure 4. Source data 5. This file contains the source data that was quantified to make the graph presented in Figure 4. Source data 6. This file contains the source data that was quantified to make the graph presented in Figure 4. Source data 7. This file contains the source data that was quantified to make the graph presented in Figure 4. Source data 8. This file contains the source data that was quantified to make the graph presented in Figure 4. Source data 9. This file contains the source data that was quantified to make the graph presented in Figure 4. Source data 10. This file contains the source data that was quantified to make the graph presented in Figure 4. Source data 11. This file contains the annotated source data that was quantified to make the graph presented in Figure 4. Figure supplement 1. Uncropped immunoblots, as in Figure 4C. Herneisen et al. eLife 2022;11:e80336. DOI: https://doi.org/10.7554/eLife.80336 14 of 41 Research article Cell Biology | Microbiology and Infectious Disease Figure 5. PP1 is a Ca2+- responsive enzyme involved in T. gondii egress and invasion. (A) The Ca2+- dependent stabilization of PP1 (TGGT1_310700) in each mass spectrometry experiment. (B) Immunoblotting for endogenously tagged PP1- mNG- Ty at different Ca2+ concentrations and thermal challenge at 58 °C. Abundance is calculated relative to the band intensity at 0 µM Ca2+ and scaled. Points in shades of gray represent different replicates. A dose-r esponse curve was calculated for the mean abundances. (C) Parasites expressing endogenously tagged PP1-m NG egress after treatment with 500 µM zaprinast or 4 µM A23187. Arrows show examples of PP1 enrichment at the apical end. Time after treatment is indicated in m:ss. (D) Rapid regulation of PP1 by endogenous tagging with mAID- HA. IAA, Indole- 3- acetic acid (IAA). (E) PP1- mAID-H A visualized in fixed intracellular parasites by immunofluorescence after 3 hr of 500 µM IAA or vehicle treatment. Hoechst and anti- CDPK1 are used as counterstains (Waldman et al., 2020). (F) PP1- mAID-H A depletion, as described in (E), monitored by immunoblotting. The expected molecular weights of PP1-m AID-H A and CDPK1 are 48 and 65 kDa, respectively. (G) Plaque assays of 1,000 TIR1 and PP1- mAID- HA parasites infected onto a host cell monolayer and allowed to undergo repeated cycles of invasion, replication, and lysis for 7 days in media with or without IAA. (H) The number of parasites per vacuole measured for PP1-m AID- HA and the TIR1 parental strain after 24 hr of 500 µM IAA treatment. Mean counts (n=3) are expressed as a percentage of all vacuoles counted. (I) Invasion assays PP1-m AID-H A or TIR1 parental strains treated with IAA or vehicle for 3 hr. Parasites were incubated on host cells for 60 min prior to differential staining of intracellular and extracellular parasites. Parasite numbers were normalized to host cell nuclei for each field. Means graphed for n=5 biological replicates (different shapes), Welch’s t- test. (J) Parasite egress stimulated with 500 µM zaprinast or 8 µM A23187 following 3 h of treatment with vehicle or IAA. Egress was monitored by the number of host cell nuclei stained with DAPI over time and was normalized to egress in the vehicle-t reated strain. Mean ±S.D. graphed for n=3 biological replicates. The online version of this article includes the following source data for figure 5: Source data 1. This file contains the source data that was quantified to make the graph presented in Figure 4. Source data 2. This file contains the source data that was quantified to make the graph presented in Figure 4. Source data 3. This file contains the annotated source data that was quantified to make the graph presented in Figure 5. Source data 4. This file contains the source data that was presented in Figure 5. Source data 5. This file contains the annotate source data that was presented in Figure 5. Herneisen et al. eLife 2022;11:e80336. DOI: https://doi.org/10.7554/eLife.80336 15 of 41 Research article Cell Biology | Microbiology and Infectious Disease treatment failed to induce the same relocalization patterns (Figure 5C and Video 2). The dynamics of PP1 activity thus appear enhanced by cGMP signaling within the parasites, which is upstream of Ca2+ (Brown et al., 2016; Sidik et al., 2016b) and was similarly observed in Plasmodium (Paul et al., 2020). The parasite apex and pellicle are hotspots for the signaling that potentiates motility and invasion, so the relocalization of PP1 suggests it may function in these processes. We created strains expressing the PP1 catalytic subunit with an endogenous C- terminal auxin- inducible degron (AID) for rapid conditional knockdown (Figure  5D). Conditional degradation of PP1, following treatment with indole-3 -a cetic acid Video 1. Representative image series of parasites (IAA) was confirmed by immunofluorescence expressing endogenously tagged PP1-m NG following and immunoblotting (Figure  5E and F). Para- treatment with 500 µM zaprinast. sites depleted of PP1 failed to form plaques https://elifesciences.org/articles/80336/figures#video1 (Figure 5G), implicating the phosphatase in one or more essential functions during the lytic cycle (Sidik et al., 2016a). Even in the absence of IAA, the PP1- AID strain formed small plaques, indicating substantial hypomorphism. PP1- AID parasites exhibited slower replication than untagged parasites, and this effect was exacerbated with IAA treat- ment (Figure 5H). Parasites depleted of PP1 had a reduced invasion efficiency (Figure 5I), although this effect was modest and subject to technical variation, likely arising from the hypomorphism of the C- terminally tagged strain. Parasites treated with IAA egressed more slowly than untreated parasites when stimulated with zaprinast; however, egress kinetics were indistinguishable with a Ca2+ iono- phore agonist (Figure 5J). Together, these results suggest that PP1 holoenzymes function at multiple steps in the lytic cycle. At least one holoenzyme relocalizes when parasite cGMP/Ca2+ pathways are stimulated. Because parasites lacking PP1 exhibit delays specifically in zaprinast-i nduced egress, we hypothesize that the peripheral holoenzyme enhances Ca2+ signaling. A PP1 holoenzyme dephosphorylates signaling enzymes at the parasite periphery PP1 is one of the major serine/threonine phos- phatases in eukaryotic cells (Brautigan and Shenolikar, 2018). As the catalytic subunit relo- calizes during cGMP/Ca2+- stimulated transitions in apicomplexans, we hypothesized that PP1 dephosphorylates crucial targets during egress and motility. To identify putative PP1 holoen- zyme targets, we first treated intracellular PP1- AID parasites with IAA or vehicle for 3  hr prior to mechanically releasing parasites for analysis. We performed sub- minute phosphoproteomics by resuspending the extracellular parasites in a zaprinast solution followed by denaturation in SDS to stop enzymatic activity after 0, 10, 30, and 60  s. Multiplexing with TMTpro reagents followed by phosphopeptide enrichment allowed us to compare the zaprinast time courses with Video 2. Representative image series of parasites or without PP1- depletion in biological duplicate expressing endogenously tagged PP1-m NG following (Figure 6A). Analysis of an unenriched fraction of treatment with 4 µM A23187. the proteome revealed significant depletion of https://elifesciences.org/articles/80336/figures#video2 Herneisen et al. eLife 2022;11:e80336. DOI: https://doi.org/10.7554/eLife.80336 16 of 41 Research article Cell Biology | Microbiology and Infectious Disease Figure 6. The PP1- dependent phosphoproteome. (A) Schematic of the phosphoproteomics time course. PP1- AID parasites were treated with IAA or vehicle for 3 hr. Extracellular parasites were then treated with zaprinast, and samples were collected during the first minute after stimulation. The experiment was performed in biological replicate (R1 and R2). (B) Five clusters were identified with respect to phosphopeptide dynamics and PP1- dependence. (C) Examples of phosphopeptides dynamically regulated by zaprinast and exhibiting PP1-d ependent dephosphorylation. (D) GO terms enriched among phosphopeptides, grouped by cluster. Gene ratio is the proportion of proteins with the indicated GO term divided by the total number of proteins belonging to each cluster. Significance was determined with a hypergeometric test; only GO terms with p<0.05 and represented by more than one protein are shown. Redundant GO terms were removed. Cluster 1 lacked enough peptides for enrichment analysis. The online version of this article includes the following source data and figure supplement(s) for figure 6: Figure supplement 1. Extended analysis of PP1 phosphoproteomics experiment. Figure supplement 1—source data 1. This file contains the source data that was presented in Figure 6—figure supplement 1. Figure supplement 1—source data 2. This file contains the source data that was presented in Figure 6—figure supplement 1. Figure supplement 1—source data 3. This file contains the annotated source data that was presented in Figure 6—figure supplement 1. PP1, which we confirmed in parallel by immunoblot (Figure 6—figure supplement 1). The hypo- morphism of the PP1-A ID strain and reduced parasite yield resulted in a phosphoproteome with lower coverage: 6,916 phosphopeptides with quantification values. To identify peptides exhibiting PP1-d ependent regulation, we selected those exhibiting a twofold difference in abundance between vehicle- and IAA-t reated samples in both replicates for at least one time point. In total, 757 peptides passed this threshold. To focus on zaprinast- dependent changes, we clustered these peptides on the basis of their abun- dance relative to the earliest time point of stimulation. The peptides fit into five clusters with respect Herneisen et al. eLife 2022;11:e80336. DOI: https://doi.org/10.7554/eLife.80336 17 of 41 Research article Cell Biology | Microbiology and Infectious Disease to PP1- dependence and kinetics (Figure 6B). Clusters 1, 2, and 3 contained peptides increasing in abundance upon stimulation in PP1-d epleted parasites. The peptides belonging to cluster 1 gener- ally increase in abundance with zaprinast treatment, which occurs more rapidly with PP1 depletion. By contrast, the abundances of the 614 peptides belonging to clusters 2 and 3 were elevated in the absence of PP1, suggesting that PP1 antagonizes these phosphorylation events. Under normal conditions, peptides in cluster 4 decreased sharply in abundance between 10 and 30 s of zaprinast treatment and recovered by 60 s; however, these peptides did not change in abundance when PP1 was depleted. Peptides in cluster 5 increased gradually between 10 and 30 s and decreased to basal levels by 60 s; when PP1 is depleted, these peptides exhibit a delay. The effect of PP1 disruption on the phosphoproteome was pervasive, which may reflect disrup- tion of numerous PP1 holoenzymes. Many of the examined phosphopeptides exhibited substantial basal hyperphosphorylation in the absence of PP1 (Figure 6C). Given the likely pleiotropy of cata- lytic subunit depletion, we focused our analysis on perturbed pathways rather than individual targets (Figure 6D). Cluster 1 did not possess enough peptides for enrichment analysis. Cluster 2 was enriched in phosphoproteins functioning in transmembrane transport, including P-t ype ATPases and ABC trans- porters. Cluster 3 contained both the guanylyl and adenylyl cyclases and was further enriched in putative sodium- hydrogen exchangers and tubulin- tyrosine ligases. Cluster 4 was solely enriched for proteins involved in ubiquitin transfer, including TGGT1_280660, an uncharacterized HECT domain- containing protein (Figure  6C). Cluster 5 phosphoproteins were notably involved in cytoskeletal motor activity, actin binding, and RNA binding. Numerous apical proteins exhibited PP1-d ependent phosphorylation (Table 1), including AC7, CIP1, and two hypothetical proteins—TGGT1_258090 and TGGT1_320030, localized to the conoid base and the second apical polar ring, respectively (Koreny et al., 2021). The protein abundances did not vary between vehicle and IAA treatment, indicating that the phosphopeptide abundance changes arose from dynamic covalent modifications. The PP1- dependent phosphoproteome therefore supports the existence of apical and peripheral PP1 holoen- zymes, as observed by live microscopy. PP1 activity is important for Ca2+ entry to enhance the kinetics of egress The phosphoproteomics data pointed to ion homeostasis dysregulation in the absence of PP1 (Figure 6D). Numerous transporters are hyperphosphorylated when PP1 is depleted, including Ca2+ ATPases, proton transporters, and MFS and ABC- family transporters. We hypothesized that aberrant Ca2+ mobilization may underlie the phenotypes observed in the motile stages of the lytic cycle upon PP1 depletion (Figure 5), as Ca2+ release and signaling are required for egress and invasion. To monitor Ca2+ mobilization in parasites, we introduced the genetically encoded Ca2+ indicator GCaMP6f into a defined chromosomal locus (Chen et al., 2013; Herneisen et al., 2020). By tracking the vacuole fluorescence of PP1- depleted parasites, we observed delayed Ca2+ mobilization and egress following zaprinast treatment (Figure 7A–C, Figure 7—figure supplement 1, and Video 3). Despite the delay, IAA- treated parasites attained the same increase in Ca2+ levels prior to egress as vehicle-t reated para- sites (Figure 7—figure supplement 1). By contrast, Ca2+ mobilization upon A23187 treatment was unperturbed in PP1- depleted parasites (Figure 7D–F, Figure 7—figure supplement 1 and Video 4). Motivated by recent experiments suggesting that T. gondii exhibits Ca2+- activated Ca2+ entry (Vella et al., 2021), we loaded PP1- AID parasites with the ratiometric Ca2+ indicator Fura2- AM and observed both reduced Ca2+ entry and lowered resting cytoplasmic [Ca2+] (Figure 7G–I and Figure 7—figure supplement 1). The integration of our global datasets therefore identified PP1 as a Ca2+- responsive enzyme and precipitated discovery of its role in Ca2+ entry. Discussion Apicomplexan signaling pathways have largely been characterized through a combination of genetic manipulation, pharmacological perturbation, and physiological observation (Bisio and Soldati-F avre, 2019; Brown et al., 2020; Lourido and Moreno, 2015). Such approaches have been sufficient to docu- ment feedback loops between cyclic nucleotide signaling and Ca2+- store release and influx. However, the molecular mechanisms linking these pathways have remained obscure, as genomic annotation or bioinformatic analysis are complicated by the evolutionary divergence of apicomplexans. Unbiased, Herneisen et al. eLife 2022;11:e80336. DOI: https://doi.org/10.7554/eLife.80336 18 of 41 Research article Cell Biology | Microbiology and Infectious Disease A zaprinast B C G 0 s 70 s 110 s 1.8 300 300 TIR1 +IAA TIR1 +IAA PP1-AID +veh. p=0.018 PP1-AID +IAA 1.6 PP1-AID +IAA p=0.013 p=0.12 10 µm 200 200 1.4 1.2 100 100 1.0 1.8 mM Ca2+ 0.8 0 0 0 50 100 150 200 TIR1 PP1−AID PP1−AID 0 100 200 300 400 time (s) + IAA + vehicle + IAA time (s) D A23187 E F H I 0 s 70 s 110 s 1.8 300 120 10 TIR1 +IAA p=0.31 p=0.0002 p=0.017PP1-AID +veh. p=0.76 1.6 PP1-AID +IAA 100p=0.047 8 10 µm 200 80 1.4 6 60 1.2 4 100 40 1.0 20 2 0.8 0 0 0 0 50 100 150 200 TIR1 PP1−AID PP1−AID TIR1 PP1-AID TIR1 PP1-AID time (s) + IAA + vehicle + IAA +IAA Figure 7. PP1 ensures rapid Ca2+ mobilization prior to zaprinast- induced egress. (A) Selected frames of time- lapse images of PP1-A ID parasites expressing the genetically encoded Ca2+ indicator GCaMP following treatment with 500 µM zaprinast. (B) Normalized GCaMP fluorescence of individual vacuoles was tracked after zaprinast treatment and prior to egress (opaque lines). The solid line represents the mean normalized fluorescence of all vacuoles across n=3 biological replicates. (C) The time to maximum normalized fluorescence of individual vacuoles after zaprinast treatment. Different replicates are shown in different shades of gray. Small points correspond to individual vacuoles; large points are the mean for each replicate. p Values were calculated from a two- tailed t- test. (D) Selected frames of time-l apse images of PP1-A ID parasites expressing the genetically encoded Ca2+ indicator GCaMP following treatment with 4 µM A23187. (E) Normalized GCaMP fluorescence of individual vacuoles was tracked after A23187 treatment and prior to egress (opaque lines). The solid line represents the mean normalized fluorescence of all vacuoles across n=3 biological replicates. (F) The time to maximum normalized fluorescence of individual vacuoles after A23187 treatment. Small points correspond to individual vacuoles; large points are the mean for each replicate. p Values were calculated from a two- tailed t- test. (G) Fluorescence intensity of Fura2/AM- loaded TIR1 or PP1- AID parasites treated with IAA for 5 hr before and after addition of the 1.8 mM Ca2+. Representative traces from three biological replicates. (H) Resting cytoplasmic [Ca2+] prior to incubation in buffers with elevated [Ca2+]. p Values were calculated from an ANOVA. (I) The rate of Ca2+ entry in the first 20 s after addition of 1.8 mM Ca2+ to parasites p values were calculated from an ANOVA. Entry rates following addition of other concentrations are shown in Figure 7—figure supplement 1. The online version of this article includes the following figure supplement(s) for figure 7: Figure supplement 1. PP1-A ID parasites egress, as quantified by video microscopy. Video 3. Representative image series of PP1- AID Video 4. Representative image series of PP1- AID parasites expressing the genetically encoded Ca2+ parasites expressing the genetically encoded Ca2+ indicator GCaMP following treatment with 500 µM indicator GCaMP following treatment with 4 µM zaprinast. A23187. https://elifesciences.org/articles/80336/figures#video3 https://elifesciences.org/articles/80336/figures#video4 Herneisen et al. eLife 2022;11:e80336. DOI: https://doi.org/10.7554/eLife.80336 19 of 41 PP1–AID TIR1 PP1–AID TIR1 +IAA +vehicle +IAA +IAA +vehicle +IAA normalized fluorescence normalized fluorescence time to max (s) time to max (s) 2+ resting cytosolic Ca2+ (nM) Δ[Ca ]cyt (nM) Ca2+ entry rate (nM/s) Research article Cell Biology | Microbiology and Infectious Disease functional proteomics approaches have allowed us to map a Ca2+-r esponsive proteome de novo. Our resource places proteins already under study within Ca2+ signaling pathways and identifies hubs serving as pressure points for the perturbation of multiple signaling processes. Our time- resolved phosphoproteomic data adds molecular depth to known changes in parasite physiology. Through the use of zaprinast we stimulate the endogenous regulation of Ca2+ stores by PKG. Our phosphoproteomic analysis therefore includes Ca2+-r esponsive pathways, as well as any other pathways downstream of cGMP or PKG, which may include Ca2+-i ndependent processes ( Bullen et al., 2019). The first cluster of phosphoproteins to respond to zaprinast stimulation includes second messenger amplifiers: cyclic nucleotide phosphodiesterases (PDE1 and PDE2); lipid kinases, phosphatases, and phospholipases; and the guanylyl cyclase (Bisio et al., 2019; Brown and Sibley, 2018). This cluster also contains early Ca2+ sensors and transducers, such as TgCDPK2A and TgCDPK7 (Bansal et al., 2021), which additionally responded to Ca2+ in our thermal profiling experiments. The next two clusters reveal rewiring of small molecule homeostasis, with non-C a2+ channels and trans- porters phosphorylated first. Parasites have long been reported to respond to extracellular H+ and K+ levels (Moudy et al., 2001; Roiko and Carruthers, 2013), although the proteins responsible and sequence of events transducing signals across the parasite plasma membrane remained obscure. Our proteomic datasets identified dynamic phosphorylation of ATPase2, NHE1, and two uncharac- terized putative Ca2+- activated K+ channels. ATPase2 was recently reported to mediate essential lipid homeostasis at the parasite plasma membrane (Chen et al., 2021). Mutants lacking NHE1 exhibited elevated cytosolic Ca2+, yet were less susceptible to Ca2+ ionophore-i nduced egress (Arrizabalaga et al., 2004), suggesting that monovalent cation and Ca2+ homeostasis mechanisms are intertwined. Indeed, the third wave of phosphorylation targets Ca2+ and other divalent cation ATPases, including TgA1, TgCuTP, and other P-t ype ATPases with uncharacterized ion specificity. Transporters belonging to these clusters are noted TgCDPK3 substrates, including ATP4, SCE1, and ApiAT5- 3 (Treeck et al., 2014; McCoy et al., 2017; Dominicus et al., 2021); disruption of these transporters failed to alter cytosolic Ca2+ levels (McCoy et al., 2017; Parker et al., 2019; Wallbank et al., 2019). A cluster of dynamic dephosphorylation events highlights the need for phosphatases to reverse modifications as parasites transition from one phase to another (Brautigan and Shenolikar, 2018; Yang and Arrizaba- laga, 2017). These phosphoproteins were also enriched for channel or transporter functions, and in several cases also contained sites regulated by kinases. Apicomplexans must maintain homeostasis as ion gradients are inverted when parasites transition from intracellular to extracellular, likely requiring extreme rewiring of cellular physiology. Our phosphoproteomics dataset reveals that parasite ion homeostasis pathways are hard- wired into second messenger signal transduction, perhaps in antici- pation of environmental challenge. We can infer mechanisms of crosstalk between second messenger signaling networks, different classes of enzymes, and post-t ranslational modifications from our global proteomic profiles. Proteins regulated by several kinases and phosphatases with different kinetics may function to concentrate these enzymes in local signaling hotspots acting as signaling scaffolds. Several of these dynami- cally regulated proteins themselves contain second messenger–sensing domains, such as EF hands (TGGT1_216620 and TgCDPK2A) and cyclic nucleotide–binding domains (PDE2 and AAP2). Intrigu- ingly, both subunits of PKA, the best- characterized cAMP receptor in T. gondii, are stabilized by Ca2+. PKA has been proposed to negatively regulate Ca2+ signaling in parasites (Jia et al., 2017; Uboldi et al., 2018). Here, we additionally observe Ca2+ regulating PKA, although the molecular basis of this stability change remains to be established. Finally, the stimulated phosphoproteome exhibits exten- sive crosstalk with GTPase regulation and ubiquitin transfer pathways, suggesting that there remain vast and unexplored signaling networks, in addition to phosphorylation, that remodel cellular states in T. gondii. The Ca2+ responsiveness of PP1 revealed by our TPP analysis was surprising. In most organisms, the dominant Ca2+-r egulated phosphatase is calcineurin (Shi, 2009). Although PP1 holoenzymes have been reported to modulate Ca2+ channels, the Ca2+- responsive behavior of PP1 has not been reported in any other eukaryote. The function of the phosphatase had not been directly studied in T. gondii, although PP1 inhibitors have been shown to block invasion of host cells (Delorme et al., 2002) and recent experiments in P. falciparum suggest that PP1 is required for an as-y et undetermined step in the kinetic phase of infection (Paul et al., 2020). Guided by our global experiments, we endoge- nously tagged the C terminus of PP1 and observed diverse localization patterns. Numerous regulatory Herneisen et al. eLife 2022;11:e80336. DOI: https://doi.org/10.7554/eLife.80336 20 of 41 Research article Cell Biology | Microbiology and Infectious Disease subunits control the activity, localization, and specificity of the single catalytic subunit of the phos- phatase (Yang and Arrizabalaga, 2017), giving rise to functionally distinct PP1 holoenzymes that play conserved roles in cell cycle progression (Bertolotti, 2018). Stimulating the transition from the replicative to the kinetic phase by treating parasites with the phosphodiesterase inhibitor zaprinast led to rapid enrichment of PP1 at the parasite periphery and apical end, hotspots of parasite signaling. Notably, PP1 was recently reported to relocalize to the apical complex in Plasmodium berghei ooki- netes (Zeeshan et al., 2021). Although the proposed role of PP1 in ookinetes was to set cell polarity, given our phenotypic observations in Toxoplasma, we propose instead that the apical PP1 holoen- zyme serves apicomplexan- specific, Ca2+- responsive functions supporting the parasitic lifestyle (Park et al., 2019; Paul et al., 2015; Philip and Waters, 2015). Indeed, parasites conditionally depleted of PP1 exhibited reduced invasion and delayed host cell lysis in response to zaprinast. Many parasite PP1 holoenzymes are as yet undefined. A recent phosphoproteome of P. falciparum blood- stage parasites depleted of PP1 identified ~5,000 phosphosites from ~1,000 proteins (Paul et al., 2020). Despite complications in determining direct substrates, the study proposed PP1 regu- lation of a HECT E3 ligase and guanylate cyclase upstream of PKG, providing a potential functional link between PP1 and motility. Our experimental design benefited from higher temporal resolution and the ability to directly trigger the signaling pathways in question, providing a more compre- hensive analysis of phosphatase activity during stimulated motility. We uncovered two patterns of PP1-d ependent regulation: (i) sites hyperphosphorylated only in the absence of PP1 (clusters 2 and 3), and (ii) zaprinast-d ependent sites no longer regulated in the absence of PP1 (cluster 4). The guanylyl cyclase and possible ortholog of the HECT E3 ligase identified in the P. falciparum study both included PP1- dependent phosphosites in our experiment, suggesting that these PP1- regulated cascades are important and conserved in apicomplexan parasites. The guanylyl cyclase belongs to one of the hyperphosphorylated clusters. If guanylyl cyclase function is compromised in the absence of PP1, accumulation of cGMP to the levels required for egress may be delayed following inhibition of phosphodiesterases, partially explaining the delayed zaprinast response observed in PP1-d epleted parasites (Paul et al., 2020). Additional work probing the cGMP pathway—for example by directly monitoring cyclic nucleotide levels or leveraging mutants expressing conditional depletion alleles of phosphodiesterases, cyclases, and kinases—is required to precisely place PP1 in the network that primes parasites for the kinetic phase of the infection cycle. Moreover, the set of phosphoproteins normally dephosphorylated by PP1 during zaprinast treatment was uniquely enriched in ubiquitin- transfer enzymes. Prior proteomic enrichment studies observed an abundance of ubiquitin modifi- cations in cell-c ycle- regulated proteins, including the peripheral parasite inner membrane complex (Silmon de Monerri et al., 2015). However, the ubiquitination pathway is largely unexplored in the context of apicomplexan motility. HECT E3 ligase inhibitors reduced merozoite egress in P. falciparum (Paul et al., 2020). Our proteomic datasets reveal that the phosphorylation states of ubiquitin trans- ferases and hydrolases change within 60 s of cGMP elevation, suggesting that dynamic ubiquitination intersects with dynamic phosphorylation to regulate transitions between intracellular and extracellular states. The enrichment of transporters in the basally dysregulated sites motivated us to explore PP1 pheno- types associated with Ca2+ mobilization and entry. We observed that PP1 activity enhances Ca2+ entry, which itself was reported to enhance parasite motility (Pace et al., 2014; Vella et al., 2021). Synthe- sizing the phosphoproteomic, phenotypic, and thermal profiling data, which reports a micromolar- range EC50 of PP1 for Ca2+, we propose that the Ca2+- regulated PP1 holoenzyme assembles following an initial wave of Ca2+ store release, and then inhibits Ca2+ uptake, thereby contributing to the overall elevation of cytoplasmic Ca2+. The identification of specific PP1 holoenzymes and genetic manipulation of regulatory subunits is required to test this model; however, such negative feedback loops between Ca2+- stimulated kinases and Ca2+ uptake have been reported for TgCDPK1 and TgCDPK3 (Herneisen et al., 2020; Stewart et al., 2017). Alternatively, PP1 itself may change membrane association or structure. The Ca2+ signaling field has benefited from decades of research in other model organisms, in which pathway connections were incrementally extended based on prior knowledge about the system (Berridge et al., 2000; Clapham, 2007; Luan and Wang, 2021). By contrast, we developed a global mass spectrometry- based approach to systematically identify the protein components of Ca2+ signaling pathways on the basis of dynamic phosphoregulation and biochemical thermal stability. We applied this technology to apicomplexans, which are widespread eukaryotic human pathogens Herneisen et al. eLife 2022;11:e80336. DOI: https://doi.org/10.7554/eLife.80336 21 of 41 Research article Cell Biology | Microbiology and Infectious Disease whose signaling pathways remain largely unmapped due to their evolutionary divergence from model organisms (Lourido and Moreno, 2015). Prior proteomic surveys investigating the contribution of individual kinases (Alam et al., 2015; Balestra et al., 2021; Brochet et al., 2014; Jia et al., 2017; Patel et al., 2019; Treeck et al., 2014) or phosphatases (Paul et al., 2020; Yang et al., 2019) to global phosphoregulation lacked the temporal resolution to map transitions between the replicative and kinetic phases of the apicomplexan infection cycle, which occur in a matter of seconds. Advances in proteomics technologies, such as enhanced sensitivity and improved multiplexing methods, are now being leveraged to monitor the sub-m inute processes of exflagellation and egress (Dominicus et al., 2021; Invergo et al., 2017). Second- messenger signaling pathways underpin these transitions at distinct developmental stages in diverse apicomplexans, with many kinases and phosphatases func- tioning at multiple steps in the parasitic life cycle. Meta-a nalyses of orthogonal proteomic profiles, such as a phosphoproteome tracking ionophore- induced egress (Dominicus et  al., 2021), could reveal signatures of signaling subnetworks. A recent study points to extensive regulation of enzymes in the cyclic nucleotide pathways following treatment with calcium ionophore, which likely bypasses some of the regulated steps of calcium release and uptake from parasite stores—PDE1, PDE2, and the guanylyl and adenylyl cyclases were all heavily regulated (Dominicus et al., 2021). We identified some of these enzymes (PDE2, GC, and the CNB protein AAP2) as candidate signaling platforms, as in our study they belonged to several kinetically resolved clusters, suggesting regulation by multiple kinases or phosphatases. In principle, our approach can be applied to other post-t ranslational modifi- cations, natural ligands (Lim et al., 2018; Sridharan et al., 2019), and organisms (Jarzab et al., 2020) to establish—and in some cases re- evaluate—the topology of complex signaling pathways. Materials and methods Key resources table Reagent type (species) or resource Designation Source or reference Identifiers Additional information Strain, strain background (T. gondii) TIR1 PMID:28465425 RH/TIR1/∆KU80/∆HXGPRT Strain, strain RH/TIR1/pTUB1-G CaMP6f- background (T. gondii) TIR1/GCaMP6f PMID:35484233 3'DHFR/∆KU80/∆HXGPRT Strain, strain background (T. gondii) DiCre PMID:31577230 RH/∆KU80::DiCre_T2A/∆HXGPRT Strain, strain RH/TIR1/∆KU80/∆HXGPRT/pMIC2- MIC2- GLuc- background (T. gondii) TIR1/MIC2- GLuc-P 2A-G CaMP6f This paper myc- P2A-G CaMP6f Strain, strain RH/TIR1/pTUB1-G CaMP6f- 3'DHFR/∆KU80/PKA background (T. gondii) PKA C1- Ty (Figure 4) This paper TGGT1_226030 C1- mCherry-V 5-m AID-T y/PKA R- V5-3 HA Strain, strain RH/TIR1/pTUB1- GCaMP6f- 3'DHFR/∆KU80/ background (T. gondii) Eps15- HA (Figure 4) This paper TGGT1_227800 Eps15-V 5- mCherry- mAID- HA Strain, strain RH/TIR1/∆KU80/∆HXGPRT/TGGT1_286470- background (T. gondii) zinc finger- HA (Figure 4) This paper TGGT1_286710 V5- 3HA RH/TIR1/pTUB1- GCaMP6f- Strain, strain 3'DHFR/∆KU80/∆HXGPRT/TGGT1_309290-V 5- background (T. gondii) hypothetical-V 5 (Figure 4) This paper TGGT1_309290 mNeonGreen- mAID-T y Strain, strain RH/∆KU80::DiCre_T2A/∆HXGPRT/RON13- background (T. gondii) RON13- HA (Figure 4) This paper TGGT1_321650 HA- U1 Strain, strain RH/TIR1/∆KU80/∆HXGPRT/pMIC2-M IC2- GLuc- background (T. gondii) PP1-A ID (Figure 5) This paper TGGT1_310700 myc-P 2A- GCaMP6f/PP1- V5- mAID- HA Strain, strain PP1-A ID/TIR1/MIC2-G Luc- P2A- RH/TIR1/∆KU80/∆HXGPRT/pMIC2- MIC2-G Luc- background (T. gondii) GCaMP6f (Figure 7) This paper TGGT1_310700 myc- P2A- GCaMP6f/PP1- V5- mAID- HA Human Foreskin Fibroblasts Cell line (Homo sapiens) (HFFs) ATCC SCRC- 1041 Guinea pig monoclonal anti- Antibody CDPK1 Covance Custom antibody IF (1/10000), WB (1/40000) Continued on next page Herneisen et al. eLife 2022;11:e80336. DOI: https://doi.org/10.7554/eLife.80336 22 of 41 Research article Cell Biology | Microbiology and Infectious Disease Continued Reagent type (species) or resource Designation Source or reference Identifiers Additional information Developmental Studies Mouse monoclonal anti- TUB1 Hybridoma Bank at the Antibody (clone 12G10) University of Iowa RRID:AB_1157911 WB (1/5000) Mouse monoclonal anti- Ty1 Antibody (clone BB2) PMID:8813669 IF (1/1000), WB(1/2000) Antibody Rabbit polyclonal anti- HA Invitrogen Invitrogen:71–5500 WB (1/1000) Mouse monoclonal anti- MIC2 Antibody (6D10) PMID:10799515 WB (1:2000) Antibody Mouse monoclonal anti- V5 Invitrogen Invitrogen:R960- 25 IF (1:1000), WB (1:2000) Antibody Rabbit polyclonal PMID:18312842 WB (1:2000) Mouse monoclonal anti- HA Antibody (16B12) Biolegend Biolegend:901533 IF (1:1000) Antibody Mouse polyclonal anti- SAG1 PMID:3183382 WB (1/1000) Alexa Fluor 594 polyclonal goat Antibody anti-g uinea pig Life Technologies Life Technologies:A11076 IF (1/1000) Alexa Fluor 488 polyclonal goat Antibody anti- mouse Life Technologies Life Technologies:A11029 IF (1/1000) IRDye 800CW polyclonal Goat anti- Mouse IgG1- Specific Antibody Secondary Antibody LICOR LICOR:926–32350 WB (1/10000) IRDye 680LT polyclonal Goat anti- Antibody Mouse IgG Secondary Antibody LICOR LICOR:926–68020 WB (1/10000) IRDye 800CW polyclonal Donkey anti-G uinea Pig IgG Secondary Antibody Antibody LICOR LICOR:926–32411 WB (1/10000) IRDye 680RD polyclonal Donkey anti- Guinea Pig IgG Secondary Antibody Antibody LICOR LICOR:926–68077 WB (1/10000) IRDye 800CW polyclonal Goat anti- Rabbit IgG Secondary Antibody Antibody LICOR LICOR:926–32211 WB (1/10000) IRDye 680LT polyclonal Goat anti- Antibody Rabbit IgG Secondary Antibody LICOR LICOR:926–68021 WB (1/10000) Chemical compound, drug Hoechst Santa Cruz Santa Cruz:sc-3 94039 IF (1/20000) Chemical compound, drug Prolong Diamond Thermo Fisher Thermo Fisher:P36965 Chemical compound, drug zaprinast Calbiochem Calbiochem:684500 Chemical compound, drug A23187 Calbiochem Calbiochem:100105 Commercial assay or kit S- trap micro Protifi Protific:C02- micro-8 0 TMT10plex Isobaric Label Commercial assay or kit Reagent Set Thermo Fisher Scientific Thermo Fisher Scientific:90111 TMTpro 16plex Label Reagent Commercial assay or kit Set Thermo Fisher Scientific Thermo Fisher Scientific:A44522 High-S elect TiO2 Commercial assay or kit Phosphopeptide Enrichment Kit Thermo Fisher Scientific Thermo Fisher Scientific:A32993 High- Select Fe- NTA Commercial assay or kit Phosphopeptide Enrichment Kit Thermo Fisher Scientific Thermo Fisher Scientific:A32992 Continued on next page Herneisen et al. eLife 2022;11:e80336. DOI: https://doi.org/10.7554/eLife.80336 23 of 41 Research article Cell Biology | Microbiology and Infectious Disease Continued Reagent type (species) or resource Designation Source or reference Identifiers Additional information Pierce High pH Reversed- Phase Commercial assay or kit Peptide Fractionation Kit Thermo Fisher Scientific Thermo Fisher Scientific:84868 Calcium Calibration Buffer Kit #1, Commercial assay or kit zero and 10 mM CaEGTA Life Technologies Life Technologies:C3008MP Hydrophobic Sera- Mag Speed Commercial assay or kit Beads GE Healthcare GE Healthcare:65152105050250 Hydrophilic Sera- Mag Speed Commercial assay or kit Beads GE Healthcare GE Healthcare:45152105050250 Recombinant DNA All plasmids used in this study are reagent listed in Supplementary file 7 All primers and oligonucleotides Sequence-b ased used in this study are listed in reagent Supplementary file 7 Software, algorithm Proteome Discoverer 4.2 Thermo Fisher R Foundation for Statistical Software, algorithm R version 4.0 Computing Software, algorithm mineCETSA version 1.1.1 Dziekan et al., 2020 https://github.com/nkdailingyun/mineCETSA Software, algorithm Prism 8 GraphPad Software, algorithm HHPRED PMID:29258817 Other Halt protease inhibitor Thermo Fisher Thermo Fisher:87786 Materials and Methods: lysis buffer Halt protease and phosphatase Other inhibitor Thermo Fisher Thermo Fisher:PI78440 Materials and Methods: lysis buffer Other Benzonase Sigma Aldrich Sigma Aldrich:E1014 Materials and Methods: lysis buffer Reagents and cell culture T. gondii parasites of the type I RH strain ∆ku80/∆hxgprt genetic background (ATCC PRA- 319) were grown in human foreskin fibroblasts (HFFs, ATCC SCRC-1 041) maintained in DMEM (GIBCO) supple- mented with 3% inactivated fetal calf serum (IFS) and 10 μg/mL gentamicin (Thermo Fisher Scientific). When noted, DMEM was supplemented with 10% IFS and 10 μg/mL gentamicin. Parasite transfection and strain construction Genetic background of parasite strains Existing T. gondii RH strains were used as genetic backgrounds for this study. All strains contain the ∆ku80∆hxgprt mutations to facilitate homologous recombination (Huynh and Carruthers, 2009). TIR1/RH expresses the TIR1 ubiquitin ligase (Brown et al., 2017). TIR1/pTUB1-G CaMP additionally expresses a GCaMP6f transgene (Smith et  al., 2022). DiCre_T2A/RH uses the same promoter to express both subunits of Cre, separated by T2A skip peptides (Hunt et al., 2019). TIR1/pMIC2-MIC2Gluc-P2A-GCaMP6f A new TIR1/RH strain expressing GCaMP6f under different regulatory elements was constructed as a part of this study. The sequence pMIC2_MIC2Gluc_myc_P2A_GCaMP6f_DHFR3′UTR was amplified with primers P1 and P2 from plasmid Genbank ON312866 to yield a repair template with homology to the 5′ and 3′ ends of a defined, neutral genomic locus (Markus et al., 2019). Approximately 1×107 extracellular TIR1 parasites were transfected with 20  µg gRNA/Cas9-e xpression plasmid pBM041 (GenBank MN019116.1) and 10 µg repair template. GFP-p ositive clones were isolated by limiting dilution following fluorescence- activated cell sorting. Endogenous tagging of candidate genes Genes were endogenously tagged using the previously described high- throughput (HiT) tagging strategy (Smith et al., 2022). Cutting units specific to each candidate were purchased as IDT gBlocks Herneisen et al. eLife 2022;11:e80336. DOI: https://doi.org/10.7554/eLife.80336 24 of 41 Research article Cell Biology | Microbiology and Infectious Disease (P3- P7) and assembled with the following empty HiT vector backbones: pALH083 (V5-m Cherry- mAID- Ty; GenBank: ON312867), pALH052 (V5- mCherry-m AID- HA; GenBank: OM863784), pALH047 (V5- 3HA; GenBank: ON312868), pGL015 (V5- mNG-m AID- Ty; GenBank OM640005), pEC341 (HA- loxP- U1; GenBank: OM863788), and pALH086 (V5-m AID- HA; GenBank ON312869). Between 20 and 50 µg of each vector was BsaI- linearized and co- transfected with 50 μg of the pSS014 Cas9- expression plasmid (GenBank: OM640002). Vectors targeting TGGT1_226030 (PKA-C 1) and TGGT1_227800 (Eps15) were transfected into TIR1/pTUB1-G CaMP or TIR1/pMIC2-M IC2Gluc- P2A- GCaMP6f. Vectors targeting TGGT1_286710 and TGGT1_309290 were transfected into TIR1/RH. Vectors targeting TGGT1_321650 were transfected into DiCre/RH. Parasite populations were selected for 1 week in standard media with 1 µM pyrimethamine or 25 µg/mL mycophenolic acid and 50 µg/mL xanthine, followed by subcloning into 96-w ell plates. Single clones were screened for tag expression by immunofluorescence and immunoblot. Endogenous tagging of PP1 A cutting unit specific to the C-t erminus of PP1 (TGGT1_310700; P8) was assembled into the pALH086 HiT vector backbone. Approximately 50  µg of this vector was BsaI- linearized and co-t ransfected with 50 μg of the pSS014 Cas9-e xpression plasmid into TIR1/RH or Parasite popu- lations were selected for 1 week in standard media with 25 µg/mL mycophenolic acid and 50 µg/ mL xanthine and were then subcloned by limiting dilution. Single clones were screened for tag expression by immunofluorescence, immunoblot and sequencing of the junction spanning the 3′ CDS and the tag. Sub-minute phosphoproteomics experiments Parasite harvest and treatment T. gondii tachyzoites from the RH strain were infected onto confluent HFF monolayers in 4–8 15 cm dishes. After the parasites had completely lysed the host cell monolayer (40–48 hr post- infection), extracellular parasites were passed through 5 µm filters into 50 ml conical vials. The samples were spun at 1000 x g for 7 min. The supernatant was decanted, and the parasite pellet was resuspended in 1 ml FluoroBrite DMEM lacking serum and transferred to a 1.5 ml protein low- bind tube. The sample was spun in a mini centrifuge at 1000 x g for 7 min. The supernatant was aspirated, and the parasite pellet was resuspended in 800 µl FluoroBrite DMEM. The sample was split into aliquots of 300 and 500 µl, which were spun at 1000 x g for 7 min followed by aspiration of the supernatant. The pellet containing ⅝ of the parasite harvest was resuspended in 250 µl 500 µM zaprinast in FluoroBrite DMEM. Aliquots of 50 µl were removed and combined with 50 µl 2 X lysis buffer (10% SDS, 4 mM MgCl2, 100 mM TEAB pH 7.55 with 2 X Halt Protease and Phosphatase Inhibitors and 500 U/ml benzonase) at 0, 5, 10, 30, and 60 s post- stimulation. The final composition of the lysate was 5% SDS, 2 mM MgCl2, 50 mM TEAB pH 8.5 with 1 X Halt Protease and Phosphatase Inhibitors and 250 U/ml benzonase. The pellet containing ⅜ of the parasite harvest was resuspended in 150 µl 0.5% DMSO in FluoroBrite DMEM. Aliquots of 50 µl were removed and combined with 50 µl 2 X lysis buffer at 0, 10, and 30 s post- stimulation. Protein cleanup and digestion Proteins were cleaned up and digested using a modified version of the S-t rap protocol (Protifi). Proteins in the lysates were reduced with the addition of 5 mM TCEP and heating at 55 °C for 20 min. Alkylation was performed for 10 min at room temperature with 20 mM MMTS. The lysates were acid- ified to a final concentration of 2.5% v/v phosphoric acid. A 6 X volume of S- trap binding buffer (90% methanol, 100 mM TEAB pH 7.55) was added to each sample to precipitate proteins. The solution was loaded onto S- trap mini columns (Protifi) and spun at 4000 x g until all of the solution had been passed over the column. The columns were washed four times with 150 µl S-t rap binding buffer, followed by a 30 s spin at 4000 x g. Proteins were digested overnight at 37 °C in a humidified incubator with 20 µl of 50 mM TEAB pH 8.5 containing 2 µg of trypsin/LysC mix (Thermo Fisher Scientific). Peptides were eluted in three 40 µl washes with 50 mM TEAB, 0.2% formic acid, and 50% acetonitrile/0.2% formic acid. The eluted peptides snap- frozen and lyophilized. Herneisen et al. eLife 2022;11:e80336. DOI: https://doi.org/10.7554/eLife.80336 25 of 41 Research article Cell Biology | Microbiology and Infectious Disease TMTpro labeling The dried peptides were resuspended in 50 µl 100 mM TEAB 8.5. The peptide concentrations of 1/50 dilutions of the samples were quantified using the Pierce Fluorometric Peptide Assay according to manufacturer’s instructions. Sample abundances were normalized to 50 µg peptides in 50 µl 100 mM TEAB pH 8.5. Each sample was combined with TMTpro reagents at a 5:1 label:peptide weight/weight ratio. Labeling reactions proceeded for 1 hr at room temperature with shaking at 350 rpm. Unreacted TMT reagent was quenched with 0.2% hydroxylamine. The samples were pooled, acidified to 3% with formic acid, and were loaded onto an EasyPep Maxi Sample Prep column. The samples were washed and eluted according to the manufacturer’s instructions. 5% by volume of the elution was reserved as an unenriched proteome sample. The eluted peptides were snap- frozen and lyophilized until dry. Phosphopeptide enrichment Phosphopeptides were enriched using the SMOAC protocol (Tsai et al., 2014). Resuspended TMTpro- labeled samples were enriched with the High- Select TiO2 Phosphopeptide Enrichment Kit (Thermo Fisher Scientific A32993) according to the manufacturer’s instructions. The flow-t hrough and the eluate containing phosphopeptides were immediately snap-f rozen and lyophilized. The flow-t hrough was resuspended and enriched with the High-S elect Fe-N TA Phosphopeptide Enrichment Kit (Thermo Fisher Scientific A32992) according to the manufacturer’s instructions. The eluted phosphopeptides were immediately snap- frozen and lyophilized. Fractionation Unenriched and enriched proteome samples were fractionated with the Pierce High pH Reversed- Phase Peptide Fractionation Kit according to the manufacturer’s instructions for TMT-l abeled peptides. The phosphopeptides enriched with the TiO2 and Fe- NTA kits were combined prior to fractionation. MS data acquisition The fractions were lyophilized and resuspended in 10–20 µl of 0.1% formic acid for MS analysis and were analyzed on an Exploris 480 Orbitrap mass spectrometer equipped with a FAIMS Pro source (Bekker- Jensen et al., 2020) connected to an EASY-n LC chromatography system as described above. Peptides were separated at 300 nl/min on a gradient of 5–20% B for 110 min, 20–28% B for 10 min, 28–95% B for 10 min, 95% B for 10 min, and a seesaw gradient of 95–2% B for 2 min, 2% B for 2 min, 2–98% B for 2 min, 98% B for 2 min, 98–2% B for 2 min, and 2% B for 2 min. The orbitrap and FAIMS were operated in positive ion mode with a positive ion voltage of 1800 V; with an ion transfer tube temperature of 270 °C; using standard FAIMS resolution and compensation voltages of –50 and –65 V, an inner electrode temperature of 100 °C, and outer electrode temperature of 80 °C with 4.6 ml/min carrier gas. Full scan spectra were acquired in profile mode at a resolution of 60,000, with a scan range of 400–1400 m/z, automatically determined maximum fill time, 300% AGC target, intensity threshold of 5×104, 2–5 charge state, and dynamic exclusion of 30 s with a cycle time of 1.5 s between master scans. MS2 spectra were generated with a HCD collision energy of 32 at a resolution of 45,000 using TurboTMT settings with a first mass at 110 m/z, an isolation window of 0.7 m/z, 200% AGC target, and 120 ms injection time. Phosphoproteomics time course analysis Raw files were analyzed in Proteome Discoverer 2.4 (Thermo Fisher Scientific) to generate peak lists and protein and peptide IDs using Sequest HT (Thermo Fisher Scientific) and the ToxoDB release49 GT1 protein database. The unenriched sample search included the following post-t ranslational modifications: dynamic oxidation (+15.995  Da; M), dynamic acetylation (+42.011  Da; N- terminus), static TMTpro (+304.207 Da; any N- terminus), static TMTpro (+304.207 Da; K), and static methylthio (+45.988 Da; C). The enriched sample search included the same post-t ranslational modifications, but with the addition of dynamic phosphorylation (+79.966 Da; S, T, Y). The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE partner repository (Perez- Riverol et al., 2022) with the dataset identifier PXD033765 and 10.6019/PXD033765. Protein and peptide abundances are reported in Supplementary file 1 and Supplementary file 2, respectively. Exported peptide and protein abundance files from Proteome Discoverer 2.4 were loaded into R (version 4.0.4). Summed abundances were calculated for both DMSO and zaprinast-t reated samples, Herneisen et al. eLife 2022;11:e80336. DOI: https://doi.org/10.7554/eLife.80336 26 of 41 Research article Cell Biology | Microbiology and Infectious Disease incorporating the quantification values from both replicates. The ratio of zaprinast and DMSO peptide abundances was transformed into a modified Z score. Peptides with modified Z scores of 3 or above (839 peptides) or –1.8 or below (154 peptides) were considered dynamically changing. The normal- ized abundance values of peptides passing these thresholds were used as input for a principal compo- nent analysis using the R stats package (version 3.6.2). Clustering analysis was performed with the mclust package (Scrucca et al., 2016; version 5.4.7) using the log2 ratios of zaprinast-t reated samples relative to the DMSO-t reated t=0 samples. Cluster assignments are reported in Supplementary file 3. PP1-AID phosphoproteomics Samples were prepared as described above, with the following modifications. PP1- AID parasites were infected onto confluent HFF monolayers. IAA or vehicle (PBS) was added to the cell culture media 29 hr post- infection. At 32 hr post-i nfection, the monolayers were scraped, and parasites were mechan- ically released from host cells by passage through a 27-g auge needle. The parasites were harvested and concentrated as described above. The time-c ourse was initiated when the parasite pellet was resuspended in 250 µl of 500 µM zaprinast solution. Aliquots of 50 µl were removed and combined with 2 X lysis buffer at 0, 10, 30, and 60 s post-s timulation. The experiment was performed in biolog- ical duplicate, with vehicle- and IAA- treated parasites, for a total of 16 samples. The samples were prepared for mass spectrometry, labeled with TMTpro reagents, and enriched as described above. The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE partner repository (Perez- Riverol et al., 2022) with the dataset identifier PXD033765 and 10.6019/PXD033765. Exported peptide and protein abundance files from Proteome Discoverer 2.4 were loaded into R (version 4.0.4). Peptides were clustered if they exhibited 2-f old difference in abundance between vehicle- and IAA- treated samples in both replicates for at least one time point. To focus on peptides dynamically changing during stimulation, we used for each peptide the grouped log2 ratios relative to the t=0 abundance for each treatment as input to the mclust package (Scrucca et al., 2016; version 5.4.7). Thermal profiling temperature range experiments Parasite harvest and treatment Temperature range thermal profiling experiments were performed in biological duplicate on different days. We have optimized thermal profiling for T. gondii and have published detailed descriptions of the protocol (Herneisen and Lourido, 2021). Confluent HFF cells in 15 cm dishes were infected with 2–5 ✕ 107 RH tachyzoites each. When the parasites had fully lysed the host cell monolayer (40–48 hr later), the extracellular parasites passed through a 5 µm filter. The parasite solution was concentrated by centrifugation for 10 min at 1000 x g. Parasites were resuspended in 1 ml of wash buffer (5 mM NaCl, 142 mM KCl, 1 mM MgCl2, 5.6 mM glucose, 25 mM HEPES pH 7.2) and spun again for 10 min at 1000 x g. The parasite pellet was resuspended in 1200 µl of lysis buffer (5 mM NaCl, 142 mM KCl, 1 mM MgCl2, 5.6 mM glucose, 25 mM HEPES pH 7.2 with 0.8% IGEPAL CA-6 30, 1 X Halt Protease Inhibitors, and 250 U/ml benzonase) and subjected to three freeze-t haw cycles. Half of the lysate was combined with an equivalent amount of 2 X EGTA or 10 µM Ca2+ solution. The parasite extracts were incubated at 37 °C/5% CO2 for 5 min. The parasite suspension was aliquoted into PCR tubes and was heated at the following temperatures on two 48-w ell heat blocks: 37, 41, 43, 47, 50, 53, 56, 59, 63, and 67 °C. After 3 min, the tubes were removed from the heat blocks and were chilled on ice for 5 min. The lysates were transferred to a TLA- 100 rotor and were spun at 100,000 x g for 20 min at 4 °C in a Beckman Ultra MAX benchtop ultracentrifuge. The solution, containing the soluble protein fraction, was removed for further processing. Protein cleanup and digestion The concentrations of the 37 °C samples were determined with a DC Protein Assay (BioRad). A volume corresponding to 50 µg of the 37 °C sample was used for further analysis, and equivalent volumes were used from the remaining temperature range samples. Samples were prepared using a modified version of the SP3 protocol (Hughes et al., 2019). A total of 500 µg of a 1:1 mix of hydrophobic and hydrophilic beads (GE Healthcare 65152105050250 and 45152105050250) was added to each Herneisen et al. eLife 2022;11:e80336. DOI: https://doi.org/10.7554/eLife.80336 27 of 41 Research article Cell Biology | Microbiology and Infectious Disease sample, followed by a 6 X volume of 100% ethanol to induce aggregation. The samples were incu- bated at room temperature and 1000 rpm on a thermomixer for 30 min. The beads were magnetically separated, and the supernatant was removed. The beads were washed three times with 80% ethanol followed by magnetic separation. Each sample was resuspended in 88 µl reduction buffer (10 mM TCEP and 50 mM TEAB, pH 8.5) and heated at 55 °C for 1 hr. Samples were alkylated with 25 mM iodoacetamide shielded from light for 1 hr. 2 µg trypsin was added to each sample, and digestion proceeded overnight at 37 °C on a thermomixer shaking at 1000 rpm. The beads were magnetically separated, and the peptide eluatse were removed from further analysis. TMT10plex labeling Peptides were labeled according to the manufacturer’s protocol (Thermo Scientific 90111), with the following modifications. The eluates were combined with 200 µg of TMT10plex reagent in 41 µl of acetonitrile, for an estimated 1:4 w/w peptide:tag labeling reaction. The labeling proceeded for 1 hr at room temperature and was quenched for 15 min with 5% hydroxylamine. The samples were then pooled, flash- frozen, and lyophilized to dryness. The peptides were resuspended in 10% glacial acetic acid and were desalted using a Waters SepPak Light C18 cartridge according to the manufacturer’s instructions. Fractionation Samples were fractionated offline via reversed-p hase high performance liquid chromatography. Samples were applied to a 10 cm ×2.1 mm column packed with 2.6 µm Aeris PEPTIDE XB-C 18 media (Phenomenex) using 20 mM ammonium formate pH 10 in water as Buffer A, 100% acetonitrile as Buffer B, and Shimadzu LC- 20AD pumps. The gradient was isocratic 1% Buffer A for 1 min at 150 µl/ min with increasing Buffer B concentrations to 16.7% B at 20.5 min, 30% B at 31 min and 45% B at 36 min. Fractions were collected with a FRC-1 0A fraction collector and were pooled to eight samples per condition, followed by lyophilization. MS data acquisition The fractions were lyophilized and resuspended in 10–20 µl of 0.1% formic acid for MS analysis and were analyzed on a Q- Exactive HF-X Orbitrap mass spectrometer connected to an EASY- nLC chroma- tography system using 0.1% formic acid as Buffer A and 80% acetonitrile/0.1% formic acid as Buffer B. Peptides were separated at 300 nl/min on a gradient of 6–9% B for 3 min, 9–31% B for 100 min, 31–75% B for 20 min, and 75 to 100% B over 15 min. The orbitrap was operated in positive ion mode. Full scan spectra were acquired in profile mode with a scan range of 375–1400 m/z, resolution of 120,000, maximum fill time of 50ms, and AGC target of 3×106 with a 15  s dynamic exclusion window. Precursors were isolated with a 0.8  m/z window and fragmented with a NCE of 32. The top 20  MS2 spectra were acquired over a scan range of 350–1500 m/z with a resolution of 45,000, AGC target of 8×103, and maximum fill time of 100ms, and first fixed mass of 100 m/z. Thermal profiling temperature range analysis Raw files were analyzed in Proteome Discoverer 2.4 (Thermo Fisher Scientific) to generate peak lists and protein and peptide IDs using Sequest HT (Thermo Fisher Scientific) and the ToxoDB release49 GT1 protein database. The search included the following post-t ranslational modifications: dynamic phosphorylation (+79.966 Da; S, T, Y), dynamic oxidation (+15.995 Da; M), dynamic acetyl- ation (+42.011  Da; N-t erminus), static TMT6plex (+229.163  Da; any N- terminus), static TMT6plex (+229.163 Da; K), and static carbamidomethyl (+57.021 Da; C). Normalization was turned off. The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE partner repository (Perez-R iverol et al., 2022) with the dataset identifier PXD033713 and 10.6019/PXD033713. Protein abundances are reported in Supplementary file 3. Protein abundances were loaded into the R environment (version 4.0.4) and were analyzed using the mineCETSA package (Dziekan et  al., 2020), which performed normalization, generated log- logistic fits of the temperature profiles, and calculated a Euclidean distance score (Supplementary file 3). Individual curve plots were generated using the drc package (Ritz et al., 2015) to fit relative abundances at each temperature or concentration. Herneisen et al. eLife 2022;11:e80336. DOI: https://doi.org/10.7554/eLife.80336 28 of 41 Research article Cell Biology | Microbiology and Infectious Disease Thermal profiling concentration range: Experiment 1 Parasite harvest and treatment Concentration range thermal profiling experiments were performed in biological duplicate on different days. Confluent HFF cells in 15 cm dishes were infected with 2–5×107 RH tachyzoites each. When the parasites had fully lysed the host cell monolayer (40–48 hr later), the extracellular parasites passed through a 5 µm filter. The parasite solution was concentrated by centrifugation for 10 min at 1000 x g. Parasites were resuspended in 1 ml of wash buffer (5 mM NaCl, 142 mM KCl, 1 mM MgCl2, 5.6 mM glucose, 25 mM HEPES pH 7.2) and spun again for 10 min at 1000 x g. The parasite pellet was resuspended in 1200 µl of lysis buffer (5 mM NaCl, 142 mM KCl, 1 mM MgCl2, 5.6 mM glucose, 25 mM HEPES pH 7.2 with 0.8% IGEPAL CA-6 30, 1 X Halt Protease Inhibitors, and 250 U/ml benzonase) and subjected to three freeze-t haw cycles. An equivalent volume of parasite lysate was combined with 2 X [Ca2+]free buffers to attain the final concentrations: 0 nM, 10 nM, 100 nM, 250 nM, 500 nM, 750 nM, 1 µM, 10 µM, 100 µM, and 1 mM. The solutions were aliquoted into two PCR tubes and were incubated for 5 min at 5% CO2/37 °C. The tubes were placed on heat blocks pre- warmed at 54 °C or 58 °C for 3 min and were then immediately placed on ice. The lysates were transferred to a TLA- 100 rotor and were spun at 100,000 x g for 20 min at 4 °C in a Beckman Ultra MAX benchtop ultracentrifuge. The solution, containing the soluble protein fraction, was removed for further processing. MS sample preparation and data acquisition Samples were prepared for mass spectrometry as described in the Temperature Range methods. MS data was acquired using the same instrument and methods as described above. Raw files were searched in Proteome Discoverer 2.4 using the search parameters described in the Temperature Range methods. Normalization was performed in the Proteome Discoverer software. The mass spec- trometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE partner repository (Perez- Riverol et al., 2022) with the dataset identifier PXD033642 and 10.6019/ PXD033642. Protein abundances are reported in Supplementary file 4. Thermal profiling concentration range: Experiment 2 Parasite harvest and treatment Parasites and host cells were grown for one week in SILAC media and were considered separate biological replicates. SILAC media was prepared from DMEM for SILAC (Thermo Fisher 88364), 10% dialyzed FBS (Thermo Fisher A3382001) and L- Arginine HCl/L- Lysine-2 - HCl (Thermo Fisher 89989 and 89987) or 13C 156 N4 L-A rginine HCl/13C 156 N2 L-L ysine-2 HCl (Thermo Fisher 89990 and 88209). The parasites were harvested, lysed, and treated as described in the Concentration range: Experiment 1 procedure. The lysates were incubated on heat blocks pre-w armed at 50 °C, 54 °C, or 58 °C for 3 min and were then immediately placed on ice. Insoluble aggregates were removed by filtration as described in Herneisen and Lourido, 2021. In brief, the lysates were applied to a pre- equilibrated 96- well filter plate (Millipore MSHVN4510) and were spun at 500 x g for 5 min. The concentrations of the filtrates containing soluble fractions were quantified with a DC assay. The heavy and light samples were combined at 1:1 wt/wt, yielding an estimated 50 µg total per concentration. Protein cleanup and digestion Proteins were reduced with 5 mM TCEP for 20 min at 50 °C. Alkylation of cysteines was performed with 15 mM MMTS for 10 min at room temperature. The samples were precipitated in 80% ethanol and were washed using the SP3 protocol as described in the Temperature range procedure. After the final wash, the samples were resuspended in 35 µl of digest buffer (50 mM TEAB) and 1 µg Trypsin. Digestion proceeded overnight at 37 °C on a thermomixer shaking at 1000  rpm. The beads were magnetically separated, and the peptide eluatse were removed from further analysis. TMT10plex labeling Peptides were labeled according to the manufacturer’s protocol (Thermo Scientific 90111), with the following modifications. The eluates were combined with 100 µg of TMT10plex reagent in 15 µl of Herneisen et al. eLife 2022;11:e80336. DOI: https://doi.org/10.7554/eLife.80336 29 of 41 Research article Cell Biology | Microbiology and Infectious Disease acetonitrile, for an estimated 1:2 w/w peptide:tag labeling reaction. The labeling proceeded for 1 hr at room temperature and was quenched for 15 min with 5% hydroxylamine. The samples were then pooled, flash- frozen, and lyophilized to dryness. Fractionation Samples were fractionated with the Pierce High pH Reversed- Phase Peptide Fractionation Kit according to the manufacturer’s instructions for TMT- labeled peptides. MS data acquisition The fractions were lyophilized and resuspended in 10–20 µl of 0.1% formic acid for MS analysis and were analyzed on an Exploris 480 Orbitrap mass spectrometer equipped with a FAIMS Pro source (Bekker- Jensen et al., 2020) connected to an EASY-n LC chromatography system as described above. Peptides were separated at 300 nl/min on a gradient of 6–21% B for 41 min, 21–36% B for 20 min, 36–50% B for 10 min, and 50 to 100% B over 15 min. The orbitrap and FAIMS were operated in posi- tive ion mode with a positive ion voltage of 1800 V; with an ion transfer tube temperature of 270 °C; using standard FAIMS resolution and compensation voltages of –50 and –65 V (injection 1) or –40 and –60 (injection 2). Full scan spectra were acquired in profile mode at a resolution of 120,000, with a scan range of 350–1200 m/z, automatically determined maximum fill time, standard AGC target, intensity threshold of 5×103, 2–5 charge state, and dynamic exclusion of 30 s with a cycle time of 2 s between master scans. MS2 spectra were generated with a HCD collision energy of 36 at a resolution of 30,000 using TurboTMT settings with a first mass at 110 m/z, an isolation window of 0.7 m/z, standard AGC target, and auto injection time. Thermal profiling concentration range analysis Raw files were analyzed in Proteome Discoverer 2.4 (Thermo Fisher Scientific) to generate peak lists and protein and peptide IDs using Sequest HT (Thermo Fisher Scientific) and the ToxoDB release49 GT1 protein database. The search included the following post- translational modifications for the first biological replicate, which was collected from parasites grown in light SILAC medium: dynamic phosphorylation (+79.966 Da; S, T, Y), dynamic oxidation (+15.995 Da; M), dynamic acetyl- ation (+42.011  Da; N-t erminus), static TMT6plex (+229.163  Da; any N-t erminus), static TMT6plex (+229.163 Da; K), and static methylthio (+45.988 Da; C). The heavy samples were searched for the additional modifications Lys8- TMT6plex (+237.177 Da; K) and static Label:13C 156 N4 (+10.008 Da; R). The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE partner repository (Perez- Riverol et al., 2022) with the dataset identifier PXD033650 and 10.6019/PXD033650. Protein abundances are reported in Supplementary file 4. Protein abundances from Proteome Discoverer 2.4 were loaded into the R environment (version 4.0.4) and were analyzed using the mineCETSA package (Dziekan et al., 2020). The data were not further normalized using the package, as normalization had already been performed by the Proteome Discoverer software. Log- logistic fitting of the relative abundance profiles was performed using default settings. The AUC reflects the mean stability change when proteins were detected in both replicates of an experiment. If proteins were not detected, the AUC was calculated from one replicate. Proteins were considered Ca2+-r esponsive if the curve- fit parameters had an R2 >0.8 and AUC two modified Z scores from the median. The two concentration range experiments were analyzed separately. Indi- vidual curve plots were generated using the drc package (Ritz et al., 2015) to fit relative abundances at each temperature or concentration. Enrichment analysis Sets of gene ontology terms from the differentially regulated (or Ca2+- responsive) and background proteome (all proteins with quantification values in each mass spectrometry experiment) were down- loaded from ToxoDB.org (Molecular Function, Computed evidence, P-V alue cutoff set to 1). Gene ontology terms were tested for enrichment across all gene ontology terms identified in the back- ground proteome. A p value for the likelihood of a given enrichment to have occurred by chance was obtained using a hypergeometric test. Herneisen et al. eLife 2022;11:e80336. DOI: https://doi.org/10.7554/eLife.80336 30 of 41 Research article Cell Biology | Microbiology and Infectious Disease Immunoblotting Samples were prepared as described in the thermal profile concentration range or phosphopro- teomics procedures prior to proteomics sample preparation. The samples, which had already been treated with benzonase, were combined with 5 X laemmli sample buffer (10% SDS, 50% glycerol, 300 mM Tris HCl pH 6.8, 0.05% bromophenol blue, 5% beta- mercaptoethanol) and were incubated at 37 °C for 10 min. The samples were then run on precast 4–15% SDS gels (BioRad) and were trans- ferred overnight onto nitrocellulose membranes at 4 °C and 30 mA in 25 mM TrisHCl, 192 mM glycine, and 20% methanol. Blocking and antibody incubations were performed in 5% milk in TBS- T for 1 hr at room temperature. The membrane was washed three times with TBS-T between antibody incuba- tions. Imaging was performed with the LICOR Odyssey CLx. Immunofluorescence assays Confluent HFFs seeded onto coverslips were infected with extracellular parasites and were grown at 37 °C/5% CO2. Approximately 21 hr later, IAA or a vehicle solution of PBS were added to the wells to a final concentration of 500 µM where indicated. At 24 hr post-i nfection, the media was aspirated, and coverslips were fixed in 4% formaldehyde in PBS. Following three washes in PBS, the fixed cells were permeabilized with 0.25% triton for 10  min at room temperature. Residual permeabilization solution was removed with three washes of PBS. The coverslips were incubated in blocking solution (5% IFS/5% NGS in PBS) for 10 min at room temperature, followed by a 60-m in incubation in primary antibody solution. An anti- CDPK1 antibody (Covance) was used as a parasite counterstain (Waldman et al., 2020). After three washes with PBS, the coverslips were incubated in blocking solution at room temperature for 5 min, followed by a 30- min incubation in secondary antibody solution. The coverslips were washed three times in PBS and once in water. Coverslips were mounted with Prolong Diamond and were set for 30 min at 37 °C. Imaging was performed with the Nikon Ti Eclipse and NIS Elements software package. Invasion assays Confluent HFFs seeded onto coverslips were incubated with 5× 106 extracellular parasites for 60 min at 37 °C/5% CO2. The coverslips were washed four six times with PBS and were fixed for 10 min at room temperature with 4% formaldehyde in PBS. The coverslips were incubated in blocking solution (1% BSA in PBS) for 10 min. Extracellular parasites were stained with mouse anti-S AG1 for 30 min at room temperature. Following permeabilization with 0.25% triton- X100 for 10 min, the coverslips were incubated with guinea- pig anti- CDPK1 as a parasite counterstain for 30 min at room temperature. The coverslips were incubated with a secondary antibody solution containing Hoechst and were mounted on coverglass with Prolong Diamond. The number of parasites invaded was calculated by normalizing the number of intracellular, invaded parasites to host cell nuclei in a field of view. Five random fields of view were imaged per coverslip. Each experiment was performed in technical duplicate. Egress assays Automated, plate- based egress assays were performed as previously described (Shortt and Lourido, 2020). In brief, HFF monolayers in a clear-b ottomed 96- well plate were infected with 7.5× 104 or 1× 105 parasites of the TIR1 or PP1-A ID strains, respectively. IAA or PBS were added to a final concen- tration of 500 µM 20 hr later. After 3 hr, the media was exchanged for FluoroBrite supplemented with 3% calf serum. Three images were taken before zaprinast (final concentration 500 μM) or A23187 (final concentration 8 µM) and DAPI (final concentration 5 ng /mL) were injected. Imaging of DAPI-s tained host cell nuclei continued for 9 additional minutes before 1% Triton X- 100 was injected into all wells to determine the total number of host cell nuclei. Imaging was performed at 37 °C and 5% CO2 using a Biotek Cytation 3. Results are the mean of two wells per condition and are representative of three independent experiments. Replication assays Parasites were inoculated onto coverslips containing HFFs. After 1 hr, the media was aspirated and replaced with media containing 500 µM IAA or PBS vehicle. At 24 hr post- IAA addition, the cover- slips with intracellular parasites were fixed, permeabilized, and stained with CDPK1 antibody and Hoechst as described under ‘Immunofluorescence assays’. For each sample, multiple fields of view Herneisen et al. eLife 2022;11:e80336. DOI: https://doi.org/10.7554/eLife.80336 31 of 41 Research article Cell Biology | Microbiology and Infectious Disease were acquired with an Eclipse Ti microscope (Nikon) and the number of nuclei per vacuole were calcu- lated from the full field of view (at least 100 vacuoles). Results are the mean of three independent experiments. Plaque assays A total of 500 parasites were inoculated into 12-w ell plates of HFFs maintained in D10 and allowed to grow undisturbed for 7 days. IAA or vehicle (PBS) was added to a final concentration of 100 µM. Plates were washed with PBS and fixed for 10 min at room temperature with 100% ethanol. Staining was performed for 5 min at room temperature with crystal violet solution, followed by two washes with PBS, one wash with water, and drying. Live microscopy PP1- mNG parasites were grown in HFFs in glass- bottom 35 mm dishes (Ibidi) for 24 hours. The media was decanted and the dish was washed once with 1 ml Ringer’s buffer (155 mM NaCl, 2 mM CaCl2, 3 mM KCl, 1 mM MgCl2, 3 mM NaH2PO4, 10 mM HEPES, 10 mM glucose). Parasites were stimulated to egress with 500 μM zaprinast or 4 µM A23187 in Ringer’s buffer (155 mM NaCl, 2 mM CaCl2, 3 mM KCl, 1 mM MgCl2, 3 mM NaH2PO4, 10 mM HEPES, 10 mM glucose) supplemented with 1% FBS (v/v) and recorded every 2 s for 300 s using an Eclipse Ti microscope (Nikon) with an enclosure maintained at 37 °C. Cytosolic Ca2+ measurements with FURA-2AM Fura-2 AM loading of T. gondii tachyzoites was done as described previously (Moreno and Zhong, 1996). Briefly, freshly collected parasites (TIR1 or PP1-A ID parasites treated with IAA for 5 hr) were washed twice with buffer A plus glucose or BAG (116 mM NaCL, 5.4 mM KCL, 0.8 mM MgSO4.7H2O, 50 mM Hepes pH 7.3, 5 mM Glucose) by centrifugation (706 x g for 10 min) and re- suspended to a final density of 1 x l09 parasites/ml in loading buffer (BAG plus 1.5% sucrose, and 5 μM Fura-2 AM). The suspension was incubated for 26 min at 26 °C with mild agitation. Subsequently, parasites were washed twice by centrifugation (2000 x g for 2 min) with BAG to remove extracellular dye, re-s us- pended to a final density of 1x109 parasites per ml in BAG and kept on ice. Parasites are viable for a few hours under these conditions. For fluorescence measurements, 2x107 parasites/mL were placed in a cuvette with 2.5 mL of BAG. Fluorescence measurements were done in a Hitachi F-7 100 fluorescence spectrometer using the Fura 2 conditions for excitation (340 and 380 nm) and emission (510  nm). The Fura- 2 fluorescence response to Ca2+ was calibrated from the ratio of 340/380  nm fluorescence values after subtraction of background fluorescence at 340 and 380 nm as described previously (Grynkiewicz et al., 1985). Ca2+ release rate is the change in Ca2+ concentration during the initial 20 s after compound addition. Delta [Ca2+] was calculated by the difference between the higher Ca2+ peak and basal Ca2+. Acknowledgements We thank Emily Shortt for assistance with cell culture, Eric Spooner of the Whitehead Proteomics Core facility for assistance with sample preparation, and Tyler Smith for providing assistance with tagging vector design. We thank L David Sibley for the MIC2 and SAG1 antibodies, Dominique Soldati-F avre for the GAP45 antibody, and Marc-J an Gubbels for the TUB1 antibody. This research was supported by funds from National Institutes of Health grants to SL (R01AI144369) and SNJM (R01AI128356 and R21AI15493), and a National Science Foundation Graduate Research Fellowship to ALH (174530). Additional information Competing interests Sebastian Lourido: Reviewing editor, eLife. The other authors declare that no competing interests exist. Herneisen et al. eLife 2022;11:e80336. DOI: https://doi.org/10.7554/eLife.80336 32 of 41 Research article Cell Biology | Microbiology and Infectious Disease Funding Funder Grant reference number Author National Institutes of R01AI144369 Sebastian Lourido Health National Science 174530 Alice L Herneisen Foundation National Institutes of R01AI128356 Silvia NJ Moreno Health National Institutes of R21AI15493 Silvia NJ Moreno Health The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication. Author contributions Alice L Herneisen, Conceptualization, Data curation, Formal analysis, Investigation, Visualization, Writing - original draft, Writing - review and editing; Zhu- Hong Li, Investigation; Alex W Chan, Meth- odology; Silvia NJ Moreno, Resources, Supervision, Writing - review and editing; Sebastian Lourido, Conceptualization, Resources, Supervision, Funding acquisition, Writing - review and editing Author ORCIDs Alice L Herneisen http://orcid.org/0000-0003-3368-0893 Silvia NJ Moreno http://orcid.org/0000-0002-2041-6295 Sebastian Lourido http://orcid.org/0000-0002-5237-1095 Decision letter and Author response Decision letter https://doi.org/10.7554/eLife.80336.sa1 Author response https://doi.org/10.7554/eLife.80336.sa2 Additional files Supplementary files • Supplementary file 1. Sub-m inute phosphoproteomics time course protein and abundance assignments from Proteome Discoverer 2.4. • Supplementary file 2. Sub-m inute phosphoproteomics time course peptide and abundance assignments from Proteome Discoverer 2.4. Mclust cluster assignments (column 118) of phosphopeptides dynamically changing during zaprinast treatment. • Supplementary file 3. Data pertaining to the temperature range thermal profiling experiment. 1. Protein and abundance assignments from Proteome Discoverer 2.4 for samples with 0 µM Ca2+, replicate 1. 2. Protein and abundance assignments from Proteome Discoverer 2.4 for samples with 0 µM Ca2+, replicate 2. 3. Protein and abundance assignments from Proteome Discoverer 2.4 for samples with 10 µM Ca2+, replicate 1. 4. Protein and abundance assignments from Proteome Discoverer 2.4 for samples with 10 µM Ca2+, replicate 2. 5. Curve fit output from the mineCETSA package. 6. Area under the euclidean distance score calculations from the mineCETSA package. • Supplementary file 4. Data pertaining to the concentration range thermal profiling experiments. 1. Protein and abundance assignments from Proteome Discoverer 2.4 for Experiment 1 samples with 54 °C, replicate 1. 2. Protein and abundance assignments from Proteome Discoverer 2.4 for Experiment 1 samples with 54 °C, replicate 2. 3. Protein and abundance assignments from Proteome Discoverer 2.4 for Experiment 1 samples with 58 °C, replicate 1. 4. Protein and abundance assignments from Proteome Discoverer 2.4 for Experiment 1 samples with 58 °C, replicate 2. 5. Curve fit output for concentration range Experiment 1 from the mineCETSA package. 6. Area under the curve score calculations from the mineCETSA package for concentration range Experiment 2. 7. Protein and abundance assignments from Proteome Discoverer 2.4 for Experiment 2 samples with 50 °C, replicate 1. 8. Protein and abundance assignments from Proteome Discoverer 2.4 for Experiment 2 samples with 50 °C, replicate 2. 8. Protein and abundance assignments from Proteome Discoverer 2.4 for Experiment 2 samples with 54 °C, replicate 1. 9. Protein and abundance assignments from Proteome Discoverer 2.4 for Experiment 2 samples with 54 °C, replicate 2. 10. Protein and abundance assignments from Proteome Discoverer 2.4 for Experiment 2 samples Herneisen et al. eLife 2022;11:e80336. DOI: https://doi.org/10.7554/eLife.80336 33 of 41 Research article Cell Biology | Microbiology and Infectious Disease with 58 °C, replicate 1. 11. Protein and abundance assignments from Proteome Discoverer 2.4 for Experiment 2 samples with 58 °C, replicate 2. 12. Curve fit output for concentration range Experiment 2 from the mineCETSA package. 13. Area under the curve score calculations from the mineCETSA package for concentration range Experiment 2. • Supplementary file 5. PP1 depletion zaprinast phosphoproteomics time course protein and abundance assignments from Proteome Discoverer 2.4. • Supplementary file 6. PP1 depletion zaprinast phosphoproteomics time course peptide and abundance assignments from Proteome Discoverer 2.4. Mclust cluster assignments (column 2) of phosphopeptides dynamically changing during zaprinast treatment when PP1 is depleted. • Supplementary file 7. Sequences and accessions of oligonucleotides and plasmids used in this study. • MDAR checklist Data availability All mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE partner repository with the dataset identifier PXD033765 and https://doi.org/10.6019/ PXD033765. All other information is provided in the supplementary files. The following dataset was generated: Author(s) Year Dataset title Dataset URL Database and Identifier Herneisen AL 2022 De novo mapping https://www. ebi.a c. PRIDE, PXD033765 of the apicomplexan uk/ pride/ archive/ Ca2+-responsive projects/ PXD033765 phosphoproteome References Alam MM, Solyakov L, Bottrill AR, Flueck C, Siddiqui FA, Singh S, Mistry S, Viskaduraki M, Lee K, Hopp CS, Chitnis CE, Doerig C, Moon RW, Green JL, Holder AA, Baker DA, Tobin AB. 2015. Phosphoproteomics reveals malaria parasite protein kinase G as a signalling hub regulating egress and invasion. Nature Communications 6:1–15. DOI: https://doi.org/10.1038/ncomms8285, PMID: 26149123 Alves E, Benns HJ, Magnus L, Dominicus C, Dobai T, Blight J, Wincott CJ, Child MA. 2021. An extracellular redox signal triggers calcium release and impacts the asexual development of Toxoplasma gondii Frontiers in Cellular and Infection Microbiology 11:728425. DOI: https://doi.org/10.3389/fcimb.2021.728425, PMID: 34447699 Arrizabalaga G, Ruiz F, Moreno S, Boothroyd JC. 2004. Ionophore- resistant mutant of Toxoplasma gondii reveals involvement of a sodium/hydrogen exchanger in calcium regulation. The Journal of Cell Biology 165:653–662. DOI: https://doi.org/10.1083/jcb.200309097, PMID: 15173192 Back PS, O’Shaughnessy WJ, Moon AS, Dewangan PS, Hu X, Sha J, Wohlschlegel JA, Bradley PJ, Reese ML. 2020. Ancient MAPK ERK7 is regulated by an unusual inhibitory scaffold required for toxoplasma apical complex biogenesis. PNAS 117:12164–12173. DOI: https://doi.org/10.1073/pnas.1921245117, PMID: 32409604 Balestra AC, Koussis K, Klages N, Howell SA, Flynn HR, Bantscheff M, Pasquarello C, Perrin AJ, Brusini L, Arboit P, Sanz O, Castaño LPB, Withers- Martinez C, Hainard A, Ghidelli- Disse S, Snijders AP, Baker DA, Blackman MJ, Brochet M. 2021. Ca2+ signals critical for egress and gametogenesis in malaria parasites depend on a multipass membrane protein that interacts with PKG. Science Advances 7:eabe5396. DOI: https://doi.org/ 10.1126/sciadv.abe5396, PMID: 33762339 Bansal P, Antil N, Kumar M, Yamaryo-B otté Y, Rawat RS, Pinto S, Datta KK, Katris NJ, Botté CY, Prasad TSK, Sharma P. 2021. Protein kinase tgcdpk7 regulates vesicular trafficking and phospholipid synthesis in Toxoplasma gondii. PLOS Pathogens 17:e1009325. DOI: https://doi.org/10.1371/journal.ppat.1009325, PMID: 33635921 Barylyuk K, Koreny L, Ke H, Butterworth S, Crook OM, Lassadi I, Gupta V, Tromer E, Mourier T, Stevens TJ, Breckels LM, Pain A, Lilley KS, Waller RF. 2020. A comprehensive subcellular atlas of the toxoplasma proteome via hyperlopit provides spatial context for protein functions. Cell Host & Microbe 28:752–766.. DOI: https:// doi.org/10.1016/j.chom.2020.09.011, PMID: 33053376 Bekker- Jensen DB, Martínez- Val A, Steigerwald S, Rüther P, Fort KL, Arrey TN, Harder A, Makarov A, Olsen JV. 2020. A compact quadrupole- orbitrap mass spectrometer with FAIMS interface improves proteome coverage in short LC gradients. Molecular & Cellular Proteomics 19:716–729. DOI: https://doi.org/10.1074/mcp.TIR119. 001906, PMID: 32051234 Berridge MJ, Lipp P, Bootman MD. 2000. The versatility and universality of calcium signalling. Nature Reviews. Molecular Cell Biology 1:11–21. DOI: https://doi.org/10.1038/35036035, PMID: 11413485 Herneisen et al. eLife 2022;11:e80336. DOI: https://doi.org/10.7554/eLife.80336 34 of 41 Research article Cell Biology | Microbiology and Infectious Disease Bertolotti A. 2018. The split protein phosphatase system. The Biochemical Journal 475:3707–3723. DOI: https:// doi.org/10.1042/BCJ20170726, PMID: 30523060 Billker O, Lourido S, Sibley LD. 2009. Calcium- dependent signaling and kinases in apicomplexan parasites. Cell Host & Microbe 5:612–622. DOI: https://doi.org/10.1016/j.chom.2009.05.017, PMID: 19527888 Birnbaum J, Scharf S, Schmidt S, Jonscher E, Hoeijmakers WAM, Flemming S, Toenhake CG, Schmitt M, Sabitzki R, Bergmann B, Fröhlke U, Mesén- Ramírez P, Blancke Soares A, Herrmann H, Bártfai R, Spielmann T. 2020. A kelch13- defined endocytosis pathway mediates artemisinin resistance in malaria parasites. Science 367:51–59. DOI: https://doi.org/10.1126/science.aax4735, PMID: 31896710 Bisio H, Lunghi M, Brochet M, Soldati-F avre D. 2019. Phosphatidic acid governs natural egress in Toxoplasma gondii via a guanylate cyclase receptor platform. Nature Microbiology 4:420–428. DOI: https://doi.org/10. 1038/s41564-018-0339-8, PMID: 30742070 Bisio Hugo, Soldati- Favre D. 2019. Signaling cascades governing entry into and exit from host cells by Toxoplasma gondii Annual Review of Microbiology 73:579–599. DOI: https://doi.org/10.1146/annurev-micro- 020518-120235, PMID: 31500539 Blader IJ, Coleman BI, Chen CT, Gubbels MJ. 2015. Lytic cycle of Toxoplasma gondii: 15 years later. Annual Review of Microbiology 69:463–485. DOI: https://doi.org/10.1146/annurev-micro-091014-104100, PMID: 26332089 Brautigan DL, Shenolikar S. 2018. Protein serine/threonine phosphatases: keys to unlocking regulators and substrates. Annual Review of Biochemistry 87:921–964. DOI: https://doi.org/10.1146/annurev-biochem- 062917-012332, PMID: 29925267 Brochet M, Collins MO, Smith TK, Thompson E, Sebastian S, Volkmann K, Schwach F, Chappell L, Gomes AR, Berriman M, Rayner JC, Baker DA, Choudhary J, Billker O. 2014. Phosphoinositide metabolism links cgmp- dependent protein kinase G to essential ca2+. PLOS Biology 12:e1001806. DOI: https://doi.org/10.1371/ journal.pbio.1001806, PMID: 24594931 Brown KM, Lourido S, Sibley LD. 2016. Serum albumin stimulates protein kinase G- dependent microneme secretion in Toxoplasma gondii. The Journal of Biological Chemistry 291:9554–9565. DOI: https://doi.org/10. 1074/jbc.M115.700518, PMID: 26933037 Brown KM, Long S, Sibley LD. 2017. Plasma membrane association by N-a cylation governs PKG function in Toxoplasma gondii. MBio 8:e00375-1 7. DOI: https://doi.org/10.1128/mBio.00375-17, PMID: 28465425 Brown KM, Sibley LD. 2018. Essential cgmp signaling in toxoplasma is initiated by a hybrid P-t ype atpase- guanylate cyclase. Cell Host & Microbe 24:804–816.. DOI: https://doi.org/10.1016/j.chom.2018.10.015, PMID: 30449726 Brown KM, Tonkin CJ, Billker O, Sibley LD. 2020. Calcium and cyclic nucleotide signaling networks in Toxoplasma gondii. Weiss LM, Kim K (Eds). Toxoplasma gondii. Third Edition. Academic Press. p. 577–605. DOI: https://doi.org/10.1016/B978-0-12-815041-2.00013-X Bullen HE, Jia Y, Yamaryo- Botté Y, Bisio H, Zhang O, Jemelin NK, Marq JB, Carruthers V, Botté CY, Soldati- Favre D. 2016. Phosphatidic acid-m ediated signaling regulates microneme secretion in toxoplasma. Cell Host & Microbe 19:349–360. DOI: https://doi.org/10.1016/j.chom.2016.02.006, PMID: 26962945 Bullen HE, Bisio H, Soldati-F avre D. 2019. The triumvirate of signaling molecules controlling toxoplasma microneme exocytosis: cyclic GMP, calcium, and phosphatidic acid. PLOS Pathogens 15:e1007670. DOI: https://doi.org/10.1371/journal.ppat.1007670, PMID: 31121005 Chang L, Dykes EJ, Li J, Moreno SNJ, Hortua Triana MA. 2019. Characterization of two EF- hand domain- containing proteins from Toxoplasma gondii. The Journal of Eukaryotic Microbiology 66:343–353. DOI: https:// doi.org/10.1111/jeu.12675, PMID: 30063275 Charoenwattanasatien R, Zinzius K, Scholz M, Wicke S, Tanaka H, Brandenburg JS, Marchetti GM, Ikegami T, Matsumoto T, Oda T, Sato M, Hippler M, Kurisu G. 2020. Calcium sensing via EF- hand 4 enables thioredoxin activity in the sensor- responder protein calredoxin in the green alga Chlamydomonas reinhardtii The Journal of Biological Chemistry 295:170–180. DOI: https://doi.org/10.1074/jbc.RA119.008735, PMID: 31776187 Chen T-W , Wardill TJ, Sun Y, Pulver SR, Renninger SL, Baohan A, Schreiter ER, Kerr RA, Orger MB, Jayaraman V, Looger LL, Svoboda K, Kim DS. 2013. Ultrasensitive fluorescent proteins for imaging neuronal activity. Nature 499:295–300. DOI: https://doi.org/10.1038/nature12354, PMID: 23868258 Chen AL, Kim EW, Toh JY, Vashisht AA, Rashoff AQ, Van C, Huang AS, Moon AS, Bell HN, Bentolila LA, Wohlschlegel JA, Bradley PJ. 2015. Novel components of the toxoplasma inner membrane complex revealed by bioid. MBio 6:e02357- 14. DOI: https://doi.org/10.1128/mBio.02357-14, PMID: 25691595 Chen K, Günay- Esiyok Ö, Klingeberg M, Marquardt S, Pomorski TG, Gupta N. 2021. Aminoglycerophospholipid flipping and P4- atpases in Toxoplasma gondii. The Journal of Biological Chemistry 296:100315. DOI: https:// doi.org/10.1016/j.jbc.2021.100315, PMID: 33485966 Chern JH, Pasquarelli RR, Moon AS, Chen AL, Sha J, Wohlschlegel JA, Bradley PJ. 2021. A novel toxoplasma inner membrane complex suture- associated protein regulates suture protein targeting and colocalizes with membrane trafficking machinery. MBio 12:e0245521. DOI: https://doi.org/10.1128/mBio.02455-21, PMID: 34634933 Clapham DE. 2007. Calcium signaling. Cell 131:1047–1058. DOI: https://doi.org/10.1016/j.cell.2007.11.028, PMID: 18083096 Dai L, Prabhu N, Yu LY, Bacanu S, Ramos AD, Nordlund P. 2019. Horizontal cell biology: monitoring global changes of protein interaction states with the proteome- wide cellular thermal shift assay (CETSA). Annual Review of Biochemistry 88:383–408. DOI: https://doi.org/10.1146/annurev-biochem-062917-012837, PMID: 30939043 Herneisen et al. eLife 2022;11:e80336. DOI: https://doi.org/10.7554/eLife.80336 35 of 41 Research article Cell Biology | Microbiology and Infectious Disease Delorme V, Garcia A, Cayla X, Tardieux I. 2002. A role for Toxoplasma gondii type 1 ser/thr protein phosphatase in host cell invasion. Microbes and Infection 4:271–278. DOI: https://doi.org/10.1016/s1286-4579(02)01538-1, PMID: 11909736 Dominicus C, Nofal SD, Broncel M, Katris NJ, Flynn HR, Arrizabalaga G, Botté CY, Invergo BM, Treeck M. 2021. A positive feedback loop mediates crosstalk between calcium, cyclic nucleotide and lipid signalling in Toxoplasma gondii. bioRxiv. DOI: https://doi.org/10.1101/2021.11.29.470317 Dos Santos Pacheco N, Tosetti N, Krishnan A, Haase R, Maco B, Suarez C, Ren B, Soldati- Favre D. 2021. Revisiting the role of Toxoplasma gondii ERK7 in the maintenance and stability of the apical complex. MBio 12:e0205721. DOI: https://doi.org/10.1128/mBio.02057-21, PMID: 34607461 Dziekan JM, Wirjanata G, Dai L, Go KD, Yu H, Lim YT, Chen L, Wang LC, Puspita B, Prabhu N, Sobota RM, Nordlund P, Bozdech Z. 2020. Cellular thermal shift assay for the identification of drug-t arget interactions in the Plasmodium falciparum proteome. Nature Protocols 15:1881–1921. DOI: https://doi.org/10.1038/s41596-020- 0310-z, PMID: 32341577 Engelberg K, Chen CT, Bechtel T, Sánchez Guzmán V, Drozda AA, Chavan S, Weerapana E, Gubbels MJ. 2020. The apical annuli of Toxoplasma gondii are composed of coiled-c oil and signalling proteins embedded in the inner membrane complex sutures. Cellular Microbiology 22:e13112. DOI: https://doi.org/10.1111/cmi.13112, PMID: 31470470 Erler H, Ren B, Gupta N, Beitz E. 2018. The intracellular parasite Toxoplasma gondii harbors three druggable FNT- type formate and l- lactate transporters in the plasma membrane. The Journal of Biological Chemistry 293:17622–17630. DOI: https://doi.org/10.1074/jbc.RA118.003801, PMID: 30237165 Fang J, Marchesini N, Moreno SNJ. 2006. A Toxoplasma gondii phosphoinositide phospholipase C (tgpi-P LC) with high affinity for phosphatidylinositol. The Biochemical Journal 394:417–425. DOI: https://doi.org/10.1042/ BJ20051393, PMID: 16288600 Farrell A, Thirugnanam S, Lorestani A, Dvorin JD, Eidell KP, Ferguson DJP, Anderson- White BR, Duraisingh MT, Marth GT, Gubbels MJ. 2012. A DOC2 protein identified by mutational profiling is essential for apicomplexan parasite exocytosis. Science 335:218–221. DOI: https://doi.org/10.1126/science.1210829, PMID: 22246776 Feng B, Stemmer PM. 1999. Interactions of calcineurin A, calcineurin B, and ca2+. Biochemistry 38:12481– 12489. DOI: https://doi.org/10.1021/bi990492w, PMID: 10493818 Feng B, Stemmer PM. 2001. Ca2+ binding site 2 in calcineurin-B modulates calmodulin-d ependent calcineurin phosphatase activity. Biochemistry 40:8808–8814. DOI: https://doi.org/10.1021/bi0025161, PMID: 11467941 Francia ME, Wicher S, Pace DA, Sullivan J, Moreno SNJ, Arrizabalaga G. 2011. A Toxoplasma gondii protein with homology to intracellular type na+/H+ exchangers is important for osmoregulation and invasion. Experimental Cell Research 317:1382–1396. DOI: https://doi.org/10.1016/j.yexcr.2011.03.020, PMID: 21501607 Frénal K, Marq JB, Jacot D, Polonais V, Soldati- Favre D, Frischknecht F. 2014. Plasticity between myoc- and myoa- glideosomes: an example of functional compensation in Toxoplasma gondii invasion. PLOS Pathogens 10:e1004504. DOI: https://doi.org/10.1371/journal.ppat.1004504, PMID: 25393004 Frénal K, Jacot D, Hammoudi PM, Graindorge A, Maco B, Soldati- Favre D. 2017. Myosin- dependent cell-c ell communication controls synchronicity of division in acute and chronic stages of Toxoplasma gondii. Nature Communications 8:15710. DOI: https://doi.org/10.1038/ncomms15710, PMID: 28593938 Garcia CRS, Alves E, Pereira PHS, Bartlett PJ, Thomas AP, Mikoshiba K, Plattner H, Sibley LD. 2017. InsP3 signaling in apicomplexan parasites. Current Topics in Medicinal Chemistry 17:2158–2165. DOI: https://doi.o rg /10.2174/1568026617666170130121042, PMID: 28137231 Gaskins E, Gilk S, DeVore N, Mann T, Ward G, Beckers C. 2004. Identification of the membrane receptor of a class XIV myosin in Toxoplasma gondii. The Journal of Cell Biology 165:383–393. DOI: https://doi.org/10.1083/ jcb.200311137, PMID: 15123738 Gazarini ML, Garcia CRS. 2004. The malaria parasite mitochondrion senses cytosolic ca2+ fluctuations. Biochemical and Biophysical Research Communications 321:138–144. DOI: https://doi.org/10.1016/j.bbrc. 2004.06.141, PMID: 15358226 Giorgi C, Marchi S, Pinton P. 2018. The machineries, regulation and cellular functions of mitochondrial calcium. Nature Reviews. Molecular Cell Biology 19:713–730. DOI: https://doi.org/10.1038/s41580-018-0052-8, PMID: 30143745 Graindorge A, Frénal K, Jacot D, Salamun J, Marq JB, Soldati- Favre D. 2016. The conoid associated motor myoh is indispensable for Toxoplasma gondii entry and exit from host cells. PLOS Pathogens 12:e1005388. DOI: https://doi.org/10.1371/journal.ppat.1005388, PMID: 26760042 Grynkiewicz G, Poenie M, Tsien RY. 1985. A new generation of ca2+ indicators with greatly improved fluorescence properties. The Journal of Biological Chemistry 260:3440–3450. DOI: https://doi.org/10.1016/ S0021-9258(19)83641-4, PMID: 3838314 Guttery DS, Pittman JK, Frénal K, Poulin B, McFarlane LR, Slavic K, Wheatley SP, Soldati- Favre D, Krishna S, Tewari R, Staines HM. 2013. The plasmodium berghei ca(2+)/H(+) exchanger, pbcax, is essential for tolerance to environmental ca(2+) during sexual development. PLOS Pathogens 9:e1003191. DOI: https://doi.org/10. 1371/journal.ppat.1003191 Harding CR, Gow M, Kang JH, Shortt E, Manalis SR, Meissner M, Lourido S. 2019. Alveolar proteins stabilize cortical microtubules in Toxoplasma gondii. Nature Communications 10:401. DOI: https://doi.org/10.1038/ s41467-019-08318-7, PMID: 30674885 Harding CR, Sidik SM, Petrova B, Gnädig NF, Okombo J, Herneisen AL, Ward KE, Markus BM, Boydston EA, Fidock DA, Lourido S. 2020. Genetic screens reveal a central role for heme metabolism in artemisinin Herneisen et al. eLife 2022;11:e80336. DOI: https://doi.org/10.7554/eLife.80336 36 of 41 Research article Cell Biology | Microbiology and Infectious Disease susceptibility. Nature Communications 11:4813. DOI: https://doi.org/10.1038/s41467-020-18624-0, PMID: 32968076 Havelaar AH, Kirk MD, Torgerson PR, Gibb HJ, Hald T, Lake RJ, Praet N, Bellinger DC, Silva NR, Gargouri N, Speybroeck N, Cawthorne A, Mathers C, Stein C, Angulo FJ, Devleesschauwer B. 2015. World health organization global estimates and regional comparisons of the burden of foodborne disease in 2010. PLOS Medicine 12:1001923. DOI: https://doi.org/10.1371/journal.pmed.1001923 Hayward JA, Rajendran E, Zwahlen SM, Faou P, van Dooren GG, Soldati- Favre D. 2021. Divergent features of the coenzyme Q:cytochrome c oxidoreductase complex in Toxoplasma gondii parasites. PLOS Pathogens 17:e1009211. DOI: https://doi.org/10.1371/journal.ppat.1009211, PMID: 33524071 Heaslip AT, Nishi M, Stein B, Hu K. 2011. The motility of a human parasite, Toxoplasma gondii, is regulated by a novel lysine methyltransferase. PLOS Pathogens 7:e1002201. DOI: https://doi.org/10.1371/journal.ppat. 1002201, PMID: 21909263 Heaslip AT, Nelson SR, Warshaw DM, Steinberg G. 2016. Dense granule trafficking in Toxoplasma gondii requires a unique class 27 myosin and actin filaments. Molecular Biology of the Cell 27:2080–2089. DOI: https://doi.org/10.1091/mbc.E15-12-0824, PMID: 27146112 Heredero- Bermejo I, Varberg JM, Charvat R, Jacobs K, Garbuz T, Sullivan WJ, Arrizabalaga G. 2019. TgDrpC, an atypical dynamin- related protein in Toxoplasma gondii, is associated with vesicular transport factors and parasite division. Molecular Microbiology 111:46–64. DOI: https://doi.org/10.1111/mmi.14138, PMID: 30362624 Herneisen AL, Sidik SM, Markus BM, Drewry DH, Zuercher WJ, Lourido S. 2020. Identifying the target of an antiparasitic compound in toxoplasma using thermal proteome profiling. ACS Chemical Biology 15:1801–1807. DOI: https://doi.org/10.1021/acschembio.0c00369, PMID: 32597628 Herneisen AL, Lourido S. 2021. Thermal proteome profiling to identify protein- ligand interactions in the apicomplexan parasite Toxoplasma gondii. Bio- Protocol 11:e4207. DOI: https://doi.org/10.21769/BioProtoc. 4207, PMID: 34859122 Hochmal AK, Zinzius K, Charoenwattanasatien R, Gäbelein P, Mutoh R, Tanaka H, Schulze S, Liu G, Scholz M, Nordhues A, Offenborn JN, Petroutsos D, Finazzi G, Fufezan C, Huang K, Kurisu G, Hippler M. 2016. Calredoxin represents a novel type of calcium- dependent sensor-r esponder connected to redox regulation in the chloroplast. Nature Communications 7:11847. DOI: https://doi.org/10.1038/ncomms11847, PMID: 27297041 Hortua Triana MA, Márquez-N ogueras KM, Chang L, Stasic AJ, Li C, Spiegel KA, Sharma A, Li ZH, Moreno SNJ. 2018. Tagging of weakly expressed Toxoplasma gondii calcium- related genes with high-a ffinity tags. The Journal of Eukaryotic Microbiology 65:709–721. DOI: https://doi.org/10.1111/jeu.12626, PMID: 29672999 Huet D, Rajendran E, van Dooren GG, Lourido S. 2018. Identification of cryptic subunits from an apicomplexan ATP synthase. eLife 7:e38097. DOI: https://doi.org/10.7554/eLife.38097, PMID: 30204085 Hughes CS, Moggridge S, Müller T, Sorensen PH, Morin GB, Krijgsveld J. 2019. Single- pot, solid-p hase- enhanced sample preparation for proteomics experiments. Nature Protocols 14:68–85. DOI: https://doi.org/ 10.1038/s41596-018-0082-x, PMID: 30464214 Hunt A, Russell MRG, Wagener J, Kent R, Carmeille R, Peddie CJ, Collinson L, Heaslip A, Ward GE, Treeck M. 2019. Differential requirements for cyclase- associated protein (CAP) in actin- dependent processes of Toxoplasma gondii eLife 8:e50598. DOI: https://doi.org/10.7554/eLife.50598, PMID: 31577230 Huynh MH, Carruthers VB. 2009. Tagging of endogenous genes in a Toxoplasma gondii strain lacking ku80. Eukaryotic Cell 8:530–539. DOI: https://doi.org/10.1128/EC.00358-08, PMID: 19218426 Ingram JR, Knockenhauer KE, Markus BM, Mandelbaum J, Ramek A, Shan Y, Shaw DE, Schwartz TU, Ploegh HL, Lourido S. 2015. Allosteric activation of apicomplexan calcium-d ependent protein kinases. PNAS 112:E4975– E4984. DOI: https://doi.org/10.1073/pnas.1505914112, PMID: 26305940 Invergo BM, Brochet M, Yu L, Choudhary J, Beltrao P, Billker O. 2017. Sub- minute phosphoregulation of cell cycle systems during plasmodium gamete formation. Cell Reports 21:2017–2029. DOI: https://doi.org/10. 1016/j.celrep.2017.10.071, PMID: 29141230 Jacot D, Daher W, Soldati-F avre D. 2013. Toxoplasma gondii myosin F, an essential motor for centrosomes positioning and apicoplast inheritance. The EMBO Journal 32:1702–1716. DOI: https://doi.org/10.1038/emboj. 2013.113, PMID: 23695356 Jarzab A, Kurzawa N, Hopf T, Moerch M, Zecha J, Leijten N, Bian Y, Musiol E, Maschberger M, Stoehr G, Becher I, Daly C, Samaras P, Mergner J, Spanier B, Angelov A, Werner T, Bantscheff M, Wilhelm M, Klingenspor M, et al. 2020. Meltome atlas-t hermal proteome stability across the tree of life. Nature Methods 17:495–503. DOI: https://doi.org/10.1038/s41592-020-0801-4, PMID: 32284610 Jia Y, Marq J- B, Bisio H, Jacot D, Mueller C, Yu L, Choudhary J, Brochet M, Soldati- Favre D. 2017. Crosstalk between PKA and PKG controls ph- dependent host cell egress of Toxoplasma gondii The EMBO Journal 36:3250–3267. DOI: https://doi.org/10.15252/embj.201796794, PMID: 29030485 Katris NJ, van Dooren GG, McMillan PJ, Hanssen E, Tilley L, Waller RF. 2014. The apical complex provides a regulated gateway for secretion of invasion factors in toxoplasma. PLOS Pathogens 10:e1004074. DOI: https:// doi.org/10.1371/journal.ppat.1004074, PMID: 24743791 Kelley LA, Mezulis S, Yates CM, Wass MN, Sternberg MJE. 2015. The phyre2 web portal for protein modeling, prediction and analysis. Nature Protocols 10:845–858. DOI: https://doi.org/10.1038/nprot.2015.053, PMID: 25950237 Herneisen et al. eLife 2022;11:e80336. DOI: https://doi.org/10.7554/eLife.80336 37 of 41 Research article Cell Biology | Microbiology and Infectious Disease Kenthirapalan S, Waters AP, Matuschewski K, Kooij TWA. 2014. Copper- transporting atpase is important for malaria parasite fertility. Molecular Microbiology 91:315–325. DOI: https://doi.org/10.1111/mmi.12461, PMID: 24237419 Koreny L, Zeeshan M, Barylyuk K, Tromer EC, van Hooff JJE, Brady D, Ke H, Chelaghma S, Ferguson DJP, Eme L, Tewari R, Waller RF. 2021. Molecular characterization of the conoid complex in toxoplasma reveals its conservation in all apicomplexans, including plasmodium species. PLOS Biology 19:e3001081. DOI: https://doi. org/10.1371/journal.pbio.3001081, PMID: 33705380 Lehane AM, Dennis ASM, Bray KO, Li D, Rajendran E, McCoy JM, McArthur HM, Winterberg M, Rahimi F, Tonkin CJ, Kirk K, van Dooren GG. 2019. Characterization of the ATP4 ion pump in Toxoplasma gondii The Journal of Biological Chemistry 294:5720–5734. DOI: https://doi.org/10.1074/jbc.RA118.006706, PMID: 30723156 Lentini G, Ben Chaabene R, Vadas O, Ramakrishnan C, Mukherjee B, Mehta V, Lunghi M, Grossmann J, Maco B, Visentin R, Hehl AB, Korkhov VM, Soldati-F avre D. 2021. Structural insights into an atypical secretory pathway kinase crucial for Toxoplasma gondii invasion. Nature Communications 12:3788. DOI: https://doi.org/10.1038/ s41467-021-24083-y, PMID: 34145271 Leung JM, He Y, Zhang F, Hwang YC, Nagayasu E, Liu J, Murray JM, Hu K. 2017. Stability and function of a putative microtubule- organizing center in the human parasite Toxoplasma gondii Molecular Biology of the Cell 28:1361–1378. DOI: https://doi.org/10.1091/mbc.E17-01-0045, PMID: 28331073 Li M, Wang H, Liu J, Hao P, Ma L, Liu Q. 2015. The apoptotic role of metacaspase in Toxoplasma gondii. Frontiers in Microbiology 6:1560. DOI: https://doi.org/10.3389/fmicb.2015.01560, PMID: 26834715 Li J, Van Vranken JG, Pontano Vaites L, Schweppe DK, Huttlin EL, Etienne C, Nandhikonda P, Viner R, Robitaille AM, Thompson AH, Kuhn K, Pike I, Bomgarden RD, Rogers JC, Gygi SP, Paulo JA. 2020. TMTpro reagents: a set of isobaric labeling mass tags enables simultaneous proteome- wide measurements across 16 samples. Nature Methods 17:399–404. DOI: https://doi.org/10.1038/s41592-020-0781-4, PMID: 32203386 Li ZH, King TP, Ayong L, Asady B, Cai X, Rahman T, Vella SA, Coppens I, Patel S, Moreno SNJ. 2021. A plastid two- pore channel essential for inter- organelle communication and growth of Toxoplasma gondii. Nature Communications 12:5802. DOI: https://doi.org/10.1038/s41467-021-25987-5, PMID: 34608145 Li W, Grech J, Stortz JF, Gow M, Periz J, Meissner M, Jimenez- Ruiz E. 2022. A splitcas9 phenotypic screen in Toxoplasma gondii identifies proteins involved in host cell egress and invasion. Nature Microbiology 7:882– 895. DOI: https://doi.org/10.1038/s41564-022-01114-y, PMID: 35538310 Lim YT, Prabhu N, Dai L, Go KD, Chen D, Sreekumar L, Egeblad L, Eriksson S, Chen L, Veerappan S, Teo HL, Tan CSH, Lengqvist J, Larsson A, Sobota RM, Nordlund P. 2018. An efficient proteome- wide strategy for discovery and characterization of cellular nucleotide-p rotein interactions. PLOS ONE 13:e0208273. DOI: https://doi.org/10.1371/journal.pone.0208273, PMID: 30521565 Long S, Anthony B, Drewry LL, Sibley LD. 2017a. A conserved ankyrin repeat- containing protein regulates conoid stability, motility and cell invasion in Toxoplasma gondii. Nature Communications 8:2236. DOI: https://doi.org/ 10.1038/s41467-017-02341-2, PMID: 29269729 Long S, Brown KM, Drewry LL, Anthony B, Phan IQH, Sibley LD. 2017b. Calmodulin- like proteins localized to the conoid regulate motility and cell invasion by Toxoplasma gondii. PLOS Pathogens 13:e1006379. DOI: https:// doi.org/10.1371/journal.ppat.1006379, PMID: 28475612 Lourido S, Shuman J, Zhang C, Shokat KM, Hui R, Sibley LD. 2010. Calcium-d ependent protein kinase 1 is an essential regulator of exocytosis in toxoplasma. Nature 465:359–362. DOI: https://doi.org/10.1038/ nature09022, PMID: 20485436 Lourido S, Tang K, Sibley LD. 2012. Distinct signalling pathways control toxoplasma egress and host- cell invasion. The EMBO Journal 31:4524–4534. DOI: https://doi.org/10.1038/emboj.2012.299, PMID: 23149386 Lourido S, Moreno SNJ. 2015. The calcium signaling toolkit of the apicomplexan parasites Toxoplasma gondii and plasmodium spp. Cell Calcium 57:186–193. DOI: https://doi.org/10.1016/j.ceca.2014.12.010, PMID: 25605521 Luan S, Wang C. 2021. Calcium signaling mechanisms across kingdoms. Annual Review of Cell and Developmental Biology 37:311–340. DOI: https://doi.org/10.1146/annurev-cellbio-120219-035210, PMID: 34375534 Luo S, Vieira M, Graves J, Zhong L, Moreno SN. 2001. A plasma membrane- type ca(2+)- atpase co-l ocalizes with A vacuolar H(+)- pyrophosphatase to acidocalcisomes of Toxoplasma gondii. The EMBO Journal 20:55–64. DOI: https://doi.org/10.1093/emboj/20.1.55, PMID: 11226155 Luo S, Ruiz FA, Moreno SNJ. 2005. The acidocalcisome ca2+-atpase (tga1) of Toxoplasma gondii is required for polyphosphate storage, intracellular calcium homeostasis and virulence. Molecular Microbiology 55:1034– 1045. DOI: https://doi.org/10.1111/j.1365-2958.2004.04464.x, PMID: 15686552 Markus BM, Bell GW, Lorenzi HA, Lourido S. 2019. Optimizing systems for cas9 expression in Toxoplasma gondii. MSphere 4:e00386- 19. DOI: https://doi.org/10.1128/mSphere.00386-19, PMID: 31243081 Márquez- Nogueras KM, Hortua Triana MA, Chasen NM, Kuo IY, Moreno SN. 2021. Calcium signaling through a transient receptor channel is important for Toxoplasma gondii growth. eLife 10:e63417. DOI: https://doi.org/ 10.7554/eLife.63417, PMID: 34106044 Mateus A, Kurzawa N, Becher I, Sridharan S, Helm D, Stein F, Typas A, Savitski MM. 2020. Thermal proteome profiling for interrogating protein interactions. Molecular Systems Biology 16:e9232. DOI: https://doi.org/10. 15252/msb.20199232, PMID: 32133759 Herneisen et al. eLife 2022;11:e80336. DOI: https://doi.org/10.7554/eLife.80336 38 of 41 Research article Cell Biology | Microbiology and Infectious Disease Matthiesen SH, Shenoy SM, Kim K, Singer RH, Satir BH. 2001. A parafusin-r elated toxoplasma protein in ca2+-regulated secretory organelles. European Journal of Cell Biology 80:775–783. DOI: https://doi.org/10. 1078/0171-9335-00214, PMID: 11831391 McCoy JM, Whitehead L, van Dooren GG, Tonkin CJ. 2012. TgCDPK3 regulates calcium- dependent egress of Toxoplasma gondii from host cells. PLOS Pathogens 8:e1003066. DOI: https://doi.org/10.1371/journal.ppat. 1003066, PMID: 23226109 McCoy JM, Stewart RJ, Uboldi AD, Li D, Schröder J, Scott NE, Papenfuss AT, Lehane AM, Foster LJ, Tonkin CJ. 2017. A forward genetic screen identifies A negative regulator of rapid ca2+- dependent cell egress (MS1) in the intracellular parasite Toxoplasma gondii The Journal of Biological Chemistry 292:7662–7674. DOI: https://doi. org/10.1074/jbc.M117.775114, PMID: 28258212 Meier A, Söding J. 2015. Automatic prediction of protein 3D structures by probabilistic multi- template homology modeling. PLOS Computational Biology 11:e1004343. DOI: https://doi.org/10.1371/journal.pcbi. 1004343, PMID: 26496371 Melatti C, Pieperhoff M, Lemgruber L, Pohl E, Sheiner L, Meissner M. 2019. A unique dynamin-r elated protein is essential for mitochondrial fission in Toxoplasma gondii. PLOS Pathogens 15:e1007512. DOI: https://doi.org/ 10.1371/journal.ppat.1007512, PMID: 30947298 Møller AB, Asp T, Holm PB, Palmgren MG. 2008. Phylogenetic analysis of P5 P- type atpases, a eukaryotic lineage of secretory pathway pumps. Molecular Phylogenetics and Evolution 46:619–634. DOI: https://doi.org/ 10.1016/j.ympev.2007.10.023, PMID: 18155930 Moreno SN, Zhong L. 1996. Acidocalcisomes in Toxoplasma gondii tachyzoites. The Biochemical Journal 313 (Pt 2):655–659. DOI: https://doi.org/10.1042/bj3130655, PMID: 8573106 Morlon- Guyot J, Berry L, Chen CT, Gubbels MJ, Lebrun M, Daher W. 2014. The Toxoplasma gondii calcium- dependent protein kinase 7 is involved in early steps of parasite division and is crucial for parasite survival. Cellular Microbiology 16:95–114. DOI: https://doi.org/10.1111/cmi.12186, PMID: 24011186 Moss WJ, Patterson CE, Jochmans AK, Brown KM. 2022. Functional analysis of the expanded phosphodiesterase gene family in Toxoplasma gondii tachyzoites. MSphere 7:e0079321. DOI: https://doi.org/10.1128/msphere. 00793-21, PMID: 35107337 Moudy R, Manning TJ, Beckers CJ. 2001. The loss of cytoplasmic potassium upon host cell breakdown triggers egress of Toxoplasma gondii. The Journal of Biological Chemistry 276:41492–41501. DOI: https://doi.org/10. 1074/jbc.M106154200, PMID: 11526113 Mueller C, Klages N, Jacot D, Santos JM, Cabrera A, Gilberger TW, Dubremetz JF, Soldati- Favre D. 2013. The toxoplasma protein ARO mediates the apical positioning of rhoptry organelles, a prerequisite for host cell invasion. Cell Host & Microbe 13:289–301. DOI: https://doi.org/10.1016/j.chom.2013.02.001, PMID: 23498954 Mueller C, Samoo A, Hammoudi PM, Klages N, Kallio JP, Kursula I, Soldati- Favre D. 2016. Structural and functional dissection of Toxoplasma gondii armadillo repeats only protein. Journal of Cell Science 129:1031– 1045. DOI: https://doi.org/10.1242/jcs.177386, PMID: 26769898 Mühleip A, Kock Flygaard R, Ovciarikova J, Lacombe A, Fernandes P, Sheiner L, Amunts A. 2021. ATP synthase hexamer assemblies shape cristae of toxoplasma mitochondria. Nature Communications 12:120. DOI: https:// doi.org/10.1038/s41467-020-20381-z, PMID: 33402698 Nagamune K, Sibley LD. 2006. Comparative genomic and phylogenetic analyses of calcium atpases and calcium- regulated proteins in the apicomplexa. Molecular Biology and Evolution 23:1613–1627. DOI: https:// doi.org/10.1093/molbev/msl026, PMID: 16751258 Nebl T, Prieto JH, Kapp E, Smith BJ, Williams MJ, Yates JR, Cowman AF, Tonkin CJ. 2011. Quantitative in vivo analyses reveal calcium- dependent phosphorylation sites and identifies a novel component of the toxoplasma invasion motor complex. PLOS Pathogens 7:e1002222. DOI: https://doi.org/10.1371/journal.ppat.1002222, PMID: 21980283 O’Shaughnessy WJ, Hu X, Beraki T, McDougal M, Reese ML. 2020. Loss of a conserved MAPK causes catastrophic failure in assembly of a specialized cilium- like structure in Toxoplasma gondii Molecular Biology of the Cell 31:881–888. DOI: https://doi.org/10.1091/mbc.E19-11-0607, PMID: 32073987 Pace DA, McKnight CA, Liu J, Jimenez V, Moreno SNJ. 2014. Calcium entry in Toxoplasma gondii and its enhancing effect of invasion- linked traits. The Journal of Biological Chemistry 289:19637–19647. DOI: https:// doi.org/10.1074/jbc.M114.565390, PMID: 24867952 Pace DA, Moreno SNJ, Lourido S. 2020. Chapter 12 - Calcium storage and homeostasis in Toxoplasma gondii. Weiss LM, Kim K (Eds). Toxoplasma gondii. Elsevier. p. 547–575. DOI: https://doi.org/10.1016/b978-0-12- 815041-2.00012-8 Park HS, Lee SC, Cardenas ME, Heitman J. 2019. Calcium-c almodulin-c alcineurin signaling: A globally conserved virulence cascade in eukaryotic microbial pathogens. Cell Host & Microbe 26:453–462. DOI: https://doi.org/10. 1016/j.chom.2019.08.004, PMID: 31600499 Parker KER, Fairweather SJ, Rajendran E, Blume M, McConville MJ, Bröer S, Kirk K, van Dooren GG. 2019. The tyrosine transporter of Toxoplasma gondii is a member of the newly defined apicomplexan amino acid transporter (apiat) family. PLOS Pathogens 15:e1007577. DOI: https://doi.org/10.1371/journal.ppat.1007577, PMID: 30742695 Patel A, Perrin AJ, Flynn HR, Bisson C, Withers- Martinez C, Treeck M, Flueck C, Nicastro G, Martin SR, Ramos A, Gilberger TW, Snijders AP, Blackman MJ, Baker DA. 2019. Cyclic AMP signalling controls key components of malaria parasite host cell invasion machinery. PLOS Biology 17:e3000264. DOI: https://doi.org/10.1371/journal. pbio.3000264, PMID: 31075098 Herneisen et al. eLife 2022;11:e80336. DOI: https://doi.org/10.7554/eLife.80336 39 of 41 Research article Cell Biology | Microbiology and Infectious Disease Paul AS, Saha S, Engelberg K, Jiang RHY, Coleman BI, Kosber AL, Chen CT, Ganter M, Espy N, Gilberger TW, Gubbels MJ, Duraisingh MT. 2015. Parasite calcineurin regulates host cell recognition and attachment by apicomplexans. Cell Host & Microbe 18:49–60. DOI: https://doi.org/10.1016/j.chom.2015.06.003, PMID: 26118996 Paul AS, Miliu A, Paulo JA, Goldberg JM, Bonilla AM, Berry L, Seveno M, Braun- Breton C, Kosber AL, Elsworth B, Arriola JSN, Lebrun M, Gygi SP, Lamarque MH, Duraisingh MT. 2020. Co- option of Plasmodium falciparum PP1 for egress from host erythrocytes. Nature Communications 11:3532. DOI: https://doi.org/10. 1038/s41467-020-17306-1, PMID: 32669539 Perez- Riverol Y, Bai J, Bandla C, García- Seisdedos D, Hewapathirana S, Kamatchinathan S, Kundu DJ, Prakash A, Frericks- Zipper A, Eisenacher M, Walzer M, Wang S, Brazma A, Vizcaíno JA. 2022. The pride database resources in 2022: a hub for mass spectrometry- based proteomics evidences. Nucleic Acids Research 50:D543–D552. DOI: https://doi.org/10.1093/nar/gkab1038, PMID: 34723319 Philip N, Waters AP. 2015. Conditional degradation of plasmodium calcineurin reveals functions in parasite colonization of both host and vector. Cell Host & Microbe 18:122–131. DOI: https://doi.org/10.1016/j.chom. 2015.05.018, PMID: 26118994 Reinhard FBM, Eberhard D, Werner T, Franken H, Childs D, Doce C, Savitski MF, Huber W, Bantscheff M, Savitski MM, Drewes G. 2015. Thermal proteome profiling monitors ligand interactions with cellular membrane proteins. Nature Methods 12:1129–1131. DOI: https://doi.org/10.1038/nmeth.3652, PMID: 26524241 Ritz C, Baty F, Streibig JC, Gerhard D. 2015. Dose- response analysis using R. PLOS ONE 10:e0146021. DOI: https://doi.org/10.1371/journal.pone.0146021, PMID: 26717316 Roiko MS, Carruthers VB. 2013. Functional dissection of Toxoplasma gondii perforin-l ike protein 1 reveals a dual domain mode of membrane binding for cytolysis and parasite egress. The Journal of Biological Chemistry 288:8712–8725. DOI: https://doi.org/10.1074/jbc.M113.450932, PMID: 23376275 Salunke R, Mourier T, Banerjee M, Pain A, Shanmugam D, Striepen B. 2018. Highly diverged novel subunit composition of apicomplexan F- type ATP synthase identified from Toxoplasma gondii. PLOS Biology 16:e2006128. DOI: https://doi.org/10.1371/journal.pbio.2006128, PMID: 30005062 Savitski MM, Reinhard FBM, Franken H, Werner T, Savitski MF, Eberhard D, Martinez Molina D, Jafari R, Dovega RB, Klaeger S, Kuster B, Nordlund P, Bantscheff M, Drewes G. 2014. Tracking cancer drugs in living cells by thermal profiling of the proteome. Science 346:1255784. DOI: https://doi.org/10.1126/science. 1255784, PMID: 25278616 Scrucca L, Fop M, Murphy T, Raftery A. 2016. Mclust 5: clustering, classification and density estimation using gaussian finite mixture models. The R Journal 8:289. DOI: https://doi.org/10.32614/RJ-2016-021 Seidi A, Muellner- Wong LS, Rajendran E, Tjhin ET, Dagley LF, Aw VY, Faou P, Webb AI, Tonkin CJ, van Dooren GG. 2018. Elucidating the mitochondrial proteome of Toxoplasma gondii reveals the presence of a divergent cytochrome c oxidase. eLife 7:e38131. DOI: https://doi.org/10.7554/eLife.38131, PMID: 30204084 Shi Y. 2009. Serine/threonine phosphatases: mechanism through structure. Cell 139:468–484. DOI: https://doi. org/10.1016/j.cell.2009.10.006, PMID: 19879837 Shortt E, Lourido S. 2020. Plate- based quantification of stimulated toxoplasma egress. Methods in Molecular Biology 2071:171–186. DOI: https://doi.org/10.1007/978-1-4939-9857-9_10, PMID: 31758453 Sidik SM, Huet D, Ganesan SM, Huynh MH, Wang T, Nasamu AS, Thiru P, Saeij JPJ, Carruthers VB, Niles JC, Lourido S. 2016a. A genome- wide CRISPR screen in toxoplasma identifies essential apicomplexan genes. Cell 166:1423–1435. DOI: https://doi.org/10.1016/j.cell.2016.08.019, PMID: 27594426 Sidik S, Hortua Triana MA, Paul AS, El Bakkouri M, Hackett CG, Tran F, Westwood NJ, Hui R, Zuercher WJ, Duraisingh MT, Moreno SNJ, Lourido S. 2016b. Using a genetically encoded sensor to identify inhibitors of Toxoplasma gondii ca2+ signaling. The Journal of Biological Chemistry 291:9566–9580. DOI: https://doi.org/ 10.1074/jbc.M115.703546, PMID: 26933036 Silmon de Monerri NC, Yakubu RR, Chen AL, Bradley PJ, Nieves E, Weiss LM, Kim K. 2015. The ubiquitin proteome of Toxoplasma gondii reveals roles for protein ubiquitination in cell-c ycle transitions. Cell Host & Microbe 18:621–633. DOI: https://doi.org/10.1016/j.chom.2015.10.014, PMID: 26567513 Smith TA, Lopez- Perez GS, Herneisen AL, Shortt E, Lourido S. 2022. Screening the toxoplasma kinome with high- throughput tagging identifies a regulator of invasion and egress. Nature Microbiology 7:868–881. DOI: https://doi.org/10.1038/s41564-022-01104-0, PMID: 35484233 Sridharan S, Kurzawa N, Werner T, Günthner I, Helm D, Huber W, Bantscheff M, Savitski MM. 2019. Proteome- wide solubility and thermal stability profiling reveals distinct regulatory roles for ATP. Nature Communications 10:1155. DOI: https://doi.org/10.1038/s41467-019-09107-y, PMID: 30858367 Stewart RJ, Whitehead L, Nijagal B, Sleebs BE, Lessene G, McConville MJ, Rogers KL, Tonkin CJ. 2017. Analysis of ca2+ mediated signaling regulating toxoplasma infectivity reveals complex relationships between key molecules. Cellular Microbiology 19:12685. DOI: https://doi.org/10.1111/cmi.12685 Stommel EW, Ely KH, Schwartzman JD, Kasper LH. 1997. Toxoplasma gondii:dithiol-i nduced ca2+flux causes egress of parasites from the parasitophorous vacuole. Experimental Parasitology 87:88–97. DOI: https://doi. org/10.1006/expr.1997.4187, PMID: 9326884 Subramanian SV, Satir BH. 1992. Carbohydrate cycling in signal transduction: parafusin, a phosphoglycoprotein and possible ca(2+)- dependent transducer molecule in exocytosis in paramecium. PNAS 89:11297–11301. DOI: https://doi.org/10.1073/pnas.89.23.11297, PMID: 1333606 Subramanian SV, Wyroba E, Andersen AP, Satir BH. 1994. Cloning and sequencing of parafusin, a calcium- dependent exocytosis- related phosphoglycoprotein. PNAS 91:9832–9836. DOI: https://doi.org/10.1073/pnas. 91.21.9832, PMID: 7937900 Herneisen et al. eLife 2022;11:e80336. DOI: https://doi.org/10.7554/eLife.80336 40 of 41 Research article Cell Biology | Microbiology and Infectious Disease Tosetti N, Dos Santos Pacheco N, Bertiaux E, Maco B, Bournonville L, Hamel V, Guichard P, Soldati- Favre D. 2020. Essential function of the alveolin network in the subpellicular microtubules and conoid assembly in Toxoplasma gondii eLife 9:e56635. DOI: https://doi.org/10.7554/eLife.56635, PMID: 32379047 Treeck M, Sanders JL, Gaji RY, LaFavers KA, Child MA, Arrizabalaga G, Elias JE, Boothroyd JC. 2014. The calcium- dependent protein kinase 3 of toxoplasma influences basal calcium levels and functions beyond egress as revealed by quantitative phosphoproteome analysis. PLOS Pathogens 10:e1004197. DOI: https://doi.org/ 10.1371/journal.ppat.1004197, PMID: 24945436 Tsai CF, Hsu CC, Hung JN, Wang YT, Choong WK, Zeng MY, Lin PY, Hong RW, Sung TY, Chen YJ. 2014. Sequential phosphoproteomic enrichment through complementary metal- directed immobilized metal ion affinity chromatography. Analytical Chemistry 86:685–693. DOI: https://doi.org/10.1021/ac4031175, PMID: 24313913 Uboldi AD, McCoy JM, Blume M, Gerlic M, Ferguson DJP, Dagley LF, Beahan CT, Stapleton DI, Gooley PR, Bacic A, Masters SL, Webb AI, McConville MJ, Tonkin CJ. 2015. Regulation of starch stores by a ca(2+)- dependent protein kinase is essential for viable cyst development in Toxoplasma gondii. Cell Host & Microbe 18:670–681. DOI: https://doi.org/10.1016/j.chom.2015.11.004, PMID: 26651943 Uboldi AD, Wilde ML, McRae EA, Stewart RJ, Dagley LF, Yang L, Katris NJ, Hapuarachchi SV, Coffey MJ, Lehane AM, Botte CY, Waller RF, Webb AI, McConville MJ, Tonkin CJ. 2018. Protein kinase A negatively regulates ca2+ signalling in Toxoplasma gondii. PLOS Biology 16:e2005642. DOI: https://doi.org/10.1371/ journal.pbio.2005642, PMID: 30208022 Varberg JM, Coppens I, Arrizabalaga G, Gaji RY. 2018. TgTKL1 is a unique plant- like nuclear kinase that plays an essential role in acute toxoplasmosis. MBio 9:e00301-1 8. DOI: https://doi.org/10.1128/mBio.00301-18, PMID: 29559568 Vella SA, Moore CA, Li ZH, Hortua Triana MA, Potapenko E, Moreno SNJ. 2021. The role of potassium and host calcium signaling in Toxoplasma gondii egress. Cell Calcium 94:102337. DOI: https://doi.org/10.1016/j.ceca. 2020.102337, PMID: 33524795 Villalobo A, González-M uñoz M, Berchtold MW. 2019. Proteins with calmodulin- like domains: structures and functional roles. Cellular and Molecular Life Sciences 76:2299–2328. DOI: https://doi.org/10.1007/s00018-019- 03062-z, PMID: 30877334 Vo KC, Günay- Esiyok Ö, Liem N, Gupta N. 2020. The protozoan parasite Toxoplasma gondii encodes a gamut of phosphodiesterases during its lytic cycle in human cells. Computational and Structural Biotechnology Journal 18:3861–3876. DOI: https://doi.org/10.1016/j.csbj.2020.11.024, PMID: 33335684 Waldman BS, Schwarz D, Wadsworth MH, Saeij JP, Shalek AK, Lourido S. 2020. Identification of a master regulator of differentiation in toxoplasma. Cell 180:359–372. DOI: https://doi.org/10.1016/j.cell.2019.12.013, PMID: 31955846 Wallbank BA, Dominicus CS, Broncel M, Legrave N, Kelly G, MacRae JI, Staines HM, Treeck M. 2019. Characterisation of the Toxoplasma gondii tyrosine transporter and its phosphorylation by the calcium- dependent protein kinase 3. Molecular Microbiology 111:1167–1181. DOI: https://doi.org/10.1111/mmi. 14156, PMID: 30402958 Weiss LM, Wastling J, Jeffers V, Sullivan WJ, Kim K. 2020. Chapter 22 - proteomics and posttranslational protein modifications in Toxoplasma gondii. Weiss LM, Kim K (Eds). Toxoplasma gondii. Third Edition. Academic Press. p. 983–1020. DOI: https://doi.org/10.1016/B978-0-12-815041-2.00022-0 Wernimont AK, Artz JD, Finerty P, Lin YH, Amani M, Allali- Hassani A, Senisterra G, Vedadi M, Tempel W, Mackenzie F, Chau I, Lourido S, Sibley LD, Hui R. 2010. Structures of apicomplexan calcium- dependent protein kinases reveal mechanism of activation by calcium. Nature Structural & Molecular Biology 17:596–601. DOI: https://doi.org/10.1038/nsmb.1795, PMID: 20436473 Wickstead B, Gull K, Richards TA. 2010. Patterns of kinesin evolution reveal a complex ancestral eukaryote with a multifunctional cytoskeleton. BMC Evolutionary Biology 10:110. DOI: https://doi.org/10.1186/1471-2148-10- 110, PMID: 20423470 Yang C, Arrizabalaga G. 2017. The serine/threonine phosphatases of apicomplexan parasites. Molecular Microbiology 106:1–21. DOI: https://doi.org/10.1111/mmi.13715, PMID: 28556455 Yang C, Broncel M, Dominicus C, Sampson E, Blakely WJ, Treeck M, Arrizabalaga G. 2019. A plasma membrane localized protein phosphatase in Toxoplasma gondii, PPM5C, regulates attachment to host cells. Scientific Reports 9:5924. DOI: https://doi.org/10.1038/s41598-019-42441-1, PMID: 30976120 Zeeshan M, Pandey R, Subudhi AK, Ferguson DJP, Kaur G, Rashpa R, Nugmanova R, Brady D, Bottrill AR, Vaughan S, Brochet M, Bollen M, Pain A, Holder AA, Guttery DS, Tewari R. 2021. Protein phosphatase 1 regulates atypical mitotic and meiotic division in plasmodium sexual stages. Communications Biology 4:760. DOI: https://doi.org/10.1038/s42003-021-02273-0, PMID: 34145386 Zeng JM, Hapuarachchi SV, Shafik SH, Martin RE, Kirk K, van Dooren GG, Lehane AM. 2021. Identifying the major lactate transporter of Toxoplasma gondii tachyzoites. Scientific Reports 11:6787. DOI: https://doi.org/10. 1038/s41598-021-86204-3, PMID: 33762657 Herneisen et al. eLife 2022;11:e80336. DOI: https://doi.org/10.7554/eLife.80336 41 of 41