Introduction

Apicomplexans parasitize a majority of warm-blooded species, including an estimated quarter of the human population (Jones et al., 2014). Central to parasite success is the ability to transition between different life cycle stages to balance proliferation and transmission (Lourido, 2019). The acute phase of infection by T. gondii, a model apicomplexan, induces the pathology of disease. Parasites invade host cells, establish an intracellular niche that subverts innate immune responses, replicate, and exit the host via lysis to establish new sites of infection (Blader et al., 2015). These pathogens transcend the inverted biomolecular gradients within and without the host cell in minutes. Such rapid changes in cellular state depend on signal transduction in T. gondii and related apicomplexans. Second messenger cascades within the parasite cytoplasm are especially well-suited to rapidly transduce and amplify signals from a changing environment. Accordingly, apicomplexans—and T. gondii in particular—have an expanded repertoire of second messenger signaling components, from kinases to cyclases and phosphodiesterases (Brown et al., 2020; Pace et al., 2020).

The roles of second messenger kinases—collectively referred to as AGC kinases for the founding members protein kinases A, G, and C—in apicomplexan life cycles have been elaborated in candidate-by-candidate approaches. cGMP-dependent protein kinase (or protein kinase G, PKG) is necessary for the secretion of parasite-specific adhesins and perforins enabling efficient escape and invasion of host cells (Brown et al., 2017, 2016; Wiersma et al., 2004). In apicomplexan cells, PKG operates upstream of calcium release and thus potentiates an orthogonal second messenger signaling network. Consequently, PKG inhibitors have been pursued as anti-parasitic compounds (Baker et al., 2017; Donald et al., 2006; Sidik et al., 2016b). Another cyclic nucleotide-dependent protein kinase, PKA, similarly promotes parasite spread during infection; in Plasmodium spp. merozoites, PKA function enables parasite invasion of host cells (Flueck et al., 2019; Patel et al., 2019; Wilde et al., 2019), and in T. gondii balanced PKA C1 function is required to ensure that parasites do not exit the host cell prematurely (Jia et al., 2017; Uboldi et al., 2018). T. gondii has an additional ortholog of the PKA catalytic subunit, PKA C3, which maintains parasites in the acute phase of the infection and reduces conversion to a slowly proliferating, chronic infection termed the bradyzoite stage (Sugi et al., 2016).

Several recent studies have sought to characterize the downstream targets of second messenger kinases in apicomplexans, including PKA, PKG, and calcium-dependent protein kinases (Alam et al., 2015; Balestra et al., 2021; Brochet et al., 2014; Chan et al., 2023; Herneisen et al., 2022; Jia et al., 2017; Nofal et al., 2022). The second messenger networks in these parasites are intricately interwoven, with documented crosstalk between parasite cAMP, cGMP, and calcium signaling cascades. Comparatively little work has addressed upstream signal integration in apicomplexans. In metazoans, growth factor kinases such as mammalian target of rapamycin (mTOR) and phosphoinositide-dependent protein kinase 1 (PDK1) prime numerous AGC kinases for further second messenger-based activation (Laplante and Sabatini, 2012; Mora et al., 2004).

We recently identified the Store Potentiating/Activating Regulatory Kinase (SPARK), an apicomplexan ortholog of metazoan PDK1, as a fitness-conferring candidate in pooled screens of the T. gondii kinome. Parasites depleted of SPARK for multiple replication cycles failed to enter and exit host cells upon stimulation with zaprinast, a chemical that elevates cGMP in parasites (Smith et al., 2022). In a concurrent study, chemical mutagenesis implicated the Plasmodium falciparum ortholog of SPARK in invasion and proliferation through the regulation of the PKA pathway (Hitz et al., 2021). Although parasites with disrupted SPARK alleles were characterized phenotypically, evidence for the PDK1 activity of SPARK was indirect. Here, we map the topology of the signaling network regulated by SPARK in T. gondii. SPARK complexes with an elongin-like protein we name SPARKEL. Phenotypic and proteomic experiments show that SPARK dysregulation lessens PKA C1, PKG, and PKA C3 function. Phenotypes attributed to SPARK depletion—reduced invasion, host cell lysis, calcium signaling, and elevated differentiation—can be explained by the attenuated activities of SPARK’s client kinases. The proper functioning of SPARK thus ensures that T. gondii progresses through the acute stage of its infection cycle.

Results

SPARK complexes with an elongin-like protein, SPARKEL

To identify proteins interacting with SPARK, we immunopurified the mNG epitope of the previously described SPARK-mNG-AID strain (Smith et al., 2022) with anti-mNG nanobodies and quantified protein abundance using label-free quantitative proteomics (Figure 1—figure supplement 1 and Supplementary File 1). The immunoprecipitation (IP) successfully enriched SPARK as well as a hypothetical protein, TGGT1_291010. IPs with lysates of parasites expressing SPARK-mNG lacking AID demonstrated that the interaction was not due to the presence of the degron; TGGT1_291010 was once again highly enriched (Figure 1A), along with a putative AGC kinase (TGGT1_205550). The AGC kinase is restricted to coccidians but is dispensable during acute T. gondii infection (Sidik et al., 2016a); we therefore focused our efforts on the hypothetical protein, which contributes to parasite fitness.

Figure 1.

TGGT1_291010 is a 23 kDa protein with a C-terminal SKP1/BTB/POZ domain (Figure 1B). This domain shares homology with that seen in metazoan Elongin C. For this reason, we named the gene TGGT1_291010 “SPARK Elongin-Like protein” (SPARKEL, Table 1). SPARKEL homologs are found in the free-living alveolates Vitrella brassicaformis and Chromera velia but are absent in Aconoidasida (Figure 1C), suggesting that the gene was present in the ancestor of apicomplexans but was selectively lost in some lineages.

Table 1.

Immunoprecipitation and mass spectrometry of SPARKEL enriched SPARK (Figure 1D and Supplementary File 1) and several putative components of daughter cell replication, including TGGT1_212770 (Dos Santos Pacheco et al., 2021), BCC10, and EFPP/BCC5. To confirm that SPARKEL and SPARK interact in situ, we introduced a C-terminal TurboID fusion at the SPARKEL endogenous locus (Branon et al., 2018). These transgenic parasites were treated with biotin while intracellular, and biotinylated proteins were enriched from lysates via streptavidin precipitation. As a control, we conducted a parallel experiment with transgenic parasites expressing TurboID fused to a cytosolic mNG fluorophore. SPARK and the AGC kinase TGGT1_205550 were highly enriched in the SPARKEL proximity-labeling experiment (Figure 1E and Supplementary File 2). A proximity labeling experiment with parasites expressing TurboID at the N terminus of SPARKEL similarly enriched SPARK and TGGT1_205550 (Figure 1—figure supplement 1 and Supplementary File 2). Thus, the interactions between SPARK, SPARKEL, and TGGT1_205550 were consistently observed across a variety of approaches.

Endogenous tagging of SPARK and SPARKEL in the same parasite strain revealed punctate cytosolic staining for each protein (Figure 1F); however, the low expression level of SPARKEL precluded robust colocalization. Nevertheless, strong mass spectrometry evidence for the physical interaction between SPARK and SPARKEL motivated further analysis of their functional relationship.

SPARK and SPARKEL depletion phenocopies at multiple steps in the lytic cycle

To determine the regulatory interaction between SPARK and SPARKEL, we first generated a SPARKEL conditional knockdown allele by inserting V5, HaloTag, mini auxin-inducible degron (mAID), and Ty epitopes at the endogenous SPARKEL C terminus. In this strain, we tagged SPARK with a C-terminal V5-mCherry-HA. Parasites treated with IAA to deplete SPARKEL failed to form plaques in host cell monolayers (Figure 2A). As indicated by immunoblot, IAA treatment led to a reduction in detectable SPARKEL signal within one hour, as well as co-depletion of SPARK on a similar time scale (Figure 2B). We next generated a strain expressing SPARK-V5-mAID-HA. Similarly, SPARK levels were reduced within one hour of IAA addition (Figure 2—figure supplement 1), as reported for a similar strain previously (Smith et al., 2022). In this genetic background, we tagged SPARKEL endogenously with a C-terminal V5-mNG-Ty. IAA treatment reduced SPARK levels and led to a reduction in SPARKEL abundance (Figure 2C and Figure 2—figure supplement 1). Together, these results suggest that the AID system leads to co-depletion of SPARK and SPARKEL.

Figure 2.

To assess the impact of complex knockdown on parasite replication and transitions between intracellular and extracellular environments, we treated SPARKEL-AID and SPARK-AID parasites with IAA and performed replication, invasion, and egress assays. SPARKEL-AID parasites exhibited normal replication kinetics after 24h of IAA treatment (Figure 2D). Invasion and egress efficiency of the AID- tagged strains were reduced after parasites were treated with IAA for 24 hours (Figure 2E and 2F), as previously reported for SPARK (Smith et al., 2022). Intracellular calcium measurements of GCaMP- expressing parasites revealed a similar trend, with SPARKEL- and SPARK-depleted cells exhibiting reduced calcium mobilization after 24 hours of IAA and treatment when stimulated with zaprinast, a compound that stimulates calcium release from parasite intracellular stores (Brown et al., 2016; Sidik et al., 2016b) (Figure 2G, 2H, and Figure 2—figure supplement 1).

SPARK was previously identified in a genetic screen that distinguished acute and delayed death phenotypes in a pooled population of parasites (Smith et al., 2022). As parasites with conditional SPARK alleles were assessed only in the context of delayed death phenotypes—i.e. after several replication cycles with IAA treatment—we characterized the ability of SPARK- and SPARKEL-AID parasites to invade and egress after 3 hours of IAA, which is shorter than the T. gondii cell cycle. In marked contrast to the 24 hour treatment regime, SPARKEL- and SPARK-AID parasites exhibited no deficiencies in their ability to invade, egress, or mobilize calcium after 3 hours of depletion (Figure 2E– J). Thus, in all measured assays, mutants lacking SPARK and SPARKEL are phenotypically identical, likely due to reciprocal co-depletion, but such phenotypes develop over time.

SPARK and SPARKEL depletion leads to AGC kinase down-regulation and up-regulation of bradyzoite-stage proteins

Having verified the interaction between SPARK and SPARKEL, we next investigated pathways regulated by the complex. Initially, we performed a quantitative proteomics experiment following 24 hours of SPARK depletion—the previously described depletion window (Smith et al., 2022). Extended SPARK depletion led to down-regulation of PKG and both the catalytic and regulatory subunits of PKA (Figure 3A and Table 1). We also observed up-regulation of bradyzoite-specific genes upon SPARK depletion (Figure 3A and Table 1), including MAG1, BRP1, TGME49_208740, CST1, CST10, MIC13, and BFD2.

Figure 3.

To determine the critical window of SPARK and SPARKEL function, we performed quantitative proteomics with tandem mass tag multiplexing that included samples with 0, 3, 8, and 24 hours of SPARK or SPARKEL depletion. The experiments included internal TIR1 parental strain controls to account for basal downregulation arising from the degron (Figure 3B). We detected 3,333 and 3,880 proteins with quantification values in the SPARK and SPARKEL depletion time courses, respectively (Supplementary File 3). PKA C1, PKA R, and PKG protein abundances only began to drop after 8 hours of IAA treatment (Figure 3C). Globally, these kinase subunits were the most down-regulated proteins following 24 hours of either SPARK or SPARKEL depletion (Figure 3—figure supplement 1). Bradyzoite proteins up-regulated upon SPARK depletion primarily increased between 8 and 24 hours of IAA treatment (Figure 3D); the same proteins were initially down-regulated in SPARKEL-AID parasites relative to the parental strain and increased most in abundance before 24 hours of IAA treatment (Figure 3D). Thus, depletion of the SPARK-SPARKEL complex consistently leads to up-regulation of proteins associated with the bradyzoite stage of development, as well as down-regulation of PKG and the PKA C1 complex.

A depletion phosphoproteome implicates SPARK in the signaling, gene regulation, and metabolic states of the lytic cycle

To determine SPARK-dependent phosphoregulation, we enriched and analyzed phosphopeptides from the depletion proteome. We identified 18,518 phosphopeptides with quantification values, of which 8,867 were quantified with more than one peptide sequence match (PSM). Principal component analysis (PCA) separated the SPARK-AID samples treated for 24 h with IAA from all other samples, while closely clustering biological replicates (Figure 3—figure supplement 1). Despite the low abundance of SPARK and SPARKEL, individual phosphopeptides were detected for each protein. Likely due to basal degradation, the abundance of the SPARK peptide was ten-fold lower in SPARK-AID samples compared to the TIR1 control such that further depletion could not be measured (Figure 3— figure supplement 1). The SPARKEL peptide abundance was also reduced relative to the TIR1 control sample; however, it decreased further between 3 and 8 hours of IAA treatment (Figure 3—figure supplement 1). Because no other SPARKEL peptides were detected, it is not possible to determine whether this decrease is phosphosite-specific or due to a general decrease in SPARKEL protein abundance.

SPARK-regulated phosphopeptides may manifest in a strain-dependent or IAA-dependent manner. To account for basal downregulation, we ranked peptides by their log₂-ratios of altered abundances in the SPARK-AID strain relative to the TIR1 parental strain (Figure 3E). When a cutoff of 3.5 modified Z- scores was used, 91 peptides (44 with >1 PSMs) and 69 peptides (19 with >1 PSMs) were down- or up-regulated, respectively. The aforementioned phosphosite belonging to SPARK was the most down-regulated peptide by this metric. A phosphatidylinositol 3- and 4-kinase (PI3,4K) was also strongly down-regulated (Table 1); however, the protein abundance of this enzyme was not quantified. To identify peptides most altered by IAA treatment, and hence SPARK depletion, we summed and ranked peptides by log₂-ratios of abundances in the SPARK-AID strain relative to the untreated samples (Figure 3F). By this metric, 138 peptides (51 with >1 PSMs) and 93 peptides (36 with >1 PSMs) were down- or up-regulated by more than 3.5 modified Z-scores, respectively. Phosphoproteins belonging to this category are discussed in greater depth below.

To uncover the molecular pathways regulated by SPARK, we performed enrichment analysis of the peptides exhibiting the highest degree of altered abundance upon SPARK depletion (Figure 3G). Changes in SPARK function altered PTM cascades both directly and indirectly. SPARK depletion reduced the abundance of phosphopeptides belonging to cyclic nucleotide-regulated kinases, which was in part driven by a reduction in kinase abundance. SPARK depletion also coincided with a reduction in phosphosites on several phosphatases (Table 1), including BCC6, PPM1, and an uncharacterized dual-specificity phosphatase. By contrast, phosphorylation of other kinases and phosphatases increased upon SPARK depletion, which may arise as an indirect consequence of SPARK function. Upregulated phosphosites were found on pyruvate dehydrogenase and STE kinase, CDPK3, cyclin2-related protein, and the kinases TGGT1_204280 and TGGT1_239885; as well as the protein phosphatases PPM2B, calcineurin, and EFPP. Phosphosites belonging to proteins involved in ubiquitin transfer and hydrolysis also increased in abundance upon SPARK depletion. Examples of this class included UBP1, USP39, an RBR E3 ligase, CCR4-NOT subunit, and UBR box E3 ligase. To a lesser degree, phosphosites on ubiquitin ligases were also down-regulated following SPARK depletion—for example, a HECT E3 and putative E2 enzyme. Thus, SPARK function has large implications for PTM cascades.

Additional functional categories enriched for SPARK-dependent regulation included nucleic acid– binding proteins (Figure 3G and Table 1). A number of down-regulated phosphosites belong to proteins involved in DNA accessibility and transcription—for example, TFIIS, SET1, NOC3p, and TgSRCAP. RNA modification enzymes also exhibited phosphoregulation, such as two PUS proteins, PRP8, SLU7, a FHA protein, three KH proteins, and a DCP2 homolog. Other regulated phosphoproteins may function in mRNA cap or untranslated region (UTR) binding, such as TgPUF1, a MIF4G protein, CBP80, eIF3, and GCN1. Some of these proteins have been linked to expression changes during stage conversion, with PUF1 and SRCAP levels having been noted to increase during stress (Joyce et al., 2013; Liu et al., 2014; Sullivan et al., 2003). Moreover, CBP80, eIF3, and the KH protein TGGT1_216670 co-purified with the Alba complex (Gissot et al., 2013), which is necessary for bradyzoite differentiation. Regulation of these proteins thus may be linked to the elevated expression of bradyzoite-stage proteins upon SPARK depletion.

Several putative transporters also exhibited decreased phospho-abundance upon SPARK knockdown (Figure 3G). Some of these phosphoproteins channel inorganic cations, such as TgZnT (Chasen et al., 2019). Others likely transport small organic molecules—for example, a nucleoside and two ABC transporters. In aggregate, these analyses implicate SPARK function in signaling, gene regulation, and metabolic exchange.

Clustering of phosphopeptide kinetics identifies seven response signatures

Mammalian PDK1 phosphorylates numerous AGC kinases at activation loop threonines (Mora et al., 2004; Pearce et al., 2010), as well as non-kinase substrates. We reasoned that SPARK-dependent phosphopeptides depleted within 3 hours of IAA treatment may represent sites constitutively regulated by the kinase. By contrast, phosphosites down-regulated only after 8 hours of SPARK depletion—in excess of the length of the T. gondii cell cycle—may represent newly synthesized proteins dependent on SPARK for activation. We sought to resolve the phosphoproteome into kinetically distinct clusters. Gaussian mixture-model algorithms heuristically clustered peptides identified by more than one PSM into seven classes on the basis of their response signatures following SPARK depletion (Figure 3H). Two phosphopeptides in Cluster 1 were distinguished by their rapid and sustained depletion kinetics and corresponded to CBP80 and eIF3, both of which function in mRNA and translational regulation. The second cluster of phosphopeptides was largely down-regulated by 3 hours of IAA treatment and lacked dynamics following this window. PKA and PKG phosphopeptides belonged to Cluster 3, which decreased only after 8 hours of IAA treatment. The final cluster of depleted phosphopeptides, Cluster 4, only exhibits down-regulation at 8h of IAA treatment and may reflect phosphoregulation specific to the cell cycle.

Up-regulated phosphopeptides mirrored the dynamics of the down-regulated clusters. Cluster 5 phosphopeptides were elevated within 3 hours of IAA treatment and continued to increase thereafter. Cluster 6 was up-regulated by 3 hours of IAA treatment and lacked dynamics following this window. Cluster 7 increased gradually and predominantly after 8 hours of IAA treatment. Upregulated phosphosites may arise from a number of mechanisms following SPARK depletion. For example, SPARK activity may inhibit downstream kinases or activate downstream phosphatases; sterically block access of serines/threonines to other enzymes; or give rise to crosstalk between PTMs. Approaches that take into account subcellular proximity were required to distinguish direct substrates from pathways downstream of SPARK function.

Proximity labeling identifies putative SPARK targets in situ

To gain insight into proximal SPARK interactors in intracellular parasites, we tagged the endogenous C terminus of SPARK with a TurboID domain and carried out proximity labeling experiments, as previously described for SPARKEL. We observed overlap in the phosphoproteomics and proximity labeling experiments. PKG was enriched (Figure 3I), while PKA C1 and PI4P5K only narrowly missed our cutoffs. A putative HECT E3 ubiquitin ligase with two phosphosites down-regulated in the SPARK- dependent phosphoproteome was a candidate interactor. We also detected enrichment of PKA C3. PKA C3 is a coccidian-specific ortholog of the apicomplexan PKA catalytic subunit (Sugi et al., 2016); however, no PKA C3 peptides were quantified in our depletion phosphoproteome experiments, likely due to the low abundance of the protein. These integrated proteomics approaches hint that the phenotypes associated with SPARK depletion may arise from the altered activity of AGC kinases regulated by SPARK.

SPARK functionally interacts with PKA and PKG

Our proteomic data suggest SPARK controls the lytic cycle via activation of PKA C1 and PKG. SPARK depletion reduced levels of PKA C1 T190—which corresponds to the kinase activation loop—more rapidly than the corresponding protein levels (Figure 4A). Similarly, SPARK depletion triggered the down-regulation of several phosphosites residing in the N terminus of PKA R (e.g. S17, S27, and S94) more robustly than the changes in protein levels (Figure 4B). To confirm the depletion kinetics of PKA, we tagged the kinase subunits with fluorophores in the SPARK-AID parasite strain (Figure 4C). Prolonged SPARK depletion (24 hours) reduced PKA C1 and PKA R immunofluorescence signals in intracellular parasites (Figure 4D). To monitor kinase levels kinetically, we performed flow cytometry at several time points following the addition of IAA (Figure 4E). PKA C1 and PKA R fluorophore intensity decreased predominantly after 8 hours of IAA treatment, consistent with proteomic results. In parallel experiments, we tagged PKA C1 and PKA R in the SPARKEL-AID strain, revealing a comparable loss of the two markers following SPARKEL depletion (Figure 4F).

Figure 4.

We considered the possibility that the SPARK depletion phosphoproteome represents a convolution of phosphoproteomes arising from broader AGC kinase down-regulation. We systematically tested this possibility by generating phosphoproteomes of parasites depleted of PKA R and PKG (Supplementary File 4). As association with PKA R inhibits PKA C1 activity (Jia et al., 2017), we reasoned that phosphosites up-regulated after a short window of PKA R depletion represent candidate PKA C1 targets. Similarly, phosphosites down-regulated after PKG depletion represent candidate PKG substrates (Brown et al., 2017). Peptides above two modified Z-scores were considered candidate targets of each kinase.

Compared to the bulk phosphoproteome, the PKA C1-dependent phosphopeptides significantly decreased in abundance between 8 and 24 hours of SPARK depletion (Figure 4G). When we applied the same analysis to a similarly acquired PP1-dependent phosphoproteome (Herneisen et al., 2022), we found no significant relationship to the SPARK phosphoproteome (Figure 4–figure supplement 1), suggesting that the association with PKA C1 is specific. Several of the overlapping SPARK and PKA C1– dependent phosphoproteins localized to the parasite periphery and apex (Figure 4H), for example, IMC25, AC13, CIP1, CIP2, and AC3. Others are putatively involved in RNA-binding functions, such as eIF4E and two KH proteins, including a component of the METTL3/METTL14 core (Farhat et al., 2021). We also observed regulation of phosphoproteins involved in ubiquitin transfer, including the aforementioned E2 and HECT E3 proteins (Table 1). Regulation of these proteins may account for both the kinetic- and replicative-phase phenotypes attributed to perturbed PKA C1 activity (Jia et al., 2017; Uboldi et al., 2018).

We carried out related experiments to validate that SPARK regulates PKG function. SPARK depletion reduced levels of several PKG phosphopeptides more rapidly than the corresponding protein levels (Figure 5A). Notably, T838 corresponds to the PKG activation loop. Attempts to tag PKG with a fluorophore in a functional SPARK-AID background were unsuccessful. Global genetic screens have indicated that parasites with hypomorphic PKG—as can arise through endogenous tagging—are rapidly outcompeted by wildtype parasites (Fang et al., 2018; Smith et al., 2022). Therefore, we instead relied on chemical-genetic interactions to validate the relationship between PKG and SPARK. We reasoned that parasites with down-regulated PKG activity might be sensitized to the specific PKG inhibitor Compound 1 (Donald et al., 2002). Previous studies had shown that SPARK was dispensable for A23187- induced egress (Smith et al., 2022); however, PKG function is required for ionophore-induced egress (Brown et al., 2017; Lourido et al., 2012). We performed egress assays using parasites depleted of SPARK for 24 hours and treated with different concentrations of Compound 1. As anticipated, SPARK depletion sensitized parasites to A23187-induced egress inhibition with Compound 1 (Figure 5B–C), suggesting that PKG activity is decreased in the absence of SPARK.

Figure 5.

As with PKA, we quantified the behavior of putative PKG substrates in the SPARK-AID phosphoproteome. Compared to the bulk phosphoproteome, the PKG-dependent phosphopeptides significantly decreased in abundance between 8 and 24 hours of SPARK depletion (Figure 5D). Several overlapping SPARK and PKG-dependent phosphoproteins localized to the parasite periphery, including an MFS transporter, TgA1, MyoG, IMC18, IMC28, and IMC41. Other PKG/SPARK-dependent sites were found in the apical proteins AAP3, AC5/TLAP3, AC13, and SPM1; while several other proteins localized to the ER, including two TB2/DP1, HVA22 family proteins, a TPR protein, and an ABC transporter (Figure 5E). Others function in the secretory or protein trafficking pathways, such as Sec7, COG6, TgTBC10, and TgDHHC1. We also observed down-regulation of phosphosites belonging to proteins involved in ubiquitin transfer, such as an E2 enzyme and a ubiquitin family protein with a C-terminal extension. The precise pathways regulated by TgPKG remain to be defined; however, recent studies triggering PKG activity in T. gondii (Chan et al., 2023; Herneisen et al., 2022) or related Plasmodium spp. (Alam et al., 2015; Balestra et al., 2021; Brochet et al., 2014) suggested that PKG substrates relevant for egress reside at the parasite periphery and ER, as observed for the overlapping PKG and SPARK targets here.

Depletion of SPARK, SPARKEL, or PKA C3 promotes chronic differentiation

Several pieces of evidence suggested that SPARK and SPARKEL depletion may prompt differentiation from the tachyzoite to bradyzoite state through dysregulation of another AGC kinase, PKA C3. First, down-regulation of SPARK and SPARKEL coincided with up-regulation of bradyzoite markers (Figure 3A and 3D). Second, SPARK proximity labeling resulted in enrichment of PKA C3 (Figure 3I)—a negative regulator of the bradyzoite stage (Sugi et al., 2016) and, to date, the only AGC kinase characterized with a function in differentiation (Augusto et al., 2020). Prolonged SPARK and SPARKEL depletion indeed elevated differentiation under normal culture conditions, as measured by staining of the cyst walls with dolichos binding lectin, a characteristic of chronic-stage bradyzoite vacuoles (Figure 6A–B). N-terminal tagging of SPARKEL with AID similarly elevated differentiation upon IAA treatment (Figure 6–figure supplement 1). Alkaline media in combination with serum starvation and low carbon dioxide is commonly used to elevate bradyzoite transformation in cell culture, and SPARK and SPARKEL down-regulation enhanced differentiation under such conditions (Figure 6A and 6B).

Figure 6.

To dissect the hypothesized PKA C3–dependent arm of SPARK regulation, we altered the PKA C3 genomic locus with a V5-mAID-mNG-Ty tagging payload. Addition of auxin to the culture medium resulted in down-regulation of PKA C3 levels to below the detection limit within one hour (Figure 6C– D). PKA C3 exhibited low expression levels in cytoplasmic puncta (Figure 6C), as previously reported (Sugi et al., 2016). We confirmed that prolonged PKA C3 depletion led to elevated differentiation rates under normal and alkaline-stress culture conditions (Figure 6A–B) and observed that the effect was similar to SPARK or SPARKEL knock-down. Given that SPARK also regulates PKG and PKA C1, we considered the possibility that differentiation may be linked to a block in motility. However, the inability to egress is not sufficient for transformation, as conditional knockdown of PKG or CDPK1 (Brown et al., 2017; Chan et al., 2023; Shortt et al., 2022) failed to significantly enhance the frequency of differentiation markers (Figure 6A–B). We did not assess differentiation in PKA C1 knockdown parasites, as these mutants spontaneously egress from host cells (Jia et al., 2017; Uboldi et al., 2018).

The RH strain in which we generated mutants is reported to be refractory to differentiation into the bradyzoite form (Dubey et al., 1999), like due to defects downstream of the pathways that lead to the initial synthesis of the cyst wall. We generated SPARKEL- and PKA C3-AID strains in the ME49 background (Figure 6–figure supplement 1), which readily undergoes differentiation in tissue culture and mice, and confirmed that depletion prompts spontaneous differentiation in tissue culture (Figure 6E–F). We did not assess alkaline-induced differentiation, as this stress treatment results in near-complete differentiation in the ME49 strain (Waldman et al., 2020), and thus does not provide a differential for mutants that enhance differentiation.

To assess whether differentiation in the RH-type strain is dependent on the master regulator of differentiation, BFD1 (Waldman et al., 2020), we knocked out BFD1 CDS with a sortable dTomato cassette in the SPARK- and PKA C3-AID strains (Figure 6–figure supplement 1). The resulting Δbfd1 mutants failed to undergo differentiation as measured by cyst wall staining (Figure 6G–H). We additionally profiled transcriptomic remodeling in PKA C3-AID parasites that had been depleted of the kinase for 24 hours and compared these changes to the chronic-stage transcriptome of parasites overexpressing BFD1 (Figure 6I and (Waldman et al., 2020)). We observed up-regulation of several chronic-stage proteins upon PKA C3 depletion, including the bradyzoite antigen BAG1, CST1, CST10, MIC13; and TgSPT2, which functions in ceramide synthesis and was recently reported to exhibit BFD2- dependent abundance (Licon et al., 2023). PKA C3 mRNA abundance was down-regulated, suggesting that the tagging payload may alter steady-state transcript levels. Together, these results suggest that the bradyzoite stage conversion induced by knockdown of PKA C3 and SPARK proceeds through the canonical transcriptional pathway.

PKA C3 interacts with the SPARK complex

Several AGC kinases, including PKA C3, share the activation loop motif TLC/VGTxxY, which displayed SPARK-dependent phosphorylation in PKA C1 and PKG. PKA C3 was not detected in our SPARK depletion proteomics experiments (Figure 3). To determine whether SPARK and PKA C3 interact, we immunopurified PKA C3 lysates, using the mNG epitope as a handle. The immunoprecipitated PKA C3 was highly enriched for SPARK (Figure 7A), and to a lesser extent other protein kinases, including TGGT1_205550, SRPK, and AMPK subunit beta. The latter has recently been characterized as a safeguard of parasite cell cycle regulation that enhances differentiation when knocked out (Yang et al., 2022). To a lesser extent, SPARKEL also significantly co-purified with PKA C3 (Supplementary File 1). Our results are consistent with other reports that PKA C3 does not interact with the canonical PKA regulatory subunit in T. gondii (Jia et al., 2017; Uboldi et al., 2018) and, in conjunction with TurboID data presented in Figure 3I, provide support for a strong physical interaction between SPARK and PKA C3.

Figure 7.

SPARK down-regulation does not reduce PKA C3 protein levels

To determine whether SPARK activity affects PKA C3 abundance, as for PKA C1 and PKG, we tagged the endogenous PKA C3 C terminus with an mNG fluorophore in the SPARK-AID background. After 24 hours of IAA treatment, PKA C3 levels were unaltered or slightly elevated, as measured by immunoblot (Figure 7B–C). Immunofluorescence and flow cytometry measuring mNG intensity after 0, 3, 8, and 24 hours of IAA treatment revealed a trend towards increasing PKA C3 abundance, although this difference was not significant (Figure 7D–F). Thus, unlike PKA C1, PKA C3 does not exhibit down-regulation following SPARK depletion. We considered the possibility that SPARK function could instead alter PKA C3 activity levels.

The PKA C3 phosphoproteome identifies candidates regulating the transition to the bradyzoite state

To pursue the functional consequences of PKA C3 inhibition, we performed a depletion phosphoproteomics time course experiment with the PKA C3-AID strain following 1, 3, 8, and 24 hours of IAA treatment. PKA C3-dependent phosphopeptides were identified as for the SPARK depletion proteome (Figure 7G). The resulting proteome identified 4,792 proteins and 12,721 phosphopeptides, of which 9,576 were quantified with more than one PSM. PCA separated the 24 hour IAA samples from all others (Figure 7–figure supplement 1). Similar to the SPARK depletion proteome, several of the up-regulated proteins were reported to be transcriptional signatures of the bradyzoite stage, including BRP1, TGGT1_208740, two secreted cAMP-dependent protein kinases, and BFD2 (Figure 7–figure supplement 1).

To identify peptides most altered by IAA treatment, we summed and ranked peptides by log₂-ratios of abundances in the PKA C3-AID strain relative to the untreated (0 h) samples (Figure 7G). By this metric, 129 peptides (94 with >1 PSMs) and 49 peptides (31 with >1 PSMs) were down- or up-regulated by more than 3.5 in the modified Z-scores, respectively. Gaussian mixture-model algorithms heuristically clustered peptides identified by more than one PSM into three classes (Figure 7H). Cluster 1 contained 12 peptides that rapidly decreased within 1–3 hours of PKA C3 depletion, and continued to decrease thereafter. This cluster includes CBP80, a DEP domain protein, TgZnT, TgApiAT3-1, AC13, and FtsH1 (Table 1).

Cluster 2 comprised 82 peptides that also decreased within 1–3 hours of IAA treatment, after which their abundances remained stable and decreased overall to a lesser extent than cluster 1 (Figure 7H). Cluster 2 included genes putatively involved in gene regulation, including eIF4G1, NOC3p, a zinc finger protein, a DNA repair protein, SET1, and a nucleotidyltransferase (Table 1). Several phosphosites belonging to proteins involved in daughter cell biogenesis were also down-regulated, including AC13, IAP2, BCC8, and condensin 2. We also observed down-regulation of transporters belonging to this cluster–for example, TgApiAT2, ATP4, and a nucleoside transporter. The early events in bradyzoite differentiation have not been extensively characterized but involve translational regulation of BFD1 enhanced by a BFD1/BFD2 feed-forward loop (Licon et al., 2023; Waldman et al., 2020). Down-regulation of PKA C3 may enhance conditions that activate this loop, for example due to global down-regulation of cap-dependent translation (CBP80, eIF4G1) or alteration of the replicative cycle (Radke et al., 2003).

Cluster 2 also contains numerous proteins trafficking through secretory organelles, including the micronemes, rhoptries, and dense granules. Modified phosphosites belonging to micronemal proteins included CRMP-b, TGGT1_221180, and TGGT1_304490 (Table 1). The rhoptry proteins ROP1, ROP13, ROP17, ROP40, and RON9 were also differentially regulated. Modified dense granule/parasitophorous vacuole proteins included GRA3, GRA4, GRA6, GRA7, GRA8, GRA31, GRA57, GRA62, GRA70, SFP1, and MAG1. Many of these proteins are lumenal within their respective organelles, and regulation by PKA C3 is therefore likely indirect. Nonetheless, each of the secretory organelles contain stage-specific subproteomes with implications for metabolism and protein complex assembly in bradyzoites (Sinai et al., 2020); the regulation observed here may reflect such early stages of proteomic remodeling.

The 31 phosphopeptides belonging to cluster 3 increased in abundance following PKA C3 depletion and therefore are likely indirect targets of the kinase. This cluster contained phosphoproteins involved in gene regulation (Table 1), such as AP2XII-9, an ISWI protein, CAF1, PLU-1, GCFC, and two zinc finger proteins: TGGT1_223880 and BFD2. We also identified several candidates involved in proteostasis, including HSP90, p23, CSN3, PSME4. Regulation of such proteins may be related to the transition and maintenance of the bradyzoite state, albeit not as a direct consequence of PKA activity.

The protein phosphatase PP2A has recently been characterized as a regulator of bradyzoite differentiation (Wang et al., 2022). We therefore analyzed proteins that exhibited both PKA C3- and PP2A-dependent phosphoregulation. Several candidates in this category were identified as differentially important in interferon-stimulated cells, including TGGT1_209500, GRA57, and GRA70. Numerous phosphosites belonging to AC13 were differentially regulated in both phosphoproteomes. Some proteins involved in endo- and exocytosis also exhibited dynamic phosphoregulation, including CRMP-b and Kelch13. PKA R S²⁷ was up-regulated upon PKA C3 depletion and was also up-regulated in Δpp2a mutants. Finally, nucleic acid-binding proteins exhibited differential phosphoregulation, including BFD2, a KH protein, eIF4G2, and AP2XII-1. Given the stage conversion phenotypes of parasites lacking PKA C3 or PP2A, the phosphoregulation of overlapping targets may be functionally important in the context of bradyzoite differentiation.

The SPARK phosphoproteome shares signatures with the PKA C3–depletion phosphoproteome

We anecdotally observed overlap between SPARK- and PKA C3-dependent phosphopeptides, and in particular those with putative nucleic acid binding functions. We therefore systematically compared the abundances of candidate PKA C3 targets (as defined above) to overlapping phosphosites in the SPARK phosphoproteome. Compared to the bulk phosphoproteome, the PKA C3–dependent phosphopeptides significantly decreased in abundance within three hours of SPARK depletion (Figure 7I). Several of phosphoproteins found to be altered by both PKA C3 and SPARK depletion belonged to proteins involved in gene regulation (Figure 7J and Table 1), including CBP80, eIF4G2, SET1, NOC3p, Spt5, a KH protein, a nucleotidyltransferase, and PRP4. Some of these proteins associate with large transcription factor assemblies, including the AP2IX4/MORC and GCN5-B complexes. Transporters also overlapped between the two phosphoproteomes, including ZnT, TgApiAT2, and TgApiAT3-1. Several phosphosites belonging to apical cap proteins were also shared between the phosphoproteomes, including numerous AC13 peptides.

The SPARK and PKA C3 phosphoproteomes notably diverged with regards to regulation of dense granule proteins, which was observed as a subpopulation of PKA C3–dependent peptides that increased in abundance upon SPARK depletion (Figure 7I and 7K). This effect was specific to the indicated phosphopeptides, as corresponding protein levels of the dense granule proteins were not changed (Supplementary File 4). The differentially regulated dense granule proteins include many of those belonging to cluster 2 of the PKA C3 depletion proteome (Figure 7H and 7K). The phosphoproteome of SPARK shares signatures of reduced PKA C3 activity; however, due to concomitant down-regulation of other AGC kinases, the phosphoproteomes of SPARK and PKA C3 depletion diverge in ways that remain to be explored.

Discussion

PDK1, the closest mammalian homolog to SPARK, is considered a key regulator of AGC kinases (Mora et al., 2004; Pearce et al., 2010). These kinases are united by general structural features required for activity and integrity, such as the activation segment (also known as the activation loop or T-loop), and the C-terminal hydrophobic and turn motifs, which are common loci of regulation that can position the kinase into an active conformation. In mammals, AGC kinases depend on PDK1 in various ways. For example, phosphorylation of the Akt (PKB) activation loop by PDK1 leads to partial activation; additional phosphorylation by the target of rapamycin (mTOR) complex produces maximal activity (Pearce et al., 2010). Some protein kinase C (PKC) isoforms undergo ordered phosphorylation by upstream enzymes—including PDK1—that prime PKC for activity and prevent kinase degradation (Newton, 2010). PKA catalytic isoforms require T-loop phosphorylation for activity; however, whether this modification is installed autocatalytically or by PDK1 depends on the cellular context (Cheng et al., 1998; Moore et al., 2002; Williams et al., 2000). Consequently, PDK1 dysregulation provokes diverse cellular consequences, including enhanced growth, migration, and cell survival—processes associated with cancer in animals (Gagliardi et al., 2018; Mora et al., 2004; Raimondi and Falasca, 2011). In this study, we characterize the diverse consequences of PDK1 dysregulation in the ubiquitous and divergent eukaryotic parasite of such animals, T. gondii.

SPARK down-regulation coincides with reduced PKG and PKA C1 activity due to reduced protein levels. From our own studies, multiple mass spectrometric and cell biological experiments support this observation. Furthermore, PDK1 mutations in the related apicomplexan Plasmodium falciparum were observed to act as suppressors of toxicity related to overexpression of PfPKA (Hitz et al., 2021). Both TgPKG and TgPKA C1 possess features of AGC kinases involved in PDK1 interaction, such as the PDK1- interacting fragment FXXF. PKA C1 terminates in the sequence FTSW, whereas PKG has two non-terminal FXXF motifs, FGDF and FLYF. Although these motifs are not strictly C-terminal, as in animals, their position resembles that observed in some plant AGC kinases, up to 50 amino acids upstream of the C-terminus (Rademacher and Offringa, 2012). The T-loop phosphorylation of each kinase—PKA C1 T190 and PKG T838—are down-regulated to a greater extent than corresponding protein levels following SPARK knockdown (Supplementary File 3). Heterologous expression and activity assays of PKA C1 and PKG may resolve the importance of activation loop phosphorylation for kinase function; however, such biochemical assays are not currently available for the T. gondii kinases. The requirement of T-loop threonine for PKG activity is unclear, as recombinant Plasmodium spp. PKG was still active when this residue was mutated (Bakkouri et al., 2019). Conformational changes due to cGMP binding— in particular, rearrangement of the N-terminal autoinhibitory segment otherwise adjacent to the activation loop—appeared particularly important for PKG activity. Upon SPARK depletion, multiple phosphorylation sites in the N-terminal and cyclic nucleotide-binding domains of PKG (S20, S22, T59, S62, T96, T99, and S105) similarly decrease in abundance to a greater extent than protein levels, and may function in the regulation of these domains. The cause for reduction of PKA C1 and PKG levels requires further study; in principle, active degradation or dilution mechanisms due to cell division are possible. In support of the latter, down-regulation of both kinases was only observed after 8 hours of SPARK depletion, a period in excess of the duration of the T. gondii cell cycle.

Independent of biochemical mechanisms, genetic evidence supports the hypothesis that SPARK phenotypes associated with kinetic transitions in the lytic cycle—such as motility, egress, and invasion—arise from dysregulated PKA C1 and PKG activity. In T. gondii, PKA negatively regulates egress and positively regulates the transition between invasion and the establishment of a replicative niche (Jia et al., 2017; Uboldi et al., 2018); by contrast, in P. falciparum, PKA mediates parasite invasion (Flueck et al., 2019; Patel et al., 2019; Wilde et al., 2019). PKG is more broadly required for apicomplexan motility, invasion, and egress and functions by mobilizing parasite intracellular calcium stores. Genetic pathway analysis of T. gondii PKA C1 and PKG inhibition supports a model in which PKG lies downstream of PKA (Jia et al., 2017). Increasingly, proteomic, genetic, and physiological studies implicate negative feedback mechanisms between cAMP and cGMP signaling pathways in apicomplexans (Alam et al., 2015; Bisio et al., 2019; Brochet et al., 2014; Moss et al., 2022; Nofal et al., 2022). In T. gondii, SPARK inhibition releases the breaks that PKA C1 holds on egress and jams the engine of PKG activity, resulting in parasite stasis. The PKA signaling circuit differs in Plasmodium spp., where a sole PKA C subunit is indispensable for invasion but not egress (Patel et al., 2019; Wilde et al., 2019). However, SPARK inhibition would lead to a similarly inhibitory effect on parasite proliferation, as has been suggested recently in studies employing altered P. falciparum PDK1 and PKA alleles (Hitz et al., 2021).

1. T. gondii has two additional paralogs of the PKA catalytic subunit, relative to Plasmodium spp., PKA C2 and PKA C3. PKA C2 expression levels are low outside of the parasite sexual stages (Ramakrishnan et al., 2019; Sugi et al., 2016). PKA C3 disruption was previously reported to increase conversion to the bradyzoite stage (Sugi et al., 2016). Our results using a PKA C3 conditional knockdown support this finding. Furthermore, PKA C3 and SPARK associate, as revealed through approaches that probe reciprocal interactions. PKA C3 possesses a canonical PDK1-interacting fragment (FDNF) and shares an activation loop motif with PKA C1. Although PKA C3 peptide levels were below the detection limit in our depletion phosphoproteome experiments, we propose that SPARK depletion leads to down-regulation of PKA C3 activity separately from a reduction in PKA C3 protein abundance. SPARK and PKA C3 depletion phenocopy each other with respect to in vitro differentiation. Candidates for the molecular targets responsible for differentiation may reside in the overlapping phosphosites in the SPARK and PKA C3 depletion phosphoproteomes.

The regulation of SPARK remains an open inquiry. Mammalian PDK1 is constitutively active due to auto-phosphorylation in trans. The PH domain triggers mammalian PDK1 homodimerization upon binding of phosphatidylinositol (3,4,5)-trisphosphate (Levina et al., 2022; Masters et al., 2010; Pearce et al., 2010). SPARK lacks the PH domain (Smith et al., 2022); therefore, apicomplexans have likely evolved alternative mechanisms for regulating PDK1. The complex of proteins interacting with SPARK—the elongin-like protein SPARKEL and the putative AGC kinase TGME49_205550—may represent a means of modulating SPARK activity. Elongin C forms a component of multiprotein complexes. E3 ubiquitin ligases are well-known examples; however, elongin C can affect steady-state protein levels through mechanisms distinct from ubiquitination (Hyman et al., 2002). In metazoans, PDK1 has been reported to associate with E3 ligases and exhibit monobiquitination (Jiang et al., 2021; Uras et al., 2012). The anti-correlated cell cycle expression profiles of SPARK and SPARKEL suggest a specific window of activity for the complex. However, given the co-depletion of the complex—as has been reported for other tightly associated proteins studied with the AID system (Tosetti et al., 2020)— further biochemical studies are required to discern the regulatory interactions between SPARK and SPARKEL. Plasmodium spp. appear to have lost SPARKEL, which is otherwise conserved in the apicomplexans, possibly requiring tight regulation of SPARK expression or other mechanism to control its activity.

Our unbiased and integrative proteomic approaches uncovered ancient features of signaling networks de novo and revealed how such relationships can be repurposed for parasitism. Apicomplexans are well-adapted to transition between life stages in their hosts. Synthesizing genetic, proteomic, and phenotypic data, we propose that SPARK fundamentally integrates the topology of cyclic nucleotide signaling networks in apicomplexan parasites. In T. gondii, proper SPARK function preserves tachyzoite physiology by activating PKG, PKA C1, and PKA C3, thus promoting cycling within the acute stages of the infection in favor of differentiation into the chronic stage wherein AGC kinase activity is lower (Fu et al., 2021). In principle, SPARK may serve as a central node for regulation of the asexual stages.

Key resources table

Materials & methods

Reagents and cell culture

1. T. gondii parasites of the type I RH strain Δku80/Δhxgprt genetic background (ATCC PRA-319, (Huynh and Carruthers, 2009)) or the type II ME49 strain Δku80/Δhxgprt background (Waldman et al., 2020) were grown in human foreskin fibroblasts (HFFs, ATCC SCRC-1041) maintained in DMEM (GIBCO) supplemented with 3% calf serum and 10 μg/mL gentamicin (Thermo Fisher Scientific). When noted, DMEM was supplemented with 10% inactivated fetal calf serum and 10 μg/mL gentamicin. HFFs and parasites were tested routinely for mycoplasma using the ATCC Universal Mycoplasma Detection Kit (30-1012K).

Oligonucleotide sequences are listed in Supplementary File 5.

Parasite transfection and strain construction

Genetic background of established parasite strains

This study used the following strains that have been characterized in other publications: RHΔku80Δhxgprt/TIR1 (Brown et al., 2017), ME49Δku80Δhxgprt/TIR1 (Licon et al., 2023) RHΔku80Δhxgprt/TIR1/pTUB1-GCaMP6f (Smith et al., 2022), RHΔku80Δhxgprt/TIR1/pMIC2-MIC2- Gluc-myc-P2A-GCaMP6f (Herneisen et al., 2022), PKG-mAID-HA/RHΔku80/TIR1 (Brown et al., 2017), CDPK1-mNG-mAID-Ty/RHΔku80Δhxgprt/TIR1 (Shortt et al., 2022; Smith et al., 2022), SPARK-mNG- mAID-Ty/RHΔku80Δhxgprt/TIR1 (Smith et al., 2022), SPARK-mCherry-mAID- HA/RHΔku80Δhxgprt/TIR1/pTUB1-GCaMP6f (Smith et al., 2022), and pMIC2-mNG-TurboID- Ty/RHΔku80Δhxgprt/TIR1 (Chan et al., 2023). The strains described below were derived from these parental lines.

SPARK-V5-mNG-Ty/RHΔku80Δhxgprt

The HiT vector cutting unit gBlock described for SPARK in (Smith et al., 2022) was cloned into the pALH193 HiT empty vector, which is based on the pGL015 cloning vector (GenBank: OM640005) but possesses the XTEN-V5-mNG-Ty-3’CDPK3/DHFR tagging unit. The vector was linearized with BsaI and co-transfected with the pSS014 Cas9 expression plasmid (GenBank: OM640002) into RHΔku80Δhxgprt parasites. Clones were selected with 1 µM pyrimethamine and isolated via limiting dilution. Clones were verified by immunofluorescence microscopy and immunoprecipitation and mass spectrometry.

SPARKEL-V5-mNG-Ty/SPARK-V5-mCherry-HA/RHΔku80Δhxgprt/TIR1

The HiT vector cutting unit gBlock for SPARKEL (P1) was cloned into the pALH193 HiT empty vector. The vector was linearized with BsaI and co-transfected with the pSS014 Cas9 expression plasmid into RHΔku80Δhxgprt/TIR1 parasites. Clones were selected with 1 µM pyrimethamine and isolated via limiting dilution to generate the SPARKEL-V5-mNG-Ty/RHΔku80Δhxgprt/TIR1 strain. The HiT vector cutting unit gBlock described for SPARK in (Smith et al., 2022) was cloned into the pALH142 HiT empty vector, which is based on the pALH052 cloning vector (GenBank: OM863784) but possesses the XTEN- V5-mCherry-HA-3’CDPK3/HXGPRT tagging unit. The BsaI-linearized template was co-transfected with pSS014 into the SPARKEL-V5-mNG-Ty/RHΔku80Δhxgprt/TIR1 strain. Clones were selected with 25 µg ml⁻¹ mycophenolic acid and 50 µg ml⁻¹ xanthine and isolated via limiting dilution. Clones were verified by immunofluorescence microscopy.

SPARKEL-V5-HaloTag-mAID-Ty/RHΔku80Δhxgprt/TIR1, SPARKEL-V5-HaloTag-mAID- Ty/ME49Δku80Δhxgprt/TIR1, SPARKEL-V5-HaloTag-mAID-Ty/RHΔku80Δhxgprt/TIR1/pMIC2-MIC2-Gluc-myc-P2A-GCaMP6f and SPARKEL-V5-HaloTag-mAID-Ty/SPARK-V5-mCherry-HA/RHΔku80Δhxgprt/TIR1

The HiT vector cutting unit gBlock for SPARKEL was cloned into the pALH076 HiT empty vector, which is based on the pGL015 cloning vector (GenBank: OM640005) but possesses the XTEN-V5-HaloTag-mAID-Ty-3’CDPK3/DHFR tagging unit. The vector was linearized with BsaI and co-transfected with the pSS014 Cas9 expression plasmid into RHΔku80Δhxgprt/TIR1, RHΔku80Δhxgprt/TIR1/pMIC2-MIC2- Gluc-myc-P2A-GCaMP6f, or ME49Δku80Δhxgprt/TIR1 parasites. Clones were selected with 1 µM pyrimethamine and isolated via limiting dilution to generate the SPARKEL-V5-HaloTag-mAID- Ty/RHΔku80Δhxgprt/TIR1 strain. The linearized SPARK-V5-mCherry-HA HiT vector described above was co-transfected with pSS014 into the SPARKEL-V5-HaloTag-mAID-Ty/RHΔku80Δhxgprt/TIR1 strain. Clones were selected with 25 µg ml⁻¹ mycophenolic acid and 50 µg ml⁻¹ xanthine and isolated via limiting dilution. Clones were verified by PCR amplification and sequencing of the junction between the 3′ end of SPARKEL (5’-GGGAGGCCACAACGGCGC-3’) and 5′ end of the protein tag (5’- gggggtcggtcatgttacgt-3’), immunoblotting, and immunoprecipitation and mass spectrometry.

SPARK-V5-mAID-HA and PKA R-V5-mAID-HA/RHΔku80Δhxgprt/TIR1

HiT cutting units previously described for SPARK (Smith et al., 2022) and PKA R (P2) were cloned into HiT empty vector pALH086 (XTEN-V5-mAID-HA-3’CDPK3/HXGPRT; GenBank ON312869). In the case of the PKA R vector, the CDPK3 3’UTR was exchanged for the sequence 1092 base pairs downstream of the PKA R stop codon. The vectors were linearized and co-transfected with the Cas9 plasmid pSS014 into the RHΔku80Δhxgprt/TIR1 strain. Clones were selected with 25 µg ml⁻¹ mycophenolic acid and 50 µg ml⁻¹ xanthine and isolated via limiting dilution. In-frame tagging of PKA R was notably rare. Clones were verified by PCR amplification and sequencing of the junction between the 5′ end of the protein tag (5’-gggggtcggtcatgttacgt-3’) and 3′ end of SPARK (5’-GACGCAAAACTGGCAAGACG-3’) or the 5’UTR of the DHFR selection cassette (5’-tcgacaacgaatgacacaca-3’) and 3’ end of PKA R (5’- GTGCGCATCTTGGAAGACATGGATC-3’).

3HA-mAID-SPARKEL/RHΔku80Δhxgprt/TIR1

An N-terminal HiT cutting unit designed for SPARKEL (P3) was cloned into HiT empty vector pALH460, which is based on the pTUB1-DD HiT vector described in (Licon et al., 2023) but with the 5’DHFR-DHFR- T2A-3HA-mAID-XTEN tagging unit. The vector was linearized and co-transfected with the Cas9 plasmid pSS014 into the RHΔku80Δhxgprt/TIR1 strain. Clones were selected with 1 µm pyrimethamine and isolated via limiting dilution. Clones were verified for in-frame tagging via immunofluorescence microscopy.

PKA C1-mNG-3myc/PKA R-mCherry-Ty/RHΔku80Δhxgprt/TIR1 and SPARKEL-V5-HaloTag-mAID-Ty/PKA C1-mNG-3myc/PKA R-mCherry-Ty

The HiT cutting unit for PKA C1 (P4) was cloned into the HiT empty vector pALH4052, which is based on the pGL015 cloning vector (GenBank: OM640005) but possesses the XTEN-mNG-3myc-3’CDPK3 tagging unit. The HiT cutting unit for PKA R (P3) was cloned into the HiT empty vector pALH4082, which is based on the pGL015 cloning vector (GenBank: OM640005) but possesses the XTEN-mCherry-Ty-3’CDPK3 tagging unit. The vectors were linearized and co-transfected with the Cas9 plasmid pSS014 into the RHΔku80Δhxgprt/TIR1 strain. Clones were sorted on the basis of mNG and mCherry fluorescence and were isolated via limiting dilution. Clones were verified by live microscopy and flow cytometry. The linearized SPARKEL-V5-HaloTag-mAID-Ty HiT vector described above was transfected into the PKA C1-mNG-3myc/PKA R-mCherry-Ty/RHΔku80Δhxgprt/TIR1 strain with pSS014. Clones were selected with 1 µM pyrimethamine and isolated via limiting dilution.

SPARK-V5-mAID-HA/PKA C1-mNG-3myc/PKA R-V5-mCherry-2A-DHFR

The PKA C1-mNG-3myc HiT vector described above was co-transfected with pSS014 into the SPARK- V5-mAID-HA/RHΔku80/TIR1 strain. Clones were sorted on the basis of mNG fluorescence and were isolated via limiting dilution. Clones were verified by live microscopy. The HiT cutting unit described for PKA R above was cloned into the HiT empty vector pALH5292, which is based on the pGL015 cloning vector (GenBank: OM640005) but possesses the XTEN-V5-mCherry-P2A-DHFR tagging unit. This vector was linearized with BsaI and was transfected with pSS014 into the SPARK-V5-mAID-HA/PKA C1- mNG-3myc/RHΔku80/TIR1 strain. Clones were selected with 1 µM pyrimethamine and isolated via limiting dilution. Clones were verified by live microscopy.

PKA C3-V5-mNG-mAID-Ty/RHΔku80Δhxgprt/TIR1 and PKA C3-V5-mNG-mAID-Ty/ME49Δku80Δhxgprt/TIR1

The HiT vector cutting unit gBlock for PKA C3 (P5) was cloned into the pGL015 HiT empty vector (XTEN-V5-mNG-mAID-Ty-3’CDPK3/DHFR; GenBank OM640005). The vector was linearized with BsaI and co-transfected with the pSS014 Cas9 expression plasmid into RHΔku80Δhxgprt/TIR1 or ME49Δku80Δhxgprt/TIR1 parasites. Clones were selected with 1 µM pyrimethamine and isolated via limiting dilution. Clones were verified by PCR amplification and sequencing of the junction between the 3′ end of PKA C3 (5’-CGGGGTCATGGGCTACCTG-3’) and 5′ end of the protein tag (5’- gggggtcggtcatgttacgt-3’), immunofluorescence microscopy, immunoblotting, and immunoprecipitation and mass spectrometry.

PKA C3-V5-mNG-mAID-Ty/Δbfd1::dTomato/RHΔku80Δhxgprt/TIR1 and SPARK-V5-mAID- HA/Δbfd1::dTomato/RHΔku80/TIR1

The sequences 5’GGGGGGCTTGATGTAACAGA-3’ and 5’-cacacacttaagtacgggga-3’, which target upstream of the BFD1 transcription start site (Markus et al., 2020; Waldman et al., 2020) and downstream of its last exon, respectively, were cloned as small guide RNAs into the Cas9 plasmid pSS013 pU6-Universal (Sidik et al., 2016a). A repair template sequence with 40 bases of homology upstream to each of these guides was amplified from a pTUB1-dTomato-3’DHFR plasmid using the oligos 5’-caacctgcaggggtcacctctacagtgtttcgcaccatccgtgtcatgtagcctgccaga-3’ and 5’- cgactctacacggggaggaaggacgtcaacagaccctcttCATGCATGTCCCGCGTTCGT-3’. This repair template was co-transfected with the two sgRNA plasmids described above into PKA C3-V5-mNG-mAID- Ty/RHΔku80Δhxgprt/TIR1 and SPARK-V5-mAID-HA/RHΔku80Δhxgprt/TIR1. Clones were sorted on the basis of dTomato fluorescence and were isolated via limiting dilution. Clones were verified by live microscopy and PCR amplification of the junctions spanning upstream of the BFD1 transcription start site and in the first exon of BFD1 or the 3’ sequence of dTomato. The oligos used for this amplification are included in Supplementary File 5.

PKA C3-mNG-P2A-DHFR/SPARK-V5-mAID-HA/RHΔku80/TIR1

The HiT vector cutting unit gBlock for PKA C3 described above was cloned into the pALH459 HiT empty vector, which is based on the pGL015 cloning vector (GenBank: OM640005) but possesses the XTEN- mNG-P2A-DHFR tagging unit. The vector was linearized with BsaI and co-transfected with the pSS014 Cas9 expression plasmid into SPARK-V5-mAID-HA/RHΔku80/TIR1 parasites (described above). Clones were selected with 1 µM pyrimethamine and isolated via limiting dilution. Clones were verified by PCR amplification and sequencing of the junction between the 3′ end of PKA C3 (5’- CGGGGTCATGGGCTACCTG-3’) and 3′ end of the DHFR selection cassette (5’-gggcagcttctgtatttccg-3’), immunofluorescence microscopy, immunoblotting, and flow cytometry.

SPARK-TurboID-Ty and SPARKEL-TurboID-Ty/RHΔku80Δhxgprt/TIR1

The HiT vector cutting unit gBlocks for SPARK and SPARKEL described above were cloned into the pALH173 HiT empty vector with the XTEN-TurboID-Ty-3’CDPK3/DHFR tagging payload (Chan et al., 2023). The vectors were linearized with BsaI and co-transfected with the pSS014 Cas9 expression plasmid into RHΔku80Δhxgprt/TIR1 parasites. Clones were selected with 1 µM pyrimethamine and isolated via limiting dilution. Clones were verified by PCR amplification and sequencing of the junction between the 3′ end of SPARK and SPARKEL (using oligos described above) and the 5′ end of the protein tag (5’-gggggtcggtcatgttacgt-3’) as well as mass spectrometry.

3HA-TurboID-SPARKEL/RHΔku80Δhxgprt/TIR1

An N-terminal HiT cutting unit designed for SPARKEL (P3) was cloned into HiT empty vector pALH461, which is based on the pTUB1-DD HiT vector described in (Licon et al., 2023) but with the 5’DHFR- DHFR-T2A-3HA-TurboID-XTEN tagging unit. The vector was linearized and co-transfected with the Cas9 plasmid pSS014 into the RHΔku80Δhxgprt/TIR1 strain. Parasites selected with 1 µm pyrimethamine. Integration of the construct was verified by immunofluorescence microscopy.

Immunoprecipitation and mass spectrometry

SPARK-mNG-mAID, SPARK-mNG, and PKA C3-mNG IPs

Approximately 2 ✕ 10⁸ extracellular SPARK-mNG/SPARK-mNG-mAID or Δku80Δhxgprt/ TIR1 tachyzoites were concentrated by spinning and resuspended in a volume of 400 µl lysis buffer (5 mM NaCl, 142 mM KCl, 1 mM MgCl₂, 5.6 mM glucose, 25 mM HEPES pH 7.2 with 0.8% IGEPAL CA-630, 1X Halt Protease Inhibitors, and 250 U/ml benzonase). After 20 minutes on ice, the parasite lysates were combined with mNeonGreen-Trap Magnetic Agarose (ChromoTek) and were incubated with rotation at 4°C for 1 hour. The beads were washed three times with lysis buffer. Bound proteins were eluted with 5% SDS in 50 mM TEAB pH 8.5 for 10 minutes at 70°C. The eluates were reduced with 5 mM TCEP for 20 minutes at 50°C. Alkylation was performed with 20 mM MMTS at room temperature for 10 minutes. The samples were then precipitated, washed, and digested with 1 µg trypsin/LysC using the S-trap mini columns (Protifi) according to the manufacturer’s protocol. Eluted peptides were lyophilized and stored at -80°C until analysis. PKA C3-mNG samples were harvested and processed in the same way, except the corresponding control strain was mNG-TurboID-Ty/Δku80Δhxgprt.

SPARKEL-Ty IPs

Approximately 2 ✕10⁸ extracellular SPARKEL-mNG-Ty or Δku80Δhxgprt tachyzoites were concentrated by spinning and resuspended in a volume of 400 µl lysis buffer (10 mM Tris-HCl pH 7.5, 140 mM NaCl, 1% IGEPAL CA-630, 0.1% sodium deoxycholate, 0.1% SDS, Halt protease inhibitors, and Pierce universal nuclease). After 20 minutes on ice, the parasite lysates were combined with anti-Ty (BB2) conjugated Protein G beads and were incubated with rotation at 4°C for 1 hour (Huet et al., 2018). The beads were washed three times with lysis buffer. Bound proteins were eluted with 5% SDS in 50 mM TEAB pH 8.5 for 10 minutes at 70°C. The eluates were reduced with 5 mM TCEP for 20 minutes at 50°C. Alkylation was performed with 20 mM MMTS at room temperature for 10 minutes. The samples were then precipitated, washed, and digested with 1 µg trypsin/LysC using the S-trap mini columns (Protifi) according to the manufacturer’s protocol. Eluted peptides were lyophilized and stored at -80°C until analysis.

Mass spectrometry data acquisition and analysis

The lyophilized peptides were 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-nLC chromatography system using 0.1% formic acid in water as Buffer A and 0.1% formic acid in 80% acetonitrile as Buffer B. SPARK-mNG and PKA C3-mNG IP samples were separated at 300 nl/min on a 60-minute gradient of 1-2% B for 1 minute, 2-25% B for 41 minutes, 25-40% B for 6 minutes, and 40-100% B for 12 minutes. SPARKEL-Ty IP samples were separated at 300 nl/min on a 90-minute gradient of 1-6% B for 1 minute, 6-21% B for 41 minutes, 21- 36% B for 21 minutes, 36-50% B for 10 minutes, and 50-100% B for 14 minutes. Both methods ended with a seesaw gradient of 2% B for 6 minutes and 98% B for 6 minutes.

The orbitrap and FAIMS were operated in positive ion mode with a positive ion voltage of 1800V; with an ion transfer tube temperature of 270°C; using standard FAIMS resolution and compensation voltages of -50 and -65V. Peptides were filtered for a charge state of 2–6, and a dynamic exclusion of 20 seconds was applied. Full scan spectra were acquired in profile mode with 1 microscan at a resolution of 60,000, with a scan range of 350-1400 m/z, fill time of 25 milliseconds, normalized AGC target of 300%, intensity threshold of 5 × 10³. The top 15 MS2 spectra were acquired. MS2 spectra were generated with a HCD collision energy of 30 at a resolution of 15,000 with a first mass at 110 m/z, an isolation window of 1.3 m/z, AGC target of 200%, and auto injection time.

Raw mass spectrometry files were processed with the Proteome Discoverer 2.4 software (Thermo). Spectra were searched against the ToxoDB release49 GT1 protein database using Sequest HT (Thermo Fisher Scientific). The search included the following post-translational modifications: dynamic oxidation (+15.995 Da; M), dynamic acetylation (+42.011 Da; N-Terminus), dynamic phosphorylation (+79.966 Da; S, T, Y), and static methylthio (+45.988 Da; C) with a mass tolerance of 10 ppm for precursor ions and 0.2 Da for fragment ions. Enzymatic cleavage was set to trypsin, with 2 allowed missed cleavages per peptide. False discovery was assessed using Percolator with a concatenated target/decoy strategy using a strict FDR of 0.01, relaxed FDR of 0.05, and maximum Delta CN of 0.05. Protein abundance values were calculated using default workflows. Only unique peptide quantification values were used. Protein abundances were calculated from summation of peptide abundances. Sample ratios were calculated based on protein abundances. Hypothesis testing was performed by t- test on the background population of proteins and peptides. Protein abundances and statistical tests are reported in Supplementary File 1. The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE partner repository (Perez-Riverol et al., 2022) with the following dataset identifiers: SPARK-mNG-mAID IPs, PXD039922; SPARK-mNG IPs, PXD039919; SPARKEL-Ty IPs, PXD039979; PKA C3 IPs, PXD039896.

Proximity labeling of intracellular TurboID strains and mass spectrometry

Parasite treatment and harvest

Parasites freshly egressed from host cell monolayers and expressing SPARK-TurboID-Ty, SPARKEL- TurboID-Ty, or mNG-TurboID-Ty were infected onto a 15-cm dish confluent with HFFs (approximately 1-4 × 10⁷ parasites per dish). The parasites invaded and replicated in host cells for approximately 32 hours, at which point the media was aspirated and replaced with 10 ml of Ringer’s solution (155 mM NaCl, 2 mM CaCl₂, 3 mM KCl, 1mM MgCl₂, 3 mM NaH₂PO₄, 10 mM HEPES, 10 mM glucose) with 1% IFS and 500 µM biotin. The dishes were incubated at 37°C and ambient CO₂ for 30 minutes to allow equilibration of biotin across the host and parasite plasma membranes. The monolayers were then scraped, and parasites were mechanically released by passage through a 27-gauge syringe into a conical tube. The parasite suspension was diluted with phosphate-buffered saline (PBS) and was spun at 1,500 x g for 5 minutes at ambient temperature. The supernatant was decanted. The parasite pellet was resuspended in 1 ml PBS and transferred to a 1.5 ml protein low-bind tube. The tubes were spun at 4°C and 21,000 x g for 1 minute, after which the supernatant was discarded and the parasite pellet was resuspended in 1 ml PBS. This spin and wash cycle was repeated twice more, for three washes total to remove excess biotin and serum. After the final wash, the parasite pellet was resuspended in 500 µl lysis buffer (140 mM NaCl, 1% IGEPAL CA-630, 0.1% sodium deoxycholate, 0.1% SDS, 10 mM Tris-HCl, pH 7.5 supplemented with Pierce protease inhibitors and benzonase), was snap-frozen in liquid nitrogen, and was stored at –80°C until sample processing.

Biotinylated protein enrichment and MS

The frozen lysates were thawed and combined with 30 µl of Pierce magnetic streptavidin beads. The samples were rotated for 1 hour at room temperature. The beads were magnetically separated and were washed with 1 ml lysis buffer, 1 ml 1M KCl, 0.1 M NaCO3, 2M urea/10 mM Tris pH 8.0, and 1 ml lysis buffer. Elution of bound proteins was performed with two incubations of 10 minutes at 37 or 70°C in 5% SDS, 50 mM TEAB pH 8.5 with 2.5 mM biotin. The eluates were immediately reduced, alkylated, and digested as described in the immunoprecipitation methods. TurboID samples were analyzed by MS as were IP samples. The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE partner repository (Perez-Riverol et al., 2022) with the following dataset identifiers: SPARK-TurboID, PXD039983; SPARKEL-TurboID, PXD039985; TurboID- SPARKEL, PXD039986.

Immunofluorescence microscopy

Where indicated, parasites were infected onto confluent HFFs and were allowed to grow for 24 hours. Coverslips or coverglass were fixed in 4% formaldehyde in PBS unless otherwise noted. Following three washes in PBS, the fixed cells were permeabilized with 0.25% triton for 10 minutes at room temperature followed by three washes with PBS. The coverslips were incubated in blocking solution (2% BSA in PBS) for 10 minutes at room temperature, followed by a 60-minute incubation in primary antibody solution. After three washes with PBS, the coverslips were incubated in blocking solution at room temperature for 5 minutes, followed by a 60-minute 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 minutes at 37°C. Imaging was performed with the Nikon Ti Eclipse and NIS Elements software package.

Immunoblotting

SPARK/SPARKEL co-depletion immunoblotting

HFFs were infected with parasites of the TIR1, SPARKEL-HaloTag-mAID-Ty/TIR1, or SPARK-V5- mCherry-HA/SPARKEL-V5-HaloTag-mAID-Ty/TIR1 strains. Approximately 16 hours later, IAA was added to a final concentration of 500 µM to initiate the 24-hour depletion period. The following day, IAA was added to the 3- and 1-hour treatment samples. The host cell monolayers were scraped and harvested with spinning for 7 minutes at 1,500 x g. The parasite pellet was resuspended in 1 ml of HHE and was spun again for 10 minutes at 2,000 x g and 4°C. The pellet was resuspended to 2 × 10⁷ parasite equivalents per 10 µl in lysis buffer (5 mM NaCl, 142 mM KCl, 1 mM MgCl₂, 5.6 mM glucose, 25 mM HEPES pH 7.2 with 0.8% IGEPAL CA-630, 1X Halt Protease Inhibitors, 1X PMSF, and 250 U/ml benzonase). The lysates were diluted with 5x laemmli buffer (10% SDS, 50% glycerol, 300 mM Tris HCl pH 6.8, 0.05% bromophenol blue) and were boiled at 95°C for 5 minutes. Samples were run with 1 × 10⁷ parasite equivalents per lane on a precast 4-15% SDS-PAGE gel. The gel was incubated in 20% ethanol for 10 minutes and subsequently transfer buffer (25 mM TrisHCl, 192 mM glycine, 20% methanol) for 10 minutes. Samples were transferred onto a nitrocellulose membrane at 4°C for 1h at 100V. The membrane was blocked for 1h in 5% milk in TBS-T and incubated overnight in primary antibody solution (mouse anti-V5 at 1:1,000 dilution or mouse anti-TUB1 at 1:5,000 dilution) in 5% milk/TBS-T at 4°C. The membrane was incubated in secondary antibody solution (1:10,000 Goat anti-Mouse IgG IRDye 800 or 680) in 5% milk/TBS-T for 1h at room temperature and was visualized by LI-COR Odyssey CLx. Three 5-minute washes were performed between incubations in antibodies.

PKA C3-AID immunoblotting

HFFs were infected with parasites of the TIR1 and PKA C3-V5-mNG-mAID-Ty/TIR1 strains. Approximately 16 hours later, IAA was added to a final concentration of 500 µM to initiate the 24-hour depletion period. The following day, IAA was added to the 3- and 1-hour treatment samples. The host cell monolayers were scraped and harvested with spinning for 7 minutes at 1,500 x g. The parasite pellet was resuspended in 1 ml of HHE and was spun again for 10 minutes at 2,000 x g and 4°C. The pellet was resuspended to 1 × 10⁷ parasite equivalents per 10 µl in lysis buffer (5 mM NaCl, 142 mM KCl, 1 mM MgCl₂, 5.6 mM glucose, 25 mM HEPES pH 7.2 with 0.8% IGEPAL CA-630, 1X Halt Protease Inhibitors, 1X PMSF, and 250 U/ml benzonase). The lysates were diluted with 5x laemmli buffer (10% SDS, 50% glycerol, 300 mM Tris HCl pH 6.8, 0.05% bromophenol blue) and were boiled at 95°C for 5 minutes. Samples were run with 2 × 10⁷ parasite equivalents per lane on a precast 4-15% SDS-PAGE gel. The gel was incubated in 20% ethanol for 10 minutes and subsequently transfer buffer (25 mM TrisHCl, 192 mM glycine, 20% methanol) for 10 minutes. Samples were transferred onto a nitrocellulose membrane at 4°C for 1h at 100V. The membrane was blocked for 1h in 5% milk in TBS-T and incubated overnight in primary antibody solution (mouse anti-V5 at 1:1,000 dilution or mouse anti-TUB1 at 1:5,000 dilution) in 5% milk/TBS-T at 4°C. The membrane was incubated in secondary antibody solution (1:10,000 Goat anti-Mouse IgG IRDye 800 or 680) in 5% milk/PBS for 1h at room temperature and was visualized by LI- COR Odyssey CLx. Three 5-minute washes with TBS-T or PBS were performed between incubations in primary and secondary antibodies, respectively. The expected molecular weight of PKA C3-V5-mNG- mAID-Ty is 97.4 kDa, from the 40.2 kDa tagging payload (Smith et al., 2022) and 57.2 kDa protein sequence.

Plaque assays

500 parasites of the indicated strain were inoculated into 12-well 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 500 µM. Plates were washed with PBS and fixed for 10 min at room temperature with 70% 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.

Replication assays

Parasites recently egressed from host cells were allowed to invade coverslips confluent with HFFs for 1 hour, after which the media was aspirated and replaced with media containing 500 µM IAA or vehicle. After 24 hours, coverslips were fixed with 4% formaldehyde for 10 minutes at room temperature. The samples were prepared for immunofluorescence microscopy as described above, using anti-GAP45 as a parasite counterstain. The number of parasites per vacuole was enumerated from three fields of view acquired using a 40x objective. Significance was assessed using ANOVA and Tukey’s test in R (4.0.4).

Flow cytometry

Parasites were infected into 6-well dishes. Eight hours later, the media was replaced. IAA was added to wells at the appropriate time periods prior to sample harvest. The monolayers were scraped, and parasites were mechanically released from host cells by passage through a 27-gauge syringe. Parasites were then passed through a 5 µm filter. The samples were concentrated for 5 minutes via spinning at 1,000 x g. The media was replaced with 1 ml PBS. The parasite solutions were analyzed by flow cytometry with a Miltenyi MACSQuant VYB. Histograms are representative of two biological replicates.

Invasion assays

Parasites expressing SPARKEL-AID and SPARK-AID (Smith et al., 2022) were grown in HFFs for 40h with 24, 3, or 0h of treatment with 500 µM IAA. Extracellular parasites or parasites mechanically released from the host cells were spun at 1,500 x g for 5 minutes and were resuspended in invasion media (HEPES-buffered DMEM without phenol red) supplemented with 1% IFS. The parasite suspension was normalized to 1 ✕10⁶ cells/ml, and 2 ✕10⁵ parasites were added to three wells of a 96- well clear-bottom plate with confluent HFFs. The plate was spun at 290 x g for 5 minutes to distribute parasites over the host cells. Parasites were allowed to invade for 20 minutes at 37°C and 5% CO₂. The plate was washed once with PBS to remove unattached parasites and was fixed with 4% formaldehyde in PBS for 15 minutes. Extracellular parasites were stained with mouse anti-SAG1 for 30 minutes at room temperature. Following permeabilization with 0.25% triton-X100 for 10 minutes, all parasites were stained with guinea-pig anti-CDPK1 for 30 minutes at room temperature. The wells were incubated with a secondary antibody solution containing Hoechst for 30 minutes. The plate was imaged at 20x magnification using a Nikon Ti Eclipse epifluorescence scope. The number of parasites invaded was calculated by normalizing the number of intracellular, invaded parasites to host cell nuclei. Assays were performed in biological triplicate. Mean invasion efficiencies for each biological replicate were calculated from technical triplicates. Plotting and statistical tests were performed with ggplot2 and the compare_means() function with a one-sided t-test.

Zaprinast egress assays and GCaMP6f measurements

Parasites expressing SPARKEL-AID and pMIC2-MIC2-Gluc-P2A-GCaMP6f at a defined genomic locus (Herneisen et al., 2022) were grown in HFFs in glass-bottom dishes (35 mm, 1.5, Ibidi) for 25h with 24, 3, or 0h of treatment with 500 µM IAA. The dishes were washed once with one volume of Ringer’s buffer (155 mM NaCl, 2 mM CaCl₂, 3 mM KCl, 1mM MgCl₂, 3 mM NaH₂PO₄, 10 mM HEPES, 10 mM glucose), and the media was replaced with 500 µl of Ringer’s solution. Imaging was recorded every 2s for 300 using an Eclipse Ti microscope (Nikon) with an enclosure maintained at 37 °C. After approximately 10s of imaging, 1 ml of 750 µM zaprinast was added to the dish, for a final compound concentration of 500 µM. The same procedure was used to measure calcium mobilization in parasites expressing SPARK-AID and GCaMP6f (Smith et al., 2022). GCaMP6f intensity was measured in regions of interest around vacuoles containing two or more parasites using FIJI. Assays were performed on different days in biological triplicate. Vacuoles were manually scored for time to egress. Plotting and statistical tests were performed with ggplot2 and the scale_compare_means() function with a one-sided t-test.

Compound 1 treatment egress assays

HFF monolayers in a clear-bottomed 96-well plate were infected with 1 × 10⁵ parasites of the TIR1 or SPARK-AID strains in fluorobrite media supplemented with 10% calf serum. IAA or PBS were added to a final concentration of 500 µM 24 hours before the start of the assay. The media on the wells was aspirated, and wells were washed once with PBS to remove extracellular parasites. A fluorobrite solution containing a 3x concentration of compound 1 was added to each well, and the plate was incubated at 37°C and 5% CO₂ for 10 minutes. Ringer’s solution with 2 µM A23187 was added to each well, resulting in dilution of compound 1 to a 1x concentration. The plate was incubated at 37°C and 5% CO₂ for 10 minutes, after which 16% formaldehyde was added to each well, resulting in a final concentration of 4% formaldehyde in each well. Fixation occurred for 10 minutes. The plate was washed with PBS three times. Wells were permeabilized with 0.25% Triton X-100 followed by three washes with PBS. Parasites were stained with guinea pig anti-CDPK1 (Waldman et al., 2020) for 1 hour at room temperature. After three washes with PBS, the wells were incubated in a secondary antibody solution (Alexa Fluor 594 goat anti-guinea pig) for 30 minutes. The wells were washed three times with PBS and were imaged using a BioTek Cytation 3; four fields were acquired at 20x magnification per well. Percent egress was calculated as the ratio of liberated parasites to intact vacuoles. Curves were fit to the percent egress data in Prism using a sigmoidal dose-response model and difference between the TIR1 and SPARK-AID response curves was assessed with an extra sum of squares F-test. IC₅₀ values from individual replicates were assessed for difference using a paired t-test in Prism.

Differentiation assays

Plate-based differentiation assays

HFFs in black glass-bottom 96-well plates were infected with three dilutions ranging from 5 × 10³–2 × 10⁴ tachyzoites. The plates were spun for 5 minutes at 290 x g. After a four-hour invasion period at 37°C and 5% CO₂, the media was exchanged for DMEM supplemented with 10% FBS for spontaneous differentiation assays or RPMI with 1% FBS for alkaline stress differentiation assays. Spontaneous and alkaline stress differentiation assay samples were incubated at 37°C/5%CO₂ or 37°C/0% CO₂ for 48 hours, respectively, before fixation in 100% cold methanol for 2 minutes. The assays were incubated for 10 minutes at room temperature in blocking solution (1% bovine serum albumin in PBS). The wells were then incubated in a primary antibody blocking solution containing rabbit anti-GAP45 or biotinylated dolichos (Vector labs B-1035) for 30 minutes at room temperature. Secondary staining was performed with streptavidin-APC and goat anti-rabbit Alexa fluor 488 in blocking solution for 30 minutes at room temperature. Three PBS washes were performed between antibody incubations. Images were acquired with a widefield Nikon Ti epifluorescence scope. The percentage of DBL+ vacuoles was determined by manual quantification of vacuoles with DBL staining as a proportion of vacuoles with GAP45 staining. Images included in figure panels were false-colored in FIJI.

Coverglass-based differentiation assays

HFFs confluent on coverglass were infected with 1 × 10⁴–5 × 10⁴ tachyzoites. After a four-hour invasion period at 37°C and 5% CO₂, the media was exchanged for DMEM supplemented with 10% FBS for spontaneous differentiation assays or RPMI with 1% FBS for alkaline stress differentiation assays. Spontaneous and alkaline stress differentiation assay samples were incubated at 37°C/5%CO₂ or 37°C/0% CO₂ for 48 hours, respectively, before fixation in 4% formaldehyde in PBS for 10 minutes. The assays were incubated for 10 minutes at room temperature in blocking solution (1% bovine serum albumin in PBS). The wells were then incubated in a primary antibody blocking solution containing rabbit anti-GAP45 or guinea pig anti-CDPK1 or biotinylated dolichos (Vector labs B-1035) for 60 minutes at room temperature. Secondary staining was performed with streptavidin-APC and goat anti-rabbit Alexa fluor 488 or goat anti-guinea pig Alexa fluor 488 in blocking solution for 60 minutes at room temperature. Three PBS washes were performed between antibody incubations. Images were acquired with a widefield Nikon Ti epifluorescence scope with 20x objective, with 3 images per coverslip. The percentage of DBL+ vacuoles was determined by manual quantification of vacuoles with DBL staining as a proportion of vacuoles with GAP45 staining. Images included in figure panels were false-colored in FIJI.

SPARK 24h depletion proteome

Parasite treatment and harvest

HFF monolayers were infected with parasites expressing SPARK-mNG-mAID. Eight hours later, uninvaded parasites were washed away. The replacement media contained 500 µM IAA or PBS vehicle. After 24 hours, the host cell monolayers were scraped, and parasites were mechanically released by passage through syringes with 27-gauge needles. The suspension was passed through 5 µm filters and collected into conical tubes, which were spun for 7 minutes at 1,000 x g. The supernatant was decanted, and the parasite pellet was resuspended in 1 ml Fluorobrite DMEM lacking serum. The suspension was spun for 1,000 x g for 5 minutes. The supernatant was aspirated, and an additional wash in Fluorobrite DMEM followed by a spin at 1,000 x g for 5 minutes was performed. The parasite pellet was resupended in 250 µl, and four technical replicates of 50 µl were combined with 50 µl of 2x lysis buffer (10% SDS, 4 mM MgCl₂, 100 mM TEAB pH 7.55 with 2X Halt Protease and Phosphatase Inhibitors and 500 U/ml benzonase). The experiment was performed in biological duplicate.

Proteomics sample preparation and TMT labeling

Samples were prepared as previously described (Herneisen et al., 2022). Lysates were reduced 5 mM TCEP at 55°C for 15 minutes. Alkylation was performed with 20 mM MMTS for 10 minutes at room temperature. After the lysates were acidified to a final concentration of 2.5% v/v phosphoric acid, a 6X volume of S-trap binding buffer (90% methanol, 100 mM TEAB pH 7.55) was added. The solution was loaded onto S-trap mini columns (Protifi) and spun at 4,000 x g until all of the solution had been passed through the column. The columns were washed four times with 150 µl S-trap binding buffer, followed by a 30 second spin at 4,000 x g. Proteins were digested overnight in 20 µl of 50 mM TEAB pH 8.5 containing 2 µg of trypsin/LysC mix (Thermo Fisher Scientific) at 37°C in a humidified incubator. 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 were snap-frozen and lyophilized.

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 hour 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. The eluted peptides were snap-frozen and lyophilized until dry. Samples were fractionated with the Pierce High pH Reversed-Phase Peptide Fractionation Kit according to the manufacturer’s instructions for TMT-labeled peptides.

LC-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-nLC chromatography system. Peptides were separated at 300 nl/min on a gradient of 5–20% B for 110 minutes, 20–28% B for 10 minutes, 28–95% B for 10 minutes, 95% B for 10 minutes, and a seesaw gradient of 95–2% B for 2 minutes, 2% B for 2 minutes, 2–98% B for 2 minutes, 98% B for 2 minutes, 98–2% B for 2 minutes, and 2% B for 2 minutes. The orbitrap and FAIMS were operated in positive ion mode with a positive ion voltage of 1800V; with an ion transfer tube temperature of 270°C; using FAIMS user defined mode with compensation voltages of -50 and - 65V, an inner electrode temperature of 100°C, an outer electrode temperature of 85°C, with 4.6 ml/min carrier gas and a default charge state of 2. 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 × 10⁴, 2-5 charge state, and dynamic exclusion of 30 seconds with a cycle time of 1.5 seconds between master scans. MS2 spectra were generated with a HCD collision energy of 32 at a resolution of 45,000 with a first mass at 110 m/z, an isolation window of 0.7 m/z, 200% AGC target, and 120 ms injection time.

Data analysis

Raw files were analyzed by Proteome Discoverer 4.2. Peak lists were generated using the Sequest HT search engine and ToxoDB GT1 version 49 sequence database. Trypsin was specified as the digestion enzyme, with a maximum of two missed cleavages. Modifications included dynamic oxidation (+15.995 Da; M), dynamic acetylation (+42.011 Da; N-terminus), static TMT6plex (+229.163 Da; N-terminus), static TMTpro (+304.207 Da; K), and static methylthio (+45.988 Da; C). The enriched sample analysis additionally searched for dynamic phosphorylation (+79.966 Da; S, T, Y). The allowed mass tolerance for precursor and fragment ions was 10 ppm and 0.02 Da, respectively. False discovery was assessed using Percolator with a concatenated target/decoy strategy using a strict FDR of 0.01, relaxed FDR of 0.05, and maximum Delta CN of 0.05. Protein abundance values were calculated using default workflows. Only unique peptide quantification values were used. Co-isolation and signal-to-noise thresholds were set to 50% and 10, respectively. Protein abundances were calculated from summation of peptide abundances. Normalization was performed according to total peptide amount. Sample ratios were calculated based on protein abundances. Hypothesis testing was performed by ANOVA against individual proteins. Protein abundances are reported in Supplementary File 3. 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 PXD040598.

SPARK, SPARKEL, and PKA C3 depletion time course proteomes

Parasite treatment and harvest

HFF monolayers were infected with parasites expressing SPARK-mNG-mAID, SPARKEL-HaloTag-mAID, or PKA C3-mNG-mAID. Eight hours later, uninvaded parasites were washed away. IAA was added to 500 µM for the indicated periods of time before the samples were harvested together 32 hours post-infection. The host cell monolayers were scraped, and parasites were mechanically released by passage through syringes with 27-gauge needles. The suspension was passed through 5 µm filters and collected into conical tubes, which were spun for 7 minutes at 1,000 x g. The supernatant was decanted, and the parasite pellet was resuspended in 1 ml Fluorobrite DMEM lacking serum. The suspension was spun for 1,000 x g for 5 minutes. The supernatant was aspirated, and an additional wash in Fluorobrite DMEM followed by a spin at 1,000 x g for 5 minutes was performed. The parasite pellet was resupended in 50 µl of 1x lysis buffer (5% SDS, 2 mM MgCl₂, 100 mM TEAB pH 7.55 with 1X Halt Protease and Phosphatase Inhibitors and 500 U/ml benzonase). The experiment was performed in biological duplicate.

Proteomics sample preparation and TMT labeling

Samples were reduced, alkylated, and digested as described for the SPARK 24h depletion proteome. The lyophilized peptides were quantified using a Pierce fluorometric peptide assay. Samples destined for phosphopeptide enrichment were normalized to 100 µg of peptides per channel and were labeled with 200 µg of TMT10plex reagents according to the manufacturer’s protocol, with the following modifications. The eluates were combined with 200 µg of TMT10plex reagent in 15 µl of acetonitrile, for an estimated 1:2 w/w peptide:tag labeling reaction. The labeling proceeded for 1 hour at room temperature and was quenched for 15 minutes with 5% hydroxylamine. The samples were then pooled, flash-frozen, and lyophilized to dryness. The SPARKEL-AID proteome, which was not phospho-enriched, contained only 7 µg per channel and was labeled with 15 µg of TMT10plex reagent.

Phosphopeptide enrichment

Phosphopeptide enrichment was performed as previously described (Herneisen et al., 2022) using the SMOAC protocol (Tsai et al., 2014). Resuspended TMT10plex-labeled samples were enriched with the High-Select™ TiO2 Phosphopeptide Enrichment Kit (Thermo Fisher Scientific A32993). The flow-through and the eluate from IMAC enrichment were immediately snap-frozen and lyophilized. The flow-through was resuspended and enriched with the High-Select™ Fe-NTA Phosphopeptide Enrichment Kit (Thermo Fisher Scientific A32992) according to the manufacturer’s instructions. The eluted phosphopeptides were immediately snap-frozen and lyophilized. 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-labeled peptides. The phosphopeptides enriched with the TiO2 and Fe-NTA methods were combined prior to fractionation.

LC-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-nLC chromatography system as described above. Peptides were separated at 300 nl/min on a gradient of 6–21% B for 41 minutes, 21–36% B for 20 minutes, 36– 50% B for 10 minutes, and 50 to 100% B over 15 minutes. The orbitrap and FAIMS were operated in positive ion mode with a positive ion voltage of 1800V; with an ion transfer tube temperature of 270°C; using standard FAIMS resolution and compensation voltages of -50 and -65V (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 × 10³, 2-5 charge state, and dynamic exclusion of 30 seconds with a cycle time of 2 seconds 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.

Data 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 GT1 version 49 sequence database. The unenriched sample search included the following post-translational modifications: dynamic oxidation (+15.995 Da; M), dynamic acetylation (+42.011 Da; N-terminus), static TMT6plex (+229.163 Da; any N-terminus), static TMT6plex (+229.163 Da; K), and static methylthio (+45.988 Da; C). The enriched sample search included the same post-translational modifications, but with the addition of dynamic phosphorylation (+79.966 Da; S, T, Y). The mass spectrometry proteomics datasets have been deposited to the ProteomeXchange Consortium via the PRIDE partner repository (Perez-Riverol et al., 2022) with the following dataset identifiers: SPARK-AID phosphoproteomics, PXD040602; SPARKEL-AID proteomics, PXD040635; PKA C3-AID phosphoproteomics, PXD042989.

Exported peptide and protein abundance files from Proteome Discoverer 2.4 were loaded into R (version 4.0.4). The normalized abundance values of all peptides were used as input for a principal component analysis using the R stats package (version 3.6.2). The strain basal downregulation score was determined for the SPARK-AID phosphoproteome by calculating the modified z-score of the log₂- ratios of altered abundances in the SPARK-AID strain with no IAA treatment relative to the TIR1 parental strain. The IAA score was calculated by summing the log₂-ratios of altered abundances in the SPARK-AID strain relative to the SPARK-AID untreated samples, followed by transformation into a modified z-score. Clustering analysis was performed with the mclust package (Scrucca et al., 2016) (version 5.4.7) using the log₂ ratios of relative to the vehicle-treated SPARK-AID samples. Cluster assignments are reported in Supplementary File X. The PKA C3-AID phosphoproteome IAA score was calculated as described above; a strain basal downregulation score was not calculated because the experiment did not include a TIR1 parental sample.

Enrichment Analysis

Gene ontology enrichment analysis was performed as described (Herneisen et al., 2022). Sets of gene ontology terms from the differentially regulated (SPARK-dependent) and background proteome (all proteins with quantification values in the SPARK-AID mass spectrometry experiment) were downloaded from ToxoDB.org (Molecular Function, Computed evidence, P-Value cutoff set to 1). Gene ontology terms were tested for enrichment across all gene ontology terms identified in the background proteome. A p-value for the likelihood of a given enrichment to have occurred by chance was obtained using a hypergeometric test.

PKA R depletion phosphoproteome

Parasite treatment and harvest

HFF monolayers were infected with parasites expressing PKA R-V5-mAID-3HA. The next morning, IAA was added to 500 µM one hour prior to sample harvest. The host cell monolayers were scraped, and parasites were mechanically released by passage through syringes with 27-gauge needles. The suspension was passed through 5 µm filters and collected into conical tubes, which were spun for 7 minutes at 1,000 x g. The supernatant was decanted, and the parasite pellet was resuspended in 1 ml Fluorobrite DMEM lacking serum. The suspension was spun for 1,000 x g for 5 minutes. The supernatant was aspirated, and an additional wash in Fluorobrite DMEM followed by a spin at 1,000 x g for 5 minutes was performed. The parasite pellet was resuspended in 250 µl Fluorobrite DMEM containing 500 µM zaprinast, and 50 µl of this suspension was combined at 0, 10, 30, and 60 seconds post-treatment with 50 µl of 2x lysis buffer for a final composition of 5% SDS, 2 mM MgCl₂, 100 mM TEAB pH 7.55 with 1X Halt Protease and Phosphatase Inhibitors and 500 U/ml benzonase. The experiment was performed in biological duplicate.

Proteomics sample preparation and TMT labeling

Samples were reduced, alkylated, and digested as described previously (Herneisen et al., 2022). In brief, the eluates were reduced with 5 mM TCEP for 20 minutes at 50°C. Alkylation was performed with 20 mM MMTS at room temperature for 10 minutes. The samples were then precipitated, washed, and digested with 1 µg trypsin/LysC using the S-trap mini columns (Protifi) according to the manufacturer’s protocol. The lyophilized peptides were quantified using a Pierce fluorometric peptide assay. Samples destined for phosphopeptide enrichment were normalized to 50 µg of peptides per channel at a concentration of 1 µg/µl and were labeled with 250 µg of TMTpro reagents according to the manufacturer’s protocol, with volumes adjusted accordingly. The labeling proceeded for 1 hour at room temperature and was quenched for 15 minutes with 5% 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. The eluted peptides were snap-frozen and lyophilized until dry.

Phosphopeptide enrichment

Phosphopeptide enrichment was performed as previously described (Herneisen et al., 2022) 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). The flow-through and the eluate from IMAC enrichment were immediately snap-frozen and lyophilized. The flow-through was resuspended and enriched with the High-Select™ Fe-NTA Phosphopeptide Enrichment Kit (Thermo Fisher Scientific A32992) according to the manufacturer’s instructions. The eluted phosphopeptides were immediately snap-frozen and lyophilized. Unenriched proteome samples were fractionated with the Pierce High pH Reversed-Phase Peptide Fractionation Kit according to the manufacturer’s instructions for TMT-labeled peptides. The phosphopeptides enriched with the TiO2 and Fe-NTA methods were combined prior to fractionation. Enriched proteome samples were fractionated with the Pierce High pH Reversed-Phase Peptide Fractionation Kit using the following gradients (% acetonitrile/0.1%triethylamine): 1.5, 3, 5, 7.5, 10, 15, 50. Samples were immediately flash-frozen and lyophilized to dryness.

LC-MS data acquisition

The fractions were lyophilized and resuspended in 5-15 µ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-nLC chromatography system with a 15 cm column with 75 µm diameter (Thermo Fisher ES900). Peptides were separated at 300 nl/min on a gradient of 5–20% B for 110 minutes, 20–28% B for 10 minutes, 28–95% B for 10 minutes, 95% B for 10 minutes, and a seesaw gradient of 95–2% B for 2 minutes, 2% B for 2 minutes, 2–98% B for 2 minutes, 98% B for 2 minutes, 98–2% B for 2 minutes, and 2% B for 2 minutes. The orbitrap and FAIMS were operated in positive ion mode with a positive ion voltage of 1800V; with an ion transfer tube temperature of 270°C; using standard FAIMS resolution and compensation voltages of -50 and -65V, 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 × 10⁴, 2-5 charge state, and dynamic exclusion of 30 seconds with a cycle time of 1.5 seconds 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.

Data 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 as described for the SPARK 24h depletion proteome. The mass spectrometry proteomics datasets have been deposited to the ProteomeXchange Consortium via the PRIDE partner repository (Perez-Riverol et al., 2022) with the dataset identifier PXD044398.

Exported peptide and protein abundance files from Proteome Discoverer 2.4 were loaded into R (version 4.0.4). Ratios relative to the vehicle-treated timepoint were standardized by a modified Z- score. Peptides with the same site phosphorylations were matched across the PKA R and SPARK depletion phosphoproteomes. We defined basally PKA C1-dependent phosphorylation sites as phosphopeptides that were up-regulated by more than two modified Z-scores when PKA R was depleted for one hour. The distributions of all phosphopeptides identified in both the SPARK and PKA R phosphoproteomes were compared with the distributions of PKA C1-dependent sites using a Kolmogorov-Smirnov test.

PKG depletion phosphoproteome

Parasite treatment and harvest

HFF monolayers were infected with parasites expressing PKG-mAID-3HA or the parental strain expressing TIR1. The next morning, IAA was added to 500 µM four hours prior to sample harvest. The host cell monolayers were scraped, and parasites were mechanically released by passage through syringes with 27-gauge needles. The suspension was passed through 5 µm filters and collected into conical tubes, which were spun for 7 minutes at 1,000 x g. The supernatant was decanted, and the parasite pellet was resuspended in 1 ml Fluorobrite DMEM lacking serum. The suspension was spun for 1,000 x g for 5 minutes. The supernatant was aspirated, and an additional wash in Fluorobrite DMEM followed by a spin at 1,000 x g for 5 minutes was performed. The parasite pellet was resuspended in 250 µl Fluorobrite DMEM containing 500 µM zaprinast, and 50 µl of this suspension was combined at 0, 10, 30, and 60 seconds post-treatment with 50 µl of 2x lysis buffer for a final composition of 5% SDS, 2 mM MgCl₂, 100 mM TEAB pH 7.55 with 1X Halt Protease and Phosphatase Inhibitors and 500 U/ml benzonase. The experiment was performed in biological duplicate.

Proteomics sample preparation and TMT labeling

Samples were reduced, alkylated, and digested as described previously (Herneisen et al., 2022). In brief, the eluates were reduced with 5 mM TCEP for 20 minutes at 50°C. Alkylation was performed with 20 mM MMTS at room temperature for 10 minutes. The samples were then precipitated, washed, and digested with 1 µg trypsin/LysC using the S-trap mini columns (Protifi) according to the manufacturer’s protocol. The lyophilized peptides were quantified using a Pierce fluorometric peptide assay. Samples destined for phosphopeptide enrichment were normalized to 50 µg of peptides per channel at a concentration of 1 µg/µl and were labeled with 250 µg of TMTpro reagents according to the manufacturer’s protocol, with volumes adjusted accordingly. The labeling proceeded for 1 hour at room temperature and was quenched for 15 minutes with 5% hydroxylamine. The samples were then pooled, flash-frozen, and lyophilized to dryness. The dry pooled sample was resuspended in 1 ml of 5% formic acid and was desalted using the EasyPep Maxi Kit (Thermo Scientific) according to the manufacturer’s instructions. The eluted peptides were flash-frozen and lyophilized to dryness.

Phosphopeptide enrichment

Phosphopeptide enrichment was performed as previously described (Herneisen et al., 2022) 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). The flow-through and the eluate from IMAC enrichment were immediately snap-frozen and lyophilized. The flow-through was resuspended and enriched with the High-Select™ Fe-NTA Phosphopeptide Enrichment Kit (Thermo Fisher Scientific A32992) according to the manufacturer’s instructions. The eluted phosphopeptides were immediately snap-frozen and lyophilized. Unenriched proteome samples were fractionated with the Pierce High pH Reversed-Phase Peptide Fractionation Kit according to the manufacturer’s instructions for TMT-labeled peptides. The phosphopeptides enriched with the TiO2 and Fe-NTA methods were combined prior to fractionation. Enriched proteome samples were fractionated with the Pierce High pH Reversed-Phase Peptide Fractionation Kit using the following gradients (% acetonitrile/0.1%triethylamine): 1.5, 3, 5, 7.5, 10, 15, 50. Samples were immediately flash-frozen and lyophilized to dryness.

LC-MS data acquisition

The fractions were lyophilized and resuspended in 5-15 µ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-nLC chromatography system with a 25 cm column with 75 µm diameter (Thermo Fisher ES902). Peptides were separated at 300 nl/min on a gradient of 5–20% B for 110 minutes, 20–28% B for 10 minutes, 28–95% B for 10 minutes, 95% B for 10 minutes, and a seesaw gradient of 95–2% B for 2 minutes, 2% B for 2 minutes, 2–98% B for 2 minutes, 98% B for 2 minutes, 98–2% B for 2 minutes, and 2% B for 2 minutes. The orbitrap and FAIMS were operated in positive ion mode with a positive ion voltage of 1800V; with an ion transfer tube temperature of 270°C; using standard FAIMS resolution and compensation voltages of -50 and -65V, an inner electrode temperature of 100°C, and outer electrode temperature of 85°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 × 10⁴, 2-5 charge state, and dynamic exclusion of 30 seconds with a cycle time of 1.5 seconds 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.

Data 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 as described for the SPARK 24h depletion proteome. The mass spectrometry proteomics datasets have been deposited to the ProteomeXchange Consortium via the PRIDE partner repository (Perez-Riverol et al., 2022) with the dataset identifier PXD044361.

Exported peptide and protein abundance files from Proteome Discoverer 2.4 were loaded into R (version 4.0.4). Ratios relative to the vehicle-treated timepoint were standardized by a modified Z- score. Peptides with the same site phosphorylations were matched across the PKG and SPARK depletion phosphoproteomes. We defined basally PKG-dependent phosphorylation sites as phosphopeptides that were down-regulated by more than two modified Z-scores when PKG was depleted for four hours. The distributions of all phosphopeptides identified in both the SPARK and PKG phosphoproteomes were compared with the distributions of PKG-dependent sites using a Kolmogorov-Smirnov test.

PKA C3 depletion transcriptomics

Parasite Harvest

TIR1/RH and PKA C3-mAID parasites were inoculated onto 15 cm dishes of HFF monolayers. After ∼6-8 hours, auxin was added to a final concentration of 500 uM. 24 hours after addition of auxin, dishes were placed on ice and washed twice with PBS containing 100 μg/mL cycloheximide. Lysis buffer (20 mM HEPES pH 7.5, 100 mM KCl, 5 mM MgCl₂, 1% Triton X-100, 2 mM DTT, 100 ug/mL cycloheximide, 500 U/mL RNasIn Plus, and cOmplete, Mini Protease Inhibitor Cocktail) was added to the dish, the monolayer was scraped, and lysate was transferred to an eppendorf tube to incubate for 10 minutes on ice. Lysate was passed through a 27-gauge needle 6 times to facilitate lysis. Lysate was then centrifuged for 10 minutes at 1300xg, and the supernatant was flash-frozen and stored at -80℃. 3 biological replicates of each strain were collected, and each sample was split into 2 technical replicates for downstream library preparation.

RNA-sequencing and Analysis

Lysates were thawed on ice and a spike-in of D. melanogaster mRNA was added to each sample before Trizol extraction of RNA. RNA was polyA selected using the NEXTFLEX® Poly(A) Beads 2.0 Kit, and sequencing libraries were prepared according to the manufacturer’s protocol with the NEXTFLEX® Rapid Directional RNA-Seq Kit 2.0 and NEXTFLEX® RNA-Seq 2.0 Unique Dual Index barcodes for multiplexing. Libraries were sequenced on one lane of a NovaSeq S4 with 150×150 paired-end reads at the Whitehead Institute Genome Technology Core.

Reads adapters were trimmed with TrimGalore (version 0.6.7). A metagenome of T. gondii, H. sapiens, and D. melanogaster was prepared using STAR (version 2.7.1a) from genome FASTA files and gff annotation files downloaded from VEuPathDb Release 59. Trimmed reads were aligned to the metagenome using STAR, and mapped reads were quantified using featureCounts (version 1.6.2). Differential expression analysis of T. gondii genes with at least 10 reads across all samples was performed using DESeq2 (version 1.30.1) in R.

SPARKEL phylogenetic analysis

SPARKEL (TGGT1_291010) homologs were identified from representative apicomplexan and metazoan (outgroup) genomes. The Skp1 domain sequences were extracted based on annotation with Interpro domain IPR011333 or PFam domain PF03931. Domain alignment was performed in MEGA using ClustalW with a gap opening penalty of 3, gap extension penalty of 1.8, and remaining default parameters. A Neighbor-Joining Tree was generated and visualized with MEGA using the bootstrap method with 1000 replications as a test of phylogeny and a Jones-Taylor-Thornton model and remaining default parameters.

Acknowledgements

We thank Faye Harling for assistance with cell culture, Emily Shortt for the CDPK1-mNG-mAID strain, L. David Sibley for ALD1 antibody and RH TIR1 strain, Dominique Soldati-Favre and Drew Etheridge for the GAP45 antibodies, Marc-Jan Gubbels for the TUB1 antibody, Kevin Brown for the PKG-AID strain, and members of the Lourido laboratory for helpful discussions. Aditi Shukla developed the BFD1 knockout strategy. The Whitehead Institute Genome Technology Core assisted in nucleic acid sequencing, and the Whitehead Flow Cytometry Core assisted with sorting of fluorescent parasite strains. This work relied extensively on VEuPathDB.org and we thank all contributors to this resource. This research was supported by funds from R01AI144369 to S.L. and a National Science Foundation Graduate Research Fellowship to A.L.H. (174530) and T.A.S. (2018259980).

Author details

Alice L Herneisen

Whitehead Institute for Biomedical Research, Cambridge, MA, USA

Biology Department, Massachusetts Institute of Technology, Cambridge, MA, USA

Contribution: Conceptualization, Methodology, Validation, Formal analysis, Investigation, Writing - original draft, Visualization

Competing interests: No competing interests declared ORCID: 0000-0003-3368-0893

Michelle L Peters

Whitehead Institute for Biomedical Research, Cambridge, MA, USA

Biology Department, Massachusetts Institute of Technology, Cambridge, MA, USA

Contribution: Methodology, Formal analysis, Investigation

Competing interests: No competing interests declared ORCID: 0000-0002-2919-2523

Tyler A Smith

Whitehead Institute for Biomedical Research, Cambridge, MA, USA

Biology Department, Massachusetts Institute of Technology, Cambridge, MA, USA Contribution: Resources, Investigation

Competing interests: No competing interests declared ORCID: 0000-0002-3399-4043

Sebastian Lourido

Whitehead Institute for Biomedical Research, Cambridge, MA, USA

Biology Department, Massachusetts Institute of Technology, Cambridge, MA, USA

Contribution: Conceptualization, Resources, Supervision, Funding acquisition, Methodology, Writing - review and editing

For correspondence: lourido@wi.mit.edu Competing interests: Reviewing editor, eLife ORCID: 0000-0002-5237-1095

Competing interests

The authors declare no additional competing interests.

Supplementary files

Supplementary File 1. Protein quantification and statistical tests for immunoprecipitation-mass spectrometry experiments (SPARK, SPARKEL, and PKA C3-tagged strains), exported from the Proteome Discoverer 2.4 software.

Supplementary File 2. Protein quantification and statistical tests for TurboID mass spectrometry experiments (SPARKEL and SPARK TurboID-tagged strains), exported from the Proteome Discoverer 2.4 software.

Supplementary File 3. Protein and peptide quantification and statistical tests for SPARK, SPARKEL, and PKA C3 depletion proteomes, exported from the Proteome Discoverer 2.4 software.

Supplementary File 4. Peptide ratios quantified in the PKA R-AID, and PKG-AID, and SPARK-AID depletion phosphoproteomes. Ratios are reported relative to the vehicle-treated sample within each experiment. Modified Z-scores were calculated by standardizing each value with respect to the median and median absolute deviation within each ratio.

Supplementary File 5. Oligonucleotides and DNA sequences used in this study.

Source data

Figure 2—source data 1. This file contains source data that was used to generate the blot in Figure 2B. V5, LICOR.

Figure 2—source data2. This file contains source data that was used to generate the blot in Figure 2B. TUB1, LICOR.

Figure 2—source data 3. This file contains source data that was used to generate the blot in Figure 2C. V5, LICOR.

Figure 2—source data 4. This file contains source data that was used to generate the blot in Figure 2C. CDPK1, LICOR.

Figure 2—figure supplement 1 source data 1. This file contains source data that was used to generate the blot in Figure 2—figure supplement 1B. V5, LICOR.

Figure 2—figure supplement 1 source data 2. This file contains source data that was used to generate the blot in Figure 2—figure supplement 1B. ALD1, LICOR.

Figure 6—source data 1. This file contains source data that was used to generate the blot in Figure 6D. V5, LICOR.

Figure 6—source data 2. This file contains source data that was used to generate the blot in Figure 6D. ALD1, LICOR.

Figure 7—source data 1. This file contains source data that was used to generate the blot in Figure 7D. mNG, LICOR.

Figure 7—source data 2. This file contains source data that was used to generate the blot in Figure 7D. V5, LICOR.

Figure 7—source data 3. This file contains source data that was used to generate the blot in Figure 7D. ALD1, LICOR.

Materials availability statement

All mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE partner repository with the accession numbers PXD039896, PXD039919, PXD039922, PXD039979, PXD039983, PXD039985, PXD039986, PXD040598, PXD040602, PXD040635, PXD042989, PXD044361, PXD044398. Sequences of cloning vectors generated for this study are listed in Supplementary File 5. Strains and plasmids generated for this study can be obtained by emailing the corresponding author.

Figure legends

Figure 1—figure supplement 1.

Figure 2—figure supplement 1.

Figure 3—figure supplement 1.

Figure 4—figure supplement 1.

Figure 6—figure supplement 1.

Figure 7—figure supplement 1.