Comparative Proteomics Analysis for Elucidating the Interaction Between Host Cells and Toxoplasma gondii

Toxoplasma gondii, a representative model organism belonging to the phylum Apicomplexa, can infect almost all warm-blooded organisms, including humans. The invasion of host cells via host–parasite interaction is the key step for T. gondii to complete its life cycle. Herein we performed tandem mass tag analysis to investigate global proteomic changes in host cells (human foreskin fibroblasts, HFFs) [HFFs infected with T. gondii (HT) vs. HFFs (H)] and T. gondii [HT vs. T. gondii (T)] during intracellular infection. Overall, 3477 and 1434 proteins were quantified, of which 375 and 1099 proteins were differentially expressed (adjusted p-value < 0.05 and >1.5 or <0.67-fold change) in host cells and T. gondii, respectively. T. gondii invasion relies on the secretion of numerous secretory proteins, which originate from three secretory organelles: micronemes, rhoptries, and dense granules. In the HT vs. T group, few secretory proteins were upregulated, such as microneme proteins (MICs: MIC6, MIC10), rhoptry bulb proteins (ROPs: ROP5, ROP17), and dense granule proteins (GRAs: GRA4, GRA5, GRA12). In contrast, dozens of known secretory proteins were significantly downregulated in T. gondii-infected HFFs. In HFFs, gene ontology and Kyoto Encyclopedia of Genes and Genomes pathway analyses revealed a large number of differentially expressed proteins (DEPs) enriched in metabolic processes and immune-associated signaling pathways, such as NF-κB, cAMP, and Rap1 signaling pathways. Further, in case of T. gondii, DEPs were involved in ribosome biogenesis, citrate cycle, and galactose metabolism, indicating that cell biosynthesis and metabolism of T. gondii were altered after host cell invasion. These findings reveal novel modifications in the proteome of host cells as well as T. gondii, helping us better understand the mechanisms underlying host–parasite interaction.


INTRODUCTION
Toxoplasmosis is a life-threatening opportunistic infection caused by the protozoan parasite Toxoplasma gondii, an obligate intracellular organism. Around one-third of the global population is considered to have a latent infection of T. gondii (Montoya and Liesenfeld, 2004). T. gondii is a representative model organism of the phylum Apicomplexa that contains many important parasites of humans and animals, including Plasmodium spp., Eimeria spp., Neospora, and Cryptosporidium (Zhang et al., 2012;Janousǩovec et al., 2015). T. gondii can infect almost all warm-blooded organisms, including humans, and can exist in various tissues and organs of infected hosts. Moreover, toxoplasmosis can be fatal in immunocompromised patients, such as those with AIDS and those undergoing immunosuppressive chemotherapy (Kim and Weiss, 2008;Flegr, 2013;Torrey and Yolken, 2013). In addition, T. gondii can often cause adverse pregnancy outcomes, such as abortion and preterm labor, and even if pregnant women can successfully deliver, surviving offspring often show serious neurological and ocular sequelae (Robert-Gangneux et al., 2011;Jiao et al., 2017).
T. gondii actively invades host cells and can infect and replicate within almost any nucleated cell, which suggests that the pertinent mechanism is precise as well as effective. Host cell invasion is pivotal for T. gondii to complete its life cycle and is a consequence of host-parasite interaction (Hakimi et al., 2017). T. gondii participates in the signal transduction of host cells, making them prone to invasion. The invasion of host cells by T. gondii can be achieved in merely 15-20 s and relies on the actin-myosin motor complex and multiprotein complexes secreted by specialized organelles at its apical end, such as micronemes, rhoptries, and dense granules (Rabenau et al., 2001;Freńal et al., 2010;Swierzy et al., 2017). These secretory organelles have a crucial function in host cell invasion and parasitophorous vacuole establishment . The microneme protein (MIC) plays an important role in the recognition, adhesion, and attachment of T. gondii to host cells during the early stage of invasion. Rhoptry proteins, comprising rhoptry neck proteins (RONs) and rhoptry bulb proteins (ROPs), are secreted during invasion and are believed to contribute to parasitophorous vacuole formation. Further, dense granule proteins (GRAs) participate in the modification of parasitophorous vacuole and its membrane, facilitating survival and replication of the organism.
The host-parasite response to an infection and the identification of host proteins that interact with T. gondii proteins should help in elucidating the mechanism of infection and virulence at the tachyzoite stage. Previous studies have majorly used two-dimensional electrophoresisto detect proteomic changes in human foreskin fibroblasts (HFFs) or T. gondii (Cohen et al., 2002;Nelson et al., 2008); however, this method has a relatively low throughput for protein identification or quantification. Herein we performed tandem mass tag (TMT) analysis to investigate global proteomic changes in both host cells and T. gondii during intracellular infection. We report novel insights into the proteomics landscape of T. gondii before and after infection of host cells and also report data pertaining to modifications in the host cell proteome during T. gondii infection, with the aim of elucidating the complex relationship between host cells and parasites.
For cell invasion experiments, tachyzoites were added to HFF monolayers (approx. 85% confluency) at an infection ratio of 3:1. Tachyzoites that had not invaded the cells were washed away with phosphate-buffered saline after 2 h. The cells were harvested at 24 h when approximately 80% of them appeared infected. The control groups included an equal number of non-infected cells and an equal number of T. gondii tachyzoites. All cells were washed using cold phosphate-buffered saline, centrifuged at 400 g for 10 min, and stored at −80°C until needed. Each sample was replicated three times to yield nine samples (three groups × three biological replicates) and the three groups were designated H (HFFs), T (T. gondii), and HT (HFFs infected with T. gondii).

Protein Extraction and Quantification
The samples were suspended in a lysis buffer (8 M urea and 1% protease inhibitor) and disrupted three times on ice by sonication. The supernatant was then separated by centrifugation at 4°C and 12,000 g for 10 min, and the pellet was removed. Finally, protein concentration was determined using a BCA protein assay kit (catalog no. P0011-1, Beyotime, China), according to manufacturer instructions.

Trypsin Digestion and TMT Labeling
Dithiothreitol (5 mM, catalog no. D9163, Sigma-Aldrich, USA) was added to the protein solution, and reduction was allowed to proceed at 56°C for 30 min. The protein solution was then incubated with 11 mM iodoacetamide (catalog no. V900335, Sigma-Aldrich, USA) at room temperature in darkness for 15 min. The urea concentration of the solution was then diluted to <2 M using 100 mM tetraethylammonium bromide (catalog no. 17902, Sigma-Aldrich, USA). Finally, trypsin (catalog no. V5111, Promega, USA) was added at a trypsin-toprotein mass ratio of 1:50 for the first round of digestion for overnight and 1:100 for the second round of digestion for 4 h. Peptides were desalted on a Strata X C18 SPE column (Phenomenex, USA) and vacuum dried after trypsin digestion. They were then labeled using a TMT 10-plex isobaric label kit (catalog no. 90111, Thermo, USA), according to manufacturer instructions.

High-Performance Liquid Chromatography Fractionation and LC-MS/MS
The labeled peptides were fractionated on a 300 Extend C18 column (4.6 mm × 250 mm, 5 mm, Agilent, USA) and separated with a gradient of 8%-32% acetonitrile (pH 9.0) over 60 min into 60 fractions. The peptides samples were then combined into 14 fractions and dried in a vacuum centrifuge. Subsequently, they were dissolved in 0.1% formic acid solution and loaded onto an EASY-nLC 1000 UPLC system (Thermo, USA). Solvent A was 0.1% formic acid in 90% acetonitrile solution. The gradient comprised an increase from 9% to 26% of solvent B (0.1% formic acid in 90% acetonitrile) over 40 min, a rise from 26% to 35% in 14 min, and then from 35% to 80% in 3 min, followed by holding at 80% for the last 4 min; the flow rate was kept constant at 350 nL/min at all times. The peptides were subjected to NSI source, followed by tandem mass spectrometry (MS/MS) on a Q Exactive Plus Orbitrap LC-MS/MS System (Thermo, USA). The electrospray voltage was 2.0 kV. The mass spectrometer scan range was 350-1800 m/z and the orbitrap resolution was set to 70,000. The MS/MS scan range and orbitrap resolution were 100 m/z and 35,000, respectively. Datadependent acquisition was used for data collection. Automatic gain control was set to 5E4, and the dynamic exclusion duration was 30 s.

Database Search
The raw files from the same batch were processed together with MaxQuant against the SwissProt Homo sapiens and UniPort T. gondii ATCC_50861 databases, concatenated with the reverse decoy database. The cleavage enzyme was trypsin/P, with no more than two missing cleavages and the minimum peptide length was seven. For precursor ions, the mass tolerance of the first search and main search was set to 20 ppm and 5 ppm, respectively. The mass tolerance of fragment ions was 0.02 Da. Carbamidomethylation of cysteine residues and methionine oxidation were specified as a fixed and variable modification, respectively. The protein quantification method was set with TMT 10-plex and the reporter ion isotopic impurity distribution of TMT 10-plex was set as instructed by to the product manual. The results were filtered based on the peptide false discovery rate (FDR) of < 1%.

Data Analysis and Bioinformatics
We chose unique peptides for protein quantification; relative quantification of protein expression levels were calculated as the median of unique peptides. There biological replicates were assessed for each group. For each protein, the average ratio calculated using the three biological replicates indicated final protein expression. The protein quantification results were statistically analyzed using two-sample t-test and the p-value was corrected by the method of Benjamin and Hochberg FDR analysis. Proteins with a fold change > 1.5 or < 0.67 and adjusted p-value < 0.05 were considered to be differentially abundant. Differentially expressed proteins (DEPs) were selected for subcellular localization, clusters of orthologous groups of protein, gene ontology (GO) classification, and enrichment analyses. WoLF PSORT (v0.2 http://www.genscript.com/psort/ wolf_psort.html) was used to predict the subcellular localization of DEPs. For GO analysis, all annotations were derived from the UniProt-GOA database (http://www.ebi.ac.uk/GOA/) and complemented with InterProScan soft (v5.14-53.0 http://www. ebi.ac.uk/InterProScan/). All DEPs were subjected to analysis of eukaryotic orthologous group (KOG) by aligning their sequences with the KOG protein sequence database (http://genome.jgi.doe. gov/help/kogbrowser.jsf). The Kyoto Encyclopedia of Genes and Genomes (KEGG) database was used for pathway enrichment analyses of DEPs. For GO and KEGG enrichment analyses, twotailed Fisher's exact test was applied to test DEPs against all identified proteins. A corrected p value of <0.05 indicated statistical significance. MS data have been deposited in the ProteomeXchange Consortium via the PRIDE (Perez-Riverol et al., 2019) partner repository, with the dataset identifier PXD021736.

Parallel Reaction Monitoring (PRM) Validation
To confirm the authenticity and accuracy of results of TMT analysis, we performed PRM assay that included 14 peptides. Protein extraction and digestion were achieved using the aforementioned methods. The peptides were subjected to MS/ MS on a Q Exactive Plus Orbitrap LC-MS/MS System (Thermo, USA). The electrospray voltage was 2.1 kV, MS scan range was 350-1050 m/z, and orbitrap resolution was 70,000. The MS/MS orbitrap resolution was 17,500. Data-independent acquisition was used for data collection. Automatic gain control of MS and MS/MS was set to 3E6 and 1E5, respectively. The maximum injection time was 50 ms for MS and 180 ms for MS/MS. The isolation window for MS/MS was 1.6 m/z. Three biological replicates were assessed.

Data Analysis of PRM Validation
The resulting MS data were processed using Skyline (v3.6). The peptide settings were as follows: the enzyme was set as trypsin [KR/P] and the max missed cleavage was set as 0. The peptide length was 7-25, and carbamidomethylation of cysteine residues was the fixed modification. The transition settings were as follows: precursor charges were set as 2, 3; ion charges were set as 1; and ion types were set as b, y. The product ions were set as from ion 3 to the last ion, and the ion match tolerance was 0.02 Da.

Protein Identification and Expression
From all nine samples, 240,537 spectra were generated, from which 26,251 peptides and 25,528 unique peptides were obtained. A total of 4111 and 1708 proteins were identified as human and T. gondii proteins, of which 3477 and 1434 proteins were quantified, respectively (Supplementary Table S1). The standard criteria were adjusted p-value of <0.05 and > 1.5 or < 0.67 -fold change accordingly, 150 up-and 225 downregulated proteins were identified in the HT group in comparison with the H group ( Figure 1A and Supplementary Table S2A). Further, 107 up-and 992 downregulated proteins were found in the HT group in comparison with the T group ( Figure 1B and Supplementary Table S2B). All DEPs were subsequently subjected to comprehensive bioinformatics analysis.

Comparative Proteomics Analyses of Host Cells in the HT vs. H Group
In the HT vs. H group, the subcellular location of DEPs was as follows: the nucleus (25.33%), cytoplasm (24.26%), and extracellular (21.06%; Figure 2A). As evident from Figure   We performed GO analysis to study biological changes in infected HFFs in the cellular component, biological process, and molecular function categories. Among the up-and downregulated proteins, 8, 14, and 8 terms in the cellular component, biological process, and molecular function categories were significantly enriched ( Figure 4 and Supplementary Table S3), respectively. In the biological process category, the enriched GO terms for immunity/defense-related functions were upregulated, such as regulation of lymphocyte-mediated immunity, regulation of leukocyte-mediated immunity, and regulation of adaptive immune response ( Figure 4A), whereas those for cell metabolism-related functions were downregulated, such as oligosaccharide catabolic process, glycosaminoglycan catabolic process, and aminoglycan catabolic process ( Figure 4B). According to KEGG pathway analysis, up-and downregulated proteins were significantly enriched in 19 and 25 pathways, respectively; some immune-related signaling pathways were present among them, such as the NF-kB signaling pathway, cAMP signaling pathway, cell adhesion molecules (CAMs), Rap1 signaling pathway, and calcium signaling pathway ( Figure 5A and Supplementary Table S4).
Host Proteins Associated With the Immune Response Were Up-and Downregulated Upon Infection With T. gondii T. gondii stimulated the host immune system. The expression of proteins associated with the immune response was upregulated, such as apoptotic proteins (Cytochrome c and TNF receptorassociated factor 1), adhesion proteins [intercellular adhesion molecule 1 (ICAM1), programmed cell death 1 ligand 1 (CD274), receptor-type tyrosine-protein phosphatase C (PTPRC), and integrin beta-2 (ITGB2)], and transcriptional regulatory factors (nuclear factor NF-kappa-B p100 subunit and prostaglandin G/H synthase 2). In contrast, the expression  of ras-related proteins (rap-1b, m-ras, rab-32, and r-ras2) was downregulated (Supplementary Table S2A).

Signaling Pathways in Host Cells Were Altered During T. gondii Infection
In the HT vs. H group, all DEPs could be mapped to 91 pathways, of which 25 were significantly enriched ( Figure 6 and Supplementary Table S4). These significantly enriched pathways could be classified into cell metabolism and catabolism pathways (glycan degradation, phagosome and lysosome), immune-related signaling pathways (leukocyte transendothelial migration, calcium signaling pathway, Ras signaling pathway, cell adhesion molecules (CAMs), Rap1 signaling pathway, p53 signaling pathway, cAMP signaling pathway, complement and coagulation cascades) and signaling molecules and interaction pathways (cytokine-cytokine receptor interaction and ECM-receptor interaction).

Comparative Proteomics Analysis of T. gondii in the HT vs. T Group
In the HT vs. T group, the subcellular location of DEPs was as follows: the nucleus (26.56%), cytoplasm (22.02%), and mitochondria (16.47%; Figure 2B).  Figure 3B).
In the HT vs. T group, 1099 DEPs were classified into eight terms of biological process, five terms of cellular component, and five terms of molecular function. Among the upregulated proteins, eight, six, and eight terms for cell component, biological process, and molecular function were significantly enriched, respectively, whereas the downregulated proteins were significantly enriched in one term of biological process and eight terms of molecular function (Figure 7 and Supplementary Table S5). In the biological process category, cell biosynthesis and metabolism were upregulated on infection with T. gondii, including cellular protein metabolic, peptide metabolic, and peptide biosynthetic processes ( Figure 7A). Further, KEGG pathway analysis of the HT vs. T group revealed that two significantly enriched pathways were upregulated, including ribosome and ribosome biogenesis in eukaryotes, and six enriched pathways were downregulated, including the pentose phosphate pathway, citrate cycle (TCA cycle), and galactose metabolism ( Figure 5B and Supplementary  Table S6).

T. gondii Infection Altered the Expression of Numerous Secretory Proteins
T. gondii invasion relies on the secretion of numerous secretory proteins, which mainly originate from three secretory organelles: micronemes, rhoptries, and dense granules. In the HT vs. T group, the expression of few secretory proteins was upregulated, such as MICs (MIC6, MIC10), ROPs (ROP5, ROP17), and GRAs (GRA4, GRA5, and GRA12). In contrast, the expression of dozens of known secretory proteins was significantly downregulated in T. gondii-infected HFFs ( Table  1). As shown in Table 1, 10 SAG1-related sequence proteins (SRSs) exhibited a 0.67-fold or lower decrease in their expression levels in the HT vs. T group, whereas four SRSs exhibited a 1.5-fold or greater increase in their expression levels, including SRS29A (2.938 fold), SRS36D (1.919 fold), SRS57 (1.769 fold), and SRS20A (1.524 fold).

Validation of DEPs
We performed PRM assay for an independent comparison and validation of the results of TMT analysis. In the HT vs. H group, six proteins-desmoplakin protein (DSP), plexin-B2 (PLXNB2), annexin A1 (ANXA1), early endosome antigen 1 (EEA1), DNA topoisomerase 1 (TOP1), and probable ATP-dependent RNA helicase (DDX5)-were selected for analysis. As evident from Table 2, the results of PRM assay were consistent with those of TMT analysis. In case of the HT vs. T group, we selected MIC1, MIC2, ROP5, GRA12, and two heat shock proteins (HSP28 and HSP90) for PRM assay. The expression level of MIC1, MIC2, and HSP90 was downregulated and that of ROP5, GRA12, and HSP28 was upregulated, validating the consistency of results obtained using the two methods ( Table 2 and Supplementary  Table S7).

DISCUSSION
With technological advancements, mass spectrometry-based proteomic approaches have rapidly developed; compared to two-dimensional electrophoresis and difference gel electrophoresis, they are high throughput methods for protein identification and quantification. Nelson et al. (2008) identified 157 DEPs in HFFs before and after infection with T. gondii using two-dimensional electrophoresis and mass spectrometry, whereas we herein identified 1400 DEPs (adjusted p-value < 0.05 and >1.2 or <0.83 fold change, Supplementary Table S1A) in response to infection with T. gondii using TMT analysis. The host demonstrated global reprogramming of cell metabolism upon T. gondii invasion, as evident from the differential expression levels of proteins involved in key metabolic pathways (glycolysis, lipid and sterol metabolism, mitosis, apoptosis, and structural protein expression). In the present study, we focused on host cell signaling pathways related to cell survival, innate recognition of parasites, and immune response.
In the process of infecting hosts, parasites secrete various proteins to regulate the host signal transduction mechanism and modulate the immune response of host cells (Shapira et al., 2005;Muniz-Feliciano et al., 2013). In this study, immune-related signaling pathways were significantly enriched, such as the NF-kB, cAMP, Rap1, Ras, cytokine-cytokine receptor interaction, p53, platelet activation, and MAPK signaling pathways (Supplementary Table S4), suggesting that the innate and adaptive immunity of the host was majorly modulated in response to T. gondii infection.
As obligate intracellular parasites, T. gondii must traverse the host cell plasma membrane to initiate an infection and to ensure survival and replication. T. gondii invasion and infection establishment relies on secretory proteins, including MICs, ROPs, RONs, GRAs, and SRSs (Manger et al., 1998;Carey et al., 2000;Lekutis et al., 2001;Tomley and Soldati, 2001;Mercier et al., 2002;Bradley et al., 2005). The majority of MICs comprise multiple copies of a conserved adhesive domain, allowing them to bind to host cell surface receptors or carbohydrates during pathogen invasion. Carruthers and Tomley (2008) indicated that MIC secretion is rapidly upregulated when parasites make contact with host cells. We herein observed that the expression levels of MIC6 and MIC10 were upregulated, but those of most other MICs (MIC1-5, MIC7, MIC8, MIC11, MIC15, MIC17A, and M2AP) were downregulated; this could be because we analyzed the proteome at a single timepoint. In addition, the adhesionrelated domains of MICs not only adhere to host cells but also promote complex formation between different MICs, making them more stable on the T. gondii surface and increasing its invasive capability. Multiple MICs have been previously reported, including MIC1-MIC4-MIC6, MIC2-M2AP, and MIC3-MIC8; all complex components were downregulated, except MIC6, in the current study. GRA6 can evidently activate the host transcription factor NFAT4, which promotes the synthesis of the chemokines Cxcl2 and Ccl2, and chemokines are known to control the spread of parasites (Ma et al., 2014;Hakimi et al., 2017). In this study, GRA6 expression level was downregulated after host cell invasion, consequently facilitating T. gondii survival. It has been reported that many SRSs are stage-specific; they play a key role in mediating the attachment to host cells and activate host immunity to regulate T. gondii virulence (Jung et al., 2004). For example, SRS57 serves as an adhesion protein by binding to the host cell surface, SRS29B is a virulence factor, and SRS34A and SRS16B play a crucial role in host cell immune modulation (Dzierszinski et al., 2000;Van et al., 2007;Tomita et al., 2018).
In the HT group, the expression level of many SRS proteins, including, SRS29A, SRS36D, SRS57, and SRS20A, was upregulated and that of others (SRS38B, SRS30A, SRS29C, SRS67, SRS34A, SRS43, SRS30C, SRS44, SRS25, and SRS52A) was downregulated in comparison with the expression levels in the T group. So far, few studies exist on SRSs, and the understanding of their various biological functions needs to be enhanced in future.
Parasite invasion is an active process that relies on the cytoskeleton for motility and also on the host cytoskeleton (Dobrowolski and Sibley, 1996). Host microtubules have been observed to play a role in interactions with invading T. gondii tachyzoites, and the host cytoskeleton apparently contributes to the anchoring of the moving junction (Sweeney et al., 2010;Takemae et al., 2013). Microtubules are composed of heterodimers of alpha-and beta-tubulin. Takemae et al. (2013) reported that TUBB2C and TUBB5 interact with TgRON4. In our study, TUBB2C and TUBB5 were not identified as DEPs in the HT vs. H group, but other alpha-and beta-tubulin proteins were detected, including TUBA1C (0.625 fold), TUBA4A (0.595 fold), TUBB2A (0.801 fold), TUBB3 (0.688 fold), TUBB4B (0.646 fold), and TUBB6 (0.753 fold) (Supplementary Table S1A), indicating that infection with T. gondii causes the remodeling of the host cytoskeleton. CAMs are distributed throughout regions of the plasma membrane and come in contact with other cells and extracellular matrix (Figliuolo da Paz et al., 2019). Herein the pathways of CAMs were significantly enriched and included neuronal growth regulator 1, ICAM1, syndecan-2, CD274, neural cell adhesion molecule L1, PTPRC, junctional adhesion molecule C, and ITGB2. ICAM1 is expressed on various cell types and participates in diverse cellular processes, including host-pathogen interactions (Bhalla et al., 2015). ICAM1 has been reported to interact with PfEMP1 in Plasmodium falciparum infection (Lennartz et al., 2019). Further, ICAM1 expression level has been noted to be upregulated upon host cell invasion by T. gondii in vitro and also during human toxoplasma uveitis in vivo, with ICAM1 interacting with TgMIC2 (Klok et al., 1999;Nagineni et al., 2000;Barragan et al., 2005). In this study, compared with uninfected HFFs, ICAM1 expression in infected cells was upregulated (1.96 fold). Blader et al. (2001) demonstrated that the ICAM1 gene was upregulated by at least 2-fold in HFFs after T. gondii invasion, and Nagineni et al. (2000) reported that ICAM1 secretion in human retinal pigment epithelial cells infected with T. gondii was augmented 5-fold over the uninfected control, and no significant changes in ICAM1 levels were found in human retinal pigment epithelial cells incubated with heat-killed T. gondii. These findings suggest that ICAM1 expression is derived from the interaction between host cells and the parasite.
Based on their roles in host-pathogen interaction, Spear et al. (2006) clustered host genes into three distinct classes: "pro-host", host genes required for host defense; "pro-parasite", host genes required for parasite growth; and "bystander", host genes incidentally regulated as a consequence of modulating the first two classes. In present study, we performed TMT-based proteomic analysis of T. gondii as well as host cells at 24 h of infection. The infection induced significant changes in the proteome of T. gondii and also led to the significant enrichment of immune-related signaling pathways in HFFs. Various DEPs, including cytoskeleton and immune-related proteins, were induced in response to host cell invasion, indicating that infection with T. gondii leads to considerable remodeling of the host cytoskeleton and also the crucial role played by the host immune system upon encountering T. gondii. Future studies are warranted to validate these proteins and to assess their clinical value.

DATA AVAILABILITY STATEMENT
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.

AUTHOR CONTRIBUTIONS
HS and JC designed this study and made the final revision. HS, JL and LW performed the experiments. KY, CX, GL and TX helped to create the tables and figures. BH and QW contributed to the reagents and materials. MG supervised the process. HS wrote the draft manuscript and JC revised the manuscript. All authors contributed to the article and approved the submitted version. to JC). The funders had no role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript.