Next Article in Journal
Medical Management versus PACK-CXL in Dogs with Infectious Keratitis: A Randomized Controlled Trial Protocol
Next Article in Special Issue
Molecular Detection of Babesia gibsoni in Cats in China
Previous Article in Journal
Effects of a Combined Geothermal and Solar Heating System as a Renewable Energy Source in a Pig House and Estimation of Energy Consumption Using Artificial Intelligence-Based Prediction Model
Previous Article in Special Issue
Comparison of a Commercial Enzyme-Linked Immunosorbent Assay (ELISA) with the Modified Agglutination Test (MAT) for the Detection of Antibodies against Toxoplasma gondii in a Cohort of Hunting Dogs
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Animals
Volume 12
Issue 20
10.3390/ani12202858
animals-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Glycosylation Analysis of Feline Small Intestine Following Toxoplasma gondii Infection

by
Bintao Zhai
1,2,
Shichen Xie
3,4,
Junjie Peng
2,
Yanhua Qiu
1,5,
Yang Liu
6,
Xingquan Zhu
3,4,7,
Junjun He
7,* and
Jiyu Zhang
1,*
1
Key Laboratory of Veterinary Pharmaceutical Development, Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Ministry of Agriculture, Lanzhou 730050, China
2
State Key Laboratory of Veterinary Etiological Biology, Key Laboratory of Veterinary Parasitology of Gansu Province, Lanzhou Veterinary Research Institute, Chinese Academy of Agricultural Sciences, Lanzhou 730046, China
3
College of Veterinary Medicine, Shanxi Agricultural University, Jinzhong 030801, China
4
Research Center for Parasites & Vectors, College of Veterinary Medicine, Hunan Agricultural University, Changsha 410128, China
5
College of Veterinary Medicine, Northwest A&F University, Xianyang 712100, China
6
College of Life Science, Ningxia University, Yinchuan 750021, China
7
Key Laboratory of Veterinary Public Health of Yunnan Province, College of Veterinary Medicine, Yunnan Agricultural University, Kunming 650201, China
*
Authors to whom correspondence should be addressed.
Animals 2022, 12(20), 2858; https://doi.org/10.3390/ani12202858
Submission received: 3 September 2022 / Revised: 12 October 2022 / Accepted: 17 October 2022 / Published: 20 October 2022
(This article belongs to the Special Issue Parasitic Diseases in Canines and Felines and Their Vectors)
Browse Figures
Review Reports Versions Notes

Abstract

:

Simple Summary

Toxoplasma gondii has a serious impact on public health and the economic development of animal husbandry. Glycosylation, especially N-glycosylation, the pattern modification of proteins, is closely related to the biological functions of proteins, and our study used it to analyze glycosylation alterations in the small intestine of cats infected with T. gondii. The results of the present study showed that 56 glycosylated peptides were upregulated and 37 glycosylated peptides were downregulated. Additionally, we also identified eight N-glycosylated proteins of T. gondii including eight N-glycopeptides and eight N-glycosylation sites. Moreover, the protein eEF2 and its corresponding peptide sequence were identified, with GO terms (i.e., cellular process and metabolic process, cell and cell part, and catalytic activity) that were significantly enriched in the T. gondii MAPK pathway. In addition, the Clusters of Orthologous Groups of proteins (COG) function prediction results showed that posttranslational modification, protein turnover, and chaperones (11%) had the highest enrichment for T. gondii. The host proteins ICAM-1 and PPT1 and the endoplasmic reticulum stress pathway may play an important role in the glycosylation of T. gondii-infected hosts. Our study may provide a new target for T. gondii detection to prevent the spread of T. gondii oocysts in the future.

Abstract

Toxoplasma gondii (T. gondii) is responsible for severe human and livestock diseases, huge economic losses, and adversely affects the health of the public and the development of animal husbandry. Glycosylation is a common posttranslational modification of proteins in eukaryotes, and N-glycosylation is closely related to the biological functions of proteins. However, glycosylation alterations in the feline small intestine following T. gondii infection have not been reported. In this study, the experimental group was intragastrically challenged with 600 brain cysts of the Prugniuad (Pru) strain that were collected from infected mice. The cats’ intestinal epithelial tissues were harvested at 10 days post-infection and then sent for protein glycosylation analysis. High-performance liquid chromatography coupled to tandem mass spectrometry was used to analyze the glycosylation alterations in the small intestine of cats infected with T. gondii. The results of the present study showed that 56 glycosylated peptides were upregulated and 37 glycosylated peptides were downregulated in the feline small intestine infected by T. gondii. Additionally, we also identified eight N-glycosylated proteins of T. gondii including eight N-glycopeptides and eight N-glycosylation sites. The protein A0A086JND6_TOXGO (eEF2) and its corresponding peptide sequence were identified in T. gondii infection. Some special GO terms (i.e., cellular process and metabolic process, cell and cell part, and catalytic activity) were significantly enriched, and the Clusters of Orthologous Groups of proteins (COG) function prediction results showed that posttranslational modification, protein turnover, and chaperones (11%) had the highest enrichment for T. gondii. Interestingly, eEF2, a protein of T. gondii, is also involved in the significantly enriched T. gondii MAPK pathway. The host proteins ICAM-1 and PPT1 and the endoplasmic reticulum stress pathway may play an important role in the glycosylation of Toxoplasma-infected hosts. This is the first report showing that T. gondii oocysts can undergo N-glycosylation in the definitive host and that eEF2 is involved, which may provide a new target for T. gondii detection to prevent the spread of T. gondii oocysts in the future.

1. Introduction

Toxoplasma gondii is a strict intracellular parasite that can infect all warm-blooded vertebrates including humans and birds. Generally, infection with T. gondii is asymptomatic in healthy individuals but fatal for individuals with immunodeficiency [1]. For AIDS patients, T. gondii infection causes secondary encephalitis, pneumonia, and disseminated infection [2]; for pregnant women, the situation is more serious, and T. gondii infection can lead to premature birth, miscarriage, deformity, and even stillbirth [3,4]. Although T. gondii has a wide range of intermediate hosts, it only carries out sexual reproduction in felids [5]. Once cats are infected with T. gondii, bradyzoites settle within enterocytes, undergo sexual development, and form unsporulated oocysts [6]. The unsporulated oocysts are shed through the feces and then sporulated into infectious oocysts [7]. The discharged oocysts can be scattered on land, grass, water, or anywhere [8]. Contaminated water/grass/land, if eaten by livestock, can cause T. gondii infection in humans or companion animals [9]. T. gondii oocysts are a major public health problem.
Glycosylation is the most extensive and functional posttranslational modification (PTM) of proteins. Approximately 50% of proteins in eukaryotes undergo N-glycosylation modification and O-glycan [10]. The N-glycosylation of proteins is closely related to the occurrence and development of diseases such as diabetes, immune system diseases, nervous system diseases, and cancer [11,12]. Protein N-glycosylation plays a key role in numerous eukaryotic biological processes including molecular interactions, signal transduction, cell adhesion, endocytosis, protein stability, and other relevant biological functions [13,14]. In the biosynthesis of N-glycans, under the action of glycosidase, N-acetylglucosamine (GlcNAc) and mannose (Man) molecules are associated with asparagine (Asn) residues to form the basic structure of N-glycans. Galactose (Gal), N-acetylneuraminic acid (Neu5Ac), GalNAc, and other molecules participate in the process of sugar chain formation [10,15]. Asparagine-linked glycan (N-glycan)-dependent quality control (QC) systems for protein folding and endoplasmic reticulum (ER)-related degradation exist in most eukaryotes; however, in Giardia and some other protists, truncated N-glycans have been lost [16,17]. Plasmodium falciparum and T. gondii have also been found in lack proteins, which are involved in the N-glycan-dependent QC of protein folding [16]. Because the N-glycan structure of Plasmodium falciparum and T. gondii is controversial [18], and the T. gondii N-glycosylation function is critical for both parasite motility and host cell invasion [19], we were therefore very interested in whether N-glycosylation also occurs in the definitive host infection in T. gondii.
Although some articles highlight that the N-glycosylation of proteins plays an important role in the biological process of eukaryotes, this kind of PTM is rare and controversial in apicomplexan parasites [14]. However, some reports indicate that N-glycosylation is very common in the tachyzoite and cyst stages in T. gondii [14,20]. Studies have shown that the N-glycans of apicomplexans transmitted by arthropods are very short or absent, while the apicomplexans (i.e., Toxoplasma, Neosporidia, Cryptosporidium, and Eimeria) transmitted orally via oocysts have longer N-glycan precursors [18]. We speculated that the protein N-glycosylation of T. gondii is significant and indispensable when it invades its final host intestine. In this manuscript, we identified many known and novel candidate glycoproteins in the feline small intestine following T. gondii infection by using LC-MS/MS and further analyzed the pathway involvement of these glycoproteins.

2. Materials and Methods

2.1. Ethics Approval

The study design was reviewed and approved by the Animal Ethics Committee of Lanzhou Veterinary Research Institute (LVRI), Chinese Academy of Agricultural Sciences (Permit No. LVRIAEC2018-06). The procedures involving animals were carried out in accordance with the Animal Ethics Procedures and Guidelines of the People’s Republic of China. All efforts were made to minimize suffering and to reduce the number of cats used in the experiment.

2.2. Laboratory Animals

Kunming mice (8 weeks, 19–20 g, female, Lanzhou, China) were purchased from the laboratory animal center of LVRI, and had free access to sterile food and clean bottled water and were kept within a spacious cage at room temperature (25 ± 2 °C). The mice were intraperitoneally inoculated with T. gondii type II Prugniuad tachyzoites, and tissue cysts of the Pru strain were collected from the infected mice brain; then, it was crushed in a mortar, ground, diluted with phosphate buffered saline (PBS), and stored at 4 °C until use.
Six domestic cats (Felis catus, Chinese Li Hua breed, 7–9 months, Lanzhou, China) were purchased from a local breeder, and were randomly assigned to the experimental group (10_DPI_1, 10_DPI_2, and 10_DPI_3) and control group (Control), with 3 cats in each group. All of the experimental cats were confirmed to be negative for T. gondii by using the modified agglutination test (MAT; cut-off: 1:25) [21], and no other major feline diseases (i.e., feline immunodeficiency virus, feline leukemia virus, feline calicivirus, and feline parvovirus) were detected by commercial ELISA kits (Enzo, EastCoast Bio, Beijing, China) according to the manufacturer’s instructions. Each cat in the experimental group was intragastrically challenged with 600 Pru tissue cysts diluted in PBS, while the cats in the control groups were treated with the same amount of PBS. The cat’s intestinal epithelial tissue was scraped by a cell scraper at 10 days post-infection (DPI), flash frozen in liquid nitrogen, and then sent to the company BGI-Shenzhen (Shenzhen, China) for protein glycosylation analysis. Finally, the remaining laboratory animal tissues and T. gondii were subjected to sterilizing treatments such as autoclaving.

2.3. Confirmation of Infection Model Establishment

Genomic DNA was extracted from the harvested samples (30 mg) using a commercial kit (QIAamp DNA Mini Kit 50T, Cat. No. 1304, Dusseldorf, Germany) according to the manufacturer’s instructions. The PCR assay detection of T. gondii was performed using the B1 gene as described previously [22]. Parasites for histological observation were stained with H&E and observed by light microscopy (Olympus, BX41, Tokyo, Japan) [23].

2.4. Total Protein Extraction and Digestion

The sample (15–20 mg) was mixed with 5 mm magnetic beads (BeaverBeads™ his-tag protein purification, Solarbio, Beijing, China) and lysis buffer 3 (containing 1 mM PMSF, 2 mM EDTA, and 10 mM dithiothreitol) and centrifuged to obtain the supernatant, and then 10 mM dithiothreitol was added and the sample was incubated in a water bath at 56 °C for 1 h. After adding 55 mM IAM at room temperature (avoiding light incubation for 45 min), 4 times the volume of acetone (−20 °C, 2 h) was added, and the same volume of acetone (−20 °C, 2 h) was added again two to three times until the supernatant was colorless. The precipitate was collected after centrifugation, mixed with 5 mm magnetic beads with the lysis buffer 3, shaken for 2 min, centrifuged, and used for quantification. Quantification of the protein was carried out using the Bradford Kit (Bradford Assay Kit, Bio-Rad, Shanghai, China) [24] and analyzed by SDS–PAGE [25] as described in the manufacturer’s protocol.
A total of 500 μg of the protein solution was taken from each sample, and the enzyme hydrolysis activity was determined according to the protein solution: trypsin enzyme = 40:1 (37 °C, 4 h). A Strata X column was utilized to remove salt and then the enzymatic peptides were subjected to vacuum drying. A small portion of the peptides was used to detect the enzymatic hydrolysis effect using mass spectrometry.

2.5. Enrichment of N-glycosylated Peptides and Deglycosylation

Three hundred microliters of qualified peptides were dissolved in 60% ACN and 0.1% FA, and the optimized HILIC method was then used for enrichment and grouping fractionation by HPLC (5 µm, 150 × 4.6 mm, Merck, Beijing, China) [26]. Fractionated product collection started at the 30th min and ended at the 54th min, with one tube per minute. A total of 24 tubes were collected and labeled and then pumped by a freezer pump. Each of the 4 adjacent tubes were reconstituted with a total of 50 µL of 50 mM NH4HCO3, combined into one tube, and finally divided into 6 fragments. Then, 2.5 µL of PNGase F was added to each fragment, vortexed, mixed well, centrifuged immediately, and incubated overnight at 37 °C. Then, it was frozen and drained.

2.6. High-Performance Liquid Chromatography Coupled to Tandem Mass Spectrometry

The extracted peptide samples were redissolved in mobile phase A (2% ACN, 0.1% FA) and centrifuged at 20,000× g for 10 min, and then the supernatant was used for gradient separation by UHPLC (UltiMate 3000, Thermo, Shanghai, China). The separated peptides were ionized by the nanoESI source and transferred to the tandem mass spectrometer Q-Exactive HF (Thermo Fisher Scientific, San Jose, CA, USA) for data-dependent acquisition (DDA) mode detection [27]. Main parameter settings: The MS1 scan range was 350~1500 m/z, and the resolution was set to 120,000. The MS2 started at 100 m/z, and the resolution was 30,000. The ion fragmentation mode was HCD, and the dynamic exclusion time was set to 30 s. The AGC settings: 1E6 for the first level, and 5E4 for the second level.

2.7. Enrichment of N-glycosylated Peptides and Deglycosylation

The host proteins/peptides analysis process: The FFPE sample data generated by the high-resolution mass spectrometer were identified using MaxQuant (Andromeda, v.1.5.3.30, http://www.maxquant.org, 7 February 2019) [28]; the fixed modifications were carbamidomethyl (C) and the variable modifications were oxidation (M), acetyl (protein N-term), and Asn- > Asp (N). The filtration was completed with PSM-level FDR <= 1% at the spectrum level and filtered with protein level FDR <= 1% at the protein level, and the database was searched by the UniProt Felis catus and Toxoplasma gondii (111,780 sequences). At the same time, the modified sites were filtered with 1% FDR to obtain significant modification results. MaxQuant was used to perform the quantitative analysis based on the peak strength, peak area, retention time, and other information regarding the peptide segment related to the primary mass spectrometry, extracting the peptide segment related to the target modification, and conducting a series of statistical analyses and quality controls. The Welch’s t-test (p-value) was used to test the significance for MaxQuant, and the modification quantification was screened according to a fold change ≥ 2 and a p-value ≤ 0.05 as the criteria for a significant difference.
The parasite proteins/peptides analysis process: The identification of N-glycosylated proteins was mainly based on the matching of experimental tandem mass spectrometry data with theoretical mass spectrometry data obtained by the database simulation to obtain protein identification results. The protein identification was performed using Mascot software (Mascot 2.3.02, https://mascot-distiller.software.informer.com/, 1 November 2021). The fixed modifications were carbamidomethyl (C), and the variable modifications were oxidation (M), acetyl (protein N-term), and Asn -> Asp (N). The filtration was completed with the peptide mass tolerance set at 20 ppm and filtered with the fragment mass tolerance set at 0.05 Da. The max missed cleavages was 2, and the database was searched using the UniProt Toxoplasma gondii library (77,541 sequences, https://www.uniprot.org/, 8 January 2019). This process used Percolator (http://noble.gs.washington.edu/proj/percolator, 1 November 2017) to preprocess and re-score the results generated by the search engine to improve the accuracy of the correct matching and random matching. The output results were then filtered with PSM-level FDR <= 0.01 to obtain a list of significantly identified spectra and peptides. Then, based on the “parsimony principle”, this process performed protein inference on the peptides and generated a series of proteomes.
In this study, the Universal Protein Resource (UniProt) database (https://www.uniprot.org/, 1 April 2022) was used to annotate the identified proteins. The Kyoto Encyclopedia of Genes and Genomes (KEGG) Orthology-Based Annotation System 3.0 (KOBAS) (http://bioinfo.org/kobas, 20 December 2020) is a database resource for pathway annotation. KEGG is the main public database of pathways, and pathway analysis can identify the major biochemical metabolic pathways and signal transduction pathways in which proteins are involved [29]. The Gene Ontology (GO, http://geneontology.org/, 7 October 2022) database was used for the gene enrichment analyses, and GO function annotation analysis was performed on all of the identified proteins. The results include protein2go (providing a list of IDs of all corresponding GO functions for each protein) and go2protein (i.e., biological process, BP; cellular component, CC; molecular function, MF). The Clusters of Orthologous Groups of proteins (COGs) database (https://www.ncbi.nlm.nih.gov/research/cog/, 1 March 2022) was used for the orthologous classification of the proteins. In addition, the Minimal Common Oncology Data Elements (MCODE) plugin for Cytoscape software (https://cytoscape.org/, 21 April 2021) was applied to analyze the highly interconnected clusters in the protein–protein interaction (PPI) network using the default parameters.

3. Results

3.1. Differential Quantitative Analysis of the Modified Peptides/Proteins

After challenge with T. gondii (TOX) in the small intestine of cats, the DNA detection of the B1 gene showed that T. gondii-infected cat models were successfully established (Supplementary Materials Figure S1). The pathological tissue section of the cats’ small intestines also illustrated the successful establishment of infection. As shown in Supplementary Materials Figure S2, compared with the control group, the experimental group showed swelling of the muscle fibers in the muscle layer (Supplementary Materials Figure S2A). Supplementary Materials Figure S2B shows that the number of intestinal microvilli in the experimental group decreased sharply due to the stromal hyperplasia squeezing the intestinal microvilli. In addition, the experimental group also had basal intestinal hemorrhage accompanied by muscle tissue hyperplasia. The diffuse hemorrhage of the intestinal villi was also found in the infected small intestine (Supplementary Materials Figure S2C). All of these results indicate that the model of T. gondii infection in the small intestine of cats was successfully established and could be applied to further study.
In this study, a total of 12,202 proteins, 3122 (Felis: 2656, TOX: 466) modified peptides, and 1813 (Felis: 1411, TOX: 402) modified proteins were identified by analyzing the glycosylation results of the T. gondii-infected small intestine of cats (the details are listed in Supplementary Materials Table S1). The mass distribution of the identified modified proteins and the number of unique peptide are shown in Figure 1A,B, and the motif distribution of the posttranslational modification sites of the identified modified peptides is shown in Figure 1C. The Welch’s t-test (fold change ≥ 2 and q-value < 0.05) was utilized to identify the differentially glycosylated peptides between the experimental group and the control group. Fifty-eight upregulated and 37 downregulated modified peptides with significant differences were detected in the infected small intestine (Figure 2 and Supplementary Materials Table S2). The two most significantly upregulated genes (M3XF94_FELCA and M3W4C1_FELCA) and the two most significantly downregulated genes (A0A2I2U7X3_FELCA and A0A2I2UHR2_FELCA) are shown in Figure 2. In addition, we identified 8 N-glycosylated T. gondii proteins including 8 N-glycosylated peptides and 8 N-glycosylation sites (Table 1 and Supplementary Materials Table S3). The identification map of N-glycosylation peptides is shown in Supplementary Materials Figure S3, and the distribution of N-glycosylation motifs is shown in Supplementary Materials Figure S4. An interesting finding is that the protein A0A086JND6_TOXGO (putative translation elongation factor 2 family protein, eEF2) was identified in groups 10_DPI_1, 10_DPI_2, and 10_DPI_3, and the protein A0A086K8H6_TOXGO (ribosomal protein RPL11) was identified in 10_DPI_2 and 10_DPI_3. Another meaningful discovery was that the peptide sequence NMSVIAHVDHGK (A0A086JND6_TOXGO) was identified in 10_DPI_1 (once), 10_DPI_2 (twice), and 10_DPI_3 (twice), and the peptide sequence KKNFSDSGNFGFGIQEHIDLGIK (A0A086K8H6_TOXGO) was identified once in 10_DPI_2 and 10_DPI_3 (Additional Table S5). Even more coincidentally, we found that peptide_seq NMSVIAHVDHGK was not only a part of T. gondii protein A0A086JND6_TOXGO but also a part of protein A0A139Y075_TOXGO and protein A0A2G8Y563_TOXGO.

3.2. Functional Annotation and Enrichment Analysis of the Proteins with Differentially Modified Peptides

A total of 42 GO terms, including 19 biological process (BP) terms, 15 cellular component (CC) terms, and 8 molecular function (MF) terms, were significantly enriched with the proteins corresponding to 95 (58 up- and 37 downregulated) differentially modified proteins in the small intestine of cats (Figure 3). In the BP category, the top two upregulated enriched GO terms were single-organism process and cellular process, and the top two downregulated enriched GO terms were cellular process and metabolic process. In the CC category, membrane and cell/cell part were the top two upregulated GO terms, and the top two downregulated enriched GO terms were cell and cell part. In the MF category, the top two upregulated/downregulated GO terms were catalytic activity and binding. In addition to these common GO terms, some special BP terms (biological regulation and response to stimulus), CC terms (cell junction), and MF terms (signal transducer activity, structural molecule activity, and molecular function regulator) were also enriched. Interestingly, there were also many T. gondii-modified proteins involved in these pathways, for example, cellular process and metabolic process, cell and cell part, and catalytic activity.
The 1211 proteins were divided into 24 functional classes, among which “general function prediction only” accounted for 14% (166/1211), “posttranslational modification, protein turnover, and chaperones” (11%, 129/1211), and “cell wall/membrane/envelope biogenesis” (10%, 117/1211) accounted for the top three. We also found that in these host-related COG classifications, some of the T. gondii proteins were also involved in eight classifications, but the proteins involved are relatively rare (Figure 4). “General function prediction only” included the WD40 repeat, secreted protein containing bacterial Ig-like domain, and vWFA domain, predicted pyrophosphatase or phosphodiesterase, alkaline phosphatase (AlkP) superfamily, uncharacterized conserved protein YfaS, alpha-2-macroglobulin family, tetratricopeptide repeat, uncharacterized conserved protein RhaS, contains 28 RHS repeats, GTPase SAR1 family domain, DNA-binding beta-propeller fold protein YncE, predicted ATPase, etc. The function of “cell wall/membrane/envelope biogenesis” is numerous and mainly involves membrane proteins involved in colicin uptake, murein tripeptide amidase MpaA, O-glycosyl hydrolase, and the PASTA domain, which binds beta-lactams. The details are presented in Supplementary Materials Table S4. The “Posttranslational modification, protein turnover, and chaperones”-specific functions included the molecular chaperone DnaK (HSP70 and A0A2G8YAT4_TOXGO), molecular chaperone (HSP90 family and B9QJS4_TOXGV), peptidyl-prolyl cis-trans isomerase–rotamase (cyclophilin family and Q26994_TOXGO), chaperonin groEL (HSP60 family), zn-dependent amino- or carboxypeptidase (M28 family), serine protease inhibitor, asparagine N-glycosylation enzyme (membrane subunit Stt3), cysteine protease (C1A family), serine protease (subtilase family), and glutathione s-transferase.
The identified modified peptides in this study were mainly annotated into six major subsystems (i.e., metabolism, genetic information processing, environmental information processing, cellular processes, organismal systems, and human diseases), as shown in Figure 5. The most commonly involved pathways were signal transduction, transport and catabolism, immune system, and infectious diseases. Figure 6 shows the top 30 enriched KEGG pathways of the host, and the top three significantly enriched pathways were African trypanosomiasis, fatty acid elongation, and primary immunodeficiency. Meanwhile, the KEGG enrichment results of the T. gondii-related protein are shown in Table 2; among them, A0A086JND6_TOXGO participates in the AMP-activated protein kinase (AMPK) signaling pathway and the oxytocin signaling pathway.

3.3. PPI Networks for Differentially Modified Peptides/Proteins

The PPI networks were established of the differentially modified peptides/proteins using the STRING database. We found that a variety of interactions in the host cells were up-/downregulated at the glycosylation level at 10 days after T. gondii infection. Then, we applied cytotype MCODE for further analysis of the highly relevant subnetworks and identified three MCODE clusters (i.e., viral myocarditis, fatty acid elongation, and protein processing in endoplasmic reticulum) (Figure 7). These analysis results were consistent with those of the KEGG enrichment analysis. The proteins ICAM1 (M3WBF1_FELCA) and CDH17 (M3WMR7_FELCA) were upregulated, while ITGAL (M3 WC52_FELCA) and ITGAE (M3W8N0_FELCA) were downregulated in the viral myocarditis pathway. In the fatty acid elongation pathway, we found that RPL3 (M3W4R9_ FELCA) was upregulated, while PPT1 (M3X4P7_ FELCA) and PPT2 (M3W2C6_ FELCA) were downregulated. For protein processing in the endoplasmic reticulum pathway, in addition to the downregulation of EDEM3 (M3W7T1_FELCA), the other involved proteins, COPB2 (M3VWP9_FELCA), SEL1L (M3WAN5_FELCA), LMAN2 (M3XF94_FELCA), and ERGIC3 (M3VZZ0_FELCA), were all upregulated.

4. Discussion

Glycosylation is an important posttranslational modification. Most membrane proteins and secretory proteins are glycoproteins [30]. Glycosylation modification not only affects the spatial conformation, activity, transportation, and localization of proteins but also plays a critical role in signal transduction, molecular recognition, and immunity [31,32]. In this study, we performed a global glycosylation modification analysis of T. gondii-infected small intestines of cats after 10 DPI by using LC-MS/MS technology. A total of 1813 modified proteins (Felis: 1411, TOX: 402) were identified in the glycosylation analysis after T. gondii infected the terminal host, indicating that glycosylation of T. gondii may occur frequently in the process of infection with definitive hosts, which is consistent with previous reports of the glycosylation of T. gondii [14]. A total of 3122 (Felis: 2656, TOX: 466) modified peptides were identified. Meanwhile, under strict screening conditions, we found that there were more upregulated differentially modified peptides (56) than downregulated ones (37), which may mean that most of the cat-modified peptides were upregulated after T. gondii infection.
Among these findings, eight N-glycosylated proteins and eight N-glycopeptides of T. gondii were screened. We infer that the eEF2 protein (A0A086JND6_TOXGO) and its corresponding peptide sequence (NMSVIAHVDHGK) may be meaningful in future research. We noted that eEF2 plays a major role in the protein synthesis and survival of Eimeria tenella [33], and low levels of total eEF2 protein in T. gondii reduced TgPHYb expression [34]. eEF2 has been identified or confirmed as an antimalarial and anticoccidial drug target of Plasmodium falciparum and Eimeria tenella, respectively [33,34]. It is clear that the eEF2 of T. gondii localizes to the nucleus [35]. Based on these findings, we speculate that the eEF2 protein and its corresponding peptide sequences (NMSVIAHVDHGK) may be quite unexpected candidate glycoproteins for T. gondii control, and its glycosylation status and contribution to T. gondii and host will await further experimental investigations.
The results of the GO analysis of the glycosylated proteins after T. gondii infection in cats suggested that several special GO terms, such as biological regulation and response to stimulus, that play important roles in the NF-κB signaling pathway and the regulation of MIC protein were enriched in T. gondii and the cats [36,37]. The cell junction term plays an important role in maintaining the invasive force of T. gondii [38]. Signal transducer activity, structural molecule activity, and molecular function regulator terms are linked to the process by which ROP16 is released from rhoptries into the host cell [39,40]. The GO terms cellular process and metabolic process, cell and cell part, and catalytic activity were not only involved in the pathway enrichment of the host but were also a significant pathway of T. gondii, which indicates that after T. gondii infects the host, some pathways changed significantly and synchronously between the host and T. gondii.
The COG potential function prediction showed that the identified glycosylated proteins were mainly enriched in three parts: “general function prediction only” (14%); “posttranslational modification, protein turnover, and chaperones” (11%); “cell wall/membrane/envelope biogenesis” (10%); moreover, each of them contained many functional proteins. For example, in the “posttranslational modification, protein turnover, and chaperones” section, the HSP70 and HSP90 families are not only molecular chaperones [41] but also glycosylated proteins of T. gondii identified in this study. They are related to T. gondii differentiation, which is an important process in T. gondii pathology [42]. Taken together, these findings indicate that the glycosylation of Hsp90/Hsp70 could act as a key process in the life of T. gondii. The cyclophilins of peptidyl-prolyl cis-trans isomerases (PPIases) are believed to be involved in protein folding [43], and their role in T. gondii may be regulated by P glycoprotein [44]. Heat shock protein 60 (TgHSP60) plays an important role in intracellular survival and in the differentiation of T. gondii, and it has been isolated from cats [45]. In our study, HSP60 was identified just after the cat was infected with T. gondii, which further shows that HSP60 does play an important role in the differentiation of oocysts after cat infection. Serine protease inhibitors display immunomodulatory properties [46], and T. gondii serine protease inhibitor 1 (TgPI1) not only affects the virulence of T. gondii but also affects the differentiation of bradyzoites [47]. It has been speculated that this may occur through inhibition of T. gondii and/or host serine protease [47]. In our study, serine protease inhibitors were detected in the host glycosylated protein after T. gondii infection in cats, which also confirmed this inference from another aspect.
The KEGG enrichment results for the glycosylated proteins showed that the signal transduction, transport, and catabolism; immune system; infectious diseases were significantly enriched, indicating that these pathways play an important role in the glycosylation process of the host. In addition, fatty acid elongation and primary immunodeficiency pathways were the top significantly enriched pathways, and these two pathways are closely related to cat glycosylation. For example, acetyl-CoA synthetase is essential for the fatty elongation pathway to generate fatty acids used for T. gondii membrane biogenesis [48], and endoplasmic reticulum-localized fatty acid elongation and very long-chain unsaturated fatty acids are essential for T. gondii growth but are not supplied by the host cell [49]. In addition, the AMPK signaling pathway and oxytocin signaling pathway were significantly enriched during the process of T. gondii infection, and both are involved in the T. gondii A0A086K8H6_TOXGO protein. AMPK signaling is required for CD40-induced autophagic killing of T. gondii and plays a pivotal role in regulating cellular energy homeostasis [49]. Based on this, we can speculate that the AMPK signaling pathway may significantly change the cell’s metabolism after T. gondii infects the definitive host.
Three MCODE_Clusters (i.e., viral myocarditis, fatty acid elongation, and protein processing in endoplasmic reticulum) of Cytoscape by the differentially expressed peptides/proteins are shown in Figure 7. Regarding the viral myocarditis pathway, T. gondii can cause severe myocarditis and even death when people have immune deficiency or immune insufficiency [50,51]. Intercellular adhesion molecule 1 (ICAM-1) is not only upregulated on cellular barriers when the host is infected but is also related to T. gondii, which exploits the host to cross cellular barriers and disseminate to deep tissues [52]. In this study, ICAM-1 was also upregulated, and we speculated that this upregulation may be linked to the host migration of oocysts. T. gondii possesses a single cytosolic acetyl-CoA synthetase (TgACS) that is involved in providing acetyl-CoA for the essential fatty elongation pathway to generate fatty acids used for membrane biogenesis [48], and it was found that T. gondii oocyst walls were acid fast and contained long fatty acyl chains that might be synthesized by an abundant polyketide synthase [53]. We assume that this may be because the host deliberately coats the oocyst wall with acid-fast lipids to make them resistant to the environment. T. gondii could activate the unfolded protein response upon endoplasmic reticulum stress, and this pathway could not only cause activation of autophagy but also induce phosphorylation of Toxoplasma gondii eIF2α and inhibit translation initiation [54,55]. PPT1 of T. gondii is based on amino acid residues distal to the palmitoyl cysteine, which has distinct specificity and could be used to detect the activity of a specific depalmitoylase in complex proteomes [56]. All of these highlighted pathways were found in the differentially glycosylated proteins, whether the alteration contributed to the migration or formation of T. gondii oocysts in the cat’s small intestine remains to be determined.

5. Conclusions

T. gondii oocysts cause great harm to humans and the surrounding environment. In this study, we found that N-glycosylation occurs in infected cats and T. gondii, and the N-glycosylation of eEF2 in T. gondii may be an important clue for studying the sexual reproduction of T. gondii. The results of the COG and GO analyses showed that the “posttranslational modification, protein turnover, and chaperones” pathways were significantly enriched during the glycosylation process of T. gondii-infected cats. This proteome data provides a landscape of N-glycosylation in the protein of T. gondii and the infected cat.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/ani12202858/s1, Figure S1. The results of the PCR for the samples. Identification of T. gondii-infected cat small intestine samples at 10 DPI by the T. gondii B1 gene. The order of the sample holes was M: Takara DNA marker (500 bp)—E-group (lanes 2–4) and C-group (lanes 6–8); P: T. gondii Pru strain PCR positive control; N: PCR negative control. E-group: experimental group; C-group: control group. Figure S2. H&E staining of the samples. The pathological results of the cats’ small intestine infected with T. gondii at 10 DPI. The arrow points to the lesion. The experimental group showed swelling of the muscle fibers in the muscle layer (A). (B) The number of intestinal microvilli in the experimental group decreased sharply due to the stromal hyperplasia squeezing the intestinal microvilli. In addition, the experimental group (C) also had basal intestinal hemorrhage accompanied by muscle tissue hyperplasia, and diffuse hemorrhage of the intestinal villi was also found in the infected small intestine. Figure S3. Identification map of the N-glycosylated peptides. Distribution map of the number of identifications of each peptide map (A: 10_DPI_1, B: 10_DPI_2, and C: 10_DPI_3). The x-axis is the number of spectrum supports corresponding to each N-glycosylated peptide, and the y-axis is the number of N-glycosylation peptides. The identification of the N-glycosylated peptides was performed by searching and identifying each spectrum generated by mass spectrometry. For a peptide, the more spectrums identified, the higher the confidence of the peptide. Figure S4. Distribution map of the N-glycosylation motifs—A: 10_DPI_1, B: 10_DPI_2, and C: 10_DPI_3. The x-axis represents the base number, and the y-axis represents the corrected score. The higher the base height, the higher the probability of the base appearing in the motif. Table S1: List of all modified peptides/protein results. Table S2: List of all modified peptides with significant differences. Table S3: List of 8 N-glycosylated peptides and 8 N-glycosylation sites. Table S4: List of the COG classification results of all modified proteins.

Author Contributions

J.Z., J.H. and X.Z. conceived and designed the study and critically revised the manuscript; B.Z. performed the experiment, analyzed the data, and drafted the manuscript; S.X., J.P. and Y.Q. helped in the study implementation and data analysis; Y.L. analyzed the histopathological results. All authors have read and agreed to the published version of the manuscript.

Funding

Project support was kindly provided by the National Natural Science Foundation of China (Grant No. 32102701), the Central Public-Interest Scientific Institution Basal Research Fund (Grant No. 1610322022010), the Innovation Project of the Chinese Academy of Agricultural Sciences (Grant No. 25-LZIHPS-05), the Agricultural Science and Technology Innovation Program (ASTIP) (Grant No. CAAS-ASTIP-2016-LVRI-03), the Yunnan Expert Workstation (Grant No. 202005AF150041), the Fund for Shanxi “1331 Project” (Grant No. 20211331-13), and the Youth Science and Technology Fund Program of Gansu Province (Grant No. 20JR5RA581).

Institutional Review Board Statement

The study design was reviewed and approved by the Animal Ethics Committee of Lanzhou Veterinary Research Institute (LVRI), Chinese Academy of Agricultural Sciences (Permit No. LVRIAEC2018-06).

Informed Consent Statement

Not applicable.

Data Availability Statement

The mass spectrometry data obtained in this study were deposited in the National Center for Biotechnology Information (NCBI) Sequence Read Archive (SRA) database (https://www.ncbi.nlm.nih.gov/sra), under accession number PRJNA673972, or in Mendeley Data (https://data.mendeley.com/datasets/rn2m56vdky/1).

Acknowledgments

We thank BGI-Shenzhen for the technical assistance.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Yarovinsky, F. Innate immunity to Toxoplasma gondii infection. Nat. Rev. Immunol. 2014, 14, 109–121. [Google Scholar] [CrossRef] [PubMed]
  2. Wang, Z.D.; Wang, S.C.; Liu, H.H.; Ma, H.Y.; Li, Z.Y.; Wei, F.; Zhu, X.Q.; Liu, Q. Prevalence and burden of Toxoplasma gondii infection in HIV-infected people: A systematic review and meta-analysis. Lancet HIV 2017, 4, e177–e188. [Google Scholar] [CrossRef]
  3. Moncada, P.A.; Montoya, J.G. Toxoplasmosis in the fetus and newborn: An update on prevalence, diagnosis and treatment. Expert. Rev. Anti-Infect. Ther. 2014, 10, 815–828. [Google Scholar] [CrossRef] [PubMed]
  4. Smith, N.C.; Goulart, C.; Hayward, J.A.; Kupz, A.; Miller, C.M.; van Dooren, G.G. Control of human toxoplasmosis. Int. J. Parasitol. 2021, 51, 95–121. [Google Scholar] [CrossRef]
  5. Dubey, J.P.; Cerqueira-Cezar, C.K.; Murata, F.; Kwok, O.; Yang, Y.R.; Su, C. All about toxoplasmosis in cats: The last decade. Vet. Parasitol. 2020, 283, 109145. [Google Scholar] [CrossRef]
  6. Di Genova, B.M.; Wilson, S.K.; Dubey, J.P.; Knoll, L.J. Intestinal delta-6-desaturase activity determines host range for Toxoplasma sexual reproduction. PLoS Biol. 2019, 17, e3000364. [Google Scholar]
  7. Zulpo, D.L.; Sammi, A.S.; Dos, S.J.R.; Sasse, J.P.; Martins, T.A.; Minutti, A.F.; Cardim, S.T.; de Barros, L.D.; Navarro, I.T.; Garcia, J.L. Toxoplasma gondii: A study of oocyst re-shedding in domestic cats. Vet. Parasitol. 2018, 249, 17–20. [Google Scholar] [CrossRef]
  8. Shapiro, K.; Bahia-Oliveira, L.; Dixon, B.; Dumetre, A.; de Wit, L.A.; VanWormer, E.; Villena, I. Environmental transmission of Toxoplasma gondii: Oocysts in water, soil and food. Food Waterborne Parasitol. 2019, 15, e00049. [Google Scholar] [CrossRef]
  9. Alizadeh, A.M.; Jazaeri, S.; Shemshadi, B.; Hashempour-Baltork, F.; Sarlak, Z.; Pilevar, Z.; Hosseini, H. A review on inactivation methods of Toxoplasma gondii in foods. Pathog. Glob. Health 2018, 112, 306–319. [Google Scholar] [CrossRef]
  10. Gomord, V.; Fitchette, A.C.; Menu-Bouaouiche, L.; Saint-Jore-Dupas, C.; Plasson, C.; Michaud, D.; Faye, L. Plant-specific glycosylation patterns in the context of therapeutic protein production. Plant Biotechnol. J. 2010, 8, 564–587. [Google Scholar] [CrossRef]
  11. Lauc, G.; Huffman, J.E.; Pučić, M.; Zgaga, L.; Adamczyk, B.; Mužinić, A.; Novokmet, M.; Polašek, O.; Gornik, O.; Krištić, J.; et al. Loci associated with N-glycosylation of human immunoglobulin G show pleiotropy with autoimmune diseases and haematological cancers. PLoS Genet. 2013, 9, e1003225. [Google Scholar] [CrossRef]
  12. Hwang, H.; Zhang, J.P.; Chung, K.A.; Leverenz, J.B.; Zabetian, C.P.; Peskind, E.R.; Jankovic, J.; Su, Z.; Hancock, A.M.; Pan, C.; et al. Glycoproteomics in neurodegenerative diseases. Mass Spectrum. Rev. 2010, 29, 79–125. [Google Scholar] [CrossRef] [Green Version]
  13. Zaia, J. Mass spectrometry and the emerging field of glycomics. Chem. Biol. 2008, 15, 881–892. [Google Scholar] [CrossRef] [Green Version]
  14. Luk, F.C.; Johnson, T.M.; Beckers, C.J. N-linked glycosylation of proteins in the protozoan parasite Toxoplasma gondii. Mol. Biochem. Parasitol. 2008, 157, 169–178. [Google Scholar] [CrossRef] [Green Version]
  15. Nazarova, L.A.; Ochoa, R.J.; Jones, K.A.; Morrissette, N.S.; Prescher, J.A. Extracellular Toxoplasma gondii tachyzoites metabolize and incorporate unnatural sugars into cellular proteins. Microbes Infect. 2016, 18, 199–210. [Google Scholar] [CrossRef] [Green Version]
  16. Banerjee, S.; Vishwanath, P.; Cui, J.; Kelleher, D.J.; Gilmore, R.; Robbins, P.W.; Samuelson, J. The evolution of N-glycan-dependent endoplasmic reticulum quality control factors for glycoprotein folding and degradation. Proc. Natl. Acad. Sci. USA 2007, 104, 11676–11681. [Google Scholar] [CrossRef] [Green Version]
  17. Helenius, A.; Aebi, M. Roles of N-linked glycans in the endoplasmic reticulum. Annu. Rev. Biochem. 2004, 73, 1019–1049. [Google Scholar] [CrossRef]
  18. Bushkin, G.G.; Ratner, D.M.; Cui, J.; Banerjee, S.; Duraisingh, M.T.; Jennings, C.V.; Dvorin, J.D.; Gubbels, M.J.; Robertson, S.D.; Steffen, M.; et al. Suggestive evidence for darwinian selection against asparagine-linked glycans of Plasmodium falciparum and Toxoplasma gondii. Eukaryot. Cell 2010, 9, 228–241. [Google Scholar] [CrossRef] [Green Version]
  19. Fauquenoy, S.; Hovasse, A.; Sloves, P.J.; Morelle, W.; Alayi, T.D.; Slomianny, C.; Werkmeister, E.; Schaeffer, C.; Dorsselaer, A.V.; Tomavo, S. Unusual N-glycan structures required for trafficking Toxoplasma gondii GAP50 to the inner membrane complex regulate host cell entry through parasite motility. Mol. Cell. Proteom. 2011, 10, M111.008953. [Google Scholar]
  20. Gaskins, E.; Gilk, S.; DeVore, N.; Mann, T.; Ward, G.; Beckers, C. Identification of the membrane receptor of a class XIV myosin in Toxoplasma gondii. J. Cell Biol. 2004, 165, 383–393. [Google Scholar] [CrossRef] [Green Version]
  21. Rani, S.; Cerqueira-Cezar, C.K.; Murata, F.; Sadler, M.; Kwok, O.; Pradhan, A.K.; Hill, D.E.; Urban, J.J.; Dubey, J.P. Toxoplasma gondii tissue cyst formation and density of tissue cysts in shoulders of pigs 7 and 14 days after feeding infected mice tissues. Vet. Parasitol. 2019, 269, 13–15. [Google Scholar] [CrossRef]
  22. Ma, J.; He, J.J.; Hou, J.L.; Zhou, C.X.; Zhang, F.K.; Elsheikha, H.M.; Zhu, X.Q. Metabolomic signature of mouse cerebral cortex following Toxoplasma gondii infection. Parasit. Vectors 2019, 12, 373. [Google Scholar] [CrossRef]
  23. Zhou, C.X.; Zhou, D.H.; Elsheikha, H.M.; Liu, G.X.; Suo, X.; Zhu, X.Q. Global metabolomic profiling of mice brains following experimental infection with the cyst-forming Toxoplasma gondii. PLoS ONE 2015, 10, e139635. [Google Scholar]
  24. Wu, D.; Xue, X.W.; Gao, C.C.; Liu, Y.H.; Wang, T.T.; Li, L.S.; Tong, X.H.; Li, F.; Xu, J.D. Rhubarb-Evoke mucus secretion through aggregation and degranulation of mast cell in the colon of rat: In vivo and ex vivo studies. Sci. Rep. 2019, 9, 19375. [Google Scholar] [CrossRef] [Green Version]
  25. Haider, S.R.; Reid, H.J.; Sharp, B.L. Tricine-SDS-PAGE. Methods Mol. Biol. 2019, 1855, 151–160. [Google Scholar]
  26. Zappacosta, F.; Scott, G.F.; Huddleston, M.J.; Annan, R.S. An optimized platform for hydrophilic interaction chromatography-immobilized metal affinity chromatography enables deep coverage of the rat liver phosphoproteome. J. Proteome Res. 2015, 14, 997–1009. [Google Scholar] [CrossRef]
  27. Cole, J.; Hanson, E.J.; James, D.C.; Dockrell, D.H.; Dickman, M.J. Comparison of data-acquisition methods for the identification and quantification of histone post-translational modifications on a Q Exactive HF hybrid quadrupole Orbitrap mass spectrometer. Rapid Commun. Mass Spectrom. 2019, 33, 897–906. [Google Scholar] [CrossRef] [Green Version]
  28. Cox, J.; Mann, M. MaxQuant enables high peptide identification rates, individualized p.p.b.-range mass accuracies and proteome-wide protein quantification. Nat. Biotechnol. 2008, 26, 1367–1372. [Google Scholar] [CrossRef]
  29. Kanehisa, M.; Araki, M.; Goto, S.; Hattori, M.; Hirakawa, M.; Itoh, M.; Katayama, K.; Kawashima, S.; Okuda, S.; Tokimatsu, T.; et al. KEGG for linking genomes to life and the environment. Nucleic Acids Res. 2008, 36, D480–D484. [Google Scholar] [CrossRef]
  30. Pedrazzini, E.; Caprera, A.; Fojadelli, I.; Stella, A.; Rocchetti, A.; Bassin, B.; Martinoia, E.; Vitale, A. The Arabidopsis tonoplast is almost devoid of glycoproteins with complex N-glycans, unlike the rat lysosomal membrane. J. Exp. Bot. 2016, 67, 1769–1781. [Google Scholar] [CrossRef] [Green Version]
  31. Wiederschain, G.Y. Essentials of Glycobiology; Cold Spring Harbor Laboratory Press: New York, NY, USA, 2009. [Google Scholar]
  32. Moremen, K.W.; Tiemeyer, M.; Nairn, A.V. Vertebrate protein glycosylation: Diversity, synthesis and function. Nat. Rev. Mol. Cell Biol. 2012, 13, 448–462. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Zhou, B.H.; Jia, L.S.; Guo, H.W.; Ding, H.Y.; Yang, J.Y.; Wang, H.W. Eukaryotic elongation factor 2 is involved in the anticoccidial action of diclazuril in the second-generation merozoites of Eimeria tenella. Vet. Parasitol. 2019, 276, 108991. [Google Scholar] [CrossRef] [PubMed]
  34. Florimond, C.; Cordonnier, C.; Taujale, R.; van der Wel, H.; Kannan, N.; West, C.M.; Blader, I.J. A toxoplasma prolyl hydroxylase mediates oxygen stress responses by regulating translation elongation. mBio 2019, 10, e00234-19. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Barylyuk, K.; Koreny, L.; Ke, H.; Butterworth, S.; Crook, O.M.; Lassadi, I.; Gupta, V.; Tromer, E.; Mourier, T.; Stevens, T.J.; et al. Comprehensive subcellular atlas of the Toxoplasma proteome via hyperLOPIT provides spatial context for protein functions. Cell Host Microbe 2020, 28, 752–766.e9. [Google Scholar] [CrossRef]
  36. Liu, Y.M.; Zou, X.Q.; Ou, M.L.; Ye, X.H.; Zhang, B.H.; Wu, T.Y.; Dong, S.R.; Chen, X.J.; Liu, H.X.; Zheng, Z.; et al. Toxoplasma gondii cathepsin C1 inhibits NF-kappaB signalling through the positive regulation of the HIF-1alpha/EPO axis. Acta Trop. 2019, 195, 35–43. [Google Scholar] [CrossRef]
  37. Mendonça-Natividade, F.C.; Duque, L.C.; Ricci-Azevedo, R.; Sardinha-Silva, A.; Figueiredo, P.C.; Paiva, A.A.; Alegre-Maller, A.C.P.; Nohara, L.L.; Carneiro, A.B.; Panunto-Castelo, A.; et al. Receptor heterodimerization and co-receptor engagement in TLR2 activation induced by MIC1 and MIC4 from Toxoplasma gondii. Int. J. Mol. Sci. 2019, 20, 5001. [Google Scholar] [CrossRef] [Green Version]
  38. Bichet, M.; Joly, C.; Henni, A.H.; Guilbert, T.; Xemard, M.; Tafani, V.; Lagal, V.; Charras, G.; Tardieux, I. The toxoplasma-host cell junction is anchored to the cell cortex to sustain parasite invasive force. BMC Biol. 2014, 12, 773. [Google Scholar] [CrossRef] [Green Version]
  39. Denkers, E.Y.; Bzik, D.J.; Fox, B.A.; Butcher, B.A. An inside job: Hacking into Janus kinase/signal transducer and activator of transcription signaling cascades by the intracellular protozoan Toxoplasma gondii. Infect. Immun. 2012, 80, 476–482. [Google Scholar] [CrossRef] [Green Version]
  40. Nast, R.; Staab, J.; Meyer, T.; Luder, C. Toxoplasma gondii stabilises tetrameric complexes of tyrosine-phosphorylated signal transducer and activator of transcription-1 and leads to its sustained and promiscuous DNA binding. Cell Microbiol. 2018, 20, e12887. [Google Scholar] [CrossRef]
  41. Figueras, M.J.; Martin, O.A.; Echeverria, P.C.; de Miguel, N.; Naguleswaran, A.; Sullivan, W.J.; Corvi, M.M.; Angel, S.O. Toxoplasma gondii Sis1-like J-domain protein is a cytosolic chaperone associated to HSP90/HSP70 complex. Int. J. Biol. Macromol. 2012, 50, 725–733. [Google Scholar] [CrossRef] [Green Version]
  42. Weiss, L.M.; Ma, Y.F.; Takvorian, P.M.; Tanowitz, H.B.; Wittner, M. Bradyzoite development in Toxoplasma gondii and the hsp70 stress response. Infect. Immun. 1998, 66, 3295–3302. [Google Scholar] [CrossRef] [Green Version]
  43. Bell, A.; Monaghan, P.; Page, A.P. Peptidyl-prolyl cis-trans isomerases (immunophilins) and their roles in parasite biochemistry, host-parasite interaction and antiparasitic drug action. Int. J. Parasitol. 2006, 36, 261–276. [Google Scholar] [CrossRef] [Green Version]
  44. Silverman, J.A.; Hayes, M.L.; Luft, B.J.; Joiner, K.A. Characterization of anti-Toxoplasma activity of SDZ 215–918, a cyclosporin derivative lacking immunosuppressive and peptidyl-prolyl-isomerase-inhibiting activity: Possible role of a P glycoprotein in Toxoplasma physiology. Antimicrob. Agents Chemother. 1997, 41, 1859–1866. [Google Scholar] [CrossRef] [Green Version]
  45. Lu, J.; Zhou, D.H.; Chen, J.; Zhang, N.Z.; Wang, R.A.; Weng, Y.B.; Zhu, X.Q. Characterization of the Toxoplasma gondii hsp60 gene sequences from different hosts and geographical locations. Genet. Mol. Res. 2014, 13, 6906–6911. [Google Scholar] [CrossRef]
  46. Soto, A.S.; Fenoy, I.M.; Sanchez, V.R.; March, F.; Sibilia, M.D.P.; Aldirico, M.D.L.A.; Picchio, M.S.; Arcon, N.; Acosta, P.; Polack, F.P.; et al. Toxoplasma gondii serine-protease inhibitor-1: A new adjuvant candidate for asthma therapy. PLoS ONE 2017, 12, e0187002. [Google Scholar] [CrossRef] [Green Version]
  47. Pszenny, V.; Davis, P.H.; Zhou, X.W.; Hunter, C.A.; Carruthers, V.B.; Roos, D.S. Targeted disruption of Toxoplasma gondii serine protease inhibitor 1 increases bradyzoite cyst formation in vitro and parasite tissue burden in mice. Infect. Immun. 2012, 80, 1156–1165. [Google Scholar] [CrossRef] [Green Version]
  48. Dubois, D.; Fernandes, S.; Amiar, S.; Dass, S.; Katris, N.J.; Botte, C.Y.; Yamaryo-Botte, Y. Toxoplasma gondii acetyl-CoA synthetase is involved in fatty acid elongation (of long fatty acid chains) during tachyzoite life stages. J. Lipid Res. 2018, 59, 994–1004. [Google Scholar] [CrossRef] [Green Version]
  49. Ramakrishnan, S.; Docampo, M.D.; MacRae, J.I.; Ralton, J.E.; Rupasinghe, T.; McConville, M.J.; Striepen, B. The intracellular parasite Toxoplasma gondii depends on the synthesis of long-chain and very long-chain unsaturated fatty acids not supplied by the host cell. Mol. Microbiol. 2015, 97, 64–76. [Google Scholar] [CrossRef] [Green Version]
  50. Filipowicz, A.; Coca, M.N.; Blair, B.M.; Chang, P.Y. Acute myocarditis with cardiogenic shock and multiple organ failure, followed by bilateral panuveitis masquerading as endogenous endophthalmitis, due to Toxoplasma gondii in an immunocompetent patient. Retin. Cases Brief Rep. 2021, 15, 575–580. [Google Scholar] [CrossRef]
  51. Sauvant, G.; Bossart, W.; Kurrer, M.O.; Follath, F. Diagnosis and course of myocarditis: A survey in the medical clinics of Zurich university hospital 1980 to 1998. Schweiz. Med. Wochenschr. 2000, 130, 1265–1271. [Google Scholar]
  52. Barragan, A.; Brossier, F.; Sibley, L.D. Transepithelial migration of Toxoplasma gondii involves an interaction of intercellular adhesion molecule 1 (ICAM-1) with the parasite adhesin MIC2. Cell. Microbiol. 2005, 7, 561–568. [Google Scholar] [CrossRef]
  53. Bushkin, G.G.; Motari, E.; Carpentieri, A.; Dubey, J.P.; Costello, C.E.; Robbins, P.W.; Samuelson, J. Evidence for a structural role for acid-fast lipids in oocyst walls of Cryptosporidium, Toxoplasma, and Eimeria. mBio 2013, 4, e00387-13. [Google Scholar] [CrossRef] [Green Version]
  54. Nguyen, H.M.; Berry, L.; Sullivan, W.J.; Besteiro, S. Autophagy participates in the unfolded protein response in Toxoplasma gondii. FEMS Microbiol. Lett. 2017, 364, fnx153. [Google Scholar] [CrossRef] [Green Version]
  55. Joyce, B.R.; Tampaki, Z.; Kim, K.; Wek, R.C.; Sullivan, W.J. The unfolded protein response in the protozoan parasite Toxoplasma gondii features translational and transcriptional control. Eukaryot. Cell 2013, 12, 979–989. [Google Scholar] [CrossRef] [Green Version]
  56. Amara, N.; Foe, I.T.; Onguka, O.; Garland, M.; Bogyo, M. Synthetic fluorogenic peptides reveal dynamic substrate specificity of depalmitoylases. Cell Chem. Biol. 2019, 26, 35–47. [Google Scholar] [CrossRef]
Animals 12 02858 g001 550
Figure 1. Differential quantitative analysis of the modified peptides/proteins. (A) The mass distribution of the modified proteins. The x-axis represents the protein amount interval (kilodalton), and the y-axis is the number of modified proteins. (B) The unique peptide distribution of the target-modified proteins. The x-axis represents the number of unique matching peptides per protein, and the y-axis represents the number of modified proteins. (C) The motif distribution of the posttranslational modification sites of all the modified peptides. The x-axis represents the base number, and the y-axis represents the corrected score. The higher the base height, the higher the probability of the base appearing in the motif.
Figure 1. Differential quantitative analysis of the modified peptides/proteins. (A) The mass distribution of the modified proteins. The x-axis represents the protein amount interval (kilodalton), and the y-axis is the number of modified proteins. (B) The unique peptide distribution of the target-modified proteins. The x-axis represents the number of unique matching peptides per protein, and the y-axis represents the number of modified proteins. (C) The motif distribution of the posttranslational modification sites of all the modified peptides. The x-axis represents the base number, and the y-axis represents the corrected score. The higher the base height, the higher the probability of the base appearing in the motif.
Animals 12 02858 g001
Animals 12 02858 g002 550
Figure 2. Volcano of the up/downregulated differentially modified peptides. Red indicates upregulation, blue indicates downregulation, and gray indicates no significant difference. The first two significantly enriched upregulated and downregulated modified peptides are marked with horizontal lines. The Toxoplasma gondii-modified proteins A0A086JND6_TOXGO (putative translation elongation factor 2 family) and A0A086L1S1_TOXGO (condensin complex subunit 1) were not significantly different.
Figure 2. Volcano of the up/downregulated differentially modified peptides. Red indicates upregulation, blue indicates downregulation, and gray indicates no significant difference. The first two significantly enriched upregulated and downregulated modified peptides are marked with horizontal lines. The Toxoplasma gondii-modified proteins A0A086JND6_TOXGO (putative translation elongation factor 2 family) and A0A086L1S1_TOXGO (condensin complex subunit 1) were not significantly different.
Animals 12 02858 g002
Animals 12 02858 g003 550
Figure 3. GO term distribution of the proteins sequence corresponding to the modified peptides in the small intestine of cats. Red indicates upregulation, and green indicates downregulation. The bar graph shows the number of modified peptides enriched in GO terms belonging to the three GO categories, biological process (BP), cellular component (CC), and molecular function (MF), at 10 DPI. The x-axis represents the GO terms, and the y-axis represents the number of upregulated and downregulated modified peptides in the different GO terms. The five-pointed star (★) represents that some Toxoplasma-related modified peptides were also involved in the host’s GO pathway. One five-pointed star (★) indicates that less than 4 were involved, and two five-pointed stars (★★) indicate that more than 7 peptides were involved.
Figure 3. GO term distribution of the proteins sequence corresponding to the modified peptides in the small intestine of cats. Red indicates upregulation, and green indicates downregulation. The bar graph shows the number of modified peptides enriched in GO terms belonging to the three GO categories, biological process (BP), cellular component (CC), and molecular function (MF), at 10 DPI. The x-axis represents the GO terms, and the y-axis represents the number of upregulated and downregulated modified peptides in the different GO terms. The five-pointed star (★) represents that some Toxoplasma-related modified peptides were also involved in the host’s GO pathway. One five-pointed star (★) indicates that less than 4 were involved, and two five-pointed stars (★★) indicate that more than 7 peptides were involved.
Animals 12 02858 g003
Animals 12 02858 g004 550
Figure 4. COG term distribution of the modified peptides for feline small intestine. The x-axis is the number of proteins corresponding to functional categories, and the y-axis is the COG category entries. The graph represents the statistical number of proteins with different functions in the sample. In the COG classification of the 24 host-related modified peptides, 8 (O, G, S, T, P, P, J, L, C)-related modified peptides of T. gondii were also involved.
Figure 4. COG term distribution of the modified peptides for feline small intestine. The x-axis is the number of proteins corresponding to functional categories, and the y-axis is the COG category entries. The graph represents the statistical number of proteins with different functions in the sample. In the COG classification of the 24 host-related modified peptides, 8 (O, G, S, T, P, P, J, L, C)-related modified peptides of T. gondii were also involved.
Animals 12 02858 g004
Animals 12 02858 g005 550
Figure 5. KEGG pathway annotation for the small intestine of cats. The pathways that were significantly enriched for proteins corresponding to the differentially modified peptides. The x-axis represents the number of differentially modified peptides in the corresponding KEGG pathways in each KEGG subsystem. The y-axis represents the main clusters of the KEGG pathways.
Figure 5. KEGG pathway annotation for the small intestine of cats. The pathways that were significantly enriched for proteins corresponding to the differentially modified peptides. The x-axis represents the number of differentially modified peptides in the corresponding KEGG pathways in each KEGG subsystem. The y-axis represents the main clusters of the KEGG pathways.
Animals 12 02858 g005
Animals 12 02858 g006 550
Figure 6. KEGG pathway enrichment of the host. The top 30 enriched KEGG pathways of the identified modified peptides in this study. The y-axis represents the distinct KEGG pathways, and the x-axis represents the Rich factor. The Rich factor refers to the ratio of modified peptides annotated in the pathway to the total number of modified peptides annotated in the pathway. The greater the Rich factor, the greater the degree of pathway enrichment. The dot size represents the number of modified peptides (bigger dots denote a large number of modified peptides and vice versa). The colors of the dots represent the p-values of enrichment: red indicates high enrichment, blue indicates low enrichment.
Figure 6. KEGG pathway enrichment of the host. The top 30 enriched KEGG pathways of the identified modified peptides in this study. The y-axis represents the distinct KEGG pathways, and the x-axis represents the Rich factor. The Rich factor refers to the ratio of modified peptides annotated in the pathway to the total number of modified peptides annotated in the pathway. The greater the Rich factor, the greater the degree of pathway enrichment. The dot size represents the number of modified peptides (bigger dots denote a large number of modified peptides and vice versa). The colors of the dots represent the p-values of enrichment: red indicates high enrichment, blue indicates low enrichment.
Animals 12 02858 g006
Animals 12 02858 g007 550
Figure 7. Protein–protein interaction (PPI) networks for the differentially modified peptides/proteins of the host. Red indicates upregulation, and green indicates downregulation. The filled circle represents protein, and the direction of the arrow indicates regulation. The dotted circles represent three MCODE clusters (MCODE clusters 1: viral myocarditis, MCODE clusters 2: fatty acid elongation, and MCODE clusters 3: protein processing in endoplasmic reticulum).
Figure 7. Protein–protein interaction (PPI) networks for the differentially modified peptides/proteins of the host. Red indicates upregulation, and green indicates downregulation. The filled circle represents protein, and the direction of the arrow indicates regulation. The dotted circles represent three MCODE clusters (MCODE clusters 1: viral myocarditis, MCODE clusters 2: fatty acid elongation, and MCODE clusters 3: protein processing in endoplasmic reticulum).
Animals 12 02858 g007
Table
Table 1. Results of the N-glycosylation identification by Toxoplasma gondii.
Table 1. Results of the N-glycosylation identification by Toxoplasma gondii.
SampleTotal_SpectraIdentified_SpectraProteinsGlycoproteins aPeptidesGlycopeptides bGlycosites
10_DPI_1287,4611510209429044
10_DPI_2288,6861622171425954
10_DPI_3281,3731491165324043
a The protein A0A086JND6_TOXGO was identified in 10_DPI_1, 10_DPI_2, and 10_DPI_3, and the protein A0A086K8H6_TOXGO was identified in 10_DPI_2 and 10_DPI_3. b The peptide_seq NMSVIAHVDHGK was identified in 10_DPI_1 (once), 10_DPI_2 (twice), and 10_DPI_3 (twice), and the peptide_seq KKNFSDSGNFGFGIQEHIDLGIK was identified once in both 10_DPI_2 and 10_DPI_3.
Table
Table 2. The KEGG enrichment results of the T. gondii-related protein.
Table 2. The KEGG enrichment results of the T. gondii-related protein.
No.PathwayDifferent Proteins with Pathway AnnotationProteinsPathway ID
1Ribosome1 (50%)tr|A0A086K8H6|A0A086K8H6_TOXGOko03010
2AMPK signaling pathway1 (50%)tr|A0A086JND6|A0A086JND6_TOXGOko04152
3Oxytocin signaling pathway1 (50%)tr|A0A086JND6|A0A086JND6_TOXGOko04921
4Amino sugar and nucleotide sugar metabolism1 (25%)tr|A0A151HN22|A0A151HN22_TOXGOko00520
5Aminoacyl-tRNA biosynthesis1 (25%)tr|A0A125YXZ8|A0A125YXZ8_TOXGVko00970
6Ribosome biogenesis in eukaryotes1 (25%)tr|A0A086KYF0|A0A086KYF0_TOXGOko03008
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhai, B.; Xie, S.; Peng, J.; Qiu, Y.; Liu, Y.; Zhu, X.; He, J.; Zhang, J. Glycosylation Analysis of Feline Small Intestine Following Toxoplasma gondii Infection. Animals 2022, 12, 2858. https://doi.org/10.3390/ani12202858

AMA Style

Zhai B, Xie S, Peng J, Qiu Y, Liu Y, Zhu X, He J, Zhang J. Glycosylation Analysis of Feline Small Intestine Following Toxoplasma gondii Infection. Animals. 2022; 12(20):2858. https://doi.org/10.3390/ani12202858

Chicago/Turabian Style

Zhai, Bintao, Shichen Xie, Junjie Peng, Yanhua Qiu, Yang Liu, Xingquan Zhu, Junjun He, and Jiyu Zhang. 2022. "Glycosylation Analysis of Feline Small Intestine Following Toxoplasma gondii Infection" Animals 12, no. 20: 2858. https://doi.org/10.3390/ani12202858

APA Style

Zhai, B., Xie, S., Peng, J., Qiu, Y., Liu, Y., Zhu, X., He, J., & Zhang, J. (2022). Glycosylation Analysis of Feline Small Intestine Following Toxoplasma gondii Infection. Animals, 12(20), 2858. https://doi.org/10.3390/ani12202858

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop