Issue |
Parasite
Volume 31, 2024
|
|
---|---|---|
Article Number | 23 | |
Number of page(s) | 10 | |
DOI | https://doi.org/10.1051/parasite/2024027 | |
Published online | 17 May 2024 |
Research Article
Quantitative phosphoproteomic analysis of chicken DF-1 cells infected with Eimeria tenella, using tandem mass tag (TMT) and parallel reaction monitoring (PRM) mass spectrometry
Analyse phosphoprotéomique quantitative de cellules DF-1 de poulet infectées par Eimeria tenella, par spectrométrie de masse avec marqueur de masse en tandem (TMT) et surveillance des réactions parallèles (PRM)
Shanghai Veterinary Research Institute, Chinese Academy of Agricultural Sciences, Key Laboratory of Animal Parasitology of Ministry of Agriculture, Minhang, Shanghai 200241, PR China
* Corresponding author: donghui@shvri.ac.cn
Received:
13
November
2023
Accepted:
1
May
2024
Eimeria tenella is an obligate intracellular parasite which causes great harm to the poultry breeding industry. Protein phosphorylation plays a vital role in host cell–E. tenella interactions. However, no comprehensive phosphoproteomic analyses of host cells at various phases of E. tenella infection have been published. In this study, quantitative phosphoproteomic analysis of chicken embryo DF-1 fibroblasts that were uninfected (UI) or infected with E. tenella for 6 h (PI6, the early invasion phase) or 36 h (PI36, the trophozoite development phase) was conducted. A total of 10,122 phosphopeptides matched to 3,398 host cell phosphoproteins were identified and 13,437 phosphorylation sites were identified. Of these, 491, 1,253, and 275 differentially expressed phosphorylated proteins were identified in the PI6/UI, PI36/UI, and PI36/PI6 comparisons, respectively. KEGG pathway enrichment analysis showed that E. tenella modulated host cell processes through phosphorylation, including focal adhesion, regulation of the actin cytoskeleton, and FoxO signaling to support its early invasion phase, and modulating adherens junctions and the ErbB signaling pathway to favor its trophozoite development. These results enrich the data on the interaction between E. tenella and host cells and facilitate a better understanding of the molecular mechanisms underlying host–parasite relationships.
Résumé
Eimeria tenella est un parasite intracellulaire obligatoire qui cause de graves dommages à l’industrie de l’élevage de volailles. La phosphorylation des protéines joue un rôle essentiel dans les interactions entre la cellule hôte et E. tenella. Cependant, aucune analyse phosphoprotéomique complète des cellules hôtes à différentes phases de l’infection par E. tenella n’a été publiée. Dans cette étude, une analyse phosphoprotéomique quantitative de fibroblastes DF-1 d’embryon de poulet non infectés (NI) ou infectés par E. tenella pendant 6 h (PI6, la phase d’invasion précoce) ou 36 h (PI36, la phase de développement des trophozoïtes) a été réalisée. Un total de 10 122 phosphopeptides correspondant à 3 398 phosphoprotéines de cellules hôtes ont été identifiés et 13 437 sites de phosphorylation ont été identifiés. Parmi celles-ci, 491, 1 253 et 275 protéines différentiellement phosphorylées exprimées ont été identifiées respectivement dans les comparaisons PI6/NI, PI36/NI et PI36/PI6. L’analyse d’enrichissement de la voie KEGG a montré qu’E. tenella modulait les processus de la cellule hôte par phosphorylation, y compris l’adhésion focale, la régulation du cytosquelette d’actine et la signalisation FoxO, pour aider sa phase d’invasion précoce, et la modulation des jonctions adhérentes et de la voie de signalisation ErbB pour favoriser le développement de son trophozoïte. Ces résultats enrichissent les données sur l’interaction entre E. tenella et les cellules hôtes et facilitent une meilleure compréhension des mécanismes moléculaires sous-jacents aux relations hôtes–parasites.
Key words: Coccidiosis / Eimeira tenella / DF-1 cells / TMT / Phosphoproteomics
© L.-S. Jia et al., published by EDP Sciences, 2024
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Introduction
Chicken coccidiosis is caused by the genus Eimeria and is one of the most common and devastating parasitic diseases in the poultry industry [3]. It can damage the intestinal epithelial cells of chickens, hinder the absorption of nutrients, and reduce the production of meat and eggs, causing economic losses of approximately £10.4 billion annually worldwide [4]. Currently, the prevention and control methods of coccidiosis mainly include anticoccidial drugs and live vaccines. However, these methods raise several issues ranging from drug resistance to food security and production costs [6, 39]. Therefore, it is necessary to find novel measures to control coccidiosis.
Eimeria relies on host cells to complete its life cycle. In the process of Eimeria infection, some physiological and biochemical processes of host cells are regulated, supporting the invasion and development of the parasite. For instance, invasion of host cells by Eimeria tenella causes host F-actin aggregation on the sporozoite surface. This aggregation is closely related to the invasion efficiency of E. tenella, and inhibition of F-actin aggregation by cytochalasin D reduces the invasion rate of E. tenella [47]. The expression of host fatty acid binding protein 4 increases significantly at 72 h after invasion of E. tenella sporozoites, and overexpression of this protein significantly inhibits the invasion of E. tenella sporozoites [45]. Incubation of host cells with polyclonal antibodies targeting receptors for activated C kinase 1 (RACK1) can promote the invasion of E. tenella sporozoites [51]. The vimentin transcript and protein levels increase continually at 6–72 h after E. tenella infection, and treating DF-1 cells with polyclonal anti-vimentin or knocking down vimentin by small interfering RNA significantly improves the invasion efficiency of sporozoites, indicating that vimentin plays an inhibitory role during the invasion of sporozoites [31]. These results demonstrate that host proteins play an indispensable role in Eimeria invasion and intracellular development. However, the extent and relevance of Eimeria-induced host responses remain poorly studied.
Many studies have confirmed that protein phosphorylation is involved in the interaction between pathogenic microorganisms and host cells [11, 18]. Phosphorylation is one of the most common post-translational modifications in organisms and plays an important role in cell growth, differentiation, apoptosis, and cell signaling in healthy conditions [41]. It occurs mainly on serine, threonine, or tyrosine residues and can affect the function of one third of the proteins [9]. Studies have found that the PI3K–PKB/Akt pathway may be one of the major routes through which Toxoplasma gondii prevents host cell apoptosis and that T. gondii phosphorylates the pro-apoptotic protein bad to prevent apoptosis [36]. So far, no studies on host cell phosphoproteomics after E. tenella infection have been published.
In recent years, tandem mass tag (TMT)-labeled quantitative proteomics technology coupled with high-resolution mass spectrometry (MS) has been proven to be an efficient technique for proteome studies simultaneously identifying and comparing the relative contents of proteins in up to 10 different samples. It has been widely used in the study of phosphoproteomics [29]. In the present study, the phosphoproteomes of DF-1 cells before and after E. tenella infection were comparatively analyzed by TMT-labeled quantitative proteomics. The data in this study are helpful to identify key host proteins associated with E. tenella invasion and provide an important basis to further analyze E. tenella–host cell interactions.
Materials and methods
Parasite, cell culture and ethics statement
The E. tenella Shanghai strain used in the current study has been maintained in our laboratory since 1993. E. tenella was propagated by passage through 14-day-old coccidia-free chickens [43]. Sporozoites were collected and purified from cleaned sporulated oocysts with standard procedures [17, 50]. DF-1 cells were cultured in Dulbecco’s modified Eagle’s medium (DMEM) (Life Technology, Carlsbad, CA, USA) supplemented with 10% fetal bovine serum at 37 °C in a 5% CO2 incubator.
The animal protocols were approved by the Institutional Animal Care and Use Committee of Shanghai Veterinary Research Institute, Chinese Academy of Agricultural Sciences (Permit Number: SHVRI-SZ-20230323-4).
Parasite infection and sample collection
DF-1 cells were seeded into T25 culture flasks (1 × 106 cells per flask) and allowed to reach 80–90% confluence before infection. Freshly excysted sporozoites were incubated with DMEM supplemented with 2% fetal bovine serum and 5% penicillin/streptomycin at 37 °C for 2 h. DF-1 cells were infected with pre-treated sporozoites at a multiplicity of infection of 3 for 6 h (6 h post-infection, PI6) or 36 h (36 h post-infection, PI36). Uninfected (UI) cells were used as control. Each experiment was performed with three independent biological replicates.
At 6 h after infection, the cell culture supernatant containing the sporozoites was removed and the non-invasive sporozoites were removed by washing gently with PBS three times. Cells in the UI and PI6 groups were harvested with a cell scraper; cells in the PI36 group were cultured with complete DMEM for another 30 h and then harvested with a cell scraper. All cells were pelleted by centrifugation and stored at −80 °C for subsequent analysis.
Protein extraction and peptide enzymatic hydrolysis
The cell samples were lysed with SDT buffer (4% SDS, 100 mM Tris/HCl pH 7.6, and 0.1 M DTT), incubated for 3 min in boiling water, sonicated on ice twice, and clarified by centrifugation at 16,000 ×g at 25 °C for 10 min. Protein content was quantified with a BCA Protein Assay Kit (Bio-Rad, Hercules, CA, USA).
A moderate amount of protein was extracted from each sample, and trypsin enzymatic hydrolysis was performed using the filter-aided sample preparation method [49], followed by desalting enzymolysis of peptides using a C18 Cartridge. Lyophilized peptides were dissolved with 40 μL dissolution buffer and the peptides were quantified by measuring the optical density at 280 nm.
TMT labeling and SCX chromatographic classification
First, 100 μg peptide mixture was taken from each sample and labeled using a TMT labeling kit, according to the manufacturer’s instructions (Thermo Fisher, Waltham, MA, USA). The labeled peptides were fractionated by strong cation exchange (SCX) liquid chromatography (LC) using an AKTA Purifier 100 system. Gradient elution was performed using buffer A (10 mM KH2PO4, 25% acetonitrile, pH 3.0) and buffer B (10 mM KH2PO4, 500 mM KCl, 25% acetonitrile, pH 3.0). The chromatographic column was balanced using buffer A and the sample was applied to the column for separation. The flow rate was 1 mL/min and the gradient elution program was as follows: 0–25 min, 0–10% buffer B; 25–32 min, 10–20% buffer B; 32–42 min, 20–45% buffer B; 42–47 min, 45–100% buffer B; 47–60 min, 100% buffer B; 60 min, 0% buffer B. The absorbance at 214 nm was monitored during elution, fractions were collected every 1 min, and the collected fractions were desalted with a C18 Cartridge [23].
TiO2 enrichment of phosphopeptides
The labeled peptide mixtures were vacuum-freeze dried and dissolved with 1× DHB buffer [obtained by diluting 5× DHB buffer (3% dihydrobenzoic acid, 80% acetonitrile, 0.1% trifluoroacetic acid) with water at a 1:4 ratio]. Next, TiO2 beads were added to the solution, samples were oscillated for 40 min and centrifuged, and the supernatant was removed. The precipitates were washed three times with wash buffer I (30% acetonitrile and 3% trifluoroacetic acid) and three times with washing buffer II (80% acetonitrile and 0.3% trifluoroacetic acid). The enriched phosphorylated peptides were eluted with elution buffer (40% acetonitrile and 15% NH4OH) and then vacuum concentrated. The phosphorylated peptides were redissolved in 10 μL of 0.1% formic acid, and 5 μL was taken for MS analysis.
LC-MS/MS analysis
Each sample was analyzed by a nano-scale high-performance liquid phase system (EASY-nLC) using buffer A (0.1% aqueous formic acid solution) and buffer B (0.1% aqueous acetonitrile solution). The chromatographic column was balanced using 95% of buffer A, and then sample was added to the sample column (Thermo Scientific EASY column, 100 μm × 2 cm, 5 μm, C18) and separated by analytical columns (Thermo Scientific EASY column, 75 μm × 10 cm, 3 μm, C18). The flow rate was 250 nL/min, and the gradient program was as follows: 0–50 min, 0–35% buffer B; 50–58 min, 35–100% buffer B; 58–60 min, 100% buffer B.
After chromatographic separation, the samples were analyzed by LC-MS/MS using a Q-Exactive HF-X mass spectrometer. The analysis time was 120 min, detection was performed in positive ion mode, and the scanning range of parent ions was 300–1,800 m/z. The resolution of the first-order mass spectrum was 70,000 at m/z 200, the automatic control target was set at 3e6, and the maximum injection time was set to 10 ms. The dynamic exclusion time was 40.0 s. Ten MS2 scan debris maps were collected after each full scan. The activation type of MS2 was high collision dissociation, and the isolation window was 2 m/z. The resolution of secondary MS was 17,500 at m/z 200. The normalized collision energy was 30 eV and the underfill ratio was set to 0.1%.
Phosphoprotein identification and quantitative analysis
Raw data were processed by Proteome Discoverer (Version 1.4, Thermo Electron, San Jose, CA, USA) and matched against the UniProt Gallus gallus database (20230902, 35,124 sequences) by scoring through the Mascot server (Version 2.3, Matrix Science, London, UK). Data were searched using the following parameters: peptide mass tolerance = 20 ppm; enzyme = trypsin; fragment mass tolerance = 0.1 Da; fixed modification = TMT 10plex (N-term), carbamidomethyl (C), TMT 10plex (K); variable modification = oxidation (M), TMT 10plex (Y), deamidated (NQ), phosphorylation (S/T/Y); max missed cleavage = 2. The incorporated Target Decoy PSM Validator in Proteome Discoverer was used to validate the search results with only the hits with a false discovery rate of ≤0.01. The protein ratios were calculated as the median of only unique peptides of the protein. All peptide ratios were normalized by the median protein ratio. The median protein ratio should be 1 after normalization. A phosphoprotein showing a median fold change of >1.5 or <0.67 in phosphorylation level within the group at p ≤ 0.05 was considered to be a significantly differentially regulated phosphoprotein.
Bioinformatics analysis
All significantly differentially regulated phosphoproteins in the three comparisons (PI6/UI, PI36/UI, and PI36/PI6) were subjected to bioinformatics analysis. Blast2GO (https://www.blast2go.com/) was used to for Gene Ontology (GO) annotation and enrichment analysis in the cell component, biological process, and molecular function categories. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis was performed using the online KEGG Automatic Annotation Server, and the annotation results were mapped to the KEGG pathway database using KEGG mapper. A two-tailed Fisher exact test was employed to assess the enrichment of proteins with differentially regulated phosphosites (GO, KEGG pathway, and protein domain; p < 0.05 was considered significant).
PRM validation of proteins with differentially regulated phosphosites
Validation of proteomics data by parallel reaction monitoring (PRM) is the latest development in targeted MS, serving as the “golden standard” for protein MS analysis. The differentially expressed phosphorylated proteins (DEPPs) from the results of the TMT analysis were further validated by using PRM analysis. Briefly, the protein samples were prepared by the same method as used for LC-MS/MS analysis. PRM analyses were performed on a Q-Exactive mass spectrometer (Thermo Fisher Scientific). The reverse-phase liquid chromatography (RPLC) fractions of each sample were analyzed using an acquisition method. For each sample, 2 μg of peptides were taken and mixed with 20 fMol of stable isotope internal standard peptide. A total of nine samples were obtained. PRM data analysis was performed using Skyline software (MacCoss Lab, University of Washington, Seattle, WA, USA).
Results
Identification and quantification of phosphopeptides
The phosphorylated peptides of each group were enriched with TiO2 and identified by LC-MS/MS. Through the statistical analysis of raw data, a total of 39,960 peptide spectrum matches (PSMs) were identified and matched to 3,516 host cell proteins (Supplementary Table S1). Among them, 10,122 phosphopeptides (Supplementary Table S2) corresponded to 3,398 host cell phosphoproteins, and 13,437 phosphorylation sites were detected. Among the phosphorylation sites, phospho-serine (pS) was the most abundant, accounting for 84.81% of all the phosphorylated amino acids, followed by phospho-threonine (pT, 14.46%) and phospho-tyrosine (pY, 0.73%) (Figure 1).
Figure 1 Proportion of serine, threonine, and tyrosine in phosphorylation sites. |
Identification and quantification of the significantly differentially regulated phosphoproteins
To speculate on host proteins that play a role during parasite infection, we analyzed the DEPPs at different time points of infection. The significantly differentially expressed phosphorylated peptides were identified based on a fold change value of >1.5 or <0.67 and a p-value of <0.05. A comparative phosphorylation level analysis was performed among the UI, PI6, and PI36 groups.
The comparative analysis led to the identification of 489 upregulated and 2 downregulated phosphoproteins in the PI6/UI comparison, 1,228 upregulated and 25 downregulated phosphoproteins in the PI36/UI comparison, and 141 upregulated and 134 downregulated phosphoproteins in the PI36/PI6 comparison (Table 1, Figure 2, Supplementary Table S3). The differentially expressed peptides in the clustering heatmap (Figure 3) show that the protein level trends are consistent between biological replicates in the infected and uninfected groups.
Figure 2 Volcano plots from different group comparisons. The abscissa indicates difference multiple (logarithmic transformation based on 2), the ordinate indicates the significant of difference (logarithmic transformation based on 10). The red point is significantly upregulated phosphorylated peptide segment, the blue point is significantly downregulated phosphorylated peptide segment and the gray point is a phosphorylated peptide segment with no significant difference. |
Figure 3 Clustering heatmap of different expression phosphorylated peptides. Each row represents a phosphorylated peptide segment, each column represents a group of samples. The logarithmic value (logarithmic transformation based on 2) of the significantly differentially expressed phosphorylated peptides in different samples is displayed in the clustering heatmap in different colors. Red represents significant up-regulation of phosphorylated peptides; blue represents significant down-regulation of phosphorylated peptides. |
The number of significantly regulated phosphopeptides/proteins among the three comparison groups.
GO analysis of differentially expressed phosphoproteins
Fisher’s exact test was used for GO functional enrichment analysis of the DEPPs corresponding to differentially expressed phosphorylated peptides in the three comparisons (PI6/UI, PI36/UI, and PI36/PI6) to identify the top five enriched terms for each comparison.
Interestingly, in the cellular component category, the main DEPPs of the three comparisons (PI6/UI, PI36/UI, and PI36/PI6) were all concentrated in the intracellular part. In addition, in the cellular component category, the GO terms were similar between the PI6/UI and PI36/PI6 comparisons, including intracellular, organelle, and intracellular organelle part. In the PI36/UI comparison, the GO terms were mostly related to the nucleus, including nuclear part and nuclear lumen (Figure 4A).
Figure 4 Gene ontology annotations of the differentially expressed phosphorylated proteins. The abscissa indicates the enriched GO functional classification, including biological process (A), cellular component (B), and molecular function (C). The ordinate indicates the size of the significance of corresponding to each entry, the more left, the more significant. |
In the biological process category, the DEPPs of PI6/UI and PI36/PI6 were found to be mainly involved in cellular component organization, cellular component organization or biogenesis, and organelle organization. However, it is worth mentioning that we found that metabolic-related biological processes were enriched in the PI36/UI comparison, including regulation of metabolic process, regulation of macromolecule metabolic process, and regulation of cellular metabolic process. Meanwhile, some DEPPs in the PI6/UI comparison were also related to regulation of metabolic process (Figure 4B). These results suggest that sporozoite invasion and intracellular development may depend on the energy provided by host cells.
Additionally, the GO terms in the molecular function category of all three comparisons were associated with binding, including poly(A) RNA binding, RNA binding, binding, and protein binding (Figure 4C). This suggests that both transcription and translation levels were quite active in cells at the parasite infection stage.
KEGG pathway analysis of differentially expressed phosphoproteins
The functionally annotated DEPPs were classified using KEGG pathway enrichment analysis (Figure 5). In both the PI6/UI and PI36/UI comparisons, the top three enriched signaling pathways were focal adhesion, regulation of actin cytoskeleton, and FoxO signaling pathways, while in the PI36/PI6 comparison, adherens junctions, the ErbB signaling pathway, and herpes simplex infection were enriched.
Figure 5 KEGG pathways of the differentially expressed phosphorylated proteins. The abscissa indicates the first 10 significantly enriched KEGG pathways and the ordinate indicates the significance of enriched KEGG pathways, the more left, the more significant. |
Moreover, some signaling pathways associated with cytoskeleton regulation were enriched in all three comparisons, including regulation of actin cytoskeleton and focal adhesion, indicating that cytoskeletal changes run through the whole process of E. tenella sporozoite infection. The mitogen-activated protein kinase (MAPK) and ErbB signaling pathways were also enriched in the PI36/UI comparison, which play important roles in cell signal transduction, transcription, inflammation, and apoptosis.
PRM results
We selected eight DEPPs for PRM validation from the TMT screening results. Skyline was used to analyze the target peptides, which included the information of chromatographic peaks, original peak area, and the contrast histogram of the original peak area. The results showed that the changes of phosphorylation levels of Zfp36l1, TGFB1I1, MKI67, GIGYF2, and ANLN in all comparisons were consistent with the TMT results, while the changes of phosphorylation levels of VIM, PDCD4, and MAP1B in some comparisons were slightly different from the TMT results. This is because TMT is large-scale screening of non-target proteins; there is a certain false positive rate, so there may be some differences with the PRM results. In conclusion, our TMT results are reliable and reproducible (Table 2).
PRM analysis and comparison of the quantitative results for candidate phosphoproteins.
Discussion
Several studies have shown that a large number of host phosphoproteins are involved in the interaction between parasites and host cells. The analysis and identification of key molecules related to the invasion of host cells by E. tenella is of positive significance for elucidating the molecular mechanism of invasion and infection of E. tenella and effectively controlling coccidiosis. In this study, we screened DEPPs of DF-1 cells infected with E. tenella sporozoites at 6 h and 36 h by TMT-labeled quantitative proteomic technology combined with LC-MS/MS. In total, 10,122 phosphopeptides were identified, matching to 3,398 host phosphoproteins. These DEPPs are mainly involved in cytoskeleton regulation, apoptosis regulation, the immune response, and substance metabolism.
It has been found that many apicomplexans, such as T. gondii, Plasmodium, and Cryptosporidium, remodel the host cytoskeleton in the process of host cell invasion [13, 14, 19, 33, 41, 48]. In the present study, we found that a large number of DEPPs are involved in regulation of the actin cytoskeleton. For example, the phosphorylation level of Rho GTPase activating protein was upregulated at 6 h after infection, and actin-depolymerizing factor, focal adhesion kinase 1, and Gelsolin were only upregulated at 36 h after infection. However, the phosphorylation levels of vimentin, vinculin, paxillin, and LIM domain kinase 1 were upregulated at 6 h and 36 h after infection. Vimentin is an important component of the cytoskeleton; phosphorylation can induce disintegration and spatial orientation alterations of vimentin, thereby regulating cell contraction and local adhesion kinetics [40]. Vinculin can be activated by phosphorylation and plays a role in binding and rearrangement of the actin cytoskeleton [1]. Rho GTPase is a key regulator of cytoskeletal dynamics and affects many cellular processes [20]. Phosphorylation inhibits the binding of GTP to Rac1 and reduces its overall activity [28]. Therefore, we speculate that sporozoites may regulate the cytoskeleton in many ways to facilitate its invasion and intracellular development.
Previous studies have shown that parasites can regulate apoptosis of host cells and create favorable conditions for their parasitism in cells [7, 12, 21, 27]. We found that the phosphorylation level of apoptosis inhibitor 5 (API-5) after 36 h of sporozoite infection was significantly higher than after 6 h and in the uninfected group. Studies have shown that the API-5 phosphorylation level is negatively correlated with the apoptosis rate [37, 38]. Furthermore, the phosphorylation level of apoptosis-antagonizing transcription factor (AATF) after 36 h of sporozoite infection was significantly higher than that of the uninfected group. AATF is a key regulatory factor in the p53 response, and phosphorylation causes AATF to be released from the cytoplasm of MRLC3 and bind to the promoter regions of PUMA, Bax, and BAK, while inhibiting the expression of p53 driven pro apoptotic genes [22]. At 36 h after sporozoite infection, the parasite is at the trophozoite developmental stage and needs the host to provide a suitable development site. Therefore, we speculate that the sporozoite may inhibit apoptosis by regulating the activity of apoptosis-related proteins to facilitate its intracellular development.
Parasites can escape host immune attacks through a variety of immune escape mechanisms. Some studies have found that ATG5 expression was significantly downregulated after E. tenella sporozoites infect CEF cells and that ATG5 plays an important role in the formation of autophagy vesicles and has proapoptotic activity [44, 46]. In our study, the phosphorylation level of Atg13 after 36 h of sporozoite infection was 1.685 times that of the uninfected group and 1.309 times that at 6 h after infection. Phosphorylation of Atg13 can prevent its binding to Atg1, while tight Atg1–Atg13 binding is required for autophagic activity [26]. It is suggested that sporozoites may escape immune clearance by interfering with host autophagy to facilitate their invasion and intracellular development.
In addition, host cells also fight parasitic infections through a variety of immune regulatory mechanisms and immune-related proteins [30, 35]. In the present study, the phosphorylation levels of proteins associated with MAPK pathways, such as MAPK1, MAPK11, and MAPK14, were significantly increased. MAPK proteins are mainly transmitters of signals from the cell surface to the nucleus and have important regulatory roles in the production of inflammatory mediators [34]. At 6 h after infection, the expression of Akirin2, a nuclear protein that plays an important role in the innate immune response and embryonic development [15, 32, 42], was 1.6 times higher than in the uninfected group. The transcription factor AP-1 was upregulated at 6 h and 36 h after infection. It is suggested that host cells activate multiple transcription factors such as AP-1 to mediate inflammation initiated by E. tenella infection.
All obligate intracellular parasites face the great challenge of getting nutrients from host cells. E. falciformis [24] and T. gondii [2] can widely affect host cell metabolic pathways or change gene expression patterns. Plasmodium are completely dependent on host cell glucose to obtain energy and actively internalize phospholipids from the erythrocyte membrane and the extracellular medium [16, 25]. We found that at 6 h and 36 h after sporozoite infection, the phosphorylation level of 6-phosphofructokinase was 1.4 and 1.5 times higher than in the uninfected group. This enzyme catalyzes the conversion of fructose 6-phosphate to fructose 1,6-diphosphate, which is one of the rate-limiting steps of glycolysis. Phosphorylation changes the tetramer conformation of 6-phosphofructokinase, making its structure stable and preventing its activity from being inhibited by ATP [5]. Furthermore, the increased affinity of phosphorylated 6-phosphofructokinase to F-actin prevents the inhibitory effect of lactic acid on fructose 6-phosphofructokinase and improves glycolysis efficiency [8, 10]. The changes in phosphorylation levels of metabolic-related enzymes suggest that sporozoites may ensure adequate nutritional supply during parasitism by regulating the activity of metabolic pathway-related enzymes.
Conclusion
In this study, we analyzed the DEPPs between DF-1 cells infected with E. tenella sporozoites and uninfected cells, identifying a total of 3,398 differentially expressed phosphoproteins. Our results suggest that the host cytoskeleton may be remodeled during E. tenella sporozoite infection, while parasites also evade immune clearance by regulating host apoptosis and interfering with host autophagy. This study provides relevant data and new directions for further exploring the interactions between parasites and host cells.
Acknowledgments
This research was supported by the National Natural Science Foundation of China (No. 32373038) and the National Parasitic Resources Center (NPRC-2019-194-30).
Conflicts of Interest
The authors declare that there is no conflict of interest.
Author contributions statement
Liu-Shu Jia and Zhan Liu designed the experiment, prepared test samples, performed the experiments, and wrote the original draft. Shun-Hai Zhu, Qi-Ping Zhao, and Hong-Yu Han designed animal experiments and participated in experiments. Huan-Zhi Zhao and Yu Yu assisted with the experiment and analyzed data. Hui Dong designed the experiments and revised the manuscript. All authors have read and agreed to the published version of the manuscript.
Supplementary Material
Table S1. Details list of proteins of DF-1 cells infected with Eimeria tenella.
Table S2. Details list of phosphorylation peptides of DF-1 cells infected with Eimeria tenella.
Table S3. Details list of significantly regulated phosphopeptides of DF-1 cells infected with Eimeria tenella.
Access hereReferences
- Bays JL, DeMali KA. 2017. Vinculin in cell-cell and cell-matrix adhesions. Cellular and Molecular Life Sciences, 74(16), 2999–3009. [CrossRef] [PubMed] [Google Scholar]
- Blader IJ, Manger ID, Boothroyd JC. 2001. Microarray analysis reveals previously unknown changes in Toxoplasma gondii-infected human cells. Journal of Biological Chemistry, 276(26), 24223–24231. [CrossRef] [Google Scholar]
- Blake DP, Knox J, Dehaeck B, Huntington B, Rathinam T, Ravipati V, Ayoade S, Gilbert W, Adebambo AO, Jatau ID, Raman M, Parker D, Rushton J, Tomley FM. 2020. Re-calculating the cost of coccidiosis in chickens. Veterinary Research, 51(1), 115. [CrossRef] [PubMed] [Google Scholar]
- Blake DP, Tomley FM. 2014. Securing poultry production from the ever-present Eimeria challenge. Trends in Parasitology, 30(1), 12–19. [CrossRef] [PubMed] [Google Scholar]
- Cai GZ, Callaci TP, Luther MA, Lee JC. 1997. Regulation of rabbit muscle phosphofructokinase by phosphorylation. Biophysical Chemistry, 64(1–3), 199–209. [CrossRef] [PubMed] [Google Scholar]
- Chapman HD, Blake DP. 2022. Genetic selection of Eimeria parasites in the chicken for improvement of poultry health: Implications for drug resistance and live vaccine development. Avian Pathology, 51(6), 521–534. [CrossRef] [PubMed] [Google Scholar]
- Chen XM, Levine SA, Splinter PL, Tietz PS, Ganong AL, Jobin C, Gores GJ, Paya CV, LaRusso NF. 2001. Cryptosporidium parvum activates nuclear factor kappaB in biliary epithelia preventing epithelial cell apoptosis. Gastroenterology, 120(7), 1774–1783. [CrossRef] [PubMed] [Google Scholar]
- Coelho WS, Costa KC, Sola-Penna M. 2007. Serotonin stimulates mouse skeletal muscle 6-phosphofructo-1-kinase through tyrosine-phosphorylation of the enzyme altering its intracellular localization. Molecular Genetics and Metabolism, 92(4), 364–370. [CrossRef] [PubMed] [Google Scholar]
- Cohen P. 2001. The role of protein phosphorylation in human health and disease. The Sir Hans Krebs Medal Lecture. European Journal of Biochemistry, 268(19), 5001–5010. [CrossRef] [PubMed] [Google Scholar]
- Costa Leite T, Da Silva D, Guimarães Coelho R, Zancan P, Sola-Penna M. 2007. Lactate favours the dissociation of skeletal muscle 6-phosphofructo-1-kinase tetramers down-regulating the enzyme and muscle glycolysis. Biochemical Journal, 408(1), 123–130. [CrossRef] [PubMed] [Google Scholar]
- Dapat C, Saito R, Suzuki H, Horigome T. 2014. Quantitative phosphoproteomic analysis of host responses in human lung epithelial (A549) cells during influenza virus infection. Virus Research, 179, 53–63. [CrossRef] [PubMed] [Google Scholar]
- del Cacho E, Gallego M, López-Bernad F, Quílez J, Sánchez-Acedo C. 2004. Expression of anti-apoptotic factors in cells parasitized by second-generation schizonts of Eimeria tenella and Eimeria necatrix. Veterinary Parasitology, 125(3–4), 287–300. [CrossRef] [PubMed] [Google Scholar]
- Delorme-Walker V, Abrivard M, Lagal V, Anderson K, Perazzi A, Gonzalez V, Page C, Chauvet J, Ochoa W, Volkmann N, Hanein D, Tardieux I. 2012. Toxofilin upregulates the host cortical actin cytoskeleton dynamics, facilitating Toxoplasma invasion. Journal of Cell Science, 125(Pt 18), 4333–4342. [PubMed] [Google Scholar]
- Elliott DA, Coleman DJ, Lane MA, May RC, Machesky LM, Clark DP. 2001. Cryptosporidium parvum infection requires host cell actin polymerization. Infection and Immunity, 69(9), 5940–5942. [CrossRef] [PubMed] [Google Scholar]
- Goto A, Matsushita K, Gesellchen V, El Chamy L, Kuttenkeuler D, Takeuchi O, Hoffmann JA, Akira S, Boutros M, Reichhart JM. 2008. Akirins are highly conserved nuclear proteins required for NF-kappaB-dependent gene expression in drosophila and mice. Nature Immunology, 9(1), 97–104. [CrossRef] [PubMed] [Google Scholar]
- Haldar K. 1992. Lipid transport in Plasmodium. Infectious Agents and Disease, 1(5), 254–262. [PubMed] [Google Scholar]
- Han HY, Lin JJ, Zhao QP, Dong H, Jiang LL, Xu MQ, Zhu SH, Huang B. 2010. Identification of differentially expressed genes in early stages of Eimeria tenella by suppression subtractive hybridization and cDNA microarray. Journal of Parasitology, 96(1), 95–102. [CrossRef] [PubMed] [Google Scholar]
- He C, Kong L, Puthiyakunnon S, Wei HX, Zhou LJ, Peng HJ. 2019. iTRAQ-based phosphoproteomic analysis reveals host cell’s specific responses to Toxoplasma gondii at the phases of invasion and prior to egress. Biochimica et Biophysica Acta-Proteins and Proteomics, 1867(3), 202–212. [CrossRef] [Google Scholar]
- He C, Kong L, Zhou L, Xia J, Wei H, Liu M, Peng H. 2017. Host cell vimentin restrains Toxoplasma gondii invasion and phosphorylation of vimentin is partially regulated by interaction with TgROP18. International Journal of Biological Sciences, 13(9), 1126–1137. [CrossRef] [PubMed] [Google Scholar]
- Heasman SJ, Ridley AJ. 2008. Mammalian Rho GTPases: New insights into their functions from in vivo studies. Nature Reviews Molecular Cell Biology, 9(9), 690–701. [CrossRef] [PubMed] [Google Scholar]
- Heussler VT, Machado J Jr, Fernandez PC, Botteron C, Chen CG, Pearse MJ, Dobbelaere DA. 1999. The intracellular parasite Theileria parva protects infected T cells from apoptosis. Proceedings of the National Academy of Sciences, 96(13), 7312–7317. [CrossRef] [PubMed] [Google Scholar]
- Höpker K, Hagmann H, Khurshid S, Chen S, Hasskamp P, Seeger-Nukpezah T, Schilberg K, Heukamp L, Lamkemeyer T, Sos ML, Thomas RK, Lowery D, Roels F, Fischer M, Liebau MC, Resch U, Kisner T, Röther F, Bartram MP, Müller RU, Fabretti F, Kurschat P, Schumacher B, Gaestel M, Medema RH, Yaffe MB, Schermer B, Reinhardt HC, Benzing T. 2012. AATF/Che-1 acts as a phosphorylation-dependent molecular modulator to repress p53-driven apoptosis. EMBO Journal, 31(20), 3961–3975. [CrossRef] [Google Scholar]
- Huang A, Zhang M, Li T, Qin X. 2018. Serum proteomic analysis by Tandem Mass Tags (TMT) based quantitative proteomics in gastric cancer patients. Clinical Laboratory, 64(5), 855–866. [PubMed] [Google Scholar]
- Huang G, Zhang S, Zhou C, Tang X, Li C, Wang C, Tang X, Suo J, Jia Y, El-Ashram S, Yu Z, Cai J, Gupta N, Suo X, Liu X. 2018. Influence of Eimeria falciformis infection on gut microbiota and metabolic pathways in mice. Infection and Immunity, 86(5), e00073–18. [CrossRef] [PubMed] [Google Scholar]
- Joët T, Morin C, Fischbarg J, Louw AI, Eckstein-Ludwig U, Woodrow C, Krishna S. 2003. Why is the Plasmodium falciparum hexose transporter a promising new drug target? Expert Opinion on Therapeutic Targets, 7(5), 593–602. [CrossRef] [PubMed] [Google Scholar]
- Kamada Y, Funakoshi T, Shintani T, Nagano K, Ohsumi M, Ohsumi Y. 2000. Tor-mediated induction of autophagy via an Apg1 protein kinase complex. Yonsei Medical Journal, 150(6), 1507–1513. [Google Scholar]
- Kim JY, Ahn MH, Jun HS, Jung JW, Ryu JS, Min DY. 2006. Toxoplasma gondii inhibits apoptosis in infected cells by caspase inactivation and NF-kappaB activation. Yonsei Medical Journal, 47(6), 862–869. [CrossRef] [PubMed] [Google Scholar]
- Kwon T, Kwon DY, Chun J, Kim JH, Kang SS. 2000. Akt protein kinase inhibits Rac1-GTP binding through phosphorylation at serine 71 of Rac1. Journal of Biological Chemistry, 275(1), 423–428. [CrossRef] [Google Scholar]
- Li Z, Li M, Li X, Xin J, Wang Y, Shen QW, Zhang D. 2018. Quantitative phosphoproteomic analysis among muscles of different color stability using tandem mass tag labeling. Food Chemistry, 249, 8–15. [CrossRef] [PubMed] [Google Scholar]
- Lima TS, Lodoen MB. 2019. Mechanisms of human innate immune evasion by Toxoplasma gondii. Frontiers in Cellular and Infection Microbiology, 9, 103. [CrossRef] [PubMed] [Google Scholar]
- Liu Z, Geng X, Zhao Q, Zhu S, Han H, Yu Y, Huang W, Yao Y, Huang B, Dong H. 2022. Effects of host vimentin on Eimeria tenella sporozoite invasion. Parasites & Vectors, 15(1), 8. [CrossRef] [PubMed] [Google Scholar]
- Man CL, Mu WT, Chang Y, Zhao DX. 2015. Akirin2 expression in response to vaccine-induced immunity in chicken. Genetics and Molecular Research, 14(4), 17489–17495. [CrossRef] [PubMed] [Google Scholar]
- Melo EJ, Carvalho TM, De Souza W. 2001. Behaviour of microtubules in cells infected with Toxoplasma gondii. Biocell. 25(1), 53–59. [PubMed] [Google Scholar]
- Pearson G, Robinson F, Beers Gibson T, Xu BE, Karandikar M, Berman K, Cobb MH. 2001. Mitogen-activated protein (MAP) kinase pathways: regulation and physiological functions. Endocrine Reviews, 22(2), 153–183. [PubMed] [Google Scholar]
- Petry F, Jakobi V, Tessema TS. 2010. Host immune response to Cryptosporidium parvum infection. Experimental Parasitology. 126(3), 304–309. [CrossRef] [PubMed] [Google Scholar]
- Quan JH, Cha GH, Zhou W, Chu JQ, Nishikawa Y, Lee YH. 2013. Involvement of PI 3 kinase/Akt-dependent Bad phosphorylation in Toxoplasma gondii-mediated inhibition of host cell apoptosis. Experimental Parasitology, 133(4), 462–471. [CrossRef] [PubMed] [Google Scholar]
- Ren K, Zhang W, Shi Y, Gong J. 2010. Pim-2 activates API-5 to inhibit the apoptosis of hepatocellular carcinoma cells through NF-kappaB pathway. Pathology & Oncology Research, 16(2), 229–237. [CrossRef] [PubMed] [Google Scholar]
- Rigou P, Piddubnyak V, Faye A, Rain JC, Michel L, Calvo F, Poyet JL. 2009. The antiapoptotic protein AAC-11 interacts with and regulates Acinus-mediated DNA fragmentation. Embo Journal. 28(11), 1576–1588. [CrossRef] [PubMed] [Google Scholar]
- Sharman PA, Smith NC, Wallach MG, Katrib M. 2010. Chasing the golden egg: Vaccination against poultry coccidiosis. Parasite Immunology, 32(8), 590–598. [CrossRef] [PubMed] [Google Scholar]
- Sihag RK, Inagaki M, Yamaguchi T, Shea TB, Pant HC. 2007. Role of phosphorylation on the structural dynamics and function of types III and IV intermediate filaments. ExperimentaL Cell Research, 313(10), 2098–2109. [CrossRef] [PubMed] [Google Scholar]
- Singh V, Ram M, Kumar R, Prasad R, Roy BK, Singh KK. 2017. Phosphorylation: Implications in cancer. Protein Journal, 36(1), 1–6. [CrossRef] [PubMed] [Google Scholar]
- Tartey S, Takeuchi O. 2015. Chromatin remodeling and transcriptional control in innate immunity: Emergence of Akirin2 as a novel player. Biomolecules, 5(3), 1618–1633. [CrossRef] [PubMed] [Google Scholar]
- Tomley F. 1997. Techniques for isolation and characterization of apical organelles from Eimeria tenella sporozoites. Methods, 13(2), 171–176. [CrossRef] [PubMed] [Google Scholar]
- Wang L, Zhu S, Zhao Q, Huang B, Lv L, Liu G, Li Z, Zhao H, Han H, Dong H. 2019. Effects of host fatty acid-binding protein 4 on Eimeria tenella sporozoites invasion of cells. Parasitology Research, 118(6), 1919–1926. [CrossRef] [PubMed] [Google Scholar]
- Wang L, Zhu S, Zhao Q, Huang B, Han H, Liu G, Li Z, Zhao H, Dong H. 2020. Analysis of differentially expressed proteins in CEF cells infected with Eimeria tenella by iTRAQ technology. Chinese Journal of Animal Infectious Diseases, 28, 8–16. [Google Scholar]
- Wang Y, Zhang H. 2019. Regulation of autophagy by mTOR signaling pathway. Advances in Experimental Medicine and Biology, 1206, 67–83. [CrossRef] [PubMed] [Google Scholar]
- Wang Z. 2016. Research on the relationship between the host cell F-actin aggregation and chicken coccidian invasion. South China Agricultural University. [Google Scholar]
- Warncke JD, Beck HP. 2019. Host cytoskeleton remodeling throughout the blood stages of Plasmodium falciparum. Microbiology and Molecular Biology Reviews, 83(4), e00013–e00019. [CrossRef] [PubMed] [Google Scholar]
- Wiśniewski JR, Zougman A, Nagaraj N, Mann M. 2009. Universal sample preparation method for proteome analysis. Nature Methods, 6(5), 359–362. [CrossRef] [PubMed] [Google Scholar]
- Xie MQ, Gilbert JM, Fuller AL, McDougald LR. 1990. A new method for purification of Eimeria tenella merozoites. Parasitology Research, 76(7), 566–569. [CrossRef] [Google Scholar]
- Zhao ZP, Zhu S, Huang B, Han H, Zhao Q, Lv L, Cheng T, Yan M, Dong H. 2018. Preliminary study of RACK1 effects on Eimeria tenella sporozoites invasion of cells. Chinese Journal of Animal Infectious Diseases, 26, 62–68. [Google Scholar]
Cite this article as: Jia L-S, Liu Z, Zhu S-H, Zhao Q-P, Han H-Y, Zhao HZ, Yu Y & Dong H. 2024. Quantitative phosphoproteomic analysis of chicken DF-1 cells infected with Eimeria tenella, using tandem mass tag (TMT) and parallel reaction monitoring (PRM) mass spectrometry. Parasite 31, 23.
All Tables
The number of significantly regulated phosphopeptides/proteins among the three comparison groups.
PRM analysis and comparison of the quantitative results for candidate phosphoproteins.
All Figures
Figure 1 Proportion of serine, threonine, and tyrosine in phosphorylation sites. |
|
In the text |
Figure 2 Volcano plots from different group comparisons. The abscissa indicates difference multiple (logarithmic transformation based on 2), the ordinate indicates the significant of difference (logarithmic transformation based on 10). The red point is significantly upregulated phosphorylated peptide segment, the blue point is significantly downregulated phosphorylated peptide segment and the gray point is a phosphorylated peptide segment with no significant difference. |
|
In the text |
Figure 3 Clustering heatmap of different expression phosphorylated peptides. Each row represents a phosphorylated peptide segment, each column represents a group of samples. The logarithmic value (logarithmic transformation based on 2) of the significantly differentially expressed phosphorylated peptides in different samples is displayed in the clustering heatmap in different colors. Red represents significant up-regulation of phosphorylated peptides; blue represents significant down-regulation of phosphorylated peptides. |
|
In the text |
Figure 4 Gene ontology annotations of the differentially expressed phosphorylated proteins. The abscissa indicates the enriched GO functional classification, including biological process (A), cellular component (B), and molecular function (C). The ordinate indicates the size of the significance of corresponding to each entry, the more left, the more significant. |
|
In the text |
Figure 5 KEGG pathways of the differentially expressed phosphorylated proteins. The abscissa indicates the first 10 significantly enriched KEGG pathways and the ordinate indicates the significance of enriched KEGG pathways, the more left, the more significant. |
|
In the text |
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.